Shock compression of deuterium near 100 GPa pressures

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.

Measurements of the equations of state and spectrum of nonideal xenon plasma under shock compression

Experimental equations of state on generation of nonideal xenon plasma by intense shock wave compression was presented in the ranges of pressure of 2-16 GPa and temperature of 31-50 kK, and the xenon plasma with the nonideal coupling parameter Γ range from 0.6-2.1 was generated. The shock wave was produced using the flyer plate impact and accelerated up to ∼6 km/s with a two-stage light gas gun. Gaseous specimens were shocked from two initial pressures of 0.80 and 4.72 MPa at room temperature. Time-resolved spectral radiation histories were recorded by using a multiwavelength channel pyrometer. The transient spectra with the wavelength range of 460-700 nm were recorded by using a spectrometer to evaluate the shock temperature. Shock velocity was measured and particle velocity was determined by the impedance matching methods. The equations of state of xenon plasma and ionization degree have been discussed in terms of the self-consistent fluid variational theory.

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.

Excited electron dynamics modeling of warm dense matter

We present a model (the electron force field, or eFF) based on a simplified solution to the time-dependent Schrödinger equation that with a single approximate potential between nuclei and electrons correctly describes many phases relevant for warm dense hydrogen. Over a temperature range of 0 to 100,000 K and densities up to 1 g/cm(3), we find excellent agreement with experimental, path integral Monte Carlo, and linear mixing equations of state, as well as single-shock Hugoniot curves from shock compression experiments. In principle eFF should be applicable to other warm dense systems as well.

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.

Monte Carlo results for the hydrogen Hugoniot

We propose a theoretical Hugoniot relation obtained by combining results for the equation of state from the direct path integral Monte Carlo technique (DPIMC) and those from reaction ensemble Monte Carlo (REMC) simulations. The main idea of this proposal is based on the fact that the DPMIC technique provides first-principle results for a wide range of densities and temperatures including the region of partially ionized plasmas. On the other hand, for lower temperatures where the formation of molecules becomes dominant, DPIMC simulations become cumbersome and inefficient. For this region it is possible to use accurate REMC simulations where bound states (molecules) are treated on the Born-Oppenheimer level. The remaining interaction is then reduced to the scattering between neutral particles which is reliably treated classically by applying effective potentials. The resulting Hugoniot is located between the experimental values of Knudson et al. [Phys. Rev. Lett. 87, 225501 (2001)] and Collins et al. [Science 281, 1178 (1998)].

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.

Equation of state for weakly coupled quantum plasmas

We calculate thermodynamic properties for a dense hydrogen plasma and a quantum electron gas using thermodynamic Green's function techniques. Our perturbation approach is appropriate to give reliable results in the weak coupling regime. In particular, the contribution of the exchange term of the order e(4) is fully included for the nondegenerate case as well as for the dense highly degenerate quantum region. We compare our results for the equation of state with data obtained by different numerical simulations.

Insulator to metal transition in fluid deuterium

We have investigated the insulator to metal transition in fluid deuterium using first principles simulations. Both density functional and quantum Monte Carlo calculations indicate that the electronic energy gap of the liquid vanishes at about ninefold compression and 3000 K. At these conditions the computed conductivity values are characteristic of a poor metal. These findings are consistent with those of recent shock wave experiments but the computed conductivity is larger than the measured value. From our ab initio results we conclude that the transition is driven by molecular dissociation rather than disorder and that both temperature and pressure play a key role in determining structural changes in the fluid.

Metallization of fluid nitrogen and the mott transition in highly compressed low-Z fluids

Electrical conductivities are reported for degenerate fluid nitrogen at pressures up to 180 GPa (1.8 Mbar) and temperatures of approximately 7000 K. These extreme quasi-isentropic conditions were achieved with multiple-shock compression generated with a two-stage light-gas gun. Nitrogen undergoes a nonmetal-metal transition at 120 GPa, probably in the monatomic state. These N data and previous conductivity data for H, O, Cs, and Rb are used to develop a general picture of the systematics of the nonmetal-metal transition in these fluids. Specifically, the density dependences of electrical conductivities in the semiconducting fluid are well correlated with the radial extent of the electronic charge-density distributions of H, N, O, Cs, and Rb atoms. These new data for N scale with previous data for O, as expected from their similar charge-density distributions.

Raman spectroscopy of hot dense hydrogen

High P-T Raman measurements of solid and fluid hydrogen to above 1100 K at 70 GPa and to above 650 K in 150 GPa range, conditions previously inaccessible by static compression experiments, provide new insight into the behavior of the material under extreme conditions. The data give a direct measure of the melting curve that extends previous optical investigations by up to a factor of 4 in pressure. The magnitude of the vibron frequency temperature derivative (dnu/dT)(P) increases by a factor of approximately 30 over the measured pressure range, indicating an increase in intrinsic anharmonicity and weakening of the molecular bond.

Equation-of-state measurements of polyimide at pressures up to 5.8 TPa using low-density foam with laser-driven shock waves

The laser-driven equation-of-state (EOS) experiments for polyimide are presented. The experiments were performed with emission measurements from the rear sides of shocked targets at up to a laser intensity of 10(14) W/cm(2) or higher with 351 nm wavelength and 2.5 ns duration. Polyimide Hugoniot data were obtained up to 0.6 TPa with good accuracy. Applying low-density foam ablator to the EOS unknown material, we also obtained the data at a highest pressure of 5.8 TPa in the nonmetal materials. Those data were in agreement with the theoretical curves.