Phase field crystal study of deformation and plasticity in nanocrystalline materials

Atomic-Scale Insights into the Deformation Mechanism of the Microstructures in Precipitation-Strengthening Alloys

Clarifying the deformation behaviors of microstructures could greatly help us understand the precipitation-strengthening mechanism in alloys. However, it is still a formidable challenge to study the slow plastic deformation of alloys at the atomic scale. In this work, the phase-field crystal method was used to investigate the interactions between precipitates, grain boundary, and dislocation during the deformation processes at different degrees of lattice misfits and strain rates. The results demonstrate that the pinning effect of precipitates becomes increasingly strong with the increase of lattice misfit at relatively slow deformation with a strain rate of 10^-4. The cut regimen prevails under the interaction between coherent precipitates and dislocations. In the case of a large lattice misfit of 19.3%, the dislocations tend to move toward the incoherent phase interface and are absorbed. The deformation behavior of the precipitate-matrix phase interface was also investigated. Collaborative deformation is observed in coherent and semi-coherent interfaces, while incoherent precipitate deforms independently of the matrix grains. The faster deformations (strain rate is 10^-2) with different lattice misfits all are characterized by the generation of a large number of dislocations and vacancies. The results contribute to important insights into the fundamental issue about how the microstructures of precipitation-strengthening alloys deform collaboratively or independently under different lattice misfits and deformation rates.

A Second-Order Crank-Nicolson Leap-Frog Scheme for the Modified Phase Field Crystal Model with Long-Range Interaction

In this paper, we construct a fully discrete and decoupled Crank-Nicolson Leap-Frog (CNLF) scheme for solving the modified phase field crystal model (MPFC) with long-range interaction. The idea of CNLF is to treat stiff terms implicity with Crank-Nicolson and to treat non-stiff terms explicitly with Leap-Frog. In addition, the scalar auxiliary variable (SAV) method is used to allow explicit treatment of the nonlinear potential, then, these technique combines with CNLF can lead to the highly efficient, fully decoupled and linear numerical scheme with constant coefficients at each time step. Furthermore, the Fourier spectral method is used for the spatial discretization. Finally, we show that the CNLF scheme is fully discrete, second-order decoupled and unconditionally stable. Ample numerical experiments in 2D and 3D are provided to demonstrate the accuracy, efficiency, and stability of the proposed method.

Control of phase ordering and elastic properties in phase field crystals through three-point direct correlation

Effects of three-point direct correlation on properties of the phase field crystal (PFC) modeling are examined for the control of various ordered and disordered phases and their coexistence in both three-dimensional and two-dimensional systems. Such effects are manifested via the corresponding gradient nonlinearity in the PFC free-energy functional that is derived from classical density functional theory. Their significant impacts on the stability regimes of ordered phases, phase diagrams, and elastic properties of the system, as compared to those of the original PFC model, are revealed through systematic analyses and simulations. The nontrivial contribution from three-point direct correlation leads to the variation of the critical point of order-disorder transition to which all the phase boundaries in the temperature-density phase diagram converge. It also enables the variation and control of system elastic constants over a substantial range as needed in modeling different types of materials with the same crystalline structure but different elastic properties. The capability of this PFC approach in modeling both solid and soft matter systems is further demonstrated through the effect of three-point correlation on controlling the vapor-liquid-solid coexistence and transitions for body-centered cubic phase and on achieving the liquid-stripe or liquid-lamellar phase coexistence. All these provide a valuable and efficient method for the study of structural ordering and evolution in various types of material systems.

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.

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

Mechanical relaxation and fracture of phase field crystals

