All-Electron Path Integral Monte Carlo Simulations of Warm Dense Matter: Application to Water and Carbon Plasmas

We develop an all-electron path integral Monte Carlo method with free-particle nodes for warm dense matter and apply it to water and carbon plasmas. We thereby extend path integral Monte Carlo studies beyond hydrogen and helium to elements with core electrons. Path integral Monte Carlo results for pressures, internal energies, and pair-correlation functions compare well with density functional theory molecular dynamics calculations at temperatures of (2.5--7.5)×10⁵ K, and both methods together form a coherent equation of state over a density-temperature range of 3--12 g/cm³ and 10⁴--10⁹ K.