Quantum effects and dynamics in hydrogen-bonded systems: a first-principles approach to spectroscopic experiments

Computer simulations have become an important tool in physics. Especially systems in the solid state have been investigated extensively with the help of modern computational methods. This thesis focuses on the simulation of hydrogen-bonded systems, using quantum chemical methods combined with molecular dynamics (MD) simulations. MD simulations are carried out for investigating the energetics and structure of a system under conditions that include physical parameters such as temperature and pressure. Ab initio quantum chemical methods have proven to be capable of predicting spectroscopic quantities. The combination of these two features still represents a methodological challenge. Furthermore, conventional MD simulations consider the nuclei as classical particles. Not only motional effects, but also the quantum nature of the nuclei are expected to influence the properties of a molecular system. This work aims at a more realistic description of properties that are accessible via NMR experiments. With the help of the path integral formalism the quantum nature of the nuclei has been incorporated and its influence on the NMR parameters explored. The effect on both the NMR chemical shift and the Nuclear Quadrupole Coupling Constants (NQCC) is presented for intra- and intermolecular hydrogen bonds. The second part of this thesis presents the computation of electric field gradients within the Gaussian and Augmented Plane Waves (GAPW) framework, that allows for all-electron calculations in periodic systems. This recent development improves the accuracy of many calculations compared to the pseudopotential approximation, which treats the core electrons as part of an effective potential. In combination with MD simulations of water, the NMR longitudinal relaxation times for 17O and 2H have been obtained. The results show a considerable agreement with the experiment. Finally, an implementation of the calculation of the stress tensor into the quantum chemical program suite CP2K is presented. This enables MD simulations under constant pressure conditions, which is demonstrated with a series of liquid water simulations, that sheds light on the influence of the exchange-correlation functional used on the density of the simulated liquid.

Highly accurate quantum chemistry : spin-orbit splittings via multireference coupled-cluster methods and applications in heavy-atom main-group chemistry

Diese Dissertation demonstriert und verbessert die Vorhersagekraft der Coupled-Cluster-Theorie im Hinblick auf die hochgenaue Berechnung von Moleküleigenschaften. Die Demonstration erfolgt mittels Extrapolations- und Additivitätstechniken in der Single-Referenz-Coupled-Cluster-Theorie, mit deren Hilfe die Existenz und Struktur von bisher unbekannten Molekülen mit schweren Hauptgruppenelementen vorhergesagt wird. Vor allem am Beispiel von cyclischem SiS_2, einem dreiatomigen Molekül mit 16 Valenzelektronen, wird deutlich, dass die Vorhersagekraft der Theorie sich heutzutage auf Augenhöhe mit dem Experiment befindet: Theoretische Überlegungen initiierten eine experimentelle Suche nach diesem Molekül, was schließlich zu dessen Detektion und Charakterisierung mittels Rotationsspektroskopie führte. Die Vorhersagekraft der Coupled-Cluster-Theorie wird verbessert, indem eine Multireferenz-Coupled-Cluster-Methode für die Berechnung von Spin-Bahn-Aufspaltungen erster Ordnung in 2^Pi-Zuständen entwickelt wird. Der Fokus hierbei liegt auf Mukherjee's Variante der Multireferenz-Coupled-Cluster-Theorie, aber prinzipiell ist das vorgeschlagene Berechnungsschema auf alle Varianten anwendbar. Die erwünschte Genauigkeit beträgt 10 cm^-1. Sie wird mit der neuen Methode erreicht, wenn Ein- und Zweielektroneneffekte und bei schweren Elementen auch skalarrelativistische Effekte berücksichtigt werden. Die Methode eignet sich daher in Kombination mit Coupled-Cluster-basierten Extrapolations-und Additivitätsschemata dafür, hochgenaue thermochemische Daten zu berechnen.

Higher-order molecular properties and excitation energies in single-reference and multireference coupled-cluster theory

