Phase field crystal models with applications to laser deposition: A review.

In this article, we address the application of phase field crystal (PFC) theory, a hybrid atomistic-continuum approach, for modeling nanostructure kinetics encountered in laser deposition. We first provide an overview of the PFC methodology, highlighting recent advances to incorporate phononic and heat transport mechanisms. To simulate laser heating, energy is deposited onto a number of polycrystalline, two-dimensional samples through the application of initial stochastic fluctuations. We first demonstrate the ability of the model to simulate plasticity and recrystallization events that follow laser heating in the isothermal limit. Importantly, we also show that sufficient kinetic energy can cause voiding, which serves to suppress shock propagation. We subsequently employ a newly developed thermo-density PFC theory, coined thermal field crystal (TFC), to investigate laser heating of polycrystalline samples under non-isothermal conditions. We observe that the latent heat of transition associated with ordering can lead to long lasting metastable structures and defects, with a healing rate linked to the thermal diffusion. Finally, we illustrate that the lattice temperature simulated by the TFC model is in qualitative agreement with predictions of conventional electron-phonon two-temperature models. We expect that our new TFC formalism can be useful for predicting transient structures that result from rapid laser heating and re-solidification processes.

In situ study of Au nanoparticle formation in a mechanochemical-aging-based method.

As we strive to perform chemical transformations in a more sustainable fashion, enabling solid-state reactions through mechanochemistry has emerged as a highly successful approach. Due to the wide-ranging applications of gold nanoparticles (AuNPs), mechanochemical strategies have already been employed for their synthesis. However, the underlying processes surrounding gold salt reduction, nucleation and growth of AuNPs in the solid state are yet to be understood. Here, we present a mechanically activated aging synthesis of AuNPs, through a solid-state Turkevich reaction. Solid reactants are only briefly exposed to input of mechanical energy before being aged statically over a period of six weeks at different temperatures. This system offers an excellent opportunity for an in situ analysis of both reduction and nanoparticle formation processes. During the aging period, the reaction was monitored using a combination of X-ray photoelectron spectroscopy, diffuse reflectance spectroscopy, powder X-ray diffraction and transmission electron microscopy, to gain meaningful insights into the mechanisms of solid-state formation of gold nanoparticles. The acquired data allowed for the establishment of the first kinetic model for solid-state nanoparticle formation.

Direct observation by high resolution transmission electron microscopy of gold(III) particle transformation during aging reduction reaction.

We use a gold nanoparticle synthesis as a model system to study the morphological and compositional changes in gold(III) precursor particles, while reduction is taking place during aging after mechanical activation. Scanning transmission electron microscopy coupled with a high-angle annular dark field detector revealed the nanoscale changes in particle morphology, while electron energy loss spectroscopy mapped the changes in the chemical landscape during the reduction process. Tracking a specific region of interest on the sample grid allowed for comparisons to be made of the same particles across a two day monitoring period. High-angle annular dark field images permitted the visualization of particle size reduction of the gold salt while electron energy loss spectroscopy captured the surprising mobility of the lighter chlorine and sodium ions in a solid matrix during the reduction process. This system offers unique insight into precursor particle reactivity in the solid phase, which is relevant for many mechanochemical and aging-based reactions.

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.

Quantitative phase field simulations of polycrystalline solidification using a vector order parameter.

A vector order parameter phase field model derived from a grand potential functional is presented as an alternative approach for modeling polycrystalline solidification of alloys. In this approach, the grand potential density is designed to contain a discrete set of finite wells, a feature that naturally allows for the growth and controlled interaction of multiple grains using a single vector field. We verify that dendritic solidification in binary alloys follows the well-established quantitative behavior in the thin interface limit. In addition, it is shown that grain boundary energy and solute back-diffusion are quantitatively consistent with earlier theoretical work, with grain boundary energy being controlled through a simple solid-solid interaction parameter. Moreover, when considering polycrystalline aggregates and their coarsening, we show that the kinetics follow the expected parabolic growth law. Finally, we demonstrate how this vector order parameter model can be used to describe nucleation in polycrystalline systems via thermal fluctuations of the vector order parameter, a process that cannot be treated consistently with multiphase or multi-order-parameter based phase field models. The presented vector order parameter model serves as a practical and efficient computational tool for simulating polycrystalline materials. We also discuss the extension of the order parameter to higher dimensions as a simple method for modeling multiple solid phases.

