Microscopic dynamics in the liquid Li-Na alloy: an ab initio molecular dynamics study

Orbital-Free Density Functional Theory: An Attractive Electronic Structure Method for Large-Scale First-Principles Simulations

Kohn-Sham Density Functional Theory (KSDFT) is the most widely used electronic structure method in chemistry, physics, and materials science, with thousands of calculations cited annually. This ubiquity is rooted in the favorable accuracy vs cost balance of KSDFT. Nonetheless, the ambitions and expectations of researchers for use of KSDFT in predictive simulations of large, complicated molecular systems are confronted with an intrinsic computational cost-scaling challenge. Particularly evident in the context of first-principles molecular dynamics, the challenge is the high cost-scaling associated with the computation of the Kohn-Sham orbitals. Orbital-free DFT (OFDFT), as the name suggests, circumvents entirely the explicit use of those orbitals. Without them, the structural and algorithmic complexity of KSDFT simplifies dramatically and near-linear scaling with system size irrespective of system state is achievable. Thus, much larger system sizes and longer simulation time scales (compared to conventional KSDFT) become accessible; hence, new chemical phenomena and new materials can be explored. In this review, we introduce the historical contexts of OFDFT, its theoretical basis, and the challenge of realizing its promise via approximate kinetic energy density functionals (KEDFs). We review recent progress on that challenge for an array of KEDFs, such as one-point, two-point, and machine-learnt, as well as some less explored forms. We emphasize use of exact constraints and the inevitability of design choices. Then, we survey the associated numerical techniques and implemented algorithms specific to OFDFT. We conclude with an illustrative sample of applications to showcase the power of OFDFT in materials science, chemistry, and physics.

Static, dynamic and electronic properties of expanded fluid mercury in the metal-nonmetal transition range. An ab initio study

Fluid Hg undergoes a metal-nonmetal (M-NM) transition when expanded toward a density of around 9 g cm(-3). We have performed ab initio molecular dynamics simulations for several thermodynamic states around the M-NM transition range and the associated static, dynamic and electronic properties have been analyzed. The calculated static structure shows a good agreement with the available experimental data. It is found that the volume expansion decreases the number of nearest neighbors from 10 (near the triple point) to around 8 at the M-NM transition region. Moreover, these neighbors are arranged into two subshells and the decrease in the number of neighbors occurs in the inner subshell. The calculated dynamic structure factors agree fairly well with their experimental counterparts obtained by inelastic x-ray scattering experiments, which display inelastic side peaks. The derived dispersion relation exhibits some positive dispersion for all the states, although its value around the M-NM transition region is not as marked as suggested by the experiment. We have also calculated the electronic density of states, which shows the appearance of a gap at a density of around 8.3 g cm(-3).

Structural and dynamical properties of liquid Mg. An orbital-free molecular dynamics study

Several static and dynamic properties of liquid magnesium near melting have been evaluated by the orbital-free ab initio molecular dynamics method. The calculated static structure shows good agreement with recent experimental data, including an asymmetric second peak in the structure factor which has been linked to the existence of an important icosahedral short-range order in the liquid. As for the dynamic structure, we obtain collective density excitations with an associated dispersion relation which closely follows recent experimental results. Accurate estimates have also been obtained for several transport coefficients.

Competing Thermal Activation Mechanisms in the Meltinglike Transition of NaN (N = 135−147) Clusters

The meltinglike transition in unsupported icosahedral Na(N)() clusters, with N = 135-147, has been studied by isokinetic molecular dynamics simulations based on an orbital-free version of density functional theory. A maximum in the melting temperature, T(m), is obtained for Na141, while the latent heat, deltaE, and entropy of melting, deltaS, are maximal for Na147. These observations are in close agreement with calorimetric experiments on N clusters. The size evolution of deltaS is rationalized by the emergence of important premelting effects associated with the diffusive motion of atomic vacancies at the cluster surface. The precise location of the maximum in T(m) is explained in terms of two different thermally activated structural instability mechanisms which trigger the meltinglike transition in the size ranges N = 135-141 and N = 141-147, respectively.

Anomalous size dependence in the melting temperatures of free sodium clusters: an explanation for the calorimetry experiments

The meltinglike transition in unsupported Na(N) clusters (N = 55, 92, 147, 181, 189, 215, 249, 271, 281 and 299) is studied by first-principles isokinetic molecular dynamics simulations. The irregular size dependence of the melting temperatures Tm observed in the calorimetry experiments of Schmidt et al. [Nature (London) 393, 238 (1998)] is quantitatively reproduced. We demonstrate that structural effects alone can explain all broad features of experimental observations. Specifically, maxima in Tm(N) correlate with high surface stability and with structural features such as a high compactness degree.

Liquid-vapor interface in liquid binary alloys: an ab initio molecular dynamics study

We report the results of an ab initio molecular dynamics simulation of the liquid-vapor interface of two binary liquid alloys, Na(0.3)K(0.7) and Li(0.4)Na(0.6), whose bulk behavior exhibits rather differing ordering tendencies. The study has been performed using samples of 2000 and 3000 particles, respectively, in a slab geometry with periodic boundary conditions. In both cases, the total ionic density distributions along the normal to the interface display some layering with a virtually pure monolayer of the lower surface tension component located outermost at the interface. However, the two systems behave very differently below the interface which can be accounted for by their different ordering tendencies in the bulk.