Coupled-cluster (CC) theory is one of the most successful approaches in high-accuracy quantum chemistry. The present thesis makes a number of contributions to the determination of molecular properties and excitation energies within the CC framework. The multireference CC (MRCC) method proposed by Mukherjee and coworkers (Mk-MRCC) has been benchmarked within the singles and doubles approximation (Mk-MRCCSD) for molecular equilibrium structures. It is demonstrated that Mk-MRCCSD yields reliable results for multireference cases where single-reference CC methods fail. At the same time, the present work also illustrates that Mk-MRCC still suffers from a number of theoretical problems and sometimes gives rise to results of unsatisfactory accuracy. To determine polarizability tensors and excitation spectra in the MRCC framework, the Mk-MRCC linear-response function has been derived together with the corresponding linear-response equations. Pilot applications show that Mk-MRCC linear-response theory suffers from a severe problem when applied to the calculation of dynamic properties and excitation energies: The Mk-MRCC sufficiency conditions give rise to a redundancy in the Mk-MRCC Jacobian matrix, which entails an artificial splitting of certain excited states. This finding has established a new paradigm in MRCC theory, namely that a convincing method should not only yield accurate energies, but ought to allow for the reliable calculation of dynamic properties as well. In the context of single-reference CC theory, an analytic expression for the dipole Hessian matrix, a third-order quantity relevant to infrared spectroscopy, has been derived and implemented within the CC singles and doubles approximation. The advantages of analytic derivatives over numerical differentiation schemes are demonstrated in some pilot applications.

Experimental investigation into the effect of stress on dissolution and growth of very soluble brittle salts in aqueous solution

Das Studium der Auflösungs- und Wachstumsprozesse an Feststoff-Flüssigkeits-Grenzflächen unter nicht-hydrostatischen Beanspruchungen ist wesentlich für das Verständnis von Defor-mationsprozessen, die in der Erde ablaufen. Unter diesen genannten Prozessen gehört die Drucklösung zu den wichtigsten duktilen Deformationsprozessen, von der Diagenese bishin zur niedrig- bis mittelgradigen metamorphen Bedingungen. Bisher ist allerdings wenig darüber bekannt, welche mechanischen, physikalischen oder chemischen Potentialenergie-Gradienten die Drucklösung steuern. I.a. wird angenommen, daß die Drucklösung durch Un-terschiede kristallplastischer Verformungsenergien oder aber durch Unterschiede der Normal-beanspruchung an Korngrenzen gesteuert wird. Unterschiede der elastischen Verformungs-energien werden dabei allerdings als zu gering erachtet, um einen signifikanten Beitrag zu leisten. Aus diesem Grund werden sie als mögliche treibende Kräfte für die Drucklösung vernachlässigt. Andererseits haben neue experimentelle und theoretische Untersuchungen gezeigt, daß die elastische Verformung in der Tat einen starken Einfluß auf Lösungs- und Wachstumsmechanismen von Kristallen in einer Lösung haben kann. Da die in der Erdkruste vorherrschenden Deformationsmechanismen überwiegend im elastischen Verformungsbereich der Gesteine ablaufen, ist es sehr wichtig, das Verständnis für die Effekte, die die elastische Verformung verursacht, zu erweitern, und ihre Rolle während der Deformation durch Drucklösung zu definieren. Die vorliegende Arbeit beschäftigt sich mit Experimenten, bei denen der Effekt der mechanisch kompressiven Beanspruchung auf Lösungs- und Wachstumsprozesse von Einzelkristallen unterschiedlicher, sehr gut löslicher, elastisch/spröder Salze untersucht wurde. Diese Salze wurden als Analoga gesteinsbildender Minerale wie Quarz und Calcit ausgewählt. Der Einfluß von Stress auf die Ausbildung der Oberflächenmikrostrukturen in einer untersättigten Lösung wurde an Kaliumalaun untersucht.Lösungsrillen (20 – 40 µm breit, 10 – 40 µm tief und 20 – 80 µm Abstand) entwickelten sich in den Bereichen, in denen die Beanspruchung im Kristall am größten war. Sie verschwanden wieder, sobald der Kristall entlastet wurde. Diese Rillen entwickelten sich parallel zu niedrig indizierten kristallographischen Richtungen und sub-perpendikular zu den Trajektorien, die der maximalen, lokalen kompressiven Beanspruchung entsprachen. Die Größe der Lösungsrillen hing von der lokalen Oberflächenbeanspruchung, der Oberflächenenergie und dem Untersättigungsgrad der wässrigen Lösung ab. Die mikrostrukturelle Entwicklung der Kristalloberflächen stimmte gut mit den theoretischen Vorhersagen überein, die auf den Modellen von Heidug & Leroy (1994) und Leroy & Heidug (1994) basieren. Der Einfluß der Beanspruchung auf die Auflösungsrate wurde an Natriumchlorat-Einzelkristallen untersucht. Dabei wurde herausgefunden, daß sich gestresste Kristalle schneller lösen als Kristalle, auf die keine Beanspruchung einwirkt. Der experimentell beobachtete Anstieg der Auflösungsrate der gestressten Kristalle war ein bis zwei Größenordnungen höher als theoretisch erwartet. Die Auflösungsrate stieg linear mit dem Stress an, und der Anstieg war um so größer, je stärker die Lösung untersättigt war. Außerdem wurde der Effekt der Bean-spruchung auf das Kristallwachstum an Kaliumalaun- und Kaliumdihydrogenphosphat-Ein-zelkristallen untersucht. Die Wachstumsrate der Flächen {100} und {110} von Kalium-alaun war bei Beanspruchung stark reduziert. Für all diese Ergebnisse spielte die Oberflächenrauhigkeit der Kristalle eine Schlüsselrolle, indem sie eine nicht-homogene Stressverteilung auf der Kristalloberfläche verursachte. Die Resultate zeigen, daß die elastische Verformung eine signifikante Rolle während der Drucklösung spielen kann, und eine signifikante Deformation in der oberen Kruste verursachen kann, bei Beanspruchungen, die geringer sind, als gemeinhin angenommen wird. Somit folgt, daß die elastische Bean-spruchung berücksichtigt werden muß, wenn mikrophysikalische Deformationsmodelle entwickelt werden sollen.

