Effect of three-body interactions on the zero-temperature equation of state of HCP solid 4He

Abstract

Previous studies have pointed to the importance of three-body interactions in high density 4He solids. However the computational cost often makes it unfeasible to incorporate these interactions into the simulation of large systems. We report the implementation and evaluation of a computationally efficient perturbative treatment of three-body interactions in hexagonal close packed solid 4He utilizing the recently developed nonadditive three-body potential of Cencek et al. This study represents the first application of the Cencek three-body potential to condensed phase 4He systems. Ground state energies from quantum Monte Carlo simulations, with either fully incorporated or perturbatively treated three-body interactions, are calculated in systems with molar volumes ranging from 21.3 cm3/mol down to 2.5 cm3/mol. These energies are used to derive the zero-temperature equation of state for comparison against existing experimental and theoretical data. The equations of state derived from both perturbative and fully incorporated three-body interactions are found to be in very good agreement with one another, and reproduce the experimental pressure-volume data with significantly better accuracy than is obtained when only two-body interactions are considered. At molar volumes below approximately 4.0 cm3/mol, neither two-body nor three-body equations of state are able to accurately reproduce the experimental pressure-volume data, suggesting that below this molar volume four-body and higher many-body interactions are becoming important.