Charge transport in amorphous organic semiconductors

Organic semiconductors with the unique combination of electronic and mechanical properties may offer cost-effective ways of realizing many electronic applications, e.g. large-area flexible displays, printed integrated circuits and plastic solar cells. In order to facilitate the rational compound design of organic semiconductors, it is essential to understand relevant physical properties e.g. charge transport. This, however, is not straightforward, since physical models operating on different time and length scales need to be combined. First, the material morphology has to be known at an atomistic scale. For this atomistic molecular dynamics simulations can be employed, provided that an atomistic force field is available. Otherwise it has to be developed based on the existing force fields and first principle calculations. However, atomistic simulations are typically limited to the nanometer length- and nanosecond time-scales. To overcome these limitations, systematic coarse-graining techniques can be used. In the first part of this thesis, it is demonstrated how a force field can be parameterized for a typical organic molecule. Then different coarse-graining approaches are introduced together with the analysis of their advantages and problems. When atomistic morphology is available, charge transport can be studied by combining the high-temperature Marcus theory with kinetic Monte Carlo simulations. The approach is applied to the hole transport in amorphous films of tris(8-hydroxyquinoline)aluminium (Alq3). First the influence of the force field parameters and the corresponding morphological changes on charge transport is studied. It is shown that the energetic disorder plays an important role for amorphous Alq3, defining charge carrier dynamics. Its spatial correlations govern the Poole-Frenkel behavior of the charge carrier mobility. It is found that hole transport is dispersive for system sizes accessible to simulations, meaning that calculated mobilities depend strongly on the system size. A method for extrapolating calculated mobilities to the infinite system size is proposed, allowing direct comparison of simulation results and time-of-flight experiments. The extracted value of the nondispersive hole mobility and its electric field dependence for amorphous Alq3 agree well with the experimental results.

Ionic liquids in a multiscale study

Trotz des hohen Interesse an Ionischen Flüssigkeiten wird das zielgerichtete Design und die Anwendung Ionischer Flüssigkeiten durch fehlendes grundlegendes Verständnis erschwert. Deshalb wurde die Balance der molekularen Wechselwirkungen in Ionischen Flüssigkeiten studiert, um die Eigenschaften dieser zu verstehen und die Kraftfeldentwicklung im Rahmen des Multiskalenansatzes zu systematisieren. Es wurden reine Imidazolium-basierte Ionische Flüssigkeiten, Mischungen mit kleinen Molekülen und eine protische Ionische Flüssigkeit mit ab-initio-Methoden, hauptsächlich Car-Parrinello-Molekulardynamik, untersucht. Weiterhin wurden Eigenschaften der Flüssigphase mit denen von Ionenpaaren verglichen.rnIm Fokus standen die molekularen elektrostatischen Eigenschaften und es wurde gezeigt, dass Coulomb-Wechselwirkungen zu einzigartigen Charakteristika führten. So waren die Ionen-Nettoladungen stets reduziert, die molekularen Dipolmomentverteilungen sehr breit, elektronische Polarisation war entscheidend. Die elektrostatischen Eigenschaften waren allgemein lokal auf molekularen Größen- und Zeitskalen und hingen stark von Phasenzustand und Zusammensetzung ab. Für andere molekulare Eigenschaften, wie der Neigung zu dispersiven Kontakten oder Wasserstoffbrücken, wurde gezeigt, dass sie einen entscheidenden Einfluss auf die Feinstruktur Ionischer Flüssigkeiten hatten. Das Gleichgewicht der Wechselwirkungen zeigte sich auch in Leistungsspektren, die sich aus den ab-initio-Molekulardynamiksimulationen ergaben. Diese boten einen neuen Weg für den Vergleich zum Experiment und für einen Einblick in die schnelle Dynamik Ionischer Flüssigkeiten.

Calibration study of bivalve shells : implications for paleoenvironmental reconstruction