Ternäre polymerhaltige Lösungen - Phasenverhalten und Grenzflächenspannung

Im Rahmen dieser Arbeit wurden experimentelle und theoretische Untersuchungen zum Phasen- und Grenzflächenverhalten von ternären Systemen des Typs Lösungsmittel/Fällungsmittel/Polymer durchgeführt. Diese Art der Mischungen ist vor allem für die Planung und Durchführung der Membranherstellung von Bedeutung, bei der die genaue Kenntnis des Phasendiagramms und der Grenzflächenspannung unabdingbar ist. Als Polymere dienten Polystyrol sowie Polydimethylsiloxan. Im Fall des Polystyrols kam Butanon-2 als Lösungsmittel zum Einsatz, wobei drei niedrigmolekulare lineare Alkohole als Fällungsmittel verwendet wurden. Für Polydimethylsiloxan eignen sich Toluol als Lösungsmittel und Ethanol als Fällungsmittel. Durch Lichtstreumessungen, Dampfdruckbestimmungen mittels Headspace-Gaschromatographie (VLE-Gleichgewichte) sowie Quellungsgleichgewichten lassen sich die thermodynamischen Eigenschaften der binären Subsysteme charakterisieren. Auf Grundlage der Flory-Huggins-Theorie kann das experimentell bestimmte Phasenverhalten (LLE-Gleichgewichte) in guter Übereinstimmung nach der Methode der Direktminimierung der Gibbs'schen Energie modelliert werden. Zieht man die Ergebnisse der Aktivitätsbestimmung von Dreikomponenten-Mischungen mit in Betracht, so ergeben sich systematische Abweichungen zwischen Experiment und Theorie. Sie können auf die Notwendigkeit ternärer Wechselwirkungsparameter zurückgeführt werden, die ebenfalls durch Modellierung zugänglich sind.Durch die aus den VLE- und LLE-Untersuchungen gewonnenen Ergebnissen kann die sog. Hump-Energie berechnet werden, die ein Maß für die Entmischungstendenz darstellt. Diese Größe eignet sich gut zur Beschreibung von Grenzflächenphänomenen mittels Skalengesetzen. Die für binäre Systeme gefundenen theoretisch fundierten Skalenparameter gelten jedoch nur teilweise. Ein neues Skalengesetz lässt erstmals eine Beschreibung über die gesamte Mischungslücke zu, wobei ein Parameter durch eine gemessene Grenzflächenspannung (zwischen Fällungsmittel/Polymer) ersetzt werden kann.

Development of a comprehensive on-line multidimensional high performance column liquid chromatography system for protein and peptide mapping with integrated size selective sample fractionation

