Equation of State and First Principles Prediction of the Vibrational Matrix Shift of Solid Parahydrogen

Abstract

We generate the equation of state (EOS) of solid parahydrogen using a path-integral Monte Carlo (PIMC) simulation based on a highly accurate first-principles adiabatic hindered rotor (AHR) potential energy curve for the parahydrogen dimer. The EOS curves for the fcc and hcp structures of solid parahydrogen near the equilibrium density show that the hcp structure is the more stable of the two, in agreement with experiment. To accurately reproduce the structural and energy properties of solid parahydrogen, we eliminated by extrapolation the systematic errors associated with the choice of simulation parameters used in the PIMC calculation. We also investigate the temperature dependence of the EOS curves, and the invariance of the equilibrium density with temperature is satisfyingly reproduced. The pressure as a function of density, and the compressibility as a function of pressure, are both calculated using the obtained EOS and are compared with previous simulation results and experiments. We also report the first ever a priori prediction of a vibrational matrix shift from first-principles two-body potential functions, and its result for the equilibrium state agrees well with experiment.

References in corpus (6)

• A six-dimensional H2-H2 potential energy surface for bound state spectroscopy
• A molecular superfluid: non-classical rotations in doped para-hydrogen clusters
• Microscopic molecular superfluid response: theory and simulations
• Path integral Monte Carlo ground state approach: Formalism, implementation, and applications
• Diffusion Monte Carlo study of the equation of state of solid ortho-D${}_2$
• The Physics of Liquid Para-Hydrogen

Cited by in corpus (5)

• Equation of state of solid parahydrogen using ab initio two-body and three-body interaction potentials
• Three-body potential energy surface for parahydrogen
• 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