Study of phase-separation dynamics by use of cell dynamical systems. II. Two-dimensional demonstrations

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.

Hybrid Time-Dependent Ginzburg–Landau Simulations of Block Copolymer Nanocomposites: Nanoparticle Anisotropy

Block copolymer melts are perfect candidates to template the position of colloidal nanoparticles in the nanoscale, on top of their well-known suitability for lithography applications. This is due to their ability to self-assemble into periodic ordered structures, in which nanoparticles can segregate depending on the polymer–particle interactions, size and shape. The resulting coassembled structure can be highly ordered as a combination of both the polymeric and colloidal properties. The time-dependent Ginzburg–Landau model for the block copolymer was combined with Brownian dynamics for nanoparticles, resulting in an efficient mesoscopic model to study the complex behaviour of block copolymer nanocomposites. This review covers recent developments of the time-dependent Ginzburg–Landau/Brownian dynamics scheme. This includes efforts to parallelise the numerical scheme and applications of the model. The validity of the model is studied by comparing simulation and experimental results for isotropic nanoparticles. Extensions to simulate nonspherical and inhomogeneous nanoparticles are discussed and simulation results are discussed. The time-dependent Ginzburg–Landau/Brownian dynamics scheme is shown to be a flexible method which can account for the relatively large system sizes required to study block copolymer nanocomposite systems, while being easily extensible to simulate nonspherical nanoparticles.

Two-dimensional Cahn-Hilliard simulations for coarsening kinetics of spinodal decomposition in binary mixtures

The evolution of the microstructure due to spinodal decomposition in phase separated mixtures has a strong impact on the final material properties. In the late stage of coarsening, the system is characterized by the growth of a single characteristic length scale L ∼ Ct^α. To understand the structure-property relationship, the knowledge of the coarsening exponent α and the coarsening rate constant C is mandatory. Since the existing literature is not entirely consistent, we perform phase field simulations based on the Cahn-Hilliard equation. We restrict ourselves to binary mixtures using a symmetric Flory-Huggins free energy and a constant composition-independent mobility term and show that the coarsening for off-critical mixtures is slower than the expected t^1/3-growth. Instead, we find α to be dependent on the mixture composition and associate this with the observed morphologies. Finally, we propose a model to describe the complete coarsening kinetics including the rate constant C.

Ordering kinetics in the active model B

We undertake a detailed numerical study of the Active Model B proposed by Wittkowski et al., [Nature Commun. 5, 4351 (2014)]2041-172310.1038/ncomms5351. We find that the introduction of activity has a drastic effect on the ordering kinetics. First, the domain growth law shows a crossover from the usual Lifshitz-Slyozov growth law for phase separation (L∼t^{1/3}, where t is the time) to a novel growth law (L∼t^{1/4}) at late times. Second, the equal-time correlation function of the density field exhibits dynamical scaling for a given activity strength λ, but the scaling function depends on λ.

Cell dynamics simulations of coupled charge and magnetic phase transformation in correlated oxides

We present a comprehensive numerical study on the kinetics of phase transition that is characterized by two nonconserved scalar order parameters coupled by a special linear-quadratic interaction. This particular Ginzburg-Landau theory has been proposed to describe the coupled charge and magnetic transition in nickelates and the collinear stripe phase in cuprates. The inhomogeneous state of such systems at low temperatures consists of magnetic domains separated by quasimetallic domain walls where the charge order is reduced. By performing large-scale cell dynamics simulations, we find a two-stage phase-ordering process in which a short period of independent evolution of the two order parameters is followed by a correlated coarsening process. The long-time growth and coarsening of magnetic domains is shown to follow the Allen-Cahn power law. We further show that the nucleation-and-growth dynamics during phase transformation to the ordered states is well described by the Kolmogorov-Johnson-Mehl-Avrami theory in two dimensions. On the other hand, the presence of quasimetallic magnetic domain walls in the ordered states gives rise to a very different kinetics for transformation to the high-temperature paramagnetic phase. In this scenario, the phase transformation is initiated by the decay of magnetic domain walls into two insulator-metal boundaries, which subsequently move away from each other. Implications of our findings to recent nano-imaging experiments on nickelates are also discussed.

Domain growth in ferronematics: slaved coarsening, emergent morphologies and growth laws

