Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Log inSign up
Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Blog
Job Board
Q3AS

© 2026 Aqora Quantum S.A.S.

TermsPrivacyLegal Notice

arXiv:2010.08459 · Physical Review Letters

New and Practical Formulation for Overlaps of Bogoliubov Vacua

B. GILLIS CARLSSON, J. Rotureau

Published Apr 26, 2021• 25 citations
arXivHTMLPDF

Abstract

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.

DOI
10.1103/physrevlett.126.172501

Related papers

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

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

Luis Miguel Robledo

Several technical aspects concerning the evaluation of the overlap between two mean-field wave functions of the Hartree-Fock-Bogoliubov type are discussed. The limit when several orbitals become fully occupied is derived as well as the formula to reduce the dimensionality of the problem when exactly empty orbitals are present. The formalism is also extended to deal with different bases for each of the wave functions. Several practical results concerning the evaluation of pfaffians, as well as the canonical decomposition of norm overlaps, are also discussed in the Appendices.

arXiv:1102.2151•Physical Review C•Jul 8, 2011•27 citations

Norm overlap between many-body states: Uncorrelated overlap between arbitrary Bogoliubov product states

B. Bally, T. Duguet

Background: State-of-the-art multi-reference energy density functional calculations require the computation of norm overlaps between different Bogoliubov quasiparticle many-body states. It is only recently that the efficient and unambiguous calculation of such norm kernels has become available under the form of Pfaffians [L. M. Robledo, Phys. Rev. C 79, 021302 (2009)]. Recently developed particle-number-restored Bogoliubov coupled-cluster (PNR-BCC) and particle-number-restored Bogoliubov many-body perturbation (PNR-BMBPT) ab initio theories [T. Duguet and A. Signoracci, J. Phys. G 44, 015103 (2017)] make use of generalized norm kernels incorporating explicit many-body correlations. In PNR-BCC and PNR-BMBPT, the Bogoliubov states involved in the norm kernels differ specifically via a global gauge rotation.Purpose: The goal of this work is threefold. We wish (i) to propose and implement an alternative to the Pfaffian method to compute unambiguously the norm overlap between arbitrary Bogoliubov quasiparticle states, (ii) to extend the first point to explicitly correlated norm kernels, and (iii) to scrutinize the analytical content of the correlated norm kernels employed in PNR-BMBPT. Point (i) constitutes the purpose of the present paper while points (ii) and (iii) are addressed in a forthcoming paper.Methods: We generalize the method used in another work [T. Duguet and A. Signoracci, J. Phys. G 44, 015103 (2017)] in such a way that it is applicable to kernels involving arbitrary pairs of Bogoliubov states. The formalism is presently explicated in detail in the case of the uncorrelated overlap between arbitrary Bogoliubov states. The power of the method is numerically illustrated and benchmarked against known results on the basis of toy models of increasing complexity.Results: The norm overlap between arbitrary Bogoliubov product states is obtained under a closed-form expression allowing its computation without any phase ambiguity. The formula is physically intuitive, accurate, and versatile. It equally applies to norm overlaps between Bogoliubov states of even or odd number parity. Numerical applications illustrate these features and provide a transparent representation of the content of the norm overlaps.Conclusions: The complex norm overlap between arbitrary Bogoliubov states is computed, without any phase ambiguity, via elementary linear algebra operations. The method can be used in any configuration mixing of orthogonal and non-orthogonal product states. Furthermore, the closed-form expression extends naturally to correlated overlaps at play in PNR-BCC and PNR-BMBPT. As such, the straight overlap between Bogoliubov states is the zero-order reduction of more involved norm kernels to be studied in a forthcoming paper.

arXiv:1704.05324•Physical Review C•Feb 5, 2018•23 citations

Symmetry Restoration in Hartree-Fock-Bogoliubov Based Theories

G. F. Bertsch, Luis Miguel Robledo

We present a Pfaffian formula for projection and symmetry restoration for wave functions of the general Bogoliubov form, including quasiparticle excited states and linear combinations of them. This solves a long-standing problem in calculating states of good symmetry, arising from the sign ambiguity of the commonly used determinant formula. A simple example is given of projecting a good particle number and angular momentum from a Bogoliubov wave function in the Fock space of a single j-shell.

arXiv:1108.5479•Physical Review Letters•Jan 26, 2012•94 citations

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

Guo P. Chen, Gustavo E. Scuseria

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.

arXiv:2304.13780•The Journal of Chemical Physics•Jun 15, 2023•18 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
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.
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−1V^{-1}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 (Nmesh=16N_\text{mesh}=16Nmesh​=16) is enough to perform the particle number projection even for the heavy nuclear system.
Several technical aspects concerning the evaluation of the overlap between two mean-field wave functions of the Hartree-Fock-Bogoliubov type are discussed. The limit when several orbitals become fully occupied is derived as well as the formula to reduce the dimensionality of the problem when exactly empty orbitals are present. The formalism is also extended to deal with different bases for each of the wave functions. Several practical results concerning the evaluation of pfaffians, as well as the canonical decomposition of norm overlaps, are also discussed in the Appendices.
Background: State-of-the-art multi-reference energy density functional calculations require the computation of norm overlaps between different Bogoliubov quasiparticle many-body states. It is only recently that the efficient and unambiguous calculation of such norm kernels has become available under the form of Pfaffians [L. M. Robledo, Phys. Rev. C 79, 021302 (2009)]. Recently developed particle-number-restored Bogoliubov coupled-cluster (PNR-BCC) and particle-number-restored Bogoliubov many-body perturbation (PNR-BMBPT) ab initio theories [T. Duguet and A. Signoracci, J. Phys. G 44, 015103 (2017)] make use of generalized norm kernels incorporating explicit many-body correlations. In PNR-BCC and PNR-BMBPT, the Bogoliubov states involved in the norm kernels differ specifically via a global gauge rotation.Purpose: The goal of this work is threefold. We wish (i) to propose and implement an alternative to the Pfaffian method to compute unambiguously the norm overlap between arbitrary Bogoliubov quasiparticle states, (ii) to extend the first point to explicitly correlated norm kernels, and (iii) to scrutinize the analytical content of the correlated norm kernels employed in PNR-BMBPT. Point (i) constitutes the purpose of the present paper while points (ii) and (iii) are addressed in a forthcoming paper.Methods: We generalize the method used in another work [T. Duguet and A. Signoracci, J. Phys. G 44, 015103 (2017)] in such a way that it is applicable to kernels involving arbitrary pairs of Bogoliubov states. The formalism is presently explicated in detail in the case of the uncorrelated overlap between arbitrary Bogoliubov states. The power of the method is numerically illustrated and benchmarked against known results on the basis of toy models of increasing complexity.Results: The norm overlap between arbitrary Bogoliubov product states is obtained under a closed-form expression allowing its computation without any phase ambiguity. The formula is physically intuitive, accurate, and versatile. It equally applies to norm overlaps between Bogoliubov states of even or odd number parity. Numerical applications illustrate these features and provide a transparent representation of the content of the norm overlaps.Conclusions: The complex norm overlap between arbitrary Bogoliubov states is computed, without any phase ambiguity, via elementary linear algebra operations. The method can be used in any configuration mixing of orthogonal and non-orthogonal product states. Furthermore, the closed-form expression extends naturally to correlated overlaps at play in PNR-BCC and PNR-BMBPT. As such, the straight overlap between Bogoliubov states is the zero-order reduction of more involved norm kernels to be studied in a forthcoming paper.
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.
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.