Aufbau einer kontinuierlichen, mehrdimensionalen Hochleistungs-flüssigchromatographie-Anlage für die Trennung von Proteinen und Peptiden mit integrierter größenselektiver ProbenfraktionierungEs wurde eine mehrdimensionale HPLC-Trennmethode für Proteine und Peptide mit einem Molekulargewicht von <15 kDa entwickelt.Im ersten Schritt werden die Zielanalyte von höhermolekularen sowie nicht ionischen Bestandteilen mit Hilfe von 'Restricted Access Materialien' (RAM) mit Ionenaustauscher-Funktionalität getrennt. Anschließend werden die Proteine auf einer analytischen Ionenaustauscher-Säule sowie auf Reversed-Phase-Säulen getrennt. Zur Vermeidung von Probenverlusten wurde ein kontinuierlich arbeitendes, voll automatisiertes System auf Basis unterschiedlicher Trenngeschwindigkeiten und vier parallelen RP-Säulen aufgebaut.Es werden jeweils zwei RP-Säulen gleichzeitig, jedoch mit zeitlich versetztem Beginn eluiert, um durch flache Gradienten ausreichende Trennleistungen zu erhalten. Während die dritte Säule regeneriert wird, erfolgt das Beladen der vierte Säule durch Anreicherung der Proteine und Peptide am Säulenkopf. Während der Gesamtanalysenzeit von 96 Minuten werden in Intervallen von 4 Minuten Fraktionen aus der 1. Dimension auf die RP-Säulen überführt und innerhalb von 8 Minuten getrennt, wobei 24 RP-Chromatogramme resultieren.Als Testsubstanzen wurden u.a. Standardproteine, Proteine und Peptide aus humanem Hämofiltrat sowie aus Lungenfibroblast-Zellkulturüberständen eingesetzt. Weiterhin wurden Fraktionen gesammelt und mittels MALDI-TOF Massenspektrometrie untersucht. Bei einer Injektion wurden in den 24 RP-Chromatogrammen mehr als 1000 Peaks aufgelöst. Der theoretische Wert der Peakkapazität liegt bei ungefähr 3000.

Supramolekulare Chemie mit Cyclodextrinen

Zusammenfassung Die vorliegende Arbeit gliedert sich in zwei Teilgebiete, die beide unter dem Gesichtspunkt der Supramolekularen Chemie mit Cyclodextrinen bearbeitet wurden. Der erste Teil befasste sich mit der Synthese und Charakterisierung von chemisch modifizierten Cyclodextrinen, die flüssigkristalline Eigenschaften besitzen sollten. Die verwendeten Cyclodextrine wurden mit verschieden substituierten Benzoesäurederivaten umgesetzt und durch polarisationsmikroskopische Untersuchungen auf flüssigkristalline Phasenübergänge hin untersucht. Weitere Untersuchungen erfolgten mittels DSC und Röntgendiffraktometrie.In wie weit die flüssigkristallinen Eigenschaften eine Selbstorganisation zu Überstrukturen induzierten, sollten Untersuchungen als Monoschichten auf der Luft/Wasser-Grenzfläche mittels einer Filmwaage zeigen. Daraus ergaben sich weitergehende Untersuchungen mit der AFM-Technik. Es zeigte sich, dass besonders die quasi symmetrischen, über Esterverknüpfungen modifizierten Cyclodextrine flüssigkristalline Eigenschaften zeigten. Es konnte auch gezeigt werden, dass sie sich in Mono- und Multischichten zu Lipid-Strukturen organisieren. Der zweite Teil dieser Arbeit befasste sich mit der Wechselwirkung von Cyclodextrinen mit LCST-Copolymeren, die supramolekulare Komplexe bilden. Als Schlüsselgruppe fungierte das polymergebundene Adamantanmolekül. Mittels Trübungsphotometrie wurden verschiedene Additive auf ihren Einfluss auf den LCST-Phasenübergang der synthetisierten Copolymer untersucht. Hier konnte gezeigt werden, dass die adamantanhaltigen LCST-Copolymere durch Cyclodextrinzusätze in ihrem LCST-Verhalten drastisch zu beinflussen sind. Damit konnte ein System gefunden werden, das mit einfachen analytischen Mitteln eine supramolekulare Erkennung ermöglicht.

Beiträge zur Chemie der schwersten Elemente in den Bereichen Elektrochemie und Chromatographie

