Ionic liquids in a multiscale study

Trotz des hohen Interesse an Ionischen Flüssigkeiten wird das zielgerichtete Design und die Anwendung Ionischer Flüssigkeiten durch fehlendes grundlegendes Verständnis erschwert. Deshalb wurde die Balance der molekularen Wechselwirkungen in Ionischen Flüssigkeiten studiert, um die Eigenschaften dieser zu verstehen und die Kraftfeldentwicklung im Rahmen des Multiskalenansatzes zu systematisieren. Es wurden reine Imidazolium-basierte Ionische Flüssigkeiten, Mischungen mit kleinen Molekülen und eine protische Ionische Flüssigkeit mit ab-initio-Methoden, hauptsächlich Car-Parrinello-Molekulardynamik, untersucht. Weiterhin wurden Eigenschaften der Flüssigphase mit denen von Ionenpaaren verglichen.rnIm Fokus standen die molekularen elektrostatischen Eigenschaften und es wurde gezeigt, dass Coulomb-Wechselwirkungen zu einzigartigen Charakteristika führten. So waren die Ionen-Nettoladungen stets reduziert, die molekularen Dipolmomentverteilungen sehr breit, elektronische Polarisation war entscheidend. Die elektrostatischen Eigenschaften waren allgemein lokal auf molekularen Größen- und Zeitskalen und hingen stark von Phasenzustand und Zusammensetzung ab. Für andere molekulare Eigenschaften, wie der Neigung zu dispersiven Kontakten oder Wasserstoffbrücken, wurde gezeigt, dass sie einen entscheidenden Einfluss auf die Feinstruktur Ionischer Flüssigkeiten hatten. Das Gleichgewicht der Wechselwirkungen zeigte sich auch in Leistungsspektren, die sich aus den ab-initio-Molekulardynamiksimulationen ergaben. Diese boten einen neuen Weg für den Vergleich zum Experiment und für einen Einblick in die schnelle Dynamik Ionischer Flüssigkeiten.

Dynamik von Polymerschmelzen und Quellverhalten ungeordneter Netzwerke

Deutsch:Mit Hilfe eines parallelen Molekulardynamik-Programmswurden einfach Modelle von Homopolymerschmelzen simuliert.Langkettige Schmelzen zeigten eine sehr gute Übereinstimmungmit den Vorhersagendes Reptationsmodells. Die intermediären Reptationsbereichemit den vorhergesagtenExponenten konnten wesentlich klarer als bisher verifiziertwerden. Es stellteVerschiedene, gebräuchliche Analyse-Methoden führten jedochzuunterschiedlichen Aussagen für die Verhakungslänge. Fürkurze Kettenbzw. kurze Unterketten wurden in dichten Schmelzendie Abweichungen vom Rouse-Modell aufgezeigt. DieseAbweichungenkönnen als Korrelationslocheffekt interpretiert werden undsteht in teilweiser Übereinstimmung zu den Vorhersagenrenormierter Rouse-Modelle. Aus den Schmelzen wurden Netzwerke in einem speziellenZufallsvernetzungsprozeßhergestellt, der die Bildung von Defektstrukturenunterbindet. Diesewurden bzgl. ihres Quell- und Deformationsverhaltenuntersucht. Der maximaleQuellgrad eines Netzwerkes war bereits bei verhaltnismäßigkurzen Kettenlängenverhakungslimitiert. Die Struktur der Ketten in einem biszum osmotischenGleichgewicht gequollenen Netzwerk unterhalb derMaschengröße ist die überstreckter,selbstvermeidender Ketten mit einer Fraktaldimension von D =1.4,jenseits der Maschengröße (bzw. Verhakungslänge) nehmen siedie Struktur einesIrrfluges an (D = 2). Gequollene Netzwerke zeigten, wie auch in Experimenten, einestarkeZunahme von Dichtefluktuationen, welche unter Verstreckunganisotropwurde und in der Streufunktion zu sogenanntenButterfly-Mustern führt.Diese Fluktuationen sind statischer Natur als Folge desEinfrierenseines ungeordneten Zustandes während der Vernetzung.

Computational modeling of surface interactions

