Supramolecular order and dynamics of discotic materials studied by solid-state NMR recoupling methods

The topic of this thesis is the investigation of structure,order and dynamics in discotic mesogens by advancedsolid-state NMR spectroscopy. Most of the discotic mesogensunder investigation are hexa-peri-hexabenzocoronene (HBC)derivatives which are of particular interest for potentialdevice applications due to their high one-dimensional chargecarrier mobilities. The supramolecular stacking arrangement of the discoticcores was investigated by 2D 1H-1H double-quantum (DQ)methods, which were modified by incorporating the WATERGATEsuppression technique into the experiments in order toovercome severe phase problems arising from the strongsignal of the long alkyl sidechains. Molecular dynamics and sample orientation was probed throughthe generation of sideband patterns by reconversion rotorencoding in 2D recoupling experiments. These experimentswere extended by new recoupling schemes to enable thedistinction of motion and orientation effects. The solid-state NMR studies presented in this work aim tothe understanding of structure-property relationships in theinvestigated discotic materials, while the experimentsapplied to these materials include new recoupling schemeswhich make the desired information on molecular orientationand dynamics accessible without isotope labelling.

Dynamics, rheology and critical properties of colloidal fluid mixtures

Liquids under the influence of external fields exhibit a wide range of intriguing phenomena that can be markedly different from the behaviour of a quiescent system. This work considers two different systems — a glassforming Yukawa system and a colloid-polymer mixture — by Molecular Dynamics (MD) computer simulations coupled to dissipative particle dynamics. The former consists of a 50-50 binary mixture of differently-sized, like-charged colloids interacting via a screened Coulomb (Yukawa) potential. Near the glass transition the influence of an external shear field is studied. In particular, the transition from elastic response to plastic flow is of interest. At first, this model is characterised in equilibrium. Upon decreasing temperature it exhibits the typical dynamics of glassforming liquids, i.e. the structural relaxation time τα grows strongly in a rather small temperature range. This is discussed with respect to the mode-coupling theory of the glass transition (MCT). For the simulation of bulk systems under shear, Lees-Edwards boundary conditions are applied. At constant shear rates γ˙ ≫ 1/τα the relevant time scale is given by 1/γ˙ and the system shows shear thinning behaviour. In order to understand the pronounced differences between a quiescent system and a system under shear, the response to a suddenly commencing or terminating shear flow is studied. After the switch-on of the shear field the shear stress shows an overshoot, marking the transition from elastic to plastic deformation, which is connected to a super-diffusive increase of the mean squared displacement. Since the average static structure only depends on the value of the shear stress, it does not discriminate between those two regimes. The distribution of local stresses, in contrast, becomes broader as soon as the system starts flowing. After a switch-off of the shear field, these additional fluctuations are responsible for the fast decay of stresses, which occurs on a time scale 1/γ˙ . The stress decay after a switch-off in the elastic regime, on the other hand, happens on the much larger time scale of structural relaxation τα. While stresses decrease to zero after a switch-off for temperatures above the glass transition, they decay to a finite value for lower temperatures. The obtained results are important for advancing new theoretical approaches in the framework of mode-coupling theory. Furthermore, they suggest new experimental investigations on colloidal systems. The colloid-polymer mixture is studied in the context of the behaviour near the critical point of phase separation. For the MD simulations a new effective model with soft interaction potentials is introduced and its phase diagram is presented. Here, mainly the equilibrium properties of this model are characterised. While the self-diffusion constants of colloids and polymers do not change strongly when the critical point is approached, critical slowing down of interdiffusion is observed. The order parameter fluctuations can be determined through the long-wavelength limit of static structure factors. For this strongly asymmetric mixture it is shown how the relevant structure factor can be extracted by a diagonalisation of a matrix that contains the partial static structure factors. By presenting first results of this model under shear it is demonstrated that it is suitable for non-equilibrium simulations as well.

Coarse-graining methods and polymer dynamics

