Path integral monte carlo calculation of the deuterium hugoniot

Effect of zero-point motion on properties of quantum particles adsorbed on a substrate

We qualitatively investigate the effect of zero-point motion (ZPM) on the structure and properties of a film composed of quantum particles adsorbed on a graphite substrate. The amplitude of ZPM is controlled by a change of the particle mass while keeping the interactions fixed. In that sense it is assumed that the interactions can be controlled by future doping methods. The worm-algorithm path integral Monte Carlo (WAPIMC) method is applied to simulate this system in the grand-canonical ensemble, where particles can be exchanged with the external particle reservoir. Another method, namely the multiconfigurational time-dependent Hartree method for bosons is additionally applied to verify some of the WAPIMC results and to provide further information on the entropy and the condensate fraction. Several important findings are reported. It is found that ZPM plays an important role in defining order and disorder in the crystalline structure of the adsorbed film. The total energy of the film drops with a reduction in the amplitude of ZPM, that is, it becomes more negative which is an indication to stronger adsorption. For a few particle numbers, a significant condensate fraction is detected that however drops sharply at critical values of the ZPM amplitude. Most importantly, a connection is established between chaos, in coordinate as well as momentum space, and the Heisenberg uncertainty principle. The importance of the present study lies in the fact that adsorbed two-dimensional films serve as an excellent experimental testbed for demonstrating low-dimensional quantum phenomena in the ground state. The present examination contributes also to a further understanding of the properties of heavy quantum particles adsorbed on substrates.

Equation of state of partially ionized hydrogen and deuterium plasma revisited

We present improved first-principle fermionic path integral Monte Carlo (PIMC) simulation results for a dense partially ionized hydrogen (deuterium) plasma, for temperatures in the range 15000K≤T≤400000K and densities 7×10^{-7}g/cm^{3}≤ρ_{H}≤0.085g/cm^{3} (1.4×10^{-6}g/cm^{3}≤ρ_{D}≤0.17g/cm^{3}), corresponding to 100≥r_{s}≥2, where r_{s}=r[over ¯]/a_{B} is the ratio of the mean interparticle distance to the Bohr radius. These simulations are based on the fermionic propagator PIMC (FP-PIMC) approach in the grand canonical ensemble [Filinov et al., Contrib. Plasma Phys. 61, e202100112 (2021)0863-104210.1002/ctpp.202100112] and fully account for correlation and quantum degeneracy and spin effects. For the application to hydrogen and deuterium, we develop a combination of the fourth-order factorization and the pair product ansatz for the density matrix. Moreover, we avoid the fixed node approximation that may lead to uncontrolled errors in restricted PIMC (RPIMC). Our results allow us to critically reevaluate the accuracy of the RPIMC simulations for hydrogen by Hu et al. [Phys. Rev. B 84, 224109 (2011)1098-012110.1103/PhysRevB.84.224109] and of various chemical models. The deviations are generally found to be small, but for the lowest temperature, T=15640 K they reach several percent. We present detailed tables with our first principles results for the pressure and energy isotherms. We expect our updated results will serve as a valuable benchmark for comparison with other methods.

Development of a new quantum trajectory molecular dynamics framework

An extension to the wave packet description of quantum plasmas is presented, where the wave packet can be elongated in arbitrary directions. A generalized Ewald summation is constructed for the wave packet models accounting for long-range Coulomb interactions and fermionic effects are approximated by purpose-built Pauli potentials, self-consistent with the wave packets used. We demonstrate its numerical implementation with good parallel support and close to linear scaling in particle number, used for comparisons with the more common wave packet employing isotropic states. Ground state and thermal properties are compared between the models with differences occurring primarily in the electronic subsystem. Especially, the electrical conductivity of dense hydrogen is investigated where a 15% increase in DC conductivity can be seen in our wave packet model compared with other models. This article is part of the theme issue 'Dynamic and transient processes in warm dense matter'.

Equilibrium properties of warm dense deuterium calculated by the wave packet molecular dynamics and density functional theory method

