New path integral simulation algorithms and their application to creep in the quantum sine-Gordon chain

A path integral simulation algorithm which includes a higher-order Trotter approximation (HOA)is analyzed and compared to an approach which includes the correct quantum mechanical pair interaction (effective Propagator (EPr)). It is found that the HOA algorithmconverges to the quantum limit with increasing Trotter number P as P^{-4}, while the EPr algorithm converges as P^{-2}.The convergence rate of the HOA algorithm is analyzed for various physical systemssuch as a harmonic chain,a particle in a double-well potential, gaseous argon, gaseous helium and crystalline argon. A new expression for the estimator for the pair correlation function in the HOA algorithm is derived. A new path integral algorithm, the hybrid algorithm, is developed.It combines an exact treatment of the quadratic part of the Hamiltonian and thehigher-order Trotter expansion techniques.For the discrete quantum sine-Gordon chain (DQSGC), it is shown that this algorithm works more efficiently than all other improved path integral algorithms discussed in this work. The new simulation techniques developed in this work allow the analysis of theDQSGC and disordered model systems in the highly quantum mechanical regime using path integral molecular dynamics (PIMD)and adiabatic centroid path integral molecular dynamics (ACPIMD).The ground state phonon dispersion relation is calculated for the DQSGC by the ACPIMD method.It is found that the excitation gap at zero wave vector is reduced by quantum fluctuations. Two different phases exist: One phase with a finite excitation gap at zero wave vector, and a gapless phase where the excitation gap vanishes.The reaction of the DQSGC to an external driving force is analyzed at T=0.In the gapless phase the system creeps if a small force is applied, and in the phase with a gap the system is pinned. At a critical force, the systems undergo a depinning transition in both phases and flow is induced. The analysis of the DQSGC is extended to models with disordered substrate potentials. Three different cases are analyzed: Disordered substrate potentials with roughness exponent H=0, H=1/2,and a model with disordered bond length. For all models, the ground state phonon dispersion relation is calculated.

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.

Molekulardynamik-Computersimulation einer amorph-kristallinen SiO 2 Grenzschicht

In dieser Arbeit werden Molekulardynamik-Computersimulationen zur Untersuchung der statischen und dynamischen Eigenschaften einer amorph/kristallinen Siliziumdioxid(SiO2)-Grenzschicht durchgefuehrt.Die Grenzflaeche wird von der [100]-Ebene des beta-Kristobalit-Kristalls und der fluessigen SiO2-Phase gebildet und in einem Temperaturbereich zwischen 2900K und 3100K im Zustand eines metastabilen Gleichgewichts untersucht. Als Modellpotential zur Beschreibung der mikroskopischen Wechselwirkungen zwischen den Teilchen wird ein einfaches Paarpotential aus der Literatur verwendet, das sowohl die Struktur der kristallinen Phase als auch die der fluessigen Phase gut reproduziert. Bezogen auf die Dichte und die potentielle Energie der Teilchen erstreckt sich der Uebergang von der fluessigen in die kristalline Phase ueber 3-5 Atomlagen. Ein Layering-Effekt der Dichte in der fluessigen Phase in der Naehe der Grenzschicht wird nicht beobachtet. Der Einfluss der Grenzschicht auf statische Groessen, welche das System auf einer mittelreichweitigen Laengenskala beschreiben (z. B. Koordinationszahlverteilung und Ringverteilung) reicht im Vergleich dazu weiter in die fluessige Phase hinein und manifestiert sich in Defektstrukturen, wie z. B. der Erhoehung der Wahrscheinlichkeit fuer das Auftreten von 5-fach koordiniertem Silizium und der vermehrten Bildung von 2er-Ringen in der Fluessigkeit. Dies beguenstigt das Aufbrechen und Umklappen von Si-O-Bindungen und fuehrt zu einer Beschleunigung der Dynamik und einer Erhoehung der Diffusionsgeschwindigkeit in der Fluessigkeit. Im weiteren wird die Hochfrequenzdynamik der reinen SiO2-Fluessigkeit untersucht. Dazu berechnen wir die vibratorische Zustandsdichte in harmonischer Naeherung aus der inhaerenten Struktur. Wir finden einen stark ausgepraegten Peak bei einer Frequenz von 0.6 THz. Dieser Peak kann der niederenergetischsten transversalen akustischen Mode zugeordnet werden, die auch als Scherschwingung des Systems direkt sichtbar ist.