Controlling solid-liquid interfacial energy anisotropy through the isotropic liquid.

Although the anisotropy of the solid-liquid interfacial free energy for most alloy systems is very small, it plays a crucial role in the growth rate, morphology and crystallographic growth direction of dendrites. Previous work posited a dendrite orientation transition via compositional additions. In this work we examine experimentally the change in dendrite growth behaviour in the Al-Sm (Samarium) system as a function of solute concentration and study its interfacial properties using molecular dynamics simulations. We observe a dendrite growth direction which changes from [Formula: see text] to [Formula: see text] as Sm content increases. The observed change in dendrite orientation is consistent with the simulation results for the variation of the interfacial free energy anisotropy and thus provides definitive confirmation of a conjecture in previous works. In addition, our results provide physical insight into the atomic structural origin of the concentration dependent anisotropy, and deepen our fundamental understanding of solid-liquid interfaces in binary alloys.

Classical nucleation theory in the phase-field crystal model.

A full understanding of polycrystalline materials requires studying the process of nucleation, a thermally activated phase transition that typically occurs at atomistic scales. The numerical modeling of this process is problematic for traditional numerical techniques: commonly used phase-field methods' resolution does not extend to the atomic scales at which nucleation takes places, while atomistic methods such as molecular dynamics are incapable of scaling to the mesoscale regime where late-stage growth and structure formation takes place following earlier nucleation. Consequently, it is of interest to examine nucleation in the more recently proposed phase-field crystal (PFC) model, which attempts to bridge the atomic and mesoscale regimes in microstructure simulations. In this work, we numerically calculate homogeneous liquid-to-solid nucleation rates and incubation times in the simplest version of the PFC model, for various parameter choices. We show that the model naturally exhibits qualitative agreement with the predictions of classical nucleation theory (CNT) despite a lack of some explicit atomistic features presumed in CNT. We also examine the early appearance of lattice structure in nucleating grains, finding disagreement with some basic assumptions of CNT. We then argue that a quantitatively correct nucleation theory for the PFC model would require extending CNT to a multivariable theory.

A phase-field-crystal alloy model for late-stage solidification studies involving the interaction of solid, liquid and gas phases.

We present a multiphase binary alloy phase-field-crystal model. By introducing density difference between solid and liquid into a previous alloy model, this new fusion leads to a practical tool that can be used to investigate formation of defects in late-stage alloy solidification. It is shown that this model can qualitatively capture the liquid pressure drop due to solidification shrinkage in confined geometry. With an inherited gas phase from a previous multiphase model, cavitation of liquid from shrinkage-induced pressure is also included in this framework. As a unique model that has both solute concentration and pressure-induced liquid cavitation, it also captures a modified Scheil-Gulliver-type segregation behaviour due to cavitation. Simulation of inter-dendritic channel solidification using this model demonstrates a strong cooling rate dependence of the resulting microstructure.This article is part of the theme issue 'From atomistic interfaces to dendritic patterns'.

Nanoscale interfacial defect shedding in a growing nematic droplet.

Interfacial defect shedding is the most recent known mechanism for defect formation in a thermally driven isotropic-to-nematic phase transition. It manifests in nematic-isotropic interfaces going through an anchoring switch. Numerical computations in planar geometry established that a growing nematic droplet can undergo interfacial defect shedding, nucleating interfacial defect structures that shed into the bulk as +1/2 point defects. By extending the study of interfacial defect shedding in a growing nematic droplet to larger length and time scales, and to three dimensions, we unveil an oscillatory growth mode involving shape and anchoring transitions that results in a controllable regular distributions of point defects in planar geometry, and complex structures of disclination lines in three dimensions.

