Patterns in and kinetics of phase separation in binary mixtures

This work focused mainly on two aspects of kinetics of phase separation in binary mixtures. In the first part, we studied the interplay of hydrodynamics and the phase separation of binary mixtures. A considerably flat container (a laterally extended geometry), at an aspect ratio of 14:1 (diameter: height) was chosen, so that any hydrodynamic instabilities, if they arise, could be tracked. Two binary mixtures were studied. One was a mixture of methanol and hexane, doped with 5% ethanol, which phase separated under cooling. The second was a mixture of butoxyethanol and water, doped with 2% decane, which phase separated under heating. The dopants were added to bring down the phase transition temperature around room temperature.rnrnAlthough much work has been done already on classical hydrodynamic instabilities, not much has been done in the understanding of the coupling between phase separation and hydrodynamic instabilities. This work aimed at understanding the influence of phase separation in initiating any hydrodynamic instability, and also vice versa. Another aim was to understand the influence of the applied temperature protocol on the emergence of patterns characteristic to hydrodynamic instabilities. rnrnOn slowly cooling the system continuously, at specific cooling rates, patterns were observed in the first mixture, at the start of phase separation. They resembled the patterns observed in classical Rayleigh-Bénard instability, which arises when a liquid continuously is heated from below. To suppress this classical convection, the cooling setup was tuned such that the lower side of the sample always remained cooler by a few millikelvins, relative to the top. We found that the nature of patterns changed with different cooling rates, with stable patterns appearing for a specific cooling rate (1K/h). On the basis of the cooling protocol, we estimated a modified Rayleigh number for our system. We found that the estimated modified Rayleigh number is near the critical value for instability, for cooling rates between 0.5K/h and 1K/h. This is consistent with our experimental findings. rnrnThe origin of the patterns, in spite of the lower side being relatively colder with respect to the top, points to two possible reasons. 1) During phase separation droplets of either phases are formed, which releases a latent heat. Our microcalorimetry measurements show that the rise in temperature during the first phase separation is in the order of 10-20millikelvins, which in some cases is enough to reverse the applied temperature bias. Thus phase separation in itself initiates a hydrodynamic instability. 2) The second reason comes from the cooling protocol itself. The sample was cooled from above and below. At sufficiently high cooling rates, there are situations where the interior of the sample is relatively hotter than both top and bottom of the sample. This is sufficient to create an instability within the cell. Our experiments at higher cooling rates (5K/h and above) show complex patterns, which hints that there is enough convection even before phase separation occurs. Infact, theoretical work done by Dr.Hayase show that patterns could arise in a system without latent heat, with symmetrical cooling from top and bottom. The simulations also show that the patterns do not span the entire height of the sample cell. This is again consistent with the cell sizes measured in our experiment.rnrnThe second mixture also showed patterns at specific heating rates, when it was continuously heated inducing phase separation. In this case though, the sample was turbid for a long time until patterns appeared. A meniscus was most probably formed before the patterns emerged. We attribute the reason of patterns in this case to Marangoni convection, which is present in systems with an interface, where local differences in surface tension give rise to an instability. Our estimates for the Rayleigh number also show a significantly lower number than that's required for RB-type instability.rnrnIn the first part of the work, therefore, we identify two different kinds of hydrodynamic instabilities in two different mixtures. Both are observed during, or after the first phase separation. Our patterns compare with the classical convection patterns, but here the origins are from phase separation and the cooling protocol.rnrnIn the second part of the work, we focused on the kinetics of phase separation in a polymer solution (polystyrene and methylcyclohexane), which is cooled continuously far down into the two phase region. Oscillations in turbidity, denoting material exchange between the phases are seen. Three processes contribute to the phase separation: Nucleation of droplets, their growth and coalescence, and their subsequent sedimentation. Experiments in low molecular binary mixtures had led to models of oscillation [43] which considered sedimentation time scales much faster than the time scales of nucleation and growth. The size and shape of the sample therefore did not matter in such situations. The oscillations in turbidity were volume-dominated. The present work aimed at understanding the influence of sedimentation time scales for polymer mixtures. Three heights of the sample with same composition were studied side by side. We found that periods increased with the sample height, thus showing that sedimentation time determines the period of oscillations in the polymer solutions. We experimented with different cooling rates and different compositions of the mixture, and we found that periods are still determined by the sample height, and therefore by sedimentation time. rnrnWe also see that turbidity emerges in two ways; either from the interface, or throughout the sample. We suggest that oscillations starting from the interface are due to satellite droplets that are formed on droplet coalescence at the interface. These satellite droplets are then advected to the top of the sample, and they grow, coalesce and sediment. This type of an oscillation wouldn't require the system to pass the energy barrier required for homogenous nucleation throughout the sample. This mechanism would work best in sample where the droplets could be effectively advected throughout the sample. In our experiments, we see more interface dominated oscillations in the smaller cells and lower cooling rates, where droplet advection is favourable. In larger samples and higher cooling rates, we mostly see that the whole sample becomes turbid homogenously, which requires the system to pass the energy barrier for homogenous nucleation.rnrnOscillations, in principle, occur since the system needs to pass an energy barrier for nucleation. The height of the barrier decreases with increasing supersaturation, which in turn is from the temperature ramp applied. This gives rise to a period where the system is clear, in between the turbid periods. At certain specific cooling rates, the system can follow a path such that the start of a turbid period coincides with the vanishing of the last turbid period, thus eliminating the clear periods. This means suppressions of oscillations altogether. In fact we experimentally present a case where, at a certain cooling rate, oscillations indeed vanish. rnrnThus we find through this work that the kinetics of phase separation in polymer solution is different from that of a low molecular system; sedimentation time scales become relevant, and therefore so does the shape and size of the sample. The role of interface in initiating turbid periods also become much more prominent in this system compared to that in low molecular mixtures.rnrnIn summary, some fundamental properties in the kinetics of phase separation in binary mixtures were studied. While the first part of the work described the close interplay of the first phase separation with hydrodynamic instabilities, the second part investigated the nature and determining factors of oscillations, when the system was cooled deep into the two phase region. Both cases show how the geometry of the cell can affect the kinetics of phase separation. This study leads to further fundamental understandings of the factors contributing to the kinetics of phase separation, and to the understandings of what can be controlled and tuned in practical cases. rn

