Derivation of the phase-field-crystal model for colloidal solidification

Scale-coupling and interface-pinning effects in the phase-field-crystal model

Effects of scale coupling between mesoscopic slowly varying envelopes of liquid-solid profile and the underlying microscopic crystalline structure are studied in the phase-field-crystal (PFC) model. Such scale coupling leads to nonadiabatic corrections to the PFC amplitude equations, the effect of which increases strongly with decreasing system temperature below the melting point. This nonadiabatic amplitude representation is further coarse-grained for the derivation of effective sharp-interface equations of motion in the limit of small but finite interface thickness. We identify a generalized form of the Gibbs-Thomson relation with the incorporation of coupling and pinning effects of the crystalline lattice structure. This generalized interface equation can be reduced to the form of a driven sine-Gordon equation with Kardar-Parisi-Zhang (KPZ) nonlinearity, and can be combined with two other dynamic equations in the sharp interface limit obeying the conservation condition of atomic number density in a liquid-solid system. A sample application to the study of crystal layer growth is given, and the corresponding analytic solutions showing lattice pinning and depinning effects and two distinct modes of continuous vs nucleated growth are presented. We also identify the universal scaling behaviors governing the properties of pinning strength, surface tension, interface kinetic coefficient, and activation energy of atomic layer growth, which accommodate all range of liquid-solid interface thicknesses and different material elastic moduli.

Particles at fluid-fluid interfaces: A new Navier-Stokes-Cahn-Hilliard surface- phase-field-crystal model

Colloid particles that are partially wetted by two immiscible fluids can become confined to fluid-fluid interfaces. At sufficiently high volume fractions, the colloids may jam and the interface may crystallize. The fluids together with the interfacial colloids form an emulsion with interesting material properties and offer an important route to new soft materials. A promising approach to simulate these emulsions was presented in Aland et al. [Phys. Fluids 23, 062103 (2011)], where a Navier-Stokes-Cahn-Hilliard model for the macroscopic two-phase fluid system was combined with a surface phase-field-crystal model for the microscopic colloidal particles along the interface. Unfortunately this model leads to spurious velocities which require very fine spatial and temporal resolutions to accurately and stably simulate. In this paper we develop an improved Navier-Stokes-Cahn-Hilliard-surface phase-field-crystal model based on the principles of mass conservation and thermodynamic consistency. To validate our approach, we derive a sharp interface model and show agreement with the improved diffuse interface model. Using simple flow configurations, we show that the new model has much better properties and does not lead to spurious velocities. Finally, we demonstrate the solid-like behavior of the crystallized interface by simulating the fall of a solid ball through a colloid-laden multiphase fluid.

Solidification fronts in supercooled liquids: how rapid fronts can lead to disordered glassy solids

We determine the speed of a crystallization (or, more generally, a solidification) front as it advances into the uniform liquid phase after the system has been quenched into the crystalline region of the phase diagram. We calculate the front speed by assuming a dynamical density functional theory (DDFT) model for the system and applying a marginal stability criterion. Our results also apply to phase field crystal (PFC) models of solidification. As the solidification front advances into the unstable liquid phase, the density profile behind the advancing front develops density modulations and the wavelength of these modulations is a dynamically chosen quantity. For shallow quenches, the selected wavelength is precisely that of the crystalline phase and so well-ordered crystalline states are formed. However, when the system is deeply quenched, we find that this wavelength can be quite different from that of the crystal, so the solidification front naturally generates disorder in the system. Significant rearrangement and aging must subsequently occur for the system to form the regular well-ordered crystal that corresponds to the free energy minimum. Additional disorder is introduced whenever a front develops from random initial conditions. We illustrate these findings with simulation results obtained using the PFC model.

Modeling the structure of liquids and crystals using one- and two-component modified phase-field crystal models

