Many-body computations by stochastic sampling in Hartree-Fock-Bogoliubov space

Abstract

We describe the computational ingredients for an approach to treat interacting fermion systems in the presence of 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 (QMC) techniques. The properties of HFB states, written in the form of either product states or Thouless states, are discussed. The states are shown to be propagated by generalized one-body operators. The states can be stabilized for numerical iteration, and overlaps between two such states and one-body Green's functions 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.

11 pages, 5 figures

References in corpus (3)

• Variational Monte Carlo Method Combined with Quantum-Number Projection and Multi-Variable Optimization
• The sign of the overlap of HFB wave functions
• Quasiparticle Coupled Cluster Theory for Pairing Interactions

Cited by in corpus (12)

• Coexistence of superconductivity with partially filled stripes in the Hubbard model
• Constrained-Path Auxiliary-Field Quantum Monte Carlo for Coupled Electrons and Phonons
• Auxiliary Field Quantum Monte Carlo for Multiband Hubbard Models: Controlling the Sign and Phase Problems to Capture Hund's Physics
• Many-body Majorana braiding without an exponential Hilbert space
• Robust formulation of Wick's theorem for computing matrix elements between Hartree-Fock-Bogoliubov wavefunctions
• Interfacing branching random walks with Metropolis sampling: constraint release in auxiliary-field quantum Monte Carlo
• Ising phase transitions and thermodynamics of correlated fermions in a two-dimensional spin-dependent lattice potential
• Charge Density Waves in a Quantum Plasma
• Exchangeless braiding of Majorana zero modes in weakly coupled Kitaev chains
• Accelerating ground-state auxiliary-field quantum Monte Carlo simulations by delayed update and block force-bias update
• NNQS-AFQMC: Neural network quantum states enhanced fermionic quantum Monte Carlo
• Numerical solution of large scale Hartree-Fock-Bogoliubov equations