Atomistic simulations of cosolvent/water mixtures

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

Monte Carlo-Simulationen zum kritischen Verhalten dünner Ising-Filme

Monte Carlo-Simulationen zum kritischen Verhalten dünnerIsing-Filme Dünne Ising-Filme können als vereinfachtes Modell zurBeschreibung von binären Mischungen oder von Flüssigkeitenin schlitzartigen Kapillaren dienen. Infolge dereingeschränkten Geometrie unterscheidet sich das kritischeVerhalten dieser Systeme signifikant von dem einesBulk-Systems, es kommt zu einem Crossover von zwei- zudreidimensionalem kritischen Verhalten. Zusätzlichverschiebt sich der Phasenübergang in den ungesättigtenBereich, ein Effekt, der als 'capillary condensation'bezeichnet wird. In der vorliegenden Arbeit wurden die kritischenEigenschaften von Ising-Filmen im Rahmen einer MonteCarlo-Simulation untersucht. Zur Verbesserung der Effizienzwurde ein Cluster-Algorithmus verwendet, der um einenGhost-Spin-Term zur Behandlung der Magnetfelder erweitertwar. Bei der Datenanalyse kamen moderneMulti-Histogramm-Techniken zur Anwendung. Für alle untersuchten Schichtdicken konnten kritischeTemperatur und Magnetfeld sehr präzise bestimmt werden. DieSkalenhypothese von Fisher und Nakanishi, die dieVerschiebung des kritischen Punktes gegenüber seinesBulk-Wertes beschreibt, wurde sowohl für Systeme mit freienOberflächen als auch für Systeme mit schwachemOberflächenfeld bestätigt. Der Wert des Gap-Exponenten derOberfläche wurde mit $Delta_1$=0.459(13) in Übereinstimmungmit den Literaturwerten abgeschätzt. Die Observablen Magnetisierung und magnetischeSuszeptibilität sowie deren auf die Oberfläche bezogenenEntsprechungen zeigen kein reines zweidimensionaleskritisches Verhalten. Zu ihrer Beschreibung in der Nähe deskritischen Punktes wurden effektive Exponenten für dieeinzelnen Schichtdicken bestimmt.