This thesis studies molecular dynamics simulations on two levels of resolution: the detailed level of atomistic simulations, where the motion of explicit atoms in a many-particle system is considered, and the coarse-grained level, where the motion of superatoms composed of up to 10 atoms is modeled. While atomistic models are capable of describing material specific effects on small scales, the time and length scales they can cover are limited due to their computational costs. Polymer systems are typically characterized by effects on a broad range of length and time scales. Therefore it is often impossible to atomistically simulate processes, which determine macroscopic properties in polymer systems. Coarse-grained (CG) simulations extend the range of accessible time and length scales by three to four orders of magnitude. However, no standardized coarse-graining procedure has been established yet. Following the ideas of structure-based coarse-graining, a coarse-grained model for polystyrene is presented. Structure-based methods parameterize CG models to reproduce static properties of atomistic melts such as radial distribution functions between superatoms or other probability distributions for coarse-grained degrees of freedom. Two enhancements of the coarse-graining methodology are suggested. Correlations between local degrees of freedom are implicitly taken into account by additional potentials acting between neighboring superatoms in the polymer chain. This improves the reproduction of local chain conformations and allows the study of different tacticities of polystyrene. It also gives better control of the chain stiffness, which agrees perfectly with the atomistic model, and leads to a reproduction of experimental results for overall chain dimensions, such as the characteristic ratio, for all different tacticities. The second new aspect is the computationally cheap development of nonbonded CG potentials based on the sampling of pairs of oligomers in vacuum. Static properties of polymer melts are obtained as predictions of the CG model in contrast to other structure-based CG models, which are iteratively refined to reproduce reference melt structures. The dynamics of simulations at the two levels of resolution are compared. The time scales of dynamical processes in atomistic and coarse-grained simulations can be connected by a time scaling factor, which depends on several specific system properties as molecular weight, density, temperature, and other components in mixtures. In this thesis the influence of molecular weight in systems of oligomers and the situation in two-component mixtures is studied. For a system of small additives in a melt of long polymer chains the temperature dependence of the additive diffusion is predicted and compared to experiments.

A soft-tetramer model for diblock copolymer melts: structure formation and geometric characterization

The ability of block copolymers to spontaneously self-assemble into a variety of ordered nano-structures not only makes them a scientifically interesting system for the investigation of order-disorder phase transitions, but also offers a wide range of nano-technological applications. The architecture of a diblock is the most simple among the block copolymer systems, hence it is often used as a model system in both experiment and theory. We introduce a new soft-tetramer model for efficient computer simulations of diblock copolymer melts. The instantaneous non-spherical shape of polymer chains in molten state is incorporated by modeling each of the two blocks as two soft spheres. The interactions between the spheres are modeled in a way that the diblock melt tends to microphase separate with decreasing temperature. Using Monte Carlo simulations, we determine the equilibrium structures at variable values of the two relevant control parameters, the diblock composition and the incompatibility of unlike components. The simplicity of the model allows us to scan the control parameter space in a completeness that has not been reached in previous molecular simulations.The resulting phase diagram shows clear similarities with the phase diagram found in experiments. Moreover, we show that structural details of block copolymer chains can be reproduced by our simple model.We develop a novel method for the identification of the observed diblock copolymer mesophases that formalizes the usual approach of direct visual observation,using the characteristic geometry of the structures. A cluster analysis algorithm is used to determine clusters of each component of the diblock, and the number and shape of the clusters can be used to determine the mesophase.We also employ methods from integral geometry for the identification of mesophases and compare their usefulness to the cluster analysis approach.To probe the properties of our model in confinement, we perform molecular dynamics simulations of atomistic polyethylene melts confined between graphite surfaces. The results from these simulations are used as an input for an iterative coarse-graining procedure that yields a surface interaction potential for the soft-tetramer model. Using the interaction potential derived in that way, we perform an initial study on the behavior of the soft-tetramer model in confinement. Comparing with experimental studies, we find that our model can reflect basic features of confined diblock copolymer melts.

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.

Mode-coupling theory: generalizations, high dimensions and microscopic dynamics

This work contains several applications of the mode-coupling theory (MCT) and is separated into three parts. In the first part we investigate the liquid-glass transition of hard spheres for dimensions d→∞ analytically and numerically up to d=800 in the framework of MCT. We find that the critical packing fraction ϕc(d) scales as d²2^(-d), which is larger than the Kauzmann packing fraction ϕK(d) found by a small-cage expansion by Parisi and Zamponi [J. Stat. Mech.: Theory Exp. 2006, P03017 (2006)]. The scaling of the critical packing fraction is different from the relation ϕc(d)∼d2^(-d) found earlier by Kirkpatrick and Wolynes [Phys. Rev. A 35, 3072 (1987)]. This is due to the fact that the k dependence of the critical collective and self nonergodicity parameters fc(k;d) and fcs(k;d) was assumed to be Gaussian in the previous theories. We show that in MCT this is not the case. Instead fc(k;d) and fcs(k;d), which become identical in the limit d→∞, converge to a non-Gaussian master function on the scale k∼d^(3/2). We find that the numerically determined value for the exponent parameter λ and therefore also the critical exponents a and b depend on the dimension d, even at the largest evaluated dimension d=800. In the second part we compare the results of a molecular-dynamics simulation of liquid Lennard-Jones argon far away from the glass transition [D. Levesque, L. Verlet, and J. Kurkijärvi, Phys. Rev. A 7, 1690 (1973)] with MCT. We show that the agreement between theory and computer simulation can be improved by taking binary collisions into account [L. Sjögren, Phys. Rev. A 22, 2866 (1980)]. We find that an empiric prefactor of the memory function of the original MCT equations leads to similar results. In the third part we derive the equations for a mode-coupling theory for the spherical components of the stress tensor. Unfortunately it turns out that they are too complex to be solved numerically.

