Abstract:Numerical difficulties associated with computing matrix elements of operators between Hartree-Fock-Bogoliubov (HFB) wavefunctions have plagued the development of HFB-based many-body theories for decades. The problem arises from divisions by zero in the standard formulation of the nonorthogonal Wick's theorem in the limit of vanishing HFB overlap. In this paper, we present a robust formulation of Wick's theorem that stays well-behaved regardless of whether the HFB states are orthogonal or not. This new formulation ensures cancellation between the zeros of the overlap and the poles of the Pfaffian, which appears naturally in fermionic systems. Our formula explicitly eliminates self-interaction, which otherwise causes additional numerical challenges. A computationally efficient version of our formalism enables robust symmetry-projected HFB calculations with the same computational cost as mean-field theories. Moreover, we avoid potentially diverging normalization factors by introducing a robust normalization procedure. The resulting formalism treats even and odd number of particles on equal footing and reduces to Hartree-Fock as a natural limit. As proof of concept, we present a numerically stable and accurate solution to a Jordan-Wigner-transformed Hamiltonian, whose singularities motivated the present work. Our robust formulation of Wick's theorem is a most promising development for methods using quasiparticle vacuum states.

Generalization of Cramer's rule and its application to the projection of Hartree-Fock wave function

A Generalized Family Of Discrete Pt-Symmetric Square Wells

Symmetry restoration in Hartree-Fock-Bogoliubov based theories

Non-simultaneity in two-photon coincidence spectroscopy

Symmetry-Projected Nuclear-Electronic Hartree-Fock: Eliminating Rotational Energy Contamination

Generalized Hartree Fock Theory for Dispersion Relations of Interacting Fermionic Lattice Systems

Robust formulation of Wick's theorem for computing matrix elements between Hartree-Fock-Bogoliubov wavefunctions

Grassmann Integral and Balian–Brézin Decomposition in Hartree–Fock–Bogoliubov Matrix Elements

General technique for analytical derivatives of post-projected Hartree-Fock

Generalization of the Kohn-Sham system that can represent arbitary one electron density matrices

Technical aspects of the evaluation of the overlap of Hartree- Fock- Bogoliubov wave functions

A new generalization of Fielder's lemma with applications

Projections in Toeplitz algebra

Mehta's eigenvectors for the finite Hartely transform

Analytical Hartree-Fock gradients for periodic systems

Application of Structured Matrices for Solving Hartree-Fock Equations

A Simple One-Electron Expression for Electron Rotational Factors

Szegö projection and matrix Hilbert transform in Hermitean Clifford analysis

Symbolic, numeric and quantum computation of Hartree-Fock equation

Weyl Correspondence and P-Representation As Operator Fredholm Equations and Their Solutions

On A Generalization of Motohashi's Formula