Simulations of physics and chemistry of polar stratospheric clouds with a general circulation model

A polar stratospheric cloud submodel has been developed and incorporated in a general circulation model including atmospheric chemistry (ECHAM5/MESSy). The formation and sedimentation of polar stratospheric cloud (PSC) particles can thus be simulated as well as heterogeneous chemical reactions that take place on the PSC particles. For solid PSC particle sedimentation, the need for a tailor-made algorithm has been elucidated. A sedimentation scheme based on first order approximations of vertical mixing ratio profiles has been developed. It produces relatively little numerical diffusion and can deal well with divergent or convergent sedimentation velocity fields. For the determination of solid PSC particle sizes, an efficient algorithm has been adapted. It assumes a monodisperse radii distribution and thermodynamic equilibrium between the gas phase and the solid particle phase. This scheme, though relatively simple, is shown to produce particle number densities and radii within the observed range. The combined effects of the representations of sedimentation and solid PSC particles on vertical H2O and HNO3 redistribution are investigated in a series of tests. The formation of solid PSC particles, especially of those consisting of nitric acid trihydrate, has been discussed extensively in recent years. Three particle formation schemes in accordance with the most widely used approaches have been identified and implemented. For the evaluation of PSC occurrence a new data set with unprecedented spatial and temporal coverage was available. A quantitative method for the comparison of simulation results and observations is developed and applied. It reveals that the relative PSC sighting frequency can be reproduced well with the PSC submodel whereas the detailed modelling of PSC events is beyond the scope of coarse global scale models. In addition to the development and evaluation of new PSC submodel components, parts of existing simulation programs have been improved, e.g. a method for the assimilation of meteorological analysis data in the general circulation model, the liquid PSC particle composition scheme, and the calculation of heterogeneous reaction rate coefficients. The interplay of these model components is demonstrated in a simulation of stratospheric chemistry with the coupled general circulation model. Tests against recent satellite data show that the model successfully reproduces the Antarctic ozone hole.

Non-equilibrium molecular dynamics study of an amphiphilic model system

Diese Doktorarbeit untersucht das Verhalten von komplexenFluidenunter Scherung, insbesondere den Einfluss von Scherflüssenauf dieStrukturbildung.Dazu wird ein Modell dieser entworfen, welches imRahmen von Molekulardynamiksimulationen verwendet wird.Zunächst werden Gleichgewichtseigenschaften dieses Modellsuntersucht.Hierbei wird unter anderem die Lage desOrdnungs--Unordnungsübergangs von derisotropen zur lamellaren Phase der Dimere bestimmt.Der Einfluss von Scherflüssen auf diese lamellare Phase wirdnununtersucht und mit analytischen Theorien verglichen. Die Scherung einer parallelen lamellaren Phase ruft eineNeuausrichtung des Direktors in Flussrichtung hervor.Das verursacht eine Verminderung der Schichtdicke mitsteigender Scherrateund führt oberhalb eines Schwellwertes zu Ondulationen.Ein vergleichbares Verhalten wird auch in lamellarenSystemengefunden, an denen in Richtung des Direktors gezogen wird.Allerdings wird festgestellt, dass die Art der Bifurkationenin beidenFällen unterschiedlich ist.Unter Scherung wird ein Übergang von Lamellen parallelerAusrichtung zu senkrechter gefunden.Dabei wird beoachtet, dass die Scherspannung in senkrechterOrientierungniedriger als in der parallelen ist.Dies führt unter bestimmten Bedingungen zum Auftreten vonScherbändern, was auch in Simulationen beobachtet wird. Es ist gelungen mit einem einfachen Modell viele Apsekte desVerhalten vonkomplexen Fluiden wiederzugeben. Die Strukturbildung hängt offensichtlich nurbedingt von lokalen Eigenschaften der Moleküle ab.

Thermodynamic concepts in adaptive resolution simulations

This thesis work is devoted to the conceptual and technical development of the Adaptive Resolution Scheme (AdResS), a molecular dynamics method that allows the simulation of a system with different levels of resolution simultaneously. The simulation domain is divided into high and low resolution zones and a transition region that links them, through which molecules can freely diffuse.rnThe first issue of this work regards the thermodynamic consistency of the method, which is tested and verified in a model liquid of tetrahedral molecules. The results allow the introduction of the concept of the Thermodynamic Force, an external field able to correct spurious density fluctuations present in the transition region in usual AdResS simulations.rnThe AdResS is also applied to a system where two different representations with the same degree of resolution are confronted. This simple test extends the method from an Adaptive Resolution Scheme to an Adaptive Representation Scheme, providing a way of coupling different force fields based on thermodynamic consistency arguments. The Thermodynamic Force is successfully applied to the example described in this work as well.rnAn alternative approach of deducing the Thermodynamic Force from pressure consistency considerations allows the interpretation of AdResS as a first step towards a molecular dynamics simulation in the Grand Canonical ensemble. Additionally, such a definition leads to a practical way of determining the Thermodynamic Force, tested in the well studied tetrahedral liquid. The effects of AdResS and this correction on the atomistic domain are analyzed by inspecting the local distribution of velocities, radial distribution functions, pressure and particle number fluctuation. Their comparison with analogous results coming from purely atomistic simulations shows good agreement, which is greatly improved under the effect of the external field.rnA further step in the development of AdResS, necessary for several applications in biophysics and material science, consists of its application to multicomponent systems. To this aim, the high-resolution representation of a model binary mixture is confronted with its coarse-grained representation systematically parametrized. The Thermodynamic Force, whose development requires a more delicate treatment, also gives satisfactory results.rnFinally, AdResS is tested in systems including two-body bonded forces, through the simulation of a model polymer allowed to adaptively change its representation. It is shown that the distribution functions that characterize the polymer structure are in practice not affected by the change of resolution.rnThe technical details of the implementation of AdResS in the ESPResSo package conclude this thesis work.