Ternäre polymerhaltige Lösungen - Phasenverhalten und Grenzflächenspannung

Im Rahmen dieser Arbeit wurden experimentelle und theoretische Untersuchungen zum Phasen- und Grenzflächenverhalten von ternären Systemen des Typs Lösungsmittel/Fällungsmittel/Polymer durchgeführt. Diese Art der Mischungen ist vor allem für die Planung und Durchführung der Membranherstellung von Bedeutung, bei der die genaue Kenntnis des Phasendiagramms und der Grenzflächenspannung unabdingbar ist. Als Polymere dienten Polystyrol sowie Polydimethylsiloxan. Im Fall des Polystyrols kam Butanon-2 als Lösungsmittel zum Einsatz, wobei drei niedrigmolekulare lineare Alkohole als Fällungsmittel verwendet wurden. Für Polydimethylsiloxan eignen sich Toluol als Lösungsmittel und Ethanol als Fällungsmittel. Durch Lichtstreumessungen, Dampfdruckbestimmungen mittels Headspace-Gaschromatographie (VLE-Gleichgewichte) sowie Quellungsgleichgewichten lassen sich die thermodynamischen Eigenschaften der binären Subsysteme charakterisieren. Auf Grundlage der Flory-Huggins-Theorie kann das experimentell bestimmte Phasenverhalten (LLE-Gleichgewichte) in guter Übereinstimmung nach der Methode der Direktminimierung der Gibbs'schen Energie modelliert werden. Zieht man die Ergebnisse der Aktivitätsbestimmung von Dreikomponenten-Mischungen mit in Betracht, so ergeben sich systematische Abweichungen zwischen Experiment und Theorie. Sie können auf die Notwendigkeit ternärer Wechselwirkungsparameter zurückgeführt werden, die ebenfalls durch Modellierung zugänglich sind.Durch die aus den VLE- und LLE-Untersuchungen gewonnenen Ergebnissen kann die sog. Hump-Energie berechnet werden, die ein Maß für die Entmischungstendenz darstellt. Diese Größe eignet sich gut zur Beschreibung von Grenzflächenphänomenen mittels Skalengesetzen. Die für binäre Systeme gefundenen theoretisch fundierten Skalenparameter gelten jedoch nur teilweise. Ein neues Skalengesetz lässt erstmals eine Beschreibung über die gesamte Mischungslücke zu, wobei ein Parameter durch eine gemessene Grenzflächenspannung (zwischen Fällungsmittel/Polymer) ersetzt werden kann.

Untersuchungen zur Verfestigungskinetik in Suspensionen kolloidaler Partikel

Kolloidale Suspensionen, bei denen man die kolloidalen Teilchen als "Makroatome" in einem Kontinuum aus Lösungsmittelmolekülen auffaßt, stellen ein geeignetes Modellsystem zur Untersuchung von Verfestigungsvorgängen dar. Auf Grund der typischen beteiligten Längen- und Zeitskalen können Phasenübergänge bequem mit optischen Verfahren studiert werden. In der vorliegenden Arbeit wurde die Kinetik der Kristallisation in drei kolloidalen Systemen unterschiedlicher Teilchen-Teilchen-Wechselwirkung mit Lichtstreu- und mikroskopischen Methoden untersucht. Zur Untersuchung von Suspensionen aus sterisch stabilisierten PMMA-Teilchen, die in guter Näherung wie harte Kugeln wechselwirken, wurde ein neuartiges Laserlichtstreuexperiment aufgebaut, das die gleichzeitige Detektion von Bragg- und Kleinwinkelstreuung an einer Probe erlaubt. Damit konnte der zeitliche Verlauf der Kristallisation verfolgt sowie u.a. Nukleationsraten und erstmals auch Wachstumsgeschwindigkeiten bestimmt werden; diese wurden mit klassischer Nukleationstheorie sowie Wilson-Frenkel-Wachstum verglichen. In beiden Fällen konnte sehr gute Übereinstimmung mit der Theorie festgestellt werden. In Systemen geladener Partikel wurden mit Bragg-Mikroskopie die Wachstumsgeschwindigkeiten heterogener, an der Wand der Probenzelle aufwachsender Kristalle untersucht. Die Anpassung eines Wilson-Frenkel-Wachstumsgesetzes gelingt auch hier, wenn man die dazu eingeführte reskalierte Energiedichte auf den Schmelzpunkt bezieht. Geeignete Reskalierung der Daten erlaubt den Vergleich mit den Hartkugelsystemen. Zum ersten Mal wurde die Kristallisationskinetik in zwei verschiedenen kolloidalen binären Mischungen bestimmt und ausgewertet: In Beimischungen einer nichtkristallisierenden Teilchensorte zu einer kristallisierenden Suspension konnten die Daten mit einem modifizierten Wilson-Frenkel-Gesetz beschrieben werden, während in Mischungen aus zwei kristallisierenden Partikelsystemen eine unerwartet hohe Abnahme der Wachstumsgeschwindigkeiten beobachtet wurde. Kolloidale Suspensionen hartkugelähnlicher Mikrogel-Partikel konnten mit Hilfe des Lichtstreuaufbaues ebenfalls zum ersten Mal untersucht werden. Es wurde eine ähnliche Kristallisationskinetik wie in den PMMA-Systemen gefunden, jedoch auch einige wichtige Unterschiede, die insbesondere den Streumechanismus im Kleinwinkelbereich betrafen. Hier wurden verschiedene Interpretationsvorschläge diskutiert.

