Thermodynamics, formation dynamics, and structural correlations in the bulk amorphous phase of the phase-field crystal model.

We investigate bulk thermodynamic and microscopic structural properties of amorphous solids in the framework of the phase-field crystal (PFC) model. These are metastable states with a non-uniform density distribution, having no long-range order. From extensive numerical simulations, we determine the distribution of free energy density values in varying size amorphous systems and also the point-to-set correlation length, which is the radius of the largest volume of amorphous one can take while still having the particle arrangements within the volume determined by the particle ordering at the surface of the chosen volume. We find that in the thermodynamic limit, the free energy density of the amorphous tends toward a value that has a slight dependence on the initial state from which it was formed-i.e., it has a formation history dependence. The amorphous phase is observed to form on both sides of the liquid linear-stability limit, showing that the liquid to amorphous transition is first order, with an associated finite free energy barrier when the liquid is metastable. In our simulations, this is demonstrated when the noise in the initial density distribution is used to induce nucleation events from the metastable liquid. Depending on the strength of the initial noise, we observe a variety of nucleation pathways, in agreement with previous results for the PFC model, which show that amorphous precursor mediated multi-step crystal nucleation can occur in colloidal systems.

Phase-Field Modeling of Biomineralization in Mollusks and Corals: Microstructure vs Formation Mechanism.

While biological crystallization processes have been studied on the microscale extensively, there is a general lack of models addressing the mesoscale aspects of such phenomena. In this work, we investigate whether the phase-field theory developed in materials' science for describing complex polycrystalline structures on the mesoscale can be meaningfully adapted to model crystallization in biological systems. We demonstrate the abilities of the phase-field technique by modeling a range of microstructures observed in mollusk shells and coral skeletons, including granular, prismatic, sheet/columnar nacre, and sprinkled spherulitic structures. We also compare two possible micromechanisms of calcification: the classical route, via ion-by-ion addition from a fluid state, and a nonclassical route, crystallization of an amorphous precursor deposited at the solidification front. We show that with an appropriate choice of the model parameters, microstructures similar to those found in biomineralized systems can be obtained along both routes, though the time-scale of the nonclassical route appears to be more realistic. The resemblance of the simulated and natural biominerals suggests that, underneath the immense biological complexity observed in living organisms, the underlying design principles for biological structures may be understood with simple math and simulated by phase-field theory.

Dynamical Density Functional Theory for the Drying and Stratification of Binary Colloidal Dispersions.

We develop a dynamical density functional theory based model for the drying of colloidal films on planar surfaces. We consider mixtures of two different sizes of hard-sphere colloids. Depending on the solvent evaporation rate and the initial concentrations of the two species, we observe varying degrees of stratification in the final dried films. Our model predicts the various structures described in the literature previously from experiments and computer simulations, in particular the small-on-top stratified films. Our model also includes the influence of adsorption of particles to the interfaces.

Phase-field modelling of the effect of density change on solidification revisited: model development and analytical solutions for single component materials.

In this paper the development of a physically consistent phase-field theory of solidification shrinkage is presented. The coarse-grained hydrodynamic equations are derived directly from the N-body Hamiltonian equations in the framework of statistical physics, while the constitutive relations are developed in the framework of the standard phase-field theory, by following the variational formalism and the principles of non-equilibrium thermodynamics. To enhance the numerical practicality of the model, quasi-incompressible hydrodynamic equations are derived, where sound waves are absent (but density change is still possible), and therefore the time scale of solidification is accessible in numerical simulations. The model development is followed by a comprehensive mathematical analysis of the equilibrium and propagating one-dimensional solid-liquid interfaces for different density-phase couplings. It is shown, that the fluid flow decelerates/accelerates the solidification front in case of shrinkage/expansion of the solid compared to the case when no density contrast is present between the phases. Furthermore, such a free energy construction is proposed, in which the shape of the equilibrium planar phase-field interface is independent from the density-phase coupling, and the equilibrium interface represents an exact propagating planar interface solution of the quasi-incompressible hydrodynamic equations. Our results are in agreement with previous theoretical predictions.

Hydrodynamic theory of freezing: Nucleation and polycrystalline growth.

