First principles molecular dynamics of dense plasmas

Spectral-partitioned Kohn-Sham density functional theory

We introduce a general, variational scheme for systematic approximation of a given Kohn-Sham free-energy functional by partitioning the density matrix into distinct spectral domains, each of which may be spanned by an independent diagonal representation without requirement of mutual orthogonality. It is shown that by generalizing the entropic contribution to the free energy to allow for independent representations in each spectral domain, the free energy becomes an upper bound to the exact (unpartitioned) Kohn-Sham free energy, attaining this limit as the representations approach Kohn-Sham eigenfunctions. A numerical procedure is devised for calculation of the generalized entropy associated with spectral partitioning of the density matrix. The result is a powerful framework for Kohn-Sham calculations of systems whose occupied subspaces span multiple energy regimes. As a case in point, we apply the proposed framework to warm- and hot-dense matter described by finite-temperature density functional theory, where at high energies the density matrix is represented by that of the free-electron gas, while at low energies it is variationally optimized. We derive expressions for the spectral-partitioned Kohn-Sham Hamiltonian, atomic forces, and macroscopic stresses within the projector-augmented wave (PAW) and the norm-conserving pseudopotential methods. It is demonstrated that at high temperatures, spectral partitioning facilitates accurate calculations at dramatically reduced computational cost. Moreover, as temperature is increased, fewer exact Kohn-Sham states are required for a given accuracy, leading to further reductions in computational cost. Finally, it is shown that standard multiprojector expansions of electronic orbitals within atomic spheres in the PAW method lack sufficient completeness at high temperatures. Spectral partitioning provides a systematic solution for this fundamental problem.

First-principles equation of state database for warm dense matter computation

We put together a first-principles equation of state (FPEOS) database for matter at extreme conditions by combining results from path integral Monte Carlo and density functional molecular dynamics simulations of the elements H, He, B, C, N, O, Ne, Na, Mg, Al, and Si as well as the compounds LiF, B_{4}C, BN, CH_{4}, CH_{2}, C_{2}H_{3}, CH, C_{2}H, MgO, and MgSiO_{3}. For all these materials, we provide the pressure and internal energy over a density-temperature range from ∼0.5 to 50 g cm^{-3} and from ∼10^{4} to 10^{9} K, which are based on ∼5000 different first-principles simulations. We compute isobars, adiabats, and shock Hugoniot curves in the regime of L- and K-shell ionization. Invoking the linear mixing approximation, we study the properties of mixtures at high density and temperature. We derive the Hugoniot curves for water and alumina as well as for carbon-oxygen, helium-neon, and CH-silicon mixtures. We predict the maximal shock compression ratios of H_{2}O, H_{2}O_{2}, Al_{2}O_{3}, CO, and CO_{2} to be 4.61, 4.64, 4.64, 4.89, and 4.83, respectively. Finally we use the FPEOS database to determine the points of maximum shock compression for all available binary mixtures. We identify mixtures that reach higher shock compression ratios than their end members. We discuss trends common to all mixtures in pressure-temperature and particle-shock velocity spaces. In the Supplemental Material, we provide all FPEOS tables as well as computer codes for interpolation, Hugoniot calculations, and plots of various thermodynamic functions.

Structure and dynamics of warm dense aluminum: a molecular dynamics study with density functional theory and deep potential

We perform a systematic study on the structure and dynamics of warm dense aluminum (Al) at temperatures ranging from 0.5 to 5.0 eV with molecular dynamics utilizing both density functional theory (DFT) and the deep potential (DP) method. On one hand, unlike the Thomas-Fermi kinetic energy density functional (KEDF), we find that the orbital-free DFT method with the Wang-Teter non-local KEDF yields properties of warm dense Al that agree well with the Kohn-Sham DFT method, enabling accurate orbital-free DFT simulations of warm dense Al at relatively low temperatures. On the other hand, the DP method constructs a deep neural network that has a high accuracy in reproducing short- and long-ranged properties of warm dense Al when compared to the DFT methods. The DP method is orders of magnitudes faster than DFT and is well-suited for simulating large systems and long trajectories to yield accurate properties of warm dense Al. Our results suggest that the combination of DFT methods and the DP model is a powerful tool for accurately and efficiently simulating warm dense matter.

Path integral Monte Carlo simulations of warm dense aluminum

