arXiv:1905.04361 · The Journal of Chemical Physics
An auxiliary-Field quantum Monte Carlo perspective on the ground state of the dense uniform electron gas: An investigation with Hartree-Fock trial wavefunctions
We assess the utility of Hartree-Fock (HF) trial wavefunctions in performing phaseless auxiliary-field quantum Monte Carlo (ph-AFQMC) on the uniform electron gas (UEG) model. The combination of ph-AFQMC with spin-restricted HF (RHF+ph-AFQMC) was found to be highly accurate and efficient for systems containing up to 114 electrons in 2109 orbitals, particularly for rs ≤ 2.0. Compared to spin-restricted coupled-cluster (RCC) methods, we found that RHF+ph-AFQMC performs better than CC with singles, doubles, and triples (RCCSDT) and similarly to or slightly worse than CC with singles, doubles, triples, and quadruples (RCCSDTQ) for rs ≤ 3.0 in the 14-electron UEG model. With the 54-electron, we found RHF+ph-AFQMC to be nearly exact for rs ≤ 2.0 and pointed out potential biases in existing benchmarks. Encouraged by these, we performed RHF+ph-AFQMC on the 114-electron UEG model for rs ≤ 2.0 and provided new benchmark data for future method development. We found that the UEG models with rs = 5.0 remain to be challenging for RHF+ph-AFQMC. Employing nonorthogonal configuration expansions or unrestricted HF states as trial wavefunctions was also found to be ineffective in the case of the 14-electron UEG model with rs = 5.0. We emphasize the need for a better trial wavefunction for ph-AFQMC in simulating strongly correlated systems. With the 54-electron and 114-electron UEG models, we stress the potential utility of RHF+ph-AFQMC for simulating dense solids.
We assess the utility of Hartree-Fock (HF) trial wavefunctions in performing phaseless auxiliary-field quantum Monte Carlo (ph-AFQMC) on the uniform electron gas (UEG) model. The combination of ph-AFQMC with spin-restricted HF (RHF+ph-AFQMC) was found to be highly accurate and efficient for systems containing up to 114 electrons in 2109 orbitals, particularly for rs ≤ 2.0. Compared to spin-restricted coupled-cluster (RCC) methods, we found that RHF+ph-AFQMC performs better than CC with singles, doubles, and triples (RCCSDT) and similarly to or slightly worse than CC with singles, doubles, triples, and quadruples (RCCSDTQ) for rs ≤ 3.0 in the 14-electron UEG model. With the 54-electron, we found RHF+ph-AFQMC to be nearly exact for rs ≤ 2.0 and pointed out potential biases in existing benchmarks. Encouraged by these, we performed RHF+ph-AFQMC on the 114-electron UEG model for rs ≤ 2.0 and provided new benchmark data for future method development. We found that the UEG models with rs = 5.0 remain to be challenging for RHF+ph-AFQMC. Employing nonorthogonal configuration expansions or unrestricted HF states as trial wavefunctions was also found to be ineffective in the case of the 14-electron UEG model with rs = 5.0. We emphasize the need for a better trial wavefunction for ph-AFQMC in simulating strongly correlated systems. With the 54-electron and 114-electron UEG models, we stress the potential utility of RHF+ph-AFQMC for simulating dense solids.
We assess the utility of spin-restricted Hartree-Fock (RHF) trial wavefunctions in performing phaseless auxiliary-field quantum Monte Carlo (ph-AFQMC) on the uniform electron gas (UEG) model. This particular combination, RHF+ph-AFQMC, was found to be highly accurate and efficient for systems containing up to 114 electrons in 2109 orbitals, particularly for rs ≤ 2.0. Compared to spin-restricted coupled-cluster (RCC) methods, we found that RHF+ph-AFQMC performs better than CC with singles, doubles, and triples (RCCSDT) and similarly to or slightly worse than CC with singles, doubles, triples, and quadruples (RCCSDTQ) for rs ≤ 3.0 in the 14-electron UEG model. With the 54-electron, we found RHF+ph-AFQMC to be nearly exact for rs ≤ 2.0 and pointed out potential biases in existing benchmarks. Encouraged by these, we performed RHF+ph-AFQMC on the 114-electron UEG model for rs ≤ 2.0 and provided new benchmark data for future method development. As the UEG models with rs = 5.0 remain to be challenging for the RHF+ph-AFQMC, we emphasize the need for a better trial wavefunction for ph-AFQMC in simulating strongly correlated systems. With the 54-electron and 114-electron UEG models, we stress the potential utility of RHF+ph-AFQMC for simulating dense solids.
Phaseless auxiliary-field quantum Monte Carlo (ph-AFQMC) has emerged as a promising electronic structure method for correlated electronic systems. However, the quality of its predictions depends critically on the choice of trial wavefunction, and it is not obvious how to make an optimal choice especially for strongly correlated states of large systems. Mean-field wavefunctions are compelling trial wavefunction candidates as they map directly to chemical concepts and can be obtained with O(N4) cost. Yet in the strongly correlated regime one faces a symmetry dilemma and the existence of multiple nearly-degenerate solutions. In this work we investigate active space models of [2Fe-2S]2+, mixed-valent [4Fe-4S]2+, and [4Fe-4S]4+ and explore the sensitivity of ph-AFQMC to the choice of unrestricted Hartree-Fock trial wavefunction. We find that chemical properties and physical symmetries, rather than the variational energy, ought to guide the choice of mean-field trial for ph-AFQMC (or reference state for coupled cluster models), and show that surprisingly accurate ground-state energies for these systems can be obtained. However, in all cases we find a rapidly vanishing overlap between the stochastic wavefunction and the UHF trial, indicating that the trials are suboptimal importance functions. By analogy to a similar situation in the stretched helium dimer cation, we show how this sampling bias pushes ph-AFQMC towards artificially negative energies, which evidently can be compensated for by the phaseless bias in certain cases.
We investigate the viability of the phaseless finite-temperature auxiliary-field quantum Monte Carlo (ph-FT-AFQMC) method for ab initio systems using the uniform electron gas as a model. Through comparisons with exact results and FT coupled cluster theory, we find that ph-FT-AFQMC is sufficiently accurate at high to intermediate electronic densities. We show, both analytically and numerically, that the phaseless constraint at FT is fundamentally different from its zero-temperature counterpart (i.e., ph-ZT-AFQMC), and generally, one should not expect ph-FT-AFQMC to agree with ph-ZT-AFQMC in the low-temperature limit. With an efficient implementation, we are able to compare exchange-correlation energies to the existing results in the thermodynamic limit and find that the existing parameterizations are highly accurate. In particular, we found that ph-FT-AFQMC exchange-correlation energies are in better agreement with a known parameterization than is restricted path-integral MC in the regime of Θ ≤ 0.5 and rs ≤ 2, which highlights the strength of ph-FT-AFQMC.
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.
The ground state electronic energy of benzene is the focus of a recent blind test by Eriksen and co-workers [arXiv:2008.02678]. In this note, we report the phaseless auxiliary-field quantum Monte Carlo (ph-AFQMC) results on the identical problem. We examine trial wavefunctions of the spin-restricted Hartree-Fock (RHF), and complete active space self-consistent field (CASSCF) with an active space of 6-electron and 6-orbital (i.e., CAS(6,6)) form for the cc-pVDZ basis set. ph-AFQMC+RHF deviates from the value where many methods agreed on (i.e., -863 mEh) by -3.1(3) mEh whereas ph-AFQMC+CAS(6,6) deviates from the same value by -1.3(4) mEh. In addition to this, we report the frozen core correlation energy of ph-AFQMC+RHF in the cc-pVTZ and cc-pVQZ basis sets as well as their complete basis set limit. Our findings highlight the accuracy, flexibility, and scalability of ph-AFQMC with simple trial wavefunctions.
The phaseless auxiliary-field quantum Monte Carlo (ph-AFQMC) method is a stochastic imaginary-time projection method for computing ground-state properties of strongly correlated quantum systems, whose accuracy depends critically on the trial wavefunction. We present a benchmarking study of trial wavefunctions prepared using parameterized quantum circuits, spanning unitary coupled cluster (UCC), Hamiltonian-informed, Jastrow-inspired, and adaptively constructed ansatzes, and evaluate their accuracy, expressibility, and resource requirements within the quantum-computing-enhanced AFQMC (QC-AFQMC) framework. Trial states are prepared with GPU-accelerated statevector simulations using CUDA-Q and tested on stretched linear hydrogen chains and nine closed-shell molecules at equilibrium with varying static correlation. We find that at similar parameter counts, ph-AFQMC largely removes differences between families for weakly and moderately correlated systems: every trial yields a ph-AFQMC energy within chemical accuracy despite trial energy errors spanning nearly three orders of magnitude, whereas ansatz choice becomes consequential as static correlation grows. Our results indicate that trial variational energy is not always a reliable indicator of ph-AFQMC quality, and that adding variational parameters can lower the trial energy while degrading the ph-AFQMC energy. Across our benchmark systems, trial fidelity ranks trials more reliably than variational energy, and for spin-contaminated trials, the singlet-projected fidelity is more reliable than unprojected fidelity. The UCCSD ansatz initialized from unrestricted CCSD amplitudes provides a robust and inexpensive-to-optimize baseline across benchmark systems, with circuit depth as its principal limitation, while adaptive and compact layered ansatzes remain competitive in accuracy under strong correlation with substantially shallower circuits.