Spinodal decomposition in 3-space

Diffusive first-order phase transition: nucleation, growth and coarsening in solids

The phenomena of nucleation and growth, which fall into the category of first-order phase transitions, are of great importance. They are present everywhere in our daily lives. They enable us to understand and model a vast number of phenomena, from the formation of raindrops, to the gelling of polymers, the evolution of a virus population and the formation of galaxies. Surprisingly, this whole range of phenomena can be described by two seemingly antagonistic approaches: classical nucleation theory, which highlights the atomistic approach of the diffusion process, and the phase-field (PF) approach, which erases the discrete nature of the diffusion process. Although there is an huge quantity of articles and review papers dealing with the problem of first-order phase transition, the subject is so important and vast that it is very difficult to provide nowadays exhaustive syntheses on the subject. The revival over the past 20 years in the condensed matter world of PF approaches such as PF crystal, or the recent development of optimization methods such as gentle ascend dynamics, as well as the emergence of atom probe tomography, have enabled us to better understand the links between these antagonistic approaches, and above all to provide new experimental results to test the limits of both. This renewal has motivated the writing of this review, both to take stock of current knowledge on these two approaches. This review has two distinct objectives: summarizing generic previous models applies to discuss the nucleation, the growth and the coarsening processes. Despite some reviews already exist on these different subject, few of them present the different logical links between these models and their limitations, unifying them within the framework of the theory of macroscopic fluctuations, which has been developed over the last 20 years. In particular, we present the extension of the Cahn-Hilliard formalism to model the nucleation and growth process and we discuss the relevance of the notion of pseudo-spinodal and discuss. Such an extension allows interpreting experiments performed fat from the solubility limit and the spinodal line. Finally, this work proposes some clues to make this unified approach more predictive.

Non-equilibrium dynamic hyperuniform states

Disordered hyperuniform structures are an exotic state of matter having suppressed density fluctuations at large length-scale similar to perfect crystals and quasicrystals but without any long range orientational order. In the past decade, an increasing number of non-equilibrium systems were found to have dynamic hyperuniform states, which have emerged as a new research direction coupling both non-equilibrium physics and hyperuniformity. Here we review the recent progress in understanding dynamic hyperuniform states found in various non-equilibrium systems, including the critical hyperuniformity in absorbing phase transitions, non-equilibrium hyperuniform fluids and the hyperuniform structures in phase separating systems via spinodal decomposition.

Pattern dynamics of density and velocity fields in segregation of fluid mixtures

We present comprehensive numerical results from a study of model H, which describes phase separation kinetics in binary fluid mixtures. We study the pattern dynamics of both density and velocity fields in d = 2, 3. The density length scales show three distinct regimes, in accordance with analytical arguments. The velocity length scale shows a diffusive behavior. We also study the scaling behavior of the morphologies for density and velocity fields and observe dynamical scaling in the relevant correlation functions and structure factors. Finally, we study the effect of quenched random field disorder on spinodal decomposition in model H.

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.

Hydrodynamic effects in kinetics of phase separation in binary fluids: Critical versus off-critical compositions

Via hydrodynamics-preserving molecular dynamics simulations we study growth phenomena in a phase-separating symmetric binary mixture model. We quench high-temperature homogeneous configurations to state points inside the miscibility gap, for various mixture compositions. For compositions at the symmetric or critical value we capture the rapid linear viscous hydrodynamic growth due to advective transport of material through tubelike interconnected domains. For state points very close to any of the branches of the coexistence curve, the growth in the system, following nucleation of disconnected droplets of the minority species, occurs via a coalescence mechanism. Using state-of-the-art techniques, we have identified that these droplets, between collisions, exhibit diffusive motion. The value of the exponent for the power-law growth, related to this diffusive coalescence mechanism, has been estimated. While the exponent nicely agrees with that for the growth via the well-known Lifshitz-Slyozov particle diffusion mechanism, the amplitude is stronger. For the intermediate compositions we observe initial rapid growth that matches the expectations for viscous or inertial hydrodynamic pictures. However, at later times these types of growth cross over to the exponent that is decided by the diffusive coalescence mechanism.

Ice needles weave patterns of stones in freezing landscapes