Im Rahmen dieser Arbeit wurden Experimente auf verschiedenen Gebieten der chemischen Charakterisierung der schwersten Elemente durchgeführt. So wurden Vorarbeiten zur Elektrochemie der superschweren Elemente geleistet. Hier konnte gezeigt werden, dass sich Po aus verdünnten Säuren auf Metallelektroden spontan abscheiden lässt. Als wichtigste Reaktionsparameter wurden Temperatur, Konvektion und Viskosität des Elektrolyten identifiziert. Für die Elektrodeposition ist es von Bedeutung, reproduzierbar saubere Elektroden einzusetzen. Auch hierzu konnten Erkenntnisse gewonnen werden. Auf dem Gebiet der Charakterisierung von Nukliden wurde der Elektroneneinfang in Db-263 untersucht. Aus diesem Zerfall geht Rf-263 hervor, welches mit einer Halbwertszeit von 15min überwiegend durch Spontanspaltung zerfällt. Im Experiment wurde Rf-263 mittels Ionenaustauschchromatographie und nachfolgender Flüssig-Flüssig-Extraktion von Db und Aktiniden getrennt. Das erhaltene Präparat wurde auf Alpha-Zerfall und Spontanspaltung untersucht. Die Zahl der Rf-263 Zerfälle lässt auf einen Zerfallszweig durch Elektroneneinfang in Db-263 von 3% schließen. Ein Großteil der Arbeit beschäftigt sich mit der Bestimmung des Kd-Wertes von Sg am Anionenaustauscher Aminex A27 in 0.1M HNO3/5E-3M HF. Beim mit dem ALOHA-System durchgeführten Experiment wurde die Mehrsäulentechnik (MCT) erstmals für die wässrige Chemie von Sg genutzt. Aufgrund zahlreicher Probleme konnte letztlich kein Kd-Wert bestimmt werden. Stattdessen wurde das Experiment einer Fehleranalyse unterzogen. Hier zeigte sich eine Anfälligkeit der MCT auf natürliche Radionuklide. Weiterhin konnten Probleme bei ALOHA aufgedeckt werden, die eine Wiederholung des Experiments ausschlossen. In der Folge wurde ein alternativer Anionenaustauscher charakterisiert, sowie die Elektrolysebedingungen, als wichtiger Schritt der chemischen Aufarbeitung bei der MCT, genauer spezifiziert.

Investigations on maar-diatreme volcanoes by inversion of magnetic and gravity data from the Eifel area, Germany

A study of maar-diatreme volcanoes has been perfomed by inversion of gravity and magnetic data. The geophysical inverse problem has been solved by means of the damped nonlinear least-squares method. To ensure stability and convergence of the solution of the inverse problem, a mathematical tool, consisting in data weighting and model scaling, has been worked out. Theoretical gravity and magnetic modeling of maar-diatreme volcanoes has been conducted in order to get information, which is used for a simple rough qualitative and/or quantitative interpretation. The information also serves as a priori information to design models for the inversion and/or to assist the interpretation of inversion results. The results of theoretical modeling have been used to roughly estimate the heights and the dip angles of the walls of eight Eifel maar-diatremes — each taken as a whole. Inversemodeling has been conducted for the Schönfeld Maar (magnetics) and the Hausten-Morswiesen Maar (gravity and magnetics). The geometrical parameters of these maars, as well as the density and magnetic properties of the rocks filling them, have been estimated. For a reliable interpretation of the inversion results, beside the knowledge from theoretical modeling, it was resorted to other tools such like field transformations and spectral analysis for complementary information. Geologic models, based on thesynthesis of the respective interpretation results, are presented for the two maars mentioned above. The results gave more insight into the genesis, physics and posteruptive development of the maar-diatreme volcanoes. A classification of the maar-diatreme volcanoes into three main types has been elaborated. Relatively high magnetic anomalies are indicative of scoria cones embeded within maar-diatremes if they are not caused by a strong remanent component of the magnetization. Smaller (weaker) secondary gravity and magnetic anomalies on the background of the main anomaly of a maar-diatreme — especially in the boundary areas — are indicative for subsidence processes, which probably occurred in the late sedimentation phase of the posteruptive development. Contrary to postulates referring to kimberlite pipes, there exists no generalized systematics between diameter and height nor between geophysical anomaly and the dimensions of the maar-diatreme volcanoes. Although both maar-diatreme volcanoes and kimberlite pipes are products of phreatomagmatism, they probably formed in different thermodynamic and hydrogeological environments. In the case of kimberlite pipes, large amounts of magma and groundwater, certainly supplied by deep and large reservoirs, interacted under high pressure and temperature conditions. This led to a long period phreatomagmatic process and hence to the formation of large structures. Concerning the maar-diatreme and tuff-ring-diatreme volcanoes, the phreatomagmatic process takes place due to an interaction between magma from small and shallow magma chambers (probably segregated magmas) and small amounts of near-surface groundwater under low pressure and temperature conditions. This leads to shorter time eruptions and consequently to structures of smaller size in comparison with kimberlite pipes. Nevertheless, the results show that the diameter to height ratio for 50% of the studied maar-diatremes is around 1, whereby the dip angle of the diatreme walls is similar to that of the kimberlite pipes and lies between 70 and 85°. Note that these numerical characteristics, especially the dip angle, hold for the maars the diatremes of which — estimated by modeling — have the shape of a truncated cone. This indicates that the diatreme can not be completely resolved by inversion.

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.