Molekulardynamik-Simulationen supramolekularer Komplexe unter Einfluss einer externen Kraft

Supramolekulare Komplexe werden durch nichtkovalente Bindungen stabilisiert. Legt man eine externe Kraft an einen solchen Komplex an, ist es möglich, diese Bindungen zu öffnen. Anhand der dafür benötigten Kraft läßt sich die Stabilität des Komplexes bestimmen. Im Rahmen dieser Arbeit wurden zwei supramolekulare Komplexe, die unterschiedliche Arten von nichtkovalenten Bindungen enthalten, mit Hilfe von Molekulardynamik (MD) Simulationen untersucht. In beiden Fällen wurden die relevanten Bindungsstrukturen und deren Stabilität ermittelt.rnZum einen wurden zwei synthetische Calix[4]aren-Catenan-Dimersysteme betrachtet, in denen die beiden Monomere über Wasserstoffbrückenbindungen aneinander gebunden sind. Die Besonderheit dieser Komplexe ist, dass die Monomere aufgrund von verschlauften Alkylketten (Catenan-Struktur) nicht vollständig voneinander getrennt werden können. In Abhängigkeit der Länge derrnAlkylketten findet man für die beiden Komplexe eine unterschiedliche Zahl von relevanten Bindungsstrukturen (Zustände). Für ein System mit relativ kurzen Alkylketten findet man zwei Zustände, eine kompakte Struktur, die auch im Gleichgewicht beobachtet wird und eine gestreckte Struktur, die nur unter dem Einfluss der externen Kraft stabil ist. Verlängert man die Alkylketten,rnbeobachtet man einen weiteren Zustand, in dem das Dimer vollständig gestreckt ist und die Monomere eine größere Separation aufweisen.rnBei dem zweiten System, das untersucht wurde, handelte es sich um einen Carbohydrat-Kation-Carbohydrat Komplex, der für die Selbstadhäsion von Meeresschwämmen eine wichtige Rolle spielt. Experimentell ist bekannt, dass sich dieser Komplex zwar mit Calciumionen, nicht aber mit Magnesiumionen bildet. Im Rahmen dieser Arbeit wurde gezeigt, dass die wesentlichen Unterschiede der beiden Kationarten in Bezug auf die Komplexbildung auf den kleineren Ionenradius des Magnesiumions zurückzuführen sind. Aufgrund des kleineren Radius bindet ein solvatisiertes Magnesiumion die Hydrathülle stärker und die Komplexbindung wird kinetisch gehemmt. Zum anderen bindet im Magnesiumkomplex nur eines der beiden Carbohydrate direkt an das Kation.rnDas andere Carbohydrat bindet nur indirekt über ein Wassermolekül an das Kation. Da diese indirekte Bindung gegenüber einer direkten Bindung schwächer ist, weist der Magensiumkomplex eine geringere Stabilität auf als ein vergleichbarer Calciumkomplex.rnDes Weiteren wurde untersucht, inwieweit die Ergebnisse von MD Simulationen vom verwendeten Modell (Kraftfeld) abhängen. Allgemein ist bekannt, dass die Ergebnisse von Gleichgewichtssimulationen kraftfeldabhängig sind. Im Rahmen diese Arbeit konnte gezeigt werden, dass sich für Zugsimulationen, in denen eine externe Kraft an das System angelegt wird, eine ähnliche Kraftfeldabhängigkeit ergibt. Da sich die Unterschiede der Ergebnisse auf Unterschiede in den Gleichgewichtssimulationen zurückführen lassen, kann man annehmen, dass die externe Kraft keine zusätzliche Einschränkung in Bezug auf die Zuverlässigkeit der Kraftfelder darstellt.rnAbgesehen von den MD Simulationen wurde eine in der Literatur beschriebene Methode zur Analyse von Zweizustandssystemen unter dem Einfluss einer konstanten externen Kraft erweitert. Ein Komplex läßt sich als Zweizustandssystem beschreiben, wenn dieser zwei relevante Bindungsstrukturen aufweist. Wird an solch einen Komplex eine konstante Kraft angelegt, lassen sich Übergänge zwischen den beiden Zuständen beobachten. Ist das System weit entfernt vom Gleichgewicht, kann es problematisch sein, einen der beiden Übergänge vollständig aufzulösen. In diesen Fällen wird nun vorgeschlagen, die beiden Übergänge zu einem sogenannten Kreisübergang zusammen zu fassen und diese zu zählen. Bestimmt man die Zahl der Übergänge pro Zeit in Abhängigkeit der angelegten Kraft, können die Übergangsraten bestimmt werden. Um die Methode zu validieren wurden kinetische Monte-Carlo Simulationen durchgeführt. Es zeigt sich, dass schon mit relativ kleinen Datensätzen gute Ergebnisse erzielt werden können.

Controlled wetting on silica-based nano- and microstructured surfaces