The accretionary shells of bivalve mollusks can provide environmental information, such as water temperature, precipitation, freshwater fluxes, primary productivity and anthropogenic activities in the form of variable growth rates and variable geochemical properties, such as stable oxygen and carbon isotopes. However, paleoenvironmental reconstructions are constrained by uncertainties about isotopic equilibrium fractionation during shell formation, which is generally acknowledged as a reasonable assumption for bivalves, but it has been disputed in several species. Furthermore, the variation in shell growth rates is accepted to rely on multiple environmental variables, such as temperature, food availability and salinity, but can differ from species to species. Therefore, it is necessary to perform species-specific calibration studies for both isotope proxies and shell growth rates before they can be used with confidence for environmental interpretations of the past. Accordingly, the principal objective of this Ph.D research is to examine the reliability of selected bivalve species, the long-lived Eurhomalea exalbida (Dillwyn), the short-lived and fast growing species Paphia undulata (Born 1778), and the freshwater mussel Margaritifera falcata (Gould 1850), as paleoenvironmental proxy archives.rnThe first part is focused on δ18Oshell and shell growth history of live-collected E. exalbida from the Falkland Islands. The most remarkable finding, however, is that E. exalbida formed its shell with an offset of -0.48‰ to -1.91‰ from the expected oxygen isotopic equilibrium with the ambient water. If this remained unnoticed, paleotemperature estimates would overestimate actual water temperatures by 2.1-8.3°C. With increasing ontogenetic age, the discrepancy between measured and reconstructed temperatures increased exponentially, irrespective of the seasonally varying shell growth rates. This study clearly demonstrates that, when the disequilibrium fractionation effect is taken into account, E. exalbida can serve as a high-resolution paleoclimate archive for the southern South America. The species therefore provides quantifiable temperature estimates, which yields new insights into long-term paleoclimate dynamics for mid to high latitudes on the southern hemisphere.rnThe stable carbon isotope of biogenic carbonates is generally considered to be useful for reconstruction of seawater dissolved inorganic carbon. The δ13Cshell composition of E. exalbida was therefore, investigated in the second part of this study. This chapter focuses on inter-annual and intra-annual variations in δ13Cshell. Environmental records in δ13Cshell are found to be strongly obscured by changes in shell growth rates, even if removing the ontogenetic decreasing trend. This suggests that δ13Cshell in E. exalbida may not be useful as an environmental proxy, but a potential tool for ecological investigations. rnIn addition to long-lived bivalve species, short-lived species that secrete their shells extremely fast, can also be useful for environmental reconstructions, especially as a high-resolution recorder. Therefore, P. undulata from Daya Bay, South China Sea was utilized in Chapter 4 to evaluate and establish a potential proxy archive for past variations of the East Asian monsoon on shorter time-scales. The δ18Oshell can provide qualitative estimates of the amount of monsoonal rain and terrestrial runoff and the δ13Cshell likely reflect the relative amount of isotopically light terrestrial carbon that reaches the ocean during the summer monsoon season. Therefore, shells of P. undulata can provide serviceable proxy archives to reconstruct the frequency of exceptional summer monsoons in the past. The relative strength of monsoon-related precipitation and associated changes in ocean salinity and the δ13C ratios of the dissolved inorganic carbon signature (δ13CDIC) can be estimated from the δ18Oshell and δ13Cshell values as well as shell growth patterns. rnIn the final part, the freshwater pearl shell M. falcata from four rivers in British Columbia, Canada was preliminarily studied concerning the lifespans and the shell growth rates. Two groups separated by the Georgia Strait can be clearly distinguished. Specimens from the western group exhibit a shorter lifespan, while the eastern group live longer. Moreover, the average lifespan seems to decrease from south to north. The computed growth equations from the eastern and western groups differ as well. The western group exhibits a lower growth rate, while bivalves from the eastern group grow faster. The land use history seems to be responsible for the differences in lifespans of the specimens from the two groups. Differences in growth rate may be induced by differences in water temperature or nutrient input also related to the land use activities.

Structure and dynamics of ionic liquid,electrode interfaces

