Crystal phases and glassy dynamics in monodisperse hard ellipsoids

We have performed Monte Carlo and molecular dynamics simulations of suspensions of monodisperse, hard ellipsoids of revolution. Hard-particle models play a key role in statistical mechanics. They are conceptually and computationally simple, and they offer insight into systems in which particle shape is important, including atomic, molecular, colloidal, and granular systems. In the high density phase diagram of prolate hard ellipsoids we have found a new crystal, which is more stable than the stretched FCC structure proposed previously . The new phase, SM2, has a simple monoclinic unit cell containing a basis of two ellipsoids with unequal orientations. The angle of inclination is very soft for length-to-width (aspect) ratio l/w=3, while the other angles are not. A symmetric state of the unit cell exists, related to the densest-known packings of ellipsoids; it is not always the stable one. Our results remove the stretched FCC structure for aspect ratio l/w=3 from the phase diagram of hard, uni-axial ellipsoids. We provide evidence that this holds between aspect ratios 3 and 6, and possibly beyond. Finally, ellipsoids in SM2 at l/w=1.55 exhibit end-over-end flipping, warranting studies of the cross-over to where this dynamics is not possible. Secondly, we studied the dynamics of nearly spherical ellipsoids. In equilibrium, they show a first-order transition from an isotropic phase to a rotator phase, where positions are crystalline but orientations are free. When over-compressing the isotropic phase into the rotator regime, we observed super-Arrhenius slowing down of diffusion and relaxation, and signatures of the cage effect. These features of glassy dynamics are sufficiently strong that asymptotic scaling laws of the Mode-Coupling Theory of the glass transition (MCT) could be tested, and were found to apply. We found strong coupling of positional and orientational degrees of freedom, leading to a common value for the MCT glass-transition volume fraction. Flipping modes were not slowed down significantly. We demonstrated that the results are independent of simulation method, as predicted by MCT. Further, we determined that even intra-cage motion is cooperative. We confirmed the presence of dynamical heterogeneities associated with the cage effect. The transit between cages was seen to occur on short time scales, compared to the time spent in cages; but the transit was shown not to involve displacements distinguishable in character from intra-cage motion. The presence of glassy dynamics was predicted by molecular MCT (MMCT). However, as MMCT disregards crystallization, a test by simulation was required. Glassy dynamics is unusual in monodisperse systems. Crystallization typically intervenes unless polydispersity, network-forming bonds or other asymmetries are introduced. We argue that particle anisometry acts as a sufficient source of disorder to prevent crystallization. This sheds new light on the question of which ingredients are required for glass formation.

Novel algorithms for electronic structure based molecular dynamics with linear system-size scaling

Die vorliegende Arbeit behandelt die Entwicklung und Verbesserung von linear skalierenden Algorithmen für Elektronenstruktur basierte Molekulardynamik. Molekulardynamik ist eine Methode zur Computersimulation des komplexen Zusammenspiels zwischen Atomen und Molekülen bei endlicher Temperatur. Ein entscheidender Vorteil dieser Methode ist ihre hohe Genauigkeit und Vorhersagekraft. Allerdings verhindert der Rechenaufwand, welcher grundsätzlich kubisch mit der Anzahl der Atome skaliert, die Anwendung auf große Systeme und lange Zeitskalen. Ausgehend von einem neuen Formalismus, basierend auf dem großkanonischen Potential und einer Faktorisierung der Dichtematrix, wird die Diagonalisierung der entsprechenden Hamiltonmatrix vermieden. Dieser nutzt aus, dass die Hamilton- und die Dichtematrix aufgrund von Lokalisierung dünn besetzt sind. Das reduziert den Rechenaufwand so, dass er linear mit der Systemgröße skaliert. Um seine Effizienz zu demonstrieren, wird der daraus entstehende Algorithmus auf ein System mit flüssigem Methan angewandt, das extremem Druck (etwa 100 GPa) und extremer Temperatur (2000 - 8000 K) ausgesetzt ist. In der Simulation dissoziiert Methan bei Temperaturen oberhalb von 4000 K. Die Bildung von sp²-gebundenem polymerischen Kohlenstoff wird beobachtet. Die Simulationen liefern keinen Hinweis auf die Entstehung von Diamant und wirken sich daher auf die bisherigen Planetenmodelle von Neptun und Uranus aus. Da das Umgehen der Diagonalisierung der Hamiltonmatrix die Inversion von Matrizen mit sich bringt, wird zusätzlich das Problem behandelt, eine (inverse) p-te Wurzel einer gegebenen Matrix zu berechnen. Dies resultiert in einer neuen Formel für symmetrisch positiv definite Matrizen. Sie verallgemeinert die Newton-Schulz Iteration, Altmans Formel für beschränkte und nicht singuläre Operatoren und Newtons Methode zur Berechnung von Nullstellen von Funktionen. Der Nachweis wird erbracht, dass die Konvergenzordnung immer mindestens quadratisch ist und adaptives Anpassen eines Parameters q in allen Fällen zu besseren Ergebnissen führt.