Der Fokus dieser Doktorarbeit liegt auf der kontrollierten Benetzung von festen Oberflächen, die in vielen Bereichen, wie zum Beispiel in der Mikrofluidik, für Beschichtungen und in biologischen Studien von Zellen oder Bakterien, von großer Bedeutung ist.rnDer erste Teil dieser Arbeit widmet sich der Frage, wie Nanorauigkeit das Benetzungsverhalten, d.h. die Kontaktwinkel und die Pinningstärke, von hydrophoben und superhydrophoben Beschichtungen beeinflusst. Hierfür wird eine neue Methode entwickelt, um eine nanoraue Silika-Beschichtung über die Gasphase auf eine superhydrophobe Oberfläche, die aus rauen Polystyrol-Silika-Kern-Schale-Partikeln besteht, aufzubringen. Es wird gezeigt, dass die Topographie und Dichte der Nanorauigkeiten bestimmt, ob sich die Superhydrophobizität verringert oder erhöht, d.h. ob sich ein Flüssigkeitstropfen im Nano-Wenzel- oder Nano-Cassie-Zustand befindet. Das verstärkte Pinning im Nano-Wenzel-Zustand beruht auf dem Eindringen von Flüssigkeitsmolekülen in die Nanoporen der Beschichtung. Im Nano-Cassie-Zustand dagegen sitzt der Tropfen auf den Nanorauigkeiten, was das Pinning vermindert. Die experimentellen Ergebnisse werden mit molekulardynamischen Simulationen in Bezug gesetzt, die den Einfluss der Oberflächenbeschichtungsdichte und der Länge von fluorinierten Silanen auf die Hydrophobizität einer Oberfläche untersuchen. rnEs wurden bereits verschiedenste Techniken zur Herstellung von transparenten superhydrophoben, d.h. extrem flüssigkeitsabweisenden, Oberflächen entwickelt. Eine aktuelle Herausforderung liegt darin, Funktionalitäten einzuführen, ohne die superhydrophoben Eigenschaften einer Oberfläche zu verändern. Dies ist extrem anspruchsvoll, da funktionelle Gruppen in der Regel hydrophil sind. In dieser Arbeit wird eine innovative Methode zur Herstellung von transparenten superhydrophoben Oberflächen aus Janus-Mikrosäulen mit variierenden Dimensionen und Topographien entwickelt. Die Janus-Säulen haben hydrophobe Seitenwände und hydrophile Silika-Oberseiten, die anschließend selektiv und ohne Verlust der superhydrophoben Eigenschaften der Oberfläche funktionalisiert werden können. Diese selektive Oberflächenfunktionalisierung wird mittels konfokaler Mikroskopie und durch das chemische Anbinden von fluoreszenten Molekülen an die Säulenoberseiten sichtbar gemacht. Außerdem wird gezeigt, dass das Benetzungsverhalten durch Wechselwirkungen zwischen Flüssigkeit und Festkörper in der Nähe der Benetzungslinie bestimmt wird. Diese Beobachtung widerlegt das allgemein akzeptierte Modell von Cassie und Baxter und beinhaltet, dass hydrophile Flächen, die durch mechanischen Abrieb freigelegt werden, nicht zu einem Verlust der Superhydrophobizität führen müssen, wie allgemein angenommen.rnBenetzung kann auch durch eine räumliche Beschränkung von Flüssigkeiten kontrolliert werden, z.B. in mikrofluidischen Systemen. Hier wird eine modifizierte Stöber-Synthese verwendet, um künstliche und natürliche Faser-Template mit einer Silika-Schicht zu ummanteln. Nach der thermischen Zersetzung des organischen Templat-Materials entstehen wohldefinierte Silika-Kanäle und Kanalkreuzungen mit gleichmäßigen Durchmessern im Nano- und Mikrometerbereich. Auf Grund ihrer Transparenz, mechanischen Stabilität und des großen Länge-zu-Durchmesser-Verhältnisses sind die Kanäle sehr gut geeignet, um die Füllgeschwindigkeiten von Flüssigkeiten mit variierenden Oberflächenspannungen und Viskositäten zu untersuchen. Konfokale Mikroskopie ermöglicht es hierbei, die Füllgeschwindigkeiten über eine Länge von mehreren Millimetern, sowie direkt am Kanaleingang zu messen. Das späte Füllstadium kann sehr gut mit der Lucas-Washburn-Gleichung beschrieben werden. Die anfänglichen Füllgeschwindigkeiten sind jedoch niedriger als theoretisch vorhergesagt. Wohingegen die vorhergehenden Abschnitte dieser Arbeit sich mit der quasistatischen Benetzung beschäftigen, spielt hier die Dynamik der Benetzung eine wichtige Rolle. Tatsächlich lassen sich die beobachteten Abweichungen durch einen geschwindigkeitsabhängigen Fortschreitkontaktwinkel erklären und durch dynamische Benetzungstheorien modellieren. Somit löst diese Arbeit das seit langem diskutierte Problem der Abweichungen von der Lucas-Washburn-Gleichung bei kleinen Füllgeschwindigkeiten.