Structural aspects of crystal nucleation in undercooled liquids are explored using a nonlinear hydrodynamic theory of crystallization proposed recently [G. I. Tóth et al., J. Phys.: Condens. Matter 26, 055001 (2014)JCOMEL0953-898410.1088/0953-8984/26/5/055001], which is based on combining fluctuating hydrodynamics with the phase-field crystal theory. We show that in this hydrodynamic approach not only homogeneous and heterogeneous nucleation processes are accessible, but also growth front nucleation, which leads to the formation of new (differently oriented) grains at the solid-liquid front in highly undercooled systems. Formation of dislocations at the solid-liquid interface and interference of density waves ahead of the crystallization front are responsible for the appearance of the new orientations at the growth front that lead to spherulite-like nanostructures.

Phase-field modeling of isothermal quasi-incompressible multicomponent liquids.

In this paper general dynamic equations describing the time evolution of isothermal quasi-incompressible multicomponent liquids are derived in the framework of the classical Ginzburg-Landau theory of first order phase transformations. Based on the fundamental equations of continuum mechanics, a general convection-diffusion dynamics is set up first for compressible liquids. The constitutive relations for the diffusion fluxes and the capillary stress are determined in the framework of gradient theories. Next the general definition of incompressibility is given, which is taken into account in the derivation by using the Lagrange multiplier method. To validate the theory, the dynamic equations are solved numerically for the quaternary quasi-incompressible Cahn-Hilliard system. It is demonstrated that variable density (i) has no effect on equilibrium (in case of a suitably constructed free energy functional) and (ii) can influence nonequilibrium pattern formation significantly.

Phase-field theory of multicomponent incompressible Cahn-Hilliard liquids.

In this paper, a generalization of the Cahn-Hilliard theory of binary liquids is presented for multicomponent incompressible liquid mixtures. First, a thermodynamically consistent convection-diffusion-type dynamics is derived on the basis of the Lagrange multiplier formalism. Next, a generalization of the binary Cahn-Hilliard free-energy functional is presented for an arbitrary number of components, offering the utilization of independent pairwise equilibrium interfacial properties. We show that the equilibrium two-component interfaces minimize the functional, and we demonstrate that the energy penalization for multicomponent states increases strictly monotonously as a function of the number of components being present. We validate the model via equilibrium contact angle calculations in ternary and quaternary (four-component) systems. Simulations addressing liquid-flow-assisted spinodal decomposition in these systems are also presented.

Phase field modelling of spinodal decomposition in the oil/water/asphaltene system.

In this paper the quantitative applicability of van der Sman/van der Graaf type Ginzburg-Landau theories of surfactant assisted phase separation [van der Sman et al., Rheol. Acta, 2006, 46, 3] is studied for real systems displaying high surfactant concentrations at the liquid-liquid interface. The model is applied for the water/heptane/asphaltene system (a model of heavy crude oil), for which recent molecular dynamics (MD) simulations provide microscopic data needed to calibrate the theory. A list of general requirements is set up first, which is then followed by analytical calculations of the equilibrium properties of the system, such as the equilibrium liquid densities, the adsorption isotherm and the interfacial tension. Based on the results of these calculations, the model parameters are then determined numerically, yielding a reasonable reproduction of the MD density profiles. The results of time-dependent simulations addressing the dynamical behaviour of the system will also be presented. It will be shown that the competition between the diffusion and hydrodynamic time scales can lead to the formation of an emulsion. We also address the main difficulties and limitations of the theory regarding quantitative modelling of surfactant assisted liquid phase separation.

Analysis of Ginzburg-Landau-type models of surfactant-assisted liquid phase separation.

In this paper diffuse interface models of surfactant-assisted liquid-liquid phase separation are addressed. We start from the generalized version of the Ginzburg-Landau free-energy-functional-based model of van der Sman and van der Graaf. First, we analyze the model in the constant surfactant approximation and show the presence of a critical point at which the interfacial tension vanishes. Then we determine the adsorption isotherms and investigate the validity range of previous results. As a key point of the work, we propose a new model of the van der Sman/van der Graaf type designed for avoiding both unwanted unphysical effects and numerical difficulties present in previous models. In order to make the model suitable for describing real systems, we determine the interfacial tension analytically more precisely and analyze it over the entire accessible surfactant load range. Emerging formulas are then validated by calculating the interfacial tension from the numerical solution of the Euler-Lagrange equations. Time-dependent simulations are also performed to illustrate the slowdown of the phase separation near the critical point and to prove that the dynamics of the phase separation is driven by the interfacial tension.