A modified phase-field crystal model in which the free energy may be minimized by an order parameter profile having isolated bumps is investigated. The phase diagram is calculated in one and two dimensions and we locate the regions where modulated and uniform phases are formed and also regions where localized states are formed. We investigate the effectiveness of the phase-field crystal model for describing fluids and crystals with defects. We further consider a two-component model and elucidate how the structure transforms from hexagonal crystalline ordering to square ordering as the concentration changes. Our conclusion contains a discussion of possible interpretations of the order parameter field.

Dynamic density functional theory of solid tumor growth: Preliminary models

Cancer is a disease that can be seen as a complex system whose dynamics and growth result from nonlinear processes coupled across wide ranges of spatio-temporal scales. The current mathematical modeling literature addresses issues at various scales but the development of theoretical methodologies capable of bridging gaps across scales needs further study. We present a new theoretical framework based on Dynamic Density Functional Theory (DDFT) extended, for the first time, to the dynamics of living tissues by accounting for cell density correlations, different cell types, phenotypes and cell birth/death processes, in order to provide a biophysically consistent description of processes across the scales. We present an application of this approach to tumor growth.

Density and concentration field description of nonperiodic structures

We propose a simple nonlocal energy functional that is suitable for the continuum field characterization of nonperiodic and localized textures. The phenomenological functional is based on the pairwise direction-dependent interaction of field gradients that are separated by a fixed distance. In an appendix, we describe the numerical minimization of our functional. On that basis, we investigate the kinetic evolution of threadlike stripe patterns that are created by the functional when we start from an initially disordered state. At later stages, we find a coarse graining that shows the same scaling behavior as was obtained for the Cahn-Hilliard equation. In fact, the Cahn-Hilliard model is contained in our characterization as a limiting case. A slight modification of our model omits this coarse graining and leads to nonperiodic stripe phases. For the latter case, we investigate the temporal evolution of the defects (end points) of the threadlike stripes. In view of actual applications of this functional, we discuss the characterization of processes observed for polymeric systems and vesicles. The statistics of the growth of the threadlike structures is compared to the case of step-growth polymerization reactions. Furthermore, we demonstrate that the functional may be applied for the study of vesicles in a continuum field description. Basic features, such as the tendency of tank treading in simple shear flows and parachute folding in pipe flows, are reproduced.

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.

Microscopic and macroscopic theories for the dynamics of polar liquid crystals

We derive and analyze the dynamic equations for polar liquid crystals in two spatial dimensions in the framework of classical dynamical density functional theory (DDFT). Translational density variations, polarization, and quadrupolar order are used as order-parameter fields. The results are critically compared with those obtained using the macroscopic approach of time-dependent Ginzburg-Landau (GL) equations for the analogous order-parameter fields. We demonstrate that, for both the microscopic DDFT and the macroscopic GL approach, the resulting dissipative dynamics can be derived from a dissipation function. We obtain microscopic expressions for all diagonal contributions and for many of the cross-coupling terms emerging from a GL approach. Thus, we establish a bridge between molecular correlations and macroscopic modeling for the dissipative dynamics of polar liquid crystals.

DDFT calibration and investigation of an anisotropic phase-field crystal model

The anisotropic phase-field crystal model recently proposed and used by Prieler et al (2009 J. Phys.: Condens. Matter 21 464110) is derived from microscopic density functional theory for anisotropic particles with fixed orientation. Its morphology diagram is also explored. In particular we have investigated the influence of anisotropy and undercooling on the process of nucleation and microstructure formation from the atomic to the microscale. To that end numerical simulations were performed varying those dimensionless parameters which represent anisotropy and undercooling in our anisotropic phase-field crystal model. The results from these numerical simulations are summarized in terms of a morphology diagram of the stable state phases. These stable phases are also investigated with respect to their kinetics and characteristic morphological features.

Polar liquid crystals in two spatial dimensions: the bridge from microscopic to macroscopic modeling

