Molecular dynamics simulations of phase separation in the presence of surfactants

Effect of amphiphilic polymers on phase separating binary mixtures: A DPD simulation study

We present the phase separation dynamics of a binary (AB), simple fluid (SF), and amphiphilic polymer (AP) mixture using dissipative particle dynamics simulation at d = 3. We study the effect of different AP topologies, including block copolymers, ring block copolymers (RCP), and miktoarm star polymers, on the evolution morphologies, dynamic scaling functions, and length scale of the AB mixture. Our results demonstrate that the presence of APs leads to significantly different evolution morphologies in SF. However, the deviation from dynamical scaling is prominent, mainly for RCP. Typically, the characteristic length scale for SF follows the power law R(t) ∼ tϕ, where ϕ is the growth exponent. In the presence of high AP, we observe diffusive growth (ϕ → 1/3) at early times, followed by saturation in length scale (ϕ → 0) at late times. The extent of saturation varies with constraints imposed on the APs, such as topology, composition ratio, chain length, and stiffness. At lower composition ratios, the system exhibits inertial hydrodynamic growth (ϕ → 2/3) at asymptotic times without clearly exhibiting the viscous hydrodynamic regime (ϕ → 1) at earlier times in our simulations. Our results firmly establish the existence of hydrodynamic growth regimes in low surfactant-influenced phase separation kinetics of binary fluids and settle the related ambiguity in d = 3 systems.

Using Microemulsions: Formulation Based on Knowledge of Their Mesostructure

Microemulsions, as thermodynamically stable mixtures of oil, water, and surfactant, are known and have been studied for more than 70 years. However, even today there are still quite a number of unclear aspects, and more recent research work has modified and extended our picture. This review gives a short overview of how the understanding of microemulsions has developed, the current view on their properties and structural features, and in particular, how they are related to applications. We also discuss more recent developments regarding nonclassical microemulsions such as surfactant-free (ultraflexible) microemulsions or ones containing uncommon solvents or amphiphiles (like antagonistic salts). These new findings challenge to some extent our previous understanding of microemulsions, which therefore has to be extended to look at the different types of microemulsions in a unified way. In particular, the flexibility of the amphiphilic film is the key property to classify different microemulsion types and their properties in this review. Such a classification of microemulsions requires a thorough determination of their structural properties, and therefore, the experimental methods to determine microemulsion structure and dynamics are reviewed briefly, with a particular emphasis on recent developments in the field of direct imaging by means of electron microscopy. Based on this classification of microemulsions, we then discuss their applications, where the application demands have to be met by the properties of the microemulsion, which in turn are controlled by the flexibility of their amphiphilic interface. Another frequently important aspect for applications is the control of the rheological properties. Normally, microemulsions are low viscous and therefore enhancing viscosity has to be achieved by either having high concentrations (often not wished for) or additives, which do not significantly interfere with the microemulsion. Accordingly, this review gives a comprehensive account of the properties of microemulsions, including most recent developments and bringing them together from a united viewpoint, with an emphasis on how this affects the way of formulating microemulsions for a given application with desired properties.

Structure and dynamics of an aqueous solution containing poly-(acrylic acid) and non-ionic surfactant octaethylene glycol n-decyl ether (C10E8) aggregates and their complexes investigated by molecular dynamics simulations

A detailed molecular dynamics simulation study of the self-assembly, intermolecular structure and thermodynamic behavior of an aqueous solution of non-ionic surfactant octa ethylene glycol n-decyl ether (C₁₀E₈) in the presence of a non-ionic polar polymer poly(acrylic acid) PAA is presented. The aggregation number N_agg and concentration of surfactant C_s in the simulation systems were varied in the range 0.01-0.32 M and 5 < N_agg < 101 (dilute to concentrated) with a dilute polymer concentration (C_p = 0.01 M). Lamellar aggregates of non-ionic surfactant in bulk aqueous solution are shown by molecular level computations for the first time. Spherical micellar aggregates and lamellar aggregates are formed at low and high N_agg, respectively. The transition from the spherical micelle phase to the lamellar phase in a binary solution is captured for the first time. A conformational transition from coiled to extended PAA chains adsorbed on the surfactant aggregate occurs at a particular value of N_agg, commensurate with the transition from spherical micelle aggregates to anisotropic lamellar aggregates. Formation of the surfactant aggregate in binary and ternary solutions and the polymer-surfactant complex in a ternary solution is enthalpically favored. Adsorption of PAA on the surfactant aggregate surface is driven by hydrogen bonds (HBs) between carboxylic acid groups of PAA and ethylene oxide groups of C₁₀E₈. A significant number of HBs occur between polar oxygens of C₁₀E₈ and hydroxyl oxygens of PAA. The results are in agreement with the limited available experimental data on this system.