We perform first-principles path integral Monte Carlo (PIMC) and density functional theory molecular dynamics (DFT-MD) calculations to explore warm dense matter states of aluminum. Our equation of state (EOS) simulations cover a wide density-temperature range of 0.1-32.4gcm^{-3} and 10^{4}-10^{8} K. Since PIMC and DFT-MD accurately treat effects of the atomic shell structure, we find two compression maxima along the principal Hugoniot curve attributed to K-shell and L-shell ionization. The results provide a benchmark for widely used EOS tables, such as SESAME, QEOS, and models based on Thomas-Fermi and average-atom techniques. A subsequent multishock analysis provides a quantitative assessment for how much heating occurs relative to an isentrope in multishock experiments. Finally, we compute heat capacity, pair-correlation functions, the electronic density of states, and 〈Z〉 to reveal the evolution of the plasma structure and ionization behavior.

Dynamic electron-ion collisions and nuclear quantum effects in quantum simulation of warm dense matter

The structural, thermodynamic and transport properties of warm dense matter (WDM) are crucial to the fields of astrophysics and planet science, as well as inertial confinement fusion. WDM refers to the states of matter in a regime of temperature and density between cold condensed matter and hot ideal plasmas, where the density is from near-solid up to ten times solid density, and the temperature between 0.1 and 100 eV. In the WDM regime, matter exhibits moderately or strongly coupled, partially degenerate properties. Therefore, the methods used to deal with condensed matter and isolated atoms need to be properly validated for WDM. It is therefore a big challenge to understand WDM within a unified theoretical description with reliable accuracy. Here, we review the progress in the theoretical study of WDM with state-of-the-art simulations, i.e. quantum Langevin molecular dynamics and first principles path integral molecular dynamics. The related applications for WDM are also included.

Importance of finite-temperature exchange correlation for warm dense matter calculations

The effects of an explicit temperature dependence in the exchange correlation (XC) free-energy functional upon calculated properties of matter in the warm dense regime are investigated. The comparison is between the Karasiev-Sjostrom-Dufty-Trickey (KSDT) finite-temperature local-density approximation (TLDA) XC functional [Karasiev et al., Phys. Rev. Lett. 112, 076403 (2014)PRLTAO0031-900710.1103/PhysRevLett.112.076403] parametrized from restricted path-integral Monte Carlo data on the homogeneous electron gas (HEG) and the conventional Monte Carlo parametrization ground-state LDA XC [Perdew-Zunger (PZ)] functional evaluated with T-dependent densities. Both Kohn-Sham (KS) and orbital-free density-functional theories are used, depending upon computational resource demands. Compared to the PZ functional, the KSDT functional generally lowers the dc electrical conductivity of low-density Al, yielding improved agreement with experiment. The greatest lowering is about 15% for T=15 kK. Correspondingly, the KS band structure of low-density fcc Al from the KSDT functional exhibits a clear increase in interband separation above the Fermi level compared to the PZ bands. In some density-temperature regimes, the deuterium equations of state obtained from the two XC functionals exhibit pressure differences as large as 4% and a 6% range of differences. However, the hydrogen principal Hugoniot is insensitive to the explicit XC T dependence because of cancellation between the energy and pressure-volume work difference terms in the Rankine-Hugoniot equation. Finally, the temperature at which the HEG becomes unstable is T≥7200 K for the T-dependent XC, a result that the ground-state XC underestimates by about 1000 K.

Kinetic theory molecular dynamics and hot dense matter: theoretical foundations

Electrons are weakly coupled in hot, dense matter that is created in high-energy-density experiments. They are also mildly quantum mechanical and the ions associated with them are classical and may be strongly coupled. In addition, the dynamical evolution of plasmas under these hot, dense matter conditions involve a variety of transport and energy exchange processes. Quantum kinetic theory is an ideal tool for treating the electrons but it is not adequate for treating the ions. Molecular dynamics is perfectly suited to describe the classical, strongly coupled ions but not the electrons. We develop a method that combines a Wigner kinetic treatment of the electrons with classical molecular dynamics for the ions. We refer to this hybrid method as "kinetic theory molecular dynamics," or KTMD. The purpose of this paper is to derive KTMD from first principles and place it on a firm theoretical foundation. The framework that KTMD provides for simulating plasmas in the hot, dense regime is particularly useful since current computational methods are generally limited by their inability to treat the dynamical quantum evolution of the electronic component. Using the N-body von Neumann equation for the electron-proton plasma, three variations of KTMD are obtained. Each variant is determined by the physical state of the plasma (e.g., collisional versus collisionless). The first variant of KTMD yields a closed set of equations consisting of a mean-field quantum kinetic equation for the electron one-particle distribution function coupled to a classical Liouville equation for the protons. The latter equation includes both proton-proton Coulombic interactions and an effective electron-proton interaction that involves the convolution of the electron density with the electron-proton Coulomb potential. The mean-field approach is then extended to incorporate equilibrium electron-proton correlations through the Singwi-Tosi-Land-Sjolander (STLS) ansatz. This is the second variant of KTMD. The STLS contribution produces an effective electron-proton interaction that involves the electron-proton structure factor, thereby extending the usual mean-field theory to correlated but near equilibrium systems. Finally, a third variant of KTMD is derived. It includes dynamical electrons and their correlations coupled to a MD description for the ions. A set of coupled equations for the one-particle electron Wigner function and the electron-electron and electron-proton correlation functions are coupled to a classical Liouville equation for the protons. This latter variation has both time and momentum dependent correlations.

