Theoretical studies of microfluidics: "walking" drops and "melting" snowflakes.

The present PhD thesis summarizes two examples of research in microfluidics. Both times water was the subject of interest, once in the liquid state (droplets adsorbed on chemically functionalized surfaces), the other time in the solid state (ice snowflakes and their fractal behaviour). The first problem deals with a slipping nano-droplet of water adsorbed on a surface with photo-switchable wettability characteristics. Main focus was on identifying the underlying driving forces and mechanical principles at the molecular level of detail. Molecular Dynamics simulation was employed as investigative tool owing to its record of successfully describing the microscopic behaviour of liquids at interfaces. To reproduce the specialized surface on which a water droplet can effectively ���walk���, a new implicit surface potential was developed. Applying this new method the experimentally observed droplet slippage could be reproduced successfully. Next the movement of the droplet was analyzed at various conditions emphasizing on the behaviour of the water molecules in contact with the surface. The main objective was to identify driving forces and molecular mechanisms underlying the slippage process. The second part of this thesis is concerned with theoretical studies of snowflake melting. In the present work snowflakes are represented by filled von Koch-like fractals of mesoscopic beads. A new algorithm has been developed from scratch to simulate the thermal collapse of fractal structures based on Monte Carlo and Random Walk Simulations (MCRWS). The developed method was applied and compared to Molecular Dynamics simulations regarding the melting of ice snowflake crystals and new parameters were derived from this comparison. Bigger snow-fractals were then studied looking at the time evolution at different temperatures again making use of the developed MCRWS method. This was accompanied by an in-depth analysis of fractal properties (border length and gyration radius) in order to shed light on the dynamics of the melting process.

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.

Atomistic simulations of cosolvent/water mixtures

In this thesis, atomistic simulations are performed to investigate hydrophobic solvation and hydrophobic interactions in cosolvent/water binary mixtures. Many cosolvent/water binary mixtures exhibit non-ideal behavior caused by aggregation at the molecular scale level although they are stable and homogenous at the macroscopic scale. Force-field based atomistic simulations provide routes to relate atomistic-scale structure and interactions to thermodynamic solution properties. The predicted solution properties are however sensitive to the parameters used to describe the molecular interactions. In this thesis, a force field for tertiary butanol (TBA) and water mixtures is parameterized by making use of the Kirkwood-Buff theory of solution. The new force field is capable of describing the alcohol-alcohol, water-water and alcohol-water clustering in the solution as well as the solution components��� chemical potential derivatives in agreement with experimental data. With the new force field, the preferential solvation and the solvation thermodynamics of a hydrophobic solute in TBA/water mixtures have been studied. First, methane solvation at various TBA/water concentrations is discussed in terms of solvation free energy-, enthalpy- and entropy- changes, which have been compared to experimental data. We observed that the methane solvation free energy varies smoothly with the alcohol/water composition while the solvation enthalpies and entropies vary nonmonotonically. The latter occurs due to structural solvent reorganization contributions which are not present in the free energy change due to exact enthalpy-entropy compensation. It is therefore concluded that the enthalpy and entropy of solvation provide more detailed information on the reorganization of solvent molecules around the inserted solute. Hydrophobic interactions in binary urea/water mixtures are next discussed. This system is particularly relevant in biology (protein folding/unfolding), however, changes in the hydrophobic interaction induced by urea molecules are not well understood. In this thesis, this interaction has been studied by calculating the free energy (potential of mean force), enthalpy and entropy changes as a function of the solute-solute distance in water and in aqueous urea (6.9 M) solution. In chapter 5, the potential of mean force in both solution systems is analyzed in terms of its enthalpic and entropic contributions. In particular, contributions of solvent reorganization in the enthalpy and entropy changes are studied separately to better understand what are the changes in interactions in the system that contribute to the free energy of association of the nonpolar solutes. We observe that in aqueous urea the association between nonpolar solutes remains thermodynamically favorable (i.e., as it is the case in pure water). This observation contrasts a long-standing belief that clusters of nonpolar molecules dissolve completely in the presence of urea molecules. The consequences of our observations for the stability of proteins in concentrated urea solutions are discussed in the chapter 6 of the thesis.

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.

