Ionization potential depression and charge state of warm dense hydrogen from ab initio path integral Monte Carlo simulations

Research into Warm Dense Matter (WDM) has become more important with recent advances in inertial confinement fusion and astrophysics. The interplay between quantum degeneracy and Coulomb interactions, and the transition of condensed and plasma phases occurring under these conditions, however, make WDM extremely difficult to describe theoretically. Several methods exist to describe matter at these conditions, with recent extensions to Path Integral Monte Carlo (PIMC) allowing the calculation of the Laplace transform of the dynamic structure factor, i.e. the imaginary time correlation function (ITCF), of warm dense hydrogen [1]. While PIMC is quasi-exact, it does not give access to other important physical quantities describing a plasma state, like the ionisation and ionisation potential depression (IPD). Moreover, both are difficult to measure experimentally. To remedy this, we instead compare the ITCF from PIMC simulations with synthetic X-Ray Thomson Scattering spectra [2], computed from a Chihara decomposition [3], to obtain a best estimate for the ionisation state. The IPD is then directly calculated using the Saha equation and compared against other commonly used models. We expect this work to be relevant for future inertial confinement energy developments, particularly in validating equation of state models, and for the refinement of astrophysical models.

[1] T. Dornheim et.al., arXiv preprint arXiv:2403.08570 (2024)

[2] T. Dornheim et.al., Phys. Plasmas 30, 042707 (2023)

[3] G. Gregori et.al., Phys. Rev. E 67, 026412 (2003)