Renormalization-group theory for the phase-field crystal equation

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.

A unified theory of free energy functionals and applications to diffusion

SignificanceThe free energy functional is a central component of continuum dynamical models used to describe phase transitions, microstructural evolution, and pattern formation. However, despite the success of these models in many areas of physics, chemistry, and biology, the standard free energy frameworks are frequently characterized by physically opaque parameters and incorporate assumptions that are difficult to assess. Here, we introduce a mathematical formalism that provides a unifying umbrella for constructing free energy functionals. We show that Ginzburg-Landau framework is a special case of this umbrella and derive a generalization of the widely employed Cahn-Hilliard equation. More broadly, we expect the framework will also be useful for generalizing higher-order theories, establishing formal connections to microscopic physics, and coarse graining.

Multiscale analysis of crystalline defect formation in rapid solidification of pure aluminium and aluminium-copper alloys

Rapid solidification leads to unique microstructural features, where a less studied topic is the formation of various crystalline defects, including high dislocation densities, as well as gradients and splitting of the crystalline orientation. As these defects critically affect the material's mechanical properties and performance features, it is important to understand the defect formation mechanisms, and how they depend on the solidification conditions and alloying. To illuminate the formation mechanisms of the rapid solidification induced crystalline defects, we conduct a multiscale modelling analysis consisting of bond-order potential-based molecular dynamics (MD), phase field crystal-based amplitude expansion simulations, and sequentially coupled phase field-crystal plasticity simulations. The resulting dislocation densities are quantified and compared to past experiments. The atomistic approaches (MD, PFC) can be used to calibrate continuum level crystal plasticity models, and the framework adds mechanistic insights arising from the multiscale analysis. This article is part of the theme issue 'Transport phenomena in complex systems (part 2)'.

Orientation Gradients in Rapidly Solidified Pure Aluminum Thin Films: Comparison of Experiments and Phase-Field Crystal Simulations

Rapid solidification experiments on thin film aluminum samples reveal the presence of lattice orientation gradients within crystallizing grains. To study this phenomenon, a single-component phase-field crystal (PFC) model that captures the properties of solid, liquid, and vapor phases is proposed to model pure aluminium quantitatively. A coarse-grained amplitude representation of this model is used to simulate solidification in samples approaching micrometer scales. The simulations reproduce the experimentally observed orientation gradients within crystallizing grains when grown at experimentally relevant rapid quenches. We propose a causal connection between defect formation and orientation gradients.

Mesoscale Defect Motion in Binary Systems: Effects of Compositional Strain and Cottrell Atmospheres

The velocity of dislocations is derived analytically to incorporate and predict the intriguing effects induced by the preferential solute segregation and Cottrell atmospheres in both two-dimensional and three-dimensional binary systems of various crystalline symmetries. The corresponding mesoscopic description of defect dynamics is constructed through the amplitude formulation of the phase-field crystal model, which has been shown to accurately capture elasticity and plasticity in a wide variety of systems. Modifications of the Peach-Koehler force as a result of solute concentration variations and compositional stresses are presented, leading to interesting new predictions of defect motion due to effects of Cottrell atmospheres. These include the deflection of dislocation glide paths, the variation of climb speed and direction, and the change or prevention of defect annihilation, all of which play an important role in determining the fundamental behaviors of complex defect network and dynamics. The analytic results are verified by numerical simulations.

Enabling simulations of grains within a full rotation range in amplitude expansion of the phase-field crystal model

This paper introduces improvements to the amplitude expansion of the phase-field crystal model that enable the simulation of grains within a full range of orientations. The unphysical grain boundary between grains, rotated by a crystal's symmetry rotation, is removed using a combination of the auxiliary rotation field described in our previous work and an algorithm that correctly matches the complex amplitudes according to the differences in local rotation.

Tunable mosaic structures in van der Waals layered materials

Intrinsic mosaic structures composed of distinctive stacking domains separated by domain walls (DWs) show the potential to regulate many outstanding properties of van der Waals layered materials. A comprehensive simulation at the atomic scale is performed to explore how the lattice/twist mismatch and the interlayer interaction influence the mosaic configuration from the incommensurate Moiré pattern to commensurate mosaic structures by adapting a complex amplitude version of the phase field crystal method. It is found that after an incommensurate-commensurate transition occurs, the topology of the mosaic structure indicated by different domain wall (DW) patterns can be drastically changed. An experimentally observed intriguing spiral domain wall (SDW) network is revealed as result of the emergent mixed dislocation driven by minimizing the elastic and interlayer energies in the presence of both lattice and twist mismatches. The transition process from a herringbone domain wall (HBDW) network to a SDW network is also simulated, elucidated by a dislocation reaction and in good agreement with the experimental observations.

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'.