A joint simulation method based on the wave packet molecular dynamics and density functional theory (WPMD-DFT) is applied to study warm dense deuterium (nonideal deuterium plasmas). This method was developed recently as an extension of the wave packet molecular dynamics (WPMD) in which the equations of motion are solved simultaneously for classical ions and semiclassical electrons represented as Gaussian wave packets. Compared to the classical molecular dynamics and WPMD simulations, the method of WPMD-DFT provides a more accurate representation of quantum effects such as electron-ion coupling and electron degeneracy. It allows studying nonadiabatic dynamics of electrons and ions in equilibrium and nonequilibrium states while being more accurate and efficient at high densities than WPMD and classical molecular dynamics. In the paper, we discuss particular features of the method such as special boundary conditions and the procedure of isentrope calculation as well as the results obtained by WPMD-DFT for the shock-compressed deuterium. The compression isentrope and principal Hugoniot curves obtained by WPMD-DFT are compared with available experimental data and other simulation approaches to validate the method. It opens up a possibility of further application of the method to study nonequilibrium states and relaxation processes.

Momentum distribution function and short-range correlations of the warm dense electron gas: Ab initio quantum Monte Carlo results

In a classical plasma the momentum distribution, n(k), decays exponentially, for large k, and the same is observed for an ideal Fermi gas. However, when quantum and correlation effects are relevant simultaneously, an algebraic decay, n_{∞}(k)∼k^{-8} has been predicted. This is of relevance for cross sections and threshold processes in dense plasmas that depend on the number of energetic particles. Here we present extensive ab initio results for the momentum distribution of the nonideal uniform electron gas at warm dense matter conditions. Our results are based on first principle fermionic path integral Monte Carlo (CPIMC) simulations and clearly confirm the k^{-8} asymptotic. This asymptotic behavior is directly linked to short-range correlations which are analyzed via the on-top pair distribution function (on-top PDF), i.e., the PDF of electrons with opposite spin. We present extensive results for the density and temperature dependence of the on-top PDF and for the momentum distribution in the entire momentum range.

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.

Screened activity expansion for the grand potential of a quantum plasma and how to derive approximate equations of state compatible with electroneutrality

We consider a quantum multicomponent plasma made with S species of point charged particles interacting via the Coulomb potential. We derive the screened activity series for the pressure in the grand-canonical ensemble within the Feynman-Kac path integral representation of the system in terms of a classical gas of loops. This series is useful for computing equations of state for it is nonperturbative with respect to the strength of the interaction and it involves relatively few diagrams at a given order. The known screened activity series for the particle densities can be recovered by differentiation. The particle densities satisfy local charge neutrality because of a Debye-dressing mechanism of the diagrams in these series. We introduce a new general neutralization prescription, based on this mechanism, for deriving approximate equations of state where consistency with electroneutrality is automatically ensured. This prescription is compared to other ones, including a neutralization scheme inspired by the Lieb-Lebowitz theorem and based on the introduction of (S-1) suitable independent combinations of the activities. Eventually, we briefly argue how the activity series for the pressure, combined with the Debye-dressing prescription, can be used for deriving approximate equations of state at moderate densities, which include the contributions of recombined entities made with three or more particles.

Shock Compression of Liquid Deuterium up to 1 TPa

We present laser-driven shock compression experiments on cryogenic liquid deuterium to 550 GPa along the principal Hugoniot and reflected-shock data up to 1 TPa. High-precision interferometric Doppler velocimetry and impedance-matching analysis were used to determine the compression accurately enough to reveal a significant difference as compared to state-of-the-art ab initio calculations and thus, no single equation of state model fully matches the principal Hugoniot of deuterium over the observed pressure range. In the molecular-to-atomic transition pressure range, models based on density functional theory calculations predict the maximum compression accurately. However, beyond 250 GPa along the principal Hugoniot, first-principles models exhibit a stiffer response than the experimental data. Similarly, above 500 GPa the reflected shock data show 5%-7% higher compression than predicted by all current models.

