Spinodal decomposition in binary fluids: Effects of hydrodynamic interactions

Segregation of fluids with polymer additives at domain interfaces: a dissipative particle dynamics study

This paper investigates the phase separation kinetics of ternary fluid mixtures composed of a polymeric component (C) and two simple fluids (A and B) using dissipative particle dynamics simulations with a system dimensionality of d = 3. We model the affinities between the components to enable the settling of the polymeric component at the interface of fluids A and B. Thus, the system evolves to form polymer coated morphologies, enabling alteration of the fluids' interfacial properties. This manipulation can be utilized across various disciplines, such as the stabilization of emulsions and foams, rheological control, biomimetic design, and surface modification. We probe the effects of various parameters, such as the polymeric concentration, chain stiffness, and length, on the phase separation kinetics of the system. The simulation results show that changes in the concentration of flexible polymers exhibit perfect dynamic scaling for coated morphologies. The growth rate decreases as the polymeric composition is increased due to reduced surface tension and restricted connectivity between A- and B-rich clusters. Variations in the polymer chain rigidity at fixed composition ratios and degrees of polymerization slow the evolution kinetics of AB fluids marginally, although the effect is more pronounced for perfectly rigid chains. Whereas flexible polymer chain lengths at fixed composition ratios slow down the segregation kinetics of AB fluids slightly, varying the chain lengths of perfectly rigid polymers leads to a significant deviation in the length scale and dynamic scaling for the evolved coated morphologies. The characteristic length scale follows a power-law growth with a growth exponent ϕ that shows a crossover from the viscous to the inertial hydrodynamic regime, where the values of ϕ depend on the constraints imposed on the system.

Designing the Morphology of Separated Phases in Multicomponent Liquid Mixtures

Phase separation of multicomponent liquid mixtures plays an integral part in many processes ranging from industry to cellular biology. In many cases the morphology of coexisting phases is crucially linked to the function of the separated mixture, yet it is unclear what determines the morphology when multiple phases are present. We developed a graph theory approach to predict the topology of coexisting phases from a given set of surface energies, enumerate all topologically distinct morphologies, and reverse engineer conditions for surface energies that produce the target morphology.

Phase behavior and morphology of multicomponent liquid mixtures

Multicomponent systems are ubiquitous in nature and industry. While the physics of few-component liquid mixtures (i.e., binary and ternary ones) is well-understood and routinely taught in undergraduate courses, the thermodynamic and kinetic properties of N-component mixtures with N > 3 have remained relatively unexplored. An example of such a mixture is provided by the intracellular fluid, in which protein-rich droplets phase separate into distinct membraneless organelles. In this work, we investigate equilibrium phase behavior and morphology of N-component liquid mixtures within the Flory-Huggins theory of regular solutions. In order to determine the number of coexisting phases and their compositions, we developed a new algorithm for constructing complete phase diagrams, based on numerical convexification of the discretized free energy landscape. Together with a Cahn-Hilliard approach for kinetics, we employ this method to study mixtures with N = 4 and 5 components. We report on both the coarsening behavior of such systems, as well as the resulting morphologies in three spatial dimensions. We discuss how the number of coexisting phases and their compositions can be extracted with Principal Component Analysis (PCA) and K-means clustering algorithms. Finally, we discuss how one can reverse engineer the interaction parameters and volume fractions of components in order to achieve a range of desired packing structures, such as nested "Russian dolls" and encapsulated Janus droplets.

Distinctive phase separation dynamics of polymer blends: roles of Janus nanoparticles

The present work demonstrates that Janus nanoparticles uniquely promote the phase separation of polymer blends at the early stage of spinodal decomposition, but impede it at the late stage.

Role of a polymeric component in the phase separation of ternary fluid mixtures: a dissipative particle dynamics study

We present the results from dissipative particle dynamics (DPD) simulations of phase separation dynamics in ternary (ABC) fluids mixture in d = 3 where components A and B represent the simple fluids, and component C represents a polymeric fluid. Here, we study the role of polymeric fluid (C) on domain morphology by varying composition ratio, polymer chain length, and polymer stiffness. We observe that the system under consideration lies in the same dynamical universality class as a simple ternary fluids mixture. However, the scaling functions depend upon the parameters mentioned above as they change the time scale of the evolution morphologies. In all cases, the characteristic domain size follows l(t) ∼ tφ with dynamic growth exponent φ, showing a crossover from the viscous hydrodynamic regime (φ = 1) to the inertial hydrodynamic regime (φ = 2/3) in the system at late times.