Die Wechselwirkung zwischen Proteinen und anorganischen Oberflächen fasziniert sowohl aus angewandter als auch theoretischer Sicht. Sie ist ein wichtiger Aspekt in vielen Anwendungen, unter anderem in chirugischen Implantaten oder Biosensoren. Sie ist außerdem ein Beispiel für theoretische Fragestellungen betreffend die Grenzfläche zwischen harter und weicher Materie. Fest steht, dass Kenntnis der beteiligten Mechanismen erforderlich ist um die Wechselwirkung zwischen Proteinen und Oberflächen zu verstehen, vorherzusagen und zu optimieren. Aktuelle Fortschritte im experimentellen Forschungsbereich ermöglichen die Untersuchung der direkten Peptid-Metall-Bindung. Dadurch ist die Erforschung der theoretischen Grundlagen weiter ins Blickfeld aktueller Forschung gerückt. Eine Möglichkeit die Wechselwirkung zwischen Proteinen und anorganischen Oberflächen zu erforschen ist durch Computersimulationen. Obwohl Simulationen von Metalloberflächen oder Proteinen als Einzelsysteme schon länger verbreitet sind, bringt die Simulation einer Kombination beider Systeme neue Schwierigkeiten mit sich. Diese zu überwinden erfordert ein Mehrskalen-Verfahren: Während Proteine als biologische Systeme ausreichend mit klassischer Molekulardynamik beschrieben werden können, bedarf die Beschreibung delokalisierter Elektronen metallischer Systeme eine quantenmechanische Formulierung. Die wichtigste Voraussetzung eines Mehrskalen-Verfahrens ist eine Übereinstimmung der Simulationen auf den verschiedenen Skalen. In dieser Arbeit wird dies durch die Verknüpfung von Simulationen alternierender Skalen erreicht. Diese Arbeit beginnt mit der Untersuchung der Thermodynamik der Benzol-Hydratation mittels klassischer Molekulardynamik. Dann wird die Wechselwirkung zwischen Wasser und den [111]-Metalloberflächen von Gold und Nickel mittels eines Multiskalen-Verfahrens modelliert. In einem weiteren Schritt wird die Adsorbtion des Benzols an Metalloberflächen in wässriger Umgebung studiert. Abschließend wird die Modellierung erweitert und auch die Aminosäuren Alanin und Phenylalanin einbezogen. Dies eröffnet die Möglichkeit realistische Protein- Metall-Systeme in Computersimulationen zu betrachten und auf theoretischer Basis die Wechselwirkung zwischen Peptiden und Oberflächen für jede Art Peptide und Oberfläche vorauszusagen.

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.

Computer simulations of microstructures related to the olivine spinel transition