A computational method is developed for the study of mechanical response and fracture behavior of phase field crystals (PFC), to overcome a limitation of the PFC dynamics which lacks an effective mechanism for describing fast mechanical relaxation of the material system. The method is based on a simple interpolation scheme for PFC (IPFC) making use of a condition of the displacement field to satisfy local elastic equilibration, while preserving key characteristics of the original PFC model. We conduct a systematic study on the mechanical properties of a sample nanoribbon system with honeycomb lattice symmetry subjected to uniaxial tension, for numerical validation of the IPFC scheme and the comparison with the original PFC and modified PFC methods. Results of mechanical response, in both elasticity and fracture regimes, show the advantage and efficiency of the IPFC method across different system sizes and applied strain rates, due to its effective process of mechanical equilibration. A brittle fracture behavior is obtained in IPFC calculations, where effects of system temperature and chirality on the fracture strength and Young's modulus are also identified, with results agreeing with those found in previous atomistic simulations of graphene. The IPFC scheme developed here is generic and applicable to the mechanical studies using different types of PFC free-energy functionals designed for various material systems.

A Study of Strain-Driven Nucleation and Extension of Deformed Grain: Phase Field Crystal and Continuum Modeling

The phase-field-crystal (PFC) method is used to investigate migration of grain boundary dislocation and dynamic of strain-driven nucleation and growth of deformed grain in two dimensions. The simulated results show that the deformed grain nucleates through forming a gap with higher strain energy between the two sub-grain boundaries (SGB) which is split from grain boundary (GB) under applied biaxial strain, and results in the formation of high-density ensembles of cooperative dislocation movement (CDM) that is capable of plastic flow localization (deformed band), which is related to the change of the crystal lattice orientation due to instability of the orientation. The deformed grain stores the strain energy through collective climbing of the dislocation, as well as changing the orientation of the original grain. The deformed grain growth (DGG) is such that the higher strain energy region extends to the lower strain energy region, and its area increase is proportional to the time square. The rule of the time square of the DGG can also be deduced by establishing the dynamic equation of the dislocation of the strain-driven SGB. The copper metal is taken as an example of the calculation, and the obtained result is a good agreement with that of the experiment.

Phase field crystal simulation of stress induced localized solid-state amorphization in nanocrystalline materials

In this work, the phase field crystal (PFC) method is used to study the localized solid-state amorphization (SSA) and its dynamic transformation process in polycrystalline materials under the uniaxial tensile deformation with different factors. The impacts of these factors, including strain rates, temperatures and grain sizes, are analyzed. Kinetically, the ultra-high strain rate causes the lattice to be seriously distorted and the grain to gradually collapse, so the dislocation density rises remarkably. Therefore, localized SSA occurs. Thermodynamically, as high temperature increases the activation energy, the atoms are active and prefer to leave the original position, which induce atom rearrangement. Furthermore, small grain size increases the percentage of grain boundary and the interface free energy of the system. As a result, Helmholtz free energy increases. The dislocations and Helmholtz free energy act as the seed and driving force for the process of the localized SSA. Also, the critical diffusion-time step and the percentage of amorphous region areas are calculated. Through this work, the PFC method is proved to be an effective means to study localized SSA under uniaxial tensile deformation.

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.

Phase-field crystal model for a diamond-cubic structure

We present a structural phase-field crystal model [M. Greenwood et al., Phys. Rev. Lett. 105, 045702 (2010)] that yields a stable dc structure. The stabilization of a dc structure is accomplished by constructing a two-body direct correlation function (DCF) approximated by a combination of two Gaussian functions in Fourier space. A phase diagram containing a dc-liquid phase coexistence region is calculated for this model. We examine the energies of solid-liquid interfaces with normals along the [100], [110], and [111] directions. The dependence of the interfacial energy on a temperature parameter, which controls the heights of the peaks in the two-body DCF, is described by a Gaussian function. Furthermore, the dependence of the interfacial energy on the peak widths of the two-body DCF, which controls the excess energy associated with interfaces, defects, and strain, is described by an inverse power law. These relationships can be used to parametrize the phase-field crystal model for the dc structure to match solid-liquid interfacial energies to those measured experimentally or calculated from atomistic simulations.