Effect of bond-disorder on the phase-separation kinetics of binary mixtures: A Monte Carlo simulation study

We present Monte Carlo (MC) simulation studies of phase separation in binary (AB) mixtures with bond-disorder that is introduced in two different ways: (i) at randomly selected lattice sites and (ii) at regularly selected sites. The Ising model with spin exchange (Kawasaki) dynamics represents the segregation kinetics in conserved binary mixtures. We find that the dynamical scaling changes significantly by varying the number of disordered sites in the case where bond-disorder is introduced at the randomly selected sites. On the other hand, when we introduce the bond-disorder in a regular fashion, the system follows the dynamical scaling for the modest number of disordered sites. For a higher number of disordered sites, the evolution morphology illustrates a lamellar pattern formation. Our MC results are consistent with the Lifshitz-Slyozov power-law growth in all the cases.

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.

RNA transcription modulates phase transition-driven nuclear body assembly

Nuclear bodies are RNA and protein-rich, membraneless organelles that play important roles in gene regulation. The largest and most well-known nuclear body is the nucleolus, an organelle whose primary function in ribosome biogenesis makes it key for cell growth and size homeostasis. The nucleolus and other nuclear bodies behave like liquid-phase droplets and appear to condense from the nucleoplasm by concentration-dependent phase separation. However, nucleoli actively consume chemical energy, and it is unclear how such nonequilibrium activity might impact classical liquid-liquid phase separation. Here, we combine in vivo and in vitro experiments with theory and simulation to characterize the assembly and disassembly dynamics of nucleoli in early Caenorhabditis elegans embryos. In addition to classical nucleoli that assemble at the transcriptionally active nucleolar organizing regions, we observe dozens of "extranucleolar droplets" (ENDs) that condense in the nucleoplasm in a transcription-independent manner. We show that growth of nucleoli and ENDs is consistent with a first-order phase transition in which late-stage coarsening dynamics are mediated by Brownian coalescence and, to a lesser degree, Ostwald ripening. By manipulating C. elegans cell size, we change nucleolar component concentration and confirm several key model predictions. Our results show that rRNA transcription and other nonequilibrium biological activity can modulate the effective thermodynamic parameters governing nucleolar and END assembly, but do not appear to fundamentally alter the passive phase separation mechanism.

Phase separation in ternary fluid mixtures: a molecular dynamics study

We present detailed results from molecular dynamics (MD) simulations of phase separation in ternary (ABC) fluid mixtures for d = 2 and d = 3 systems. Our MD simulations naturally incorporate hydrodynamic effects. The domain growth law is l(t) ∼ t(ϕ) with dynamic growth exponent ϕ. Our data clearly indicate that a ternary fluid mixture reaches a dynamical scaling regime at late times with a gradual crossover from ϕ = 1/3 → 1/2 → 2/3 in d = 2 and ϕ = 1/3 → 1 in d = 3 resulting from the hydrodynamic effect in the system. These MD simulations do not yet access the inertial hydrodynamic regime (with l(t) ∼ t(2/3)) of phase separation in ternary fluid mixtures in d = 3.

Modeling multiphase flow using fluctuating hydrodynamics

Fluctuating hydrodynamics provides a model for fluids at mesoscopic scales where thermal fluctuations can have a significant impact on the behavior of the system. Here we investigate a model for fluctuating hydrodynamics of a single-component, multiphase flow in the neighborhood of the critical point. The system is modeled using a compressible flow formulation with a van der Waals equation of state, incorporating a Korteweg stress term to treat interfacial tension. We present a numerical algorithm for modeling this system based on an extension of algorithms developed for fluctuating hydrodynamics for ideal fluids. The scheme is validated by comparison of measured structure factors and capillary wave spectra with equilibrium theory. We also present several nonequilibrium examples to illustrate the capability of the algorithm to model multiphase fluid phenomena in a neighborhood of the critical point. These examples include a study of the impact of fluctuations on the spinodal decomposition following a rapid quench, as well as the piston effect in a cavity with supercooled walls. The conclusion in both cases is that thermal fluctuations affect the size and growth of the domains in off-critical quenches.

