Assessing the accuracy of the Jastrow antisymmetrized geminal power in the H₄ model system

We report a quantum Monte Carlo study, on a very simple but nevertheless very instructive model system of four hydrogen atoms, recently proposed in Gasperich et al. [J. Chem. Phys. 147, 074106 (2017)]. We find that the Jastrow correlated Antisymmetrized Geminal Power (JAGP) is able to recover most of the correlation energy even when the geometry is symmetric and the hydrogens lie on the edges of a perfect square. Under such conditions, the diradical character of the molecule ground state prevents a single determinant Ansatz to achieve an acceptable accuracy, whereas the JAGP performs very well for all geometries. Remarkably, this is obtained with a similar computational effort. Moreover, we find that the Jastrow factor is fundamental in promoting the correct resonances among several configurations in the JAGP, which cannot show up in the pure Antisymmetrized Geminal Power (AGP). We also show the extremely fast convergence of this approach in the extension of the basis set. Remarkably, only the simultaneous optimization of the Jastrow and the AGP part of our variational Ansatz is able to recover an almost perfect nodal surface, yielding therefore state of the art energies, almost converged in the complete basis set limit, when the so called diffusion Monte Carlo is applied.