The authors provide a state-of-the-art overview of projection techniques employed in beyond-mean-field many-body methods. They focus on two symmetries that are of critical importance for the description of nuclear spectroscopy, particle number, and angular momentum. The presentation covers formal aspects of symmetry restoration (group theory) and practical aspects related to numerical implementation. A number of illustrative examples are given, including cranked states in odd-$A$ nuclei.
The authors provide a state-of-the-art overview of projection techniques employed in beyond-mean-field many-body methods. They focus on two symmetries that are of critical importance for the description of nuclear spectroscopy, particle number, and angular momentum. The presentation covers formal aspects of symmetry restoration (group theory) and practical aspects related to numerical implementation. A number of illustrative examples are given, including cranked states in odd-A nuclei.
Abstract In many so-called ‘beyond-mean-field’ many-body methods, one creates symmetry-breaking states and then projects out states with good quantum number(s); the most important example is angular momentum. Motivated by the computational intensity of symmetry restoration, we investigate the numerical convergence of two competing methods for angular momentum projection with rotations over Euler angles, the textbook-standard projection through quadrature, and a recently introduced projection through linear algebra. We find well-defined patterns of convergence with increasing number of mesh points (for quadrature) and cut-offs (for linear algebra). Because the method of projection through linear algebra requires inverting matrices generated on a mesh of Euler angles, we discuss two methods for robustly reducing the number of required evaluations. Reviewing the literature, we find our inversion involving rotations about the z -axis is equivalent to trapezoidal ‘quadrature’ commonly used as well as Fomenko projection used for particle-number projection. The efficiency depends upon the number of angular momentum J to be projected, but in general inversion methods, including Fomenko projection/trapezoidal ‘quadrature’ dramatically improve the efficiency.
Given elementary inter-nucleon interactions,the resolution of the A-body Schrödinger equation gives access to the properties of the nuclear eigenstates. Many-body perturbation theory (MBPT) based on a Hartree-Fock mean field allows one to treat the dynamical correlations at play in doubly-closed-shell nuclei. Beyond shell closures,breaking U(1) symmetry associated with the conservation of neutron and proton numbers further permits to capture static correlations at the meanfield level and to formulate the well-behaved Bogoliubov many-body perturbation theory (BMBPT).Nevertheless, BMBPT results show contaminations associated with the broken symmetry that is only emergent on finite quantum system such as theatomic nucleus. Thus, the restoration of U(1) symmetry beyond the mean field is necessary for acorrect description and gives rise to the projected BMBPT method (PBMBPT). The goal is to implement PBMBPT in order to perform ab initio calculations of singly-open-shell mid-mass nuclei.The present work provides systematic solutions to the formal and technical problems arising during the implementation of PBMBPT: an automatic and safe method to generate PBMBPT diagrams and associated expressions, a framework to calculate other observables than the energy, an extension of the normal ordering approximation to symmetry breaking methods as well as an automatic and safe tool for spherical reduction of tensor networks.
The mean-field approximation based on effective interactions or density functionals plays a pivotal role in the description of finite quantum many-body systems that are too large to be treated by ab initio methods. Some examples are strongly interacting medium and heavy mass atomic nuclei and mesoscopic condensed matter systems. In this approach, the linear Schrodinger equation for the exact many-body wave function is mapped onto a non-linear density-dependent one-body potential problem. This approximation, not only provides computationally very simple solutions even for systems with many particles, but due to the non-linearity, it also allows for obtaining solutions that break essential symmetries of the system, often connected with phase transitions. In this way, additional correlations are subsumed in the system. However, the mean-field approach suffers from the drawback that the corresponding wave functions do not have sharp quantum numbers and, therefore, many results cannot be compared directly with experimental data. In this article, we discuss general group-theory techniques to restore the broken symmetries, and provide detailed expressions on the restoration of translational, rotational, spin, isospin, parity and gauge symmetries, where the latter corresponds to the restoration of the particle number. In order to avoid the numerical complexity of exact projection techniques, various approximation methods available in the literature are examined. Applications of the projection methods are presented for simple nuclear models, realistic calculations in relatively small configuration spaces, nuclear energy density functional theory, as well as in other mesoscopic systems. We also discuss applications of projection techniques to quantum statistics in order to treat the averaging over restricted ensembles with fixed quantum numbers. Further, unresolved problems in the application of the symmetry restoration methods to the energy density functional theories are highlighted in the present work.
Combining several techniques, we propose an efficient and numerically reliable method to perform the quantum number projection and configuration mixing for most general mean-field states, i.e., the Hartree-Fock-Bogoliubov (HFB) type product states without symmetry restrictions. As for the example of calculations, we show the results of the simultaneous parity, number and angular-momentum projection from HFB type states generated from the cranked Woods-Saxon mean-field with a very large basis that is composed of Nmax = 20 spherical harmonic oscillator shells.
Abstract The mean-field approximation based on effective interactions or density functionals plays a pivotal role in the description of finite quantum many-body systems that are too large to be treated by ab initio methods. Some examples are strongly interacting medium and heavy mass atomic nuclei and mesoscopic condensed matter systems. In this approach, the linear Schrödinger equation for the exact many-body wave function is mapped onto a non-linear one-body potential problem. This approximation, not only provides computationally very simple solutions even for systems with many particles, but due to the non-linearity, it also allows for obtaining solutions that break essential symmetries of the system, often connected with phase transitions. In this way, additional correlations are subsumed in the system. However, the mean-field approach suffers from the drawback that the corresponding wave functions do not have sharp quantum numbers and, therefore, many results cannot be compared directly with experimental data. In this article, we discuss general group-theory techniques to restore the broken symmetries, and provide detailed expressions on the restoration of translational, rotational, spin, isospin, parity and gauge symmetries, where the latter corresponds to the restoration of the particle number. In order to avoid the numerical complexity of exact projection techniques, various approximation methods available in the literature are examined. Applications of the projection methods are presented for simple nuclear models, realistic calculations in relatively small configuration spaces, nuclear energy density functional (EDF) theory, as well as in other mesoscopic systems. We also discuss applications of projection techniques to quantum statistics in order to treat the averaging over restricted ensembles with fixed quantum numbers. Further, unresolved problems in the application of the symmetry restoration methods to the EDF theories are highlighted in the present work.