Heterogeneous nucleation of/on nanoparticles: a density functional study using the phase-field crystal model.

Crystallization of supersaturated liquids usually starts by heterogeneous nucleation. Mounting evidence shows that even homogeneous nucleation in simple liquids takes place in two steps; first a dense amorphous precursor forms, and the crystalline phase appears via heterogeneous nucleation in/on the precursor cluster. Herein, we review recent results by a simple dynamical density functional theory, the phase-field crystal model, for (precursor-mediated) homogeneous and heterogeneous nucleation of nanocrystals. It will be shown that the mismatch between the lattice constants of the nucleating crystal and the substrate plays a decisive role in determining the contact angle and nucleation barrier, which were found to be non-monotonic functions of the lattice mismatch. Time dependent studies are essential as investigations based on equilibrium properties often cannot identify the preferred nucleation pathways. Modeling of these phenomena is essential for designing materials on the basis of controlled nucleation and/or nano-patterning.

Nonlinear hydrodynamic theory of crystallization.

We present an isothermal fluctuating nonlinear hydrodynamic theory of crystallization in molecular liquids. A dynamic coarse-graining technique is used to derive the velocity field, a phenomenology which allows a direct coupling between the free energy functional of the classical density functional theory and the Navier-Stokes equation. In contrast to the Ginzburg-Landau type amplitude theories, the dynamic response to elastic deformations is described by parameter-free kinetic equations. Employing our approach to the free energy functional of the phase-field crystal model, we recover the classical spectrum for the phonons and the steady-state growth fronts. The capillary wave spectrum of the equilibrium crystal-liquid interface is in good qualitative agreement with the molecular dynamics simulations.

Heterogeneous crystal nucleation: the effect of lattice mismatch.

A simple dynamical density functional theory is used to investigate freezing of an undercooled liquid in the presence of a crystalline substrate. We find that the adsorption of the crystalline phase on the substrate, the contact angle, and the height of the nucleation barrier are nonmonotonic functions of the lattice constant of the substrate. We show that the free-growth-limited model of particle-induced freezing by Greer et al. [Acta Mater. 48, 2823 (2000)] is valid for larger nanoparticles and a small anisotropy of the interface free energy. Faceting due to the small size of the foreign particle or a high anisotropy decouples free growth from the critical size of homogeneous nuclei.

Amorphous nucleation precursor in highly nonequilibrium fluids.

Dynamical density-functional simulations reveal structural aspects of crystal nucleation in undercooled liquids: The first appearing solid is amorphous, which promotes the nucleation of bcc crystals but suppresses the appearance of the fcc and hcp phases. These findings are associated with features of the effective interaction potential deduced from the amorphous structure.

Faceting and branching in 2D crystal growth.

Using atomic scale time-dependent density functional calculations we confirm that both diffusion-controlled and diffusionless crystallization modes exist in simple 2D systems. We provide theoretical evidence that a faceted to nonfaceted transition is coupled to these crystallization modes, and faceting is governed by the local supersaturation at the fluid-crystalline interface. We also show that competing modes of crystallization have a major influence on mesopattern formation. Irregularly branched and porous structures are emerging at the crossover of the crystallization modes. The proposed branching mechanism differs essentially from dendritic fingering driven by diffusive instability.

Polymorphism, crystal nucleation and growth in the phase-field crystal model in 2D and 3D.

We apply a simple dynamical density functional theory, the phase-field crystal (PFC) model of overdamped conservative dynamics, to address polymorphism, crystal nucleation, and crystal growth in the diffusion-controlled limit. We refine the phase diagram for 3D, and determine the line free energy in 2D and the height of the nucleation barrier in 2D and 3D for homogeneous and heterogeneous nucleation by solving the respective Euler-Lagrange (EL) equations. We demonstrate that, in the PFC model, the body-centered cubic (bcc), the face-centered cubic (fcc), and the hexagonal close-packed structures (hcp) compete, while the simple cubic structure is unstable, and that phase preference can be tuned by changing the model parameters: close to the critical point the bcc structure is stable, while far from the critical point the fcc prevails, with an hcp stability domain in between. We note that with increasing distance from the critical point the equilibrium shapes vary from the sphere to specific faceted shapes: rhombic dodecahedron (bcc), truncated octahedron (fcc), and hexagonal prism (hcp). Solving the equation of motion of the PFC model supplied with conserved noise, solidification starts with the nucleation of an amorphous precursor phase, into which the stable crystalline phase nucleates. The growth rate is found to be time dependent and anisotropic; this anisotropy depends on the driving force. We show that due to the diffusion-controlled growth mechanism, which is especially relevant for crystal aggregation in colloidal systems, dendritic growth structures evolve in large-scale isothermal single-component PFC simulations. An oscillatory effective pair potential resembling those for model glass formers has been evaluated from structural data of the amorphous phase obtained by instantaneous quenching. Finally, we present results for eutectic solidification in a binary PFC model.