Effect of mixing, particle interaction, and spatial dimension on the glass transition

In this thesis, the influence of composition changes on the glass transition behavior of binary liquids in two and three spatial dimensions (2D/3D) is studied in the framework of mode-coupling theory (MCT).The well-established MCT equations are generalized to isotropic and homogeneous multicomponent liquids in arbitrary spatial dimensions. Furthermore, a new method is introduced which allows a fast and precise determination of special properties of glass transition lines. The new equations are then applied to the following model systems: binary mixtures of hard disks/spheres in 2D/3D, binary mixtures of dipolar point particles in 2D, and binary mixtures of dipolar hard disks in 2D. Some general features of the glass transition lines are also discussed. The direct comparison of the binary hard disk/sphere models in 2D/3D shows similar qualitative behavior. Particularly, for binary mixtures of hard disks in 2D the same four so-called mixing effects are identified as have been found before by Götze and Voigtmann for binary hard spheres in 3D [Phys. Rev. E 67, 021502 (2003)]. For instance, depending on the size disparity, adding a second component to a one-component liquid may lead to a stabilization of either the liquid or the glassy state. The MCT results for the 2D system are on a qualitative level in agreement with available computer simulation data. Furthermore, the glass transition diagram found for binary hard disks in 2D strongly resembles the corresponding random close packing diagram. Concerning dipolar systems, it is demonstrated that the experimental system of König et al. [Eur. Phys. J. E 18, 287 (2005)] is well described by binary point dipoles in 2D through a comparison between the experimental partial structure factors and those from computer simulations. For such mixtures of point particles it is demonstrated that MCT predicts always a plasticization effect, i.e. a stabilization of the liquid state due to mixing, in contrast to binary hard disks in 2D or binary hard spheres in 3D. It is demonstrated that the predicted plasticization effect is in qualitative agreement with experimental results. Finally, a glass transition diagram for binary mixtures of dipolar hard disks in 2D is calculated. These results demonstrate that at higher packing fractions there is a competition between the mixing effects occurring for binary hard disks in 2D and those for binary point dipoles in 2D.

Nano- und mikrostrukturierte organische Erkennungsstrukturen in Lipidmembranen

