Molecular dynamics simulations of phase separation in chemically reactive binary mixtures

Physics of active emulsions

Phase separating systems that are maintained away from thermodynamic equilibrium via molecular processes represent a class of active systems, which we call active emulsions. These systems are driven by external energy input, for example provided by an external fuel reservoir. The external energy input gives rise to novel phenomena that are not present in passive systems. For instance, concentration gradients can spatially organise emulsions and cause novel droplet size distributions. Another example are active droplets that are subject to chemical reactions such that their nucleation and size can be controlled, and they can divide spontaneously. In this review, we discuss the physics of phase separation and emulsions and show how the concepts that govern such phenomena can be extended to capture the physics of active emulsions. This physics is relevant to the spatial organisation of the biochemistry in living cells, for the development of novel applications in chemical engineering and models for the origin of life.

Molecular dynamics study of phase separation in fluids with chemical reactions

We present results from the first d=3 molecular dynamics (MD) study of phase-separating fluid mixtures (AB) with simple chemical reactions (A⇌B). We focus on the case where the rates of forward and backward reactions are equal. The chemical reactions compete with segregation, and the coarsening system settles into a steady-state mesoscale morphology. However, hydrodynamic effects destroy the lamellar morphology which characterizes the diffusive case. This has important consequences for the phase-separating structure, which we study in detail. In particular, the equilibrium length scale (ℓ(eq)) in the steady state suggests a power-law dependence on the reaction rate ε:ℓ(eq)∼ε(-θ) with θ≃1.0.

Suppression of Ostwald ripening in active emulsions

Emulsions consisting of droplets immersed in a fluid are typically unstable since they coarsen over time. One important coarsening process is Ostwald ripening, which is driven by the surface tension of the droplets. Stability of emulsions is relevant not only in complex fluids but also in biological cells, which contain liquidlike compartments, e.g., germ granules, Cajal bodies, and centrosomes. Such cellular systems are driven away from equilibrium, e.g., by chemical reactions, and thus can be called active emulsions. In this paper, we study such active emulsions by developing a coarse-grained description of the droplet dynamics, which we analyze for two different chemical reaction schemes. We first consider the simple case of first-order reactions, which leads to stable, monodisperse emulsions in which Ostwald ripening is suppressed within a range of chemical reaction rates. We then consider autocatalytic droplets, which catalyze the production of their own droplet material. Spontaneous nucleation of autocatalytic droplets is strongly suppressed and their emulsions are typically unstable. We show that autocatalytic droplets can be nucleated reliably and their emulsions stabilized by the help of chemically active cores, which catalyze the production of droplet material. In summary, different reaction schemes and catalytic cores can be used to stabilize emulsions and to control their properties.

Influence of extended dynamics on phase transitions in a driven lattice gas

Monte Carlo simulations and dynamical mean-field approximations are performed to study the phase transition in a driven lattice gas with nearest-neighbor exclusion on a square lattice. A slight extension of the microscopic dynamics with allowing the next-nearest-neighbor hops results in dramatic changes. Instead of the phase separation into high- and low-density regions in the stationary state the system exhibits a continuous transition belonging to the Ising universality class for any driving. The relevant features of the phase diagram are reproduced by an improved mean-field analysis.

Molecular dynamics simulations of isomerization kinetics in condensed fluids

Simulations of different reaction schemes for isomerization kinetics in a condensed fluid mixture between two species with small differences in the pair energies show that one of the species dominates at late reaction times. The isomerization is performed on the basis of the energy of the two states, either by choosing minimum energy or by use of Boltzmann weighted kinetics. Both kinetics are autocatalytic and establish domain decomposition with critical fluctuations, which ensure the symmetry break. The model(s) offers a possible explanation of the origin of biomolecular chirality.

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.

Stationary state in a two-temperature model with competing dynamics

A two-dimensional half-filled lattice gas model with nearest-neighbor attractive interaction is studied where particles are coupled to two thermal baths at different temperatures T1 and T2. The hopping of particles is governed by the heat bath at temperature T1 with probability p and the other heat bath (T2) with probability 1-p independently of the hopping direction. On a square lattice the vertical and horizontal interfaces become unstable while interfaces are stable in the diagonal directions. As a consequence, particles condense into a tilted square in the novel ordered state. The p dependence of the resulting nonequilibrium stationary state is studied by a Monte Carlo simulation and dynamical mean-field approximation as well.

Modes selection in polymer mixtures undergoing phase separation by photochemical reactions

Phase separation kinetics and morphology of binary polymer mixtures (A/B) in the presence of photochemical reactions were investigated by using phase-contrast optical microscopy combined with digital image analysis. The polymers were chemically designed in such a way that two types of chemical reactions, intermolecular photodimerization and intramolecular photoisomerization, of polymer segments can be induced and controled by irradiation with ultraviolet light. Unlike the conventional case, the phase separation in the presence of these reactions is spontaneously frozen due to the suppression of the long-wavelength instabilities, resulting in stationary spatial structures with intrinsic periodicities. These characteristic length scales are determined by the competition between the two antagonistic interactions: phase separation as a relatively short-range activation and the photochemical reaction as a long-range inhibition. Furthermore, it was found that the spatial symmetry breaking of concentration fluctuations can emerge from the elastic stress associated with the nonhomogeneous kinetics of the reactions. Experimental data obtained with three types of reactions: A-A only cross-link, A-A and B-B simultaneous cross-links and the reversible A<-->B photoisomerization are described. These results do not only indicate that combination of chemical reactions and phase separation could provide a novel method to control the morphology of multiphase polymer materials, but also suggest that photoreactive polymers can be used as a chemical system to study the mode-selection process in polymers far from thermodynamic equilibrium. (c) 1999 American Institute of Physics.

Steady-state compartmentalization of lipid membranes by active proteins

Using a simple microscopic model of lipid-protein interactions, based on the hydrophobic matching principle, we study some generic aspects of lipid-membrane compartmentalization controlled by a dispersion of active integral membrane proteins. The activity of the proteins is simulated by conformational excitations governed by an external drive, and the deexcitation is controlled by interaction of the protein with its lipid surroundings. In response to the flux of energy into the proteins from the environment and the subsequent dissipation of energy into the lipid bilayer, the lipid-protein assembly reorganizes into a steady-state structure with a typical length scale determined by the strength of the external drive. In the specific case of a mixed dimyristoylphosphatidylcholine-distearoylphosphatidylcholine bilayer in the gel-fluid coexistence region, it is shown explicitly by computer simulation that the activity of an integral membrane protein can lead to a compartmentalization of the lipid-bilayer membrane. The compartmentalization is related to the dynamical process of phase separation and lipid domain formation.