Two-dimensional polar liquid crystals have been discovered recently in monolayers of anisotropic molecules. Here, we provide a systematic theoretical description of liquid-crystalline phases for polar particles in two spatial dimensions. Starting from microscopic density functional theory, we derive a phase-field-crystal expression for the free-energy density that involves three local order-parameter fields, namely the translational density, the polarization, and the nematic order parameter. Various coupling terms between the order-parameter fields are obtained, which are in line with macroscopic considerations. Since the coupling constants are brought into connection with the molecular correlations, we establish a bridge from microscopic to macroscopic modeling. Our theory provides a starting point for further numerical calculations of the stability of polar liquid-crystalline phases and is also relevant for modeling of microswimmers, which are intrinsically polar.

Stability of liquid crystalline phases in the phase-field-crystal model

The phase-field-crystal model for liquid crystals is solved numerically in two spatial dimensions. This model is formulated with three position-dependent order parameters, namely the reduced translational density, the local nematic order parameter, and the mean local direction of the orientations. The equilibrium free-energy functional involves local powers of the order parameters up to fourth order, gradients of the order parameters up to fourth order, and different couplings between the order parameters. The stable phases of the equilibrium free-energy functional are calculated for various coupling parameters. Among the stable liquid crystalline states are the isotropic, nematic, columnar, smectic-A, and plastic crystalline phases. The plastic crystals can have triangular, square, and honeycomb lattices and exhibit orientational patterns with a complex topology involving a sublattice with topological defects. Phase diagrams were obtained by numerical minimization of the free-energy functional. Their main features are qualitatively in line with much simpler one-mode approximations for the order parameters.

Eighth-order phase-field-crystal model for two-dimensional crystallization

We present a derivation of the recently proposed eighth-order phase-field crystal model [A. Jaatinen, Phys. Rev. E 80, 031602 (2009)] for the crystallization of a solid from an undercooled melt. The model is used to study the planar growth of a two-dimensional hexagonal crystal, and the results are compared against similar results from dynamical density functional theory of Marconi and Tarazona, as well as other phase-field crystal models. We find that among the phase-field crystal models studied, the eighth-order fitting scheme gives results in good agreement with the density functional theory for both static and dynamic properties, suggesting it is an accurate and computationally efficient approximation to the density functional theory.

A phase-field-crystal approach to critical nuclei

We investigate a phase-field-crystal model for homogeneous nucleation. Instead of obtaining the time evolution of a density field towards equilibrium, we use a string method to identify saddle points in phase space. The saddle points allow us to obtain the nucleation barrier and the critical nucleus. The advantage of using the phase-field-crystal model for this task is that it can be used to resolve atomistic effects. The results obtained indicate different properties of the critical nucleus as compared with those for bulk crystals and provide a detailed description of the nucleation process.

Controlling crystal symmetries in phase-field crystal models

We investigate the possibility to control the symmetry of ordered states in phase-field crystal models by tuning nonlinear resonances. In two dimensions, we find that a state of square symmetry as well as the coexistence between squares and hexagons can be easily obtained. In contrast, it is delicate to obtain the coexistence of squares and liquid. We develop a general method for constructing free energy functionals that exhibit solid-liquid coexistence with desired crystal symmetries. As an example, we develop a free energy functional for square-liquid coexistence in two dimensions. A systematic analysis for determining the parameters of the necessary nonlinear terms is provided. The implications of our findings for simulations of materials with simple cubic symmetry are discussed.

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.

A phase-field-crystal model for liquid crystals

On the basis of static and dynamical density functional theory, a phase-field-crystal model is derived which involves both the translational density and the orientational degree of ordering as well as a local director field. The model exhibits stable isotropic, nematic, smectic A, columnar, plastic-crystalline and orientationally ordered crystalline phases. As far as the dynamics is concerned, the translational density is a conserved order parameter while the orientational ordering is non-conserved. The derived phase-field-crystal model can serve for use in efficient numerical investigations of various nonequilibrium situations in liquid crystals.