Im Rahmen dieser Arbeit wurden zunächst die ausgewählten Lipide anhand ihrer p-ArnIsothermen charakterisiert. Dabei zeigte sich, dass die flüssiganaloge Phase umso ausgeprägter ist, je größer die hydrophile Kopfgruppe des jeweiligen Lipids ist, so dass diese schon bei sehr großen Flächen miteinander wechselwirken. Außerdem wurde das Mischungsverhalten und die Domänenbildung binärer und ternärer Mischungen untersucht.rnDabei konnte im Fall der binären Mischungen der polymerisierbaren Lipide 1 und 2 mit dem perfluorierten Lipid 3 eine Phasenseparation beobachtet werden. Hierbei zeigte sich, dass sowohl die Form als auch die Größe der Domänen durch zwei Faktoren kontrolliert werden konnte, nämlich zum einen durch den Lateraldruck und zum anderen durch die molare Zusammensetzung der binären Mischung. Bei niedrigen Lateraldrücken (10 mN/m) warenrnfluorierte und nichtfluorierte Lipide in der homogenen flüssiganalogen Phase vollständig mischbar. Bei sehr hohen Lateraldrücken (30 mN/m) lag eine nahezu vollständige Entmischung vor, da das Lipid 1 hier schon vollständig auskristallisiert war. Im Fall der binären Mischung 2/3 konnte ein solches Verhalten nicht beobachtet werden, was an dem Phasenverhalten des Lipids 2 liegt. Aber auch hier lies sich die Domänenbildung sowohl überrnden Lateraldruck als auch über die Komposition der binären Mischung kontrollieren. Des Weiteren konnte gezeigt werden, dass die Monoschichten dieser beiden binären Mischungen durch photochemische Polymerisation stabilisiert werden können.rnDer Einfluss des Fluorierungsgrades auf das Phasenseparationsverhalten wurde anhand der binären Mischungen 1/6 und 1/8 untersucht. Es zeigte sich hier, dass sich auch hier in beiden Fällen die Domänenbildung sowohl über den Lateraldruck als auch über die Komposiion der Mischung steuern lies. Ein Vergleich der gefundenen Strukturen belegte darüber hinaus klare Unterschiede zwischen dem fluorierten und teilfluorierten Lipid. Denn in der binären Mischung 1/8 (perfluoriertes Lipid) konnte die hexagonale Grundstruktur, die zuvor schon inrnder binären Mischung 1/3 gefunden wurde, beobachtet werden. Demgegenüber hatten die kristallinen Domänen, die in der binären Mischung 1/6 (teilfluoriertes Lipid) beobachtet wurden, die hexagonale Grundform verloren. Hier wurden rundliche unförmige Domänen gefunden, die zudem deutlich kleiner waren als die bei der Mischung mit dem perfluorierten Lipid gefundenen Domänen. In den ternären Mischungen zeigte sich, dass diese ternären Mischungen eine ähnliche Phasenseparation zeigten wie die binäre Mischung 1/3. Die hier auftretenden kristallanalogen Domänen bestanden aus Lipid 1, wobei die umgebendernflüssiganaloge Matrix die beiden fluorierten Lipide enthielt. Deswegen war es von besonderer Bedeutung wie sich die beiden fluorierten Verbindungen untereinander verhalten würden.rnDazu wurden die p-A Isothermen der binären Mischung 3/5 ausgewertet und dabei zeigte sich, dass sich das Vakzinlipid, sofern der Anteil nicht höher ist als der des fluorierten Matrixlipides, in idealer Weise mit diesem mischt. Eine Untersuchung der Antigen-Antikörper Erkennung der binären Mischungen 3/5 zeigte, dass die Reaktion mit dem Antikörper umso effektiver war, je größer der Anteil des Vakzinlipids 5 war. Allerdings konnte hier kein linearer Zusammenhang gefunden werden. Stattdessen handelte es sich umrneine Sättigungskurve, da bereits ein Anteil an Lipid 5 von 42,2 mol% eine nur geringfügig schwächere Antikörperanbindung zeigte, als das reine Vakzinlipid.

Dynamics, rheology and critical properties of colloidal fluid mixtures