Development of empirical potentials for amorphous silica

Development of empirical potentials for amorphous silica Amorphous silica (SiO2) is of great importance in geoscience and mineralogy as well as a raw material in glass industry. Its structure is characterized as a disordered continuous network of SiO4 tetrahedra. Many efforts have been undertaken to understand the microscopic properties of silica by classical molecular dynamics (MD) simulations. In this method the interatomic interactions are modeled by an effective potential that does not take explicitely into account the electronic degrees of freedom. In this work, we propose a new methodology to parameterize such a potential for silica using ab initio simulations, namely Car-Parrinello (CP) method [Phys. Rev. Lett. 55, 2471 (1985)]. The new potential proposed is compared to the BKS potential [Phys. Rev. Lett. 64, 1955 (1990)] that is considered as the benchmark potential for silica. First, CP simulations have been performed on a liquid silica sample at 3600 K. The structural features so obtained have been compared to the ones predicted by the classical BKS potential. Regarding the bond lengths the BKS tends to underestimate the Si-O bond whereas the Si-Si bond is overestimated. The inter-tetrahedral angular distribution functions are also not well described by the BKS potential. The corresponding mean value of theSiOSi angle is found to be ��� 147���, while the CP yields to aSiOSi angle centered around 135���. Our aim is to fit a classical Born-Mayer/Coulomb pair potential using ab initio calculations. To this end, we use the force-matching method proposed by Ercolessi and Adams [Europhys. Lett. 26, 583 (1994)]. The CP configurations and their corresponding interatomic forces have been considered for a least square fitting procedure. The classical MD simulations with the resulting potential have lead to a structure that is very different from the CP one. Therefore, a different fitting criterion based on the CP partial pair correlation functions was applied. Using this approach the resulting potential shows a better agreement with the CP data than the BKS ones: pair correlation functions, angular distribution functions, structure factors, density of states and pressure/density were improved. At low temperature, the diffusion coefficients appear to be three times higher than those predicted by the BKS model, however showing a similar temperature dependence. Calculations have also been carried out on crystalline samples in order to check the transferability of the potential. The equilibrium geometry as well as the elastic constants of ��-quartz at 0 K are well described by our new potential although the crystalline phases have not been considered for the parameterization. We have developed a new potential for silica which represents an improvement over the pair potentials class proposed so far. Furthermore, the fitting methodology that has been developed in this work can be applied to other network forming systems such as germania as well as mixtures of SiO2 with other oxides (e.g. Al2O3, K2O, Na2O).

Synthesis of Modified Amino Acids and Insertion in Peptides and Mimetics. Structural Aspects and Impact on Biological Activity.

I studied the effects exerted by the modifications on structures and biological activities of the compounds so obtained. I prepared peptide analogues containing unusual amino acids such as halogenated, alkylated (S)- or (R)-tryptophans, useful for the synthesis of mimetics of the endogenous opioid peptide endomorphin-1, or 2-oxo-1,3-oxazolidine-4-carboxylic acids, utilized as pseudo-prolines having a clear all-trans configuration of the preceding peptide bond. The latter gave access to a series of constrained peptidomimetics with potential interest in medicinal chemistry and in the field of the foldamers. In particular, I have dedicated much efforts to the preparation of cyclopentapeptides containing D-configured, alfa-, or beta-aminoacids, and also of cyclotetrapeptides including the retro-inverso modification. The conformational analyses confirmed that these cyclic compounds can be utilized as rigid scaffolds mimicking gamma- or beta-turns, allowing to generate new molecular and 3D diversity. Much work has been dedicated to the structural analysis in solution and in the receptor-bound state, fundamental for giving a rationale to the experimentally determined bioactivity, as well as for predicting the activity of virtual compounds (in silico pre-screen). The conformational analyses in solution has been done mostly by NMR (2D gCosy, Roesy, VT, molecular dynamics, etc.). A special section is dedicated to the prediction of plausible poses of the ligands when bound to the receptors by Molecular Docking. This computational method proved to be a powerful tool for the investigation of ligand-receptor interactions, and for the design of selective agonists and antagonists. Another practical use of cyclic peptidomimetics was the synthesis and biological evaluation of cyclic analogues of endomorphin-1 lacking in a protonable amino group. The studies revealed that a inverse type II beta-turn on D-Trp-Phe constituted the bioactive conformation.