From atomistic interfaces to dendritic patterns

Transport processes around phase interfaces, together with thermodynamic properties and kinetic phenomena, control the formation of dendritic patterns. Using the thermodynamic and kinetic data of phase interfaces obtained on the atomic scale, one can analyse the formation of a single dendrite and the growth of a dendritic ensemble. This is the result of recent progress in theoretical methods and computational algorithms calculated using powerful computer clusters. Great benefits can be attained from the development of micro-, meso- and macro-levels of analysis when investigating the dynamics of interfaces, interpreting experimental data and designing the macrostructure of samples. The review and research articles in this theme issue cover the spectrum of scales (from nano- to macro-length scales) in order to exhibit recently developing trends in the theoretical analysis and computational modelling of dendrite pattern formation. Atomistic modelling, the flow effect on interface dynamics, the transition from diffusion-limited to thermally controlled growth existing at a considerable driving force, two-phase (mushy) layer formation, the growth of eutectic dendrites, the formation of a secondary dendritic network due to coalescence, computational methods, including boundary integral and phase-field methods, and experimental tests for theoretical models-all these themes are highlighted in the present issue.This article is part of the theme issue 'From atomistic interfaces to dendritic patterns'.

Controlling the energy of defects and interfaces in the amplitude expansion of the phase-field crystal model

One of the major difficulties in employing phase-field crystal (PFC) modeling and the associated amplitude (APFC) formulation is the ability to tune model parameters to match experimental quantities. In this work, we address the problem of tuning the defect core and interface energies in the APFC formulation. We show that the addition of a single term to the free-energy functional can be used to increase the solid-liquid interface and defect energies in a well-controlled fashion, without any major change to other features. The influence of the newly added term is explored in two-dimensional triangular and honeycomb structures as well as bcc and fcc lattices in three dimensions. In addition, a finite-element method (FEM) is developed for the model that incorporates a mesh refinement scheme. The combination of the FEM and mesh refinement to simulate amplitude expansion with a new energy term provides a method of controlling microscopic features such as defect and interface energies while simultaneously delivering a coarse-grained examination of the system.

Honeycomb and triangular domain wall networks in heteroepitaxial systems

A comprehensive study is presented for the influence of misfit strain, adhesion strength, and lattice symmetry on the complex Moiré patterns that form in ultrathin films of honeycomb symmetry adsorbed on compact triangular or honeycomb substrates. The method used is based on a complex Ginzburg-Landau model of the film that incorporates elastic strain energy and dislocations. The results indicate that different symmetries of the heteroepitaxial systems lead to distinct types of domain wall networks and phase transitions among various surface Moiré patterns and superstructures. More specifically, the results show a dramatic difference between the phase diagrams that emerge when a honeycomb film is adsorbed on substrates of honeycomb versus triangular symmetry. It is also shown that in the small deformation limit, the complex Ginzburg-Landau model reduces to a two-dimensional sine-Gordon free energy form. This free energy can be solved exactly for one dimensional patterns and reveals the role of domains walls and their crossings in determining the nature of the phase diagrams.

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.

Scaling of alloy interfacial properties under compositional strain

Complex morphologies and microstructures that emerge during materials growth and solidification are often determined by both equilibrium and kinetic properties of the interface and their crystalline anisotropies. However, limited knowledge is available on alloying and, particularly, compositionally generated elastic effects on these interface characteristics. Here we systematically investigate such compositional effects on the interfacial properties of an alloy model system based on a phase-field-crystal analysis, including the solid-liquid interfacial free energy, kinetic coefficient, and lattice pinning strength. Scaling relations for these interfacial quantities over various ranges of material parameters are identified and predicted. Our results indicate the important effects of couplings among mesoscopic and microscopic length scales of alloy structure and concentration, and the influence of compressive and tensile interface stresses induced by composition variations. The approach developed here provides an efficient way to systematically identify these key material properties beyond the traditional atomistic and continuum methods.

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.

Microscopic treatment of solute trapping and drag

The long wavelength limit of a recent microscopic phase-field crystal (PFC) theory of a binary alloy mixture is used to derive an analytical approximation for the segregation coefficient as a function of the interface velocity, and relate it to the two-point correlation function of the liquid and the thermodynamic properties of solid and liquid phases. Our results offer the first analytical derivation of solute segregation from a microscopic model, and support recent molecular dynamics and numerical PFC simulations. Our results also provide an independent framework, motivated from classical density functional theory, from which to elucidate the fundamental nature of solute drag, which is still highly contested in the literature.

Morphological instability of grain boundaries in two-phase coherent solids