Dynamic scaling in phase separation kinetics for quasi-two-dimensional membranes

We consider the dynamics of phase separation in lipid bilayer membranes, modeled as flat two-dimensional liquid sheets within a bulk fluid, both in the creeping flow approximation. We present scaling arguments that suggest asymptotic coarsening in these systems is characterized by a length scale R(t) ~ t(1/2) for critical (bicontinuous) phase separation and R(t) ~t(1/3) for off-critical concentrations (droplet morphology). In this limit, the bulk fluid is the primary source of dissipation. We also address these questions with continuum stochastic hydrodynamic simulations. We see evidence of scaling violation in critical phase separation, where isolated circular domains coarsen slower than elongated ones. However, we also find a region of apparent scaling where R(t) ~ t(1/2) is observed. This appears to be due to the competition of thermal and hydrodynamic effects. We argue that the diversity of scaling exponents measured in experiment and prior simulations can in part be attributed to certain measurements lying outside the asymptotic long-length-scale regime, and provide a framework to help understand these results. We also discuss a few simple generalizations to confined membranes and membranes in which inertia is relevant.

Continuum simulations of biomembrane dynamics and the importance of hydrodynamic effects

Traditional particle-based simulation strategies are impractical for the study of lipid bilayers and biological membranes over the longest length and time scales (microns, seconds and longer) relevant to cellular biology. Continuum-based models developed within the frameworks of elasticity theory, fluid dynamics and statistical mechanics provide a framework for studying membrane biophysics over a range of mesoscopic to macroscopic length and time regimes, but the application of such ideas to simulation studies has occurred only relatively recently. We review some of our efforts in this direction with emphasis on the dynamics in model membrane systems. Several examples are presented that highlight the prominent role of hydrodynamics in membrane dynamics and we argue that careful consideration of fluid dynamics is key to understanding membrane biophysics at the cellular scale.

Phase Separation Kinetics of Binary Systems: Effects of Hydrodynamic Interaction And Surfactants

Cell dynamics computer simulation method is used to investigate effects of hydro-dynamic interaction as well as effects of added surfactants on phase separation kinetics. In the former the obtained scattering structure function for 3-dimensional system reproduces the experimental results for polymer blends remarkably well. In the latter we display self-assembling process in two-dimensional system on computer.

Kinetics of phase separation in fluids: a molecular dynamics study

We present results from extensive three-dimensional molecular dynamics (MD) simulations of phase separation kinetics in fluids. A coarse-graining procedure is used to obtain state-of-the-art MD results. We observe an extended period of temporally linear growth in the viscous hydrodynamic regime. The morphological similarity of coarsening in fluids and solids is also quantified. The velocity field is characterized by the presence of monopolelike defects, which yield a generalized Porod tail in the corresponding structure factor.

Dynamic simulations of multicomponent lipid membranes over long length and time scales

We present a stochastic phase-field model for multicomponent lipid bilayers that explicitly accounts for the quasi-two-dimensional hydrodynamic environment unique to a thin fluid membrane immersed in aqueous solution. Dynamics over a wide range of length scales (from nanometers to microns) for durations up to seconds and longer are readily accessed and provide a direct comparison to fluorescence microscopy measurements in ternary lipid-cholesterol mixtures. Simulations of phase separation kinetics agree with experiment and elucidate the importance of hydrodynamics in the coarsening process.

Spinodal decomposition of polymer solutions: a parallelized molecular dynamics simulation

In simulations of phase separation kinetics, large length and time scales are involved due to the mesoscopic size of the polymer coils, and the structure formation on still larger scales of length and time. We apply a coarse-grained model of hexadecane dissolved in supercritical carbon dioxide, for which in previous work the equilibrium phase behavior has been established by Monte Carlo methods. Using parallelized simulations on a multiprocessor supercomputer, large scale molecular dynamics simulations of phase separation following pressure jumps are presented for systems containing N=435136 coarse-grained particles, which correspond to several millions of atoms in a box with linear dimension 447 A . Even for large systems the phase separation can be observed up to the final, macroscopically segregated, equilibrium state. It is shown that in the segregation process the two order parameters of the system (density and concentration) are strongly coupled. The system does not follow the predicted growth law for the characteristic domain size l(t) proportional, variant t in binary fluid mixtures for the range of times accessible in the simulation. Instead, it exhibits a distinctly slower growth, presumably due to the dynamic asymmetry of the constituents.