Structure, organization and dynamics of functional supramolecular materials studied by solid-state NMR

Functional materials have great importance due to their many important applications. The characterization of supramolecular architectures which are held together by non-covalent interactions is of most importance to understand their properties. Solid-state NMR methods have recently been proven to be able to unravel such structure-property relations with the help of fast magic-angle spinning and advanced pulse sequences. The aim of the current work is to understand the structure and dynamics of functional supramolecular materials which are potentially important for fuel-cell (proton conducting membrane materials) and solar-cell or plastic-electronic applications (photo-reactive aromatic materials). In particular, hydrogen-bonding networks, local proton mobility, molecular packing arrangements, and local dynamics will be studied by the use of advanced solid-state NMR methods. The first class of materials studied in this work is proton conducting polymers which also form hydrogen-bonding network. Different materials, which are prepared for high 1H conduction by different approaches are studied: PAA-P4VP, PVPA-ABPBI, Tz5Si, and Triazole-functional systems. The materials are examples of the following major groups; - Homopolymers with specific functional groups (Triazole functional polysiloxanes). - Acid-base polymer blends approach (PAA-P4VP, PVPA-ABPBI). - Acid-base copolymer approach (Triazole-PVPA). - Acid doped polymers (Triazole functional polymer doped with H3PO4). Perylenebisimide (PBI) derivatives, a second type of important functional supramolecular materials with potent applications in plastic electronics, were also investigated by means of solid-state NMR. The preparation of conducting nanoscopic fibers based on the self-assembling functional units is an appealing aim as they may be incorporated in molecular electronic devices. In this category, perylene derivatives have attracted great attention due to their high charge carrier mobility. A detailed knowledge about their supramolecular structure and molecular dynamics is crucial for the understanding of their electronic properties. The aim is to understand the structure, dynamics and packing arrangements which lead to high electron conductivity in PBI derivatives.

A molecular dynamics study of the influence of chain branching on the properties of polymer systems

Die vorliegende Arbeit beschäftigt sich mit dem Einfluß von Kettenverzweigungen unterschiedlicher Topologien auf die statischen Eigenschaften von Polymeren. Diese Untersuchungen werden mit Hilfe von Monte-Carlo- und Molekular-Dynamik-Simulationen durchgeführt.Zunächst werden einige theoretische Konzepte und Modelle eingeführt, welche die Beschreibung von Polymerketten auf mesoskopischen Längenskalen gestatten. Es werden wichtige Bestimmungsgrößen eingeführt und erläutert, welche zur quantitativen Charakterisierung von Verzweigungsstrukturen bei Polymeren geeignet sind. Es wird ebenso auf die verwendeten Optimierungstechniken eingegangen, die bei der Implementierung des Computerprogrammes Verwendung fanden. Untersucht werden neben linearen Polymerketten unterschiedliche Topolgien -Sternpolymere mit variabler Armzahl, Übergang von Sternpolymeren zu linearen Polymeren, Ketten mit variabler Zahl von Seitenketten, reguläre Dendrimere und hyperverzweigte Strukturen - in Abhängigkeit von der Lösungsmittelqualität. Es wird zunächst eine gründliche Analyse des verwendeten Simulationsmodells an sehr langen linearen Einzelketten vorgenommen. Die Skalierungseigenschaften der linearen Ketten werden untersucht in dem gesamten Lösungsmittelbereich vom guten Lösungsmittel bis hin zu weitgehend kollabierten Ketten im schlechten Lösungsmittel. Ein wichtiges Ergebnis dieser Arbeit ist die Bestätigung der Korrekturen zum Skalenverhalten des hydrodynamischen Radius Rh. Dieses Ergebnis war möglich aufgrund der großen gewählten Kettenlängen und der hohen Qualität der erhaltenen Daten in dieser Arbeit, insbesondere bei den linearen ketten, und es steht im Widerspruch zu vielen bisherigen Simulations-Studien und experimentellen Arbeiten. Diese Korrekturen zum Skalenverhalten wurden nicht nur für die linearen Ketten, sondern auch für Sternpolymere mit unterchiedlicher Armzahl gezeigt. Für lineare Ketten wird der Einfluß von Polydispersität untersucht.Es wird gezeigt, daß eine eindeutige Abbildung von Längenskalen zwischen Simulationsmodell und Experiment nicht möglich ist, da die zu diesem Zweck verwendete dimensionslose Größe eine zu schwache Abhängigkeit von der Polymerisation der Ketten besitzt. Ein Vergleich von Simulationsdaten mit industriellem Low-Density-Polyäthylen(LDPE) zeigt, daß LDPE in Form von stark verzweigten Ketten vorliegt.Für reguläre Dendrimere konnte ein hochgradiges Zurückfalten der Arme in die innere Kernregion nachgewiesen werden.