Extremely Low Electron-ion Temperature Relaxation Rates in Warm Dense Hydrogen: Interplay between Quantum Electrons and Coupled Ions

Theoretical and computational modeling of nonequilibrium processes in warm dense matter represents a significant challenge. The electron-ion relaxation process in warm dense hydrogen is investigated here by nonequilibrium molecular dynamics using the constrained electron force field (CEFF) method. CEFF evolves wave packets that incorporate dynamic quantum diffraction that obviates the Coulomb catastrophe. Predictions from this model reveal temperature relaxation times as much as three times longer than prior molecular dynamics results based on quantum statistical potentials. Through analyses of energy distributions and mean free paths, this result can be traced to delocalization. Finally, an improved GMS [Gericke, Murillo, and Schlanges, Phys. Rev. E 78, 025401 (2008)PRESCM1539-375510.1103/PhysRevE.78.025401] model is proposed, in which the Coulomb logarithms are in good agreement with CEFF results.

Ab Initio Studies on the Stopping Power of Warm Dense Matter with Time-Dependent Orbital-Free Density Functional Theory

Electronic transport properties of warm dense matter, such as electrical or thermal conductivities and nonadiabatic stopping power, are of particular interest to geophysics, planetary science, astrophysics, and inertial confinement fusion (ICF). One example is the α-particle stopping power of dense deuterium-tritium (DT) plasmas, which must be precisely known for current small-margin ICF target designs to ignite. We have developed a time-dependent orbital-free density functional theory (TD-OF-DFT) method for ab initio investigations of the charged-particle stopping power of warm dense matter. Our current dependent TD-OF-DFT calculations have reproduced the recently well-characterized stopping power experiment in warm dense beryllium. For α-particle stopping in warm and solid-density DT plasmas, the ab initio TD-OF-DFT simulations show a lower stopping power up to ∼25% in comparison with three stopping-power models often used in the high-energy-density physics community.

Theoretical and experimental investigation of the equation of state of boron plasmas

We report a theoretical equation of state (EOS) table for boron across a wide range of temperatures (5.1×10^{4}-5.2×10^{8} K) and densities (0.25-49 g/cm^{3}) and experimental shock Hugoniot data at unprecedented high pressures (5608±118 GPa). The calculations are performed with first-principles methods combining path-integral Monte Carlo (PIMC) at high temperatures and density-functional-theory molecular-dynamics (DFT-MD) methods at lower temperatures. PIMC and DFT-MD cross-validate each other by providing coherent EOS (difference <1.5 Hartree/boron in energy and <5% in pressure) at 5.1×10^{5} K. The Hugoniot measurement is conducted at the National Ignition Facility using a planar shock platform. The pressure-density relation found in our shock experiment is on top of the shock Hugoniot profile predicted with our first-principles EOS and a semiempirical EOS table (LEOS 50). We investigate the self-diffusivity and the effect of thermal and pressure-driven ionization on the EOS and shock compression behavior in high-pressure and -temperature conditions. We also study the sensitivity of a polar direct-drive exploding pusher platform to pressure variations based on applying pressure multipliers to LEOS 50 and by utilizing a new EOS model based on our ab initio simulations via one-dimensional radiation-hydrodynamic calculations. The results are valuable for future theoretical and experimental studies and engineering design in high-energy density research.

Quantum molecular dynamics study on the proton exchange, ionic structures, and transport properties of warm dense hydrogen-deuterium mixtures