Spinodal decomposition via surface diffusion in polymer mixtures

We present experimental results for spinodal decomposition in polymer mixtures of gelatin and dextran. The domain growth law is found to be consistent with t 1/4 growth over extended time regimes. Similar results are obtained from lattice simulations of a polymer mixture. This slow growth arises due to the suppression of the bulk mobility of polymers. In that case, spinodal decomposition is driven by the diffusive transport of material along domain interfaces, which gives rise to a t 1/4 growth law.

Spinodal decomposition in thin films: molecular-dynamics simulations of a binary Lennard-Jones fluid mixture

We use molecular dynamics (MD) to simulate an unstable homogeneous mixture of binary fluids (AB), confined in a slit pore of width D. The pore walls are assumed to be flat and structureless and attract one component of the mixture (A) with the same strength. The pairwise interactions between the particles are modeled by the Lennard-Jones potential, with symmetric parameters that lead to a miscibility gap in the bulk. In the thin-film geometry, an interesting interplay occurs between surface enrichment and phase separation. We study the evolution of a mixture with equal amounts of A and B, which is rendered unstable by a temperature quench. We find that A-rich surface enrichment layers form quickly during the early stages of the evolution, causing a depletion of A in the inner regions of the film. These surface-directed concentration profiles propagate from the walls towards the center of the film, resulting in a transient layered structure. This layered state breaks up into a columnar state, which is characterized by the lateral coarsening of cylindrical domains. The qualitative features of this process resemble results from previous studies of diffusive Ginzburg-Landau-type models [S. K. Das, S. Puri, J. Horbach, and K. Binder, Phys. Rev. E 72, 061603 (2005)], but quantitative aspects differ markedly. The relation to spinodal decomposition in a strictly two-dimensional geometry is also discussed.

Three-dimensional lattice-Boltzmann simulations of critical spinodal decomposition in binary immiscible fluids

We use a modified Shan-Chen, noiseless lattice-BGK model for binary immiscible, incompressible, athermal fluids in three dimensions to simulate the coarsening of domains following a deep quench below the spinodal point from a symmetric and homogeneous mixture into a two-phase configuration. The model is derivable from a continuous-time Boltzmann-BGK equation in the presence of an intercomponent body force. We find the average domain size grows with time as t(gamma), where gamma increases in the range 0.545+/-0.014<gamma<0.717+/-0.002, consistent with a crossover between diffusive t(1/3) and hydrodynamic viscous, t(1.0), behavior. We find good collapse onto a single scaling function, yet the domain growth exponents differ from previous results for similar values of the unique characteristic length L0 and time T0 that can be constructed out of the fluid's parameters. This rebuts claims of universality for the dynamical scaling hypothesis. For Re=2.7 and small wave numbers q we also find a q(2)<-->q(4) crossover in the scaled structure function, which disappears when the dynamical scaling reasonably improves at later stages (Re=37). This excludes noise as the cause for a q(2) behavior, as analytically derived from Yeung and proposed by Appert et al. and Love et al. on the basis of their lattice-gas simulations. We also observe exponential temporal growth of the structure function during the initial stages of the dynamics and for wave numbers less than a threshold value, in accordance with the diffusive Cahn-Hilliard Model B. However, this exponential growth is also present in regimes proscribed by that model. There is no evidence that regions of parameter space for which the scheme is numerically stable become unstable as the simulations proceed, in agreement with finite-difference relaxational models and in contradistinction with an unconditionally unstable lattice-BGK free-energy model previously reported. Those numerical instabilities that do arise in this model are the result of large intercomponent forces which turn the equilibrium distribution negative.

Viscoelastic effects in relaxation processes of concentration fluctuations in dynamically asymmetric polymer blends