Computational modeling of surface interactions

Die Wechselwirkung zwischen Proteinen und anorganischen Oberflächen fasziniert sowohl aus angewandter als auch theoretischer Sicht. Sie ist ein wichtiger Aspekt in vielen Anwendungen, unter anderem in chirugischen Implantaten oder Biosensoren. Sie ist außerdem ein Beispiel für theoretische Fragestellungen betreffend die Grenzfläche zwischen harter und weicher Materie. Fest steht, dass Kenntnis der beteiligten Mechanismen erforderlich ist um die Wechselwirkung zwischen Proteinen und Oberflächen zu verstehen, vorherzusagen und zu optimieren. Aktuelle Fortschritte im experimentellen Forschungsbereich ermöglichen die Untersuchung der direkten Peptid-Metall-Bindung. Dadurch ist die Erforschung der theoretischen Grundlagen weiter ins Blickfeld aktueller Forschung gerückt. Eine Möglichkeit die Wechselwirkung zwischen Proteinen und anorganischen Oberflächen zu erforschen ist durch Computersimulationen. Obwohl Simulationen von Metalloberflächen oder Proteinen als Einzelsysteme schon länger verbreitet sind, bringt die Simulation einer Kombination beider Systeme neue Schwierigkeiten mit sich. Diese zu überwinden erfordert ein Mehrskalen-Verfahren: Während Proteine als biologische Systeme ausreichend mit klassischer Molekulardynamik beschrieben werden können, bedarf die Beschreibung delokalisierter Elektronen metallischer Systeme eine quantenmechanische Formulierung. Die wichtigste Voraussetzung eines Mehrskalen-Verfahrens ist eine Übereinstimmung der Simulationen auf den verschiedenen Skalen. In dieser Arbeit wird dies durch die Verknüpfung von Simulationen alternierender Skalen erreicht. Diese Arbeit beginnt mit der Untersuchung der Thermodynamik der Benzol-Hydratation mittels klassischer Molekulardynamik. Dann wird die Wechselwirkung zwischen Wasser und den [111]-Metalloberflächen von Gold und Nickel mittels eines Multiskalen-Verfahrens modelliert. In einem weiteren Schritt wird die Adsorbtion des Benzols an Metalloberflächen in wässriger Umgebung studiert. Abschließend wird die Modellierung erweitert und auch die Aminosäuren Alanin und Phenylalanin einbezogen. Dies eröffnet die Möglichkeit realistische Protein- Metall-Systeme in Computersimulationen zu betrachten und auf theoretischer Basis die Wechselwirkung zwischen Peptiden und Oberflächen für jede Art Peptide und Oberfläche vorauszusagen.

Spectroscopic signatures from first principles calculations: from surface adsorbates to liquids and polymers