Polymerization-induced polymer aggregation or polymer aggregation-enhanced polymerization? A computer simulation study

In this study, using dissipative particle dynamics simulations coupled with the stochastic reaction model, we investigate the polymerization-induced polymer aggregation process and the polymer aggregation-enhanced polymerization process in a binary solution.

Numerical simulations of bijel morphology in thin films with complete surface wetting

Bijels are a relatively new class of soft materials that have many potential energy and environmental applications. In this work, simulation results of bijel evolution confined within thin films with preferential surface wetting are presented. The computational approach used is a hybrid Cahn-Hilliard/Brownian dynamics method. In the absence of suspended particles, we demonstrate that the model accurately captures the rich kinetics associated with diffusion-based surface-directed spinodal decomposition, as evidenced by comparison with previous theoretical and simulation-based studies. When chemically neutral particles are included in the films, the simulations capture surface-modified bijel formation, with stabilized domain structures comparable with the experimental observations of Composto and coworkers. Namely, two basic morphologies - bicontinuous or discrete - are seen to emerge, with direct dependence on the film thickness, particle volume fraction, and particle radius.

Microphase separation induced by interfacial segregation of isotropic, spherical nanoparticles

In a recent experiment by Chung et al. [Nano Lett. 5, 1878 (2005)] and simulation by Stratford et al. [Science 309, 2198 (2005)] on immiscible blends containing nanoscale particles, it was shown that the phase separation of the two polymers can be prevented as a result of the aggregation of the nanoparticles at the interfaces between the two polymers. Motivated by these studies, we performed large scale systematic simulations, based on the dissipative particle dynamics approach, on immiscible binary (A-B) fluids containing moderate volume fractions of isotropic nanoscale spherical particles N. The nanoparticles preferentially segregate at the interfaces between the two fluids if the pairwise interactions between the three components are such that chi(AB)>/chi(AN)-chi(BN)/. We find that at later times, the average domain size saturates to a value, L approximately R(N)/phi(N), where R(N) and phi(N) are the radius and volume fraction of the nanoparticles, respectively. For small nanoparticles, however, full phase separation is observed.

Coarsening dynamics of ternary amphiphilic fluids and the self-assembly of the gyroid and sponge mesophases: Lattice-Boltzmann simulations

By means of a three-dimensional amphiphilic lattice-Boltzmann model with short-range interactions for the description of ternary amphiphilic fluids, we study how the phase separation kinetics of a symmetric binary immiscible fluid is altered by the presence of the amphiphilic species. We find that a gradual increase in amphiphile concentration slows down domain growth, initially from algebraic to logarithmic temporal dependence, and, at higher concentrations, from logarithmic to stretched-exponential form. In growth-arrested stretched-exponential regimes, at late times we observe the self-assembly of sponge mesophases and gyroid liquid-crystalline cubic mesophases, hence confirming that (a) amphiphile-amphiphile interactions need not be long-ranged in order for periodically modulated structures to arise in a dynamics of competing interactions, and (b) a chemically specific model of the amphiphile is not required for the self-assembly of cubic mesophases, contradicting claims in the literature. We also observe a structural order-disorder transition between sponge and gyroid phases driven by amphiphile concentration alone or, independently, by the amphiphile-amphiphile and the amphiphile-binary fluid coupling parameters. For the growth-arrested mesophases, we also observe temporal oscillations in the structure function at all length scales; most of the wave numbers show slow decay, and long-term stationarity or growth for the others. We ascribe this behavior to a combination of complex amphiphile dynamics leading to Marangoni flows.

Fast growth in phase-separating A-B-copolymer ternary mixtures with a chemical reaction