Theoretical investigation of the interaction of a polymer film with a nanoparticle

The fundamental aim in our investigation of the interaction of a polymer film with a nanoparticle is the extraction of information on the dynamics of the liquid using a single tracking particle. In this work two theoretical methods were used: one passive, where the motion of the particle measures the dynamics of the liquid, one active, where perturbations in the system are introduced through the particle. In the first part of this investigation a thin polymeric film on a substrate is studied using molecular dynamics simulations. The polymer is modeled via a 'bead spring' model. The particle is spheric and non structured and is able to interact with the monomers via a Lennard Jones potential. The system is micro-canonical and simulations were performed for average temperatures between the glass transition temperature of the film and its dewetting temperature. It is shown that the stability of the nanoparticle on the polymer film in the absence of gravity depends strongly on the form of the chosen interaction potential between nanoparticle and polymer. The relative position of the tracking particle to the liquid vapor interface of the polymer film shows the glass transition of the latter. The velocity correlation function and the mean square displacement of the particle has shown that it is caged when the temperature is close to the glass transition temperature. The analysis of the dynamics at long times shows the coupling of the nanoparticle to the center of mass of the polymer chains. The use of the Stokes-Einstein formula, which relates the diffusion coefficient to the viscosity, permits to use the nanoparticle as a probe for the determination of the bulk viscosity of the melt, the so called 'microrheology'. It is shown that for low frequencies the result obtained using microrheology coincides with the results of the Rouse model applied to the polymer dynamics. In the second part of this investigation the equations of Linear Hydrodynamics are solved for a nanoparticle oscillating above the film. It is shown that compressible liquids have mechanical response to external perturbations induced with the nanoparticle. These solutions show strong velocity and pressure profiles of the liquid near the interface, as well as a mechanical response of the liquid-vapor interface. The results obtained with this calculations can be employed for the interpretation of experimental results of non contact AFM microscopy

Molekulardynamiksimulationen zur Untersuchung des Mischalkali-Effektes in silikatischen Gläsern

In der vorliegenden Arbeit wird mittels Molekulardynamik(MD)-Computersimulationen die Dynamik von verschiedenen Alkalisilikaten in der Schmelze und im Glas untersucht. Es ist bekannt, daß diese Systeme ionenleitend sind, was auf eine hohe Mobilität der Alkaliionen im Vergleich zu den glasbildenden Komponenten Si und O zurückzuführen ist. Im Mittelpunkt des Interesses steht der sog. Mischalkalieffekt (MAE), der in ternären Mischungen aus Siliziumdioxid mit zwei Alkalioxiden auftritt. Gegenüber Mischungen mit nur einer Alkaliionensorte weisen letztere Systeme eine signifikante Verlangsamung der Alkaliionendiffusion auf. Zunächst werden zwei binäre Alkalisilikate simuliert, nämlich Lithiumdisilikat (LS2) und Kaliumdisilikat (KS2). Die Simulationen zeigen, daß der Ursprung der hohen Mobilität der Alkaliionen in der Struktur begründet ist. KS2 und LS2 weisen auf intermediären Längenskalen Ordnung auf, die in partiellen statischen Strukturfaktoren durch Prepeaks reflektiert ist. Die den Prepeaks zugrundeliegende Struktur erklärt sich durch perkolierende Netzwerke aus alkalioxidreichen Kanälen, die als Diffusionskanäle für die mobilen Alkaliionen fungieren. In diesen Kanälen bewegen sich die Ionen mittels Sprüngen (Hopping) zwischen ausgezeichneten Plätzen. In der Simulation beobachtet man für die hohen Temperaturen (4000K>=1500K) eine ähnliche Aktivierungsenergie wie im Experiment. Im Experiment findet allerdings unterhalb von ca.1200K ein Crossover in ein Arrheniusverhalten mit höherer Aktivierungsenergie statt, welches von der Simulation nicht nachvollzogen wird. Das kann mit der in der Simulation nicht im Gleichgewicht befindlichen Si-O-Matrix erklärt werden, bei der Alterungseffekte beobachtet werden. Am stärksten ist der MAE für eine Alkalikomponente, wenn deren Konzentrationsanteil in einem ternären Mischalkalisystem gegen 0 geht. Daher wird ein LS2-System untersucht, in dem ein Li-Ion gegen ein K-Ion getauscht wird. Der Einfluß des K-Ions ist sowohl lokal in den charakteristischen Abständen zu den ersten nächsten Nachbarn (NN) zu sehen, als auch in der ortsaufgelösten Koordinationszahlverteilung bis zu Längenskalen von ca. 8,5 Angstrom. Die Untersuchung der Dynamik des eingesetzten K-Ions zeigt, daß die Sprungwahrscheinlichkeit nicht mit der Lokalisierung, einem Maß für die Bewegung eines Teilchens um seine Ruheposition, korreliert ist, aber daß eine chemische Umgebung mit wenig Li- und vielen O-NN oder vielen Li- und wenig O-NN ein Sprungereignis begünstigt. Zuletzt wird ein ternäres Alkalisilikat (LKS2) untersucht, dessen Struktur alle charakteristischen Längenskalen von LS2 und KS2 aufweist. Es stellt sich also eine komplexe Struktur mit zwei perkolierenden Subnetzwerken für Alkaliionen ein. Die Untersuchung der Dynamik zeigt eine geringe Wahrscheinlichkeit dafür auf, daß Ionen in ein Subnetzwerk andersnamiger Ionen springen. Auch kann gezeigt werden, daß das Modellpotential den MAE reproduzieren kann, daß also die Diffusionskonstanten in LKS2 bei bis zu einer Größenordnung langsamer sind als in KS2 bzw. LS2. Der beobachtete Effekt stellt sich zudem vom funktionalen Verlauf her so dar, wie er beim MAE erwartet wird. Es wurde auch festgestellt, daß trotz der zeitlichen Verzögerung in den dynamischen Größen die Anzahl der Sprünge pro Zeit nicht geringer ist und daß für niedrige Temperaturen (d.h.im Glas) Sprünge auf den Nachbarplatz mit anschließendem Rücksprung auf die vorherige Position deutlich wahrscheinlicher sind als bei hohen Temperaturen (also in der Schmelze). Die vorliegenden Resultate geben Aufschluß über die Details der Mechanismen mikroskopischer Ionenleitung in binären und ternären Alkalisilikaten sowie dem MAE.