Comprehensive knowledge of physical properties such as equation of state (EOS), proton exchange, dynamic structures, diffusion coefficients, and viscosities of hydrogen-deuterium mixtures with densities from 0.1 to 5 g/cm^{3} and temperatures from 1 to 50 kK has been presented via quantum molecular dynamics (QMD) simulations. The existing multi-shock experimental EOS provides an important benchmark to evaluate exchange-correlation functionals. The comparison of simulations with experiments indicates that a nonlocal van der Waals density functional (vdW-DF1) produces excellent results. Fraction analysis of molecules using a weighted integral over pair distribution functions was performed. A dissociation diagram together with a boundary where the proton exchange (H_{2}+D_{2}⇌2HD) occurs was generated, which shows evidence that the HD molecules form as the H_{2} and D_{2} molecules are almost 50% dissociated. The mechanism of proton exchange can be interpreted as a process of dissociation followed by recombination. The ionic structures at extreme conditions were analyzed by the effective coordination number model. High-order cluster, circle, and chain structures can be founded in the strongly coupled warm dense regime. The present QMD diffusion coefficient and viscosity can be used to benchmark two analytical one-component plasma (OCP) models: the Coulomb and Yukawa OCP models.

First-principles equation of state and shock compression predictions of warm dense hydrocarbons

We use path integral Monte Carlo and density functional molecular dynamics to construct a coherent set of equations of state (EOS) for a series of hydrocarbon materials with various C:H ratios (2:1, 1:1, 2:3, 1:2, and 1:4) over the range of 0.07-22.4gcm^{-3} and 6.7×10^{3}-1.29×10^{8}K. The shock Hugoniot curve derived for each material displays a single compression maximum corresponding to K-shell ionization. For C:H = 1:1, the compression maximum occurs at 4.7-fold of the initial density and we show radiation effects significantly increase the shock compression ratio above 2 Gbar, surpassing relativistic effects. The single-peaked structure of the Hugoniot curves contrasts with previous work on higher-Z plasmas, which exhibit a two-peak structure corresponding to both K- and L-shell ionization. Analysis of the electronic density of states reveals that the change in Hugoniot structure is due to merging of the L-shell eigenstates in carbon, while they remain distinct for higher-Z elements. Finally, we show that the isobaric-isothermal linear mixing rule for carbon and hydrogen EOS is a reasonable approximation with errors better than 1% for stellar-core conditions.

Molecular dynamics simulation of strong shock waves propagating in dense deuterium, taking into consideration effects of excited electrons

We present a molecular dynamics simulation of shock waves propagating in dense deuterium with the electron force field method [J. T. Su and W. A. Goddard, Phys. Rev. Lett. 99, 185003 (2007)PRLTAO0031-900710.1103/PhysRevLett.99.185003], which explicitly takes the excitation of electrons into consideration. Nonequilibrium features associated with the excitation of electrons are systematically investigated. We show that chemical bonds in D_{2} molecules lead to a more complicated shock wave structure near the shock front, compared with the results of classical molecular dynamics simulation. Charge separation can bring about accumulation of net charges on large scales, instead of the formation of a localized dipole layer, which might cause extra energy for the shock wave to propagate. In addition, the simulations also display that molecular dissociation at the shock front is the major factor that accounts for the "bump" structure in the principal Hugoniot. These results could help to build a more realistic picture of shock wave propagation in fuel materials commonly used in the inertial confinement fusion.

Equation of state of warm dense deuterium and its isotopes from density-functional theory molecular dynamics

Of the two approaches of density-functional theory molecular dynamics, quantum molecular dynamics is limited at high temperature by computational cost whereas orbital-free molecular dynamics, based on an approximation of the kinetic electronic free energy, can be implemented in this domain. In the case of deuterium, it is shown how orbital-free molecular dynamics can be regarded as the limit of quantum molecular dynamics at high temperature for the calculation of the equation of state. To this end, accurate quantum molecular dynamics calculations are performed up to 20 eV at mass densities as low as 0.5g/cm^{3} and up to 10 eV at mass densities as low as 0.2g/cm^{3}. As a result, the limitation in temperature so far attributed to quantum molecular dynamics is overcome and an approach combining quantum and orbital-free molecular dynamics is used to construct an equation of state of deuterium. The thermodynamic domain addressed is that of the fluid phase above 1 eV and 0.2g/cm^{3}. Both pressure and internal energy are calculated as functions of temperature and mass density, and various exchange-correlation contributions are compared. The generalized gradient approximation of the exchange-correlation functional, corrected to approximately include the influence of temperature, is retained and the results obtained are compared to other approaches and to experimental shock data; in parts of the thermodynamic domain addressed, these results significantly differ from those obtained in other first-principles investigations which themselves disagree. The equations of state of hydrogen and tritium above 1 eV and above, respectively, 0.1g/cm^{3} and 0.3g/cm^{3}, can be simply obtained by mass density scaling from the results found for deuterium. This ab initio approach allows one to consistently cover a very large domain of temperature on the domain of mass density outlined above.