Atomic force microscopy of biomimetic systems

This thesis was driven by the ambition to create suitable model systems that mimic complex processes in nature, like intramolecular transitions, such as unfolding and refolding of proteins, or intermolecular interactions between different cell compo-nents. Novel biophysical approaches were adopted by employing atomic force mi-croscopy (AFM) as the main measurement technique due to its broad diversity. Thus, high-resolution imaging, adhesion measurements, and single-molecule force distance experiments were performed on the verge of the instrumental capabilities. As first objective, the interaction between plasma membrane and cytoskeleton, me-diated by the linker protein ezrin, was pursued. Therefore, the adsorption process and the lateral organization of ezrin on PIP2 containing solid-supported membranes were characterized and quantified as a fundament for the establishment of a biomimetic model system. As second component of the model system, actin filaments were coated on functionalized colloidal probes attached on cantilevers, serving as sensor elements. The zealous endeavor of creating this complex biomimetic system was rewarded by successful investigation of the activation process of ezrin. As a result, it can be stated that ezrin is activated by solely binding to PIP2 without any further stimulating agents. Additional cofactors may stabilize and prolong the active conformation but are not essentially required for triggering ezrin���s transformation into an active conformation. In the second project, single-molecule force distance experiments were performed on bis-loop tetra-urea calix[4]arene-catenanes with different loading rates (increase in force per second). These macromolecules were specifically designed to investigate the rupture and rejoining mechanism of hydrogen bonds under external load. The entangled loops of capsule-like molecules locked the unbound state of intramolecular hydrogen bonds mechanically, rendering a rebinding observable on the experimental time scale. In conjunction with Molecular Dynamics simulations, a three-well potential of the bond rupture process was established and all kinetically relevant parameters of the experiments were determined by means of Monte Carlo simulations and stochastic modeling. In summary, it can be stated that atomic force microscopy is an invaluable tool to scrutinize relevant processes in nature, such as investigating activation mechanisms in proteins, as shown by analysis of the interaction between F-actin and ezrin, as well as exploring fundamental properties of single hydrogen bonds that are of paramount interest for the complete understanding of complex supramolecular structures.

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.

Synthesis of hexa-peri-hexabenzocoronenes - from building blocks to functional materials

Discotic hexa-peri-hexabenzocoronene (HBC) derivatives have attracted intensive scientific interest due to their unique optoelectronic properties, which depends, to a large extend, upon the attached functional groups. The presented work covers the synthesis of novel HBC building blocks and new HBC derivatives as functional materials. The traditional preparation of HBC derivatives requires elaborate synthetic techniques and tremendous effort. Especially, more than 10 synthetic steps are usually necessary to approach HBCs with lower symmetries. In order to simplify the synthetic work and reduce the high costs, a novel synthetic strategy involving only four steps was developed based on 2,3,5,6-tetraphenyl-1,4-diiodobenzene intermediates and palladium catalyzed Suzuki cross coupling reactions. In order to introduce various functionalities and expand the diversity of multi-functionalizations, a novel C2v-symmetric dihalo HBC building block 2-47, which contains one iodine and one bromine in para positions, was prepared following the traditional intermolecular [4+2] Diels-Alder reaction route. The outstanding chemical selectivity between iodo and bromo groups in this compound consequently leads to lots of HBC derivatives bearing different functionalities. Directly attached heteroatoms will improve the material properties. According to the application of intramolecular Scholl reaction to a para-dimethoxy HPB, which leads to a meta-dimethoxy HBC, a phenomenon of phenyl group migration was discovered. Thereby, several interesting mechanistic details involving arenium cation intermediates were discussed. With a series of dipole functionalized HBCs, the molecular dynamics of this kind of materials was studied in different phases by DSC, 2D WAXD, solid state NMR and dielectric spectroscopies. High charge carrier mobility is an important parameter for a semiconductive material and depends on the degree of intramolecular order of the discotic molecules in thin films for HBC derivatives. Dipole ��� dipole interaction and hydrogen bonds were respectively introduced in order to achieve highly ordered supramolecular structure. The self-assembly behavior of these materials were investigated both in solution and solid state. Depending upon the different functionalities, these novel materials show either gelating or non-linear optical properties, which consequently broaden their applications as functional materials. In the field of conceivable electronic devices at a molecular level, HBCs hold high promise. Differently functionalized HBCs have been used as active component in the studies of single-molecular CFET and metal-SAMs-metal junctions. The outstanding properties shown in these materials promise their exciting potential applications in molecular devices.

