arXiv:2409.11571 · The Journal of Chemical Physics

Monte Carlo methods in the manifold of Hartree–Fock–Bogoliubov wave functions

Ettore Vitali, Peter Rosenberg, Shiwei Zhang

Published Oct 7, 2024• 2 citations

Abstract

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.

DOI: 10.1063/5.0228450

Related papers

Phaseless auxiliary-field quantum Monte Carlo in the Hartree-Fock-Bogoliubov manifold

Zhaozhan Zhang

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.

arXiv:2607.21929Jul 24, 20260 citations

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

Hao Shi, Shiwei Zhang

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.

arXiv:1611.02689Physical review. B./Physical review. BJan 27, 201728 citations

Phaseless quantum Monte-Carlo approach to strongly correlated superconductors with stochastic Hartree–Fock–Bogoliubov wavefunctions

Olivier Juillet, Alexandre Leprévost, J. Bonnard +1 more

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.

arXiv:1610.08022Journal of Physics A Mathematical and TheoreticalFeb 24, 20174 citations

Auxiliary-Field Quantum Monte Carlo Calculations of Excited States and Strongly Correlated Systems

S. Zhang

to develop quantum Monte Carlo (QMC) approaches for excited states and for more systematic and accurate treatment of ground states of strongly correlated sytems and heavier elements. We use the auxiliary-field quantum Monte Carlo (AFQMC) framework, which is a many-body total energy method that looks like many coupled mean-field calculations. I will discuss this method’s relations and dierences with traditional density-functional and quantum chemistry methods, and with other forms of QMC. The AFQMC method uses importance-sampled random walks in the space of nonorthogonal Slater determinants to project out the many-body ground state, or excited state. The random walks are constrained by a trial wave function (TWF) to remove the so-called sign problem, giving approximate (but non-perturbative) results. With singledeterminant TWFs from Hartree-Fock or density-functional theory (DFT), this approach can be thought of as a post-DFT method with no parameter. In this form, the accuracy is typically comparable to the quantum chemistry coupled-cluster CCSD(T) method near equilibrium geometries, and better than CCSD(T) when bonds are stretched. The computational cost of AFQMC scales as N 3 with characteristic size of the system. Its accuracy can be further improved at the cost of increased computing (releasing the constraint) or better TWF. With AFQMC, we can systematically target selected orbitals to correlate. This allows many-body calculations in solids with a frozen core, with down-folded Hamiltonians, or possibly embeded in a DFT calculations to extend scales. Results from a few on-going applications will be presented, including band structures in simple solids, magnetic phases in MnO, and Cobalt adsoption on graphene.

Jan 1, 20120 citations

Multi-determinant generalized Hartree-Fock wave functions in Monte Carlo calculations

Chia‐Chen Chang, Miguel A. Morales

The quantum Monte Carlo algorithm is arguably one of the most powerful computational many-body methods, enabling accurate calculation of many properties in interacting quantum systems. In the presence of the so-called sign problem, the algorithm typically relies on trial wave functions to eliminate the exponential decay of signal-to-noise ratio, usually at the expense of a bias in the result. The quality of the trial state therefore is critical for accurate simulations. In this work, benchmark results of the ground state auxiliary-field quantum Monte Carlo method are reported for the Hubbard model on several geometries. We demonstrate that when multi-determinant generalized Hartree-Fock states are used as trial wave functions, the systematic errors can be systematically reduced to a low level and the results compare favorably with the exact diagonalization data.

arXiv:1711.02154Nov 6, 20174 citations

Computational Studies Of Strongly Correlated Quantum Matter

Hao Shi

The study of strongly correlated quantum many-body systems is an outstanding challenge. Highly accurate results are needed for the understanding of practical and fundamental problems in condensed-matter physics, high energy physics, material science, quantum chemistry and so on. Our familiar mean-field or perturbative methods tend to be ineffective. Numerical simulations provide a promising approach for studying such systems. The fundamental difficulty of numerical simulation is that the dimension of the Hilbert space needed to describe interacting systems increases exponentially with the system size. Quantum Monte Carlo (QMC) methods are one of the best approaches to tackle the problem of enormous Hilbert space. They have been highly successful for boson systems and unfrustrated spin models. For systems with fermions, the exchange symmetry in general causes the infamous sign problem, making the statistical noise in the computed results grow exponentially with the system size. This hinders our understanding of interesting physics such as high-temperature superconductivity, metal-insulator phase transition. In this thesis, we present a variety of new developments in the auxiliary-field quantum Monte Carlo (AFQMC) methods, including the incorporation of symmetry in both the trial wave function and the projector, developing the constraint release method, using the force-bias to drastically improve the efficiency in Metropolis framework, identifying and solving the infinite variance problem, and sampling Hartree-Fock-Bogoliubov wave function. With these developments, some of the most challenging many-electron problems are now under control. We obtain an exact numerical solution of two-dimensional strongly interacting Fermi atomic gas, determine the ground state properties of the 2D Fermi gas with Rashba spin-orbit coupling, provide benchmark results for the ground state of the two-dimensional Hubbard model, and establish that the Hubbard model has a stripe order in the underdoped region.

Jan 1, 20170 citations