© 2026 Aqora Quantum S.A.S.

Terms Privacy Legal Notice

arXiv:2304.13780 · The Journal of Chemical Physics

Robust formulation of Wick’s theorem for computing matrix elements between Hartree–Fock–Bogoliubov wavefunctions

Guo P. Chen, Gustavo E. Scuseria

Published Jun 15, 2023• 18 citations

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 Communication, 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.

DOI: 10.1063/5.0156124

Related papers

Norm-overlap formula for Hartree–Fock–Bogoliubov states with odd number parity

Makito Oi, T. Mizusaki

A formula to calculate the norm overlap between Hartree–Fock–Bogoliubov (HFB) states with odd number parity (one quasi-particle excited states) is derived with the help of Grassmann numbers and Fermion coherent states. The final formula is expressed in terms of a product of the Pfaffian for a neighboring even–even system (the zero quasi-particle state), and an extra factor consisting of the Bogoliubov transformation matrix and the anti-symmetric matrix in Thoulessʼ HFB ansatz for the even–even system.

arXiv:1110.2340•Physics Letters B•Dec 21, 2011•20 citations

Extension of Wick’s theorem for many-particle matrix elements

K. Tanabe, K. Enami, Naotaka Yoshinaga

Wick's theorem for the expectation value of an operator between many-particle states is generalized to the case of a homogeneous linear transformation operator $\^O$, which leaves invariant the $2N$-dimensional linear space composed of the annihilation and the creation operators for the N fermion or boson single-particle basis. Such a transformation generally belongs to the group $\mathrm{SO}(2N,C)$ for the case of fermion, or the group $\mathrm{Sp}(2N,C)$ for the case of boson. This theorem is applied to derive a reduction theorem for the matrix element of the type $〈\ensuremath{\alpha}|{\ensuremath{\alpha}}_{1}{\ensuremath{\alpha}}_{2}\ensuremath{\cdot}\ensuremath{\cdot}\ensuremath{\cdot}{\ensuremath{\alpha}}_{m}\^O{\ensuremath{\beta}}_{m+1}^{\ifmmode\dagger\else\textdagger\fi{}}{\ensuremath{\beta}}_{m+2}^{\ifmmode\dagger\else\textdagger\fi{}}\ensuremath{\cdot}\ensuremath{\cdot}\ensuremath{\cdot}{\ensuremath{\beta}}_{2n}^{\ifmmode\dagger\else\textdagger\fi{}}|\ensuremath{\beta}〉,$ whose many-quasiparticle states and vacua appearing in both sides are represented in different quasiparticle pictures as indicated by $\ensuremath{\alpha}$ and $\ensuremath{\beta}.$ The reduction formula is useful in calculating any nuclear transition matrix element for the quantum-number projection calculation associated with the constrained Hartree-Fock-Bogoliubov solutions, and also the generator coordinate method with quantum number projection.

Physical Review C•May 1, 1999•14 citations

Symmetry restoration in the finite-temperature Hartree-Fock-Bogoliubov approximation

P. Fanto

The finite-temperature Hartree-Fock-Bogoliubov (HFB) approximation often breaks symmetries of the underlying many-body Hamiltonian. Restricting the calculation of the HFB partition function to a subspace with good quantum numbers through projection after variation restores some of the correlations lost in breaking these symmetries, although effects of the broken symmetries such as sharp kinks at phase transitions remain. However, the most general projection after variation formula in the finite-temperature HFB approximation is limited by a sign ambiguity. Here, we extend the Pfaffian formula for the many-body traces of HFB density operators introduced by L. M. Robledo in Ref. [1] to eliminate this sign ambiguity and evaluate the more complicated many-body traces required in projection after variation in the most general HFB case. We validate our method through a proof-of-principle calculation of the particle-number-projected HFB thermal energy in a simple model.

•Jun 28, 2017•0 citations

Overlap formula between the HFB vacua for the symmetry restoration

Zao-Chun Gao, Qing-Li Hu, Y. S. Chen

Using the Pfaffian identity, the formula for the overlap between Hartree-Fock-Bogoliubov (HFB) vacua, recently proposed by Bertsch and Robledo [Phys. Rev. Lett. {\bf 108}, 042505(2012)] is deduced into a simpler form. This new formula does not explicitly depend on the Bloch-Messiah decomposition of the HFB transformation, but it requires the existence of the $V^{-1}$ matrix. It presents the significant convenience of application for the symmetry restoration in the wave functions of general Bogoliubov form. Numerical calculations show that this formula is able to reproduce the results of the Onishi formula with the same accuracy as well as having the correct sign. It is also shown that a small number of the integrating mesh points ($N_\text{mesh}=16$) is enough to perform the particle number projection even for the heavy nuclear system.

•Jun 13, 2013•0 citations

New and Practical Formulation for Overlaps of Bogoliubov Vacua

B. GILLIS CARLSSON, J. Rotureau

In this Letter, we present a new expression for the overlaps of wave functions in Hartree-Fock-Bogoliubov based theories. Starting from the Pfaffian formula by Bertsch et al. [1], an exact and computationally stable formula for overlaps is derived. We illustrate the convenience of this new formulation with a numerical application in the context of the particle-number projection method. This new formula allows for substantially increased precision and versatility in chemical, atomic, and nuclear physics applications, particularly for methods dealing with superfluidity, symmetry restoration, and uses of nonorthogonal many-body basis states.

arXiv:2010.08459•Physical Review Letters•Apr 26, 2021•25 citations

Projection after variation in the finite-temperature Hartree-Fock-Bogoliubov approximation

P. Fanto