Multishock Compression Properties of Warm Dense Argon

Warm dense argon was generated by a shock reverberation technique. The diagnostics of warm dense argon were performed by a multichannel optical pyrometer and a velocity interferometer system. The equations of state in the pressure-density range of 20–150 GPa and 1.9–5.3 g/cm³ from the first- to fourth-shock compression were presented. The single-shock temperatures in the range of 17.2–23.4 kK were obtained from the spectral radiance. Experimental results indicates that multiple shock-compression ratio (η_i = ρ_i/ρ₀) is greatly enhanced from 3.3 to 8.8, where ρ₀ is the initial density of argon and ρ_i (i = 1, 2, 3, 4) is the compressed density from first to fourth shock, respectively. For the relative compression ratio (η_i’ = ρ_i/ρ_i-1), an interesting finding is that a turning point occurs at the second shocked states under the conditions of different experiments, and η_i’ increases with pressure in lower density regime and reversely decreases with pressure in higher density regime. The evolution of the compression ratio is controlled by the excitation of internal degrees of freedom, which increase the compression, and by the interaction effects between particles that reduce it. A temperature-density plot shows that current multishock compression states of argon have distributed into warm dense regime.

Development of Path Integral Monte Carlo Simulations with Localized Nodal Surfaces for Second-Row Elements

We extend the applicability range of fermionic path integral Monte Carlo simulations to heavier elements and lower temperatures by introducing various localized nodal surfaces. Hartree-Fock nodes yield the most accurate prediction for pressure and internal energy, which we combine with the results from density functional molecular dynamics simulations to obtain a consistent equation of state for hot, dense silicon under plasma conditions and in the regime of warm dense matter (2.3-18.6  g cm(-3), 5.0×10(5)-1.3×10(8)  K). The shock Hugoniot curve is derived and the structure of the fluid is characterized with various pair correlation functions.

First-principles equation of state of polystyrene and its effect on inertial confinement fusion implosions

Obtaining an accurate equation of state (EOS) of polystyrene (CH) is crucial to reliably design inertial confinement fusion (ICF) capsules using CH/CH-based ablators. With first-principles calculations, we have investigated the extended EOS of CH over a wide range of plasma conditions (ρ=0.1to100g/cm(3) and T=1000 to 4,000,000 K). When compared with the widely used SESAME-EOS table, the first-principles equation of state (FPEOS) of CH has shown significant differences in the low-temperature regime, in which strong coupling and electron degeneracy play an essential role in determining plasma properties. Hydrodynamic simulations of cryogenic target implosions on OMEGA using the FPEOS table of CH have predicted ∼30% decrease in neutron yield in comparison with the usual SESAME simulations. This is attributed to the ∼5% reduction in implosion velocity that is caused by the ∼10% lower mass ablation rate of CH predicted by FPEOS. Simulations using CH-FPEOS show better agreement with measurements of Hugoniot temperature and scattered light from ICF implosions.

Effective potentials of interactions and thermodynamic properties of a nonideal two-temperature dense plasma

In this article a dense nonideal, nonisothermal plasma is considered. New effective screened interaction potentials taking into account quantum-mechanical diffraction and symmetry effects have been obtained. The effective potential of the ion-ion interaction in plasmas with a strongly coupled ion subsystem and semiclassical electron subsystem is presented. Based on the obtained effective potentials the analytical expressions for internal energy and the pressure of the considered plasma were obtained.

Molecular-Atomic Transition along the Deuterium Hugoniot Curve with Coupled Electron-Ion Monte Carlo Simulations

