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

Deriving phase field crystal theory from dynamical density functional theory: Consequences of the approximations

Phase field crystal (PFC) theory is extensively used for modeling the phase behavior, structure, thermodynamics, and other related properties of solids. PFC theory can be derived from dynamical density functional theory (DDFT) via a sequence of approximations. Here, we carefully identify all of these approximations and explain the consequences of each. One approximation that is made in standard derivations is to neglect a term of form ∇·[n∇Ln], where n is the scaled density profile and L is a linear operator. We show that this term makes a significant contribution to the stability of the crystal, and that dropping this term from the theory forces another approximation, that of replacing the logarithmic term from the ideal gas contribution to the free energy with its truncated Taylor expansion, to yield a polynomial in n. However, the consequences of doing this are (i) the presence of an additional spinodal in the phase diagram, so the liquid is predicted first to freeze and then to melt again as the density is increased; and (ii) other periodic structures, such as stripes, are erroneously predicted to be thermodynamic equilibrium structures. In general, L consists of a nonlocal convolution involving the pair direct correlation function. A second approximation sometimes made in deriving PFC theory is to replace L with a gradient expansion involving derivatives. We show that this leads to the possibility of the density going to zero, with its logarithm going to -∞ while being balanced by the fourth derivative of the density going to +∞. This subtle singularity leads to solutions failing to exist above a certain value of the average density. We illustrate all of these conclusions with results for a particularly simple model two-dimensional fluid, the generalized exponential model of index 4 (GEM-4), chosen because a DDFT is known to be accurate for this model. The consequences of the subsequent PFC approximations can then be examined. These include the phase diagram being both qualitatively incorrect, in that it has a stripe phase, and quantitatively incorrect (by orders of magnitude) regarding the properties of the crystal phase. Thus, although PFC models are very successful as phenomenological models of crystallization, we find it impossible to derive the PFC as a theory for the (scaled) density distribution when starting from an accurate DDFT, without introducing spurious artifacts. However, we find that making a simple one-mode approximation for the logarithm of the density distribution lnρ(x) rather than for ρ(x) is surprisingly accurate. This approach gives a tantalizing hint that accurate PFC-type theories may instead be derived as theories for the field lnρ(x), rather than for the density profile itself.

Interfacial free energy adjustable phase field crystal model for homogeneous nucleation

To describe the homogeneous nucleation process, an interfacial free energy adjustable phase-field crystal model (IPFC) was proposed by reconstructing the energy functional of the original phase field crystal (PFC) methodology. Compared with the original PFC model, the additional interface term in the IPFC model effectively can adjust the magnitude of the interfacial free energy, but does not affect the equilibrium phase diagram and the interfacial energy anisotropy. The IPFC model overcame the limitation that the interfacial free energy of the original PFC model is much less than the theoretical results. Using the IPFC model, we investigated some basic issues in homogeneous nucleation. From the viewpoint of simulation, we proceeded with an in situ observation of the process of cluster fluctuation and obtained quite similar snapshots to colloidal crystallization experiments. We also counted the size distribution of crystal-like clusters and the nucleation rate. Our simulations show that the size distribution is independent of the evolution time, and the nucleation rate remains constant after a period of relaxation, which are consistent with experimental observations. The linear relation between logarithmic nucleation rate and reciprocal driving force also conforms to the steady state nucleation theory.

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.