Accurate homogeneous electron gas exchange-correlation free energy for local spin-density calculations

An accurate analytical parametrization for the exchange-correlation free energy of the homogeneous electron gas, including interpolation for partial spin polarization, is derived via thermodynamic analysis of recent restricted path integral Monte Carlo (RPIMC) data. This parametrization constitutes the local spin density approximation (LSDA) for the exchange-correlation functional in density functional theory. The new finite-temperature LSDA reproduces the RPIMC data well, satisfies the correct high-density and low- and high-T asymptotic limits, and is well behaved beyond the range of the RPIMC data, suggestive of broad utility.

Comparison of density functional approximations and the finite-temperature Hartree-Fock approximation in warm dense lithium

We compare the behavior of the finite-temperature Hartree-Fock model with that of thermal density functional theory using both ground-state and temperature-dependent approximate exchange functionals. The test system is bcc Li in the temperature-density regime of warm dense matter (WDM). In this exchange-only case, there are significant qualitative differences in results from the three approaches. Those differences may be important for Born-Oppenheimer molecular dynamics studies of WDM with ground-state approximate density functionals and thermal occupancies. Such calculations require reliable regularized potentials over a demanding range of temperatures and densities. By comparison of pseudopotential and all-electron results at T=0 K for small Li clusters of local bcc symmetry and bond lengths equivalent to high density bulk Li, we determine the density ranges for which standard projector augmented wave (PAW) and norm-conserving pseudopotentials are reliable. Then, we construct and use all-electron PAW data sets with a small cutoff radius that are valid for lithium densities up to at least 80 g/cm{3}.

Dynamic ionic clusters with flowing electron bubbles from warm to hot dense iron along the Hugoniot curve

Complex structures of warm and hot dense matter are essential to understanding the behavior of materials in high energy density processes and provide new features of matter constitutions. Here, around a new unified first-principles determined Hugoniot curve of iron from the normal condensed condition up to 1 Gbar, the novel structures characterized by the ionic clusters with electron bubbles are found using quantum Langevin molecular dynamics. Subsistence of complex clusters can persist in the time scale of 50 fs dynamically with quantum flowing bubbles, which are produced by the interplay of Fermi electron degeneracy, the ionic coupling, and the dynamical nature. With the inclusion of those complicated features in quantum Langevin molecular dynamics, the present equation of states could serve as a first-principles based database in a wide range of temperatures and densities.

All-electron path integral Monte Carlo simulations of warm dense matter: application to water and carbon plasmas

We develop an all-electron path integral Monte Carlo method with free-particle nodes for warm dense matter and apply it to water and carbon plasmas. We thereby extend path integral Monte Carlo studies beyond hydrogen and helium to elements with core electrons. Path integral Monte Carlo results for pressures, internal energies, and pair-correlation functions compare well with density functional theory molecular dynamics calculations at temperatures of (2.5-7.5)×10(5) K, and both methods together form a coherent equation of state over a density-temperature range of 3-12 g/cm(3) and 10(4)-10(9) K.

Equation of state for partially ionized carbon and oxygen mixtures at high temperatures