We have performed simulations of the principal deuterium Hugoniot curve using coupled electron-ion Monte Carlo calculations. Using highly accurate quantum Monte Carlo methods for the electrons, we study the region of maximum compression along the Hugoniot, where the system undergoes a continuous transition from a molecular fluid to a monatomic fluid. We include all relevant physical corrections so that a direct comparison to experiment can be made. Around 50 GPa we find a maximum compression of 4.85. This compression is approximately 5.5% higher than previous theoretical predictions and 15% higher than the most accurate experimental data. Thus first-principles simulations encompassing the most advanced techniques are in disagreement with the results of the best experiments.

Average-atom model combined with the hypernetted chain approximation applied to warm dense matter

We have combined the average-atom model with the hypernetted chain approximation (AAHNC) to describe the electronic and ionic structure in the warm dense matter regime. On the basis of the electronic and ionic structures, the x-ray Thomson scattering (XRTS) spectrum is calculated using the random-phase approximation. While the electronic structure is described within the average-atom model, the effects of other ions on the electronic structure are considered using an integral equation method of the theory of liquids, namely the hypernetted chain approximation. The ion-ion pair potential is calculated using the modified Gordon-Kim model based on the electronic density distribution. Finally, the electronic and ionic structures are determined self-consistently. The XRTS spectrum is calculated according to the Chihara formula, where the scattering contributions are divided into three components: elastic, bound-free, and free-free. Comparison of the present AAHNC results with other theoretical models and experimental data shows very good agreement. Thus the AAHNC model can give a reasonable description of the electronic and ionic structure in warm dense matter.

Combining Kohn-Sham and orbital-free density-functional theory for Hugoniot calculations to extreme pressures

The shock Hugoniot for lithium 6 deuteride ((6)LiD) was calculated via first principles using Kohn-Sham density-functional theory molecular dynamics (KSMD) for temperatures of 0.5-25 eV. The upper limit of 25 eV represents a practical limit where KSMD is no longer computationally feasible due to the number of electronic bands which are required to be populated. To push the Hugoniot calculations to higher temperatures we make use of orbital-free density-functional theory molecular dynamics (OFMD). Thomas-Fermi-Dirac-based OFMD gives a poor description of the electronic structure at low temperatures so the initial state is not well defined. We propose a method of bootstrapping the Hugoniot from OFMD to the Hugoniot from KSMD between 10 and 20 eV, where the two methods are in agreement. The combination of KSMD and OFMD allows construction of a first-principles Hugoniot from the initial state to 1000 eV. Theoretical shock-compression results are in good agreement with available experimental data and exhibit the appropriate high-temperature limits. We show that a unified KSMD-OFMD Hugoniot can be used to assess the quality of the existing equation-of-state (EOS) models and inform better EOS models based on justifiable physics.

First-principles opacity table of warm dense deuterium for inertial-confinement-fusion applications

Accurate knowledge of the optical properties of a warm dense deuterium-tritium (DT) mixture is important for reliable design of inertial confinement fusion (ICF) implosions using radiation-hydrodynamics simulations. The opacity of a warm dense DT shell essentially determines how much radiation from hot coronal plasmas can be deposited in the DT fuel of an imploding capsule. Even for the simplest species of hydrogen, the accurate calculation of their opacities remains a challenge in the warm-dense matter regime because strong-coupling and quantum effects play an important role in such plasmas. With quantum-molecular-dynamics (QMD) simulations, we have derived a first-principles opacity table (FPOT) of deuterium (and the DT mixture by mass scaling) for a wide range of densities from ρ(D)=0.5 to 673.518g/cm(3) and temperatures from T=5000K up to the Fermi temperature T(F) for each density. Compared with results from the astrophysics opacity table (AOT) currently used in our hydrocodes, the FPOT of deuterium from our QMD calculations has shown a significant increase in opacity for strongly coupled and degenerate plasma conditions by a factor of 3-100 in the ICF-relevant photon-energy range. As conditions approach those of classical plasma, the opacity from the FPOT converges to the corresponding values of the AOT. By implementing the FPOT of deuterium and the DT mixture into our hydrocodes, we have performed radiation-hydrodynamics simulations for low-adiabat cryogenic DT implosions on the OMEGA laser and for direct-drive-ignition designs for the National Ignition Facility. The simulation results using the FPOT show that the target performance (in terms of neutron yield and energy gain) could vary from ∼10% up to a factor of ∼2 depending on the adiabat of the imploding DT capsule; the lower the adiabat, the more variation is seen in the prediction of target performance when compared to the AOT modeling.