Derivation of a three-dimensional phase-field-crystal model for liquid crystals from density functional theory

Using a generalized order-parameter gradient expansion within density functional theory, we derive a phase-field-crystal model for liquid crystals composed of apolar particles in three spatial dimensions. Both the translational density and the orientational direction and ordering are included as order parameters. Different terms involving gradients in the order parameters in the resulting free-energy functional are compared to the macroscopic Ginzburg-Landau approach as well as to the hydrodynamic description for liquid crystals. Our approach provides microscopic expressions for all prefactors in terms of the particle interactions. Our phase-field-crystal model generalizes the conventional phase-field-crystal model of spherical particles to orientational degrees of freedom and can be used as a starting point to explore phase transitions and interfaces for various liquid-crystalline phases.

Phase-field-crystal dynamics for binary systems: Derivation from dynamical density functional theory, amplitude equation formalism, and applications to alloy heterostructures

The dynamics of phase field crystal (PFC) modeling is derived from dynamical density functional theory (DDFT), for both single-component and binary systems. The derivation is based on a truncation up to the three-point direct correlation functions in DDFT, and the lowest order approximation using scale analysis. The complete amplitude equation formalism for binary PFC is developed to describe the coupled dynamics of slowly varying complex amplitudes of structural profile, zeroth-mode average atomic density, and system concentration field. Effects of noise (corresponding to stochastic amplitude equations) and species-dependent atomic mobilities are also incorporated in this formalism. Results of a sample application to the study of surface segregation and interface intermixing in alloy heterostructures and strained layer growth are presented, showing the effects of different atomic sizes and mobilities of alloy components. A phenomenon of composition overshooting at the interface is found, which can be connected to the surface segregation and enrichment of one of the atomic components observed in recent experiments of alloying heterostructures.

Free energy functionals for efficient phase field crystal modeling of structural phase transformations

The phase field crystal (PFC) method is a promising technique for modeling materials with atomic resolution on mesoscopic time scales. While numerically more efficient than classical density functional theory (CDFT), its single mode free energy limits the complexity of structural transformations that can be simulated. We introduce a new PFC model inspired by CDFT, which uses a systematic construction of two-particle correlation functions that allows for a broad class of structural transformations. Our approach considers planar spacings, lattice symmetries, planar atomic densities, and atomic vibrational amplitudes in the unit cell, and parameterizes temperature and anisotropic surface energies. The power of our approach is demonstrated by two examples of structural phase transformations.

Phase-field-crystal model for fcc ordering

We develop and analyze a two-mode phase-field-crystal model to describe fcc ordering. The model is formulated by coupling two different sets of crystal density waves corresponding to <111> and <200> reciprocal lattice vectors, which are chosen to form triads so as to produce a simple free-energy landscape with coexistence of crystal and liquid phases. The feasibility of the approach is demonstrated with numerical examples of polycrystalline and (111) twin growth. We use a two-mode amplitude expansion to characterize analytically the free-energy landscape of the model, identifying parameter ranges where fcc is stable or metastable with respect to bcc. In addition, we derive analytical expressions for the elastic constants for both fcc and bcc. Those expressions show that a nonvanishing amplitude of [200] density waves is essential to obtain mechanically stable fcc crystals with a nonvanishing tetragonal shear modulus (C11-C12)/2. We determine the model parameters for specific materials by fitting the peak liquid structure factor properties and solid-density wave amplitudes following the approach developed for bcc [K.-A. Wu and A. Karma, Phys. Rev. B 76, 184107 (2007)]. This procedure yields reasonable predictions of elastic constants for both bcc Fe and fcc Ni using input parameters from molecular dynamics simulations. The application of the model to two-dimensional square lattices is also briefly examined.

Extended phase diagram of the three-dimensional phase field crystal model

