Application of reactive fluid transport modeling to hydrothermal systems

A one-dimensional multi-component reactive fluid transport algorithm, 1DREACT (Steefel, 1993) was used to investigate different fluid-rock interaction systems. A major short coming of mass transport calculations which include mineral reactions is that solid solutions occurring in many minerals are not treated adequately. Since many thermodynamic models of solid solutions are highly non-linear, this can seriously impact on the stability and efficiency of the solution algorithms used. Phase petrology community saw itself faced with a similar predicament 10 years ago. To improve performance and reliability, phase equilibrium calculations have been using pseudo compounds. The same approach is used here in the first, using the complex plagioclase solid solution as an example. Thermodynamic properties of a varying number of intermediate plagioclase phases were calculated using ideal molecular, Al-avoidance, and non-ideal mixing models. These different mixing models can easily be incorporated into the simulations without modification of the transport code. Simulation results show that as few as nine intermediate compositions are sufficient to characterize the diffusional profile between albite and anorthite. Hence this approach is very efficient, and can be used with little effort. A subsequent chapter reports the results of reactive fluid transport modeling designed to constrain the hydrothermal alteration of Paleoproterozoic sediments of the Southern Lake Superior region. Field observations reveal that quartz-pyrophyllite (or kaolinite) bearing assemblages have been transformed into muscovite-pyrophyllite-diaspore bearing assemblages due to action of fluids migrating along permeable flow channels. Fluid-rock interaction modeling with an initial qtz-prl assemblage and a K-rich fluid simulates the formation of observed mineralogical transformation. The bulk composition of the system evolves from an SiO2-rich one to an Al2O3+K2O-rich one. Simulations show that the fluid flow was up-temperature (e.g. recharge) and that fluid was K-rich. Pseudo compound approach to include solid solutions in reactive transport models was tested in modeling hydrothermal alteration of Icelandic basalts. Solid solutions of chlorites, amphiboles and plagioclase were included as the secondary mineral phases. Saline and fresh water compositions of geothermal fluids were used to investigate the effect of salinity on alteration. Fluid-rock interaction simulations produce the observed mineral transformations. They show that roughly the same alteration minerals are formed due to reactions with both types of fluid which is in agreement with the field observations. A final application is directed towards the remediation of nitrate rich groundwaters. Removal of excess nitrate from groundwater by pyrite oxidation was modeled using the reactive fluid transport algorithm. Model results show that, when a pyrite-bearing, permeable zone is placed in the flow path, nitrate concentration in infiltrating water can be significantly lowered, in agreement with proposals from the literature. This is due to nitrogen reduction. Several simulations investigate the efficiency of systems with different mineral reactive surface areas, reactive barrier zone widths, and flow rates to identify the optimum setup.

Computer simulation methods to study interfacial tensions : from the Ising model to colloidal crystals

In condensed matter systems, the interfacial tension plays a central role for a multitude of phenomena. It is the driving force for nucleation processes, determines the shape and structure of crystalline structures and is important for industrial applications. Despite its importance, the interfacial tension is hard to determine in experiments and also in computer simulations. While for liquid-vapor interfacial tensions there exist sophisticated simulation methods to compute the interfacial tension, current methods for solid-liquid interfaces produce unsatisfactory results.rnrnAs a first approach to this topic, the influence of the interfacial tension on nuclei is studied within the three-dimensional Ising model. This model is well suited because despite its simplicity, one can learn much about nucleation of crystalline nuclei. Below the so-called roughening temperature, nuclei in the Ising model are not spherical anymore but become cubic because of the anisotropy of the interfacial tension. This is similar to crystalline nuclei, which are in general not spherical but more like a convex polyhedron with flat facets on the surface. In this context, the problem of distinguishing between the two bulk phases in the vicinity of the diffuse droplet surface is addressed. A new definition is found which correctly determines the volume of a droplet in a given configuration if compared to the volume predicted by simple macroscopic assumptions.rnrnTo compute the interfacial tension of solid-liquid interfaces, a new Monte Carlo method called ensemble switch method'' is presented which allows to compute the interfacial tension of liquid-vapor interfaces as well as solid-liquid interfaces with great accuracy. In the past, the dependence of the interfacial tension on the finite size and shape of the simulation box has often been neglected although there is a nontrivial dependence on the box dimensions. As a consequence, one needs to systematically increase the box size and extrapolate to infinite volume in order to accurately predict the interfacial tension. Therefore, a thorough finite-size scaling analysis is established in this thesis. Logarithmic corrections to the finite-size scaling are motivated and identified, which are of leading order and therefore must not be neglected. The astounding feature of these logarithmic corrections is that they do not depend at all on the model under consideration. Using the ensemble switch method, the validity of a finite-size scaling ansatz containing the aforementioned logarithmic corrections is carefully tested and confirmed. Combining the finite-size scaling theory with the ensemble switch method, the interfacial tension of several model systems, ranging from the Ising model to colloidal systems, is computed with great accuracy.