A density functional tight binding model with an extended basis set and three-body repulsion for hydrogen under extreme thermodynamic conditions

We present a new DFTB-p3b density functional tight binding model for hydrogen at extremely high pressures and temperatures, which includes a polarizable basis set (p) and a three-body environmentally dependent repulsive potential (3b). We find that use of an extended basis set is necessary under dissociated liquid conditions to account for the substantial p-orbital character of the electronic states around the Fermi energy. The repulsive energy is determined through comparison to cold curve pressures computed from density functional theory (DFT) for the hexagonal close-packed solid, as well as pressures from thermally equilibrated DFT-MD simulations of the liquid phase. In particular, we observe improved agreement in our DFTB-p3b model with previous theoretical and experimental results for the shock Hugoniot of hydrogen up to 100 GPa and 25000 K, compared to a standard DFTB model using pairwise interactions and an s-orbital basis set, only. The DFTB-p3b approach discussed here provides a general method to extend the DFTB method for a wide variety of materials over a significantly larger range of thermodynamic conditions than previously possible.

Properties of warm dense polystyrene plasmas along the principal Hugoniot

Polystyrene (CH) is often chosen as the ablator material for inertial confinement fusion (ICF) targets. Its static, dynamical, and optical properties in warm, dense conditions (due to shock compression) are important for ICF designs. Using the first-principles quantum molecular dynamics (QMD) method, we have investigated the equation of state (EOS) and optical reflectivity of shock-compressed CH up to an unprecedentedly high pressure of 62 Mbar along the principal Hugoniot. The QMD results are compared with existing experimental measurements as well as the SESAME EOS model. Although the Hugoniot pressure and/or temperature from QMD calculations agrees with experiments and the SESAME EOS model at low pressures below 10 Mbar, we have identified for the first time a stiffer behavior of shocked CH at higher pressures (>10 Mbar). Such a stiffer behavior of warm, dense CH can affect the ablation pressure (shock strength), shock coalescence dynamics, and nonuniformity growth in ICF implosions. In addition, we corrected the mistake made in literature for calculating the reflectivity of shocked CH and obtained good agreements with experimental measurements, which should lend credence to future opacity calculations in a first-principles fashion.

First-principles thermal conductivity of warm-dense deuterium plasmas for inertial confinement fusion applications

Thermal conductivity (κ) of both the ablator materials and deuterium-tritium (DT) fuel plays an important role in understanding and designing inertial confinement fusion (ICF) implosions. The extensively used Spitzer model for thermal conduction in ideal plasmas breaks down for high-density, low-temperature shells that are compressed by shocks and spherical convergence in imploding targets. A variety of thermal-conductivity models have been proposed for ICF hydrodynamic simulations of such coupled and degenerate plasmas. The accuracy of these κ models for DT plasmas has recently been tested against first-principles calculations using the quantum molecular-dynamics (QMD) method; although mainly for high densities (ρ > 100 g/cm3), large discrepancies in κ have been identified for the peak-compression conditions in ICF. To cover the wide range of density-temperature conditions undergone by ICF imploding fuel shells, we have performed QMD calculations of κ for a variety of deuterium densities of ρ = 1.0 to 673.518 g/cm3, at temperatures varying from T = 5 × 103 K to T = 8 × 106 K. The resulting κQMD of deuterium is fitted with a polynomial function of the coupling and degeneracy parameters Γ and θ, which can then be used in hydrodynamic simulation codes. Compared with the "hybrid" Spitzer-Lee-More model currently adopted in our hydrocode lilac, the hydrosimulations using the fitted κQMD have shown up to ∼20% variations in predicting target performance for different ICF implosions on OMEGA and direct-drive-ignition designs for the National Ignition Facility (NIF). The lower the adiabat of an imploding shell, the more variations in predicting target performance using κQMD. Moreover, the use of κQMD also modifies the shock conditions and the density-temperature profiles of the imploding shell at early implosion stage, which predominantly affects the final target performance. This is in contrast to the previous speculation that κQMD changes mainly the inside ablation process during the hot-spot formation of an ICF implosion.