Molecular dynamics simulations of silicate and borate glasses and melts : structure, diffusion dynamics and vibrational properties

Molecular dynamics simulations of silicate and borate glasses and melts: Structure, diffusion dynamics and vibrational properties. In this work computer simulations of the model glass formers SiO2 and B2O3 are presented, using the techniques of classical molecular dynamics (MD) simulations and quantum mechanical calculations, based on density functional theory (DFT). The latter limits the system size to about 100−200 atoms. SiO2 and B2O3 are the two most important network formers for industrial applications of oxide glasses. Glass samples are generated by means of a quench from the melt with classical MD simulations and a subsequent structural relaxation with DFT forces. In addition, full ab initio quenches are carried out with a significantly faster cooling rate. In principle, the structural properties are in good agreement with experimental results from neutron and X-ray scattering, in all cases. A special focus is on the study of vibrational properties, as they give access to low-temperature thermodynamic properties. The vibrational spectra are calculated by the so-called ”frozen phonon” method. In all cases, the DFT curves show an acceptable agreement with experimental results of inelastic neutron scattering. In case of the model glass former B2O3, a new classical interaction potential is parametrized, based on the liquid trajectory of an ab initio MD simulation at 2300 K. In this course, a structural fitting routine is used. The inclusion of 3-body angular interactions leads to a significantly improved agreement of the liquid properties of the classical MD and ab initio MD simulations. However, the generated glass structures, in all cases, show a significantly lower fraction of 3-membered planar boroxol rings as predicted by experimental results (f=60%-80%). The largest boroxol ring fraction of f=15±5% is observed in the full ab initio quenches from 2300 K. In case of SiO2, the glass structures after the quantum mechanical relaxation are the basis for calculations of the linear thermal expansion coefficient αL(T), employing the quasi-harmonic approximation. The striking observation is a change change of sign of αL(T) going along with a temperature range of negative αL(T) at low temperatures, which is in good agreement with experimental results.

Influence of chain stiffness on structure and dynamics of polymers in the melt

Die vorliegende Arbeit beschaeftigt sich mit der Untersuchung vonPolymeren mit intrinsischer Steifigkeit. Es werden vor allem lokale statische unddynamische Eigenschaften anhand zweier verschiedener Simulationsmodellebetrachtet: Ein generisches Polymermodell, bei dem nur dieSteifigkeit als ein das spezifische Polymer charakterisierenden Parametereingeht und ein atomistisches Modell fuer trans-Polyisopren. Mit Hilfe des ersten Modells koennen Statik und Dynamik wurmartiger Kettenbeobachtet werden. Das Blob-Konzept ist eine angemessene statischeBeschreibung. Lokale Orientierungen haengen schwach von derSteifigkeit ab. Das Reptationsmodell kann die beobachtete Dynamik fuer lange Kettennicht mehr angemessen beschreiben. Lange Ketten bewegen sich, als obsie in Roehren gezwaengt waeren; jedoch ist die Bewegung starkabhaengig von der Steifigkeit. Fuer Ketten dieser Art konntequalitativ das Verhalten reproduziert werden, das in NMR-Experimentenbeobachtet wird. Eine Verhakungslaenge laesst sich fuer solche Kettenkaum mehr definieren. Dynamische Strukturfunktionen und insbesonderedie direkte Visualisierung der Ketten verdeutlichen die effektiv aufeine Roehre beschraenkte Bewegung. Das atomistische Polyisoprenmodell wurde mit verschiedenen Experimenten,verglichen. In den Simulationen bei konnten qualitativ undsemiquantitativ experimentelle Ergebnisse reproduziert werden. Zuletzt wurden die Laengen- und Zeitskalen der beiden Modelleerfolgreich aufeinander abgebildet.