Relaxation processes of the concentration fluctuations induced by a rapid pressure change were investigated for a dynamically asymmetric polymer blend [deuterated polybutadiene (DPB)/polyisoprene (PI)] with a composition of 50-50 by weight by using time-resolved small-angle neutron scattering. The pressure change was carried out inside the single-phase of the blend with the cell designed for polymeric systems under high pressure and temperature. Time change in the scattered intensity distribution with wave number (q) during the relaxation processes was found to be approximated by Cahn-Hilliard-Cook linearized theory. The theoretical analysis yielded the q dependence of Onsager kinetic coefficient that is characterized by the q(-2) dependence at q(xi)(ve)>1 with the characteristic length xi(ve) (with xi(ve) being the viscoelastic length) being much larger than radius of gyration of DPB or PI. The estimated xi(ve) agrees well with that calculated using the Doi and Onuki theory that takes into account the viscoelastic effects arising from the dynamical asymmetry between the component polymers in the relaxation of concentration fluctuations.

Three-dimensional hydrodynamic lattice-gas simulations of domain growth and self-assembly in binary immiscible and ternary amphiphilic fluids

We simulate the dynamics of phase assembly in binary immiscible fluids and ternary microemulsions using a three-dimensional hydrodynamic lattice-gas approach. For critical spinodal decomposition we perform the scaling analysis in reduced variables introduced by Jury et al. [Phys. Rev. E 59, R2535 (1999)] and by Bladon et al. [Phys. Rev. Lett. 83, 579 (1999)]. We find a late-stage scaling exponent consistent with the R approximately t(2/3) inertial regime. However, as observed with the previous lattice-gas model of Appert et al. [J. Stat. Phys. 81, 181 (1995)] our data do not fall in the same range of reduced length and time as those of Jury et al. and Bladon et al. For off-critical binary spinodal decomposition we observe a reduction of the effective exponent with decreasing volume fraction of the minority phase. However, the n=1 / 3 Lifshitz-Slyzov-Wagner droplet coalescence exponent is not observed. Adding a sufficient number of surfactant particles to a critical quench of binary immiscible fluids produces a ternary bicontinuous microemulsion. We observe a change in scaling behavior from algebraic to logarithmic growth for amphiphilic fluids in which the domain growth is not arrested. For formation of a microemulsion where the domain growth is halted we find that a stretched exponential growth law provides the best fit to the data.

Spinodal decomposition of two-dimensional fluid mixtures: A spectral analysis of droplet growth

The spinodal decomposition of two-dimensional fluid mixture is studied by numerical simulation. For the high viscous fluid mixture it has not been evident whether the interfacial tension is relevant to the droplet growth or not. A length scale R defined by the structure function extracting the effect of the long wavelength mode justifies a rapid growth close to R approximately t, but the length scale energetically defined reveals a much slower growth R approximately t(0.5), where t is time. This discrepancy represents the violation of the dynamical scaling with single length scale. The slow gowth of the length scale is attributed to the accumulation of the number of isolated droplets in phase separating state, whereas the rapid growth represents the relevance of the surface tension as the driving force in two dimensions. For a low viscous fluid mixture the dynamical scaling is a good assumption with the growth law R approximately t(2/3) up to a very large Reynolds number Re approximately 1500, which is the limit in the present simulation.

Two-dimensional model of phase segregation in liquid binary mixtures

The hydrodynamic effects on the late stage kinetics of phase separation in liquid mixtures is studied using the model H. Mass and momentum transport are coupled via a nonequilibrium body force, which is proportional to the Peclet number alpha, i.e., the ratio between convective and diffusive molar fluxes. Numerical simulations based on this theoretical model show that phase separation in low viscosity, liquid binary mixtures is mostly driven by convection, thereby explaining the experimental findings that the process is fast, with the typical size of single-phase domains increasing linearly with time. However, as soon as sharp interfaces form, the linear growth regime reaches an end, and the process appears to be driven by diffusion, although the condition of local equilibrium is not reached. During this stage, the typical size of the nucleating drops increases like t(n), where 1/3< n <1/2, depending on the value of the Peclet number. As the Peclet number increases, the transition between convection- and diffusion-driven regimes occurs at larger times, and therefore for larger sizes of the nucleating drops.

Particle components of dark matter

Particle candidates for astrophysical dark matter are reviewed, with particular emphasis on the lightest supersymmetric particle and the axion. The former is now constrained by accelerator experiments to have a mass above about 40 GeV, and ongoing searches at accelerators, in space, and in underground experiments have a good chance to detect it. A reevaluation of the constraint on the axion from supernova 1987a leaves open an interesting window where it may be detected if it constitutes the galactic halo.