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

Structural and transport properties of ammonia along the principal Hugoniot

We investigate, via quantum molecular dynamics simulations, the structural and transport properties of ammonia along the principal Hugoniot for temperatures up to 10 eV and densities up to 2.6 g/cm³. With the analysis of the molecular dynamics trajectories by use of the bond auto-correlation function, we identify three distinct pressure-temperature regions for local chemical structures of ammonia. We derive the diffusivity and viscosity of strong correlated ammonia with high accuracy through fitting the velocity and stress-tensor autocorrelation functions with complex functional form which includes structures and multiple time scales. The statistical error of the transport properties is estimated. It is shown that the diffusivity and viscosity behave in a distinctly different manner at these three regimes and thus present complex features. In the molecular fluid regime, the hydrogen atoms have almost the similar diffusivity as nitrogen and the viscosity is dominated by the kinetic contribution. When entering into the mixture regime, the transport behavior of the system remarkably changes due to the stronger ionic coupling, and the viscosity is determined to decrease gradually and achieve minimum at about 2.0 g/cm³ on the Hugoniot. In the plasma regime, the hydrogen atoms diffuse at least twice as fast as the nitrogen atoms.

Transport properties of hydrogen-helium mixtures at extreme density and temperature conditions

We perform a systematic study of hydrogen-helium mixtures using quantum molecular dynamics (QMD) with a focus on the equations of state and structural and transport properties such as electrical conductivity, diffusion, and viscosity at conditions of giant planet interiors of 0.2∼2.3 g/cm(3) and 1000∼80,000 K for a typical helium mass fraction of 0.245. The H-He separation is found at low temperatures by analyzing the trajectories and pair distribution functions. We show that the diffusion coefficients exhibit transitions from kinetics- to potential-, and then to demixing-dominated regimes. In addition, we identify the discontinuity feature of optical absorption of a H-He mixture at low density and temperature conditions, which results from the change from an intraband to an interband transition. The Stokes-Einstein relation between the diffusion and viscosity coefficients is also discussed.

Quantum molecular dynamics study of expanded beryllium: evolution from warm dense matter to atomic fluid

By performing quantum molecular dynamics (QMD) simulations, we investigate the equation of states, electrical and optical properties of the expanded beryllium at densities two to one-hundred lower than the normal solid density, and temperatures ranging from 5000 to 30000 K. With decreasing the density of Be, the optical response evolves from the one characteristic of a simple metal to the one of an atomic fluid. By fitting the optical conductivity spectra with the Drude-Smith model, it is found that the conducting electrons become localized at lower densities. In addition, the negative derivative of the electrical resistivity on temperature at density about eight lower than the normal solid density demonstrates that the metal to nonmetal transition takes place in the expanded Be. To interpret this transition, the electronic density of states is analyzed systematically. Furthermore, a direct comparison of the Rosseland opacity obtained by using QMD and the standard opacity code demonstrates that QMD provides a powerful tool to validate plasma models used in atomic physics approaches in the warm dense matter regime.

Thomas-Fermi Z-scaling laws and coupling stabilization for plasmas

Extending the well-known Thomas-Fermi Z-scaling laws to the Coulomb coupling parameter, we investigate the stabilization of the ionic coupling in isochoric heating [Clérouin et al., Phys. Rev. E 87, 061101 (2013)]. This stabilization is restricted to a domain in atomic number Z, temperature, and density, including strong limitations on high couplings, that can only be obtained for high-Z elements. Contact is made with recent isochoric heating experiments. The consequences for corresponding states with respect to ionic coupling are also quantified via orbital free molecular dynamics simulations. This opens avenues for future isochoric heating experiments.

Ab initio molecular dynamics simulations of dense boron plasmas up to the semiclassical Thomas-Fermi regime

We build an "all-electron" norm-conserving pseudopotential for boron which extends the use of ab initio molecular dynamics simulations up to 50 times the normal density rho0. This allows us to perform ab initio simulations of dense plasmas from the regime where quantum mechanical effects are important to the regime where semiclassical simulations based on the Thomas-Fermi approach are, by default, the only simulation method currently available. This study first allows one to establish, for the case of boron, the density regime from which the semiclassical Thomas-Fermi approach is legitimate and sufficient. It further brings forward various issues pertaining to the construction of pseudopotentials aimed at high-pressure studies.

A crossover from metal to plasma in dense fluid hydrogen

Thermodynamic properties in dense fluid hydrogen are studied by using a density-functional theory for electron-proton binary mixtures that is called quantal hypernetted-chain (QHNC) integral equation. A nonlocal approximation for the exchange-correlation potential in a finite-temperature Kohn-Sham equation is presented. Results obtained from the QHNC with the nonlocal approximation are compared with those obtained from the QHNC with a local density approximation. Temperature variation of thermodynamic quantities between 10(4) and 10(6) K are investigated along an isochor specified by a dimensionless density parameter of rs=0.5. These quantities obtained from the QHNCs show that a crossover from metal to plasma occurs around a temperature of T=1.78 x 10(5) K. Electrical resistivity Re of the dense fluid hydrogen evaluated from a Ziman formula [The Properties of Liquid Metals, edited by S. Takenohi (Wiley, New York, 1973)] extended to finite temperature is about 0.7 muOmega cm at T=10(4) K. The dense fluid hydrogen at the temperature can be considered as a metallic fluid, because the value is smaller than typical values of Re in alkali metals at room temperature. The Re slightly increases with the temperature increase, and the temperature valuation of Re is monotonic. We clearly show that the contribution from the electronic excited states plays an important role for the sharp crossover from the metal to the plasma, and that the crossover is interpreted as a crossover from degenerate electron gas to nondegenerate electron gas.

Simulations of the optical properties of warm dense aluminum

Using quantum molecular dynamics simulations, we show that the optical properties of aluminum change drastically along the nonmetal metal transition observed experimentally. As the density increases and the many-body effects become important, the optical response gradually evolves from the one characteristic of an atomic fluid to the one of a simple metal. We show that quantum molecular dynamics combined with the Kubo-Greenwood formulation naturally embodies the two limits and provides a powerful tool to calculate and benchmark the optical properties of various systems as they evolve into the warm dense matter regime.

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.

Ab initio study of deuterium in the dissociating regime: sound speed and transport properties

The sound speed and the transport properties of dense hydrogen (deuterium) are computed from local spin-density approximation molecular-dynamics simulations in the dissociating regime. The sound speed c(s) is evaluated from the thermodynamical differentiation of the equation of state in the molecular phase and is in very good agreement with recent experiments. The diffusion constant D and the viscosity eta are extracted from simulations performed at V=6, 4, and 2.7 cm(3)/mole, corresponding, respectively, for deuterium at rho=0.672, 1.0, and 1.5 g/cm(3) in a range of temperatures 1000 K<T<50,000 K. In the dissociated regime, the diffusion coefficient is well predicted by one-component plasma formulas using a renormalized coupling parameter recently proposed by Murillo [M. S. Murillo, Phys. Rev. B 62, 4115 (2000)]. The behavior of the shear viscosity in the dissociated regime is more complex and exhibits a crossover between atomic and screened plasma formulation. A comparison with recent molecular-dynamics simulations of Yukawas systems shows that the inverse of the screening length must lie between 1 and 2, in nearest-neighbor radius units, as suggested by the results on the diffusion.