Equation of state of a dense plasma by orbital-free and quantum molecular dynamics: examination of two isothermal-isobaric mixing rules

Mixture of hydrogen and methane under planetary interior conditions

We employ first-principles molecular dynamics simulations to provide equation-of-state data, pair distribution functions (PDFs), diffusion coefficients, and band gaps of a mixture of hydrogen and methane under planetary interior conditions as relevant for Uranus, Neptune, and similar icy exoplanets. We test the linear mixing approximation, which is fulfilled within a few percent for the chosen P-T conditions. Evaluation of the PDFs reveals that methane molecules dissociate into carbon clusters and free hydrogen atoms at temperatures greater than 3000 K. At high temperatures, the clusters are found to be short-lived. Furthermore, we calculate the electrical conductivity from which we derive the non-metal-to-metal transition region of the mixture. We also calculate the electrical conductivity along the P-T profile of Uranus [N. Nettelmann et al., Planet. Space Sci., 2013, 77, 143-151] and observe the transition of the mixture from a molecular to an atomic fluid as a function of the radius of the planet. The density and temperature ranges chosen in our study can be achieved using dynamic shock compression experiments and seek to aid such future experiments. Our work also provides a relevant data set for a better understanding of the interior, evolution, luminosity, and magnetic field of the ice giants in our solar system and beyond.

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.

Effect of ionic disorder on the principal shock Hugoniot

The effect of ionic disorder on the principal Hugoniot is investigated using multiple scattering theory to very high pressure (Gbar). Calculations using molecular dynamics to simulate ionic disorder are compared to those with a fixed crystal lattice, for both carbon and aluminum. For the range of conditions considered here we find that ionic disorder has a relatively minor influence. It is most important at the onset of shell ionization and we find that, at higher pressures, the subtle effect of the ionic environment is overwhelmed by the larger number of ionized electrons with higher thermal energies.

Equation of state modeling with pseudoatom molecular dynamics

Using a modified version of the pseudoatom molecular-dynamics approach, the silicon and oxygen equations of state were generated and then employed to construct the equation of state of silicon dioxide. The results are supported by the close agreement with ab initio simulations of the silicon pressure and experimental shock Hugoniot of silicon dioxide. Ion thermal contributions to thermodynamic functions provided by the PAMD simulations are compared to their counterparts obtained with the one-component plasma and charged-hard-sphere approximations.

Multiple scattering theory for dense plasmas

Dense plasmas occur in stars, giant planets, and in inertial fusion experiments. Accurate modeling of the electronic structure of these plasmas allows for prediction of material properties that can in turn be used to simulate these astrophysical objects and terrestrial experiments. But modeling them remains a challenge. Here we explore the Korringa-Kohn-Rostoker Green's function (KKR-GF) method for this purpose. We find that it is able to predict equation of state in good agreement with other state-of-the-art methods, where they are accurate and viable. In addition, it is shown that the computational cost does not significantly change with temperature, in contrast with other approaches. Moreover, the method does not use pseudopotentials-core states are calculated self consistently. We conclude that KKR-GF is a very promising method for dense plasma simulation.

Correlation and transport properties for mixtures at constant pressure and temperature

Transport properties of mixtures of elements in the dense plasma regime play an important role in natural astrophysical and experimental systems, e.g., inertial confinement fusion. We present a series of orbital-free molecular dynamics simulations on dense plasma mixtures with comparison to a global pseudo ion in jellium model. Hydrogen is mixed with elements of increasingly high atomic number (lithium, carbon, aluminum, copper, and silver) at a fixed temperature of 100 eV and constant pressure set by pure hydrogen at 2g/cm^{3}, namely, 370 Mbars. We compute ionic transport coefficients, such as self-diffusion, mutual diffusion, and viscosity for various concentrations. Small concentrations of the heavy atoms significantly change the density of the plasma and decrease the transport coefficients. The structure of the mixture evidences a strong Coulomb coupling between heavy ions and the appearance of a broad correlation peak at short distances between hydrogen atoms. The concept of an effective one component plasma is used to quantify the overcorrelation of the light element induced by the admixture of a heavy element.

Dissociation along the principal Hugoniot of the Laser Mégajoule ablator material