Intermittent dislocation density fluctuations in crystal plasticity from a phase-field crystal model

Plastic deformation mediated by collective dislocation dynamics is investigated in the two-dimensional phase-field crystal model of sheared single crystals. We find that intermittent fluctuations in the dislocation population number accompany bursts in the plastic strain-rate fluctuations. Dislocation number fluctuations exhibit a power-law spectral density 1/f2 at high frequencies f. The probability distribution of number fluctuations becomes bimodal at low driving rates corresponding to a scenario where low density of defects alternates at irregular times with high populations of defects. We propose a simple stochastic model of dislocation reaction kinetics that is able to capture these statistical properties of the dislocation density fluctuations as a function of shear rate.

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.

Active crystals and their stability

A recently introduced active phase field crystal model describes the formation of ordered resting and traveling crystals in systems of self-propelled particles. Increasing the active drive, a resting crystal can be forced to perform collectively ordered migration as a single traveling object. We demonstrate here that these ordered migrating structures are linearly stable. In other words, during migration, the single-crystalline texture together with the globally ordered collective motion is preserved even on large length scales. Furthermore, we consider self-propelled particles on a substrate that are surrounded by a thin fluid film. We find that in this case the resulting hydrodynamic interactions can destabilize the order.

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.

Phase-field crystal model with a vapor phase

Phase-field crystal (PFC) models are able to resolve atomic length scale features of materials during temporal evolution over diffusive time scales. Traditional PFC models contain solid and liquid phases, however many important materials processing phenomena involve a vapor phase as well. In this work, we add a vapor phase to an existing PFC model and show realistic interfacial phenomena near the triple point temperature. For example, the PFC model exhibits density oscillations at liquid-vapor interfaces that compare favorably to data available for interfaces in metallic systems from both experiment and molecular-dynamics simulations. We also quantify the anisotropic solid-vapor surface energy for a two-dimensional PFC hexagonal crystal and find well-defined step energies from measurements on the faceted interfaces. Additionally, the strain field beneath a stepped interface is characterized and shown to qualitatively reproduce predictions from continuum models, simulations, and experimental data. Finally, we examine the dynamic case of step-flow growth of a crystal into a supersaturated vapor phase. The ability to model such a wide range of surface and bulk defects makes this PFC model a useful tool to study processing techniques such as chemical vapor deposition or vapor-liquid-solid growth of nanowires.

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.

Traveling and resting crystals in active systems

A microscopic field theory for crystallization in active systems is proposed which unifies the phase-field-crystal model of freezing with the Toner-Tu theory for self-propelled particles. A wealth of different active crystalline states are predicted and characterized. In particular, for increasing strength of self-propulsion, a transition from a resting crystal to a traveling crystalline state is found where the particles migrate collectively while keeping their crystalline order. Our predictions, which are verifiable in experiments and in particle-resolved computer simulations, provide a starting point for the design of new active materials.

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.

Classical density functional theory and the phase-field crystal method using a rational function to describe the two-body direct correlation function

We introduce a new approach to represent a two-body direct correlation function (DCF) in order to alleviate the computational demand of classical density functional theory (CDFT) and enhance the predictive capability of the phase-field crystal (PFC) method. The approach utilizes a rational function fit (RFF) to approximate the two-body DCF in Fourier space. We use the RFF to show that short-wavelength contributions of the two-body DCF play an important role in determining the thermodynamic properties of materials. We further show that using the RFF to empirically parametrize the two-body DCF allows us to obtain the thermodynamic properties of solids and liquids that agree with the results of CDFT simulations with the full two-body DCF without incurring significant computational costs. In addition, the RFF can also be used to improve the representation of the two-body DCF in the PFC method. Last, the RFF allows for a real-space reformulation of the CDFT and PFC method, which enables descriptions of nonperiodic systems and the use of nonuniform and adaptive grids.

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.

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.