Tiefherd-Beben, die im oberen Erdmantel in einer Tiefe von ca. 400 km auftreten, werden gewöhnlich mit dem in gleicher Tiefe auftretenden druckabhängigen, polymorphen Phasenübergang von Olivine (α-Phase) zu Spinel (β-Phase) in Verbindung gebracht. Es ist jedoch nach wie vor unklar, wie der Phasenübergang mit dem mechanischen Versagen des Mantelmaterials zusammenhängt. Zur Zeit werden im Wesentlichen zwei Modelle diskutiert, die entweder Mikrostrukturen, die durch den Phasenübergang entstehen, oder aber die rheologischen Veränderungen des Mantelgesteins durch den Phasenübergang dafür verantwortlich machen. Dabei sind Untersuchungen der Olivin→Spinel Umwandlung durch die Unzugänglichkeit des natürlichen Materials vollständig auf theoretische Überlegungen sowie Hochdruck-Experimente und Numerische Simulationen beschränkt. Das zentrale Thema dieser Dissertation war es, ein funktionierendes Computermodell zur Simulation der Mikrostrukturen zu entwickeln, die durch den Phasenübergang entstehen. Des Weiteren wurde das Computer Modell angewandt um die mikrostrukturelle Entwicklung von Spinelkörnern und die Kontrollparameter zu untersuchen. Die Arbeit ist daher in zwei Teile unterteilt: Der erste Teil (Kap. 2 und 3) behandelt die physikalischen Gesetzmäßigkeiten und die prinzipielle Funktionsweise des Computer Modells, das auf der Kombination von Gleichungen zur Errechnung der kinetischen Reaktionsgeschwindigkeit mit Gesetzen der Nichtgleichgewichtsthermodynamik unter nicht-hydostatischen Bedingungen beruht. Das Computermodell erweitert ein Federnetzwerk der Software latte aus dem Programmpaket elle. Der wichtigste Parameter ist dabei die Normalspannung auf der Kornoberfläche von Spinel. Darüber hinaus berücksichtigt das Programm die Latenzwärme der Reaktion, die Oberflächenenergie und die geringe Viskosität von Mantelmaterial als weitere wesentliche Parameter in der Berechnung der Reaktionskinetic. Das Wachstumsverhalten und die fraktale Dimension von errechneten Spinelkörnern ist dabei in guter Übereinstimmung mit Spinelstrukturen aus Hochdruckexperimenten. Im zweiten Teil der Arbeit wird das Computermodell angewandt, um die Entwicklung der Oberflächenstruktur von Spinelkörnern unter verschiedenen Bedigungen zu eruieren. Die sogenannte ’anticrack theory of faulting’, die den katastrophalen Verlauf der Olivine→Spinel Umwandlung in olivinhaltigem Material unter differentieller Spannung durch Spannungskonzentrationen erklärt, wurde anhand des Computermodells untersucht. Der entsprechende Mechanismus konnte dabei nicht bestätigt werden. Stattdessen können Oberflächenstrukturen, die Ähnlichkeiten zu Anticracks aufweisen, durch Unreinheiten des Materials erklärt werden (Kap. 4). Eine Reihe von Simulationen wurde der Herleitung der wichtigsten Kontrollparameter der Reaktion in monomineralischem Olivin gewidmet (Kap. 5 and Kap. 6). Als wichtigste Einflüsse auf die Kornform von Spinel stellten sich dabei die Hauptnormalspannungen auf dem System sowie Heterogenitäten im Wirtsminerals und die Viskosität heraus. Im weiteren Verlauf wurden die Nukleierung und das Wachstum von Spinel in polymineralischen Mineralparagenesen untersucht (Kap. 7). Die Reaktionsgeschwindigkeit der Olivine→Spinel Umwandlung und die Entwicklung von Spinelnetzwerken und Clustern wird durch die Gegenwart nicht-reaktiver Minerale wie Granat oder Pyroxen erheblich beschleunigt. Die Bildung von Spinelnetzwerken hat das Potential, die mechanischen Eigenschaften von Mantelgestein erheblich zu beeinflussen, sei es durch die Bildung potentieller Scherzonen oder durch Gerüstbildung. Dieser Lokalisierungprozess des Spinelwachstums in Mantelgesteinen kann daher ein neues Erklärungsmuster für Tiefbeben darstellen.

The development of stylolites, from small-scale heterogeneities to multi-scale roughness

Stylolites are rough paired surfaces, indicative of localized stress-induced dissolution under a non-hydrostatic state of stress, separated by a clay parting which is believed to be the residuum of the dissolved rock. These structures are the most frequent deformation pattern in monomineralic rocks and thus provide important information about low temperature deformation and mass transfer. The intriguing roughness of stylolites can be used to assess amount of volume loss and paleo-stress directions, and to infer the destabilizing processes during pressure solution. But there is little agreement on how stylolites form and why these localized pressure solution patterns develop their characteristic roughness.rnNatural bedding parallel and vertical stylolites were studied in this work to obtain a quantitative description of the stylolite roughness and understand the governing processes during their formation. Adapting scaling approaches based on fractal principles it is demonstrated that stylolites show two self affine scaling regimes with roughness exponents of 1.1 and 0.5 for small and large length scales separated by a crossover length at the millimeter scale. Analysis of stylolites from various depths proved that this crossover length is a function of the stress field during formation, as analytically predicted. For bedding parallel stylolites the crossover length is a function of the normal stress on the interface, but vertical stylolites show a clear in-plane anisotropy of the crossover length owing to the fact that the in-plane stresses (σ2 and σ3) are dissimilar. Therefore stylolite roughness contains a signature of the stress field during formation.rnTo address the origin of stylolite roughness a combined microstructural (SEM/EBSD) and numerical approach is employed. Microstructural investigations of natural stylolites in limestones reveal that heterogeneities initially present in the host rock (clay particles, quartz grains) are responsible for the formation of the distinctive stylolite roughness. A two-dimensional numerical model, i.e. a discrete linear elastic lattice spring model, is used to investigate the roughness evolving from an initially flat fluid filled interface induced by heterogeneities in the matrix. This model generates rough interfaces with the same scaling properties as natural stylolites. Furthermore two coinciding crossover phenomena in space and in time exist that separate length and timescales for which the roughening is either balanced by surface or elastic energies. The roughness and growth exponents are independent of the size, amount and the dissolution rate of the heterogeneities. This allows to conclude that the location of asperities is determined by a polimict multi-scale quenched noise, while the roughening process is governed by inherent processes i.e. the transition from a surface to an elastic energy dominated regime.rn