Ferronematics (FNs) are suspensions of magnetic nanoparticles in nematic liquid crystals (NLCs). They have attracted much experimental attention, and are of great interest both scientifically and technologically. There are very few theoretical studies of FNs, even in equilibrium. In this paper, we study the non-equilibrium phenomenon of domain growth after a thermal quench (or coarsening) in this coupled system. Our modeling is based on coupled time-dependent Ginzburg-Landau (TDGL) equations for two order parameters: the LC tensor order parameter Q, and the magnetization M. We consider both shallow and deep quenches from a high-temperature disordered phase. The system coarsens by the collision and annihilation of topological defects. We focus on slaved coarsening, where a disordered Q (or M) field is driven to coarsen by an ordered M (or Q) field. We present detailed results for the morphologies and growth laws, which exhibit unusual features purely due to the magneto-nematic coupling. To the best of our knowledge, this is the first study of non-equilibrium phenomena in FNs.

Reduced Model for Properties of Multiscale Porous Media with Changing Geometry

In this paper, we consider an important problem for modeling complex coupled phenomena in porous media at multiple scales. In particular, we consider flow and transport in the void space between the pores when the pore space is altered by new solid obstructions formed by microbial growth or reactive transport, and we are mostly interested in pore-coating and pore-filling type obstructions, observed in applications to biofilm in porous media and hydrate crystal formation, respectively. We consider the impact of these obstructions on the macroscopic properties of the porous medium, such as porosity, permeability and tortuosity, for which we build an experimental probability distribution with reduced models, which involves three steps: (1) generation of independent realizations of obstructions, followed by, (2) flow and transport simulations at pore-scale, and (3) upscaling. For the first step, we consider three approaches: (1A) direct numerical simulations (DNS) of the PDE model of the actual physical process called BN which forms the obstructions, and two non-DNS methods, which we call (1B) CLPS and (1C) LP. LP is a lattice Ising-type model, and CLPS is a constrained version of an Allen–Cahn model for phase separation with a localization term. Both LP and CLPS are model approximations of BN, and they seek local minima of some nonconvex energy functional, which provide plausible realizations of the obstructed geometry and are tuned heuristically to deliver either pore-coating or pore-filling obstructions. Our methods work with rock-void geometries obtained by imaging, but bypass the need for imaging in real-time, are fairly inexpensive, and can be tailored to other applications. The reduced models LP and CLPS are less computationally expensive than DNS, and can be tuned to the desired fidelity of the probability distributions of upscaled quantities.

Surface-directed spinodal decomposition on chemically patterned substrates

Surface-directed spinodal decomposition (SDSD) is the kinetic interplay of phase separation and wetting at a surface. This process is of great scientific and technological importance. In this paper, we report results from a numerical study of SDSD on a chemically patterned substrate. We consider simple surface patterns for our simulations, but most of the results apply for arbitrary patterns. In layers near the surface, we observe a dynamical crossover from a surface-registry regime to a phase-separation regime. We study this crossover using layerwise correlation functions and structure factors and domain length scales.

Dynamical self-consistent field theory captures multi-scale physics during spinodal decomposition in a symmetric binary homopolymer blend

We study the spinodal decomposition in a symmetric, binary homopolymer blend using our recently developed dynamical self-consistent field theory. By taking the extremal solution of a dynamical functional integral, the theory reduces the interacting, multi-chain dynamics to a Smoluchowski equation describing the statistical dynamics of a single, unentangled chain in a self-consistent, time-dependent, mean force-field. We numerically solve this equation by evaluating averages over a large ensemble of replica chains, each one of which obeys single-chain Langevin dynamics, subject to the mean field. Following a quench from the disordered state, an early time spinodal instability in the blend composition develops, before even one Rouse time elapses. The dominant, unstable, growing wavelength is on the order of the coil size. The blend then enters a late-time, t, scaling regime with a growing domain size that follows the expected Lifshitz-Slyozov-Wagner t^1/3 power law, a characteristic of a diffusion-driven coarsening process. These results provide a satisfying test of this new method, which correctly captures both the early and late time physics in the blend. Our simulation spans five orders-of-magnitude in time as the domains coarsen to 20 times the coil size, while remaining faithful to the dynamics of the microscopic chain model.

Hybrid patterns from directed self-assembly of diblock copolymers by chemical patterns