Mixtures of colloidal rods and spheres in bulk and in confinement

When non-adsorbing polymers are added to an isotropic suspension of rod-like colloids, the colloids effectively attract each other via depletion forces. Monte Carlo simulations were performed to study the phase diagram of such rod-polymer mixtures. The colloidal rods were modelled as hard spherocylinders; the polymers were described as spheres of the same diameter as the rods. The polymers may overlap with no energy cost, while overlap of polymers and rods is forbidden. In this thesis the emphasis was on the depletion effects caused by the addition of spheres on the isotropic phase of rod-like particles. Although most of the present experimental studies consider systems close to or beyond the isotropic-nematic transition, the isotropic phase with depletion interactions turns out to be a not less interesting topic. First, the percolation problem was studied in canonical simulations of a system of hard rods and soft spheres, where the amount of depletant was kept low to prevent phase separation of the mixture. The lowering of the percolation threshold seen in experiment is confirmed to be due to the depletion interactions. The local changes in the structure of the fluid of rods, which were measured in the simulations, indicated that the depletion forces enhance local alignment and aggregation of the rods. Then, the phase diagram of isotropic-isotropic demixing of short spherocylinders was calculated using grand canonical ensemble simulations with successive umbrella sampling. Finite size scaling analysis allowed to estimate the location of the critical point. Also, estimates for the interfacial tension between the coexisting isotropic phases and analyses of its power-law behaviour on approach of the critical point are presented. The obtained phase diagram was compared to the predictions of the free volume theory. After an analysis of the bulk, the phase behaviour in confinement was studied. The critical point of gas-liquid demixing is shifted to higher concentrations of rods and smaller concentrations of spheres due to the formation of an orientationally ordered surface film. If the separation between the walls becomes very small, the critical point is shifted back to smaller concentrations of rods because the surface film breaks up. A method to calculate the contact angle of the liquid-gas interface with the wall is introduced and the wetting behaviour on the approach to the critical point is analysed.

Modeling and simulations of some anisotropic soft-matter systems: from biaxial to chiral materials

We have modeled various soft-matter systems with molecular dynamics (MD) simulations. The first topic concerns liquid crystal (LC) biaxial nematic (Nb) phases, that can be possibly used in fast displays. We have investigated the phase organization of biaxial Gay-Berne (GB) mesogens, considering the effects of the orientation, strength and position of a molecular dipole. We have observed that for systems with a central dipole, nematic biaxial phases disappear when increasing dipole strength, while for systems characterized by an offset dipole, the Nb phase is stabilized at very low temperatures. In a second project, in view of their increasing importance as nanomaterials in LC phases, we are developing a DNA coarse-grained (CG) model, in which sugar and phosphate groups are represented with Lennard-Jones spheres, while bases with GB ellipsoids. We have obtained shape, position and orientation parameters for each bead, to best reproduce the atomistic structure of a B-DNA helix. Starting from atomistic simulations results, we have completed a first parametrization of the force field terms, accounting for bonded (bonds, angles and dihedrals) and non-bonded interactions (H-bond and stacking). We are currently validating the model, by investigating stability and melting temperature of various sequences. Finally, in a third project, we aim to explain the mechanism of enantiomeric discrimination due to the presence of a chiral helix of poly(gamma-benzyl L-glutamate) (PBLG), in solution of dimethylformamide (DMF), interacting with chiral or pro-chiral molecules (in our case heptyl butyrate, HEP), after tuning properly an atomistic force field (AMBER). We have observed that DMF and HEP molecules solvate uniformly the PBLG helix, but the pro-chiral solute is on average found closer to the helix with respect to the DMF. The solvent presents a faster isotropic diffusion, twice as HEP, also indicating a stronger interaction of the solute with the helix.

