Reflected shock experiments on the equation-of-state properties of liquid deuterium at 100-600 GPa (1-6 mbar)

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.

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.

Equation of State and Electrical Conductivity of Dense Fluid Hydrogen and Helium

The equation of state of fluid hydrogen, helium, and their mixtures is determined within fluid variational theory. Reactions between the constituents such as dissociation and ionization are considered. Results are given for densities and temperatures relevant for the interior of giant planets. Furthermore, the electrical conductivity is determined within linear response theory. Comparison is performed with available experiments and other theoretical work.

Studies of plastic-ablator compressibility for direct-drive inertial confinement fusion on OMEGA

The compression of planar plastic targets was studied with x-ray radiography in the range of laser intensities of I approximately 0.5 to 1.5x10(15) W/cm2 using square (low-compression) and shaped (high-compression) pulses. Two-dimensional simulations with the radiative hydrocode DRACO show good agreement with measurements at laser intensities up to I approximately 10(15) W/cm2. These results provide the first experimental evidence for low-entropy, adiabatic compression of plastic shells in the laser intensity regime relevant to direct-drive inertial confinement fusion. A density reduction near the end of the drive at a high intensity of I approximately 1.5x10(15) W/cm2 has been correlated with the hard x-ray signal caused by hot electrons from two-plasmon-decay instability.

Shock wave propagation in dissociating low-Z liquids: D2

We present direct molecular dynamics simulations of shock wave propagation in liquid deuterium for a wide range of impact velocities. The calculated Hugoniot is in perfect agreement with the gas-gun data as well as with the most recent experimental data. At high impact velocities we observe a smearing of the shock wave front and propagation of fast dissociated molecules well ahead of the compressed region. This smearing occurs due to the fast deuterium dissociation at the shock wave front. The experimental results are discussed in view of this effect.

Reaction ensemble Monte Carlo technique and hypernetted chain approximation study of dense hydrogen

In spite of the simple structure of hydrogen, up to now there is no unified theoretical and experimental description of hydrogen at high pressures. Recent results of Z-pinch experiments show a large deviation from those obtained by laser driven ones. Theoretical investigations including ab initio computer simulations show considerable differences at such extreme conditions from each other and from experimental values. We apply the reaction ensemble Monte Carlo technique on one hand and a combination of the hypernetted chain approximation with the mass action law on the other to study the behavior of dense hydrogen at such conditions. The agreement between both methods for the equation of state and for the Hugoniot curve is excellent. Comparison to other methods and experimental results is also performed.

Use of a wave reverberation technique to infer the density compression of shocked liquid deuterium to 75 GPa

A novel approach was developed to probe density compression of liquid deuterium (L-D2) along the principal Hugoniot. Relative transit times of shock waves reverberating within the sample are shown to be sensitive to the compression due to the first shock. This technique has proven to be more sensitive than the conventional method of inferring density from the shock and mass velocity, at least in this high-pressure regime. Results in the range of 22-75 GPa indicate an approximately fourfold density compression, and provide data to differentiate between proposed theories for hydrogen and its isotopes.

Calculation of a deuterium double shock Hugoniot from ab initio simulations

We calculate the equation of state of dense deuterium with two ab initio simulation techniques, path integral Monte Carlo and density functional theory molecular dynamics, in the density range of 0.67 < or = rho < or = 1.60 g cm(-3). We derive the double shock Hugoniot and compare with the recent laser-driven double shock wave experiments by Mostovych et al. [Phys. Rev. Lett. 85, 3870 (2000)]. We find excellent agreement between the two types of microscopic simulations, but a significant discrepancy with the laser-driven shock measurements.

Equation of state measurements in liquid deuterium to 70 GPa

Using intense magnetic pressure, a method was developed to launch flyer plates to velocities in excess of 20 km/s. This technique was used to perform plate-impact, shock wave experiments on cryogenic liquid deuterium ( L-D(2)) to examine its high-pressure equation of state. Using an impedance matching method, Hugoniot measurements were obtained in the pressure range of 30-70 GPa. The results of these experiments disagree with previously reported Hugoniot measurements of L-D(2) in the pressure range above approximately 40 GPa, but are in good agreement with first principles, ab initio models for hydrogen and its isotopes.