Liquids under the influence of external fields exhibit a wide range of intriguing phenomena that can be markedly different from the behaviour of a quiescent system. This work considers two different systems — a glassforming Yukawa system and a colloid-polymer mixture — by Molecular Dynamics (MD) computer simulations coupled to dissipative particle dynamics. The former consists of a 50-50 binary mixture of differently-sized, like-charged colloids interacting via a screened Coulomb (Yukawa) potential. Near the glass transition the influence of an external shear field is studied. In particular, the transition from elastic response to plastic flow is of interest. At first, this model is characterised in equilibrium. Upon decreasing temperature it exhibits the typical dynamics of glassforming liquids, i.e. the structural relaxation time τα grows strongly in a rather small temperature range. This is discussed with respect to the mode-coupling theory of the glass transition (MCT). For the simulation of bulk systems under shear, Lees-Edwards boundary conditions are applied. At constant shear rates γ˙ ≫ 1/τα the relevant time scale is given by 1/γ˙ and the system shows shear thinning behaviour. In order to understand the pronounced differences between a quiescent system and a system under shear, the response to a suddenly commencing or terminating shear flow is studied. After the switch-on of the shear field the shear stress shows an overshoot, marking the transition from elastic to plastic deformation, which is connected to a super-diffusive increase of the mean squared displacement. Since the average static structure only depends on the value of the shear stress, it does not discriminate between those two regimes. The distribution of local stresses, in contrast, becomes broader as soon as the system starts flowing. After a switch-off of the shear field, these additional fluctuations are responsible for the fast decay of stresses, which occurs on a time scale 1/γ˙ . The stress decay after a switch-off in the elastic regime, on the other hand, happens on the much larger time scale of structural relaxation τα. While stresses decrease to zero after a switch-off for temperatures above the glass transition, they decay to a finite value for lower temperatures. The obtained results are important for advancing new theoretical approaches in the framework of mode-coupling theory. Furthermore, they suggest new experimental investigations on colloidal systems. The colloid-polymer mixture is studied in the context of the behaviour near the critical point of phase separation. For the MD simulations a new effective model with soft interaction potentials is introduced and its phase diagram is presented. Here, mainly the equilibrium properties of this model are characterised. While the self-diffusion constants of colloids and polymers do not change strongly when the critical point is approached, critical slowing down of interdiffusion is observed. The order parameter fluctuations can be determined through the long-wavelength limit of static structure factors. For this strongly asymmetric mixture it is shown how the relevant structure factor can be extracted by a diagonalisation of a matrix that contains the partial static structure factors. By presenting first results of this model under shear it is demonstrated that it is suitable for non-equilibrium simulations as well.

Zusammenhang zwischen der Einzelkettendynamik und der Dynamik von Konzentrationsfluktuationen in mehrkomponentigen Polymersystemen

Zusammmenfassung:Um Phasenseparation in binären Polymermischungen zuuntersuchen, werden zwei dynamische Erweiterungen der selbstkonsistenten Feldtheorie (SCFT)entwickelt. Die erste Methode benutzt eine zeitliche Entwicklung der Dichten und wird dynamische selbstkonsistente Feldtheorie (DSCFT) genannt, während die zweite Methode die zeitliche Propagation der effektiven äußeren Felder der SCFT ausnutzt. Diese Methode wird mit External Potential Dynamics (EPD) bezeichnet. Für DSCFT werden kinetische Koeffizienten verwendet, die entweder die lokale Dynamik von Punktteilchen oder die nichtlokale Dynamik von Rouse'schen Polymeren nachbilden. Die EPD-Methode erzeugt mit einem konstanten kinetischen Koeffizienten die Dynamik von Rouse'schen Ketten und benötigt weniger Rechenzeit als DSCFT. Diese Methoden werden für verschiedene Systeme angewendet.Zuerst wird spinodale Entmischung im Volumen untersucht,wobei der Unterschied zwischen lokaler und nichtlokalerDynamik im Mittelpunkt steht. Um die Gültigkeit derErgebnisse zu überprüfen, werden Monte-Carlo-Simulationen durchgeführt. In Polymermischungen, die von zwei Wänden, die beide die gleiche Sorte Polymere bevorzugen, eingeschränkt werden, wird die Bildung von Anreicherungsschichten an den Wänden untersucht. Für dünne Polymerfilme zwischen antisymmetrischen Wänden, d.h. jede Wand bevorzugt eine andere Polymerspezies, wird die Spannung einer parallel zu den Wänden gebildeten Grenzfläche analysiert und der Phasenübergang von einer anfänglich homogenen Mischung zur lokalisierten Phase betrachtet. Des Weiteren wird die Dynamik von Kapillarwellenmoden untersucht.

Grenzflächenspannung und Morphologie fließender zweiphasiger Polymermischungen

