Refine
Has Fulltext
- yes (1)
Is part of the Bibliography
- yes (1)
Year of publication
- 2010 (1)
Document Type
- Doctoral Thesis (1)
Language
- English (1)
Keywords
- Atomisierungsenergien (1)
- Dichtefunktional (1)
- Quantenchemie (1)
- Reaktionsbarrieren (1)
- Validierung (1)
Institute
Für die Lösung der quantenmechanischen Bewegungsgleichungen, die komplexe, molekulare Systeme beschreiben, sind effiziente und verlässliche Näherungsverfahren erforderlich. Die Dichtefunktionaltheorie (DFT) stellt für die Behandlung der Elektronenwechselwirkung in vielen Fällen den besten Kompromiss zwischen Effizienz und Genauigkeit dar. Im Rahmen der DFT wird die gesamte nicht-klassische Elektron-Elektron-Wechselwirkung im so genannten Austausch-Korrelationsfunktional angenähert. Viele solcher Näherungen sind semi-empirischer Natur, andere wurden ausschließlich von physikalischen Überlegungen abgeleitet. In globalen Hybridfunktionale wird ein konstanter Anteil der integrierten DFT-Austauschenergiedichte durch exakten Austausch aus der Hartree-Fock Näherung ersetzt. Das populärste Funktional B3LYP enthält 20 % exakten Austausch und mehrere empirische Parameter. Der optimale Prozentsatz hängt allerdings sehr stark von den zu berechnenden Systemen und molekularen Eigenschaften ab. Eine Lösung dieses Problems sollten lokale Hybridfunktionale liefern, in denen die Beimischung der exakten Austauschenergiedichte über eine lokale Mischfunktion (LMF) gesteuert wird und daher positions- und molekülabhängig ist. In dieser Arbeit wird ein semi-empirischer Ansatz für die Entwicklung neuer lokaler Hybridfunktionale verfolgt: während die Energiedichten unverändert aus etablierten Näherungen zum Austauschkorrelationsfunktional übernommen werden, stehen parametrisierte LMFs im Zentrum der Untersuchungen. Die verschiedenen LMFs beinhalten neben mindestens einem empirischen Parameter eine Variable die vom Quotienten der von-Weizsäcker kinetischen Energiedichte und der korrelierten kinetischen Energiedichte (sogenannte t-LMFs) bzw. dem reduzierten Dichtegradienten (bezeichnet als t-LMFs) abhängt. Weitere LMFs werden durch zusätzliche Berücksichtigung der Spinpolarisation erhalten. Alle Parameter werden an Atomisierungsenergien bzw. Reaktionsbarrieren bekannter molekularer Testsätze gefittet. Durch Visualisierung der LMFs können zusätzlich Einblicke in den physikalischen Hintergrund und in Möglichkeiten der Weiterentwicklung gewonnen werden. Es wurde beispielsweise beobachtet, dass entlang einer gedehnten Bindung höhere Werte der LMF und damit größere Beimischungen exakter Austauschenergie in Übergangszuständen einhergehen. Dieser Effekt ist für t-LMFs am ausgeprägtesten und korreliert mit besseren Ergebnissen für Reaktionsbarrieren mit lokalen Hybridfunktionalen, die auf einer t-LMF basieren. Bis auf wenige Ausnahmen leiten sich die lokalen Hybridfunktionale in dieser Arbeit aus dem Austausch- und Korrelationsfunktional der lokalen Dichtenäherung (LSDA) ab und enthalten keine Gradientenkorrektur im Sinne der GGA (generalized gradient approximation). Die neuen Funktionale wurden zunächst nicht-selbstkonsistent in eine Entwicklerversion des quantenchemischen Programmpaketes Turbomole implementiert. Das bedeutet, für gegebene Molekülorbitale bzw. eine gegeben Elektronendichte kann lediglich die Gesamtenergie berechnet werden. Dies ist eine anerkannte Näherung, die vor allem für die Optimierung der Parameter eine große Zeitersparnis darstellt. Um letztlich orbitalabhängige, molekulare Eigenschaften berechnen zu können wird neben der Gesamtenergie auch noch das zugehörige lokale Hybridpotential benötigt. Für die selbstkonsistente Implementierung wird die funktionale Ableitung der Austauschkorrelationsenergie nach den Orbitalen bestimmt. Daraus resultierend müssen neben den üblichen lokalen Austauschkorrelationspotentialtermen auch Integrale berechnet werden, die das mit der LMF gewichtete nicht-lokale exakte Austauschpotential enthalten. Die entsprechenden Terme kann man, genauso wie die exakte Austauschenergiedichte an sich, nicht analytisch berechnen. Früheren Ansätzen folgend wurden sie in der vorliegenden Arbeit in einer Basissatzentwicklung angenähert, wobei der Einfachheit halber die atomaren Basisfunktionen verwendet wurden. Um die Genauigkeit dieser sogenannten RI (resolution of the identity)-Näherung validieren zu können und auch schon im Hinblick auf die Anpassung einer Hilfsbasis, wurde darüber hinaus die numerische Berechnung aller Integrale, die das exakte Austauschpotential und die entsprechende Energiedichte enthalten, implementiert. Unter Verwendung der RI-Näherung ist der Rechenaufwand lokaler Hybride vergleichbar mit dem globaler Hybridfunktionale: Während die formale Skalierung in Abhängigkeit der Systemgröße gleich ist, ergab sich ein etwas höherer Vorfaktor für die lokalen Hybride. Verschiedene Literaturbekannte Testsätze mit Atomisierungsenergien, Reaktionsbarrieren, Dissoziationsenergien oder Gleichgewichtsabständen, die teilweise einige Schwächen bisheriger Dichtefunktionalnäherungen aufdecken, wurden berücksichtigt. Für die 223 Atomisierungsenergien des G3 Testsatzes stellen alle unsere Funktionale eine signifikante Verbesserung gegenüber B3LYP dar. Atomisierungsenergien sind insofern ein sensibler Test, da alle Bindungen gebrochen werden und Fehlerkompensation eine untergeordnete Rolle spielt. Vor allem lokale Hybridfunktionale, deren LMFs neben der kinetischen Energiedichte explizit von der Spinpolarisation abhängen, lieferten hervorragende Resultate. Obwohl im Vergleich zu Atomisierungsenergien für die korrekte Berechnung von Reaktionsbarrieren im Allgemeinen mehr exakter Austausch benötigt wird, sind unsere Funktionale auch für zwei Testsätze mit jeweils 38 Reaktionsbarrieren besser als B3LYP. Zwar kann mit einem globalen Hybrid mit 50 % exaktem Austausch eine geringere Abweichung von den Richtwerten erzielt werden, aber ein solches Funktional ist für thermochemische Daten unzureichend. Hier wurde erstmals gezeigt, dass lokale Hybridfunktionale ohne Gradientenkorrektur sowohl für Thermochemie als auch für Kinetik zufrieden stellende Ergebnisse liefern können. Das Dissoziationsverhalten symmetrischer Radikalkationen stellt für die hier diskutierten Dichtefunktionale nach wie vor eine Herausforderung dar: Die Dissoziationsenergien von sieben Modellsystemen werden mit unseren Funktionalen stark überschätzt und Gleichgewichtsabstände unterschätzt. Insgesamt sind die Werte nur marginal besser als mit B3LYP. Neben Eigenschaften von Hauptgruppenverbindungen wurden zudem Übergangsmetalldimere und -monohydride untersucht. Für erstere ist eine gute Beschreibung dynamischer sowie statischer Elektronenkorrelation ausschlaggebend. In den Hydriden andererseits dominiert mit gängigen Dichtefunktionalen die unphysikalische Selbstwechselwirkung eines Elektrons mit sich selbst. Für die 3d-Übergangsmetalldimere sind die getesteten Funktionale genauso gut wie B3LYP und für die Hydride etwas besser. Atomare s-d Transferenergien von 3d Übergangsmetallen verbleiben auch für unsere lokalen Hybridfunktionale, die insgesamt schlechtere Ergebnisse erzielen als B3LYP, noch problematisch. Das hierfür geeignetste lokale Hybridfunktional basiert auf einer s-LMF und beinhaltet LYP Korrelation. Für die isotropen Hyperfeinkopplungskonstanten (HFCCs) kleiner Hauptgruppenverbindungen wurden zufriedenstellende Ergebnisse (ähnlich wie B3LYP) mit einem t-LMF basierten lokalen Hybrid erzielt. Die RI Näherung zum lokalen Hybridpotential wurde dem numerisch exakten Potential für die Berechnung von Gesamtenergien, isotrope HFCCs und Orbitalenergien für verschiedene Basissätze gegenübergestellt. Wie erwartet ist der Fehler für Gesamtenergien mit der RI-Näherungen vergleichsweise gering, vor allem relativ zu den verbleibenden Abweichungen von experimentellen Energien. Der Vergleich der mittleren absoluten Abweichung von experimentellen Werten für 26 isotrope HFCCs zeigt sogar für mittelgroße und kontrahierte IGLO Basissätze nur geringe Unterschiede zwischen dem RI-Potential und dem numerisch exakten lokalen Hybridpotential. Die Analyse der HFCCs einzelner Moleküle und der Orbitalenergien des CN Moleküls offenbart allerdings, dass Ungenauigkeiten aufgrund der RI-Näherung hier eine größere Rolle spielen, vor allem wenn zu kleine atomare Basissätze verwendet werden. Von den untersuchten lokalen Hybriden stellen sich einige als hervorragende Kandidaten für die Berechnung thermochemischer und kinetischer Eigenschaften heraus. Jeweils unterschiedliche Funktionale erzielen darüber hinaus mit den besten bekannten Funktionalen vergleichbare Ergebnisse für isotrope Hyperfeinkopplungskonstanten und ausgewählte Eigenschaften kleiner Übergangsmetallverbindungen. Die in dieser Arbeit präsentierten lokalen Hybridfunktionale stellen daher einen wichtigen Schritt in der Entwicklung universeller Näherungen zum Austauschkorrelationsfunktional dar. Zur akkuraten Beschreibung molekularer Eigenschaften von Übergangsmetallkomplexen und dem Dissoziationsverhalten von Radikal-Kation-Dimeren neben Thermochemie und Kinetik, werden in Zukunft wohl komplexere LMFs benötigt. Um konkurrenzfähige lokale Hybride mit gradientenkorrigierter Austausch- und Korrelationsenergiedichte zu entwickeln, müssen darüber hinaus weitere Studien zum Einfluss des abweichenden Eichursprungs der miteinander kombinierten Austauschenergiedichten durchgeführt werden. Eine andere Möglichkeit ist die Entwicklung speziell abgestimmter Korrelationsfunktionale für lokale Hybride. Außerdem sollte die Qualität der RI-Näherung zum lokalen Hybridpotential detaillierter untersucht werden. Hierfür könnten zum Beispiel Ionisierungsenergien und Elektronenaffinitäten herangezogen werden. Um zusätzliche Abweichungen oder sogar fälschlicherweise "zu gute" Ergebnisse bei Validierungsrechnungen zu vermeiden, sollten Hilfsbasen für die Entwicklung des nicht-lokalen exakten Austauschpotentials implementiert und optimiert werden. Einer der nächsten Implementierungsschritte sollte auch Gradienten bezüglich der Kernkoordinaten beinhalten, um die Validierung der neuen lokalen Hybridfunktionale auf Strukturoptimierungen auszuweiten.