Glow discharge polymer hydrocarbon (GDP-CH) is used as the ablator material in inertial confinement fusion (ICF) capsules for the Laser Mégajoule and National Ignition Facility. Due to its fabrication process, GDP-CH chemical composition and structure differ from commercially available plastics and detailed knowledge of its properties in the warm dense matter regime is needed to achieve accurate design of ICF capsules. First-principles ab initio simulations of the GDP-CH principal Hugoniot up to 8 Mbar were performed using the quantum molecular dynamics (QMD) code abinit and showed that atomic bond dissociation has an effect on the compressibility. Results from these simulations are used to parametrize a quantum semiempirical model in order to generate a tabulated equation of state that includes dissociation. Hugoniot measurements obtained from an experiment conducted at the LULI2000 laser facility confirm QMD simulations as well as EOS modeling. We conclude by showing the EOS model influence on shock timing in a hydrodynamic simulation.

Transport properties of an asymmetric mixture in the dense plasma regime

We study how concentration changes ionic transport properties along isobars-isotherms for a mixture of hydrogen and silver, representative of turbulent layers relevant to inertial confinement fusion and astrophysics. Hydrogen will typically be fully ionized while silver will be only partially ionized but can have a large effective charge. This will lead to very different physical conditions for the H and Ag. Large first principles orbital free molecular dynamics simulations are performed and the resulting transport properties are analyzed. Comparisons are made with transport theory in the kinetic regime and in the coupled regime. The addition of a small amount of heavy element in a light material has a dramatic effect on viscosity and diffusion of the mixture. This effect is explained through kinetic theory as a manifestation of a crossover between classical diffusion and Lorentz diffusion.

Experimental and ab initio investigations of microscopic properties of laser-shocked Ge-doped ablator

Plastic materials (CH) doped with mid-Z elements are used as ablators in inertial confinement fusion (ICF) capsules and in their surrogates. Hugoniot equation of state (EOS) and electronic properties of CH doped with germanium (at 2.5% and 13% dopant fractions) are investigated experimentally up to 7 Mbar using velocity and reflectivity measurements of shock fronts on the GEKKO laser at Osaka University. Reflectivity and temperature measurements were updated using a quartz standard. Shocked quartz reflectivity was measured at 532 and 1064 nm. Theoretical investigation of shock pressure and reflectivity was then carried out by ab initio simulations using the quantum molecular dynamics (QMD) code abinit and compared with tabulated average atom EOS models. We find that shock states calculated by QMD are in better agreement with experimental data than EOS models because of a more accurate description of ionic structure. We finally discuss electronic properties by comparing reflectivity data to a semiconductor gap closure model and to QMD simulations.

Impact of oxygen on the 300-K isotherm of Laser Megajoule ablator using ab initio simulation

The ablator material for inertial confinement fusion (ICF) capsules on the Laser Mégajoule is a glow-discharge polymer (GDP) plastic. Its equation of state (EOS) is of primary importance for the design of such capsules, since it has direct consequences on shock timing and is essential to mitigate hydrodynamic instabilities. Using ab initio molecular dynamics (AIMD), we have investigated the 300-K isotherm of amorphous CH(1.37)O(0.08) plastic, whose structure is close to GDP plastic. The 300-K isotherm, which is often used as a cold curve within tabular EOS, is an important contribution of the EOS in the multimegabar pressure range. AIMD results are compared to analytic models within tabular EOS, pointing out large discrepancies. In addition, we show that the effect of oxygen decreases 300-K isotherm pressure by 10%-15%. The implication of these observations is the ability to improve ICF target performance, which is essential to achieve fusion ignition.

Ion-ion dynamic structure factor of warm dense mixtures

The ion-ion dynamic structure factor of warm dense matter is determined using the recently developed pseudoatom molecular dynamics method [Starrett et al., Phys. Rev. E 91, 013104 (2015)]. The method uses density functional theory to determine ion-ion pair interaction potentials that have no free parameters. These potentials are used in classical molecular dynamics simulations. This constitutes a computationally efficient and realistic model of dense plasmas. Comparison with recently published simulations of the ion-ion dynamic structure factor and sound speed of warm dense aluminum finds good to reasonable agreement. Using this method, we make predictions of the ion-ion dynamical structure factor and sound speed of a warm dense mixture-equimolar carbon-hydrogen. This material is commonly used as an ablator in inertial confinement fusion capsules, and our results are amenable to direct experimental measurement.