Numerical Simulations of a Smectic Lamellar of Amphiphilic Molecules

This thesis deals with the liquid crystalline lamellar phase L? built by amphiphilic molecules in aqueous solutions, and its interaction with macromolecules. We perform molecular dynamics simulations to study thermal fluctuations and defects appearing in a stack of parallel amphiphilic bilayers separated by layers of solvent. The idealized, coarse-grained model represents the solvent with soft spheres and the amphiphiles with bead-and-spring tetramers (two solvophilic beads and two solvophobic beads). The algorithm used for this thesis describes the lamellar phase in the isobaric isothermal ensemble without surface tension (N,P, T, ? = 0). First, we verify that the model exhibits a liquid-crystalline lamellar phase, which we characterize. In a second part, we study the elasticity of this smectic lamellar phase. The position fluctuation spectra of the bilayers are computed, and compared to the predictions of the ôDiscrete Harmonicö(DH) theory for the elasticity of smectic phases. The bilayer fluctuations observed in the simulation of a stack of fifteen bilayers are well described by the DH-theory, so that the two elastic constants - the bending rigidity Kc and the smectic compressibility modulus B - can be computed. Then, we investigate the point defects appearing in the smectic because of thermal fluctuations. It turns out that transient pores spontaneously nucleate in the bilayers of the lamellar phase. On the contrary, necks and passages between the bilayers are rarely detected. The size and shape distributions of the pores are investigated. The relationship between their area a and their contour length c is well described by the scaling law a ? c2/3 - the same scaling as two dimensional closed random walks. Additionally, the surface tension is zero. Therefore we consider that the energy of a pore depends explicitly only on the contour length of the pore. The effective free energy of individual pores and the line tension of the pore edge are then estimated from the contour-length distribution. Besides, a time-dependent analysis shows that the displacement of the pores within the bilayers during their life-time is very limited. In the last chapter of this thesis, we investigate a lamellar phase doped by a solventsoluble flexible linear polymer inserted between two bilayers. Two polymer types were simulated: adsorbing or non-adsorbing. In both cases, the interactions between the bilayers are softened in the presence of a polymer. However, the conformations of the chain strongly depend on the interactions between the polymer and the bilayers: An adsorbing polymer remains aligned with the bilayers and confined in the thin solvent layer, whereas a non-adsorbing polymer condenses into a globule. Contrarily to standard hypothesis, a non-adsorbing polymer locally modifies the interlamellar spacing, and triggers the formation of pores in its vicinity. Diese Arbeit behandelt die fl³ssig-kristalline lamellare Phase (die sogenannte L? Phase), die amphiphile Molek³le in wõssriger L÷sung ausbilden. Diese lamellare Phase besteht aus mehrere Lagen paralleler amphiphiler Doppelschichten, die durch L÷sungsmittel voneinander getrennt sind. Wir studieren mittels Molekulardynamiksimulationen die thermischen Fluktuationen der Doppelschichten und die Defekte, die in der lamellaren Phase auftreten k÷nnen, und die Auswirkungen eines Makromolek³ls, das zwischen die Doppelschichten gestetzt wird. In dem zugrundeliegenden, idealisierten ôcoarse-grainedö Modell werden das L÷sungsmittel als weiche Kugeln und die Amphiphile als Tetramere (zwei hydrophile Kugeln und zwei hydrophobe Kugeln) reprõsentiert. Der Algorithmus, der in dieser Arbeit verwendet wird, beschreibt die lamellare Phase im isobaren isothermen Ensemble ohne Oberflõchenspannung (NPT, ? = 0). Zuerst verifizieren wir, da¯ das Modell tatsõchlich eine stabile L? Phase bildet, und charakterisieren ihre fl³ssig-kristalline Struktur. Um die Elastizitõt der lamellaren Phase zu untersuchen, berechnen wir die Fluktuationsspektren der Positionen von den Doppelschichten und vergleichen sie mit den Vorhersagen der ôDiscrete Harmonic Theoryö (DH) f¿ur die Elastizitõt der smektischen Phasen. Die Ergebnisse der Simulation, die mit einem Stapel von f¿unfzehn Doppelschichten durchgef³hrt wurde, stimmen mit der DH Theorie ³berein. Daher k÷nnen die elastischen Konstanten (Steifigkeitsmodul Kc und smektischer Kompressionsmodul B) berechnet werden. Nachfolgend betrachten wir lokale Defekte, die auf Grund thermischer Fluktuationen erscheinen. Kurzlebige Poren treten spontan in den Doppelschichten der lamellaren Phase auf. Im Gegensatz dazu werden Verbindungen oder Durchgõnge zwischen den Doppelschichten selten beobachtet. Die Verteilungen von Gr÷¯e und Form der Poren werden bestimmt. Das Verhõltnis zwischen Flõche der Poren a und ihrer Konturlõnge c wird durch das Skalengesetz a ? c2/3 beschrieben (ein analoges Skalengesetz gilt f³r zweidimensionale geschlossene Irrfahrten). Zusõtzlich ist die Oberflõchenspannung null. Wir nehmen dann an, da¯ die Energie einer einzelnen Pore nur von ihrer Konturlõnge abhõngt. Diese Energie und die Linienspannung des Porenrandes werden mittels der Verteilung f³r die Konturlõnge geschõtzt. Zusõtzlich zeigt eine zeitabhõngige Analyse, da¯ die Poren innerhalb ihrer Lebensdauer nicht nennenswert diffundieren. Schlie¯lich wird eine lamellare Phase untersucht, in der ein hydrophiles, flexibles, lineares Polymer zwischen zwei Doppelschichten eingesetzt wurde. Zwei Typen von Polymeren werden simuliert: adsorbierende und nicht-adsorbierende. In beiden Fõlle werden die Wechselwirkungen zwischen den Doppelschichten in Anwesenheit des Polymer verõndert. Die Konformation des Polymers hõngt allerdings stark von den Wechselwirkungen zwischen dem Polymer und den Amphiphilen ab: Ein adsorbierendes Polymer verbleibt wõhrend der gesamten Simulation in der d³nnen L÷sungsmittelschicht zwischen den beiden Doppelschichten. Ein nicht-adsorbierendes Polymer dagegen kondensiert in einen kompakten Tropfenö. Im Gegensatz zu den ³bliches Annahmen õndert ein nichtadsorbierendes Polymer lokal den interlamellaren Abstand und erzeugt Poren in der Doppelschicht in seiner Nõhe.