Phase-field-crystal and Swift-Hohenberg equations with fast dynamics

Pattern Formation under Deep Supercooling by Classical Density Functional-Based Approach

Solidification patterns during nonequilibrium crystallization are among the most important microstructures in the natural and technical realms. In this work, we investigate the crystal growth in deeply supercooled liquid using the classical density functional-based approaches. Our result shows that the complex amplitude expanded phase-field crystal (APFC) model containing the vacancy nonequilibrium effects proposed by us could naturally reproduce the growth front nucleation (GFN) and various nonequilibrium patterns, including the faceted growth, spherulite, symmetric and nonsymmetric dendrites among others, at the atom level. Moreover, an extraordinary microscopic columnar-to-equiaxed transition is uncovered, which is found to depend on the seed spacing and distribution. Such a phenomenon could be attributed to the combined effects of the long-wave and short-wave elastic interactions. Particularly, the columnar growth could also be predicted by an APFC model containing inertia effects, but the lattice defect type in the growing crystal is different due to the different types of short-wave interactions. Two stages are identified during the crystal growth under different undercooling, corresponding to diffusion-controlled growth and GFN-dominated growth, respectively. However, compared with the second stage, the first stage becomes too short to be noticed under the high undercooling. The distinct feature of the second stage is the dramatic increments of lattice defects, which explains the amorphous nucleation precursor in the supercooled liquid. The transition time between the two stages at different undercooling is investigated. Crystal growth of BCC structure further confirms our conclusions.

Structure diagram and dynamics of formation of hexagonal boron nitride in phase-field crystal model

The phase-field crystal (PFC-model) is a powerful tool for modelling of the crystallization in colloidal and metallic systems. In the present work, the modified hyperbolic phase-field crystal model for binary systems is presented. This model takes into account slow and fast dynamics of moving interfaces for both concentration and relative atomic number density (which were taken as order parameters). The model also includes specific mobilities for each dynamical field and correlated noise terms. The dynamics of chemical segregation with origination of mixed pseudo-hexagonal binary phase (the so-called 'triangle phase') is used as a benchmark in two spatial dimensions for the developing model. Using the free energy functional and specific lattice vectors for hexagonal crystal, the structure diagram of co-existence of liquid and three-dimensional hexagonal phase for the binary PFC-model was carried out. Parameters of the crystal lattice correspond to the hexagonal boron nitride (BN) crystal, the values of which have been taken from the literature. The paper shows the qualitative agreement between the developed structure diagram of the PFC model and the previously known equilibrium diagram for BN constructed using thermodynamic functions. This article is part of the theme issue 'Transport phenomena in complex systems (part 2)'.

Nucleation and Post-Nucleation Growth in Diffusion-Controlled and Hydrodynamic Theory of Solidification

Two-step nucleation and subsequent growth processes were investigated in the framework of the single mode phase-field crystal model combined with diffusive dynamics (corresponding to colloid suspensions) and hydrodynamical density relaxation (simple liquids). It is found that independently of dynamics, nucleation starts with the formation of solid precursor clusters that consist of domains with noncrystalline ordering (ringlike projections are seen from certain angles), and regions that have amorphous structure. Using the average bond order parameter q¯6, we distinguished amorphous, medium range crystallike order (MRCO), and crystalline local orders. We show that crystallization to the stable body-centered cubic phase is preceded by the formation of a mixture of amorphous and MRCO structures. We have determined the time dependence of the phase composition of the forming solid state. We also investigated the time/size dependence of the growth rate for solidification. The bond order analysis indicates similar structural transitions during solidification in the case of diffusive and hydrodynamic density relaxation.

Traveling waves of the solidification and melting of cubic crystal lattices

Using the phase field crystal model (PFC model), an analysis of slow and fast dynamics of solid-liquid interfaces in solidification and melting processes is presented. Dynamical regimes for cubic lattices invading metastable liquids (solidification) and liquids propagating into metastable crystals (melting) are described in terms of the evolving amplitudes of the density field. Dynamical equations are obtained for body-centered cubic (bcc) and face-centered cubic (fcc) crystal lattices in one- and two-mode approximations. A universal form of the amplitude equations is obtained for the three-dimensional dynamics for different crystal lattices and crystallographic directions. Dynamics of the amplitude's propagation for different lattices and PFC mode's approximations is qualitatively compared. The traveling-wave velocity is quantitatively compared with data of molecular dynamics simulation previously obtained by Mendelev et al. [Modell. Simul. Mater. Sci. Eng. 18, 074002 (2010)MSMEEU0965-039310.1088/0965-0393/18/7/074002] for solidification and melting of the aluminum fcc lattice.

Thermodynamics of rapid solidification and crystal growth kinetics in glass-forming alloys