Structure and dynamics of supramolecular assemblies studied by advanced solid-state NMR spectroscopy

Ziel der vorliegenden Arbeit ist die Aufklärung von Struktur und Dynamik komplexer supramolekularer Systeme mittels Festkörper NMR Spektroskopie. Die Untersuchung von pi-pi Wechselwirkungen, welche einen entscheidenden Einfluss auf die strukturellen und dynamischen Eigenschaften supra- molekularer Systeme haben, hilft dabei, die Selbst- organisationsprozesse dieser komplexen Materialien besser zu verstehen. Mit dipolaren 1H-1H and 1H-13C Wiedereinkopplungs NMR Methoden unter schnellem MAS können sowohl 1H chemische Verschiebungen als auch dipolare 1H-1H und 1H-13C Kopplungen untersucht werden, ohne dass eine Isotopenmarkierung erforderlich ist. So erhält man detaillierte Informationen über die Struktur und die Beweglichkeit einzelner Molekül- segmente. In Verbindung mit sogenannten nucleus independent chemical shift (NICS) maps (berechnet mit ab-initio Methoden) lassen sich Abstände von Protonen relativ zu pi-Elektronensystemen bestimmen und so Strukturvorschläge ableiten. Mit Hilfe von homo- und heteronuklearen dipolaren Rotationsseitenbandenmustern könnenaußerdem Ordnungs- parameter für verschiedene Molekülsegmente bestimmt werden. Die auf diese Weise gewonnenen Informationen über die strukturellen und dynamischen Eigenschaften supramolekularer Systeme tragen dazu bei, strukturbestimmende Molekül- einheiten und Hauptordnungsphänomene zu identifizieren sowie lokale Wechselwirkungen zu quantifizieren, um so den Vorgang der Selbstorganisation besser zu verstehen.

Structure and dynamics of ionic liquid,electrode interfaces

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

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.

Supramolecular order and dynamics of functional materials studied by solid state NMR

The goal of this thesis was the investigation of the structure, conformation, supramolecular order and molecular dynamics of different classes of functional materials (phthalocyanine, perylene and hexa-peri-hexabenzocoronene derivatives and mixtures of those), all having planar aromatic cores modified with various types of alkyl chains. The planar aromatic systems are known to stack in the solid and the liquid-crystalline state due to p-p interactions forming columnar superstructures with high one-dimensional charge carrier mobility and potential application in photovoltaic devices. The different functionalities attached to the aromatic cores significantly influence the behavior of these systems allowing the experimentalists to modify the structures to fine-tune the desired thermotropic properties or charge carrier mobility. The aim of the presented studies was to understand the interplay between the driving forces causing self-assembly by relating the structural and dynamic information about the investigated systems. The supramolecular organization is investigated by applying 1H solid state NMR recoupling techniques. The results are related with DSC and X-ray scattering data. Detailed information about the site-specific molecular dynamics is gained by recording spinning sideband patterns using 1H-1H and 13C-1H solid state NMR recoupling techniques. The determined dipole-dipole coupling constants are then related with the coupling constants of the respective rigid pairs, thus providing local dynamic order parameters for the respective moieties. The investigations presented reveal that in the crystalline state the preferred arrangement in the columnar stack of discotic molecules modified with alkyl chains is tilted. This leads to characteristic differences in the 1H chemical shifts of otherwise chemically equivalent protons. Introducing branches and increasing the length of the alkyl chains results in lower mesophase transitions and disordered columnar stacks. In the liquid-crystalline state some of the discs lose the tilted orientation, others do not, but all start a rapid rotation about the columnar axis.

Spectroscopic properties from hybrid QM, MM molecular dynamics simulation

The central aim of this thesis work is the application and further development of a hybrid quantum mechanical/molecular mechanics (QM/MM) based approach to compute spectroscopic properties of molecules in complex chemical environments from electronic structure theory. In the framework of this thesis, an existing density functional theory implementation of the QM/MM approach is first used to calculate the nuclear magnetic resonance (NMR) solvent shifts of an adenine molecule in aqueous solution. The findings show that the aqueous solvation with its strongly fluctuating hydrogen bond network leads to specific changes in the NMR resonance lines. Besides the absolute values, also the ordering of the NMR lines changes under the influence of the solvating water molecules. Without the QM/MM scheme, a quantum chemical calculation could have led to an incorrect assignment of these lines. The second part of this thesis describes a methodological improvement of the QM/MM method that is designed for cases in which a covalent chemical bond crosses the QM/MM boundary. The development consists in an automatized protocol to optimize a so-called capping potential that saturates the electronic subsystem in the QM region. The optimization scheme is capable of tuning the parameters in such a way that the deviations of the electronic orbitals between the regular and the truncated (and "capped") molecule are minimized. This in turn results in a considerable improvement of the structural and spectroscopic parameters when computed with the new optimized capping potential within the QM/MM technique. This optimization scheme is applied and benchmarked on the example of truncated carbon-carbon bonds in a set of small test molecules. It turns out that the optimized capping potentials yield an excellent agreement of NMR chemical shifts and protonation energies with respect to the corresponding full molecules. These results are very promising, so that the application to larger biological complexes will significantly improve the reliability of the prediction of the related spectroscopic properties.

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.