In der vorliegenden Arbeit wurde die Morphologie von zweiphasigen Polymermischungen unter Scherung in situ mit Hilfe einer Kombination aus optischer Scherzelle, Durchlichtmikroskop und computergestützten CCD-Kamera untersucht. Als Modellblends dienten die unverträglichen, bei Raumtemperatur flüssigen Polymersysteme Polyisobutylen (PIB)/Polydimethylsiloxan (PDMS) (I) und Poly(dimethyl-co-methylphenyl)siloxan/PDMS (II). Alle Komponenten verhalten sich bei den verwendeten Scherraten newtonisch.Eine der wichtigsten Einflussgrößen für die Blendmorphologie ist die Grenzflächenspannung gamma 12. Sie wurde für I und II mit Hilfe der Methode der Tropfenrelaxation (dynamisch) als Funktion der Zeit bestimmt. Diese Methode erlaubt die Messung von gamma 12 für Tropfen der Phase A in B sowie von Tropfen B in A. Bei der Methode des hängenden Tropfens (statisch) muss der Tropfen aus der Phase mit der höheren Dichte bestehen. Wo der Vergleich der beiden Methoden möglich ist, stimmen die Ergebnisse für beide Systeme sehr gut überein. Bei II sind die aus der Tropfenrelaxation erhaltenen gamma 12-Werte der beiden komplementären Zusammensetzungen im Rahmen des Fehlers gleich, bei I zeigt ein PIB-Tropfen in PDMS einen um 40 % niedrigeren Wert als ein PDMS-Tropfen in PIB, dies wird auf die Diffusion von kurzkettigen Anteilen des PDMS in die Grenzschicht zurückgeführt. Die Grenzflächenspannung hängt also unter Umständen auch bei binären Systemen deutlich von der Zusammensetzung ab.Für II wurde die Blendmorphologie über den gesamten Zusammensetzungsbereich untersucht. Die häufig beobachteten cokontinuierlichen Strukturen treten bei keiner Zusammensetzung auf. Die Phaseninversion erkennt man in einer sprunghaften Änderung der Tropfengröße zwischen phiPDMS <= 0,400 und 0,500; zudem lässt sich die Zeitabhängigkeit der Radien durch Auftragung gegen das Produkt aus der Deformation und dem Quadrat des Volumenbruchs der Tropfenphase für 0 <= phiPDMS <= 0,400 sowie 0,500 <= phiPDMS <= 1 normieren. Für I und II wurde die Morphologieentwicklung bei 25 °C nach Vorscherung bei 100 bzw. 50 s-1 und anschließendem Sprung der Scherrate auf deutlich niedrigere Werte als Funktion der Zeit verfolgt. Hierbei erhält man bei genügend langer Messdauer (mindestens 200 000-300 000 Schereinheiten) konstante Tropfengrößen. Zum einen handelt es sich dabei um pseudo-stationäre Werte, die nur durch Koaleszenz bestimmt sind, zum anderen um echte stationäre Radien, die durch gleichzeitig ablaufende Koaleszenz und Zerteilung entstehen. Für I liegen die stationären Mittelwerte auf der Zerteilungskurve, für II hingegen auf der Koaleszenzkurve.Der Einfluss eines grenzflächenwirksamen Additivs wurde anhand von I durch Zugabe des Blockcopolymer PIB-b-PDMS zu PIB untersucht. Der Vergleich des zeitlichen Verlaufs von gamma 12 mit der Morphologieentwicklung zeigt, dass das Additiv eine Stabilisierung der feinen Tropfen/Matrix-Struktur des Blends durch Hinderung der Koaleszenz und nicht durch Reduktion der Grenzflächenspannung bewirkt.

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

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

Eigenschaftskorrelationen in kolloidalen Festkörpern und Fluiden aus optischen Experimenten