In dieser Arbeit werden drei wasserstoffverbrückte Systeme in der kondensierten Phase mit Hilfe von first-principles-Elektronenstruktur-Rechnungen untersucht, die auf der Dichtefunktionaltheorie (DFT) unter periodischen Randbedingungen basieren. Ihre lokalen Konformationen und Wasserstoffbrückenbindungen werden mittels ab-initio Molekulardynamiksimulationen berechnet und weiterhin durch die Bestimmung ihrer spektroskopischen Parameter charakterisiert. Der Schwerpunkt liegt dabei auf lokalen Strukturen und auf schnellen Fluktuationen der Wasserstoffbrückenbindungen, welche von zentraler Bedeutung für die physikalischen und chemischen Eigenschaften der betrachteten Systeme sind. Die für die lokalen, instantanen Konformationen berechneten Spektren werden verwendet, um die physikalischen Prozesse, die hinter den untersuchten Phänomenen stehen, zu erklären: die Wasseradsorption auf metallischen Oberflächen, die Ionensolvatisierung in wässrigen Lösungen und der Protonentransport in protonleitenden Polymeren, welche Prototypen von Membranen für Brennstoffzellen sind. Die Möglichkeit der Vorhersage spektroskopischer Parameter eröffnet vielfältige Möglichkeiten des Dialogs zwischen Experimenten und numerischen Simulationen. Die in dieser Arbeit vorgestellten Ergebnisse zeigen, dass die Zuverlässigkeit dieser theoretischen Berechnungen inzwischen für viele experimentell relevante Systeme ein quantitatives Niveau erreicht hat.

Oxidierte Kohlenwasserstoffe in der oberen Troposphäre: Kinetische Untersuchungen zum Abbau und Einfluss auf HO x und Ozon-Chemie