Neue Methoden zur Computersimulation von Polymersystemen auf verschiedenen Längenskalen und ihre Anwendungen

In dieser Dissertation stellen wir einen neuen Ansatz zurModellierungvon Polymersystemen vor. Es werden (von methodischer Seiteher) zweiautomatisierte Iterationschemata dazu eingeführt,Kraftfeldparametermesoskopischer Polymersysteme systematisch zu optimieren:DasSimplex-Verfahren und das Struktur-Differenzen-Verfahren. Sowerdendiejenigen Freiheitsgrade aus Polymersystemen eliminiert,die einehohe Auflösung erfordern, was die Modellierung größerersystemeermöglicht. Nach Tests an einfachen Flüssigkeiten werdenvergröberteModelle von drei prototypischen Polymeren (Polyacrylsäure,Polyvinylalkohol und Polyisopren) in unterschiedlichenUmgebungen(gutes Lösungsmittel und Schmelze) entwickelt und ihrVerhalten aufder Mesoskala ausgiebig geprüft. Die zugehörige Abbildung(vonphysikalischer Seite her) so zu gestalten, daß sie dieunverwechselbaren Charakteristiken jedes systems auf diemesoskopischeLängenskala überträgt, stellt eine entscheidende Anforderungan dieautomatisierten Verfahren dar. Unsere Studien belegen, daß mesoskopische Kraftfeldertemperatur- unddichtespezifisch sind und daher bei geändernden Bedingungennachoptimiert werden müssen. Gleichzeitig läßt sichabschätzen, beiwelchen Umgebungsbedingungen dies noch nicht notwendig wird.In allenFällen reichen effektive Paarpotentiale aus, einrealistischesmesoskopisches Modell zu konstruieren. VergröberteSimulationenwerden im Falle der Polyacrylsäure erfolgreich gegenexperimentelleLichtstreudaten getestet. Wir erzielen für Molmassen bis zu300000g/mol eine hervorragende Übereinstimmung für denhydrodynamischenRadius. Unsere Ergebnisse erklären auch Korrekturen zudessenVerhalten als Funktion der Kettenlänge ('Skalenverhalten'). Im Fallevon Polyisopren untersuchen wir sowohl statische als auchdynamischeGrößen und stellen klare Unterschiede unserer Ergebnisse zudeneneines einfachen semi-flexiblen Mesoskalenmodells fest. InderProteinforschung werden aus Datenbanken gewonnene effektivePaarwechselwirkungen dazu verwendet, die freie Energie einesneuensystems vorherzusagen. Wir belegen in einem Exkurs mittelsGittersimulationen, daß es selbst in einfachsten Fällennicht gelingt,dies auch nur qualitativ korrekt zu bewerkstelligen.

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.

Multichromophore Systeme auf Basis von Rylencarbonsäureimiden

Die Leistung multichromophorer Systeme geht oftmals über die der einzelnen Chromophor-Einheiten hinaus. Ziel der vorliegenden Dissertation mit dem Titel „Multichromophore Systeme auf Basis von Rylencarbonsäureimiden“ war daher die Synthese und Charakterisierung multichromophorer Molekülarchitekturen. Die verwendeten Rylenfarbstoffe zeichnen sich durch hohe photochemische Stabilitäten sowie nahezu quantitative Fluoreszenzquantenausbeuten aus. Die optischen und elektronischen Eigenschaften multichromophorer Systeme hängen stark von der geometrischen Ordnung ab, in der die Farbstoffe zueinander stehen. Daher wurden für den Einbau formpersistente Gerüststrukturen gewählt. Der erste Teil der Arbeit beschäftigt sich mit dem Einbau ein und desselben Chromophortyps und hat neben dem Verständnis von Chromophor-Wechselwirkungen vor allem die Erhöhung des Absorptionsquerschnitts und der Fluoreszenzintensität zum Ziel. Als Gerüststruktur dienen dabei Polyphenylen-Dendrimere, Ethinyl-verbrückte Dendrimere sowie Übergangsmetall-vermittelte supramolekulare Strukturen. Aufgrund der hohen Farbstoffanzahl, des ortsdefinierten Einbaus und den hohen Fluoreszenzquantenausbeuten eignen sich diese multichromophoren Systeme als Fluoreszenzsonden und als Einzelphotonenemitter. Im zweiten Teil der Arbeit werden verschiedene Chromophortypen zu multichromophoren Systemen verknüpft, mit deren Hilfe ein vektorieller Energietransfer möglich ist. Mit Hinsicht auf die Verwendung in photovoltaischen Zellen wurde eine dendritische Triade dargestellt. Eine lineare Variante einer Rylen-Triade stellt einen molekularen Draht dar, deren Brückenelement durch eine geeignete Syntheseführung verlängert und der Energietransport daher abstandsabhängig untersucht werden kann.

