A computer simulation of the time‐dependent Ginzburg–Landau model for spinodal decomposition

Study of Microstructural Morphology of Ti-6Al-4V Alloy by Crystallographic Analysis and Phase Field Simulation

Formation of a habit plane during martensitic transformation is related to an invariant plane strain transformation, which involves dislocation glide and twins. In the current work, the Phenomenological Theory of Martensitic Transformation (PTMT) is employed to study the crystallographic features while the phase field simulation is used to study the microstructure evolution for martensitic transformation of Ti-6Al-4V alloy. Results show that mechanical constraints play a key role in the microstructure evolution. It is shown that a twinned structure with very small twinned variants is geometrically difficult to form due to the lattice parameters of Ti-6Al-4V alloy. It is concluded that the predicted habit plane from the PTMT is consistent with results of the micro-elastic theory. The formation of a triangular morphology is favored geometrically and elastically.

Present and Future Trends in Polymer Blends Technology

The plastics industry shows the highest growth rate of all industrial materials. Within this group the engineering resins (ER) are developing fastest. The review focuses on polymer alloys and blends (PAB) of ER. After the discussion of the economic aspect of PAB the patent literature is briefly reviewed and illustrated by a summary of selected, recent patents on polyphenylene ether blends. There is little doubt that the future trend in PAB technology will require control of morphology, especially in most common multiphase systems. For this reason the last part of the article summarizes the available information of criteria for rheological and thermodynamic control of morphology. The review ends with a brief description of the methods and equipment used by industry to compound and form the PAB.

Gelation Impairs Phase Separation and Small Molecule Migration in Polymer Mixtures

Surface segregation of the low molecular weight component of a polymeric mixture is a ubiquitous phenomenon that leads to degradation of industrial formulations. We report a simultaneous phase separation and surface migration phenomena in oligomer-polymer ( O P ) and oligomer-gel ( O G ) systems following a temperature quench that induces demixing of components. We compute equilibrium and time varying migrant (oligomer) density profiles and wetting layer thickness in these systems using coarse grained molecular dynamics (CGMD) and mesoscale hydrodynamics (MH) simulations. Such multiscale methods quantitatively describe the phenomena over a wide range of length and time scales. We show that surface migration in gel-oligomer systems is significantly reduced on account of network elasticity. Furthermore, the phase separation processes are significantly slowed in gels leading to the modification of the well known Lifshitz-Slyozov-Wagner (LSW) law ℓ ( τ ) ∼ τ 1 / 3 . Our work allows for rational design of polymer/gel-oligomer mixtures with predictable surface segregation characteristics that can be compared against experiments.

Morphology formation in binary mixtures upon gradual destabilisation

Spontaneous liquid-liquid phase separation is commonly understood in terms of phenomenological mean-field theories. These theories correctly predict the structural features of the fluid at sufficiently long time scales and wavelengths. However, these conditions are not met in various examples in biology and materials science where the mixture is slowly destabilised, and phase separation is strongly affected by critical thermal fluctuations. We propose a mechanism of pretransitional structuring of a mixture that approaches the miscibility gap and predict scaling relations that describe how the characteristic feature size of the emerging morphology decreases with an increasing quench rate. These predictions quantitatively agree with our kinetic Monte Carlo and molecular dynamics simulations of a phase-separating binary mixture, as well as with previously reported experimental observations. We discuss how these predictions are affected by non-conserved order parameters (e.g., due to chemical reactions or alignment of liquid-crystalline molecules), hydrodynamics and active transport.

Morphological-evolution pathway during phase separation in polymer solutions with highly asymmetrical miscibility gap

Microstructural evolution during thermally induced phase separation in a polymer solution with a highly asymmetrical miscibility gap is analyzed using Flory-Huggins thermodynamics and nonlinear Cahn-Hilliard kinetics. Computer simulation results demonstrate that, in contrast to systems with symmetric miscibility gaps, interesting morphological-evolution pathways (MEPs) including cluster-to-percolation and percolation-to-cluster transitions are identified. These MEPs are rationalized according to asynchronous evolution of the two product phases. For a highly asymmetric polymer system, the initial solution composition is also found to play a critical role in determining the MEPs. According to the simulation results a map of MEPs in terms of initial solution composition and aging time of phase separation is established to guide future microstructural design in asymmetrical polymer solutions.

Numerical Simulations of Directed Self-Assembly in Diblock Copolymer Films using Zone Annealing and Pattern Templating

Bulk fabrication of surface patterns with sub-20 nm feature sizes is immensely desirable for many existing and emerging technologies. Directed self-assembly (DSA) of block copolymers (BCPs) has been a recently demonstrated approach to achieve such feature resolution over large-scale areas with minimal defect populations. However, much work remains to understand and optimize DSA methods in order to move this field forward. This paper presents large-scale numerical simulations of zone annealing and chemo-epitaxy processing of BCP films to achieve long-range orientational order. The simulations utilize a Time-Dependent Ginzburg-Landau model and parallel processing to elucidate relationships between the magnitude and velocity of a moving thermal gradient and the resulting BCP domain orientations and defect densities. Additional simulations have been conducted to study to what degree orientational order can be further improved by combining zone annealing and chemo-epitaxy techniques. It is found that these two DSA methods do synergistically enhance long-range order with a particular relationship between thermal gradient velocity and chemical template spacing.