Untersuchungen zum nichtphotochemischen Lochbrennen auf Einzelmolekülebene am Modellsystem Terrylen in p-Terphenyl

Die Arbeit beschreibt Untersuchungen zum nichtphoto- chemischen Lochbrennen, das bei 1.4 Kelvin in Form von rein lichtinduzierten Frequenzsprüngen einzelner in p-Terphenyleingebetteter Terrylenmoleküle beobachtet werden kann. Dabei zeigen alle Chromophore aus der X1-Einbaulage ein exzellent reproduzierbares Verhalten, sowohl im bistabilen primären Photozyklus wie auch in dem daran angegliederten sekundärenPhotozyklus, welcher aus drei weiteren spektralen Positionen besteht. Aus den Ergebnissen der nach der genauen Charakterisierung dieser Eigenschaft des Systems durchgeführten Experimente - Fluoreszenzspektroskopie der Photoprodukte, Stark-Effekt-Messungen und Polarisationsmodulation - wird ein Modell für die den lichtinduzierten Änderungen der Absorptionsfrequenzzugrundeliegenden Konformationsänderungender Wirt/Gast- Struktur abgeleitet und diskutiert. Die mittlerweile verfügbaren Ergebnisse von diesbezüglichen molekular- dynamischen Simulationen einer Theoriegruppe ausBordeaux, die alle grundlegenden Annahmen dieses Modellsbestätigen und eine noch genauere mikroskopische Beschreibung des Systems liefern, werden zur Abrundung der Darstellung ebenfalls vorgestellt. Außerdem geht die Dissertation auf die durchgeführten Einzelmolekül- untersuchungen an Terrylen in p-Terphenyl bei Raumtemperatur ein und stellt das im Rahmen der Arbeit aufgebaute temperaturvariable laserscannende Konfokalmikroskop im Detail vor.

Dynamik unterkühlter Flüssigkeiten in Filmen und Röhren

GERMAN:Im Rahmen der vorliegenden Arbeit soll der Einfluß einerräumlichen Beschränkung auf die Dynamik einer unterkühltenFlüssigkeit charakterisiert werden. Insbesondere sollgeklärt werden, welche Rolle die Kooperativität derTeilchenbewegung bei niedrigen Temperaturen spielt. Hierzuuntersuchen wir mit Hilfe einer Molekulardynamik-Computersimulation die dynamischen Eigenschaften eineseinfachen Modellglasbildners, einer binäre Lennard-Jones-Flüssigkeit, für Systeme mit unterschiedlichen Geometrienund Wandarten. Durch geschickte Wahl der Wandpotentiale konnte erreichtwerden, daß die Struktur der Flüssigkeit mit der im Bulknahezu identisch ist.In Filmen mit glatten Wänden beobachtet man, daß dieDynamik der Flüssigkeit in der Nähe der Wand starkbeschleunigt ist und sich diese veränderte Dynamik bis weitin den Film ausbreitet. Den umgekehrten Effekt erhält man,wenn man eine strukturierte, rauhe Wand verwendet, in derenNähe die Dynamik stark verlangsamt ist.Die kontinuierliche Verlangsamung bzw. Beschleunigung derDynamik vom Verhalten an der Oberfläche zum Bulkverhaltenin genügend großem Abstand zur Wand können wirphänomenologisch beschreiben. Hieraus kann mancharakteristische dynamische Längenskalen ablesen, die mitsinkender Temperatur kontinuierlich anwachsen, d.h. derBereich, in dem die Existenz der Wand einen (indirekten)Einfluß auf die Dynamik eines Flüssigkeitsteilchens hat,breitet sich immer weiter aus. Man kann daher vonBereichen kooperativer Bewegung sprechen, die mit sinkenderTemperatur anwachsen.Unsere Untersuchungen von Röhren zeigen, daß aufgrund desstärkeren Einflusses der Wände die beobachteten Effektegrößer sind als in Filmgeometrie. Bei Reduzierung derSystemgröße zeigen sich immer größere Unterschiede zumBulkverhalten.

2 H-NMR und Viskositätsuntersuchungen zum Studium der molekularen Dynamik unterkühlter Flüssigkeiten unter Hochdruck