Optical resonances of sphere-on-plane geometries

The optical resonances of metallic nanoparticles placed at nanometer distances from a metal plane were investigated. At certain wavelengths, these “sphere-on-plane” systems become resonant with the incident electromagnetic field and huge enhancements of the field are predicted localized in the small gaps created between the nanoparticle and the plane. An experimental architecture to fabricate sphere-on-plane systems was successfully achieved in which in addition to the commonly used alkanethiols, polyphenylene dendrimers were used as molecular spacers to separate the metallic nanoparticles from the metal planes. They allow for a defined nanoparticle-plane separation and some often are functionalized with a chromophore core which is therefore positioned exactly in the gap. The metal planes used in the system architecture consisted of evaporated thin films of gold or silver. Evaporated gold or silver films have a smooth interface with their substrate and a rougher top surface. To investigate the influence of surface roughness on the optical response of such a film, two gold films were prepared with a smooth and a rough side which were as similar as possible. Surface plasmons were excited in Kretschmann configuration both on the rough and on the smooth side. Their reflectivity could be well modeled by a single gold film for each individual measurement. The film has to be modeled as two layers with significantly different optical constants. The smooth side, although polycrystalline, had an optical response that was very similar to a monocrystalline surface while for the rough side the standard response of evaporated gold is retrieved. For investigations on thin non-absorbing dielectric films though, this heterogeneity introduces only a negligible error. To determine the resonant wavelength of the sphere-on-plane systems a strategy was developed which is based on multi-wavelength surface plasmon spectroscopy experiments in Kretschmann-configuration. The resonant behavior of the system lead to characteristic changes in the surface plasmon dispersion. A quantitative analysis was performed by calculating the polarisability per unit area /A treating the sphere-on-plane systems as an effective layer. This approach completely avoids the ambiguity in the determination of thickness and optical response of thin films in surface plasmon spectroscopy. Equal area densities of polarisable units yielded identical response irrespective of the thickness of the layer they are distributed in. The parameter range where the evaluation of surface plasmon data in terms of /A is applicable was determined for a typical experimental situation. It was shown that this analysis yields reasonable quantitative agreement with a simple theoretical model of the sphere-on-plane resonators and reproduces the results from standard extinction experiments having a higher information content and significantly increased signal-to-noise ratio. With the objective to acquire a better quantitative understanding of the dependence of the resonance wavelength on the geometry of the sphere-on-plane systems, different systems were fabricated in which the gold nanoparticle size, type of spacer and ambient medium were varied and the resonance wavelength of the system was determined. The gold nanoparticle radius was varied in the range from 10 nm to 80 nm. It could be shown that the polyphenylene dendrimers can be used as molecular spacers to fabricate systems which support gap resonances. The resonance wavelength of the systems could be tuned in the optical region between 550 nm and 800 nm. Based on a simple analytical model, a quantitative analysis was developed to relate the systems’ geometry with the resonant wavelength and surprisingly good agreement of this simple model with the experiment without any adjustable parameters was found. The key feature ascribed to sphere-on-plane systems is a very large electromagnetic field localized in volumes in the nanometer range. Experiments towards a quantitative understanding of the field enhancements taking place in the gap of the sphere-on-plane systems were done by monitoring the increase in fluorescence of a metal-supported monolayer of a dye-loaded dendrimer upon decoration of the surface with nanoparticles. The metal used (gold and silver), the colloid mean size and the surface roughness were varied. Large silver crystallites on evaporated silver surfaces lead to the most pronounced fluorescence enhancements in the order of 104. They constitute a very promising sample architecture for the study of field enhancements.

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.