Dynamics of epithelial monolayers assessed by acoustic and impedimetric whole cell biosensors

Mit der Zielsetzung der vorliegenden Arbeit wurde die detailierten Analyse von Migrationsdynamiken epithelilaler Monolayer anhand zweier neuartiger in vitro Biosensoren verfolgt, der elektrischen Zell-Substrat Impedanz Spektroskopie (electrical cell-substrate impedance sensing, ECIS) sowie der Quarz Kristall Mikrowaage (quartz crystal microbalance, QCM). Beide Methoden erwiesen sich als sensitiv gegenüber der Zellmotilität und der Nanozytotoxizität.rnInnerhalb des ersten Projektes wurde ein Fingerprinting von Krebszellen anhand ihrer Motilitätsdynamiken und der daraus generierten elektrischen oder akkustischen Fluktuationen auf ECIS oder QCM Basis vorgenommen; diese Echtzeitsensoren wurdene mit Hilfe klassicher in vitro Boyden-Kammer Migrations- und Invasions-assays validiert. Fluktuationssignaturen, also Langzeitkorrelationen oder fraktale Selbstähnlichkeit aufgrund der kollektiven Zellbewegung, wurden über Varianz-, Fourier- sowie trendbereinigende Fluktuationsanalyse quantifiziert. Stochastische Langzeitgedächtnisphänomene erwiesen sich als maßgebliche Beiträge zur Antwort adhärenter Zellen auf den QCM und ECIS-Sensoren. Des weiteren wurde der Einfluss niedermolekularer Toxine auf die Zytoslelettdynamiken verfolgt: die Auswirkungen von Cytochalasin D, Phalloidin und Blebbistatin sowie Taxol, Nocodazol und Colchicin wurden dabei über die QCM und ECIS Fluktuationsanalyse erfasst.rnIn einem zweiten Projektschwerpunkt wurden Adhäsionsprozesse sowie Zell-Zell und Zell-Substrat Degradationsprozesse bei Nanopartikelgabe charackterisiert, um ein Maß für Nanozytotoxizität in Abhangigkeit der Form, Funktionalisierung Stabilität oder Ladung der Partikel zu erhalten.rnAls Schlussfolgerung ist zu nennen, dass die neuartigen Echtzeit-Biosensoren QCM und ECIS eine hohe Zellspezifität besitzen, auf Zytoskelettdynamiken reagieren sowie als sensitive Detektoren für die Zellvitalität fungieren können.

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.

Between the scales: Water from different perspectives

In der vorliegenden Arbeit werden verschiedene Wassermodelle in sogenannten Multiskalen-Computersimulationen mit zwei Auflösungen untersucht, in atomistischer Auflösung und in einer vergröberten Auflösung, die als "coarse-grained" bezeichnet wird. In der atomistischen Auflösung wird ein Wassermolekül, entsprechend seiner chemischen Struktur, durch drei Atome beschrieben, im Gegensatz dazu wird ein Molekül in der coarse-grained Auflösung durch eine Kugel dargestellt.rnrnDie coarse-grained Modelle, die in dieser Arbeit vorgestellt werden, werden mit verschiedenen coarse-graining Methoden entwickelt. Hierbei kommen hauptsächlich die "iterative Boltzmann Inversion" und die "iterative Monte Carlo Inversion" zum Einsatz. Beides sind struktur-basierte Ansätze, die darauf abzielen bestimmte strukturelle Eigenschaften, wie etwa die Paarverteilungsfunktionen, des zugrundeliegenden atomistischen Systems zu reproduzieren. Zur automatisierten Anwendung dieser Methoden wurde das Softwarepaket "Versatile Object-oriented Toolkit for Coarse-Graining Applications" (VOTCA) entwickelt.rnrnEs wird untersucht, in welchem Maße coarse-grained Modelle mehrere Eigenschaftenrndes zugrundeliegenden atomistischen Modells gleichzeitig reproduzieren können, z.B. thermodynamische Eigenschaften wie Druck und Kompressibilität oder strukturelle Eigenschaften, die nicht zur Modellbildung verwendet wurden, z.B. das tetraedrische Packungsverhalten, welches für viele spezielle Eigenschaft von Wasser verantwortlich ist.rnrnMit Hilfe des "Adaptive Resolution Schemes" werden beide Auflösungen in einer Simulation kombiniert. Dabei profitiert man von den Vorteilen beider Modelle:rnVon der detaillierten Darstellung eines räumlich kleinen Bereichs in atomistischer Auflösung und von der rechnerischen Effizienz des coarse-grained Modells, die den Bereich simulierbarer Zeit- und Längenskalen vergrössert.rnrnIn diesen Simulationen kann der Einfluss des Wasserstoffbrückenbindungsnetzwerks auf die Hydration von Fullerenen untersucht werden. Es zeigt sich, dass die Struktur der Wassermoleküle an der Oberfläche hauptsächlich von der Art der Wechselwirkung zwischen dem Fulleren und Wasser und weniger von dem Wasserstoffbrückenbindungsnetzwerk dominiert wird.rn