The surface affinity is a critical factor for controlling the formation of monolayer nanostructures in block copolymer thin films. In general, strong surface affinity tends to induce the formation of domains with low spontaneous curvature. Abiding by this principle, we propose a facile chemoepitaxial scheme for producing large-scale ordered hybrid patterns by the directed self-assembly of diblock copolymers. The guiding chemical pattern is designed as periodic stripes with alternately changing surface affinities. As a consequence, two different geometries of domains are formed on the stripes with different affinities. The self-assembly process of block copolymers guided by the stripe patterns is investigated using cell dynamics simulations based on time-dependent Ginzburg-Landau theory, and the kinetic stability diagram is estimated. Hybrid patterns are successfully achieved with both cylinder-forming and sphere-forming diblock copolymers. In the cylinder-forming system, the major hybrid pattern exhibiting a considerable stability window is composed of parallel cylinders and perforated lamellae, while it is composed of monolayer spheres and parallel cylinders in the other system. Encouragingly, the chemoepitaxial method is valid till the period of the guiding pattern is a large multiple of the domain spacing. The chemoepitaxial scheme demonstrated in this work serves as a nice supplement to the graphoepitaxial one proposed in our previous work.

Hybrid line-dot nanopatterns from directed self-assembly of diblock copolymers by trenches

We demonstrate that the directed self-assembly of AB diblock copolymers by periodic trenches can be used to fabricate large-scale ordered hybrid line-dot nanopatterns in addition to a defect-free dot nanopattern. The formation of line or dot nanopatterns in thin films with proper surface affinities is controlled by the film thickness, which is modulated by a topographic pattern consisting of steps and trenches. Two kinds of line-dot nanopatterns are achieved with cylinder-forming and sphere-forming copolymers, respectively. One kind of hybrid nanopatterns is composed of perpendicularly standing cylinders (dots) on the steps and parallel monolayer cylinders (lines) within the trenches, while the dots of the other kind are replaced by monolayer spheres on the steps. The thermodynamic stability region of target hybrid nanopatterns is identified by constructing two-dimensional phase diagrams with respect to two control parameters of step height and film thickness using self-consistent field theory. Furthermore, a process window of the line-dot nanopatterns is estimated using cell dynamics simulations based on time-dependent Ginzburg-Landau theory, confirming their feasibility in kinetics.

Chiral selection of single helix formed by diblock copolymers confined in nanopores

Chiral selection has attracted tremendous attention from the scientific communities, especially from biologists, due to the mysterious origin of homochirality in life. The self-assembly of achiral block copolymers confined in nanopores offers a simple but useful model of forming helical structures, where the helical structures possess random chirality selection, i.e. equal probability of left-handedness and right-handedness. Based on this model, we study the stimulus-response of chiral selection to external conditions by introducing a designed chiral pattern onto the inner surface of a nanopore, aiming to obtain a defect-free helix with controllable homochirality. A cell dynamics simulation based on the time-dependent Ginzburg-Landau theory is carried out to demonstrate the tuning effect of the patterned surface on the chiral selection. Our results illustrate that the chirality of the helix can be successfully controlled to be consistent with that of the tailored surface patterns. This work provides a successful example for the stimulus response of the chiral selection of self-assembled morphologies from achiral macromolecules to external conditions, and hence sheds light on the understanding of the mechanism of the stimulus response.

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.

Ordering Dynamics in Neuron Activity Pattern Model: An Insight to Brain Functionality

We study the domain ordering kinetics in d = 2 ferromagnets which corresponds to populated neuron activities with both long-ranged interactions, V(r) ∼ r⁻ⁿ and short-ranged interactions. We present the results from comprehensive Monte Carlo (MC) simulations for the nonconserved Ising model with n ≥ 2, interaction range considering near and far neighbors. Our model results could represent the long-ranged neuron kinetics (n ≤ 4) in consistent with the same dynamical behaviour of short-ranged case (n ≥ 4) at far below and near criticality. We found that emergence of fast and slow kinetics of long and short ranged case could imitate the formation of connections among near and distant neurons. The calculated characteristic length scale in long-ranged interaction is found to be n independent (L(t) ∼ t^1/(n−2)), whereas short-ranged interaction follows L(t) ∼ t^1/2 law and approximately preserve universality in domain kinetics. Further, we did the comparative study of phase ordering near the critical temperature which follows different behaviours of domain ordering near and far critical temperature but follows universal scaling law.

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.

Instability-induced pattern formation of photoactivated functional polymers