Simulations of charge transport in organic compounds

To aid the design of organic semiconductors, we study the charge transport properties of organic liquid crystals, i.e. hexabenzocoronene and carbazole macrocycle, and single crystals, i.e. rubrene, indolocarbazole and benzothiophene derivatives (BTBT, BBBT). The aim is to find structure-property relationships linking the chemical structure as well as the morphology with the bulk charge carrier mobility of the compounds. To this end, molecular dynamics (MD) simulations are performed yielding realistic equilibrated morphologies. Partial charges and molecular orbitals are calculated based on single molecules in vacuum using quantum chemical methods. The molecular orbitals are then mapped onto the molecular positions and orientations, which allows calculation of the transfer integrals between nearest neighbors using the molecular orbital overlap method. Thus we obtain realistic transfer integral distributions and their autocorrelations. In case of organic crystals the differences between two descriptions of charge transport, namely semi-classical dynamics (SCD) in the small polaron limit and kinetic Monte Carlo (KMC) based on Marcus rates, are studied. The liquid crystals are investigated solely in the hopping limit. To simulate the charge dynamics using KMC, the centers of mass of the molecules are mapped onto lattice sites and the transfer integrals are used to compute the hopping rates. In the small polaron limit, where the electronic wave function is spread over a limited number of neighboring molecules, the Schroedinger equation is solved numerically using a semi-classical approach. The results are compared for the different compounds and methods and, where available, with experimental data. The carbazole macrocycles form columnar structures arranged on a hexagonal lattice with side chains facing inwards, so columns can closely approach each other allowing inter-columnar and thus three-dimensional transport. When taking only intra-columnar transport into account, the mobility is orders of magnitude lower than in the three-dimensional case. BTBT is a promising material for solution-processed organic field-effect transistors. We are able to show that, on the time-scales of charge transport, static disorder due to slow side chain motions is the main factor determining the mobility. The resulting broad transfer integral distributions modify the connectivity of the system but sufficiently many fast percolation paths remain for the charges. Rubrene, indolocarbazole and BBBT are examples of crystals without significant static disorder. The high mobility of rubrene is explained by two main features: first, the shifted cofacial alignment of its molecules, and second, the high center of mass vibrational frequency. In comparsion to SCD, only KMC based on Marcus rates is capable of describing neighbors with low coupling and of taking static disorder into account three-dimensionally. Thus it is the method of choice for crystalline systems dominated by static disorder. However, it is inappropriate for the case of strong coupling and underestimates the mobility of well-ordered crystals. SCD, despite its one-dimensionality, is valuable for crystals with strong coupling and little disorder. It also allows correct treatment of dynamical effects, such as intermolecular vibrations of the molecules. Rate equations are incapable of this, because simulations are performed on static snapshots. We have thus shown strengths and weaknesses of two state of the art models used to study charge transport in organic compounds, partially developed a program to compute and visualize transfer integral distributions and other charge transport properties, and found structure-mobility relations for several promising organic semiconductors.

Morphology and charge transport in conjugated polymers

