Statistical mechanics model for the transit free energy of monatomic liquids

Development of a Multiphase Beryllium Equation of State and Physics-based Variations

We construct a family of beryllium (Be) multiphase equation of state (EOS) models that consists of a baseline ("optimal") EOS and variations on the baseline to account for physics-based uncertainties. The Be baseline EOS is constructed to reproduce a set of self-consistent data and theory including known phase boundaries, the principal Hugoniot, isobars, and isotherms from diamond-anvil cell experiments. Three phases are considered, including the known hexagonal closed-packed (hcp) phase, the liquid, and the theoretically predicted high-pressure body-centered cubic (bcc) phase. Since both the high-temperature liquid and high-pressure bcc phases lack any experimental data, we carry out ab initio density functional theory (DFT) calculations to obtain new information about the EOS properties for these two regions. At extremely high temperature conditions (>87 eV), DFT-based quantum molecular dynamics simulations are performed for multiple liquid densities using the state-of-the-art Spectral Quadrature methodology in order to validate our selected models for the ion- and electron-thermal free energies of the liquid. We have also performed DFT simulations of hcp and bcc with different exchange-correlation functionals to examine their impact on bcc compressibility, which bound the hcp-bcc transition pressure to within 4 ± 0.5 Mbar. Our baseline EOS predicts the first density maximum along the Hugoniot to be 4.4-fold in compression, while the hcp-bcc-liquid triple-point pressure is predicted to be at 2.25 Mbar. In addition to the baseline EOS, we have generated eight variations to accommodate multiple sources of potential uncertainties such as (1) the choice of free-energy models, (2) differences in theoretical treatments, (3) experimental uncertainties, and (4) lack of information. These variations are designed to provide a reasonable representation of nonstatistical uncertainties for the Be EOS and may be used to assess its sensitivity to different inertial-confinement fusion capsule designs.

V-T theory for the self-intermediate scattering function in a monatomic liquid

In V-T theory the atomic motion is harmonic vibrations in a liquid-specific potential energy valley, plus transits, which move the system rapidly among the multitude of such valleys. In its first application to the self intermediate scattering function (SISF), V-T theory produced an accurate account of molecular dynamics (MD) data at all wave numbers q and time t. Recently, analysis of the mean square displacement (MSD) resolved a crossover behavior that was not observed in the SISF study. Our purpose here is to apply the more accurate MSD calibration to the SISF, and assess the results. We derive and discuss the theoretical equations for vibrational and transit contributions to the SISF. The time evolution is divided into three successive intervals: the vibrational interval when the vibrational contribution alone accurately accounts for the MD data; the crossover when the vibrational contribution saturates and the transit contribution becomes resolved; and the diffusive interval when the transit contribution alone accurately accounts for the MD data. The resulting theoretical error is extremely small at all q and t. V-T theory is compared to mode-coupling theories for the MSD and SISF, and to recent developments in Brownian motion experiments and theory.

First-principles calculation of entropy for liquid metals

We demonstrate the accurate calculation of entropies and free energies for a variety of liquid metals using an extension of the two-phase thermodynamic (2PT) model based on a decomposition of the velocity autocorrelation function into gas-like (hard sphere) and solid-like (harmonic) subsystems. The hard sphere model for the gas-like component is shown to give systematically high entropies for liquid metals as a direct result of the unphysical Lorentzian high-frequency tail. Using a memory function framework we derive a generally applicable velocity autocorrelation and frequency spectrum for the diffusive component which recovers the low-frequency (long-time) behavior of the hard sphere model while providing for realistic short-time coherence and high-frequency tails to the spectrum. This approach provides a significant increase in the accuracy of the calculated entropies for liquid metals and is compared to ambient pressure data for liquid sodium, aluminum, gallium, tin, and iron. The use of this method for the determination of melt boundaries is demonstrated with a calculation of the high-pressure bcc melt boundary for sodium. With the significantly improved accuracy available with the memory function treatment for softer interatomic potentials, the 2PT model for entropy calculations should find broader application in high energy density science, warm dense matter, planetary science, geophysics, and material science.