Equation of state of solid parahydrogen using ab initio two-body and three-body interaction potentials

Abstract

We present the equation of state (EOS) of solid parahydrogen between $0.024{\mathring{A}}^{-3}$ and $0.1{\mathring{A}}^{-3}$ at $T=4.2$ K, calculated using path-integral Monte Carlo simulations, with ab initio two-body and three-body interaction potentials. We correct for finite size simulation errors using potential tail corrections. Trotter factorization errors are accounted for, either via extrapolation, or by using a suitably small imaginary time step. We incorporate the three-body interaction using two methods; the full inclusion method, where pair and three-body interactions are used in both Monte Carlo sampling and in the energy estimators, and the perturbative method, where three-body interactions are omitted from sampling but are still present in energy estimations. Both treatments of the three-body interaction return very similar total energies and pressures. The presence of three-body interactions has only minor effects on the structural properties of the solid. Whereas the pair interaction, on its own, significantly overestimates the pressure of solid parahydrogen, the additional presence of the three-body interaction causes a severe underestimation of the pressure. Our findings suggest that accurate simulations of solid parahydrogen require four-body and possibly higher-order many-body interactions. It may also be the case that static interaction potentials are entirely unsuitable for simulations of solid parahydrogen at high densities.

References in corpus (9)

• Worm Algorithm and Diagrammatic Monte Carlo: A New Approach to Continuous-Space Path Integral Monte Carlo Simulations
• A six-dimensional H2-H2 potential energy surface for bound state spectroscopy
• Microscopic molecular superfluid response: theory and simulations
• Path integral Monte Carlo ground state approach: Formalism, implementation, and applications
• Atomic displacements in quantum crystals
• High-purity solid parahydrogen
• Three-body potential energy surface for parahydrogen
• Equation of State and First Principles Prediction of the Vibrational Matrix Shift of Solid Parahydrogen
• Computer simulations of supercooled liquid hydrogen mixtures and the possible crystallization slowdown

Cited by in corpus (3)

• Phases of ${}^4$He and H${}_2$ adsorbed on a single carbon nanotube
• A neural network-based four-body potential energy surface for parahydrogen
• Path-integral Monte Carlo simulations of solid parahydrogen using two-body, three-body, and four-body ab initio interaction potential energy surfaces