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

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.5479Physical Review LettersJan 26, 201294 citations

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

P. Fanto

Published Jun 28, 2017• 0 citations

OpenAlex

Abstract

Related papers

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.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.09373Physical Review CNov 1, 20177 citations

Particle-number projection in the finite-temperature mean-field approximation

P. Fanto, Y. Alhassid, G. F. Bertsch

Finite-temperature mean-field theories, such as the Hartree--Fock (HF) and Hartree--Fock--Bogoliubov (HFB) theories, are formulated in the grand-canonical ensemble, and their applications to the calculation of statistical properties of nuclei such as level densities require a reduction to the canonical ensemble. In a previous publication [Y. Alhassid et al., Phys. Rev. C 93, 044320 (2016)], it was found that ensemble-reduction methods based on the saddle-point approximation are not reliable in cases in which rotational symmetry or particle-number conservation is broken. In particular, the calculated HFB canonical entropy can be unphysical as a result of the inherent violation of particle-number conservation. In this work, we derive a general formula for exact particle-number projection after variation in the HFB approximation, assuming that the HFB Hamiltonian preserves time-reversal symmetry. This formula reduces to simpler known expressions in the HF and Bardeen--Cooper--Schrieffer (BCS) limits of the HFB. We apply this formula to calculate the thermodynamic quantities needed for level densities in the heavy nuclei $^{162}\mathrm{Dy}, ^{148}\mathrm{Sm}$ , and $^{150}\mathrm{Sm}$ . We find that the exact particle-number projection gives better physical results and is significantly more computationally efficient than the saddle-point methods. However, the fundamental limitations caused by broken symmetries in the mean-field approximation are still present.Finite-temperature mean-field theories, such as the Hartree--Fock (HF) and Hartree--Fock--Bogoliubov (HFB) theories, are formulated in the grand-canonical ensemble, and their applications to the calculation of statistical properties of nuclei such as level densities require a reduction to the canonical ensemble. In a previous publication [Y. Alhassid et al., Phys. Rev. C 93, 044320 (2016)], it was found that ensemble-reduction methods based on the saddle-point approximation are not reliable in cases in which rotational symmetry or particle-number conservation is broken. In particular, the calculated HFB canonical entropy can be unphysical as a result of the inherent violation of particle-number conservation. In this work, we derive a general formula for exact particle-number projection after variation in the HFB approximation, assuming that the HFB Hamiltonian preserves time-reversal symmetry. This formula reduces to simpler known expressions in the HF and Bardeen--Cooper--Schrieffer (BCS) limits of the HFB. We apply this formula to calculate the thermodynamic quantities needed for level densities in the heavy nuclei $^{162}\mathrm{Dy}, ^{148}\mathrm{Sm}$, and $^{150}\mathrm{Sm}$. We find that the exact particle-number projection gives better physical results and is significantly more computationally efficient than the saddle-point methods. However, the fundamental limitations caused by broken symmetries in the mean-field approximation are still present.

arXiv:1610.08954Physical Review CJul 7, 201718 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.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, 20130 citations

Pairing in Finite Systems: Beyond the HFB Theory

L. M. Robledo, G. F. Bertsch

The Hartree-Fock-Bogoliubov approximation is very useful for treating both long- and short-range correlations in finite quantum fermion systems, but it must be extended in order to describe detailed spectroscopic properties. One problem is the symmetry-breaking character of the HFB approximation. We present a general and systematic way to restore symmetries and to extend the configuration space using pfaffian formulas for the many-body matrix elements. The advantage of those formulas is that the sign of the matrix elements is unambiguously determined. It is also helpful to extend the space of configurations by constraining the HFB solutions in some way. A powerful method for finding these constrained solutions is the gradient method, based on the generalized Thouless transformation. The gradient method also preserves the number parity of the Bogoliubov transformation, which facilitates the application of the theory to systems with odd particle number.The Hartree-Fock-Bogoliubov approximation is very useful for treating both long- and short-range correlations in finite quantum fermion systems, but it must be extended in order to describe detailed spectroscopic properties. One problem is the symmetry-breaking character of the HFB approximation. We present a general and systematic way to restore symmetries and to extend the configuration space using pfaffian formulas for the many-body matrix elements. The advantage of those formulas is that the sign of the matrix elements is unambiguously determined. It is also helpful to extend the space of configurations by constraining the HFB solutions in some way. A powerful method for finding these constrained solutions is the gradient method, based on the generalized Thouless transformation. The gradient method also preserves the number parity of the Bogoliubov transformation, which facilitates the application of the theory to systems with odd particle number.

arXiv:1205.4443WORLD SCIENTIFIC eBooksMar 1, 20133 citations

Symmetry restoration of HFB states using projection

Carl Frostenson

The restoration of particle number and angular momentum symmetry using projection operators has been theoretically investigated for mean-field Hartree-Fock-Bogoliubov (HFB) states. A computer code for the projection of particle number was then implemented. To do so efficiently and avoid the sign ambiguity of the Onishi formula, the computation of overlap between quasi-particle vaccums using Pfaffians was investigated. Different algorithms for the the Pfaffian from~\cite{GBBertschRobledo},\cite{Wimmer}, was tested for performance and accuracy. For the same reason different truncations of the model space was also investigated. The particle number projector was then implemented in \textit{HOSPHE} \cite{Carlsson} for calculations of ground states pertaining to an effective Quadrupole plus Pairing Hamiltonian calculated using the SLy4 parametrization of the Skyrme interaction.The restoration of particle number and angular momentum symmetry using projection operators has been theoretically investigated for mean-field Hartree-Fock-Bogoliubov (HFB) states. A computer code for the projection of particle number was then implemented. To do so efficiently and avoid the sign ambiguity of the Onishi formula, the computation of overlap between quasi-particle vaccums using Pfaffians was investigated. Different algorithms for the the Pfaffian from~\cite{GBBertschRobledo},\cite{Wimmer}, was tested for performance and accuracy. For the same reason different truncations of the model space was also investigated. The particle number projector was then implemented in \textit{HOSPHE} \cite{Carlsson} for calculations of ground states pertaining to an effective Quadrupole plus Pairing Hamiltonian calculated using the SLy4 parametrization of the Skyrme interaction.

Jan 1, 20190 citations

Symmetry Restoration in Hartree-Fock-Bogoliubov Based Theories

G. F. Bertsch, Luis Miguel Robledo

arXiv:1108.5479Physical Review LettersJan 26, 201294 citations

^{162}\mathrm{Dy}, ^{148}\mathrm{Sm}

, and

^{150}\mathrm{Sm}

. We find that the exact particle-number projection gives better physical results and is significantly more computationally efficient than the saddle-point methods. However, the fundamental limitations caused by broken symmetries in the mean-field approximation are still present.

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.