In der vorliegenden Arbeit wurde die Druckabhängigkeit der molekularen Dynamik mittels 2H-NMR und Viskositätsmessungen untersucht. Für die Messungen wurde der niedermolekulare organische Glasbildner ortho-Terphenyl (OTP) ausgewählt, da dieser aufgrund einer Vielzahl vorliegender Arbeiten als Modellsubstanz angesehen werden kann. Daneben wurden auch Messungen an Salol durchgeführt.Die Untersuchungen erstreckten sich über einen weiten Druck- und Temperaturbereich ausgehend von der Schmelze bis weit in die unterkühlte Flüssigkeit. Dieser Bereich wurde aufgrund experimenteller Voraussetzungen immer durch eine Druckerhöhung erreicht.Beide Substanzen zeigten druckabhängig ein Verhalten, das dem der Temperaturvariation bei Normaldruck sehr ähnelt. Auf einer Zeitskala der molekularen Dynamik von 10E-9 s bis zu 10E+2 s wurde daher am Beispiel von OTP ein Druck-Temperatur-Zeit-Superpositionsprinzip diskutiert. Zudem konnte eine Temperatur-Dichte-Skalierung mit rho T-1/4 erfolgreich durchgeführt werden. Dies entspricht einem rein repulsiven Potentialverlauf mit rho -12±3 .Zur Entscheidung, ob die Verteilungsbreiten der mittleren Rotationskorrelationszeiten durch Druckvariation beeinflußt werden, wurden auch Ergebnisse anderer experimenteller Methoden herangezogen. Unter Hinzuziehung aller Meßergebnisse kann sowohl eine Temperatur- als auch Druckabhängigkeit der Verteilungsbreite bestätigt werden. Zur Auswertung von Viskositätsdaten wurde ein Verfahren vorgestellt, das eine quantitative Aussage über den Fragilitätsindex von unterkühlten Flüssigkeiten auch dann zuläßt, wenn die Messungen nicht bis zur Glasübergangstemperatur Tg durchgeführt werden. Die Auswertung der druckabhängigen Viskositätsdaten von OTP und Salol zeigt einen sehr differenzierten druckabhängigen Verlauf des Fragilitätsindexes für beide Glasbildner. OTP zeigt zunächst eine leichte Abnahme und danach wieder eine Zunahme des Fragilitätsindexes, dieses Ergebnis wird auch von Simulationsdaten, die der Literatur entnommen wurden, unterstützt. Salol hingegen zeigt zunächst eine deutliche Zunahme und danach eine Abnahme des Fragilitätsindexes. Das unterschiedliche Verhalten der beiden Glasbildner mit ähnlichem Fragilitätsindex bei Normaldruck wird auf die Wasserstoffbrückenbindungen innerhalb von Salol zurückgeführt.

Molekulardynamik-Untersuchungen zur atomistischen Struktur & Dynamik von binären Mischgläsern

Mit Hilfe von Molekulardynamik-Computersimulationen werdenin dieser Arbeit die Struktur und Dynamik von Gläsern undSchmelzen der Systeme 'NSx'=(Na2O)(xSiO2), mit x=2,3,5, und'AS2'=(Al2O3)(2SiO2) untersucht. Zur Beschreibung dermikroskopischen Wechselwirkungen dient ein Modellpotenzial,das auf einem effektiven Paarpotenzial aus der Literaturaufbaut. Simuliert wurden Teilchentrajektorien über mehrereNanosekunden im Bereich 6100 K >= T >= 2100 K sowie dieGlasstruktur bei 100 K (NSx) bzw. 300 K (AS2). Das Aufbrechen der tetraedrischen Netzwerkstruktur durch denNetzwerkmodifikator Na2O führt zur Ausbildung einerzusätzlichen Struktur auf intermediären Längenskalen, diebei allen Systemen NSx etwa dem Abstand übernächster Na-bzw. Si-Nachbarn eines Na-Atoms entspricht. Die diffusiveDynamik ist in allen Systeme NSx bis zu drei Größenordnungenschneller als in SiO2 . Sie nimmt mit wachsenderNa-Konzentration zu. Die Na-Diffusion zeigtArrhenius-Verhalten; hierbei vollführen die Na-Atome einaktiviertes Hüpfen durch eine erstarrte Si-O-Matrix. DieZustandsdichten werden bis ca. 20 THz durch dominanteNa-Moden bestimmt. Bei hohen Frequenzen weichen die für SiO2 typischen intratetraedrischen Si-O-Schwingungsmodenauf.Im Gegensatz zu Na2O wird Al2O3 in die tetraedrischeNetzwerkstruktur eingebunden. AS2 zeigt eine überwiegend ausAlO4- und SiO4-Tetraedern verknüpfte Polyederstruktur, beider sich die AlO4-Tetraeder lokal anders anordnen als dieSiO4-Tetraeder, um Ladungsausgleich zu gewährleisten. Esbilden sich typische 3(Si,Al)O4-Bausteine ('3-Cluster'). Sie führen auf intermediären Längenskalen zur Ausbildungeines Al-reichen perkolierenden Netzwerks, das dieSiO4-Struktur durchdringt. Im Vergleich zu SiO2 erhöht sichdie diffusive Dynamik aller Komponenten in AS2 um ca. zweiGrößenordnungen. Die intratetraedrischenAl-O-Steckschwingungsmoden sind wesentlich weicher als die typischen intratetraedrischen Si-O-Moden des SiO2.

