Development of Path Integral Monte Carlo Simulations with Localized Nodal Surfaces for Second-Row Elements

Abstract

We extend the applicability range of fermionic path integral Monte Carlo simulations to heavier elements and lower temperatures by introducing various localized nodal surfaces. Hartree-Fock nodes yield the most accurate prediction for pressure and internal energy that we combine with the results from density functional molecular dynamics simulations to obtain a consistent equation of state for hot, dense silicon under plasma conditions and in the regime of warm dense matter (2.3-18.6 g/cm, 5.0*10^5 - 1.3*10^8 K). The shock Hugoniot curve is derived and the structure of the fluid is characterized with various pair correlation functions.