In the early 20th century, Gouy, Chapman, and Stern developed a theory to describe the capacitance and the spatial ion distribution of diluted electrolytes near an electrode. After a century of research, considerable progress has been made in the understanding of the electrolyte/electrode interface. However, its molecular-scale structure and its variation with an applied potential is still under debate. In particular for room-temperature ionic liquids, a new class of solventless electrolytes, the classical theories for the electrical double layer are not applicable. Recently, molecular dynamics simulations and phenomenological theories have attempted to explain the capacitance of the ionic liquid/electrode interface with the molecular-scale structure and dynamics of the ionic liquid near the electrode. rnHowever, experimental evidence is very limited. rnrnIn the presented study, the ion distribution of an ionic liquid near an electrode and its response to applied potentials was examined with sub-molecular resolution. For this purpose, a new sample chamber was constructed, allowing in situ high energy X-ray reflectivity experiments under potential control, as well as impedance spectroscopy measurements. The combination of structural information and electrochmical data provided a comprehensive picture of the electric double layer in ionic liquids. Oscillatory charge density profiles were found, consisting of alternating anion- and cation-enriched layers at both, cathodic and anodic, potentials. This structure was shown to arise from the same ion-ion correlations dominating the liquid bulk structure that were observed as a distinct X-ray diffraction peak. Therefore, existing physically motivated models were refined and verified by comparison with independent measurements. rnrnThe relaxation dynamics of the interfacial structure upon potential variation were studied by time resolved X-ray reflectivity experiments with sub-millisecond resolution. The observed relaxation times during charging/discharging are consistent with the impedance spectroscopy data revealing three processes of vastly different characteristic time-scales. Initially, the ion transport normal to the interface happens on a millisecond-scale. Another 100-millisecond-scale process is associated with molecular reorientation of electrode-adsorbed cations. Further, a minute-scale relaxation was observed, which is tentatively assigned to lateral ordering within the first layer.

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.

Chemical interaction between ocean and atmosphere

The exchange of chemical constituents between ocean and atmosphere provides potentially important feedback mechanisms in the climate system. The aim of this study is to develop and evaluate a chemically coupled global atmosphere-ocean model. For this, an atmosphere-ocean general circulation model with atmospheric chemistry has been expanded to include oceanic biogeochemistry and the process of air-sea gas exchange. The calculation of seawater concentrations in the oceanic biogeochemistry submodel has been expanded from DMS, CO₂

Austauschprozesse an Tropopausenfalten extratropischer Zyklonen

Troposphärisches Ozon ist bekannt als wichtiges Oxidationsmittel und als Vorläufergas hoch reaktiver Radikale. Es zählt zu den wichtigsten Treibhausgasen und wirkt bei hohen Konzentrationen an der Erdoberfläche giftig für alle Lebewesen. Zwar wird der Großteil des troposphärischen Ozons photochemisch produziert, ein erheblicher Anteil hat aber stratosphärischen Ursprung und wird entlang von Tropopausenfalten in Zyklonen in die Troposphäre transportiert. Dieser Transport von Luftmassen aus der Stratosphäre in diernTroposphäre (STT) kann zu einem kurzzeitigen, starken Ozonanstieg am Boden führen und langfristig die Chemie der Troposphäre beeinflussen. Die Quantifizierung des Ozoneintrages und die Identifizierung der dafür verantwortlichen Prozesse ist mit großen Unsicherheiten belastet und ein aktuelles Forschungsthema.rnAufgrund ihrer groben Auflösung ist es mit globalen Modellen nicht möglich, die Details dieser STT-Prozesse zu erfassen. Deshalb wird in dieser Arbeit das Modellsystem MECO(n) genutzt, welches das regionale Atmosphärenchemie- und Klimamodell COSMO/MESSy mit dem globalen Chemie-Klimamodell ECHAM5/MESSy (EMAC) koppelt. Eine einheitliche Prozessparametrisierung ermöglicht konsistente, simultane Simulationen in verschiedenen Modellauflösungen. Ein als Teil dieser Arbeit neu entwickeltes Submodell erlaubt die Initialisierung künstlicher, passiver Tracer in Abhängigkeit verschiedener Variablen. Mit einem auf diese Weise freigesetzten, stratosphärischen Tracer lässt sich Ozon mit stratosphärischer Herkunft von solchem, das photochemisch produziert wurde, unterscheiden.rnIm Rahmen einer Fallstudie werden die Austauschprozesse an einer Tropopausenfalte sowohl aus der Eulerischen, als auch aus der Lagrangeschen Perspektive betrachtet. Die Analyse der STT-Prozesse zeigt, dass Luftmassen aus der Stratosphäre durch turbulente und diabatische Prozesse am Rand der Tropopausenfalte in die Troposphäre gelangen und anschließend bis zum Boden transportiert werden. Diese absinkenden, stratosphärischen Luftmassen führen in den Simulationen zu Ozonanstiegen am Boden, die mit Beobachtungsdaten evaluiert werden können. Es wird gezeigt, dass die Ergebnisse der feiner auflösendenrnModellinstanz gut mit den Messungen übereinstimmen.rnIn einer Lagrangeschen Analyse lassen sich Mischungszeitskalen für STT-Prozesse bestimmen. Es wird gezeigt, dass Luftpakete, die sich länger als zehn Stunden in der Troposphäre aufhalten, diese durch den Eintrag ihrer stratosphärischen Tracereigenschaften beeinflussen und daher nicht vernachlässigbar sind. Eine weitere Studie gibt Aufschluss über die Effektivität der Mischung an Tropopausenfalten: Fast die gesamte Luftmasse, die sich zu einem bestimmten Zeitpunkt in der Tropopausenfalte befindet, gelangt innerhalb von zwei Tagen in die Troposphäre.