Patterned ground, defined by the segregation of stones in soil according to size, is one of the most strikingly self-organized characteristics of polar and high-alpine landscapes. The presence of such patterns on Mars has been proposed as evidence for the past presence of surface liquid water. Despite their ubiquity, the dearth of quantitative field data on the patterns and their slow dynamics have hindered fundamental understanding of the pattern formation mechanisms. Here, we use laboratory experiments to show that stone transport is strongly dependent on local stone concentration and the height of ice needles, leading effectively to pattern formation driven by needle ice activity. Through numerical simulations, theory, and experiments, we show that the nonlinear amplification of long wavelength instabilities leads to self-similar dynamics that resemble phase separation patterns in binary alloys, characterized by scaling laws and spatial structure formation. Our results illustrate insights to be gained into patterns in landscapes by viewing the pattern formation through the lens of phase separation. Moreover, they may help interpret spatial structures that arise on diverse planetary landscapes, including ground patterns recently examined using the rover Curiosity on Mars.

Pattern formation of skin cancers: Effects of cancer proliferation and hydrodynamic interactions

We study pattern formation of skin cancers by means of numerical simulation of a binary system consisting of cancer and healthy cells. We extend the conventional model H for macrophase separations by considering a logistic growth of cancer cells and also a mechanical friction between dermis and epidermis. Importantly, our model exhibits a microphase separation due to the proliferation of cancer cells. By numerically solving the time evolution equations of the cancer composition and its velocity, we show that the phase separation kinetics strongly depends on the cell proliferation rate as well as on the strength of hydrodynamic interactions. A steady-state diagram of cancer patterns is established in terms of these two dynamical parameters and some of the patterns correspond to clinically observed cancer patterns. Furthermore, we examine in detail the time evolution of the average composition of cancer cells and the characteristic length of the microstructures. Our results demonstrate that different sequence of cancer patterns can be obtained by changing the proliferation rate and/or hydrodynamic interactions.

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.

Transport and adsorption under liquid flow: the role of pore geometry

We study here the interplay between transport and adsorption in porous systems with complex geometries under fluid flow. Using a lattice Boltzmann scheme extended to take into account the adsorption at solid/fluid interfaces, we investigate the influence of pore geometry and internal surface roughness on the efficiency of fluid flow and the adsorption of molecular species inside the pore space. We show how the occurrence of roughness on pore walls acts effectively as a modification of the solid/fluid boundary conditions, introducing slippage at the interface. We then compare three common pore geometries, namely honeycomb pores, inverse opal, and materials produced by spinodal decomposition. Finally, we quantify the influence of those three geometries on fluid transport and tracer adsorption. This opens perspectives for the optimization of materials' geometries for applications in dynamic adsorption under fluid flow.

Phase transition of a symmetric diblock copolymer induced by nanorods with different surface chemistry

We investigate the phase transition of a symmetric diblock copolymer induced by nanorods with different surface chemistry. The results demonstrate that the system occurs the phase transition from a disordered structure to ordered parallel lamellae and then to the tilted layered structure as the number of rods increases. The dynamic evolution of the domain size and the order parameter of the microstructure are also examined. Furthermore, the influence of rod property, rod-phase interaction, rod-rod interaction, rod length, and polymerization degree on the behavior of the polymer system is also investigated systematically. Moreover, longer amphiphilic nanorods tend to make the polymer system form the hexagonal structure. It transforms into a perpendicular lamellar structure as the polymerization degree increases. Our simulations provide an efficient method for determining how to obtain the ordered structure on the nanometer scales and design the functional materials with optical, electronic, and magnetic properties.

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.

Rigorous embedding of cell dynamics simulations in the Cahn-Hilliard-Cook framework: Imposing stability and isotropy

We have rigorously analyzed the stability of the efficient cell dynamics simulations (CDS) method by making use of the special properties of the local averaging operator 〈〈*〉〉-* in matrix form. Besides resolving a theoretical issue that has puzzled many over the past three decades, this analysis has considerable practical value: It relates CDS directly to finite-difference approximations of the Cahn-Hilliard-Cook equations and provides a straightforward recipe for replacing the original two- or three-dimensional (2D or 3D) averaging operators in CDS by an equivalent (in terms of stability) discrete Laplacian with superior isotropy and scaling behavior. As such, we open up a route to suppress the unphysical reflection of the computational grid in CDS results (grid artifacts). We found that proper rescaling of discrete Laplacians, needed to employ them in CDS, is equivalent to introducing a well-chosen time step in CDS. In turn, our analysis provides stability conditions for phase-field simulations based on the Cahn-Hilliard-Cook equations. Subsequently, we have quantitatively compared the isotropy and scaling behavior of several discrete 2D or 3D Laplacians, thereby extending the significance of this work to general field-based methodology. We found that all considered discrete Laplacians have equivalent scaling behavior along the Cartesian directions. In addition, and somewhat surprisingly, known "isotropic" discrete Laplacians, i.e., isotropic up to fourth order in |k|, become quite anisotropic for larger wave vectors, whereas "less isotropic" discrete Laplacians (second order) are only slightly anisotropic on the whole |k| range. We identified a hard limit to the accuracy with which the discrete Laplacian can emulate the two important properties of the optimal (continuum) Laplacian, as an improvement of the isotropy, by introducing additional points to the stencil, will negatively affect the scaling behavior. Within this limitation, the discrete compact Laplacians in the DnQm class known from lattice hydrodynamics, D2Q9 in 2D and D3Q19 in 3D, are found to be optimal in terms of isotropy. However, by being only slightly anisotropic on the whole range and enabling larger time steps, the discrete Laplacians that relate to the local averaging operator of Oono and Puri (2D) and Shinozaki and Oono (3D) as well as the less familiar 3D discrete BvV Laplacian developed for dynamic density functional theory are valid alternatives.

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.

