Michael Greenwood

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.

MULTISCALE MODELING OF SOLIDIFICATION: PHASE-FIELD METHODS TO ADAPTIVE MESH REFINEMENT

We review the use of phase field methods in solidification modeling, describing their fundamental connection to the physics of phase transformations. The inherent challenges associated with simulating phase field models across multiple length and time scales are discussed, as well as how these challenges have been addressed in recent years. Specifically, we discuss new asymptotic analysis methods that enable phase field equations to emulate the sharp interface limit even in the case of quite diffuse phase-field interfaces, an aspect that greatly reduces computation times. We then review recent dynamic adaptive mesh refinement algorithms that have enabled a dramatic increase in the scale of microstructures that can be simulated using phase-field models, at significantly reduced simulation times. Combined with new methods of asymptotic analysis, the adaptive mesh approach provides a truly multi-scale capability for simulating solidification microstructures from nanometers up to centimeters. Finally, we present recent results on 2D and 3D dendritic growth and dendritic spacing selection, which have been made using phase-field models solved with adaptive mesh refinement.

Crossover scaling of wavelength selection in directional solidification of binary alloys

We simulate cellular and dendritic growth in directional solidification in dilute binary alloys using a phase-field model solved with adaptive-mesh refinement. The spacing of primary branches is examined for a wide range of thermal gradients and alloy compositions and is found to undergo a maximum as a function of pulling velocity, in agreement with experimental observations. We demonstrate that wavelength selection is unambiguously described by a nontrivial crossover scaling function from the emergence of cellular growth to the onset of dendritic fingers. This result is further validated using published experimental data, which obeys the same scaling function.

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.

Multicomponent phase-field crystal model for structural transformations in metal alloys

We present a new phase field crystal model for structural transformations in multi-component alloys. The formalism builds upon the two-point correlation kernel developed in Greenwood et al. for describing structural transformations in pure materials [Phys. Rev. Lett. 105, 045702 (2010)]. We introduce an effective twopoint correlation function for multi-component alloys that uses the local species concentrations to interpolate between different crystal structures. A simplified version of the model is derived for the particular case of threecomponent (ternary) alloys, and its equilibrium properties are demonstrated. Dynamical equations of motion for the density and multiple species concentration fields are derived, and the robustness of the model is illustrated with examples of complex microstructure evolution in dendritic solidification and solid-state precipitation.

Morphology of monolayer films on quasicrystalline surfaces from the phase field crystal model

We present a computational study of the morphology of adsorbed monolayers on quasicrystalline surfaces with five- and seven-fold symmetry. The phase field crystal model is employed to first simulate the growth of the quasicrystal surfaces and then to elastically couple a two-dimensional film to the substrate. We find several distinct pseudomorphic phases that depend on the position of adsorption sites as well as the strength of the monolayer/substrate interaction, and quantify them by computing local order parameters. In qualitative agreement with recent experiments using colloids in quasiperiodic light fields, we find that the formation of quasicrystalline order is greatly inhibited on seven-fold surfaces.