Einfluß von Topologie und Verscherung auf Dynamik und mechanische Parameter bei Homopolymeren

Die FT-Rheologie wird zur Unterscheidung verschiedener Kamm-Topologien in Polystyrollösungen und –schmelzen angewendet. Die Polystyrole werden in Abhängigkeit der Deborahzahl De unter LAOS-Bedingungen vermessen. Die Meßergebnisse zeigen, daß der Schritt von wohldefinierten Systemen (lineare Ketten, Sterne) zu solchen mit statistischer Verteilung wie in Kämmen zu großen Veränderungen sowohl im linearen als auch im nichtlinearen Bereich der rheologischen Messungen führt. Sowohl die Masterkurven als auch die Intenstiäten I3/1 und Phasen Phi3 der Nichtlinearitäten der einzelnen Proben weisen jeweils deutliche Unterschiede untereinander auf. Diese sind durch die bisherigen Ergebnisse noch nicht vollständig mit topologischen Merkmalen in Verbindung zu bringen. Die Messungen wurden mit dem von McLeish eingeführten Pom-pom Modell und daraus weiterentwickelten double convected-Pom-pom Modell (DCPP) simuliert und lieferten gute Übereinstimmung sowie auch Vorhersagen über den experimentell nicht mehr zugänglichen Bereich. Zur Untersuchung des Einflusses von mechanischer Scherung auf die lokale, molekulare Dynamik wird das LAOS-Experiment in situ mit dielektrischer Spektroskopie kombiniert. Dazu wurde eine Apparatur entwickelt, die das hochsensitive ARES-Rheometer mit dem hochauflösenden dielektrischen ALPHA-Analyzer verbindet. Mit dieser Apparatur wurde das Typ-A Polymer 1,4-cis-Polyisopren, mit einem Dipolmoment entlang des Rückgrats, bei oszillatorischer Scherung unter gleichzeitiger Aufnahme eines dielektrischen Spektrums vermessen. Es konnte gezeigt werden, daß die oszillatorische Verscherung weder die charakteristische Relaxationszeit noch die Form des Normal Mode Peaks beeinflußt, wohl aber die dielektrische Stärke Delta epsilon. Diese entspricht der Fläche unter dem e“-Peak und kann mit einer Debye- und einer Cole/Davidson-Funktion angepasst werden. Die Abnahme der dielektrischen Stärke mit zunehmender Scheramplitude kann mit der Orientierungsverteilung der End-zu-End-Vektoren in der Probe erklärt werden.

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.