Computer simulations of three-dimensional Turing patterns in the Lengyel-Epstein model

We investigate numerically Turing patterns in the Lengyel-Epstein model in three dimensions. In a bulk homogeneous system under periodic boundary conditions, we obtain not only lamellar, cylindrical, and spherical structures but also several interconnected periodic structures including the Schwartz P-surface structure. In order to examine Turing patterns in the conditions accessible experimentally, we consider inhomogeneous systems where a parameter in the reaction-diffusion equations depends on the space coordinate with either Dirichlet or Neumann boundary conditions. In this situation, we find that a perforated-lamellar structure and an Fddd structure, both of which have a uniaxial symmetry, appear depending on the boundary conditions.

Spinodal decomposition of polymer solutions: molecular dynamics simulations of the two-dimensional case

As a generic model system for phase separation in polymer solutions, a coarse-grained model for hexadecane/carbon dioxide mixtures has been studied in two-dimensional geometry. Both the phase diagram in equilibrium (obtained from a finite size scaling analysis of Monte Carlo data) and the kinetics of state changes caused by pressure jumps (studied by large scale molecular dynamics simulations) are presented. The results are compared to previous work where the same model was studied in three-dimensional geometry and under confinement in slit geometry. For deep quenches the characteristic length scale ℓ(t) of the formed domains grows with time t according to a power law close to [Formula: see text]. Since in this problem both the polymer density ρ(p) and the solvent density ρ(s) matter, the time evolution of the density distribution P(L)(ρ(p),ρ(s),t) in L × L subboxes of the system is also analyzed. It is found that in the first stage of phase separation the system separates locally into low density carbon dioxide regions that contain no polymers and regions of high density polymer melt that are supersaturated with this solvent. The further coarsening proceeds via the growth of domains of rather irregular shapes. A brief comparison of our findings with results of other models is given.

Crossover in growth laws for phase-separating binary fluids: molecular dynamics simulations

Pattern and dynamics during phase separation in a symmetrical binary (A+B) Lennard-Jones fluid are studied via molecular dynamics simulations after quenching homogeneously mixed critical (50:50) systems to temperatures below the critical one. The morphology of the domains, rich in A or B particles, is observed to be bicontinuous. The early-time growth of the average domain size is found to be consistent with the Lifshitz-Slyozov law for diffusive domain coarsening. After a characteristic time, dependent on the temperature, we find a clear crossover to an extended viscous hydrodynamic regime where the domains grow linearly with time. Pattern formation in the present system is compared with that in solid binary mixtures, as a function of temperature. Important results for the finite-size and temperature effects on the small-wave-vector behavior of the scattering function are also presented.

Anisotropic velocity statistics of topological defects under shear flow

We report numerical results on the velocity statistics of topological defects during the dynamics of phase ordering and nonrelaxational evolution assisted by an external shear flow. We propose a numerically efficient tracking method for finding the position and velocity of defects and apply it to vortices in a uniform field and dislocations in anisotropic stripe patterns. During relaxational dynamics, the distribution function of the velocity fluctuations is characterized by a dynamical scaling with a scaling function that has a robust algebraic tail with an inverse cube power law. This is characteristic of defects of codimension 2, e.g., point defects in two dimensions and filaments in three dimensions, regardless of whether the motion is isotropic (as for vortices) or highly anisotropic (as for dislocations). However, the anisotropic dislocation motion leads to anisotropic statistical properties when the interaction between defects and their motion is influenced by the presence of an external shear flow transverse to the stripe orientation.

Lattice-Boltzmann-Langevin simulations of binary mixtures

