Abstract The so-called phaseless quantum Monte-Carlo method currently offers one of the best performing theoretical framework to investigate interacting Fermi systems. It allows to extract an approximate ground-state wavefunction by averaging independent-particle states undergoing a Brownian motion in imaginary-time. Here, we extend the approach to a random walk in the space of Hartree–Fock–Bogoliubov (HFB) vacua that are better suited for superconducting or superfluid systems. Well-controlled statistical errors are ensured by constraining stochastic paths with the help of a trial wavefunction. It also guides the dynamics and takes the form of a linear combination of HFB ansätze. Estimates for the observables are reconstructed through an extension of Wick’s theorem to matrix elements between HFB product states. The usual combinatory complexity associated to the application of this theorem for four- and more- body operators is bypassed with a compact expression in terms of Pfaffians. The limiting case of a stochastic motion within Slater determinants but guided with HFB trial wavefunctions is also considered. Finally, exploratory results for the spin-polarized Hubbard model in the attractive regime are presented.
Abstract The so-called phaseless quantum Monte-Carlo method currently offers one of the best performing theoretical framework to investigate interacting Fermi systems. It allows to extract an approximate ground-state wavefunction by averaging independent-particle states undergoing a Brownian motion in imaginary-time. Here, we extend the approach to a random walk in the space of Hartree–Fock–Bogoliubov (HFB) vacua that are better suited for superconducting or superfluid systems. Well-controlled statistical errors are ensured by constraining stochastic paths with the help of a trial wavefunction. It also guides the dynamics and takes the form of a linear combination of HFB ansätze. Estimates for the observables are reconstructed through an extension of Wick’s theorem to matrix elements between HFB product states. The usual combinatory complexity associated to the application of this theorem for four- and more- body operators is bypassed with a compact expression in terms of Pfaffians. The limiting case of a stochastic motion within Slater determinants but guided with HFB trial wavefunctions is also considered. Finally, exploratory results for the spin-polarized Hubbard model in the attractive regime are presented.
We describe the computational ingredients for an approach to treat interacting fermion systems with pairing fields, based on path integrals in the space of Hartree-Fock-Bogoliubov (HFB) wave functions. The path integrals can be evaluated by Monte Carlo, via random walks of HFB wave functions whose orbitals evolve stochastically. The approach combines the advantage of HFB theory in paired fermion systems and many-body quantum Monte Carlo techniques. The properties of HFB states, written in the form of either product states or Thouless states, are discussed. The states preserve forms when propagated by generalized one-body operators. They can be stabilized for numerical iteration. Overlaps and one-body Green's functions between two such states can be computed. A constrained-path or phaseless approximation can be applied to the random walks of the HFB states if a sign problem or phase problem is present. The method is illustrated with an exact numerical projection in the Kitaev model, and in the Hubbard model with attractive interaction under an external pairing field.
We formulate and implement the phaseless auxiliary-field quantum Monte Carlo (AFQMC) method in the Hartree-Fock-Bogoliubov (HFB) manifold for systems with strong pairing correlations. Our formulation introduces a more flexible representation of HFB walkers with an explicit analytical expression for the normalization factor derived from the representation matrix with respect to the reference vacuum. We also extend our previously developed stochastic gauge formalism from the Slater determinant manifold to the HFB framework, providing additional flexibility for controlling stochastic fluctuations. Benchmark calculations for the Richardson model demonstrate the accuracy and numerical stability of the method.
Since the form of PN is the same as that occurring in the statistical mechanics of the classical gas (replace f 2 ( r ij ) by exp[- V ( r ij /kT ]), we can use the same integration techniques that are used in the classical problem. W. L. McMillan, 1965 Summary Building on the random walk methods developed in the previous chapter, we show how to compute properties of many-body trial wavefunctions using a random walk. This method, called variational Monte Carlo, is the simplest stochastic quantum many-body technique. Whereas mean-field methods are usually limited to single determinants, variational Monte Carlo can treat any correlated trial function, as long as its values are computable. We discuss how to optimize such trial wavefunctions, how to compute their momentum distribution, how to use non-local pseudopotentials, how to compute excited states, and how to correct for the finite size of the simulation cell. Deterministic quantum methods have difficulties. For example, the Hartree–Fock method assumes the wavefunction is a single Slater determinant, neglecting correlation. If one expands as a sum of determinants, it is very difficult to have the results size-consistent since the number of determinants needed will grow exponentially with the system size. As we have seen, the DMFT method introduced in Ch. 16 assumes locality. In Ch. 6 we discussed general properties of many-body wavefunctions. Using Monte Carlo methods, we can directly incorporate correlation into a wavefunction, without having to make any further approximations other than the form of the correlation factors. In many cases the energy and other properties are very close to the exact results. Some of the usual restrictions on the form of the many-body wavefunction are not an issue in variational Monte Carlo. The most important generalization of the HF wavefunction is to put correlation directly into the wavefunction via the “Jastrow” factor. At next order, one can use the “backflow” wavefunction, in which correlation is also built into the determinant. The variational Monte Carlo method (VMC) was first used by McMillan [44] to calculate the ground-state properties of superfluid 4 He. One of the key problems at that time was whether the observed superfluid properties were a consequence of Bose condensation.
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 explore the possibility of implementing random walks in the manifold of Hartree-Fock-Bogoliubov wave functions. The goal is to extend state-of-the-art quantum Monte Carlo approaches, in particular the constrained-path auxiliary-field quantum Monte Carlo technique, to systems where finite pairing order parameters or complex pairing mechanisms, e.g., Fulde-Ferrell-Larkin-Ovchinnikov pairing or triplet pairing, may be expected. Leveraging the flexibility to define a vacuum state tailored to the physical problem, we discuss a method to use imaginary-time evolution of Hartree-Fock-Bogoliubov states to compute ground state correlations, extending beyond situations spanned by current formalisms. Illustrative examples are provided.
Quantum many-body systems can be studied theoretically from first principles using microscopic models, but the computational cost of numerically exact calculations increases exponentially with the number of particles and the complexity of the interaction, limiting the systems that can be analysed. The Hartree-Fock-Bogoliubov (HFB) and quasiparticle random phase approximation (QRPA) methods, heavily used in nuclear physics, rely on mean-field approximations but scale more favourably with increasing system size and complexity. In this work we apply these methods to the problem of a two-dimensional harmonically trapped Fermi gas with attractive interactions, with the goal of investigating the validity of the approximations in the few-body limit. We implement the approximate methods in code and then compare them with exact calculations. We find that the HFB method can account for pairing correlations that arise due to the attractive interaction, predicting well the ground state energy even at high interaction strengths. Also, the excitation spectra predicted by the QRPA method qualitatively match the true spectra, reproducing the expected subset of excitations. Finally, we see that perturbative corrections can be used to improve the ground state energy estimate of the HFB.