Background: Recently, variational Hartree-Fock-Bogoliubov (HFB) mean-field equations were shown to possess a mathematically well-defined zero-pairing limit, independently of the closed- or open-shell character of the system under consideration. This limit is nontrivial for open-shell systems such that HFB theory does not reduce to the Hartree-Fock (HF) formalism in all cases.Purpose: The present work extends this analysis to finite-temperature HFB (FTHFB) theory by investigating the behavior of this more general formalism in the combined zero-temperature and zero-pairing limits.Methods: The zero-pairing and zero-temperature limits of the FTHFB statistical density operator constrained to carry an arbitrary (integer) number of particles A on average is worked out analytically and realized numerically using a two-nucleon interaction.Results: While the FTHFB density operator reduces to the projector corresponding to a pure HF Slater determinant for closed-shell nuclei, the FTHFB formalism does not reduce to the HF theory in all cases in the zero-temperature and zero-pairing limits, i.e., for open-shell nuclei. However, the fact that a nucleus can be of open-shell character in these joint limits is necessarily the result of some symmetry restrictions. Whenever it is the case, the nontrivial description obtained for open-shell systems is shown to depend on the order with which both limits are taken, i.e., the two limits do not commute for these systems. When the zero-temperature limit is performed first, the FTHFB density operator is demoted to a projector corresponding to a pure state made out of a linear combination of a finite number of Slater determinants with different (even) numbers of particles. When the zero-pairing limit is performed first, the FTHFB density operator remains a statistical mixture of a finite number of Slater determinants with both even and odd particle numbers. While the entropy (pairing density) is zero in the first (second) case, it does not vanish in the second (first) case in spite of the temperature (pairing) tending towards zero. The difference between both limits can have striking consequences for the (thermal) expectation values of observables. For instance, the particle-number variance does not vanish in either case and has limiting values that differ by a factor of two in both cases.Conclusions: While in the textbook situation associated with closed-shell nuclei Hartree-Fock-Bogoliubov (finite-temperature Hartree-Fock) theory reduces to Hartree-Fock theory in the zero-pairing (zero-temperature) limit, the present analysis demonstrates that a nontrivial and unexpected limit is obtained for this formalism in open-shell systems. This result sheds a new light on certain aspects of this otherwise very well-studied many-body formalism.
Background: The variational Hartree-Fock-Bogoliubov (HFB) mean-field theory is the starting point of various (ab initio) many-body methods dedicated to superfluid systems. In this context, pairing correlations may be driven towards zero either on purpose via HFB calculations constrained on, e.g., the particle-number variance or simply because internucleon interactions cannot sustain pairing correlations in the first place in, e.g., closed-shell systems. While taking this limit constitutes a text-book problem when the system is of closed-shell character, it is typically, although wrongly, thought to be ill-defined whenever the naive filling of single-particle levels corresponds to an open-shell system.Purpose: The present work demonstrates that the zero-pairing limit of an HFB state is well-defined independently of the average particle number
A it is constrained to. Still, the nature of the limit state is shown to depend on the regime, i.e., on whether the nucleus characterizes as a closed-shell or an open-shell system when taking the limit. Finally, the consequences of the zero-pairing limit on Bogoliubov many-body perturbation theory (BMBPT) calculations performed on top of the HFB reference state are illustrated.Methods: The zero-pairing limit of a HFB state constrained to carry an arbitrary (integer) number of particles
A on average is worked out analytically and realized numerically using a two-nucleon interaction derived within the frame of chiral effective field theory.Results: The zero-pairing limit of the HFB state is mathematically well-defined, independently of the closed- or open-shell character of the system in the limit. Still, the nature of the limit state strongly depends on the underlying shell structure and on the associated naive filling reached in the zero-pairing limit for the particle number
A of interest. First, the textbook situation is recovered for closed-shell systems, i.e., the limit state is reached for a finite value of the pairing strength (the well-known BardeenCooperSchrieffer (BCS) collapse) and takes the form of a single Slater determinant displaying (i) zero pairing energy, (ii) nondegenerate elementary excitations, and (iii) zero particle-number variance. Contrarily, a nonstandard situation is obtained for open-shell systems, i.e., the limit state is only reached for a zero value of the pairing interaction (no BCS collapse) and takes the form of a specific finite linear combination of Slater determinants displaying (a) a nonzero pairing energy, (b) degenerate elementary excitations, and (c) a nonzero particle-number variance for which an analytical formula is derived. This nonzero particle-number variance acts as a lower bound that depends in a specific way on the number of valence nucleons and on the degeneracy of the valence shell. All these findings are confirmed and illustrated numerically. Last but not least, BMBPT calculations of closed-shell (open-shell) nuclei are shown to be well-defined (ill-defined) in the zero-pairing limit.Conclusions: While HFB theory has been intensively scrutinized formally and numerically over the last decades, it still uncovers unknown and somewhat unexpected features. In the present paper, the zero-pairing limit of a HFB state carrying an arbitrary number of particles has been worked out and shown to lead to drastic differences and consequences depending on the closed-shell or open-shell nature of the system in that limit. From a general perspective, the present analysis demonstrates that HFB theory does not reduce to HF theory when the pairing field is driven to zero in the HFB Hamiltonian matrix.
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
162Dy,148Sm, and
150Sm. 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.
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, 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.
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.
This paper describes a general group theoretical analysis of the temperature–Hartree-Fock–Bogoliubov (HFB) equation and its solution. The action of the symmetry group G0 of the system on the HFB Hamiltonian, the HFB density matrix, and the HFB Green function are defined. It is shown that the HFB equation and its solution is classified by a subgroup G of G0, which is the invariance group of the HFB Hamiltonian, the HFB density matrix, and the HFB Green function corresponding the solution. General expression of the instability of a solution and its decomposition into R-rep (single-valued irreducible representation over the real number field) components of the invariance group of the solution are obtained. The self-consistent field (SCF) condition is decomposed into R-rep components of G0.