Phasenverhalten, Kristallisation und Morphologie ladungsstabilisierter kolloidaler Mischungen

Die vorliegende Arbeit möchte die Anwendbarkeit ladungsstabilisierter kolloidaler Systeme als Modellsysteme für fundamentale Fragen der Festkörperphysik und Thermodynamik auf binäre Mischungen erweitern. In diesem Kontext untersucht sie das Phasenverhalten und mit ihm im Zusammenhang stehende Eigenschaften von binären Mischungen ladungsstabilisierter, sphärischer kolloidaler Partikel in wässriger Suspension. Da das Verhalten hartkugelähnlicher Systeme durch hohe Fremdionenkonzentrationen bereits gut bekannt ist, konzentriert sich diese Arbeit auf sehr langreichweitig repulsive Systeme unter deionisierten Bedingungen. Neben etablierten Methoden der Mikroskopie und statischer Lichtstreuung zur Phasendiagrammsbestimmung wird auch die Beobachtung der zeitabhängigen Entwicklung des Schermoduls verwendet, um eine langsame Erstarrungskinetik zu studieren. Es werden insbesondere Mischungen aus Komponenten unterschiedlicher Größe und Ladung der Größenverhältnisse 0,9, 0,82, 0,57, 0,39 und 0,37 untersucht. Diese zeigen in dieser Reihenfolge Phasendiagramme mit spindelförmigem fluid/kristallinen Koexistenzbereich wie auch azeotrope und eutektische Phasendiagramme. Die Strukturuntersuchungen aus der statischen Lichtstreuung stehen in praktisch allen Fällen im Einklang mit ungeordneten bcc- Substitutionskristallen, was über Modelle zu Schermodulmessungen bestätigt wird. Für das spindelförmige System wird ein überraschend weiter Koexistenzbereich beobachtet, wie er nicht von der Theorie erwartet wird. Die Lage, aber nicht die Form des Solidus stimmt quantitativ mit Simulationsvorhersagen zu einkomponentigen Systemen überein. Für das eutektische System bei einem Radienverhältnis von 0,57 wird der Einfluss der Schwerkraft auf das Phasenverhalten und die Erstarrungskinetik untersucht. Die der Kristallisation der kleineren Majoritätskomponente vorgelagerte gravitativ unterstützte Entmischung begünstigt hier die Verfestigung. Beobachtet werden Morphologien, die aus anderen Systemen bekannt sind (Facetten, Dendriten), wie auch erstmals eine kollumnare eutektische Morphologie. Aus den Ergebnissen wird der erste umfassende Überblick über das Phasenverhalten deionisierter Mischungen ladungsstabilisierter, sphärischer Partikel erstellt, die eine Diskussion der Daten anderer Autoren und unserer Gruppe über fluid-fluider Phasenseparation und einem System mit oberem azeotropen Punkt mit einschließt. Die meisten metallspezifischen Phasendiagrammtypen können mit ladungsstabilisierten kolloidalen Partikeln reproduziert werden. Die langreichweitig wechselwirkenden Partikel zeigt eine wesentlich verbesserte substitutionelle Mischbarkeit im Vergleich mit Hartkugel- und Metallsystemen. Das Größenverhältnis der sphärischen Partikel nimmt dabei die bestimmende Rolle für den Phasendiagrammtyp ein.

Magnetische Wechselwirkung Spin-Austauschgekoppelter Systeme mit Nitroxid und Nitronyl-Nitroxid Radikalen

