Assessing many-body methods on the potential energy surface of the (H2)2 hydrogen dimer

The anisotropic potential energy surface of the (H2)2 dimer represents a challenging problem for many-body methods. Here, we determine the potential energy curves of five different dimer configurations (T, Z, X, H, and L) using the lattice regularized diffusion Monte Carlo method and a number of approximate functionals within density functional theory (DFT), including advanced orbital-dependent functionals based on the random phase approximation (RPA). We assess their performance in describing the potential wells, bond distances, and relative energies. The repulsive potential wall is studied by looking at the relative stability of the different dimer configurations as a function of an applied force acting along the intermolecular axis. It is shown that most functionals within DFT break down at finite compression, even those that give an accurate description around the potential well minima. Only by including exchange within RPA, a qualitatively correct description along the entire potential energy curve is obtained. Finally, we discuss these results in the context of solid molecular hydrogen at finite pressures.