Thermodynamic driving forces and growth rates in rapid solidification are analysed. Taking into account the relaxation time of the solute diffusion flux in the model equations, the present theory uses, in a first case, the deviation from local chemical equilibrium, and ergodicity breaking. The second case of ergodicity breaking may exist in crystal growth kinetics of rapidly solidifying glass-forming metals and alloys. In this case, a theoretical analysis of dendritic solidification is given for congruently melting alloys in which chemical segregation does not occur. Within this theory, a deviation from thermodynamic equilibrium is introduced for high undercoolings via gradient flow relaxation of the phase field. A comparison of the present derivations with previously verified theoretical predictions and experimental data is given. This article is part of the theme issue 'Heterogeneous materials: metastable and non- ergodic internal structures'.

A review on computational modelling of phase-transition problems

Phase-transition problems are ubiquitous in science and engineering. They have been widely studied via theory, experiments and computations. This paper reviews the main challenges associated with computational modelling of phase-transition problems, addressing both model development and numerical discretization of the resulting equations. We focus on classical phase-transition problems, including liquid-solid, gas-liquid and solid-solid transformations. Our review has a strong emphasis on the treatment of interfacial phenomena and the phase-field method. This article is part of the theme issue 'Heterogeneous materials: metastable and non-ergodic internal structures'.

Kinetic transition in the order-disorder transformation at a solid/liquid interface

Phase-field analysis for the kinetic transition in an ordered crystal structure growing from an undercooled liquid is carried out. The results are interpreted on the basis of analytical and numerical solutions of equations describing the dynamics of the phase field, the long-range order parameter as well as the atomic diffusion within the crystal/liquid interface and in the bulk crystal. As an example, the growth of a binary A₅₀B₅₀ crystal is described, and critical undercoolings at characteristic changes of growth velocity and the long-range order parameter are defined. For rapidly growing crystals, analogies and qualitative differences are found in comparison with known non-equilibrium effects, particularly solute trapping and disorder trapping. The results and model predictions are compared qualitatively with results of the theory of kinetic phase transitions (Chernov 1968 Sov. Phys. JETP26, 1182-1190) and with experimental data obtained for rapid dendritic solidification of congruently melting alloy with order-disorder transition (Hartmann et al. 2009 Europhys. Lett.87, 40007 (doi:10.1209/0295-5075/87/40007)).This article is part of the theme issue 'From atomistic interfaces to dendritic patterns'.

Travelling-wave amplitudes as solutions of the phase-field crystal equation

The dynamics of the diffuse interface between liquid and solid states is analysed. The diffuse interface is considered as an envelope of atomic density amplitudes as predicted by the phase-field crystal model (Elder et al. 2004 Phys. Rev. E70, 051605 (doi:10.1103/PhysRevE.70.051605); Elder et al. 2007 Phys. Rev. B75, 064107 (doi:10.1103/PhysRevB.75.064107)). The propagation of crystalline amplitudes into metastable liquid is described by the hyperbolic equation of an extended Allen-Cahn type (Galenko & Jou 2005 Phys. Rev. E71, 046125 (doi:10.1103/PhysRevE.71.046125)) for which the complete set of analytical travelling-wave solutions is obtained by the [Formula: see text] method (Malfliet & Hereman 1996 Phys. Scr.15, 563-568 (doi:10.1088/0031-8949/54/6/003); Wazwaz 2004 Appl. Math. Comput.154, 713-723 (doi:10.1016/S0096-3003(03)00745-8)). The general [Formula: see text] solution of travelling waves is based on the function of hyperbolic tangent. Together with its set of particular solutions, the general [Formula: see text] solution is analysed within an example of specific task about the crystal front invading metastable liquid (Galenko et al. 2015 Phys. D308, 1-10 (doi:10.1016/j.physd.2015.06.002)). The influence of the driving force on the phase-field profile, amplitude velocity and correlation length is investigated for various relaxation times of the gradient flow.This article is part of the theme issue 'From atomistic interfaces to dendritic patterns'.

Coarse-graining for fast dynamics of order parameters in the phase-field model

In standard descriptions, the master equation can be obtained by coarse-graining with the application of the hypothesis of full local thermalization that is equivalent to the local thermodynamic equilibrium. By contrast, fast transformations proceed in the absence of local equilibrium and the master equation must be obtained with the absence of thermalization. In the present work, a non-Markovian master equation leading, in specific cases of relaxation to local thermodynamic equilibrium, to hyperbolic evolution equations for a binary alloy, is derived for a system with two order parameters. One of them is a conserved order parameter related to the atomistic composition, and the other one is a non-conserved order parameter, which is related to phase field. A microscopic basis for phenomenological phase-field models of fast phase transitions, when the transition is so fast that there is not sufficient time to achieve local thermalization between two successive elementary processes in the system, is provided. In a particular case, when the relaxation to local thermalization proceeds by the exponential law, the obtained coarse-grained equations are related to the hyperbolic phase-field model. The solution of the model equations is obtained to demonstrate non-equilibrium phenomenon of solute trapping which appears in rapid growth of dendritic crystals.This article is part of the theme issue 'From atomistic interfaces to dendritic patterns'.

Long-wavelength properties of phase-field-crystal models with second-order dynamics