Bei der Untersuchung molekularer magnetischer Materialien spielen Metall-Radikal Verbindungen eine bedeutende Rolle. Ein Forschungsschwerpunkt stützt sich auf die Familie der Nitronyl-Nitroxid (NIT) Radikale, die sich durch eine hohe chemische Stabilität auszeichnen. Im sogenannten „Metall-Radikal Ansatz“ wurden die starken Austauschwechselwirkungen zwischen stabilen Radikalen und Übergangsmetallionen in mehrdimensionalen Netzwerken ausgiebig untersucht. Um diese Netzwerke mit NIT Radikalen aufzubauen, müssen zusätzliche funktionelle Gruppen, mit einem Abstand zur spintragenden Einheit, in das Molekül eingebaut werden. Dies kann zu einer zusätzlichen schwachen Spinaustauschwechselwirkung führen. Um diese Wechselwirkung zwischen Metalldimeren mit einem einzelnen Benzoat annalogen NIT-Radikal zu untersuchen, wurden dimere Mangan(II), Kobalt(II) und Zink(II) Komplexe mit dem Chelatliganden N,N,N',N'-Tetrakis(2-benzimid-azolylalkyl)-2-hydroxy-1,3-diamino-propan synthetisiert und zusätzlich über eine periphere Carboxylat Gruppe eines NIT Radikals verbrückt.rnDie Messungen der magnetischen Suszeptibilität weisen auf eine dominante antiferromagnetische Wechselwirkung in der Metall-Radikal Verbindung hin, bei der es sich um die Spin-Austauschwechselwirkung innerhalb des Metalldimers handelt. Durch den Vergleich mit analogen Nitrobenzoat- verbrückten Mangan(II) und Kobalt(II) Verbindungen konnte gezeigt werden, dass keine Metall-Radikal Wechselwirkung beobachtet wird, obwohl eine Wechselwirkung der pi*-orbitale mit den delokalisierten pi-System des Phenylrings durch Spin-Polarisation grundsätzlich möglich ist. Auch ESR - Messungen bestätigen dies, da der Spingrundzustand das anisotrope Signal des freien NIT Radikals aufweist. Das Radikal verhält sich somit wie ein isoliertes S=1/2 Spin-Zentrum, was zusätzlich durch DFT-Rechnungen bekräftigt werden konnte. Zusammenfassend führt also die Koordination eines NIT-Benzoats an ein antiferromagnetisch gekoppeltes Metalldimer nur zur Anhebung des Spingrundzustandes und hat keinen signifikanten Effekt auf die Austauschwechselwirkung. Um trotzdem eine Metall-Radikal Wechselwirkung beobachten zu können, ist es notwendig Koordinationsverbindungen zu synthetisieren in denen hohe Spingrundzustände besetzt werden. Dies trifft auf das analoge Kupferdimer zu, wofür eine ferromagnetische Wechselwirkung zu beobachten ist.rnNach den Regeln der Spin-Polarisation müsste die Verkürzung des Austauschpfades um eine Bindung zu einer Umkehrung des Vorzeichens der magnetischen Wechselwirkung führen. Diese Verkürzung kann man durch die Verwendung des alternativen stabilen NOA-Radikals (tert-Butyl Nitroxid) erreichen. Sowohl das NIT als auch das NOA-Radikal werden an ein Kupfer(II)-dimer koordiniert, das durch die Verwendung des oben erwähnten N6O-Liganden gebildet wurde. In der Modellverbindung, ohne einen paramagnetischen Substituenten am Benzoat, zeigen die Kupferionen eine ferromagnetische Wechselwirkung mit einem Triplett Grundzustand, dessen Existenz durch die Messung der magnetischen Suszeptibilität und ESR-Spektroskopie belegt werden kann. Aufgrund der nahezu identischen Koordinationsumgebung bleibt bei allen synthetisierten Verbindungen die Kupfer-Kupfer Wechselwirkung dabei gleich. Die Daten von ESR und magnetischen Messungen zeigen weiterhin auf eine signifikante zusätzliche Metall-Radikal Wechselwirkung hin. Bei der NIT-Verbindung ist diese Austauschwechselwirkung schwach antiferromagnetisch, während die NOA-Verbindung eine schwache ferromagnetische Kopplung aufzeigt. Diese Resultate können durch DFT Rechnungen bekräftigt werden. Der Vorzeichenwechsel des Kopplungsparameters kann durch die Verkürzung des Austauschpfades vom NIT zum NOA-Benzoat um eine Bindung erklärt werden. Durch die Wahl von geeigneten Radikal- Liganden und Metallionen, zeigt sich die Möglichkeit, Systeme zu erzeugen, in denen die Radikal-Metall Wechselwirkung auch über größere Distanzen den Spin-Grundzustand des gesamten Systems signifikant beeinflussen kann. die Anwendung dieses Konzeptes auf Metall-Radikal Cluster System sollte Von großem Interesse sein.rn

Computer simulations of colloidal fluids in confinement