Crystal nucleation in the hard-sphere system revisited: a critical test of theoretical approaches.

The hard-sphere system is the best known fluid that crystallizes: the solid-liquid interfacial free energy, the equations of state, and the height of the nucleation barrier are known accurately, offering a unique possibility for a quantitative validation of nucleation theories. A recent significant downward revision of the interfacial free energy from approximately 0.61kT/sigma(2) to (0.56 +/- 0.02)kT/sigma(2) [Davidchack, R.; Morris, J. R.; Laird, B. B. J. Chem. Phys. 2006, 125, 094710] necessitates a re-evaluation of theoretical approaches to crystal nucleation. This has been carried out for the droplet model of the classical nucleation theory (CNT), the self-consistent classical theory (SCCT), a phenomenological diffuse interface theory (DIT), and single- and two-field variants of the phase field theory that rely on either the usual double-well and interpolation functions (PFT/S1 and PFT/S2, respectively) or on a Ginzburg-Landau expanded free energy that reflects the crystal symmetries (PFT/GL1 and PFT/GL2). We find that the PFT/GL1, PFT/GL2, and DIT models predict fairly accurately the height of the nucleation barrier known from Monte Carlo simulations in the volume fraction range of 0.52 < varphi < 0.54, whereas the CNT, SCCT, PFT/S1, and PFT/S2 models underestimate it significantly.

Phase field theory of interfaces and crystal nucleation in a eutectic system of fcc structure: I. Transitions in the one-phase liquid region.

The phase field theory (PFT) has been applied to predict equilibrium interfacial properties and nucleation barrier in the binary eutectic system Ag-Cu using double well and interpolation functions deduced from a Ginzburg-Landau expansion that considers fcc (face centered cubic) crystal symmetries. The temperature and composition dependent free energies of the liquid and solid phases are taken from CALculation of PHAse Diagrams-type calculations. The model parameters of PFT are fixed so as to recover an interface thickness of approximately 1 nm from molecular dynamics simulations and the interfacial free energies from the experimental dihedral angles available for the pure components. A nontrivial temperature and composition dependence for the equilibrium interfacial free energy is observed. Mapping the possible nucleation pathways, we find that the Ag and Cu rich critical fluctuations compete against each other in the neighborhood of the eutectic composition. The Tolman length is positive and shows a maximum as a function of undercooling. The PFT predictions for the critical undercooling are found to be consistent with experimental results. These results support the view that heterogeneous nucleation took place in the undercooling experiments available at present. We also present calculations using the classical droplet model [classical nucleation theory (CNT)] and a phenomenological diffuse interface theory (DIT). While the predictions of the CNT with a purely entropic interfacial free energy underestimate the critical undercooling, the DIT results appear to be in a reasonable agreement with the PFT predictions.

Phase field theory of interfaces and crystal nucleation in a eutectic system of fcc structure: II. Nucleation in the metastable liquid immiscibility region.

In the second part of our paper, we address crystal nucleation in the metastable liquid miscibility region of eutectic systems that is always present, though experimentally often inaccessible. While this situation resembles the one seen in single component crystal nucleation in the presence of a metastable vapor-liquid critical point addressed in previous works, it is more complex because of the fact that here two crystal phases of significantly different compositions may nucleate. Accordingly, at a fixed temperature below the critical point, six different types of nuclei may form: two liquid-liquid nuclei: two solid-liquid nuclei; and two types of composite nuclei, in which the crystalline core has a liquid "skirt," whose composition falls in between the compositions of the solid and the initial liquid phases, in addition to nuclei with concentric alternating composition shells of prohibitively high free energy. We discuss crystalline phase selection via exploring/identifying the possible pathways for crystal nucleation.