Since the pioneering work of Turing on the formation principles of animal coat patterns [Turing AM (1952) Phil Trans R Soc Lond B 237(641):37-72], such as the stripes of a tiger, great effort has been made to understand and explain various phenomena of self-assembly and pattern formation. Prominent examples are the spontaneous demixing in emulsions, such as mixtures of water and oil [Herzig EM, et al. (2007) Nat Mater 6:966-971]; the distribution of matter in the universe [Kibble TWB (1976) J Phys A: Math Gen 9(8):1387]; surface reconstruction in ionic crystals [Clark KW, et al. (2012) Nanotechnol 23(18):185306]; and the pattern formation caused by phase transitions in metal alloys, polymer mixtures and binary Bose-Einstein condensates [Sabbatini J, et al. (2011) Phys Rev Lett 107:230402]. Photoactivated pattern formation in functional polymers has attracted major interest due to its potential applications in molecular electronics and photoresponsive systems. Here we demonstrate that photoactivated pattern formation on azobenzene-containing polymer films can be entirely explained by the physical concept of phase separation. Using experiments and simulations, we show that phase separation is caused by an instability created by the photoactivated transitions between two immiscible states of the polymer. In addition, we have shown in accordance with theory, that polarized light has a striking effect on pattern formation indicated by enhanced phase separation.

Kinetics of lamellar formation on sparsely stripped patterns

Chemical epitaxy based on the self-assembly of block copolymers is viewed as a promising technique to achieve ordered patterns on a large scale. Herein, we study the kinetics of lamellar formation of block copolymers under the direction of sparsely stripped patterns using cell dynamics simulations of the time-dependent Ginzburg-Landau theory. First, a scaling law is unveiled with the ordering time of lamellae, tp, with respect to the multiples between the periods of lamellae and stripe patterns, which is consistent with the power law evolution of the correlation length existing in the bulk phase of lamellae. Second, the tolerative windows of perfect order, with deviation from integer multiples, are also estimated from the aspect of kinetics. The results of the ordering time and tolerative windows are of great interest for relevant experiments or applications. Finally, a two-stage evolution is explored during the pattern formation of chemical epitaxy by probing into the evolution of defects, which is of fundamental interest for us to understand the coarsening kinetics of block copolymers under the direction of chemical patterns.

Phase separation in antisymmetric films: a molecular dynamics study

We have used molecular dynamics (MD) simulations to study phase-separation kinetics in a binary fluid mixture (AB) confined in an antisymmetric thin film. One surface of the film (located at z = 0) attracts the A-atoms, and the other surface (located at z = D) attracts the B-atoms. We study the kinetic processes which lead to the formation of equilibrium morphologies subsequent to a deep quench below the miscibility gap. In the initial stages, one observes the formation of a layered structure, consisting of an A-rich layer followed by a B-rich layer at z = 0; and an analogous structure at z = D. This multi-layered morphology is time-dependent and propagates into the bulk, though it may break up into a laterally inhomogeneous structure at a later stage. We characterize the evolution morphologies via laterally averaged order parameter profiles; the growth laws for wetting-layer kinetics and layer-wise length scales; and the scaling properties of layer-wise correlation functions.

Using light to guide the motion of nanorods in photoresponsive binary blends: designing hierarchically structured nanocomposites

One of the challenges in creating high-performance polymer nanocomposites is establishing effective routes for tailoring the morphology of both the polymer mixture and the dispersed nanoparticles, which contribute desirable optical, electrical, and mechanical properties. Using computational modeling, we devise an effective method for simultaneously controlling the spatial regularity of the polymer phases and the distribution of the rods within this matrix. We focus on mixtures of photosensitive AB binary blends and A-coated nanorods; in the presence of light, the binary blends undergo a reversible chemical reaction and phase separation to yield a morphology resembling that of microphase-separated diblock copolymers. We simulate the effects of illuminating this sample with a uniform background light and a higher intensity, spatially localized light, which is rastered over the sample with a velocity v. The resulting material displays a periodically ordered, essentially defect-free morphology, with the A-like nanoparticles localized in lamellar A domains. The dynamic behavior of the rods within this system can be controlled by varying the velocity v and Γ2, the reaction rate coefficient produced by the higher intensity light. Specifically, the rastering light can drive the rods to be "pushed" along the lamellar domains or oriented perpendicular to these stripes. Given these attributes, we isolate scenarios where the system encompasses a complex hierarchical structure, with rods that are simultaneously ordered along two distinct directions within the periodic matrix. Namely, the rods form long nanowires that span the length of the sample and lie perpendicular to these wires in regularly spaced A lamellae. Hence, our approach points to new routes for producing self-organized rectangular grids, which can impart remarkable optoelectronic or mechanical properties to the materials.

Morphological phase separation in unstable thin films: pattern formation and growth

