Novel algorithms for electronic structure based molecular dynamics with linear system-size scaling

Die vorliegende Arbeit behandelt die Entwicklung und Verbesserung von linear skalierenden Algorithmen für Elektronenstruktur basierte Molekulardynamik. Molekulardynamik ist eine Methode zur Computersimulation des komplexen Zusammenspiels zwischen Atomen und Molekülen bei endlicher Temperatur. Ein entscheidender Vorteil dieser Methode ist ihre hohe Genauigkeit und Vorhersagekraft. Allerdings verhindert der Rechenaufwand, welcher grundsätzlich kubisch mit der Anzahl der Atome skaliert, die Anwendung auf große Systeme und lange Zeitskalen. Ausgehend von einem neuen Formalismus, basierend auf dem großkanonischen Potential und einer Faktorisierung der Dichtematrix, wird die Diagonalisierung der entsprechenden Hamiltonmatrix vermieden. Dieser nutzt aus, dass die Hamilton- und die Dichtematrix aufgrund von Lokalisierung dünn besetzt sind. Das reduziert den Rechenaufwand so, dass er linear mit der Systemgröße skaliert. Um seine Effizienz zu demonstrieren, wird der daraus entstehende Algorithmus auf ein System mit flüssigem Methan angewandt, das extremem Druck (etwa 100 GPa) und extremer Temperatur (2000 - 8000 K) ausgesetzt ist. In der Simulation dissoziiert Methan bei Temperaturen oberhalb von 4000 K. Die Bildung von sp²-gebundenem polymerischen Kohlenstoff wird beobachtet. Die Simulationen liefern keinen Hinweis auf die Entstehung von Diamant und wirken sich daher auf die bisherigen Planetenmodelle von Neptun und Uranus aus. Da das Umgehen der Diagonalisierung der Hamiltonmatrix die Inversion von Matrizen mit sich bringt, wird zusätzlich das Problem behandelt, eine (inverse) p-te Wurzel einer gegebenen Matrix zu berechnen. Dies resultiert in einer neuen Formel für symmetrisch positiv definite Matrizen. Sie verallgemeinert die Newton-Schulz Iteration, Altmans Formel für beschränkte und nicht singuläre Operatoren und Newtons Methode zur Berechnung von Nullstellen von Funktionen. Der Nachweis wird erbracht, dass die Konvergenzordnung immer mindestens quadratisch ist und adaptives Anpassen eines Parameters q in allen Fällen zu besseren Ergebnissen führt.

Oberflächenuntersuchungen von Metalla-Tetraazaporphyrinen

