350 likes | 477 Views
Výpočet konstant hyperjemného štěpení β-protonů v EPR spektrech substituovaných kationradikálů pyrrolu. Petr Holub. Cíle diplomové práce. Navrhnout metodiku pro výpočet HJS flexibilních organických molekul Ověřit metodiku pomocí výpočtů HJS na skupině různě substituovaných pyrrolů.
E N D
Výpočet konstant hyperjemného štěpení β-protonů v EPR spektrech substituovaných kationradikálů pyrrolu Petr Holub
Cíle diplomové práce • Navrhnout metodiku pro výpočet HJS flexibilních organických molekul • Ověřit metodiku pomocí výpočtů HJS na skupině různě substituovaných pyrrolů
Interakce spinů jádra a elektronu • Magnetický moment • Hamiltonián interakce • Štěpící konstanta = izotropická část interakce
Mechanismy vzniku spinových hustot na jádrech vodíků • mechanismus spinové polarizace • mechanismus přenosu přes prostor(„through space“) • interakce je přímo úměrná překryvu orbitalů
Metodika • Výběr metody pro optimalizaci geometrie • Výběr metody pro výpočet spinových hustot • Boltzmannovsky vážený průměr spinových hustot • Průměrování štěpících konstant methylů • Test solvatačních modelů
Výběr teoretického modelu • Výběr modelu pro optimalizaci geometrie • musí být dostatečně rychlý (relaxovaný PES scan) • musí dávat dobré výsledky pro běžné organické molekuly => unrestricted AM1
Výběr metody pro výpočet spinových hustot • přiměřeně náročná metoda • DFT funkcionály dávají velmi dobré výsledky • použití vhodné báze =>unrestricted B3LYP/D95(d,p)
Problematika bází při výpočtu spinových hustot na jádrech 1 STO vs. lineárníkombinace 5 GTO 1 STO vs. 1 GTO
Pokud jsou pohyby molekul mimo časovou škálu EPR => boltzmannovsky vážené průměry
Vlastnosti funkce a cos2(f) • Možno prokázat následující • Zjednodušení výpočtů volně rotujících methylových skupin – stačí jediná konformace
Vývoj vlastního dávkového systému nad Gaussian 98 • Zjednodušení zadávání vstupních dat • Nezávislý relaxovaný PES scan • Zjednodušení a automatizace práce velkými objemy výstupních dat • Vhodná paralelizace úlohy - cluster-buster aplikace
Ověření metodiky výpočty • Výpočet štěpící konstanty H atomů v radikálu CH3 (2,304 mT) • vynikající výsledky DFT funkcionálů, obzvláště B3LYP (stabilita výsledků přes většinu bází) • vhodná Dunningova báze D95(d,p)
2,5-dimethylpyrrol • Ověření závislosti štěpící konstranty a cos2(f) na methylových H atomech(nejsou zde další interakce) • Výpočet štěpících konstant navrženou metodikou • Testy různých solvatačních modelů (nejmenší do výpočtů zahrnutý pyrrol)
2,5-diethyl-3,4-dimethylpyrrol • Při výpočtu bez substituentu na N vyšla štěpící konstanta CH2 skupiny Et více jak 2x větší a CH3 Et více jak dvakrát menší oproti experimentu • Experimentální hodnoty:CH2 Et 0,771 mTCH3 Et 0,052 mTCH3 Me 0,357 mTN 0,424 mT
Po přidání substituentu se výpočet velmi dobře shoduje s experimentem
2,3,4,5-bis-(trimethylen)pyrrol • správný IUPAC název:2,3,4,5,6,7-Hexahydro-1H-dicyclopenta[b,d]pyrrole • RTG analýza struktury • selhání metody AM1 při optimalizaci geometrie
RTG strukturní analýza (1) • elementární buňka obsahuje 4 molekuly • dvě struktury se od sebe mírně odlišují
RTG strukturní analýza (3) • Úhel „vyklopení“ trimethylenové skupiny se pohybuje od 12 do 22 stupňů, což odpovídá výsledkům metody B3LYP/6-31G(d,p)
Další výpočty • Další počítané molekuly:2,3,4,5-tetramethylpyrrol2,5-diisopropyl-3,4-dimethylpyrrol3,4-difenyl-2,5-dimethylpyrrol
Závěr • Navržená metodika • Optimalizace geometrie na úrovni AM1 • Výpočet spinových hustot na úrovni B3LYP/D95(d,p) • Použití open-shell formalismu pro všechny výpočty • Výpočet průměrné hodnoty vážené boltzmannovskými faktory • Zjednodušení výpočtu v případě volně rotujících methylových skupin
Ověření výpočty (typická chyba výpočtů je 5-10%) • Ověření metod na methylovém radikálu • Ověření závislosti a cos2(f) • Vysvětlení HJS u 2,5-diethyl-3,4-dimethylpyrrolu • Výpočet bariéry překlápění a RTG analýza struktury 2,3,4,5-bis-trimethylenpyrrolu DALŠÍ INFORMACE... Dokumenty uvedené na konci této presentace
Poděkování • Doc. Pavel Kubáček • Doc. Luděk Matyska a SCB • Mgr. Markéta Ludmila Munzarová • Prof. Zdirad Žák • a všem posluchačům za pozornost
Odpovídající dokumenty • http://cheminfo.chemi.muni.cz/~hopet/svoc/diplomka.pdf • http://cheminfo.chemi.muni.cz/~hopet/svoc/pyroly.pdf Q&A’s...