Experimentelle Untersuchung von Phasenübergängen im kolloidalen Modellsystem harter Kugeln

Die vorliegende Arbeit behandelt den fluid-kristallinen Phasenübergang sowie den Glasübergang anhand von kolloidalen Hart-Kugel(HK)-Modellsystemen. Die Untersuchungen erfolgen dabei im Wesentlichen mit unterschiedlichen Lichtstreumethoden und daher im reziproken Raum. rnDie Analyse der Kristallisationskinetik zeigt, dass es bei der Kristallisation zu signifikanten Abweichungen vom Bild der klassischen Nukleationstheorie (CNT) kommt. Diese geht von einem einstufigen Nukleationsprozess aus, wohingegen bei den hier durchgeführten Experimenten ein mehrstufiger Prozess beobachtet wird. Vor der eigentlichen Kristallisation kommt es zunächst zur Nukleation einer metastabilen Zwischenphase, sogenannter Precursor. In einer zweiten Stufe erfolgt innerhalb der Precursor die eigentliche Nukleation der Kristallite. rnDurch weitere Analyse und den Vergleich des Kristallisations- und Verglasungsszenarios konnte das Konzept der Precursornukleation auf den Vorgang der Verglasung erweitert werden. Während die Kristallnukleation oberhalb des Glasübergangspunktes zum Erliegen kommt, bleibt der Prozess der Precursornukleation auch bei verglasenden Proben erhalten. Ein Glas erstarrt somit in einem amorphen Zustand mit lokalen Precursorstrukturen. Die Korrelation der gemessenen zeitlichen Entwicklung der strukturellen sowie der dynamischen Eigenschaften zeigt darüber hinaus, dass das bisher unverstandene Ageing-Phänomen von HK-Gläsern mit der Nukleation von Precursorn zusammenhängt.rnEin solches mehrstufiges Szenario wurde bereits in früheren Veröffentlichungen beobachtet. Die im Rahmen dieser Arbeit durchgeführten Messungen ermöglichten erstmals die Bestimmung von Kristallnukleationsratendichten (Kristall-NRD) und Ratendichten für die Precursornukleation bis über den Glasübergangspunkt hinaus. Die Kristall-NRD bestätigen die Resultate aus anderen experimentellen Arbeiten. Die weiteren Analysen der Kristall-NRD belegen, dass die fluid-kristalline Grenzflächenspannung bei der Nukleation entgegen den Annahmen der CNT nicht konstant ist, sondern mit ansteigendem Volumenbruch linear zunimmt. Die Erweiterung der CNT um eine linear zunehmende Grenzflächenspannung ermöglichte eine quantitative Beschreibung der gemessenen Kristall- sowie der Precursor-NRD, was den Schluss zulässt, dass es sich in beiden Fällen um einen Boltzmann-aktivierten Prozess handelt. rnUm die beobachteten Abweichungen des Nukleationsprozesses vom Bild der CNT näher zu untersuchen, wurden die kollektiven Partikeldynamiken in stabilen Fluiden und metastabilen Schmelzen analysiert. Im klassischen Bild wird angenommen, dass die kollektive Partikeldynamik beim Vorgang der Nukleation keine Rolle spielt. Anhand der Resultate zeigen sich Abweichungen in der Dynamik stabiler Fluide und metastabiler Schmelzen. Während die kollektive Partikeldynamik in der stabilen Schmelze von der Struktur entkoppelt ist, tritt oberhalb des Phasenübergangspunktes eine Kopplung von Struktur und Dynamik auf. Dabei treten die Abweichungen zunächst in der Umgebung des ersten Strukturfaktormaximums und somit bei den am stärksten besetzten Moden auf. Mit steigender Unterkühlung nehmen die Anzahl der abweichenden Moden sowie die Stärke der Abweichungen zu. Dieses Phänomen könnte einen signifikanten Einfluss auf den Nukleationsprozess und somit auf die Kristallisationskinetik haben. Die Analyse der Dynamik im stabilen Fluid zeigt darüber hinaus Hinweise auf eine Singularität bei Annäherung an den fluid-kristallinen Phasenübergangspunkt.rnDes Weiteren wurden im Rahmen der vorliegenden Arbeit erstmals Ratendichten für die heterogene Nukleation eines HK-Systems an einer flachen Wand mittels statischer Lichtstreuung (SLS) bestimmt. Die Ergebnisse der Messung zeigen, dass die Nukleationsbarriere der heterogenen Nukleation annähernd Null ist und folglich eine vollständige Benetzung der Wand mit einer kristallinen Monolage vorliegt. Die Erweiterung der Untersuchungen auf gekrümmte Oberflächen in Form von sphärischen Partikeln (Seeds) stellt die erste experimentelle Arbeit dar, die den Einfluss eines Ensembles von Seeds auf die Kristallisationskinetik in HK-Systemen untersucht. Die Kristallisationskinetik und die Mikrostruktur werden abhängig von Größe und Anzahldichte der Seed-Partikel signifikant beeinflusst. In Übereinstimmung mit konfokalmikroskopischen Experimenten und Simulationen spielt dabei das Radienverhältnis der Majoritäts- zur Minoritätskomponente eine entscheidende Rolle.