Klassische und ab initio Molekulardynamik-Untersuchungen zu Germaniumdioxidschmelzen

Germaniumdioxid (GeO2) ist ein Glasbildner, der wie das homologe SiO2 ein ungeordnetes tetraedrisches Netzwerk ausbildet. In dieser Arbeit werden mit Hilfe von Molekulardynamik-Computersimulationen die Struktur und Dynamik von GeO2 in Abhängigkeit von der Temperatur untersucht. Dazu werden sowohl Simulationen mit einem klassischen Paarpotentialmodell von Oeffner und Elliott als auch ab initio-Simulationen gemäß der Car-Parrinello-Molekulardynamik (CPMD), bei der elektronische Freiheitsgrade mittels Dichtefunktionaltheorie beschrieben werden, durchgeführt. In der klassischen Simulation werden dazu ein Temperaturen zwischen 6100 K und 2530 K betrachtet. Darüberhinaus ermöglichen Abkühlläufe auf T=300 K das Studium der Struktur des Glases. Zum Vergleich werden CPMD-Simulationen für kleinere Systeme mit 60 bzw. 120 Teilchen bei den Temperaturen 3760 K und 3000 K durchgeführt. In den klassischen Simulationen kann die im Experiment bis 1700 K nachgewiesene, im Vergleich zu SiO2 starke, Temperaturabhängigkeit der Dichte auch bei höheren Temperaturen beobachtet werden. Gute Übereinstimmungen der Simulationen mit experimentellen Daten zeigen sich bei der Untersuchung verschiedener struktureller Größen, wie z.B. Paarkorrelationsfunktionen, Winkelverteilungen, Koordinationszahlen und Strukturfaktoren. Es können leichte strukturelle Abweichungen der CPMD-Simulationen von den klassischen Simulationen aufgezeigt werden: 1. Die Paarabstände in CPMD sind durchweg etwas kleiner. 2. Es zeigt sich, daß die Bindungen in den ab initio-Simulationen weicher sind, was sich auch in einer etwas stärkeren Temperaturabhängigkeit der strukturellen Größen im Vergleich zu den klassischen Simulationen niederschlägt. 3. Für CPMD kann ein vermehrtes Auftreten von Dreierringstrukturen gezeigt werden. 4. In der CPMD werden temperaturabhängige Defektstrukturen in Form von Sauerstoffpaaren beobachtet, die vor allem bei 3760 K, kaum jedoch bei 3000 K auftreten. Alle strukturellen Unterschiede zwischen klassischer und CPMD-Simulation sind eindeutig nicht auf Finite-Size-Effekte aufgrund der kleinen Systemgrößen in den CPMD-Simulationen zurückzuführen, d.h. sie sind tatsächlich methodisch bedingt. Bei der Dynamik von GeO2 wird in den klassischen Simulationen ebenfalls eine gute Übereinstimmung mit experimentellen Daten beobachtet, was ein Vergleich der Diffusionskonstanten mit Viskositätsmessungen bei hohen Temperaturen belegt. Die Diffusionskonstanten zeigen teilweise ein verschiedenes Verhalten zum homologen SiO2. Sie folgen in GeO2 bei Temperaturen unter 3000 K einem Arrheniusgesetz mit einer deutlich niedrigeren Aktivierungsenergie. Darüberhinaus werden die Möglichkeiten der Parametrisierung eines neuen klassischen Paarpotentials mittels der Kräfte entlang der CPMD-Trajektorien untersucht. Es zeigt sich, daß derartige Parametrisierungen sehr stark von den gewählten Startparametern abhängen. Ferner führen sämtliche an die Schmelze parametrisierten Potentiale zu zu hohen Dichten im Vergleich zum Experiment. Zum einen liegt dies sehr wahrscheinlich daran,daß für das System GeO2 Kraftdaten allein nicht ausreichen, um grundlegende strukturelle Größen, wie z.B. Paarkorrelationen und Winkelverteilungen, der CPMD-Simulationen gut reproduzieren zu können. Zum anderen ist wohl die Beschreibung mittels Paarpotentialen nicht ausreichend und es ist erforderlich, Merkörperwechselwirkungen in Betracht zu ziehen.

Ab-initio study of hydrogen technology materials for hydrogen storage and proton conduction

This dissertation deals with two specific aspects of a potential hydrogen-based energy economy, namely the problems of energy storage and energy conversion. In order to contribute to the solution of these problems, the structural and dynamical properties of two promising materials for hydrogen storage (lithium imide/amide) and proton conduction (poly[vinyl phosphonic acid]) are modeled on an atomistic scale by means of first principles molecular dynamics simulation methods.rnrnrnIn the case of the hydrogen storage system lithium amide/imide (LiNH_2/Li_2NH), the focus was on the interplay of structural features and nuclear quantum effects. For these calculations, Path-Integral Molecular Dynamics (PIMD) simulations were used. The structures of these materials at room temperature were elucidated; in collaboration with an experimental group, a very good agreement between calculated and experimental solid-state 1H-NMR chemical shifts was observed. Specifically, the structure of Li_2NH features a disordered arrangement of the Li lattice, which was not reported in previous studies. In addition, a persistent precession of the NH bonds was observed in our simulations. We provide evidence that this precession is the consequence of a toroid-shaped effective potential, in which the protons in the material are immersed. This potential is essentially flat along the torus azimuthal angle, which might lead to important quantum delocalization effects of the protons over the torus.rnrnOn the energy conversion side, the dynamics of protons in a proton conducting polymer (poly[vinyl phosphonic acid], PVPA) was studied by means of a steered ab-initio Molecular Dynamics approach applied on a simplified polymer model. The focus was put on understanding the microscopic proton transport mechanism in polymer membranes, and on characterizing the relevance of the local environment. This covers particularly the effect of water molecules, which participate in the hydrogen bonding network in the material. The results indicate that these water molecules are essential for the effectiveness of proton conduction. A water-mediated Grotthuss mechanism is identified as the main contributor to proton conduction, which agrees with the experimentally observed decay on conductivity for the same material in the absence of water molecules.rnrnThe gain in understanding the microscopic processes and structures present in this materials can help the development of new materials with improved properties, thus contributing to the solution of problems in the implementation of fuel cells.