Capturing strong correlation in molecules with phaseless auxiliary-field quantum Monte Carlo using generalized Hartree Fock trial wavefunctions

Generalized Hartree Fock (GHF) is a long-established electronic structure method that can lower the energy (compared to spin-restricted variants) by breaking physical wavefunction symmetries, namely S^2 and Sz. After an exposition of GHF theory, we assess the use of GHF trial wavefunctions in phaseless auxiliary field quantum Monte Carlo (ph-AFQMC-G) calculations of strongly correlated molecular systems including symmetrically stretched hydrogen rings, carbon dioxide, and dioxygen. We find comparable or better accuracy than CCSD(T) with unrestricted HF and GHF references, and consistently smooth dissociation curves -- a remarkable result given the relative scalability of ph-AFQMC-G to larger system sizes. The present exploration of model strongly correlated systems marks a promising starting point for future studies of more chemically-relevant molecules, and demonstrates that ph-AFQMC-G provides a highly accurate (and, in contrast to active-space-based trials, relatively black box and always size-consistent) description of challenging systems exhibiting, e.g., antiferromagnetic coupling and/or geometric spin frustration.