We study the dynamics of phase separation of a binary A-B- polymer mixture with copolymer C, which is produced by the reaction of two counterpart reactive polymers A and B at the interface via the chemical reaction A+B right harpoon over left harpoon C. For low interfacial energy between the A and B phases, where the copolymer prefers to locate at interfaces, we show that the chemical reaction accelerates the phase separation of the system dramatically, because the backward reaction always drives the creation of immiscible A and B pairs at interfaces, which speed up the phase separation of the system, while the forward reaction process becomes more and more difficult as the interfaces are gradually saturated by copolymers. We also indicate that for a fixed chemical reaction rate constant, as the initial concentration of the copolymers increases, the domain growth at the late stage is speeded up as a result of the backward chemical reaction. However, when the interfacial energy is high, both forward and backward reactions coexist due to the occurrence of unsaturated interfaces, but the relative strength of reaction rates has no appreciable effect on domain growth during spinodal decomposition, because the interfacial energy dominates phase separation.

Real-coded lattice gas model for ternary amphiphilic fluids

We have developed an amphiphilic surfactant model in the framework of the real-coded lattice gas (RLG), in order to analyze the dynamics and structure of ternary fluids. Formation of both the oil-in-water and the bicontinuous microemulsion phases, as well as the reduction of surface tension by adsorption of surfactant at interface are successfully reproduced in numerical simulations. Our model is simple in terms of description and implementation, however, complex structures and dynamic behavior of the ternary fluids emerge from the collective dynamics of the RLG and surfactant particles.

Lattice-boltzmann model for interacting amphiphilic fluids

We develop our recently proposed lattice-Boltzmann method for the nonequilibrium dynamics of amphiphilic fluids [H. Chen, B. M. Boghosian, P. V. Coveney, and M. Nekovee, Proc. R. Soc. London, Ser. A 456, 2043 (2000)]. Our method maintains an orientational degree of freedom for the amphiphilic species and models fluid interactions at a microscopic level by introducing self-consistent mean-field forces between the particles into the lattice-Boltzmann dynamics, in a way that is consistent with kinetic theory. We present the results of extensive simulations in two dimensions which demonstrate the ability of our model to capture the correct phenomenology of binary and ternary amphiphilic fluids.

Phase transitions and steady-state microstructures in a two-temperature lattice-gas model with mobile active impurities

The nonequilibrium, steady-state phase transitions and the structure of the different phases of a two-dimensional system with two thermodynamic temperatures are studied via a simple lattice-gas model with mobile active impurities ("hot/cold spots") whose activity is controlled by an external drive. The properties of the model are calculated by Monte Carlo computer-simulation techniques. The two temperatures and the external drive on the system lead to a rich phase diagram including regions of microstructured phases in addition to macroscopically ordered (phase-separated) and disordered phases. Depending on the temperatures, microstructured phases of both lamellar and droplet symmetry arise, described by a length scale that is determined by the characteristic temperature controlling the diffusive motion of the active impurities.

Phase separation dynamics and lateral organization of two-component lipid membranes

The non-equilibrium dynamic ordering process of coexisting phases has been studied for two-component lipid bilayers composed of saturated di-acyl phospholipids with different acyl chain lengths, such as DC14PC-DC18PC and DC12PC-DC18PC. By means of a microscopic interaction model and computer-simulation techniques the non-equilibrium properties of these two mixtures have been determined with particular attention paid to the effects of the non-equilibrium ordering process on membrane heterogeneity in terms of local and global lateral membrane organization. The results reveal that a sudden temperature change that takes the lipid mixture from the fluid one-phase region into the gel-fluid phase-coexistence region leads to the formation of a large number of small lipid domains which slowly are growing in time. The growth of the lipid domains, which is limited by long-range diffusion of the lipid molecules within the two-dimensional membrane plane, gives rise to the existence of a highly heterogeneous percolative-like structure with a network of interfacial regions that have properties different from those of the phase-separated gel and fluid bulk phases. The results, which are discussed in relation to recent experimental observations interpreted in terms of a percolative-like membrane structure within the two phase region (Almeida, P.F.F., Vaz, W.L.C., and T.E. Thompson. 1992. Biochemistry 31:7198-7210), suggest that non-equilibrium effects may influence lipid domain formation and membrane organization on various length and time scales. Such effects might be of importance in relation to membrane processes that require molecular mobility of the membrane components in restricted geometrical environments of the compartmentalized lipid membrane.