Thermodynamic concepts in adaptive resolution simulations

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

Coarse-graining and quantum-classical adaptive coupling in soft matter

Die vorliegende Arbeit untersucht den Zusammenhang zwischen Skalen in Systemen weicher Materie, der für Multiskalen-Simulationen eine wichtige Rolle spielt. Zu diesem Zweck wurde eine Methode entwickelt, die die Approximation der Separierbarkeit von Variablen für die Molekulardynamik und ähnliche Anwendungen bewertet. Der zweite und größere Teil dieser Arbeit beschäftigt sich mit der konzeptionellen und technischen Erweiterung des Adaptive Resolution Scheme'' (AdResS), einer Methode zur gleichzeitigen Simulation von Systemen mit mehreren Auflösungsebenen. Diese Methode wurde auf Systeme erweitert, in denen klassische und quantenmechanische Effekte eine Rolle spielen.rnrnDie oben genannte erste Methode benötigt nur die analytische Form der Potentiale, wie sie die meisten Molekulardynamik-Programme zur Verfügung stellen. Die Anwendung der Methode auf ein spezielles Problem gibt bei erfolgreichem Ausgang einen numerischen Hinweis auf die Gültigkeit der Variablenseparation. Bei nicht erfolgreichem Ausgang garantiert sie, dass keine Separation der Variablen möglich ist. Die Methode wird exemplarisch auf ein zweiatomiges Molekül auf einer Oberfläche und für die zweidimensionale Version des Rotational Isomer State (RIS) Modells einer Polymerkette angewandt.rnrnDer zweite Teil der Arbeit behandelt die Entwicklung eines Algorithmus zur adaptiven Simulation von Systemen, in denen Quanteneffekte berücksichtigt werden. Die Quantennatur von Atomen wird dabei in der Pfadintegral-Methode durch einen klassischen Polymerring repräsentiert. Die adaptive Pfadintegral-Methode wird zunächst für einatomige Flüssigkeiten und tetraedrische Moleküle unter normalen thermodynamischen Bedingungen getestet. Schließlich wird die Stabilität der Methode durch ihre Anwendung auf flüssigen para-Wasserstoff bei niedrigen Temperaturen geprüft.

Computersimulationen zum Einfluss topologischer Beschränkungen auf Polymere