Analysis of thermodiffusive cellular instabilities in continuum combustion fronts.

We explore numerically the morphological patterns of thermodiffusive instabilities in combustion fronts with a continuum fuel source, within a range of Lewis numbers and ignition temperatures, focusing on the cellular regime. For this purpose, we generalize the recent model of Brailovsky et al. to include distinct process kinetics and reactant heterogeneity. The generalized model is derived analytically and validated with other established models in the limit of infinite Lewis number for zero-order and first-order kinetics. Cellular and dendritic instabilities are found at low Lewis numbers. These are studied using a dynamic adaptive mesh refinement technique that allows very large computational domains, thus allowing us to reduce finite-size effects that can affect or even preclude the emergence of these patterns. Our numerical linear stability analysis is consistent with the analytical results of Brailovsky et al. The distinct types of dynamics found in the vicinity of the critical Lewis number, ranging from steady-state cells to continued tip splitting and cell merging, are well described within the framework of thermodiffusive instabilities and are consistent with previous numerical studies. These types of dynamics are classified as "quasilinear" and characterized by low-amplitude cells that may be strongly affected by the mode selection mechanism and growth prescribed by the linear theory. Below this range of Lewis number, highly nonlinear effects become prominent and large-amplitude, complex cellular and seaweed dendritic morphologies emerge.

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.

Phase-Field-Crystal Model for Electromigration in Metal Interconnects.

We propose an atomistic model of electromigration (EM) in metals based on a recently developed phase-field-crystal (PFC) technique. By coupling the PFC model's atomic density field with an applied electric field through the EM effective charge parameter, EM is successfully captured on diffusive time scales. Our framework reproduces the well-established EM phenomena known as Black's equation and the Blech effect, and also naturally captures commonly observed phenomena such as void nucleation and migration in bulk crystals. A resistivity dipole field arising from electron scattering on void surfaces is shown to contribute significantly to void migration velocity. With an intrinsic time scale set by atomic diffusion rather than atomic oscillations or hopping events, as in conventional atomistic methods, our theoretical approach makes it possible to investigate EM-induced circuit failure at atomic spatial resolution and experimentally relevant time scales.

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.

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.

New density functional approach for solid-liquid-vapor transitions in pure materials.

A new phase field crystal (PFC) type theory is presented, which accounts for the full spectrum of solid-liquid-vapor phase transitions within the framework of a single density order parameter. Its equilibrium properties show the most quantitative features to date in PFC modeling of pure substances, and full consistency with thermodynamics in pressure-volume-temperature space is demonstrated. A method to control either the volume or the pressure of the system is also introduced. Nonequilibrium simulations show that 2- and 3-phase growth of solid, vapor, and liquid can be achieved, while our formalism also allows for a full range of pressure-induced transformations. This model opens up a new window for the study of pressure driven interactions of condensed phases with vapor, an experimentally relevant paradigm previously missing from phase field crystal theories.

Self-assembly via branching morphologies in nematic liquid-crystal nanocomposites.

We demonstrate that the morphological diversity in liquid-crystal hybrid systems is much richer than previously anticipated. More importantly, we reveal the existence of a dual mechanism for self-assembly of nanoparticles via morphological instabilities at phase boundaries. Using numerical simulations, we study the growth of isolated nematic droplets in an isotropic liquid crystal (LC) doped with nanoparticles (NPs) and provide insight into the nature of microstructure evolution in LC hybrids. Our work expands the numerically accessible time and length scales in these systems, capturing morphologies which develop under the competition of nonequilibrium elastic interactions, diffusive instabilities mediated by NP transport, and the anisotropy of the nematic field. By mapping nematic morphologies, we also propose a methodology for estimating various important LC material parameters that are difficult to obtain experimentally.

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.

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.

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.