We show both computationally and analytically that grain boundaries that exhibit shear-coupled motion become morphologically unstable in solid alloys that phase separate into coherent domains of distinct chemical compositions. We carry out simulations of continuum models demonstrating that this instability is mediated by long-range elastic interaction between compositional domains and grain boundaries. In addition, we perform a linear stability analysis that predicts the range of unstable wavelengths in good quantitative agreement with simulations. In nonlinear stages, this pattern-forming instability leads to the breakup of low-angle grain boundaries, thereby strongly impacting microstructural evolution in a wide range of phase-separating materials.

Phase-field-crystal models and mechanical equilibrium

Phase-field-crystal (PFC) models constitute a field theoretical approach to solidification, melting, and related phenomena at atomic length and diffusive time scales. One of the advantages of these models is that they naturally contain elastic excitations associated with strain in crystalline bodies. However, instabilities that are diffusively driven towards equilibrium are often orders of magnitude slower than the dynamics of the elastic excitations, and are thus not included in the standard PFC model dynamics. We derive a method to isolate the time evolution of the elastic excitations from the diffusive dynamics in the PFC approach and set up a two-stage process, in which elastic excitations are equilibrated separately. This ensures mechanical equilibrium at all times. We show concrete examples demonstrating the necessity of the separation of the elastic and diffusive time scales. In the small-deformation limit this approach is shown to agree with the theory of linear elasticity.

Phase-field-crystal model for magnetocrystalline interactions in isotropic ferromagnetic solids

An isotropic magnetoelastic phase-field-crystal model to study the relation between morphological structure and magnetic properties of pure ferromagnetic solids is introduced. Analytic calculations in two dimensions were used to determine the phase diagram and obtain the relationship between elastic strains and magnetization. Time-dependent numerical simulations in two dimensions were used to demonstrate the effect of grain boundaries on the formation of magnetic domains. It was shown that the grain boundaries act as nucleating sites for domains of reverse magnetization. Finally, we derive a relation for coercivity versus grain misorientation in the isotropic limit.

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.

Reductive renormalization of the phase-field crystal equation

It has been known for some time that singular perturbation and reductive perturbation can be unified from the renormalization-group theoretical point of view: Reductive extraction of space-time global behavior is the essence of singular perturbation methods. Reductive renormalization was proposed to make this unification practically accessible; actually, this reductive perturbation is far simpler than most reduction methods, such as the rather standard scaling expansion. However, a rather cryptic exposition of the method seems to have been the cause of some trouble. Here, an explicit demonstration of the consistency of the reductive renormalization-group procedure is given for partial differentiation equations (of a certain type, including time-evolution semigroup type equations). Then, the procedure is applied to the reduction of a phase-field crystal equation to illustrate the streamlined reduction method. We conjecture that if the original system is structurally stable, the reductive renormalization-group result and that of the original equation are diffeomorphic.

Patterning of heteroepitaxial overlayers from nano to micron scales

Thin heteroepitaxial overlayers have been proposed as templates to generate stable, self-organized nanostructures at large length scales, with a variety of important technological applications. However, modeling strain-driven self-organization is a formidable challenge due to different length scales involved. In this Letter, we present a method for predicting the patterning of ultrathin films on micron length scales with atomic resolution. We make quantitative predictions for the type of superstructures (stripes, honeycomb, triangular) and length scale of pattern formation of two metal-metal systems, Cu on Ru(0001) and Cu on Pd(111). Our findings are in excellent agreement with previous experiments and call for future experimental investigations of such systems.

Phase-field-crystal calculation of crystal-melt surface tension in binary alloys

A phase field crystal (PFC) density functional for binary mixtures is coarse grained and a formalism for calculating the simultaneous concentration, temperature, and density dependence of the surface energy anisotropy of a solid-liquid interface is developed. The methodology systematically relates bulk free energy coefficients arising from coarse graining to thermodynamic data, while gradient energy coefficients are related to molecular properties. Our coarse-grained formalism is applied to the determination of surface energy anisotropy in two-dimensional Zn-Al films, a situation relevant for quantitative phase field simulations of dendritic solidification in zinc coatings.

Effective model hierarchies for dynamic and static classical density functional theories

The origin and methodology of deriving effective model hierarchies are presented with applications to solidification of crystalline solids. In particular, it is discussed how the form of the equations of motion and the effective parameters on larger scales can be obtained from the more microscopic models. It will be shown that tying together the dynamic structure of the projection operator formalism with static classical density functional theories can lead to incomplete (mass) transport properties even though the linearized hydrodynamics on large scales is correctly reproduced. To facilitate a more natural way of binding together the dynamics of the macrovariables and classical density functional theory, a dynamic generalization of density functional theory based on the nonequilibrium generating functional is suggested.