Topologische Beschränkungen beeinflussen die Eigenschaften von Polymeren. Im Rahmen dieser Arbeit wird mit Hilfe von Computersimulationen im Detail untersucht, inwieweit sich die statischen Eigenschaften von kollabierten Polymerringen, Polymerringen in konzentrierten Lösungen und aus Polymerringen aufgebauten Bürsten mit topologischen Beschränkungen von solchen ohne topologische Beschränkungen unterscheiden. Des Weiteren wird analysiert, welchen Einfluss geometrische Beschränkungen auf die topologischen Eigenschaften von einzelnen Polymerketten besitzen. Im ersten Teil der Arbeit geht es um den Einfluss der Topologie auf die Eigenschaften einzelner Polymerketten in verschiedenen Situationen. Da allerdings gerade die effiziente Durchführung von Monte-Carlo-Simulationen von kollabierten Polymerketten eine große Herausforderung darstellt, werden zunächst drei Bridging-Monte-Carlo-Schritte für Gitter- auf Kontinuumsmodelle übertragen. Eine Messung der Effizienz dieser Schritte ergibt einen Beschleunigungsfaktor von bis zu 100 im Vergleich zum herkömmlichen Slithering-Snake-Algorithmus. Darauf folgt die Analyse einer einzelnen, vergröberten Polystyrolkette in sphärischer Geometrie hinsichtlich Verschlaufungen und Knoten. Es wird gezeigt, dass eine signifikante Verknotung der Polystrolkette erst eintritt, wenn der Radius des umgebenden Kapsids kleiner als der Gyrationsradius der Kette ist. Des Weiteren werden sowohl Monte-Carlo- als auch Molekulardynamiksimulationen sehr großer Ringe mit bis zu einer Million Monomeren im kollabierten Zustand durchgeführt. Während die Konfigurationen aus den Monte-Carlo-Simulationen aufgrund der Verwendung der Bridging-Schritte sehr stark verknotet sind, bleiben die Konfigurationen aus den Molekulardynamiksimulationen unverknotet. Hierbei zeigen sich signifikante Unterschiede sowohl in der lokalen als auch in der globalen Struktur der Ringpolymere. Im zweiten Teil der Arbeit wird das Skalierungsverhalten des Gyrationsradius der einzelnen Polymerringe in einer konzentrierten Lösung aus völlig flexiblen Polymerringen im Kontinuum untersucht. Dabei wird der Anfang des asymptotischen Skalierungsverhaltens, welches mit dem Modell des “fractal globules“ konsistent ist, erreicht. Im abschließenden, dritten Teil dieser Arbeit wird das Verhalten von Bürsten aus linearen Polymeren mit dem von Ringpolymerbürsten verglichen. Dabei zeigt sich, dass die Struktur und das Skalierungsverhalten beider Systeme mit identischem Dichteprofil parallel zum Substrat deutlich voneinander abweichen, obwohl die Eigenschaften beider Systeme in Richtung senkrecht zum Substrat übereinstimmen. Der Vergleich des Relaxationsverhaltens einzelner Ketten in herkömmlichen Polymerbürsten und Ringbürsten liefert keine gravierenden Unterschiede. Es stellt sich aber auch heraus, dass die bisher verwendeten Erklärungen zur Relaxationsverhalten von herkömmlichen Bürsten nicht ausreichen, da diese lediglich den anfänglichen Zerfall der Korrelationsfunktion berücksichtigen. Bei der Untersuchung der Dynamik einzelner Monomere in einer herkömmlichen Bürste aus offenen Ketten vom Substrat hin zum offenen Ende zeigt sich, dass die Monomere in der Mitte der Kette die langsamste Relaxation besitzen, obwohl ihre mittlere Verrückung deutlich kleiner als die der freien Endmonomere ist.

Laboratory study on the optical properties of absorbing atmospheric aerosols and field applications