Novel simulation methods for Coulomb and hydrodynamic interactions

This thesis presents new methods to simulate systems with hydrodynamic and electrostatic interactions. Part 1 is devoted to computer simulations of Brownian particles with hydrodynamic interactions. The main influence of the solvent on the dynamics of Brownian particles is that it mediates hydrodynamic interactions. In the method, this is simulated by numerical solution of the Navier--Stokes equation on a lattice. To this end, the Lattice--Boltzmann method is used, namely its D3Q19 version. This model is capable to simulate compressible flow. It gives us the advantage to treat dense systems, in particular away from thermal equilibrium. The Lattice--Boltzmann equation is coupled to the particles via a friction force. In addition to this force, acting on {it point} particles, we construct another coupling force, which comes from the pressure tensor. The coupling is purely local, i.~e. the algorithm scales linearly with the total number of particles. In order to be able to map the physical properties of the Lattice--Boltzmann fluid onto a Molecular Dynamics (MD) fluid, the case of an almost incompressible flow is considered. The Fluctuation--Dissipation theorem for the hybrid coupling is analyzed, and a geometric interpretation of the friction coefficient in terms of a Stokes radius is given. Part 2 is devoted to the simulation of charged particles. We present a novel method for obtaining Coulomb interactions as the potential of mean force between charges which are dynamically coupled to a local electromagnetic field. This algorithm scales linearly, too. We focus on the Molecular Dynamics version of the method and show that it is intimately related to the Car--Parrinello approach, while being equivalent to solving Maxwell's equations with freely adjustable speed of light. The Lagrangian formulation of the coupled particles--fields system is derived. The quasi--Hamiltonian dynamics of the system is studied in great detail. For implementation on the computer, the equations of motion are discretized with respect to both space and time. The discretization of the electromagnetic fields on a lattice, as well as the interpolation of the particle charges on the lattice is given. The algorithm is as local as possible: Only nearest neighbors sites of the lattice are interacting with a charged particle. Unphysical self--energies arise as a result of the lattice interpolation of charges, and are corrected by a subtraction scheme based on the exact lattice Green's function. The method allows easy parallelization using standard domain decomposition. Some benchmarking results of the algorithm are presented and discussed.

Tail bridging vs. patchiness

The subject of this thesis are the interactions between nucleosome core particles (NCPs). NCPs are the primary storage units of DNA in eucaryotic cells. Each NCP consists of a core of eight histone proteins and a strand of DNA, which is wrapped around about two times. Each histone protein has a terminal tail passing over and between the superhelix of the wrapped DNA. Special emphasis was placed on the role of the histone tails, since experimental ndings suggest that the tails have a great in uence on the mutual attraction of the NCPs. In those experiments Mangenot et al. observe a dramatic change in the con guration of the tails, which is accompanied by evidence of mutual attraction between NCPs, when a certain salt concentration is reached. Existing models used in the theoretical approaches and in simulations focus on the description of the histone core and the wrapped DNA, but neglect the histone tails. We introduce the multi chain complex as a new simulation model. Here the histone core and the wrapping DNA are modelled via a charged sphere, while the histone tails are represented by oppositely charged chains grafted on the sphere surface. We start by investigating the parameter space describing a single NCP. The Debye-Huckel potential is used to model the electrostatic interactions and to determine the e ective charge of the NCP core. This value is subsequently used for a study of the pairinteraction of two NCPs via an extensive Molecular Dynamics study. The monomer distribution of the full chain model is investigated. The existence of tail bridges between the cores is demonstrated. Finally, by discriminating between bridging and non-bridging con gurations, we can show that the effect of tail bridging between the spheres does indeed account for the observed attraction. The full chain model can serve as a model to study the acetylation of the histone tails of the nucleosome. The reduction of the charge fraction of the tails, that corresponds to the process of acetylation, leads to a reduction or even the disappearance of the attraction. A recent MC study links this e ect to the unfolding of the chromatin ber in the case of acetylated histone tails. In this case the acetylation of the histone tails leads to the formation of heterochromatin, and one could understand how larger regions of the genetic information could be inactivated through this mechanism.

Liquid crystalline oligofluorenes and their derivatives