Ziel dieser Arbeit war es, ein System zu entwickeln, in dem ein durch Licht induzierter Elektronentransfer stattfinden kann. Dazu wurden ein Kupfer(II)- und ein Zink(II)Tetraazaporphyrin mit acht 4-tert-Butylphenyl-Substituenten synthetisiert (Cu4Dinit, Zn4Dinit). Die Energielücke von 1,85 eV zwischen HOMO und LUMO von Cu4Dinit in Lösung wurde mit Hilfe von Cyclovoltammetrie und UV/Vis-Messungen bestimmt. Somit ist sie größer als für Cu4Dinit Moleküle, die auf einer Oberfläche (Wolfram(100)) liegen und mit STM-, STS-Messungen untersucht wurden. Hier beträgt die Energielücke 1,35 eV, was durch eine Drehung der Phenylringe in die Ebene der Pyrrolringe des Makrozyklus und somit durch eine bessere Überlappung der Orbitale erklärt werden kann. Um die Wechselwirkung der Moleküle mit der Oberfläche zu untersuchen, wurde Cu4Dinit, wie oben beschrieben, auf Magnetit aufgedampft. Dadurch wurde ausschließlich die Wechselwirkung zwischen den Elektronenspins des Kupfer(II)-ions und den Elektronenspins des Eisens im Magnetit betrachtet. Durch Messungen der Röntgenabsorption und des XMCD-Effektes konnten das Spinmoment, Bahnmoment und das Gesamtmoment des Kupfers berechnet und eine anisotrope Kopplung des Elektronenspins des Kupferions zum Magnetit, in Abhängigkeit der Magnetisierungsrichtung des Magnetits, festgestellt werden. Wenn der Magnetit senkrecht zur Oberfläche (out-of-plane) magnetisiert ist, ist die Kopplung ferromagnetisch, während bei einer Magnetisierungsrichtung parallel zur Ebene (in-plane) des Magnetits der Elektronenspin des Kupfers antiferromagnetisch mit dem des Eisens koppelt. Dadurch muss der Hamiltonian, der die Wechselwirkung zwischen zwei Spins beschreibt, bei einer anisotropen Kopplung um einen ansiotropen Term ergänzt werden. Das Ergebnis, dass der Elektronenspin des Kupferions durch die Richtung der Magnetisierung des Magnetits beeinflusst werden kann, eröffnet neue Wege, um die Spinkonfiguration von auf der Oberfläche liegenden Molekülen mit ungepaarten Elektronen, wie die zentralen Metallionen der Makrozyklen aber auch die Elektronenspins anderer metallorganischer Komplexe oder molekulare Magnete, durch ein externes Magnetfeld zu beeinflussen. rnDurch die stöchiometrische Templatreaktion von Pyrazino[2,3-f][1,10]-phenanthrolin-2,3-di-carbonitril (Dicnq), Bis(4-tert-Butylphenyl)-fumarodinitril (Dinit) und Kupfer(II)-acetat wurde eine Koordinationsmöglichkeit für ein Ruthenium(II)-ion in einem Tetraazaporphyrin hergestellt und so die Makrozyklen Cu3Dinit1Dicnq und Zn3Dinit1Dicnq synthetisiert, mit Rutheniumionen versetzt und ebenfalls mit Hilfe von Röntgenabsorptionsmessungen und XMCD untersucht. Durch die Vergleiche mit Zn3Dinit1Dicnq und den jeweiligen Verbindungen mit koordinierten Rutheniumionen (Cu3Dinit1Dicnq-1Ru, Zn3Dinit1Dicnq-1Ru) konnte gezeigt werden, dass eine Verschiebung der Elektronendichte des Rutheniumions zu dem zentralen Kupferion des Makrozyklus stattgefunden hat und durch die Koordination eines Rutheniumions in der Peripherie des Tetraazaporphyrins die energetische Lage der Kupferorbitale beeinflusst wird.rnDer Einfluss von vier koordinierten Ruthenium(II)-ionen auf das zentrale Kupferion wurde an Hand des in dieser Arbeit hergestellten Kupfer(II)phenanthralocyanins (Cu4Dicnq) untersucht, das aus vier Dicnq-Liganden und Kupfer(II)-acetat synthetisiert wurde. Auf Grund der schlechten Löslichkeit wurde für die Koordination der Rutheniumionen der Prekursor [Ru(bipy)2Dicnq](PF6)2 hergestellt und daraus der Makrozyklus in einer Templatsynthese mit Kufper(II)-ionen gebildet. Durch diese neue Syntheseroute war es möglich, die Verbindung Cu4Dicnq-4Ru herzustellen und ebenfalls durch Röntgenabsorption und XMCD zu untersuchen und so das Spin- und Bahnmoment zu ermitteln. Ein Teil der Elektronendichte des Rutheniumions in dieser Verbindung wird auf die zusätzlich an das Rutheniumion koordinierten 2,2'-Bipyridine und nicht auf den Makrozyklus, wie in Cu3Dinit1Dicnq-1Ru, geschoben. Trotzdem konnte die Funktionsweise als Modell des Photosystems II durch eine Oxidation durch die Bestrahlung mit einer Quecksilberlampe mit para-Benzochinon beobachtet werden. Dies bestätigte die Funktionsweise des Kupfer(II)phenanthralocyanins mit koordinierten Rutheniumionen, da ein durch Licht induzierter Elektronenübergang auf das para-Benzochinon stattgefunden hat.rn