We report a hybrid numerical method for the solution of the Model H fluctuating hydrodynamic equations for binary mixtures. The momentum conservation equations with Landau-Lifshitz stresses are solved using the fluctuating lattice Boltzmann equation while the order parameter conservation equation with Langevin fluxes is solved using stochastic method of lines. Two methods, based on finite difference and finite volume, are proposed for spatial discretization of the order parameter equation. Special care is taken to ensure that the fluctuation-dissipation theorem is maintained at the lattice level in both cases. The methods are benchmarked by comparing static and dynamic correlations and excellent agreement is found between analytical and numerical results. The Galilean invariance of the model is tested and found to be satisfactory. Thermally induced capillary fluctuations of the interface are captured accurately, indicating that the model can be used to study nonlinear fluctuations.

Diffusive domain coarsening: early time dynamics and finite-size effects

We study the diffusive dynamics of phase separation in a symmetric binary (A + B) mixture with a 50:50 composition of A and B particles, following a quench below the demixing critical temperature, both in spatial dimensions d=2 and d=3. The particular focus of this work is to obtain information about the effects of system size and correction to the growth law via the appropriate application of the finite-size scaling method to the results obtained from the Kawasaki exchange Monte Carlo simulation of the Ising model. Observations of only weak size effects and a very small correction to scaling in the growth law are significant. The methods used in this work and information thus gathered will be useful in the study of the kinetics of phase separation in fluids and other problems of growing length scale. We also provide a detailed discussion of the standard methods of understanding simulation results which may lead to inappropriate conclusions.

Domain patterns in a diblock copolymer-diblock copolymer mixture with oscillatory particles

We investigate the orientational order transition of striped patterns in microphase structures of diblock copolymer-diblock copolymer mixtures in the presence of periodic oscillatory particles. Under certain conditions, although the macrophase separation of a system is almost isotropic, microphase separation of one diblock copolymer takes place and becomes anisotropic gradually. By changing the oscillatory frequency and amplitude, the orientational order transition of a striped microphase structure from the state parallel to the oscillatory direction to the state perpendicular to the oscillatory direction is observed. We also find that the order transition occurs when we change the initial composition ratio. Furthermore, we examine the domain size and the orientational order parameter of microstructure in the process of orientational order transition. The results may provide guidance for experimentalists. This model system can also give a simple way to realize orientational order transition of soft materials by changing the oscillatory field.

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.

Polydomain growth at isotropic-nematic transitions in liquid crystalline polymers

We studied the dynamics of isotropic-nematic transitions in liquid crystalline polymers by integrating time-dependent Ginzburg-Landau equations. In a concentrated solution of rodlike polymers, the rotational diffusion constant D(r) of the polymer is severely suppressed by the geometrical constraints of the surrounding polymers so that the rodlike molecules diffuse only along their rod directions. In the early stage of phase transition, the rodlike polymers with nearly parallel orientations assemble to form a nematic polydomain. This polydomain pattern, with characteristic length ℓ, grows with self-similarity in three dimensions over time with an ℓ~t(1/4) scaling law. In the late stage, the rotational diffusion becomes significant, leading to a crossover of the growth exponent from 1/4 to 1/2. This crossover time is estimated to be on the order of t~D(r)(-1). We also examined the time evolution of a pair of disclinations placed in a confined system by solving the same time-dependent Ginzburg-Landau equations in two dimensions. If the initial distance between the disclinations is shorter than some critical length, they approach and annihilate each other; however, at larger initial separations, they are stabilized.

Hydrodynamic effects on spinodal decomposition kinetics in planar lipid bilayer membranes

The formation and dynamics of spatially extended compositional domains in multicomponent lipid membranes lie at the heart of many important biological and biophysical phenomena. While the thermodynamic basis for domain formation has been explored extensively in the past, domain growth in the presence of hydrodynamic interactions both within the (effectively) two-dimensional membrane and in the three-dimensional solvent in which the membrane is immersed has received little attention. In this work, we explore the role of hydrodynamic effects on spinodal decomposition kinetics via continuum simulations of a convective Cahn-Hilliard equation for membrane composition coupled to the Stokes equation. Our approach explicitly includes hydrodynamics both within the planar membrane and in the three-dimensional solvent in the viscously dominated flow regime. Numerical simulations reveal that dynamical scaling breaks down for critical lipid mixtures due to distinct coarsening mechanisms for elongated versus more isotropic compositional lipid domains. The breakdown in scaling should be readily observable in experiments on model membrane systems.

Effects of nanoparticles on the compatibility of PEO-PMMA block copolymers