Computer-Simulationen von Kolloidalen Fluiden in Beschränkten Geometrien Kolloidale Suspensionen, die einen Phasenübergang aufweisen, zeigen eine Vielfalt an interessanten Effekten, sobald sie auf eine bestimmte Geometrie beschränkt werden, wie zum Beispiel auf zylindrische Poren, sphärische Hohlräume oder auf einen Spalt mit ebenen Wänden. Der Einfluss dieser verschiedenen Geometrietypen sowohl auf das Phasenverhalten als auch auf die Dynamik von Kolloid-Polymer-Mischungen wird mit Hilfe von Computer-Simulationen unter Verwendung des Asakura-Oosawa- Modells, für welches auf Grund der “Depletion”-Kräfte ein Phasenübergang existiert, untersucht. Im Fall von zylindrischen Poren sieht man ein interessantes Phasenverhalten, welches vom eindimensionalen Charakter des Systems hervorgerufen wird. In einer kurzen Pore findet man im Bereich des Phasendiagramms, in dem das System typischerweise entmischt, entweder eine polymerreiche oder eine kolloidreiche Phase vor. Sobald aber die Länge der zylindrischen Pore die typische Korrelationslänge entlang der Zylinderachse überschreitet, bilden sich mehrere quasi-eindimensionale Bereiche der polymerreichen und der kolloidreichen Phase, welche von nun an koexistieren. Diese Untersuchungen helfen das Verhalten von Adsorptionshysteresekurven in entsprechenden Experimenten zu erklären. Wenn das Kolloid-Polymer-Modellsystem auf einen sphärischen Hohlraum eingeschränkt wird, verschiebt sich der Punkt des Phasenübergangs von der polymerreichen zur kolloidreichen Phase. Es wird gezeigt, dass diese Verschiebung direkt von den Benetzungseigenschaften des Systems abhängt, was die Beobachtung von zwei verschiedenen Morphologien bei Phasenkoexistenz ermöglicht – Schalenstrukturen und Strukturen des Janustyps. Im Rahmen der Untersuchung von heterogener Keimbildung von Kristallen innerhalb einer Flüssigkeit wird eine neue Simulationsmethode zur Berechnung von Freien Energien der Grenzfläche zwischen Kristall- bzw. Flüssigkeitsphase undWand präsentiert. Die Resultate für ein System von harten Kugeln und ein System einer Kolloid- Polymer-Mischung werden anschließend zur Bestimmung von Kontaktwinkeln von Kristallkeimen an Wänden verwendet. Die Dynamik der Phasenseparation eines quasi-zweidimensionalen Systems, welche sich nach einem Quench des Systems aus dem homogenen Zustand in den entmischten Zustand ausbildet, wird mit Hilfe von einer mesoskaligen Simulationsmethode (“Multi Particle Collision Dynamics”) untersucht, die sich für eine detaillierte Untersuchung des Einflusses der hydrodynamischen Wechselwirkung eignet. Die Exponenten universeller Potenzgesetze, die das Wachstum der mittleren Domänengröße beschreiben, welche für rein zwei- bzw. dreidimensionale Systeme bekannt sind, können für bestimmte Parameterbereiche nachgewiesen werden. Die unterschiedliche Dynamik senkrecht bzw. parallel zu den Wänden sowie der Einfluss der Randbedingungen für das Lösungsmittel werden untersucht. Es wird gezeigt, dass die daraus resultierende Abschirmung der hydrodynamischen Wechselwirkungsreichweite starke Auswirkungen auf das Wachstum der mittleren Domänengröße hat.

Interaktions- und Aggregationsverhalten von Nanopartikeln in Blutserum

