Thermal fluctuations in a lamellar phase of a binary amphiphile–solvent mixture: A molecular-dynamics study

Analytical model of critical buckling transition for smectic liquid crystal based on the viscoelastic scaling of coarse-grained molecular dynamics

The buckling transition of smectic liquid crystals (LCs) is important not only as fundamental physics but also for the rational design of devices to make use of their optical and mechanical properties. However, there exists a huge gap between the specific knowledge and universal analytical formulation. We have conducted coarse-grained molecular dynamics (CGMD) simulations with the force field optimized for the description of buckling phenomena including topological defects to link the molecular nature and continuum formulation. The simulations reveal the viscoelastic characteristics where the critical strain and the compression modulus highly depend on the strain rate as well as the number of layers. Therefore, we formulate the scaling model whose coupling constants depend on both strain rate and domain size. The model reproduces the CGMD results as well as experimental and theoretical values in existing literature. Furthermore, we elucidate from this model that the critical buckling behavior is determined by the competition between the suppression of compression-induced flow and the undulation fluctuation of layers. The framework consisting of the CGMD simulation and the scaling model enables us to estimate the buckling characteristics of smectic LCs reflecting their molecular structures in a wide range from the low-frequency regime that can be verified by experiments to the high-frequency regime beyond the reach of it.

FLUCTUATIONS IN LIPID BILAYERS: ARE THEY UNDERSTOOD?

We review recent computer simulation studies of undulating lipid bilayers. Theoretical interpretations of such fluctuating membranes are most commonly based on generalized Helfrich-type elastic models, with additional contributions of local "protrusions" and/or density fluctuations. Such models provide an excellent basis for describing the fluctuations of tensionless bilayers in the fluid Lαphase at a quantitative level.

However, this description is found to fail for membranes in the gel phase and for membranes subject to high tensions. The fluctuations of tilted gel membranes (Lβ′phase) show a signature of the modulated ripple structure Pβ′, which is a nearby phase observed in the pretransition regime between the Lαand Lβ′state. This complicates a quantitative analysis on mesoscopic length scales. In the case of fluid membranes under tension, the large-wavelength fluctuation modes are found to be significantly softer than predicted by theory.

In the latter context, we also address the general problem of the relation between frame tension and the fluctuation tension, which has been discussed somewhat controversially in recent years. Simulations of very simple model membranes with fixed area show that the fluctuations should be controlled by the frame tension, and not by the internal tension.

The intrinsic interfacial structure of ionic surfactant monolayers at water–oil and water–vapour interfaces

Using computer simulations, we investigate the interfacial structure of sodium dodecyl sulphate (SDS) monolayers adsorbed at the water surface and water–oil interfaces. Using an algorithm that removes the averaging effect of the capillary waves, we obtain a detailed view of the solvation structure of water around the monolayer. We investigate surface concentrations between 45 and 33 Å 2 per surfactant, which are near experimental conditions corresponding to the critical micellar concentration and the formation of Newton black films. The surfactants induce a layering structure in water, which disappears at approximately 1 nm from the monolayer plane. The water molecules exhibit a preferred orientation with the dipoles pointing towards the monolayer. The orientational order decays slowly, but it does not influence the hydrogen bond structure of water, which is significantly disrupted in the interfacial region only. These structural changes are qualitatively the same in SDS–water and oil–SDS–water interfaces. In the latter case, we find a small degree of penetration of oil in the monolayer (between 0.2 and 0.25 molecules per SDS). This small penetration has a measurable effect on the monolayer, which increases its thickness by approximately 10 per cent. The bending modulus of the SDS monolayers is of the order of the thermal energy, k B T .

Coarse-grained simulations of membranes under tension

We investigate the properties of membranes under tension by Monte Carlo simulations of a generic coarse-grained model for lipid bilayers. We give a comprising overview of the behavior of several membrane characteristics, such as the area per lipid, the monolayer overlap, the nematic order, and pressure profiles. Both the low-temperature regime, where the membranes are in a gel L(beta(')) phase, and the high-temperature regime, where they are in the fluid L(alpha) phase, are considered. In the L(beta(')) state, the membrane is hardly influenced by tension. In the fluid state, high tensions lead to structural changes in the membrane, which result in different compressibility regimes. The ripple state P(beta(')), which is found at tension zero in the transition regime between L(alpha) and L(beta(')), disappears under tension and gives way to an interdigitated phase. We also study the membrane fluctuations in the fluid phase. In the low-tension regime the data can be fitted nicely to a suitably extended elastic theory. At higher tensions the elastic fit consistently underestimates the strength of long-wavelength fluctuations. Finally, we investigate the influence of tension on the effective interaction between simple transmembrane inclusions and show that tension can be used to tune the hydrophobic mismatch interaction between membrane proteins.

Membrane-protein interactions in a generic coarse-grained model for lipid bilayers