Assessing the functional structure of molecular transporters by EPR spectroscopy

In this thesis, the self-assembled functional structure of a broad range of amphiphilic molecular transporters is studied. By employing paramagnetic probe molecules and ions, continuous-wave and pulse electron paramagnetic resonance spectroscopy reveal information about the local structure of these materials from the perspective of incorporated guest molecules. First, the transport function of human serum albumin for fatty acids is in the focus. As suggested by the crystal structure, the anchor points for the fatty acids are distributed asymmetrically in the protein. In contrast to the crystallographic findings, a remarkably symmetric entry point distribution of the fatty acid binding channels is found, which may facilitate the uptake and release of the guest molecules. Further, the metal binding of 1,2,3-triazole modified star-shaped cholic acid oligomers is studied. These biomimetic molecules are able to include and transport molecules in solvents of different polarity. A pre-arrangement of the triazole groups induces a strong chelate-like binding and close contact between guest molecule and metal ion. In absence of a preordering, each triazole moiety acts as a single entity and the binding affinity for metal ions is strongly decreased. Hydrogels based on N-isopropylacrylamide phase separate from water above a certain temperature. The macroscopic thermal collapse of these hydrogels is utilized as a tool for dynamic nuclear polarization. It is shown that a radical-free hyperpolarized solution can be achieved with a spin-labeled gel as separable matrix. On the nanoscale, these hydrogels form static heterogeneities in both structure and function. Collapsed regions protect the spin probes from a chemical decay while open, water-swollen regions act as catalytic centers. Similarly, thermoresponsive dendronized polymers form structural heterogeneities, which are, however, highly dynamic. At the critical temperature, they trigger the aggregation of the polymer into mesoglobules. The dehydration of these aggregates is a molecularly controlled non-equilibrium process that is facilitated by a hydrophobic dendritic core. Further, a slow heating rate results in a kinetically entrapped non-equilibrium state due to the formation of an impermeable dense polymeric layer at the periphery of the mesoglobule.

Characterization of structure and dynamics of submicrometric systems via fluorescence correlation spectroscopy

The presented thesis revolves around the study of thermally-responsive PNIPAAm-based hydrogels in water/based environments, as studied by Fluorescence Correlation Spectroscopy (FCS).rnThe goal of the project was the engineering of PNIPAAm gels into biosensors. Specifically, a gamma of such gels were both investigated concerning their dynamics and structure at the nanometer scale, and their performance in retaining bound bodies upon thermal collapse (which PNIPAAm undergoes upon heating above 32 ºC).rnFCS’s requirements, as a technique, match the limitations imposed by the system. Namely, the need to intimately probe a system in a solvent, which was also fragile and easy to alter. FCS, on the other hand, both requires a fluid environment to work, and is based on the observation of diffusion of fluorescents at nanomolar concentrations. FCS was applied to probe the hydrogels on the nanometer size with minimal invasivity.rnVariables in the gels were addressed in the project including crosslinking degree; structural changes during thermal collapse; behavior in different buffers; the possibility of decreasing the degree of inhomogeneity; behavior of differently sized probes; and the effectiveness of antibody functionalization upon thermal collapse.rnThe evidenced results included the heightening of structural inhomogeneities during thermal collapse and under different buffer conditions; the use of annealing to decrease the inhomogeneity degree; the use of differently sized probes to address different length scale of the gel; and the successful functionalization before and after collapse.rnThe thesis also addresses two side projects, also carried forward via FCS. One, diffusion in inverse opals, produced a predictive simulation model for diffusion of bodies in confined systems as dependent on the bodies’ size versus the characteristic sizes of the system. The other was the observation of interaction of bodies of opposite charge in a water solution, resulting in a phenomenological theory and an evaluation method for both the average residence time of the different bodies together, and their attachment likelihood.