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.

Factorization Of Commutators - The Wick Theorem For Coupled Operators

Commutation relations for functions of canonical conjugate operators

Another Formulation of the Wick's Theorem. Farewell, Pairing?

Product Formulas for Exponentials of Commutators

Cartesian operator factorization method for Hydrogen

A Q-operator for open spin chains II: boundary factorization

Post-Wick theorems for symbolic manipulation of second-quantized expressions in atomic many-body perturbation theory

Operator relationship between conventional coupled cluster and unitary coupled cluster

Factorization of the transition matrix for the general Jacobi system

Commuting operations factorise

Generalized Wick Decompositions

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

Asymmetric factorization method on supersymmetry: Complex operators

Commutators Estimates on Triebel–Lizorkin Spaces

Wigner Functions for Noncommutative Quantum Mechanics: a group representation based construction

Factorization of the Indefinite Convection-Diffusion Operator

Coarse indices of twisted operators

Infinitesimal Operators and the Distribution of Anticommutators and Commutators

Quotients of graph operators by symmetry representations

A Canonical Decomposition of Complex Symmetric Operators