---
title: Angular Momentum Projection in Nuclear Theory
url: https://www.emergentmind.com/topics/angular-momentum-projection
type: topic
---

# Angular Momentum Projection in Nuclear Theory

Angular momentum projection is the symmetry-restoration procedure that extracts states with definite total angular momentum from intrinsically deformed or otherwise symmetry-breaking many-body wave functions. In nuclear structure theory, intrinsic Slater determinants, Hartree–Fock or Hartree–Fock–Bogoliubov vacua, antisymmetrized molecular dynamics states, and cluster-model states are flexible descriptions of deformation, pairing, and clustering, but they are generally not eigenstates of \(\hat J^2\). Projection restores the laboratory-frame quantum numbers required by rotational invariance and is therefore a basic ingredient of projected mean-field theory, projected generator-coordinate methods, multiconfiguration mixing, and several algebraic alternatives to integral projection [1608.08726] [1307.6497] [1709.04028].

## 1. Standard operator framework

The standard projector onto total angular momentum \(J\), laboratory projection \(M\), and intrinsic projection \(K\) is
\[
\hat{P}_{MK}^J = \frac{2J + 1}{8\pi^2}\int D_{MK}^{J\ast} (\alpha, \beta, \gamma)\, \hat{R}(\alpha, \beta, \gamma)\, d\alpha\, d\beta\, d\gamma ,
\]
with
\[
\hat{R}(\alpha, \beta, \gamma) = e^{-i\alpha \hat{J}_z} e^{-i\beta \hat{J}_y} e^{-i\gamma \hat{J}_z},
\qquad
D^J_{MK}(\alpha,\beta,\gamma)=e^{-iM\alpha}d^J_{MK}(\beta)e^{-iK\gamma}.
\]
Acting on an intrinsic state \(|\Phi\rangle\), \(\hat P^J_{MK}|\Phi\rangle\) isolates the component carrying the specified \(J,M,K\) quantum numbers. A useful factorization is
\[
\hat{P}_{MK}^J = (2J + 1)\, \hat{P}_M^{(z)}\, \hat{P}_{MK}^{(y)J}\, \hat{P}_K^{(z)},
\qquad
\hat{P}_K^{(z)} = \frac{1}{2\pi} \int_0^{2\pi} d\phi\, e^{i(K - \hat{J}_z) \phi},
\]
which makes the \(K\)-projection content explicit [1608.08726].

For a single Hartree–Fock Slater determinant \(|\Phi\rangle\), one builds projected norm and Hamiltonian kernels,
\[
n^J_{M M'} = \langle \Phi| \hat P^J_{M M'} |\Phi\rangle,
\qquad
h^J_{M M'} = \langle \Phi| \hat H \hat P^J_{M M'} |\Phi\rangle,
\]
and then solves the generalized eigenvalue problem
\[
\sum_{M'} h^J_{M M'}\, g^{J,\lambda}_{M'} = E^J_\lambda\, \sum_{M'} n^J_{M M'}\, g^{J,\lambda}_{M'}.
\]
When the intrinsic state is represented by an occupied-orbital matrix \(\Psi\), the rotated determinant is \(\tilde\Psi(\Omega)=R(\Omega)\Psi\), the overlap is
\[
\langle \Phi|\hat R(\Omega)|\Phi\rangle=\det(\Psi^\dagger R(\Omega)\Psi),
\]
and the transition density enters the Hamiltonian kernel through the usual one-body and two-body contractions. This kernel-based formulation is the basic computational pattern for projected Hartree–Fock and related symmetry-restored methods [1304.7292].

## 2. Intrinsic \(K\), \(K\)-distribution, and alignment

In an intrinsic frame, \(K\) denotes the projection of total angular momentum on the body-fixed \(z\)-axis. For a general intrinsic state \(|\Phi\rangle\), the weight of a given \(K\) is
\[
w(K) \equiv \langle \Phi | \hat{P}_K^{(z)} | \Phi \rangle.
\]
The set \(\{w(K)\}\) is the \(K\)-distribution. A narrow \(K\)-distribution is desirable in projected multiconfiguration calculations because each non-negligible \(K\) generates an additional basis state \(\hat P^J_{MK}|\Phi\rangle\), while different \(K\)-components of the same intrinsic state are generally non-orthogonal after projection. Broad \(K\)-distributions therefore enlarge the basis, worsen conditioning of the norm matrix, and can produce severe near-linear dependence [1608.08726].