Der erste Teil dieser Arbeit befasst sich mit der Kinetik der Reaktion des OH-Radikals mit Glykolaldehyd (HOCH2CHO). Die Geschwindigkeitskonstante k1 wurde für diese Reaktion temperaturabhängig bestimmt. Durch gepulste Photolyse wurden OH-Radikale erzeugt. Anschließend wurde die laserinduzierte Fluoreszenz der OH-Radikale bei 309 nm detektiert. Die ermittelte Geschwindigkeitskonstante k1 für die Reaktion von OH mit HOCH2CHO von (8,0 ± 0,8) x 10-12 cm3 Teilchen-1 s-1 erweist sich für den Temperaturbereich von 240 K < T < 362 K als temperaturunabhängig. Zwischen 60 und 250 Torr kann zudem keine Druckabhängigkeit für k1 beobachtet werden. Die unerwartet niedrigere Geschwindigkeitskonstante für die betrachtete Reaktion im Vergleich zur Reaktion von OH mit CH3CHO konnte anhand von Überlegungen zur Korrelation zwischen der C-H-Bindungsstärke und dem H-Abstraktionskanal erklärt werden. Im zweiten Teil dieser Arbeit wurde die Photochemie von Aceton (CH3C(O)CH3), Methylethylketon (C2H5C(O)CH3, MEK) und Acetylbromid (CH3C(O)Br) betrachtet. Für die Photolyse von Aceton (bei 248 nm und 266 nm), MEK (bei 248 nm) und Acetylbromid (bei 248 nm) wurden bei 298 ± 3 K druckabhängig zwischen 5 und 1600 Torr N2 Quantenausbeuten für die Methylbildung (Phi(CH3)) bestimmt. Nach gepulster Photolyse der betrachteten Moleküle wurden die transienten Absorptionssignale der Methylradikale bei 216,4 nm verfolgt. Die Quantenausbeuten wurden relativ zur Photolyse von Methyliodid (CH3I) unter gleichen Reaktionsbedingungen ermittelt. Die erhaltenen Quantenausbeuten für CH3-Radikale nehmen für die beiden Systeme Aceton / 248 nm (Phi(CH3, Aceton) = 1,42 – 0,99) und MEK / 248 nm (Phi(CH3, MEK) = 0,45 – 0,19) druckabhängig zu hohen Drücken ab. Die Druckabhängigkeit von Phi(CH3) wird auf die Konkurrenz zwischen Stoßrelaxation und Dissoziation der schwingungsangeregten Acetylradikale (CH3CO#) zurückgeführt. Für das System Aceton / 266 nm wird keine Druckabhängigkeit von Phi(CH3) = 0,93 ± 0,1 beobachtet. Dies wird damit erklärt, dass CH3CO# nicht genügend Energie besitzt, um die Barriere zur Dissoziation zu überschreiten. Bei der Photolyse von Acetylbromid bei 248 nm wird druckunabhängig Phi(CH3) = 0,92 ± 0,10 bestimmt. In diesem System dissoziieren die schwingungsangeregten Acetylradikale bei allen Drücken vollständig. Bei 266 nm wurde die Gesamtquantenausbeute für die Photodissoziation von Aceton (Phi(diss, 266nm)) bestimmt. Die nach Photolyse erhaltenen Methyl - und Acetylradikale wurden nach Titration mit Br2 durch die Resonanzfluoreszenz der Bromatome detektiert. Phi(diss, 266nm) wurde mit 0,92 ± 0,07 bestimmt.

First principle study of ions in water and ionic liquids

Die beiden in dieser Arbeit betrachteten Systeme, wässrige Lösungen von Ionen und ionische Flüssigkeiten, zeigen vielfältige Eigenschaften und Anwendungsmöglichkeiten, im Gegensatz zu anderen Systemen. Man findet sie beinahe überall im normalen Leben (Wasser), oder ihre Bedeutung wächst (ioinische Flüssigkeiten). Der elektronische Anteil und der atomare Anteil wurden getrennt voneinander untersucht und im Zusammenhang analysiert. Mittels dieser Methode konnten die in dem jeweiligen System auftretenden Mechanismen genauer untersucht werden. Diese Methode wird "Multiscale Modeling" genannt, dabei werden die Untereinheiten eines Systems genauer betrachtet, wie in diesem Fall die elektronischen and atomaren Teilsystem. Die Ergebnisse, die aus den jeweiligen Betrachtungen hervorgehen, zeigen, dass, im Falle von hydratisierten Ionen die Wasser-Wasser Wechselwirkungen wesentlich stärker sind als die elektrostatischen Wechselwirkung zwischen Wasser und dem Ion. Anhand der Ergebnisse ergibt sich, dass normale nicht-polarisierbare Modelle ausreichen, um Ionen-Wasser Lösungen zu beschreiben. Im Falle der ionischen Flüssigkeiten betrachten wir die elektronische Ebene mittels sehr genauer post-Hartree-Fock Methoden und DFT, deren Ergebnisse dann mit denen auf molekularer Ebene (mithilfe von CPMD/klassischer MD) in Beziehung gesetzt werden. Die bisherigen Ergebnisse zeigen, dass die Wasserstoff-Brückenbindungen im Fall der ionischen Flüssigkeiten nicht vernachässigt werden können. Weiterhin hat diese Studie herausgefunden, dass die klassischen Kraftfelder die Elektrostatik (Dipol- und Quadrupolmomente) nicht genau genug beschreibt. Die Kombination des mikroskopischen Mechanismus und der molekularen Eigenschaften ist besonders sinnvoll um verschiedene Anhaltspunkte von Simualtionen (z.B. mit klassische Molekular-Dynamik) oder Experimenten zu liefern oder solche zu erklären.

Chemie an funktionalen Polymer-Opalen

In dieser Arbeit wird die Herstellung und Anwendung von funktionalen Polymer-Opalen beschrieben. Für die Synthese von funktionalen monodipsersen Kolloiden, den Bausteinen der Opale, wird die emulgatorfreie Emulsionspolymerisation (SFEP) verwendet. Je nach einzubauendem funktionalem Molekül werden verschiedene Varianten der SFEP verwendet, wie z. B. Homopolymerisation, Copolymerisation, Polymerisation mit Fremdstoffen und die Herstellung von Kern-Schale-Kolloiden. Die so hergestellten monodispersen Kolloide formen durch Selbstorganisation über horizontale (Aufpipettieren, Rakeln, Sprühen) oder vertikale Kristallisation (Ziehmaschine)hochqualitative künstliche Opale. Die eingebauten Funktionalitäten öffnen den Weg zu einer Vielzahl von Anwendungen. Über die Spaltung von funktionalen Estergruppen kann eine lichtinduzierte Strukturierung durchgeführt werden. Der Einbau von Epoxidgruppen ermöglicht eine makroskopische Vernetzung wodurch die mechanische Stabilität der Struktur erhöht wird. Der Einsatz von Reaktivestern kann zur Oberflächen- funktionalisierung verwendet werden. Durch Replizierung der Struktur zum inversen Opal können weitere funktionale Materialien eingeführt werden, was die Einsatzmöglichkeiten noch erweitert.