We present results from a comprehensive numerical study of morphological phase separation (MPS) in unstable thin liquid films on a 2-dimensional substrate. We study the quantitative properties of the evolution morphology via several experimentally relevant markers, e.g., correlation function, structure factor, domain-size and defect-size probability distributions, and growth laws. Our results suggest that the late-stage morphologies exhibit dynamical scaling, and their evolution is self-similar in time. We emphasize the analogies and differences between MPS in films and segregation kinetics in unstable binary mixtures.

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.

Dynamics of dislocations in a two-dimensional block copolymer system with hexagonal symmetry

Block copolymer thin films have attracted considerable attention for their ability to self-assemble into nanometre-scale architectures. Recent advances in the use of block copolymer thin films as nano-lithographic masks have driven research efforts in order to have better control of long-range ordering in the plane of the film. Irrespective of the method of sample preparation, different quasi-two-dimensional systems with hexagonal symmetry unavoidably contain translational defects, called dislocations. Dislocations control the process of coarsening in the nano/meso-scales and provide one of the most important mechanisms of length-scale selection in hexagonal patterns. Although in the last decade the nonlinear dynamics of topological defects in quasi-two-dimensional systems has witnessed significant progress, still little is known about the role of external fields on the creation and annihilation mechanisms involved in the relaxation process towards equilibrium states. In this paper, the dynamics of dislocations in non-optimal hexagonal patterns is studied in the framework of the Ohta-Kawasaki model for a diblock copolymer. Measurements of the climb and glide velocities as a function of the wave vector deformation reveal the main mechanisms of relaxation associated with the motion of dislocations.

Asymmetric oscillations during phase separation under continuous cooling: A simple model

We investigate the phase separation of binary mixtures under continuous cooling using the Cahn-Hilliard equation including the effect of gravity. In our simple model, sedimentation is accounted for by instantaneously "removing" droplets from the supersaturated mixture into the coexisting phase once the droplets have reached a defined maximum size. Our model predicts an oscillatory variation of turbidity. Depending on the composition, either both phases oscillate (symmetric oscillations) or only one of the phases oscillates (asymmetric oscillations). In the asymmetric case, droplet sedimentation from the majority phase into the minority phase reduces supersaturation in the minority phase. This inhibits droplet formation in the minority phase. The cooling rate dependence of the period agrees with experimental results.

Modeling the phase separation in binary lipid membrane under externally imposed oscillatory shear flow

By adding external velocity terms, the two-dimensional time-dependent Ginzburg-Landau (TDGL) equations are modified. Based on this, the phase separation in binary lipid membrane under externally imposed oscillatory shear flow is numerically modeled employing the Cell Dynamical System (CDS) approach. Considering shear flows with different frequencies and amplitudes, several aspects of such a phase evolving process are studied. Firstly, visualized results are shown via snapshot figures of the membrane shape. And then, the simulated scattering patterns at typical moments are presented. Furthermore, in order to more quantitatively discuss this phase-separation process, the time growth laws of the characteristic domain sizes in both directions parallel and perpendicular to the flow are investigated for each case. Finally, the peculiar rheological properties of such binary lipid membrane system have been discussed, mainly the normal stress difference and the viscoelastic complex shear moduli.

Phase-field modeling of the formation of lamellar nanostructures in diblock copolymer thin films under inplanar electric fields

Recent experiments show that external inplanar electric field can be employed to guide the molecular self-assembly in diblock copolymer (BCP) thin films to form lamellar nanostructures with potential applications in nanotechnology. We study this self-assembly process through a detailed coarse-grained phase-separation modeling. During the process, the free energy of the BCP films is modeled as the Ginzburg-Landau free energy with nonlocal interaction and electrostatic coupling. The resulting Cahn-Hilliard (CH) equation is solved using an efficient semi-implicit Fourier-spectral algorithm. Numerical results show that the morphology of order parameter formed in either symmetric or asymmetric BCP thin films is strongly influenced by the electric field. For symmetrical BCPs, highly ordered lamellar nanostructures evolved along the direction of the electric field. Phase nucleation and dislocation climbing in the BCP films predicted by the numerical simulation are in a good agreement with those observed in recent BCP electronanolithography. For asymmetrical BCPs, numerical simulation shows that nanodots are guided to align to the electric field. Furthermore, in the case of high electric field, nanodots formed in asymmetrical BCPs may further convert into highly ordered lamellar nanostructures (sphere-to-cylinder transition) parallel to the electric field. Effects of the magnitude of electric field, BCP asymmetry, and molecular interaction of BCPs on the self-assembly process are examined in detail using the numerical scheme developed in this study. The present study can be used for the prediction of the formation of nanostructures in BCP thin films and the quality control of BCP-based nanomanufacturing through optimizing the external electric fields.