Instability in reaction-superdiffusion systems

We study the effect of superdiffusion on the instability in reaction-diffusion systems. It is shown that reaction-superdiffusion systems close to a Turing instability are equivalent to a time-dependent Ginzburg-Landau model and the corresponding free energy is introduced. This generalized free energy which depends on the superdiffusion exponent governs the stability, dynamics, and the fluctuations of reaction-superdiffusion systems near the Turing bifurcation. In addition, we show that for a general n-component reaction-superdiffusion system, a fractional complex Ginzburg-Landau equation emerges as the amplitude equation near a Hopf instability. Numerical simulations of this equation are carried out to illustrate the effect of superdiffusion on spatiotemporal patterns. Finally, the effect of superdiffusion on the instability in the Brusselator model, as a special case of reaction-diffusion systems, is studied. In general, superdiffusion introduces a new parameter that changes the behavior of the system near the instability.

Lifetime of metastable states in a Ginzburg-Landau system: Numerical simulations at large driving forces

We developed a "brute-force" simulation method and conducted numerical "experiments" on homogeneous nucleation in an isotropic system at large driving forces (not small supersaturations) using the stochastic Ginzburg-Landau approach. Interactions in the system are described by the asymmetric (no external field), athermal (temperature-independent driving force), tangential (simple phase diagram) Hamiltonian, which has two independent "drivers" of the phase transition: supersaturation and thermal noise. We obtained the probability distribution function of the lifetime of the metastable state and analyzed its mean value as a function of the supersaturation, noise strength, and volume. We also proved the nucleation theorem in the mean-field approximation. The results allowed us to find the thermodynamic properties of the barrier state and conclude that at large driving forces the fluctuating volumes are not independent.

Structuring of polymer solutions upon solvent evaporation

The morphology of solution-cast, phase-separated polymers becomes finer with increasing solvent evaporation rate. We address this observation theoretically for a model polymer where demixing is induced by steady solvent evaporation. In contrast to what is the case for a classical, thermal quench involving immiscible blends, the spinodal instability initially develops slowly and the associated length scale is not time invariant but decreases with time as t(-1/2). After a time lag, phase separation accelerates. Time lag and characteristic length exhibit power-law behavior as a function of the evaporation rate with exponents of -2/3 and -1/6. Interestingly, at later stages the spinodal structure disappears completely while a second length scale develops. The associated structure coarsens but does not follow the usual Lifshitz-Slyozov-Wagner kinetics.

Polymerization-induced spinodal decomposition of ethylene glycol∕phenolic resin solutions under electric fields

Temporal evolution of polymerization-induced spinodal decomposition (PISD) under electric fields was investigated numerically in ethylene glycol∕phenolic resin solutions with different initial composition. A model composed of the nonlinear Cahn-Hilliard-Cook equation for spinodal decomposition and a rate equation for curing reaction was utilized to describe the PISD phenomenon. As initial composition varied, deformed droplet-like and aligned bi-continuous structures were observed in the presence of an electric field. Moreover, the anisotropic parameter (D), determined from the 2D-FFT power spectrum, was employed to quantitatively characterize the degree of morphology anisotropy. The value of D increased quickly in the early stage and then decreased in the intermediate stage of spinodal decomposition, which was attributed to the resistance of coarsening process to morphology deformation and the decline of electric stress caused by polymerization reaction. The results can also provide a guidance on how to control the morphology of monolithic porous polymer and carbon materials with anisotropic structures.

Phase-field approach to implicit solvation of biomolecules with Coulomb-field approximation

A phase-field variational implicit-solvent approach is developed for the solvation of charged molecules. The starting point of such an approach is the representation of a solute-solvent interface by a phase field that takes one value in the solute region and another in the solvent region, with a smooth transition from one to the other on a small transition layer. The minimization of an effective free-energy functional of all possible phase fields determines the equilibrium conformations and free energies of an underlying molecular system. All the surface energy, the solute-solvent van der Waals interaction, and the electrostatic interaction are coupled together self-consistently through a phase field. The surface energy results from the minimization of a double-well potential and the gradient of a field. The electrostatic interaction is described by the Coulomb-field approximation. Accurate and efficient methods are designed and implemented to numerically relax an underlying charged molecular system. Applications to single ions, a two-plate system, and a two-domain protein reveal that the new theory and methods can capture capillary evaporation in hydrophobic confinement and corresponding multiple equilibrium states as found in molecular dynamics simulations. Comparisons of the phase-field and the original sharp-interface variational approaches are discussed.

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.