Deriving surface-energy anisotropy for phenomenological phase-field models of solidification

Incorporating Noise Quantitatively in the Phase Field Crystal Model via Capillary Fluctuation Theory

A tacit assumption underlying most phase field models of nonequilibrium phase transformations is that of scale separation. Stochastic order parameter field theories utilize noise to separate atomic-scale fluctuations from the slowly varying fields that describe microstructure patterns. The mesoscale distribution of such stochastic variables is generally assumed to follow Gaussian statistics, with their magnitude following fluctuation-dissipation relations. However, there is still much debate about how atomic-scale fluctuations map onto the mesoscale upon coarse graining of microscopic theories. This Letter studies interface fluctuations in the phase field crystal (PFC) model and proposes a self-consistent method for relating how the effective noise strength and spectral filtering of the noise in the PFC model, and similar types of microscopic models, should be defined so as to attain the spectrum of mesoscale capillary fluctuations quantitatively.

Calculating the role of composition in the anisotropy of solid-liquid interface energy using phase-field-crystal theory

This work uses Ginzburg-Landau theory derived from a recent structural phase-field-crystal model of binary alloys developed by the authors to study the roles of concentration, temperature, and pressure on the interfacial energy anisotropy of a solid-liquid front. It is found that the main contribution to the change in anisotropy with concentration arises from a change in preferred crystallographic orientation controlled by solute-dependent changes in the two-point density correlation function of a binary alloy, a mechanism that leads to such phenomena as solute-induced elastic strain and dislocation-assisted solute clustering. Our results are consistent with experimental observations in recent studies by Rappaz et al. [J. Fife, P. Di Napoli, and M. Rappaz, Metall. Mater. Trans. A 44, 5522 (2013)]. This is the first PFC work, to our knowledge, to incorporate temperature, pressure, and density into the thermodynamic description of alloys.

Phase-field-crystal model of phase and microstructural stability in driven nanocrystalline systems

We present a phase-field-crystal model for driven systems which describes competing effects between thermally activated diffusional processes and those driven by externally imposed ballistic events. The model demonstrates how the mesoscopic Enrique and Bellon [Phys. Rev. Lett. 84, 2885 (2000)] model of externally induced ballistic mixing can be incorporated into the atomistic phase-field-crystal formalism. The combination of the two approaches results in a model capable of describing the microstructural and compositional evolution of a driven system while incorporating elastoplastic effects. The model is applied to the study of grain growth in nanocrystalline materials subjected to an external driving.

Phase-field-crystal methodology for modeling of structural transformations

We introduce and characterize free-energy functionals for modeling of solids with different crystallographic symmetries within the phase-field-crystal methodology. The excess free energy responsible for the emergence of periodic phases is inspired by classical density-functional theory, but uses only a minimal description for the modes of the direct correlation function to preserve computational efficiency. We provide a detailed prescription for controlling the crystal structure and introduce parameters for changing temperature and surface energies, so that phase transformations between body-centered-cubic (bcc), face-centered-cubic (fcc), hexagonal-close-packed (hcp), and simple-cubic (sc) lattices can be studied. To illustrate the versatility of our free-energy functional, we compute the phase diagram for fcc-bcc-liquid coexistence in the temperature-density plane. We also demonstrate that our model can be extended to include hcp symmetry by dynamically simulating hcp-liquid coexistence from a seeded crystal nucleus. We further quantify the dependence of the elastic constants on the model control parameters in two and three dimensions, showing how the degree of elastic anisotropy can be tuned from the shape of the direct correlation functions.

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.

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.

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.

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.

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.

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.

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.

Phase field crystal study of deformation and plasticity in nanocrystalline materials

We introduce a modified phase field crystal (MPFC) technique that self-consistently incorporates rapid strain relaxation alongside the usual plastic deformation and multiple crystal orientations featured by the traditional phase field crystal (PFC) technique. Our MPFC formalism can be used to study a host of important phase transformation phenomena in material processing that require rapid strain relaxation. We apply the MPFC model to study elastic and plastic deformations in nanocrystalline materials, focusing on the "reverse" Hall-Petch effect. Finally, we introduce a multigrid algorithm for efficient numerical simulations of the MPFC model.