Das Interesse an nanopartikulären Wirkstoffsystemen steigt sowohl auf universitärer als auch auf industrieller Seite stetig an. Da diese Formulierungen meist intravenös verabreicht werden, kommt es folglich zu einem direkten Kontakt der Nanopartikel mit den Blutbestandteilen. Adsorption von Plasma Proteinen kann eine deutliche Veränderung der charakteristischen Eigenschaften des Systems induzieren, was dann Wirkungsweise sowie Toxizität stark beinflussen kann. Derzeit findet die Charakterisierung nanopartikulärer Wirkstoffsysteme vor der in vivo Applikation in Pufferlösungen mit physiologischem Salzgehalt statt, es ist jedoch kaum etwas bekannt über deren Wechselwirkungen mit komplexen Proteinmischungen wie sie im Blutserum oder –plasma vorliegen. rnMittels dynamischer Lichtstreuung (DLS) wurde eine einfache und reproduzierbare Methode entwickelt um die Aggregatbildung zwischen Nanopartikeln, Polymeren oder verschiedenen Wirkstoff-Konjugaten in humanem Blutserum zu untersuchen. Die Anwendbarkeit dieser Methode wurde durch Untersuchung verschiedener potentieller Nanotherapeutika (z.B. Polystyrol-Nanokapseln, Liposomen, amphiphile Blockcopolymere und Nanohydrogele) bezüglich ihrer Aggregation in humanem Blutserum mittels DLS gezeigt und teilweise mit aktuellen in vivo Experimenten verglichen. rnDarüber hinaus wurden größeneinheitliche Liposomen, basierend auf Disteraoylphosphatidylcholin, Cholesterol und einem Spermin-Tensid, hergestellt. Die Einkapselung von siRNA ist, je nach Präparationsmethode, mit Einkapselungseffizienzen von 40-75% möglich. Nach detaillierter Charakterisierung der Liposomen wurden diese ebenfalls bezüglich ihres Aggregationsverhaltens in humanem Blutserum untersucht. Unbeladene Liposomen aggregieren nicht mit Komponenten des Serums. Je nach Beladungsprotokoll können aggregierende sowie nicht aggregierende Liposomen-siRNA Komplexe hergestellt werden. rnWeiterhin wurden, zur Identifikation der Aggregation induzierenden Serumkomponenten, verschiedenen Serumsfraktionierungstechniken erfolgreich angewendet. Albumin, IgG, und Lipoproteine (VLDL, LDL) sowie verschiedene Proteinmischungen konnten isoliert und für weitere Aggregationsstudien mittels DLS verwendet werden. Für einige ausgesuchte Systeme konnten die Interaktionspartner identifiziert werden. rnDie Korrelation des Aggregationsverhaltens mit den strukturellen sowie funktionellen Eigenschaften der untersuchten Nanopartikel führt zu dem generellen Ergebnis, dass leicht negativ und leicht positiv bis neutrale Partikel eine geringe Tendenz zur Aggregation in Serum haben. Auch zwitterionische Substanzen zeigen eine hohe Serumstabilität. Hingegen aggregieren stark positiv und negativ geladene Partikel vermehrt.rn

Oberflächenstrukturierte amphiphile Polyphenylendendrimere zur Imitation natürlicher Transportproteine

In dieser Arbeit werden formstabile, amphiphile, oberflächenstrukturierte Polyphenylendendrimere (PPDs) mit verschiedenen Oberflächenpolaritäten beschrieben. Die physikalisch-chemischen Eigenschaften dieser Makromoleküle wurden studiert, welche ein gutes Verständnis der Nanoumgebung amphiphiler PPDs lieferten. Auch lichtinduzierte Polaritätsänderung wurde untersucht. Mit dem Konzept einer gleichmäßigen Verteilung polarer Bereiche auf der Peripherie hydrophober PPPs gelang es, Transportsysteme für Fettsäuren und Zytostatika zu erzeugen, welche charakteristische Merkmale natürlicher Transportproteine wie Albumin in sich vereinen. Hierzu zählen eine stabile dreidimensionale Form, die Ausbildung von Bindungstaschen sowie eine definierte strukturierte Oberfläche aus hydrophilen und hydrophoben Bereichen. Die Verfügbarkeit von lipophilen Bindungstaschen übertrifft sogar die des Albumins. Im Gegensatz zu Polymeren kann die Wirkstoffaufnahme bei PPDs exakt bestimmt werden. Die Anpassung der peripheren Gruppen beeinflusst den zellulären Aufnahmemechanismus. Es konnten effiziente Zellaufnahmen in A549-Zellen sowie der Transport und die intrazelluläre Freisetzung von Doxorubicin erreicht werden. Manche PPDs bieten eine Größe und Architektur, die es ermöglicht, Endothelzellen des Gehirns zu durchdringen. Es wurde auch der andere Extremfall untersucht, indem alle polaren Gruppen auf einer Hemisphäre akkumuliert wurden. Zur Darstellung solcher Janus-Dendrimere wurde ein neues Synthesekonzept herausgearbeitet und die erhaltenen Janus-Dendrimere mittels Lichtstreuung untersucht, wobei definierte perlenschnurartige Aggregate gefunden wurden. Weiterhin wurden semifluorierte Amphiphile vorgestellt, welche die Möglichkeit zur Selbstorganisation durch Nanophasenseparation bieten.

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.