Die Arbeit beschäftigt sich mit ein- und zweikomponentigen, geladenen Kolloidsystemen, die in vollentsalzten wässrigen und organischen Dispersionen kristalline Strukturen ausbilden. Im ersten Teil der Arbeit wird die Wechselwirkung der Kolloide mit verschiedenen Methoden charakterisiert. Dabei zeigten sich quantitative Übereinstimmungen zwischen den Resultaten aus Zellenmodellrechnungen und aus elektrokinetischen Messungen einerseits und Messungen des Phasenverhaltens und der Elastizität andererseits. Diese nunmehr gut gesicherten Diskrepanzen und Korrelationen bedürfen des theoretischen Verständnisses. Im zweiten Teil der Arbeit wurde das Erstarrungsverhalten kolloidaler Scherschmelzen in den kristallinen Zustand mit (zeitaufgelöster) statischer Lichtstreuung und mikroskopischen Methoden untersucht. Dies erlaubte zunächst die kritische Überprüfung klassischer Modelle zur Kristallisationskinetik (Wilson- Frenkel- Gesetz, klassische Nukleationstheorie, Kolmogorov- Johnson- Mehl- Avrami (KJMA)- Modell). Es zeigte sich, dass diese Modelle gut geeignet sind die Verfestigung auch geladener kolloidaler Schmelzen zu beschreiben, wenn die diffusive Einteilchendynamik korrekt berücksichtigt wird. Erstmals wurden Oberflächenspannungen zwischen Kristallkeim und Schmelze für geladene Systeme bestimmt, die im Gegensatz zu Hartkugel- Systemen eine lineare Zunahme mit der Partikelkonzentration aufweisen. Der Methodenpark und die Auswerteverfahren wurden sodann auf binäre kolloidale Mischungen übertragen. Entsprechend den Einzelkomponenten kristallisieren alle Mischungen in einer kubischen Struktur. Leitfähigkeitsmessungen und Elastizität stehen meist im Einklang mit der Nukleation zufallsgeordneter Substitutionskristalle. Für mehrere Proben mit unterschiedlichen Größenverhältnissen wurde mit statischer Lichtstreuung der Einfluss der Komposition und der Partikelkonzentration auf das Nukleationsverhalten untersucht. Im Allgemeinen wurde das Nukleationsszenario einkomponentiger Systeme mit einigen unerwarteten, quantitativen Unterschieden reproduziert. Für eine Probe, die eine Kompositionsordnung andeutet, wurden interessante Korrelationen zwischen der Nukleationskinetik und den Eigenschaften des resultierenden Festkörpers gefunden.

Mixtures of colloidal rods and spheres in bulk and in confinement

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

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

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

Untersuchung der Adhäsionskräfte in interaktiven Pulvermischungen zur Inhalation und deren Veränderung in Abhängigkeit von der Dichtigkeit der Verpackung

Die Untersuchung der Adhäsionskräfte mit Colloid Probe Technik, einer Weiterentwicklung der Rasterkraftmikroskopie (Atomic Force Microscopy=AFM), an erzeugten Carrier- und Wirkstoffkristallen bei Laborbedingungen und unter Einfluss der Luftfeuchte zeigte, dass die Adhäsion von Tiotropiumbromid Monohydrat an Mannitol deutlich höher ist als an Lactose Monohydrat. Die Kohäsionskräfte des Wirkstoffes sind stärker als die Adhäsionskräfte an Carriermaterialien. Auf dieser Grundlage wurde die Hypothese aufgestellt, dass eine Mischung mit Mannitol als Carrier eine kleinere Feinpartikeldosis liefert als eine Mischung mit Lactose. Diese Theorie wurde an interaktiven Pulvermischungen unter Variation von verschiedenen Einflussfaktoren überprüft. Die binare und ternäre Lactose-basierte Mischung lieferte unabhängig vom Kapselmaterial (Gelatine- und Polyethylenkapsel) eine höhere Feinpartikeldosis als die entsprechenden Mannitol-basierten Formulierungen. Die ternäre Komponente bewirkte nur bei Mannitol-basierten Mischungen eine Verbesserung der Feinpartikeldosis. Die detaillierte Untersuchung der aerodynamischen Verteilung ternärer Mischungen zeigte, dass das Kapselmaterial nur unter dem Einfluss der Luftfeuchte und Permeabilität der Blisterverpackung die interpartikulären Wechselwirkungen beeinflusst. Mischungen mit Mannitol als Carrier lieferten unabhängig vom Kapselmaterial, von Luftfeuchte/Lagerungsbedingungen und Permeabilität der Blisterverpackung eine kleinere Feinpartikeldosis als Mischungen mit Lactose als Carrier. Die Carrierart, die Permeabilität der Blisterverpackung und die Luftfeuchte wurden als Haupteinflussfaktoren auf die aerodynamischen Eigenschaften identifiziert. Es konnte gezeigt werden, dass AFM einen wertvollen Beitrag zum Verständnis der interpartikulären Wechselwirkungen leistet und aufgrund prädiktiver Eigenschaften hilfreich in der Entwicklung inhalativer Darreichungs-formen sein kann.