The equation of state (EOS) for partially ionized carbon, oxygen, and carbon-oxygen mixtures at temperatures 3×10(5)K is less than or approximately equal to T is less than or approximately equal to 3×10(6) K is calculated over a wide range of densities, using the method of free energy minimization in the framework of the chemical picture of plasmas. The free energy model is an improved extension of our model previously developed for pure carbon [Potekhin, Massacrier, and Chabrier, Phys. Rev. E 72, 046402 (2005)]. The internal partition functions of bound species are calculated by a self-consistent treatment of each ionization stage in the plasma environment taking into account pressure ionization. The long-range Coulomb interactions between ions and screening of the ions by free electrons are included using our previously published analytical model, recently improved, in particular for the case of mixtures. We also propose a simple but accurate method of calculation of the EOS of partially ionized binary mixtures based on detailed ionization balance calculations for pure substances.

Unified first principles description from warm dense matter to ideal ionized gas plasma: electron-ion collisions induced friction

Electron-ion interactions are central to numerous phenomena in the warm dense matter (WDM) regime and at higher temperature. The electron-ion collisions induced friction at high temperature is introduced in the procedure of ab initio molecular dynamics using the Langevin equation based on density functional theory. In this framework, as a test for Fe and H up to 1000 eV, the equation of state and the transition of electronic structures of the materials with very wide density and temperature can be described, which covers a full range of WDM up to high energy density physics. A unified first principles description from condensed matter to ideal ionized gas plasma is constructed.

Static and dynamic conductivity of warm dense matter within a density-functional approach: application to aluminum and gold

The conductivity sigma(omega) of dense Al and Au plasmas is considered where all the needed inputs are obtained from density-functional theory (DFT). These calculations involve a self-consistent determination of (i) the equation of state and the ionization balance, (ii) evaluation of the ion-ion and ion-electron pair-distribution functions, (iii) determination of the scattering amplitudes, and finally the conductivity. We present results for Al and Au for compressions 0.1-2.0, and in the temperature range T=0.1-10 eV. Excellent agreement with recent first-principles calculations using multi-ion density-functional molecular dynamics is obtained where the data fields overlap. We review first-principles approaches to the optical conductivity, including many-body perturbation theory, molecular-dynamics evaluations, and simplified time-dependent DFT. The modification to the Drude conductivity in the presence of shallow bound states in typical Al plasmas is examined and numerical results are given at the level of the Fermi Golden Rule and an approximate time-dependent DFT.

Very-high-temperature molecular dynamics

It is shown that a modified scheme of density functional theory, using the Thomas-Fermi kinetic energy functional for the electrons, is well suited to perform very-high-temperature molecular dynamics simulations on high-Z elements. As an example, iron on the principal Hugoniot is simulated up to 5 keV and 5 times the normal density, giving an equation of state in agreement with current models. Ionic structure is obtained and is given to an excellent level of precision by the structure of the one-component plasma computed for a coupling parameter corresponding to Thomas-Fermi ionization.

Equation of state for partially ionized carbon at high temperatures

Equation of state for partially ionized carbon at temperatures T approximately > or = 10(5) K is calculated in a wide range of densities, using the method of free energy minimization in the framework of the chemical picture of plasmas. The free energy model includes the internal partition functions of bound species. The latter are calculated by a self-consistent treatment of each ionization stage in the plasma environment taking into account pressure ionization. The long-range Coulomb interactions between ions and screening of the ions by free electrons are included using our previously published analytical model.

Electrical conductivity of hot expanded aluminum: experimental measurements and ab initio calculations

Experimental measurements and theoretical calculations of the electrical conductivity of aluminum are presented in the strongly coupled partially degenerate regime (rho=0.3 g/cm(3), 5000<T<15 000 K). The experiments were performed in an isochoric plasma closed vessel designed to confine electrical plasma discharges up to 1.5 GPa. Aluminum properties were determined theoretically by ab initio molecular dynamics simulations in the local density approximation, from which the conductivity was computed using the Kubo-Greenwood formula. The theoretical results were validated in the dense coupled regime against previously published experimental results and then applied to our experimental low density regime, showing that the theoretical results overestimate the experimental conductivities.

Combined pressure and electrical-resistivity measurements of warm dense aluminum and titanium plasmas

Electrical resistivity, pressure, and internal energy variation of warm dense correlated titanium (density 0.2 g/cm(3)) and aluminum (density 0.1 g/cm(3)) plasmas are measured using a homogeneous and thermally equilibrated media produced inside an isochoric closed-vessel plasma. These data are compared to detailed calculations based on the density functional theory. In the studied temperature range (15,000-30,000 K), it appears that both exchange-correlation and ion-ion interaction treatments are of great importance to calculate accurate theoretical values.