The interaction between aerosols and sun light plays an important role in the radiative balance of Earth’s atmosphere. This interaction is obtained by measuring the removal (extinction), redistribution (scattering), and transformation into heat (absorption) of light by the aerosols; i.e. their optical properties. Knowledge of these properties is crucial for our understanding of the atmospheric system. rn Light absorption by aerosols is a major contributor to the direct and indirect effects on our climate system, and an accurate and sensitive measurement method is crucial to further our understanding. A homebuilt photoacoustic sensor (PAS), measuring at a 532nm wavelength, was fully characterized and its functionality validated for measurements of absorbing aerosols. The optical absorption cross-sections of absorbing polystyrene latex spheres, to be used as a standard for aerosol absorption measurements, were measured and compared to literature values. Additionally, a calibration method using absorbing aerosol of known complex refractive index was presented.rn A new approach to retrieve the effective broadband refractive indices (mbroad,eff) of aerosol particles by a white light aerosol spectrometer (WELAS) optical particle counter (OPC) was achieved. Using a tandem differential mobility analyzer (DMA)-OPC system, the nbroad,eff are obtained for both laboratory and field applications. This method was tested in the laboratory using substances with a wide range of optical properties and it was used in ambient measurements to retrieve the nbroad,eff of biomass burning aerosols in a nationwide burning event in Israel. The retrieved effective broadband refractive indices for laboratory generated scattering aerosols were: ammonium sulfate (AS), glutaric acid (GA), and sodium chloride, all within 4% of literature values. For absorbing substances, nigrosine and various mixtures of nigrosine with AS and GA were measured, as well as a lightly absorbing substance, Suwannee river fulvic acid (SRFA). For the ambient measurements, the calibration curves generated from this method were to follow the optical evolution of biomass burning (BB) aerosols. A decrease in the overall aerosol absorption and scattering for aged aerosols during the day after the fires compared to the smoldering phase of the fires was found. rn The connection between light extinction of aerosols, their chemical composition and hygroscopicity for particles with different degrees of absorption was studied. The extinction cross-section (σext) at 532nm for different mobility diameters was measured at 80% and 90% relative humidity (RH), and at an RH<10%. The ratio of the humidified aerosols to the dry ones, fRHext(%RH,Dry), is presented. For purely scattering aerosols, fRHext(%RH,Dry) is inversely proportional with size; this dependence was suppressed for lightly absorbing ones. In addition, the validity of the mixing rules for water soluble absorbing aerosols is explored. The difference between the derived and calculated real parts of the complex RIs were less than 5.3% for all substances, wavelengths, and RHs. The obtained imaginary parts for the retrieved and calculated RIs were in good agreement with each other, and well within the measurement errors of retrieval from pulsed CRD spectroscopy measurements. Finally, a core-shell structure model is also used to explore the differences between the models, for substances with low growth factors, under these hydration conditions. It was found that at 80% RH and for size parameters less than 2.5, there is less than a 5 % difference between the extinction efficiencies calculated with both models. This difference is within measurement errors; hence, there is no significant difference between the models in this case. However, for greater size parameters the difference can be up to 10%. For 90% RH the differences below a size parameter of 2.5 were up to 7%.rn Finally, the fully characterized PAS together with a cavity ring down spectrometer (CRD), were used to study the optical properties of soot and secondary organic aerosol (SOA) during the SOOT-11 project in the AIDA chamber in Karlsruhe, Germany. The fresh fractal-like soot particles were allowed to coagulate for 28 hours before stepwise coating them with SOA. The single scattering albedo for fresh fractal-like soot was measured to be 0.2 (±0.03), and after allowing the soot to coagulate for 28 hours and coating it with SOA, it increased to 0.71(±0.01). An absorption enhancement of the coated soot of up to 1.71 (±0.03) times from the non-coated coagulated soot was directly measured with the PAS. Monodisperse measurements of SOA and soot coated with SOA were performed to derive the complex refractive index (m) of both aerosols. A complex refractive index of m = 1.471(±0.008) + i0.0(±0.002) for the SOA-αO3 was retrieved. For the compact coagulated soot a preliminary complex refractive index of m = 2.04(+0.21/-0.14) + i0.34(+0.18/-0.06) with 10nm(+4/-6) coating thickness was retrieved.rn These detail properties can be use by modelers to decrease uncertainties in assessing climatic impacts of the different species and to improve weather forecasting.rn

The relevance of the silica metabolizing enzyme silicatein-alpha to biomineralization and the formation of biogenic silica in siliceous sponges