The compatibility of six kinds of designed poly(ethylene oxide)-block-poly(methyl methacrylate) (PEO-b-PMMA) copolymers was studied at 270, 298 and 400 K via mesoscopic modeling. The values of the order parameters depended on both the structures of the block copolymers and the simulation temperature, while the values of the order parameters of the long chains were higher than those of the short ones; temperature had a more obvious effect on long chains than on the short ones. Plain copolymers doped with poly(ethylene oxide) (PEO) or poly(methyl methacrylate) (PMMA) homopolymer showed different order parameter values. When a triblock copolymer had the same component at both ends and was doped with one of its component polymers as a homopolymer (such as A5B6A5 doped with B6 or A5 homopolymer), the value of its order parameter depended on the simulation temperature. The highest order parameter values were observed for A5B6A5 doped with B6 at 400 K and for A5B6A5 doped with A5 at 270 K. A study of copolymers doped with nanoparticles showed that the mesoscopic phase was influenced by not only the properties of the nanoparticles, such as the size and density, but also the compositions of the copolymers. Increasing the size of the nanoparticles used as a dopant had the most significant effect on the phase morphologies of the copolymers.

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.

Domain coarsening in two dimensions: conserved dynamics and finite-size scaling

We present results from a study of finite-size effect in the kinetics of domain growth with conserved order parameter for a critical quench. Our observation of a weak size effect is a significant and surprising result. For diffusive dynamics, appropriate scaling analysis of Monte Carlo results obtained for small systems using a two-dimensional Ising model also shows that the correction to the expected Lifshitz-Slyozov law for the domain growth is very small. The methods used in this work to understand the growth dynamics should find application in other nonequilibrium systems with increasing length scales.

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.

Turing patterns in three dimensions

We investigate three-dimensional Turing patterns in two-component reaction diffusion systems. The FitzHugh-Nagumo equation, the Brusselator, and the Gray-Scott model are solved numerically in three dimensions. Several interconnected structures of domains as well as lamellar, hexagonal, and spherical domains are obtained as stable motionless equilibrium patterns. The relative stability of these structures is studied analytically based on the reduction approximation. The relation with the microphase-separated structures in block copolymers is also discussed.

Density functional simulation of spontaneous formation of vesicle in block copolymer solutions

The author carries out numerical simulations of vesicle formation based on the density functional theory for block copolymer solutions. It is shown by solving the time evolution equations for concentrations that a polymer vesicle is spontaneously formed from the homogeneous state. The vesicle formation mechanism obtained by this simulation agrees with the results of other simulations based on the particle models as well as experiments. By changing parameters such as the volume fraction of polymers or the Flory-Huggins interaction parameter between the hydrophobic subchains and solvents, the spherical micelles, cylindrical micelles, or bilayer structures can also be obtained. The author also shows that the morphological transition dynamics of the micellar structures can be reproduced by controlling the Flory-Huggins interaction parameter.

Dynamics of ternary mixtures with photosensitive chemical reactions: creating three-dimensionally ordered blends

Using computer simulations, we establish an approach for creating defect-free, periodically ordered polymeric materials. The system involves ABC ternary mixtures where the A and B components undergo a reversible photochemical reaction. In addition, all three components are mutually immiscible and undergo phase separation. Through the simulations, we model the effects of illuminating a three-dimensional (3D) sample with spatially and temporally dependent light irradiation. Experimentally, this situation can be achieved by utilizing both a uniform background light and a spatially localized, higher intensity light, and then rastering a higher-intensity light over the 3D sample. We first focus on the case where the higher-intensity light is held stationary and focused in a distinct region within the system. The C component is seen to displace the A and B within this region and replicate the pattern formed by the higher-intensity light. In effect, one can write a pattern of C onto the AB binary system by focusing the higher-intensity light in the desired arrangement. We isolate the conditions that are necessary for producing clearly written patterns of C (i.e., for obtaining sharp interfaces between the C and A/B domains). We next consider the effect of rastering a higher-intensity light over this sample and find that this light "combs out" defects in the AB blend as it moves through the system. The resulting material displays a defect-free structure that encompasses both a periodic ordering of the A and B domains and a well-defined motif of C. In this manner, one can create hierarchically patterned materials that exhibit periodicity over two distinct length scales. The approach is fully reversible, noninvasive, and points to a novel means of patterning with homopolymers, which normally do not self-assemble into periodic structures.

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.

Interconnected Turing patterns in three dimensions

We study numerically the Turing pattern in three dimensions in a FitzHugh-Nagumo-type reaction-diffusion system. We have found that interconnected periodic domain structures such as a gyroid, Fddd, and perforated lamellar structures appear in three dimensions, which never exist in lower dimensions. The stability analysis of these structures is also performed by means of a mode expansion.