We determine the phase diagram of the phase field crystal model in three dimensions by using numerical free energy minimization methods. Previously published results, based on single mode approximations, have indicated that in addition to the uniform (liquid) phase, there would be regions of stability of body-centered cubic, hexagonal and stripe phases. We find that in addition to these, there are also regions of stability of face-centered cubic and hexagonal close packed structures in this model.

Particles on curved surfaces: a dynamic approach by a phase-field-crystal model

We present a dynamic model to study ordering of particles on arbitrary curved surfaces. Thereby the particles are represented as maxima in a density field and a surface partial differential equation for the density field is solved to the minimal energy configuration. We study annihilation of dislocations within the ordered system and premelting along grain-boundary scars. The obtained minimal energy configurations on a sphere are compared with existing results and scaling laws are computed for the number of excess dislocations as a function of system size.

Amplitude expansion of the binary phase-field-crystal model

Amplitude representations of a binary phase-field-crystal model are developed for a two-dimensional triangular lattice and three-dimensional bcc and fcc crystal structures. The relationship between these amplitude equations and the standard phase-field models for binary-alloy solidification with elasticity are derived, providing an explicit connection between phase-field-crystal and phase-field models. Sample simulations of solute migration at grain boundaries, eutectic solidification, and quantum dot formation on nanomembranes are also presented.

Dynamical transitions and sliding friction of the phase-field-crystal model with pinning

We study the nonlinear driven response and sliding friction behavior of the phase-field-crystal (PFC) model with pinning including both thermal fluctuations and inertial effects. The model provides a continuous description of adsorbed layers on a substrate under the action of an external driving force at finite temperatures, allowing for both elastic and plastic deformations. We derive general stochastic dynamical equations for the particle and momentum densities including both thermal fluctuations and inertial effects. The resulting coupled equations for the PFC model are studied numerically. At sufficiently low temperatures, we find that the velocity response of an initially pinned commensurate layer shows hysteresis with dynamical melting and freezing transitions for increasing and decreasing applied forces at different critical values. The main features of the nonlinear response in the PFC model are similar to the results obtained previously with molecular dynamics simulations of particle models for adsorbed layers.

Computer simulation studies of finite-size broadening of solid-liquid interfaces: from hard spheres to nickel

Using molecular dynamics (MD) and Monte Carlo (MC) simulations interfacial properties of crystal-fluid interfaces are investigated for the hard sphere system and the one-component metallic system Ni (the latter modeled by a potential of the embedded atom type). Different local order parameters are considered to obtain order parameter profiles for systems where the crystal phase is in coexistence with the fluid phase, separated by interfaces with (100) orientation of the crystal. From these profiles, the mean-squared interfacial width w(2) is extracted as a function of system size. We rationalize the prediction of capillary wave theory that w(2) diverges logarithmically with the lateral size of the system. We show that one can estimate the interfacial stiffness [Formula: see text] from the interfacial broadening, obtaining [Formula: see text] for hard spheres and [Formula: see text] for Ni.

Solid-liquid interfacial energies and equilibrium shapes of nanocrystals

We extract the anisotropy of the solid-liquid interfacial energy of small crystals using phase field crystal simulations. The results indicate a strong dependence of the interfacial energy on the parameters in the phase field crystal model determining the position in the solid-liquid coexistence region in the phase diagram. Furthermore a size dependence of the anisotropy is shown if the crystal shape is reduced to the size of a nucleus.

Classical density functional theory: an ideal tool to study heterogeneous crystal nucleation

Density functional theory provides an ideal microscopic theory to address freezing and crystallization problems. We review the application of static density functional theory for the calculation of equilibrium phase diagrams. We also describe the dynamical extension of density functional theory for systems governed by overdamped Brownian dynamics. Applications of density functional theory to crystallization problems, in particular to heterogeneous crystal nucleation and subsequent crystal growth, are summarized. Heterogeneous nucleation at an externally imposed nucleation cluster is discussed in detail, in particular for a simple two-dimensional dipolar system. Finally the relation of dynamical density functional theory and the phase field crystal approach are outlined.