Nuclear quantum effects and nonlocal exchange-correlation functionals applied to liquid hydrogen at high pressure

Using first-principles molecular dynamics, we study the influence of nuclear quantum effects (NQEs) and nonlocal exchange-correlation density functionals (DFs) near molecular dissociation in liquid hydrogen. NQEs strongly influence intramolecular properties, such as bond stability, and are thus an essential part of the dissociation process. Moreover, by including DFs that account for either the self-interaction error or dispersion interactions, we find a much better description of molecular dissociation and metallization than previous studies based on classical protons and/or local or semilocal DFs. We obtain excellent agreement with experimentally measured optical properties along Hugoniot curves for precompressed states, and while we still find a first-order liquid-liquid transition at low temperatures, transition pressures are increased by more than 100 GPa.

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.

Fully variational average atom model with ion-ion correlations

An average atom model for dense ionized fluids that includes ion correlations is presented. The model assumes spherical symmetry and is based on density functional theory, the integral equations for uniform fluids, and a variational principle applied to the grand potential. Starting from density functional theory for a mixture of classical ions and quantum mechanical electrons, an approximate grand potential is developed, with an external field being created by a central nucleus fixed at the origin. Minimization of this grand potential with respect to electron and ion densities is carried out, resulting in equations for effective interaction potentials. A third condition resulting from minimizing the grand potential with respect to the average ion charge determines the noninteracting electron chemical potential. This system is coupled to a system of point ions and electrons with an ion fixed at the origin, and a closed set of equations is obtained. Solution of these equations results in a self-consistent electronic and ionic structure for the plasma as well as the average ionization, which is continuous as a function of temperature and density. Other average atom models are recovered by application of simplifying assumptions.

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.

Change in inertial confinement fusion implosions upon using an ab initio multiphase DT equation of state

Improving the description of the equation of state (EOS) of deuterium-tritium (DT) has recently been shown to change significantly the gain of an inertial confinement fusion target [S. X. Hu et al., Phys. Rev. Lett. 104, 235003 (2010)]. Here we use an advanced multiphase EOS, based on ab initio calculations, to perform a full optimization of the laser pulse shape with hydrodynamic simulations starting from 19 K in DT ice. The thermonuclear gain is shown to be a robust estimate over possible uncertainties of the EOS. Two different target designs are discussed, for shock ignition and self-ignition. In the first case, the areal density and thermonuclear energy can be recovered by slightly increasing the laser energy. In the second case, a lower in-flight adiabat is needed, leading to a significant delay (3 ns) in the shock timing of the implosion.

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.

Strong coupling and degeneracy effects in inertial confinement fusion implosions

Accurate knowledge about the equation of state (EOS) of deuterium is critical to inertial confinement fusion (ICF). Low-adiabat ICF implosions routinely access strongly coupled and degenerate plasma conditions. Using the path integral Monte Carlo method, we have derived a first-principles EOS (FPEOS) table of deuterium. It is the first ab initio EOS table which completely covers typical ICF implosion trajectory in the density and temperature ranges of ρ=0.002-1596  g/cm3 and T=1.35  eV-5.5  keV. Discrepancies in internal energy and pressure have been found in strongly coupled and degenerate regimes with respect to SESAME EOS. Hydrodynamics simulations of cryogenic ICF implosions using the FPEOS table have indicated significant differences in peak density, areal density (ρR), and neutron yield relative to SESAME simulations.