Die technische Silikatproduktion erfordert in der Regel hohe Temperaturen und extreme pH-Werte. In der Natur hingegen haben insbesondere Kieselschwämme die außergewöhnliche Fähigkeit, ihr Silikatskelett, das aus einzelnen sogenannten Spiculae besteht, enzymatisch mittels des Proteins Silicatein zu synthetisieren. rnIm Inneren der Spiculae, im zentralen Kanal, befindet sich das Axialfilament, welches hauptsächlich aus Silicatein-α aufgebaut ist. Mittels Antikörperfärbungen und Elektronenmikroskopischen Analysen konnte festgestellt werden, dass Silicatein in mit Kieselsäure-gefüllten Zellorganellen (silicasomes) nachzuweisen ist. Mittels dieser Vakuolen kann das Enzym und die Kieselsäure aus der Zelle zu den Spiculae im extrazellulären Raum befördert werden, wo diese ihre endgültige Länge und Dicke erreichen. Zum ersten Mal konnte nachgewiesen werden, dass rekombinant hergestelltes Silicatein-α sowohl als Siliciumdioxid-Polymerase als auch Siliciumdioxid-Esterase wirkt. Mittels Massenspektroskopie konnte die enzymatische Polymerisation von Kieselsäure nachverfolgt werden. Durch Spaltung der Esterbindung des künstlichen Substrates Bis(p-aminophenoxy)-dimethylsilan war es möglich kinetische Parameter der Siliciumdioxid-Esterase-Aktivität des rekombinanten Silicateins zu ermitteln.rnZu den größten biogenen Silikatstukuren auf der Erde gehören die Kieselnadeln der Schwammklasse Hexactinellida. Nadelextrakte aus den Schwammklassen Demospongien (S. domuncula) und Hexactinellida (M. chuni) wurden miteinander verglichen um die potentielle Existenz von Silicatein oder Silicatein-ähnliche Molekülen und die dazu gehörige proteolytischen Aktivität nachzuweisen. Biochemische Analysen zeigten, dass das 27 kDA große isolierte Polypeptid in Monoraphis mehrere gemeinsame Merkmale mit den Silicateinen der Demospongien teilt. Dazu gehören die Größe und die Proteinase-Aktivität. rnUm die Frage zu klären, ob das axiale Filament selbst zur Formbildung der Skelettelemente beiträgt, wurde ein neues mildes Extraktionsverfahren eingeführt. Dieses Verfahren ermöglichte die Solubilisierung des nativen Silicateins aus den Spiculae. Die isolierten Silicateine lagen als Monomere (24 kDa) vor, die Dimere durch nicht-kovalente Bindungen ausbildeten. Darüber hinaus konnten durch PAGE-Gelelektrophorese Tetramere (95 kDa) und Hexamere (135 kDa) nachgewiesen werden. Die Monomere zeigten eine beträchtliche proteolytische Aktivität, die sich während der Polymerisationsphase des Proteins weiter erhöhte. Mit Hilfe der Lichtmikroskopie und Elektronenmikroskopie (TEM) konnte die Assemblierung der Proteine zu filamentartigen Strukturen gezeigt werden. Die Selbstorganisation der Silicatein-α-Monomeren scheint eine Basis für Form- und Musterbildung der wachsenden Nadeln zu bilden.rn Um die Rolle des kürzlich entdeckten Proteins Silintaphin-1, ein starker Interaktionspartner des Silicatein-α, während der Biosilifizierung zu klären, wurden Assemblierungs-Experimente mit den rekombinanten Proteinen in vitro durchgeführt. Zusätzlich wurde deren Effekt auf die Biosilikatsynthese untersucht. Elektronenmikroskopische Analysen ergaben, dass rekombinantes Silicatein-α zufällig verteilte Aggregate bildet, während die Koinkubation beider Proteine (molekulares Verhältnis 4:1) über fraktal artige Strukturen zu Filamenten führt. Auch die enzymatische Aktivität der Silicatein-α-vermittelte Biosilikatsynthese erhöhte sich in Gegenwart von Silintaphin-1 um das 5,3-fache. rn

Polykation-DNA-Komplexe : Eigenschaften und Anwendung in der Gentransfektion

Polykationen bilden mit DNA spontan Komplexe. Triebkraft ist der Entropiegewinn durch Freisetzung der Gegenionen auf den Polyelektrolyten. Solche Komplexe können in der Gentechnik verwendet werden, um fremde DNA in eine Zelle einzuschleusen. Dies bezeichnet man als Gentransfektion. In dieser Arbeit werden erstmals bürstenförmige Polykationen mit wurmförmiger Topologie zur Gentransfektion verwendet. Dazu wurde die Komplexierung von DNA mit Bürstenpolymeren mit Poly-L-Lysin- und Polyvinylpyridinium-Seitenketten und linearen Polykationen untersucht. Die Komplexbildung verläuft in allen Fällen kinetisch kontrolliert, alle Polykationen bilden sphärische Komplexe, die Topologie hat keinen Einfluss auf die Komplexgröße. Komplexe aus Bürstenpolymeren transfizieren mehr als 25% der gesamten Zellpopulation bei Schweinehirnendothelzellen. Gegenüber dem kommerziellen Transfektionsmittel Lipofektamin konnte eine deutliche Steigerung um bis zu 400% erreicht werden. Komplexe, die mit linearen Analoga gebildet wurden, zeigten bei gleicher Komplexgröße Transfektionsraten unter 5%. Freisetzungsversuche zeigen, dass die Komplexe, die gut transfizieren, recht labil sind, also die DNA unter Kompetitoreinfluss freisetzen können. Stabile Komplexe haben geringe Transfektionseffizienzen. Ebenso wichtig ist der Schutz der DNA vor Abbau durch DNase. Die PVP-Bürste bietet als einziges der untersuchten Polykationen diesen Schutz und zeigt auch die besten Transfektionsraten. Zusätzlich zu der medizinischen Anwendung wurde die Kinetik der Komplexbildung untersucht. Dazu wurde ein spezieller Aufbau entwickelt, der es ermöglicht die Streuintensität der Komplexlösung bei kleinen Streuwinkeln zeitaufgelöst im Millisekundenbereich zu detektieren. Die Komplexbildung verläuft diffusionskontrolliert, im Bereich von Ladungsverhältnissen (positive zu negativen Ladungen) von 1.8 bis 4.0 schließt sich ein fraktales Wachstum an.