The orientation of the intrinsic state is decisive. The same physical intrinsic wave function may have a very narrow \(K\)-distribution in one orientation and a very broad one in another. To address this alignment problem, a general prescription is to minimize the variance of \(\hat J_z\),
\[
\sigma^2_{J_z}=\langle \Phi | \hat{J}_z^2 | \Phi \rangle-\langle \Phi | \hat{J}_z | \Phi \rangle^2,
\]
or, more generally, to diagonalize the angular-momentum covariance matrix
\[
V_{\sigma\rho}^J = \left\langle \frac{\hat{J}_\sigma \hat{J}_\rho + \hat{J}_\rho \hat{J}_\sigma}{2} \right\rangle - \langle \hat{J}_\sigma \rangle\, \langle \hat{J}_\rho \rangle.
\]
Because \(\mathsf V^J\) transforms as \(\mathsf V^{J'}=\mathsf R\,\mathsf V^J\,\mathsf R^T\), choosing the body-fixed \(z\)-axis along the eigenvector with the smallest eigenvalue minimizes the \(J_z\) variance and concentrates the \(K\)-distribution. This procedure was benchmarked for an \(\alpha\)–\(^{24}\)Mg cluster structure in \(^{28}\)Si, a triaxially superdeformed state in \(^{40}\)Ar, and Hartree–Fock states in light nuclei. It performs similarly to moment-of-inertia alignment in strongly quadrupole-deformed even-even nuclei, but it is markedly superior for cluster structures, nearly spherical odd and odd–odd nuclei, and cases where angular momentum is carried predominantly by valence nucleons rather than by collective quadrupole deformation. In \(^{17}\)O, moment-of-inertia alignment generated a very broad \(K\)-distribution whose \(K=1/2\) and \(K=3/2\) components after \(J=5/2\) projection had overlap \(\approx 0.9999\), causing numerical pathologies, whereas variance minimization yielded a pure \(K=5/2\) state [1608.08726].

## 3. Configuration mixing and symmetry-restored many-body methods

Angular momentum projection is most powerful when combined with configuration mixing. A general projected multiconfiguration state can be written as
\[
|\Psi_{J_n M}\rangle = \sum_{ij}\hat{P}^{J}_{M K_i}\, |\Phi_j\rangle\, f^{J_n}_{ij},
\]
so that each intrinsic configuration \(|\Phi_j\rangle\) contributes through its projected \(K\)-components. This structure underlies HF–GCM, AMD, FMD, and Monte Carlo shell-model–type calculations [1608.08726].

In triaxial projected GCM with the Gogny force, one combines particle-number and angular-momentum projection with shape mixing in the \((\beta,\gamma)\) plane,
\[
|IM;NZ\,\sigma\rangle =\sum_{K\beta\gamma} f^{I;NZ,\sigma}_{K\beta\gamma}\,|IMK;NZ;\beta\gamma\rangle,
\]
and solves the Hill–Wheeler–Griffin equation built from projected norm and Hamiltonian kernels. The first full triaxial implementation in the \((\beta,\gamma)\) plane with the Gogny force reproduced the low-lying spectrum and transition probabilities of \(^{24}\)Mg reasonably well and identified a ground band and a \(K=2\) \(\gamma\) band. The same work introduced the RVAMPIR approximation, which gives a good description of the ground and gamma bands in the absence of strong mixing [1004.2877].

Projected generator-coordinate methods also expose a methodological distinction between restoration before and after solving the collective secular equation. In ab initio PGCM studies of monopole resonances, performing angular-momentum restoration only after solving the secular equation contaminates the monopole response with an unphysical coupling to the rotational motion, whereas an a priori restoration is necessary to handle rotations and vibrations consistently in the same framework. This result aligns projected GCM with the broader view that symmetry restoration must enter the variational problem itself when collective motions are strongly coupled [2407.01325].

For rotational bands at high spin, the same logic leads to multicranked configuration mixing. Instead of projecting a single intrinsic state, one superposes several cranked mean-field states \(|\Phi(\omega_{\rm rot})\rangle\) after projection, implementing the Peierls–Thouless idea that the conjugate collective momentum must enter the generator set. In superdeformed yrast bands of \(^{152}\)Dy and \(^{194}\)Hg, this multicranked angular-momentum-projection scheme gives a good description of the rotational spectrum as well as \(\mathcal J^{(1)}\) and \(\mathcal J^{(2)}\), and particle-number projection on top of angular-momentum projection was performed for the first time with the multicranked configuration-mixing [1906.11418]. In weakly bound systems, DRHBc+AMP extends projection to deformed relativistic Hartree–Bogoliubov theory in the continuum; with a Dirac Woods–Saxon basis and Bogoliubov pairing, it reproduces the ground-state rotational bands of \(^{36,38,40}\)Mg reasonably well [2107.05925].

## 4. Alternative formulations beyond Euler-angle quadrature

Although the integral projector is standard, angular momentum projection is not tied to numerical integration over Euler angles. One alternative is algebraic. For \(n\) particles with spins \(j_\alpha\), one introduces collective lowering operators
\[
J_-(x_i)=\sum_{\alpha=1}^{n}\frac{1}{x_i-\epsilon_\alpha}\,J_-^{(\alpha)},
\]
and highest-weight states
\[
|\zeta;J,J\rangle = J_-(x^{(\zeta)}_1)\cdots J_-(x^{(\zeta)}_k)|{\rm h.w.}\rangle,
\qquad
J=\sum_\alpha j_\alpha-k.
\]
Requiring \(J_+|\zeta;J,J\rangle=0\) leads to Bethe ansatz equations
\[
\sum_{\alpha=1}^{n}\frac{2j_\alpha}{x_i-\epsilon_\alpha}
-\sum_{t\ne i}\frac{2}{x_i-x_t}=0.
\]
These equations are equivalent to the condition that the \(x_i\) are zeros of degree-\(k\) Heine–Stieltjes polynomials satisfying a second-order Fuchsian equation. The formalism yields a constructive projection method for systems of arbitrary spins, connects angular momentum projection to the strong-pairing limit of Richardson–Gaudin models, and simplifies substantially for identical bosons, where the relevant zeros can be expressed in terms of Jacobi polynomials. For identical bosons or fermions the projected states are generally non-orthogonal in the multiplicity label and require orthonormalization [1307.6497].

A second alternative replaces the three-dimensional Euler-angle integral by the solution of a linear system. The rotated norm and Hamiltonian kernels can be written as finite sums,
\[
\langle \Psi | \hat{R}(\Omega) | \Psi \rangle = \sum_{J,M,K} \mathcal{D}^{(J)}_{M K}(\Omega)\, N^J_{M K},
\]
\[
\langle \Psi | \hat{H}\hat{R}(\Omega) | \Psi \rangle = \sum_{J,M,K} \mathcal{D}^{(J)}_{M K}(\Omega)\, H^J_{M K},
\]
so projection becomes an inversion problem for the unknown \(N^J_{MK}\) and \(H^J_{MK}\). In the implementation tested on \(^{48}\)Cr and \(^{60}\)Fe in the \(pf\) shell, a quadrature with 40 points per Euler angle, \(40^3=64{,}000\) points, agreed with the linear-algebra projection to within about \(2\) keV for all \(J\) except \(J=12\), where the agreement was within about \(10\) keV. A “need-to-know” strategy reduced the number of Euler-angle evaluations for \(^{48}\)Cr from about \(12{,}615\) to about \(4{,}375\), while preserving keV-level agreement. A later convergence study showed that inversion methods, including Fomenko projection and trapezoidal quadrature in the \(z\)-rotation angles, dramatically improve efficiency in many cases, with the optimal choice depending on how many \(J\) values must be projected [1709.04028] [1808.05672].

## 5. Spectroscopy, performance, and known failure modes

Projected Hartree–Fock based on a single Slater determinant is a controlled but limited approximation. A direct benchmark against configuration-interaction shell-model calculations in identical model spaces showed that low-lying excitation spectra for rotational nuclides are well reproduced, whereas spectra for vibrational nuclides and, more generally, the complex spectra for odd-\(A\) and odd-odd nuclides are less well reproduced in detail. In the \(sd\)-shell nucleus \(^{24}\)Mg, the projected-HF ground state lies \(1.78\) MeV above CI, projection recovers \(73\%\) of the correlation energy, and the shifted excitation spectrum has an rms deviation of \(0.077\) MeV. In \(^{30}\)Al, by contrast, the projected-HF ground state is \(6.40\) MeV above CI, projection recovers only about \(25\%\) of the correlation energy, and the shifted excitation spectrum has an rms deviation of \(1.004\) MeV. The same survey found ground-state \(J^\pi\) correctness frequencies of \(15/15\) for even-even \(sd\)-shell nuclei, \(10/10\) for even-even \(pf\)-shell nuclei, but much lower success rates for odd-\(A\) and odd-odd systems [1304.7292].

A common simplifying approximation in deformed QRPA calculations is the needle approximation, which treats different rotations of a strongly deformed intrinsic state as effectively orthogonal. Exact angular-momentum projection in axially deformed proton–neutron QRPA modifies both ground-state energies and transition strengths: it reduces calculated beta-decay half-lives from those that use the needle approximation by up to \(60\%\), and even more when the effects of projection on the ground-state energy are included [2510.16313]. This does not imply that exact projection automatically improves agreement with experiment in an unrecalibrated EDF; rather, it shows that symmetry restoration is a quantitatively large effect that must be incorporated consistently.

The same caution applies to post hoc symmetry restoration in collective-response calculations. In deformed PGCM calculations of monopole resonances, a posteriori angular-momentum restoration can generate large artificial peaks because the unprojected excited states retain overlap with the projected ground state. This suggests that the mathematical order of projection and configuration mixing is not a technicality: when rotations couple strongly to the mode of interest, restoration after the fact can misidentify rotational admixtures as physical response [2407.01325].

## 6. Collective modes, discrete symmetries, and broader scope

Angular momentum projection is also the natural laboratory for rotational modes that depend on triaxiality or discrete point-group symmetry. In microscopic projection from cranked triaxially deformed mean-field states, wobbling bands appear as multiple projected rotational sequences. The resulting spectra reflect the dynamics of the angular-momentum vector in the intrinsic frame and distinguish transverse from longitudinal wobbling, paralleling the phenomenological triaxial-rotor and particle–rotor pictures [1802.06988]. In the same framework, chiral doublet bands arise naturally in \(^{128}\)Cs and \(^{104}\)Rh without any explicit rotor core; the projected bands display the characteristic \(B(E2)\) and \(B(M1)\) properties expected from the triaxial particle–rotor model [1802.06991].

When parity is broken together with rotational symmetry, angular momentum projection can be combined with parity projection. A microscopic study of tetrahedrally symmetric nuclei showed that for pure tetrahedral deformation the projected excitation patterns satisfy the characteristic features predicted by group-representation theory for the tetrahedral symmetry group, including the sequence \(0^+,3^-,4^+,6^\pm,\dots\), and that the energy-versus-spin behavior evolves from approximately linear to the rigid-rotor \(I(I+1)\) pattern as the tetrahedral deformation increases [1301.3279].

In microscopic fission theory, projection is used differently but no less fundamentally. TDDFT fragment states near scission are deformed and oriented, hence they are not eigenstates of \(\hat{\mathbf J}^2\). Applying the standard projector
\[
\hat P^{J}_{MK} = \frac{2J+1}{16\pi^2} \int d\Omega\, \mathcal{D}^{J*}_{MK}(\Omega)\, e^{i\alpha\hat J_z} e^{i\beta\hat J_y} e^{i\gamma\hat J_z}
\]
yields fragment spin distributions
\[
P(J_F)=\sum_{K_F}\langle\Psi|\hat P^{J_F}_{K_FK_F}|\Psi\rangle.
\]
A comparison between projected fragment-spin distributions and a spin cut-off formula derived from the uncertainty relation between orientation angle and angular momentum found that a large portion of the spin distribution obtained from projection methods can be explained by quantum uncertainty associated with fragment orientation, mainly due to quadrupole deformation and, to a lesser extent, octupole deformation [2512.02207].

Taken together, these developments define angular momentum projection as more than a post-processing step. It is the operator-theoretic mechanism by which intrinsic, symmetry-breaking descriptions are converted into laboratory-frame spectroscopy, and its practical success depends on how projection is interleaved with alignment, pairing, configuration mixing, and the numerical representation of the many-body state. The recurring lesson across projected HF, Gogny and relativistic EDF frameworks, PGCM, algebraic projection, and fission applications is that restoring rotational symmetry is both conceptually essential and computationally delicate.

Source: https://www.emergentmind.com/topics/angular-momentum-projection