Untersuchung von Mischungs- und Transportprozessen in der oberen Troposphäre - unteren Stratosphäre basierend auf in-situ Spurengasmessungen

Die obere Troposphäre / untere Stratosphäre (UTLS: Upper Troposphere / Lower Stratosphere)ist die Übergangsgregion zwischen den dynamisch, chemisch und mikrophysikalisch sehr verschiedenen untersten Atmosphärenschichten, der Troposphäre und der Stratosphäre. Strahlungsaktive Spurengase, wie zum Beispiel Wasserdampf (H2O), Ozon (O3) oder Kohlenstoffdioxid (CO2), und Wolken in der UTLS beeinflussen das Strahlungsbudget der Atmosphäre und das globale Klima. Mögliche Veränderungen in den Verteilungen und Konzentrationen dieser Spurengase modifizieren den Strahlungsantrieb der Atmosphäre und können zum beobachteten Klimawandel beitragen. Ziel dieser Arbeit ist es, Austausch- und Mischungsprozesse innerhalb der UTLS besser zu verstehen und damit Veränderungen der Spurengaszusammensetzung dieser Region genauer prognostizieren zu können. Grundlage hierfür bilden flugzeuggetragene in-situ Spurengasmessungen in der UTLS, welche während der Flugzeugmesskampagnen TACTS / ESMVal 2012 und AIRTOSS - ICE 2013 durchgeführt wurden. Hierbei wurde bei den Messungen von AIRTOSS - ICE 2013 das im Rahmen dieser Arbeit aufgebaute UMAQS (University of Mainz Airborne QCLbased Spectrometer) - Instrument zur Messung der troposphärischen Spurengase Distickstoffmonoxid (N2O) und Kohlenstoffmonoxid (CO) eingesetzt. Dieses erreicht bei einer zeitlichen Auflösung von 1 s eine Messunsicherheit von 0,39 ppbv und 1,39 ppbv der N2O bzw. CO-Mischungsverhältnisse. Die hohe Zeitauflösung und Messgenauigkeit der N2O- und CO- Daten erlaubt die Untersuchung von kleinskaligen Austauschprozessen zwischen Troposphäre und Stratosphäre im Bereich der Tropopause auf räumlichen Skalen kleiner 200 m. Anhand der N2O-Daten von AIRTOSS - ICE 2013 können in-situ detektierte Zirruspartikel in eisübersättigter Luft oberhalb der N2O-basierten chemischen Tropopause nachgewiesen werden. Mit Hilfe der N2O-CO-Korrelation sowie der Analyse von ECMWF-Modelldaten und der Berechnung von Rückwärtstrajektorien kann deren Existenz auf das irreversible Vermischen von troposphärischen und stratosphärischen Luftmassen zurückgeführt werden. Mit den in-situ Messungen von N2O, CO und CH4 (Methan) von TACTS und ESMVal 2012 werden die großräumigen Spurengasverteilungen bis zu einer potentiellen Temperatur von Theta = 410 K in der extratropischen Stratosphäre untersucht. Hierbei kann eine Verjüngung der Luftmassen in der extratropischen Stratosphäre mit Delta Theta > 30 K (relativ zur dynamischen Tropopause) über den Zeitraum der Messkampagne (28.08.2012 - 27.09.2012) nachgewiesen werden. Die Korrelation von N2O mit O3 zeigt, dass diese Verjüngung aufgrund des verstärkten Eintrages von Luftmassen aus der tropischen unteren Stratosphäre verursacht wird. Diese werden über den flachen Zweig der Brewer-Dobson-Zirkulation auf Zeitskalen von wenigen Wochen in die extratropische Stratosphäre transportiert. Anhandrnder Analyse der CO-O3-Korrelation eines Messfluges vom 30.08.2012 wird das irreversible Einmischen von Luftmassen aus der tropischen Stratosphäre in die Extratropen auf Isentropen mit Theta > 380 K identifiziert. Rückwärtstrajektorien zeigen, dass der Ursprung der eingemischten tropischen Luftmassen im Bereich der sommerlichen Antizyklone des asiatischen Monsuns liegt.

