Simulations of fluid hydrogen: comparison of a dissociation model with tight-binding molecular dynamics

Wave packet simulation of dense hydrogen

Dense hydrogen is studied in the framework of wave packet molecular dynamics. In this semiquantal many-body simulation method the electrons are represented by wave packets which are suitably parametrized. The equilibrium properties and time evolution of the system are obtained with the help of a variational principle. At room temperature the results for the isotherms are in good agreement with anvil experiments. At higher densities beyond the range of the experimental data a transition from a molecular to a metallic state is predicted. The wave packets become delocalized and the electrical conductivity increases sharply. The phase diagram is calculated in a wide range of the pressure-density-temperature space. The observed transition from the molecular to metallic state is accompanied by an increase in density in agreement with recent reverberating shock wave experiments.

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.

Local-spin-density-approximation molecular-dynamics simulations of dense deuterium

Local-spin-density-approximation molecular-dynamics simulations of deuterium in the dissociating regime are presented, with a particular emphasis on the molecular phase of two isochores corresponding for deuterium to V=6 cm(3)/mole, rho=0.670 g/cm(3) and V=4 cm(3)/mole, rho=1 g/cm(3). It is shown that the transition from the molecular regime, well described by the local-spin-density-approximation functional, to the dissociated regime where previous local-density-approximation results are recovered, comes with a negative curvature deltaP/deltaT<0 in the isochore. We show that this effect is not enough to explain the large compressibility measured in the laser experiments [L. B. DaSilva et al., Phys. Rev. Lett. 78, 483 (1997); G. W. Collins et al., Science 281, 1178 (1998); P. Celliers et al., Phys. Rev. Lett. 84, 5564 (2000)].