The phase-field-crystal (PFC) approach extends the notion of phase-field models by describing the topology of the microscopic structure of a crystalline material. One of the consequences is that local variation of the interatomic distance creates an elastic excitation. The dynamics of these excitations poses a challenge: pure diffusive dynamics cannot describe relaxation of elastic stresses that happen through phonon emission. To this end, several different models with fast dynamics have been proposed. In this article we use the amplitude expansion of the PFC model to compare the recently proposed hydrodynamic PFC amplitude model with two simpler models with fast dynamics. We compare these different models analytically and numerically. The results suggest that in order to have proper relaxation of elastic excitations, the full hydrodynamical description of the PFC amplitudes is required.

Consistent Hydrodynamics for Phase Field Crystals

We use the amplitude expansion in the phase field crystal framework to formulate an approach where the fields describing the microscopic structure of the material are coupled to a hydrodynamic velocity field. The model is shown to reduce to the well-known macroscopic theories in appropriate limits, including compressible Navier-Stokes and wave equations. Moreover, we show that the dynamics proposed allows for long wavelength phonon modes and demonstrate the theory numerically showing that the elastic excitations in the system are relaxed through phonon emission.

Coarse graining for the phase-field model of fast phase transitions

Fast phase transitions under lack of local thermalization between successive elementary steps of the physical process are treated analytically. Non-Markovian master equations are derived for fast processes, which do not have enough time to reach energy or momentum thermalization during rapid phase change or freezing of a highly nonequilibrium system. These master equations provide a further physical basis for evolution and transport equations of the phase-field model used previously in the analyses of fast phase transitions.

Unconditionally stable method and numerical solution of the hyperbolic phase-field crystal equation

The phase-field crystal model (PFC model) resolves systems on atomic length scales and diffusive time scales and lies in between standard phase-field modeling and atomistic methods. More recently a hyperbolic or modified PFC model was introduced to describe fast (propagative) and slow (diffusive) dynamics. We present a finite-element method for solving the hyperbolic PFC equation, introducing an unconditionally stable time integration algorithm. A spatial discretization is used with the traditional C^{0}-continuous Lagrange elements with quadratic shape functions. The space-time discretization of the PFC equation is second-order accurate in time and is shown analytically to be unconditionally stable. Numerical simulations are used to show a monotonic decrease of the free energy during the transition from the homogeneous state to stripes. Benchmarks on modeling patterns in two-dimensional space are carried out. The benchmarks show the applicability of the proposed algorithm for determining equilibrium states. Quantitatively, the proposed algorithm is verified for the problem of lattice parameter and velocity selection when a crystal invades a homogeneous unstable liquid.

Localized states in the conserved Swift-Hohenberg equation with cubic nonlinearity

The conserved Swift-Hohenberg equation with cubic nonlinearity provides the simplest microscopic description of the thermodynamic transition from a fluid state to a crystalline state. The resulting phase field crystal model describes a variety of spatially localized structures, in addition to different spatially extended periodic structures. The location of these structures in the temperature versus mean order parameter plane is determined using a combination of numerical continuation in one dimension and direct numerical simulation in two and three dimensions. Localized states are found in the region of thermodynamic coexistence between the homogeneous and structured phases, and may lie outside of the binodal for these states. The results are related to the phenomenon of slanted snaking but take the form of standard homoclinic snaking when the mean order parameter is plotted as a function of the chemical potential, and are expected to carry over to related models with a conserved order parameter.

Phase-field-crystal study of solute trapping

In this study we have incorporated two time scales into the phase-field-crystal model of a binary alloy to explore different solute trapping properties as a function of crystal-melt interface velocity. With only diffusive dynamics, we demonstrate that the segregation coefficient, K as a function of velocity for a binary alloy is consistent with the model of Kaplan and Aziz where K approaches unity in the limit of infinite velocity. However, with the introduction of wavelike dynamics in both the density and concentration fields, the trapping follows the kinetics proposed by Sobolev [Phys. Lett. A 199, 383 (1995)], where complete trapping occurs at a finite velocity.

Unconditionally gradient-stable computational schemes in problems of fast phase transitions

Equations of fast phase transitions, in which the phase boundaries move with velocities comparable with the atomic diffusion speed or with the speed of local structural relaxation, are analyzed. These equations have a singular perturbation due to the second derivative of the order parameter with respect to time, which appears due to phenomenologically introduced local nonequilibrium. To develop unconditionally stable computational schemes, the Eyre theorem [D. J. Eyre, unpublished] proved for the classical equations, based on hypotheses of local equilibrium, is used. An extension of the Eyre theorem for the case of equations for fast phase transitions is given. It is shown that the expansion of the free energy on contractive and expansive parts, suggested by Eyre for the classical equations of Cahn-Hilliard and Allen-Cahn, is also true for the equations of fast phase transitions. Grid approximations of these equations lead to gradient-stable algorithms with an arbitrary time step for numerical modeling, ensuring monotonic nonincrease of the free energy. Special examples demonstrating the extended Eyre theorem for fast phase transitions are considered.

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.