To assist rational compound design of organic semiconductors, two problems need to be addressed. First, the material morphology has to be known at an atomistic level. Second, with the morphology at hand, an appropriate charge transport model needs to be developed in order to link charge carrier mobility to structure.rnrnThe former can be addressed by generating atomistic morphologies using molecular dynamics simulations. However, the accessible range of time- and length-scales is limited. To overcome these limitations, systematic coarse-graining methods can be used. In the first part of the thesis, the Versatile Object-oriented Toolkit for Coarse-graining Applications is introduced, which provides a platform for the implementation of coarse-graining methods. Tools to perform Boltzmann inversion, iterative Boltzmann inversion, inverse Monte Carlo, and force-matching are available and have been tested on a set of model systems (water, methanol, propane and a single hexane chain). Advantages and problems of each specific method are discussed.rnrnIn partially disordered systems, the second issue is closely connected to constructing appropriate diabatic states between which charge transfer occurs. In the second part of the thesis, the description initially used for small conjugated molecules is extended to conjugated polymers. Here, charge transport is modeled by introducing conjugated segments on which charge carriers are localized. Inter-chain transport is then treated within a high temperature non-adiabatic Marcus theory while an adiabatic rate expression is used for intra-chain transport. The charge dynamics is simulated using the kinetic Monte Carlo method.rnrnThe entire framework is finally employed to establish a relation between the morphology and the charge mobility of the neutral and doped states of polypyrrole, a conjugated polymer. It is shown that for short oligomers, charge carrier mobility is insensitive to the orientational molecular ordering and is determined by the threshold transfer integral which connects percolating clusters of molecules that form interconnected networks. The value of this transfer integral can be related to the radial distribution function. Hence, charge mobility is mainly determined by the local molecular packing and is independent of the global morphology, at least in such a non-crystalline state of a polymer.

Structure formation as studied by EPR spectroscopy: From simple solutions to nanoparticles

In this thesis, three nitroxide based ionic systems were used to investigate structure and dynamics of their respective solutions in mixed solvents by means of electron paramagnetic resonance (EPR) and electron nuclear double resonance (ENDOR) spectroscopy at X- and W-band (9.5 and 94.5 GHz, respectively). rnFirst, the solvation of the inorganic radical Fremy���s salt (K2ON(SO3)2) in isotope substituted binary solvent mixtures (methanol/water) was investigated by means of high-field (W-band) pulse ENDOR spectroscopy and molecular dynamics (MD) simulations. From the analysis of orientation-selective 1H and 2H ENDOR spectra the principal components of the hyperfine coupling (hfc) tensor for chemically different protons (alcoholic methyl vs. exchangeable protons) were obtained. The methyl protons of the organic solvent approach with a mean distance of 3.5 �� perpendicular to the approximate plane spanned by ON(S)2 of the probe molecule. Exchangeable protons were found to be distributed isotropically, approaching closest to Fremy���s salt from the hydrogen-bonded network around the sulfonate groups. The distribution of exchangeable and methyl protons as found in MD simulations is in full agreement with the ENDOR results. The solvation was found to be similar for the studied solvent ratios between 1:2.3 and 2.3:1 and dominated by an interplay of H-bond (electrostatic) interactions and steric considerations with the NO group merely involved into H-bonds.rnFurther, the conformation of spin labeled poly(diallyldimethylammonium chloride) (PDADMAC) solutions in aqueous alcohol (methanol, ethanol, n-propanol, ethylene glycol, glycerol) mixtures in dependence of divalent sodium sulfate was investigated with double electron-electron resonance (DEER) spectroscopy. The DEER data was analyzed using the worm-like chain model which suggests that in organic-water solvent mixtures the polymer backbones are preferentially solvated by the organic solvent. We found a less serve impact on conformational changes due to salt than usually predicted in polyelectrolyte theory which stresses the importance of a delicate balance of hydrophobic and electrostatic interactions, in particular in the presence of organic solvents.rnFinally, the structure and dynamics of miniemulsions and polymerdispersions prepared with anionic surfactants, that were partially replaced by a spin labeled fatty acid in presence and absence of a lanthanide beta-diketonate complex was characterized by CW EPR spectroscopy. Such miniemulsions form multilayers with the surfactant head group bound to the lanthanide ion. Beta-diketonates were formerly used as NMR shift reagents and nowadays find application as luminescent materials in OLEDs and LCDs and as contrast agent in MRT. The embedding of the complex into a polymer matrix results in an easy processable material. It was found that the structure formation takes place in miniemulsion and is preserved during polymerization. For surfactants with carboxyl-head group a higher order of the alkyl chains and less lateral diffusion is found than for sulfat-head groups, suggesting a more uniform and stronger coordination to the metal ion. The stability of these bilayers depends on the temperature and the used surfactant which should be considered for the used polymerization temperature if a maximum output of the structured regions is wished.

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.