Homo-oligofluorenes (OFn), polyfluorenes (PF2/6) and oligofluorenes with one fluorenenone group in the center (OFnK) were synthesized. They were used as model compounds to understand of the structure-property relationships of polyfluorenes and the origin of the green emission in the photoluminescence (after photooxidation of the PFs) and the electroluminescence (EL) spectra. The electronic, electrochemical properties, thermal behavior, supramolecular self-assembly, and photophysical properties of OFn, PF2/6 and OFnK were investigated. Oligofluorenes with 2-ethylhexyl side chain (OF2-OF7) from the dimer up to the heptamer were prepared by a series of stepwise transition metal mediated Suzuki and Yamamoto coupling reactions. Polyfluorene was synthesized by Yamamoto coupling of 2,7-dibromo-9,9-bis(2-ethylhexyl)fluorene. Oligofluorenes with one fluorenone group in the center (OF3K, OF5K, OF7K) were prepared by Suzuki coupling between the monoboronic fluorenyl monomer, dimer, trimer and 2, 7-dibromofluorenone. The electrochemical and electronic properties of homo-oligofluorenes (OFn) were systematically studied by several combined techniques such as cyclic voltammetry, differential pulse voltammetry, UV-vis absorption spectroscopy, steady and time-resolved fluorescence spectroscopy. It was found that the oligofluorenes behave like classical conjugated oligomers, i.e., with the increase of the chain-length, the corresponding oxidation potential, the absorption and emission maximum, ionization potential, electron affinity, band gap and the photoluminescence lifetime displayed a very good linear relation with the reciprocal number of the fluorene units (1/n). The extrapolation of these linear relations to infinite chain length predicted the electrochemical and electronic properties of the corresponding polyfluorenes. The thermal behavior, single-crystal structure and supramolecular packing, alignment properties, and molecular dynamics of the homo-oligofluorenes (OFn) up to the polymer were studied using techniques such as TGA, DSC, WAXS, POM and DS. The OFn from tetramer to heptamer show a smectic liquid crystalline phase with clearly defined isotropization temperature. The oligomers do show a glass transition which exhibits n-1 dependence and allows extrapolation to a hypothetical glass transition of the polymer at around 64 °C. A smectic packing and helix-like conformation for the oligofluorenes from tetramer to heptamer was supported by WAXS experiments, simulation, and single-crystal structure of some oligofluorene derivatives. Oligofluorenes were aligned more easily than the corresponding polymer, and the alignability increased with the molecular length from tetramer to heptamer. The molecular dynamics in a series of oligofluorenes up to the polymer was studied using dielectric spectroscopy. The photophysical properties of OFn and PF2/6 were investigated by the steady-state spectra (UV-vis absorption and fluorescence spectra) and time-resolved fluorescence spectra both in solution and thin film. The time-resolved fluorescence spectra of the oligofluorenes were measured by streak camera and gate detection technique. The lifetime of the oligofluorenes decreased with the extension of the chain-length. No green emission was observed in CW, prompt and delayed fluorescence for oligofluorenes in m-THF and film at RT and 77K. Phosphorescence was observed for oligofluorenes in frozen dilute m-THF solution at 77K and its lifetime increased with length of oligofluorenes. A linear relation was obtained for triplet energy and singlet energy as a function of the reciprocal degree of polymerization, and the singlet-triplet energy gap (S1-T1) was found to decrease with the increase of degree of polymerization. Oligofluorenes with one fluorenone unit at the center were used as model compounds to understand the origin of the low-energy (“green”) emission band in the photoluminescence and electroluminescence spectra of polyfluorenes. Their electrochemical properties were investigated by CV, and the ionization potential (Ip) and electron affinity (Ea) were calculated from the onset of oxidation and reduction of OFnK. The photophysical properties of OFnK were studied in dilute solution and thin film by steady-state spectra and time-resolved fluorescence spectra. A strong green emission accompanied with a weak blue emission were obtained in solution and only green emission was observed on film. The strong green emission of OFnK suggested that rapid energy transfer takes place from higher energy sites (fluorene segments) to lower energy sites (fluorenone unit) prior to the radiative decay of the excited species. The fluorescence spectra of OFnK also showed solvatochromism. Monoexponential decay behaviour was observed by time-resolved fluorescence measurements. In addition, the site-selective excitation and concentration dependence of the fluorescence spectra were investigated. The ratio of green and blue emission band intensities increases with the increase of the concentration. The observed strong concentration dependence of the green emission band in solution suggests that increased interchain interactions among the fluorenone-containing oligofluorene chain enhanced the emission from the fluorenone defects at higher concentration. On the other hand, the mono-exponential decay behaviour and power dependence were not influenced significantly by the concentration. We have ruled out the possibility that the green emission band originates from aggregates or excimer formation. Energy transfer was further investigated using a model system of a polyfluorene doped by OFnK. Förster-type energy transfer took place from PF2/6 to OFnK, and the energy transfer efficiency increased with increasing of the concentration of OFnK. Efficient funneling of excitation energy from the high-energy fluorene segments to the low-energy fluorenone defects results from energy migration by hopping of excitations along a single polymer chain until they are trapped on the fluorenone defects on that chain or transferred onto neighbouring chains by Förster-type interchain energy transfer process. These results imply that the red-shifted emission in polyfluorenes can originate from (usually undesirable) keto groups at the bridging carbon atoms-especially if the samples have been subject to photo- or electro-oxidation or if fluorenone units are present due to an improper purification of the monomers prior to polymerization.