We study membrane-protein interactions and membrane-mediated protein-protein interactions by Monte Carlo simulations of a generic coarse-grained model for lipid bilayers with cylindrical hydrophobic inclusions. The strength of the hydrophobic force and the hydrophobic thickness of the proteins are systematically varied. The results are compared with analytical predictions of two popular analytical theories: The Landau-de Gennes theory and the elastic theory. The elastic theory provides an excellent description of the fluctuation spectra of pure membranes and successfully reproduces the deformation profiles of membranes around single proteins. However, its prediction for the potential of mean force between proteins is not compatible with the simulation data for large distances. The simulations show that the lipid-mediated interactions are governed by five competing factors: direct interactions; lipid-induced depletion interactions; lipid bridging; lipid packing; and a smooth long-range contribution. The mechanisms leading to hydrophobic mismatch interactions are critically analyzed.

Pores in bilayer membranes of amphiphilic molecules: coarse-grained molecular dynamics simulations compared with simple mesoscopic models

We investigate pores in fluid membranes by molecular dynamics simulations of an amphiphile-solvent mixture, using a molecular coarse-grained model. The amphiphilic membranes self-assemble into a lamellar stack of amphiphilic bilayers separated by solvent layers. We focus on the particular case of tensionless membranes, in which pores spontaneously appear because of thermal fluctuations. Their spatial distribution is similar to that of a random set of repulsive hard disks. The size and shape distribution of individual pores can be described satisfactorily by a simple mesoscopic model, which accounts only for a pore independent core energy and a line tension penalty at the pore edges. In particular, the pores are not circular: their shapes are fractal and have the same characteristics as those of two-dimensional ring polymers. Finally, we study the size-fluctuation dynamics of the pores, and compare the time evolution of their contour length to a random walk in a linear potential.

Shear-induced parallel-to-perpendicular orientation transition in the amphiphilic lamellar phase: A nonequilibrium molecular-dynamics simulation study

The present work is devoted to a study of the shear-induced parallel-to-perpendicular orientation transition in the lamellar system by the large-scale nonequilibrium molecular-dynamics (NEMD) simulation. An effective generic model-A2B2 tetramer for amphiphilies is used. The NEMD simulation produces unambiguous evidence that undulation instability along the vorticity direction sets in well above a critical shear rate and grows in magnitude as the shear rate is further increased. At a certain high shear rate, the coherent undulation instability grows so large that defects are nucleated and the global lamellar monodomain breaks into several aligned lamellar domains. Subsequently layers in these domains rotate into the perpendicular orientation with the rotation of chains towards the y direction, merge into a global perpendicular-aligned lamellar monodomain, and organize into a perfect well-aligned perpendicular lamellar phase by the migration and annihilation of edge dislocations and disclinations. The macroscopic observable viscosity as a function of time or shear rate is correlated with the structural response such as the mesoscopic domain morphology and the microscopic chain conformation. The onset of undulation instability concurs with the start-up of shear-thinning behavior. During the orientation transformation at the high shear rate, the complex time-dependent thixotropic behavior is observed. The smaller viscosity in the perpendicular lamellar phase gives an energetic reason for the shear-induced orientation transition.

Phase diagrams of self-assembled mono-tethered nanospheres from molecular simulation and comparison to surfactants

We perform Brownian dynamics simulations on model 3-D systems of mono-tethered nanospheres (TNS) to study the equilibrium morphologies formed by their self-assembly in a selective solvent. We predict that in contrast to flexible amphiphiles the nanospheres are locally ordered and there is an increase in the local order with an increase in concentration or relative nanoparticle diameter. We present the temperature vs concentration phase diagram for a system of TNS and propose a dimensionless scaling factor F(v) (headgroup volume/tether volume) that allows a comparison between the morphologies formed from TNS and traditional surfactants.

Implicit solvent simulation models for biomembranes

Fully atomic simulation strategies are infeasible for the study of many processes of interest to membrane biology, biophysics and biochemistry. We review various coarse-grained simulation methodologies with special emphasis on methods and models that do not require the explicit simulation of water. Examples from our own research demonstrate that such models have potential for simulating a variety of biologically relevant phenomena at the membrane surface.

Formation and structure of the microemulsion phase in two-dimensional ternary AB+A+B polymeric emulsions

We present an analysis of the structure of the fluctuation-induced microemulsion phase in a ternary blend of balanced AB diblock copolymers with equal amounts of A and B homopolymers. To this end, graphical analysis methods are employed to characterize two-dimensional configuration snapshots obtained with the recently introduced field-theoretic Monte Carlo method. We find that a microemulsion forms when the mean curvature diameter of the lamellar phase coincides roughly with the periodicity of the lamellar phase. Further, we provide evidence to the effect of a subclassification of the microemulsion into a genuine and a defect-driven region.