A phase field crystal study of epitaxial island formation on nanomembranes

In this paper the liquid phase heteroepitaxial growth of two-dimensional strained islands on nanomembranes is examined via an amplitude expansion of a binary phase field crystal model. The maximum size that the islands can grow to coherently is shown to be strongly dependent on the nanomembrane thickness and to a lesser extent on the flux rate. For a large membrane-island misfit of about 10%, islands were found to be able to grow coherently with the membrane to almost twice as large as those grown on thick membranes (or infinite substrates). It was also found that when islands are growing on both sides of the membrane, strain relaxation in the membrane leads to more and less favorable growth regions. For thinner membranes this effect increases the degree of ordering of the islands, a result consistent with recent experimental findings.

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.

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.

Nonlinear elasticity of the phase-field crystal model from the renormalization group

The rotationally covariant renormalization group equations of motion for the density wave amplitudes in the phase field crystal model are shown to follow from a dynamical equation driven by an effective free energy density that we derive. We show that this free energy can be written purely as a function of the strain tensor and thence derive the corresponding equations governing the nonlinear elastic response.

An anisotropic phase-field crystal model for heterogeneous nucleation of ellipsoidal colloids

We derive a generalized model for isotropic as well as anisotropic crystal lattice systems of arbitrary Poisson ratio within the framework of the continuum phase-field crystal (PFC) approach (Elder and Grant 2004 Phys. Rev. E 70 051606). To this end we extend the simplest PFC model defined by a free energy functional, which is based upon the Swift-Hohenberg model of pattern formation (Swift and Hohenberg 1993 Phys. Rev. A 15 851) to a conservative, anisotropic Langevin equation. By studying heterogeneous nucleation of ellipsoidal colloids at a wall, we demonstrate the capacity of our approach to contribute to the more precise understanding of condensed matter systems built up from non-spherical units at the atomic scale. In particular we address the question of how (a) the orientation of the ellipsoids as well as (b) the interaction potential with the wall determine the resulting contact angle.

Thermodynamics of bcc metals in phase-field-crystal models

We examine the influence of different forms of the free-energy functionals used in the phase-field-crystal (PFC) model, and compare them with the second-order density-functional theory (DFT) of freezing, by using bcc iron as an example case. We show that there are large differences between the PFC and the DFT and it is difficult to obtain reasonable parameters for existing PFC models directly from the DFT. Therefore, we propose a way of expanding the correlation function in terms of gradients that allows us to incorporate the bulk modulus of the liquid as an additional parameter in the theory. We show that this functional reproduces reasonable values for both bulk and surface properties of bcc iron, and therefore it should be useful in modeling bcc materials. As a further demonstration, we also calculate the grain boundary energy as a function of misorientation for a symmetric tilt boundary close to the melting transition.

Molecular dynamics on diffusive time scales from the phase-field-crystal equation

We extend the phase-field-crystal model to accommodate exact atomic configurations and vacancies by requiring the order parameter to be non-negative. The resulting theory dictates the number of atoms and describes the motion of each of them. By solving the dynamical equation of the model, which is a partial differential equation, we are essentially performing molecular dynamics simulations on diffusive time scales. To illustrate this approach, we calculate the two-point correlation function of a fluid.

Comment on "Renormalization-group theory for the phase-field crystal equation"

Athreya, Goldenfeld, and Dantzig [Phys. Rev. E 74, 011601 (2006)] claim that the current implementation of the renormalization-group method neglects the proper ordering of renormalization and differentiation. Their analysis is, however, based on the wrong multiple-scales method results.

Mesoscopic and microscopic modeling of island formation in strained film epitaxy

The instability of strained films for island formation is examined through an approach incorporating both discrete microscopic details and continuum mechanics. A linear relationship between the island wave number and misfit strain is found for large strains, while only in the small strain limit is a crossover to the continuum elasticity result obtained. A universal scaling relation accommodating all range of misfit strains is identified. Our results indicate that continuum mechanics may break down even at relatively small misfit stress due to the discrete nature of crystalline surfaces.

Adaptive mesh computation of polycrystalline pattern formation using a renormalization-group reduction of the phase-field crystal model

We implement an adaptive mesh algorithm for calculating the space and time dependence of the atomic density field in microscopic material processes. Our numerical approach uses the systematic renormalization-group formulation of a phase-field crystal model of a pure material to provide the underlying equations for the complex amplitude of the atomic density field--a quantity that is spatially uniform except near topological defects, grain boundaries, and other lattice imperfections. Our algorithm employs a hybrid formulation of the amplitude equations, combining Cartesian and polar decompositions of the complex amplitude. We show that this approach leads to an acceleration by three orders of magnitude in model calculations of polycrystalline grain growth in two dimensions.