Equation of state of partially ionized hydrogen and deuterium plasma revisited

Ionization calculations using classical molecular dynamics

By performing an ensemble of molecular dynamics simulations, the model-dependent ionization state is computed for strongly interacting systems self-consistently. This is accomplished through a free energy minimization framework based on the technique of thermodynamic integration. To illustrate the method, two simple models applicable to partially ionized hydrogen plasma are presented in which pair potentials are employed between ions and neutral particles. Within the models, electrons are either bound in the hydrogen ground state or distributed in a uniform charge-neutralizing background. Particular attention is given to the transition between atomic gas and ionized plasma, where the effect of neutral interactions is explored beyond commonly used models in the chemical picture. Furthermore, pressure ionization is observed when short-range repulsion effects are included between neutrals. The developed technique is general, and we discuss the applicability to a variety of molecular dynamics models for partially ionized warm dense matter.

Ab Initio Path Integral Monte Carlo Simulations of the Uniform Electron Gas on Large Length Scales

The accurate description of non-ideal quantum many-body systems is of prime importance for a host of applications within physics, quantum chemistry, materials science, and related disciplines. At finite temperatures, the gold standard is given by ab initio path integral Monte Carlo (PIMC) simulations, which do not require any empirical input but exhibit an exponential increase in the required computation time for Fermionic systems with an increase in system size N. Very recently, computing Fermionic properties without this bottleneck based on PIMC simulations of fictitious identical particles has been suggested. In our work, we use this technique to perform very large (N ≤ 1000) PIMC simulations of the warm dense electron gas and demonstrate that it is capable of providing a highly accurate description of the investigated properties, i.e., the static structure factor, the static density response function, and the local field correction, over the entire range of length scales.

Virial coefficients of the uniform electron gas from path-integral Monte Carlo simulations

The properties of plasmas in the low-density limit are described by virial expansions. Analytical expressions are known from Green's function approaches only for the first three virial coefficients. Accurate path-integral Monte Carlo (PIMC) simulations have recently been performed for the uniform electron gas, allowing the virial expansions to be analyzed and interpolation formulas to be derived. The exact expression for the second virial coefficient is used to test the accuracy of the PIMC simulations and the range of validity of the interpolation formula of Groth et al. [Phys. Rev. Lett. 119, 135001 (2017)0031-900710.1103/PhysRevLett.119.135001], and we discuss the fourth virial coefficient, which is not exactly known yet. Combining PIMC simulations with benchmarks from exact virial expansion results would allow us to obtain more accurate representations of the equation of state for parameter ranges of conditions which are of interest, e.g., for helioseismology.