The finite-temperature Hartree-Fock-Bogoliubov (HFB) approximation often breaks symmetries of the underlying many-body Hamiltonian. Restricting the calculation of the HFB partition function to a subspace with good quantum numbers through projection after variation restores some of the correlations lost in breaking these symmetries, although effects of the broken symmetries such as sharp kinks at phase transitions remain. However, the most general projection after variation formula in the finite-temperature HFB approximation is limited by a sign ambiguity. Here, I extend the Pfaffian formula for the many-body traces of HFB density operators introduced by Robledo [L. M. Robledo, Phys. Rev. C. 79, 021302(R) (2009)] to eliminate this sign ambiguity and evaluate the more complicated many-body traces required in projection after variation in the most general HFB case. The method is validated through a proof-of-principle calculation of the particle-number-projected HFB thermal energy in a simple model.

arXiv:1706.09373•Physical Review C•Nov 1, 2017•7 citations

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 Communication, 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.

A formula to calculate the norm overlap between Hartree–Fock–Bogoliubov (HFB) states with odd number parity (one quasi-particle excited states) is derived with the help of Grassmann numbers and Fermion coherent states. The final formula is expressed in terms of a product of the Pfaffian for a neighboring even–even system (the zero quasi-particle state), and an extra factor consisting of the Bogoliubov transformation matrix and the anti-symmetric matrix in Thoulessʼ HFB ansatz for the even–even system.

Wick's theorem for the expectation value of an operator between many-particle states is generalized to the case of a homogeneous linear transformation operator

\^O

, which leaves invariant the

2N

-dimensional linear space composed of the annihilation and the creation operators for the N fermion or boson single-particle basis. Such a transformation generally belongs to the group

\mathrm{SO}(2N,C)

for the case of fermion, or the group

\mathrm{Sp}(2N,C)

for the case of boson. This theorem is applied to derive a reduction theorem for the matrix element of the type $〈\ensuremath{\alpha}|{\ensuremath{\alpha}}_{1}{\ensuremath{\alpha}}_{2}\ensuremath{\cdot}\ensuremath{\cdot}\ensuremath{\cdot}{\ensuremath{\alpha}}_{m}\^O{\ensuremath{\beta}}_{m+1}^{\ifmmode\dagger\else\textdagger\fi{}}{\ensuremath{\beta}}_{m+2}^{\ifmmode\dagger\else\textdagger\fi{}}\ensuremath{\cdot}\ensuremath{\cdot}\ensuremath{\cdot}{\ensuremath{\beta}}_{2n}^{\ifmmode\dagger\else\textdagger\fi{}}|\ensuremath{\beta}〉,$ whose many-quasiparticle states and vacua appearing in both sides are represented in different quasiparticle pictures as indicated by $\ensuremath{\alpha}$ and $\ensuremath{\beta}.$ The reduction formula is useful in calculating any nuclear transition matrix element for the quantum-number projection calculation associated with the constrained Hartree-Fock-Bogoliubov solutions, and also the generator coordinate method with quantum number projection.

The finite-temperature Hartree-Fock-Bogoliubov (HFB) approximation often breaks symmetries of the underlying many-body Hamiltonian. Restricting the calculation of the HFB partition function to a subspace with good quantum numbers through projection after variation restores some of the correlations lost in breaking these symmetries, although effects of the broken symmetries such as sharp kinks at phase transitions remain. However, the most general projection after variation formula in the finite-temperature HFB approximation is limited by a sign ambiguity. Here, we extend the Pfaffian formula for the many-body traces of HFB density operators introduced by L. M. Robledo in Ref. [1] to eliminate this sign ambiguity and evaluate the more complicated many-body traces required in projection after variation in the most general HFB case. We validate our method through a proof-of-principle calculation of the particle-number-projected HFB thermal energy in a simple model.

Using the Pfaffian identity, the formula for the overlap between Hartree-Fock-Bogoliubov (HFB) vacua, recently proposed by Bertsch and Robledo [Phys. Rev. Lett. {\bf 108}, 042505(2012)] is deduced into a simpler form. This new formula does not explicitly depend on the Bloch-Messiah decomposition of the HFB transformation, but it requires the existence of the

V^{-1}

matrix. It presents the significant convenience of application for the symmetry restoration in the wave functions of general Bogoliubov form. Numerical calculations show that this formula is able to reproduce the results of the Onishi formula with the same accuracy as well as having the correct sign. It is also shown that a small number of the integrating mesh points (

N_\text{mesh}=16

) is enough to perform the particle number projection even for the heavy nuclear system.

In this Letter, we present a new expression for the overlaps of wave functions in Hartree-Fock-Bogoliubov based theories. Starting from the Pfaffian formula by Bertsch et al. [1], an exact and computationally stable formula for overlaps is derived. We illustrate the convenience of this new formulation with a numerical application in the context of the particle-number projection method. This new formula allows for substantially increased precision and versatility in chemical, atomic, and nuclear physics applications, particularly for methods dealing with superfluidity, symmetry restoration, and uses of nonorthogonal many-body basis states.

The finite-temperature Hartree-Fock-Bogoliubov (HFB) approximation often breaks symmetries of the underlying many-body Hamiltonian. Restricting the calculation of the HFB partition function to a subspace with good quantum numbers through projection after variation restores some of the correlations lost in breaking these symmetries, although effects of the broken symmetries such as sharp kinks at phase transitions remain. However, the most general projection after variation formula in the finite-temperature HFB approximation is limited by a sign ambiguity. Here, I extend the Pfaffian formula for the many-body traces of HFB density operators introduced by Robledo [L. M. Robledo, Phys. Rev. C. 79, 021302(R) (2009)] to eliminate this sign ambiguity and evaluate the more complicated many-body traces required in projection after variation in the most general HFB case. The method is validated through a proof-of-principle calculation of the particle-number-projected HFB thermal energy in a simple model.