Optische Experimente zur Untersuchung der Nichtgleichgewichtsdynamik kolloidaler Suspensionen

In der vorliegenden Forschungsarbeit wird die Konkurrenz von Kristallisation und Verglasung in metastabilen kolloidalen Hart-Kugel(HK)-Modellsystemen mit dynamischer und zeitaufgelöster statischer Lichtstreuung untersucht. Durch gleichzeitige Messungen mit beiden Methoden an derselben Probe gelang es, mit hoher Genauigkeit und aussagekräftiger Statistik nachzuweisen, dass in beiden Systemklassen eine starke Korrelation der strukturellen und dynamischen Eigenschaften vorliegt und diese Korrelation zu quantifizieren. Ein zentraler Teil der Arbeit bestand in dem Aufbau einer geeigneten Lichtstreuanlage mit der erstmalig Messungen der Dynamik und der Struktur simultan an derselben nicht-ergodischen Probe durchgeführt werden konnten. Für die dynamische Lichtstreuung wurde ein Flächendetektor (CCD-Kamera) verwendet. In Kombination mit einer speziellen Detektionsoptik ermöglicht dies, die gleichzeitige Detektion von Streulicht aus unterschiedlichen Probenbereichen (Subensembles). Damit kann gleichzeitig die Dynamik in unterschiedlichen Subensembles mit einer Auflösung von 15,8x15,8µm2 untersucht werden. Die Lichtstreuanlage wurde ausführlich charakterisiert und ihre korrekte Funktionsweise mithilfe von Vergleichsmessungen an etablierten Lichtstreuanlagen bestätigt. Die zeitliche Entwicklung der Dynamik und der Struktur von metastabilen Proben wurde unterhalb, am und oberhalb des Glasübergangs quantifiziert. Dabei zeigte das untersuchte kolloidale HK-Modellsystem alle typischen Eigenschaften eines HK-Systems. Die kristallisierenden Proben zeigten das etablierte zweistufige Kristallisationsszenario mit entsprechender Kristallisationskinetik und die Glasproben zeigten das erwartete Alterungsverhalten. Bei dem zweistufigen Kristallisationsszenario kommt es zuerst zur Nukleation einer metastabilen Zwischenphase von sogenannten Precursorn. In einer zweiten Stufe bilden sich Kristallite innerhalb dieser Precursor. Durch Vergleich zwischen kristallisierenden und verglasenden Proben konnte auch während der Verglasung die Bildung von Precursorn beobachtet werden. Die Korrelation zwischen der Anzahl an Precursorn und der Partikeldynamik legt die Vermutung nahe, dass das immer noch unverstandene Phänomen der Alterung von Gläsern mit der Bildung von Precursorn zusammenhängt. Verhinderte Kristallisation führt zu einer starken Verlangsamung der Partikeldynamik. Die Partikeldynamik einer Probe am Glasübergang zeigt, dass die Probe vor Einsetzen der Kristallisation eine glasartige Dynamik aufwies. Dies legt die Vermutung nahe, dass einkomponentige kolloidale HK-Gläser den Gleichgewichtszustand (Kristall) auf langen Zeitskalen erreichen können. Durch die Untersuchung der Partikeldynamik von metastabilen Proben in einzelnen Subensembles konnte eine heterogene Verteilung der Partikeldynamik nachgewiesen werden. Es existieren Bereiche, in denen die Partikeldynamik schneller oder langsamer ist als in anderen Bereichen. Gleichzeitig zeigen die Messungen der strukturellen Eigenschaften, dass metastabile Proben auch heterogen bezüglich ihrer Struktur sind. Mithilfe dieser Messungen konnte die zeitliche Entwicklung des Anteils an langsamen Partikeln und des Anteils an Partikeln innerhalb von Objekten höherer Ordnung (Precursor/Kristallite) bestimmt werden. Es zeigte sich eine direkte Korrelation zwischen dem Anteil an langsamen Partikeln und dem Anteil an Partikeln in Objekten höherer Ordnung. Die Untersuchung der Dynamik und der Struktur in einzelnen Subensembles lieferte einen weiteren Hinweis darauf, dass Subensembles, in denen eine stärker ausgeprägte strukturelle Ordnung vorliegt, auch bevorzugt eine langsamere Partikeldynamik aufweisen.