---
title: Projected Shell Model in Nuclear Physics
url: https://www.emergentmind.com/topics/projected-shell-model
type: topic
---

# Projected Shell Model in Nuclear Physics

Searching arXiv for recent and foundational Projected Shell Model papers to support the article.
The Projected Shell Model (PSM) is a microscopic nuclear-structure framework in which many-body states are constructed from a deformed quasiparticle basis, symmetries broken at the mean-field level are restored by projection, and the resulting projected configurations are mixed by diagonalization of an effective Hamiltonian. In the literature summarized here, the PSM and its triaxial extension (TPSM) are used for rotational spectroscopy, \(\gamma\)- and \(\gamma\gamma\)-band structure, chiral doublets, negative-parity bands, stellar weak-interaction rates, first-forbidden \(\beta\) decay, neutrinoless double-\(\beta\) decay matrix elements, and electromagnetic properties of heavy deformed nuclei such as \(^{229}\)Th [1806.06548; 2105.02649; 2508.19686].

## 1. Formal structure of the model

The common PSM starting point is a deformed Nilsson single-particle basis with BCS pairing. The intrinsic states are quasiparticle configurations \(|\Phi_\kappa\rangle\) built on a deformed quasiparticle vacuum, and rotational symmetry is restored by the three-dimensional angular-momentum projection operator
$$
\hat P^I_{MK}=\frac{2I+1}{8\pi^2}\int d\Omega\, D^{I*}_{MK}(\Omega)\,\hat R(\Omega).
$$
Projected configurations are then mixed to form laboratory-frame eigenstates with good total angular momentum:
$$
|\Psi^I_M\rangle=\sum_{\kappa K} F_{\kappa K}^I \hat P^I_{MK}|\Phi_\kappa\rangle.
$$
Equivalent forms of this ansatz are used throughout the axial and triaxial implementations, including odd-mass and odd-odd systems [2105.02649; 1806.06548].

The shell-model Hamiltonian is typically of pairing-plus-quadrupole form,
$$
\hat H = \hat H_0 - \frac{1}{2}\chi \sum_\mu \hat Q^\dagger_\mu \hat Q_\mu - G_M \hat P^\dagger \hat P - G_Q \sum_\mu \hat P^\dagger_\mu \hat P_\mu,
$$
with \(\hat H_0\) the spherical single-particle term and the residual interaction comprising quadrupole-quadrupole, monopole-pairing, and quadrupole-pairing terms. In heavy nuclei and octupole-deformed systems, the Hamiltonian has also been augmented to include \(L=2\) to \(4\) multipole channels,
$$
\hat H = \hat H_0 - \frac{1}{2}\sum_{L=2}^{4}\chi_L \sum_{\omega=-L}^{L}\hat Q_{L\omega}^\dagger \hat Q_{L\omega} - G_M \hat P^\dagger \hat P - G_Q \sum_{\omega=-2}^{2}\hat P_{2\omega}^\dagger \hat P_{2\omega},
$$
which was used in the microscopic study of \(^{229}\)Th [2508.19686].

Diagonalization in the nonorthogonal projected basis is carried out through a Hill–Wheeler-type generalized eigenvalue problem. The same machinery is used in axial PSM, TPSM, and number-projected variants. Some implementations enforce particle number only on average through BCS, whereas others employ explicit neutron and proton number projection before angular-momentum projection [2508.19686; 1806.06548].

## 2. Configuration spaces and major extensions

A central feature of the PSM is controlled truncation in a deformed quasiparticle basis. For even-even nuclei, standard TPSM calculations use projected \(0\)-qp, \(2\)-qp, and \(4\)-qp states,
$$
\hat P^I_{MK}\ket{\Phi},\quad
\hat P^I_{MK} a^\dagger_{p_1} a^\dagger_{p_2}\ket{\Phi},\quad
\hat P^I_{MK} a^\dagger_{n_1} a^\dagger_{n_2}\ket{\Phi},\quad
\hat P^I_{MK} a^\dagger_{p_1} a^\dagger_{p_2} a^\dagger_{n_1} a^\dagger_{n_2}\ket{\Phi},
$$
with the \(K=0,2,4\) projections from the triaxial vacuum generating the ground, \(\gamma\), and \(\gamma\gamma\) bands, respectively [2009.04135; 1401.1737].

Odd-mass extensions replace the \(0\)-qp vacuum bandhead by one-quasiparticle configurations and include higher broken-pair structures at high spin. For \(^{103}\)Nb, the odd-proton TPSM basis was extended to
$$
\left\{ \hat P^I_{MK} a^\dagger_{p} \lvert \Phi \rangle,\; \hat P^I_{MK} a^\dagger_{p} a^\dagger_{n1} a^\dagger_{n2} \lvert \Phi \rangle \right\},
$$
while in odd-neutron Xe isotopes the basis was enlarged further to include one-quasineutron, one-neutron-plus-two-proton, three-quasineutron, and five-quasiparticle configurations, enabling description up to and including the second band crossing [1004.1873; 2203.12851].

Odd-odd implementations use proton-neutron two-quasiparticle configurations,
$$
\{ |\phi_\kappa\rangle = a^\dagger_\nu a^\dagger_\pi |0\rangle \},
$$
which form the basis for chiral doublet studies in \(^{128}\)Cs and \(^{124,126,130,132}\)Cs, as well as the multi-configuration PSM treatment of \(^{130}\)Cs [1112.3749; 1312.6963; 1806.06548].

Several formal enlargements have broadened the scope of the method. A negative-parity TPSM for even-even nuclei allowed quasiparticles to occupy two major oscillator shells of opposite parity, thereby making negative-parity rotational bands accessible microscopically [2307.14827]. In stellar weak-interaction work, one extended PSM included up to six-quasiparticle configurations for rates in \(rp\)-process waiting-point nuclei [2308.11297], while an electron-capture study of \(^{93}\)Nb employed explicit multi-quasiparticle configurations up to seven quasiparticles [2111.14288]. In heavy odd-mass \(^{229}\)Th, the configuration space was truncated at the 3-quasiparticle level, with the authors stating that 5-qp contributions are small [2508.19686].

## 3. Collective spectroscopy, \(\gamma\) bands, and high-spin structure

In the TPSM, \(\gamma\)-band structure emerges from projection rather than from a phonon degree of freedom introduced by hand. For even-even nuclei, the staggering quantity
$$
S(I)= \frac{[E(I)-E(I-1)]-[E(I-1)-E(I-2)]}{E(2^{+}_1)}
$$
is used as a diagnostic of the \(\gamma\)-band phase. A systematic study of twenty-three nuclei found that most exhibit even-\(I\)-below-odd-\(I\) staggering after quasiparticle configuration mixing, while only \(^{76}\)Ge, \(^{112}\)Ru, \(^{170}\)Er, and \(^{232}\)Th retain odd-\(I\)-below-even-\(I\) staggering [2009.04135].

The Ge/Se region provides a particularly explicit TPSM test of rigid triaxiality against \(\gamma\)-softness. For \(^{76}\)Ge, the best description of the yrast and \(\gamma\)-vibrational bands required \(\epsilon=0.20\) and \(\epsilon'=0.16\), corresponding to \(\gamma\approx 30^\circ\), and the calculation predicted a \(\gamma\gamma\) band around \(2.5\) MeV. In neighboring Ge and Se isotopes, the opposite staggering phase emerged after configuration mixing, and the paper argued that the resulting \(\gamma\)-soft feature comes from mixing of the ground-state configuration with multi-quasiparticle states [1401.1737].

In odd-mass nuclei, the same projected-\(K\) logic yields multi-\(\gamma\) bands built on a quasiparticle configuration. In \(^{103}\)Nb, the ground, one-\(\gamma\), and two-\(\gamma\) bands were identified primarily with the \(K=5/2\), \(9/2\), and \(13/2\) projections of the same lowest one-quasiparticle configuration. The calculated bandhead ratio \(E_{\gamma\gamma}/E_\gamma=2.01\), compared with the experimental value \(1.87\), supported the assignment of simultaneous one- and two-\(\gamma\)-phonon bands while also indicating residual anharmonicity [1004.1873].

High-spin rotational behavior is likewise treated through band diagrams and configuration mixing. In the superdeformed band of \(^{40}\)Ar, an axial PSM with \(0\)-, \(2\)-, and \(4\)-qp configurations reproduced the nearly linear \(\gamma\)-ray energy sequence and showed that the high-spin yrast structure is dominated by mixed \(0\)-, \(2\)-, and \(4\)-qp configurations. The triaxial deformation extracted from the calculated observables was small, roughly \(2^\circ\)–\(6^\circ\), leading to the conclusion that triaxiality is not very important for this superdeformed band [1508.04055].

A recurring result in comparative studies is that energies alone do not exhaust the model discrimination. When a collective Bohr Hamiltonian was fitted to TPSM energies for \(\gamma\)-bands, the staggering pattern could be reproduced reasonably well, but interband \(B(E2)\) values differed strongly. This suggests that transition strengths, not only level energies, are required to distinguish fixed-triaxial configuration mixing from collective \(\gamma\)-soft dynamics [2009.04135].

## 4. Chiral rotation and triaxial angular-momentum geometry

The PSM and TPSM provide a microscopic description of nuclear chirality by projecting good angular momentum from triaxial proton-neutron quasiparticle configurations. In \(^{128}\)Cs, a TPSM calculation with \(\epsilon=0.220\) and \(\epsilon'=0.14\), corresponding to \(\gamma\sim 30^\circ\), reproduced the observed energies and electromagnetic transition probabilities of the nearly degenerate chiral dipole bands. The splitting between the two bands was shown to depend strongly on \(\epsilon'\), becoming small for \(\epsilon' \gtrsim 0.125\), which confirmed that triaxiality is essential for chirality in this case [1112.3749].

The odd-odd Cs isotopes \(^{124,126,130,132}\)Cs were studied with a multi-quasiparticle TPSM in terms of energies, staggering, \(B(E2)\), \(B(M1)\), and \(B(M1)/B(E2)\) ratios. An explicit methodological point in that work is that near degeneracy alone is not sufficient to establish chirality; the electromagnetic transition patterns are also required. For \(^{126}\)Cs in particular, the agreement with extensive transition-probability data was taken to support the chiral interpretation, while \(^{124,130,132}\)Cs were used to generate predictions for future tests [1312.6963].

The most detailed intrinsic-geometry analysis in the material summarized here is the PSM study of \(^{130}\)Cs with configuration mixing. In that work, chiral geometry was analyzed through the \(K\) plot, which displays distributions of angular-momentum components along the intrinsic axes, and the azimuthal plot, which gives distributions of orientation angles \((\theta,\phi)\) of the total angular momentum. The resulting picture evolved with spin from chiral vibration at \(I=10,11\hbar\), to static chirality at \(I=14,15\hbar\), and then to principal-axis rotation at \(I=18,19\hbar\). The authors concluded that the chiral geometry is stable against configuration mixing [1806.06548].

These studies jointly establish that, within the projected-shell-model framework, chiral doublets are not inferred from energy quasi-degeneracy alone. Rather, they are identified through a conjunction of projected triaxial geometry, band mixing, and characteristic \(E2\) and \(M1\) systematics [1112.3749; 1312.6963; 1806.06548].

## 5. Weak interactions, \(\beta\) decay, and stellar rates

The PSM has been extended beyond spectroscopy to weak-interaction observables in thermally populated nuclei. For stellar electron capture, the rate was written as
$$
\lambda^{\rm EC} = \frac{\ln 2}{K} \sum_i \frac{(2J_i+1)e^{-E_i/(k_B T)}}{G(Z,A,T)} \sum_f B_{if}\,\Phi_{if}^{\rm EC},
$$
with the partition function
$$
G(Z,A,T)=\sum_i (2J_i+1)e^{-E_i/(k_B T)}.
$$
In the \(^{93}\)Nb \(\rightarrow\) \(^{93}\)Zr example, the PSM was used precisely because excited configurations are explicitly included as multi-quasiparticle states. That study concluded that the Brink-Axel approximation is not reliable for odd-mass nuclei such as \(^{93}\)Nb, and that at \(T_9 \approx 15\) the inclusion of excited parent states can increase the electron-capture rate by more than an order of magnitude [2111.14288].

A broader stellar weak-interaction application proposed a PSM for even-even and odd-odd nuclei with an extended configuration space including up to six-quasiparticle configurations, and calculated stellar weak-interaction rates for eight \(rp\)-process waiting-point nuclei: \(^{64}\)Ge, \(^{68}\)Se, \(^{72}\)Kr, \(^{76}\)Sr, \(^{80}\)Zr, \(^{84}\)Mo, \(^{88}\)Ru, and \(^{92}\)Pd. In that work, higher-order quasiparticle configurations were found to affect the underlying Gamow-Teller strength distributions and the corresponding stellar weak-interaction rates. Under \(rp\)-process environments with high temperatures and densities, thermal population of excited parent states tends to decrease the stellar \(\beta^+\) decay rates, whereas electron capture tends to contribute increasingly with temperature and density. The effective half-lives under peak \(rp\)-process conditions were predicted to be reduced relative to the terrestrial case, especially for \(^{64}\)Ge and \(^{68}\)Se [2308.11297].

The operator content has also been expanded. A PSM development for first-forbidden nuclear \(\beta\) decay was presented as a first application of the model to that transition class, with 35 dominant first-forbidden transitions calculated and compared systematically with data. The corresponding experimental Log\(f_0 t\) values were described reasonably, and quenching factors of nuclear matrix elements were reported to affect both the Log\(f_0 t\) values and the related shape factors [2310.19523].

Within this body of work, the PSM serves not merely as a spectroscopic model but as a state-by-state framework for weak processes in nuclei where thermal population, deformation, and explicit excited configurations are all structurally important [2111.14288; 2308.11297].

## 6. Double-\(\beta\) decay, electromagnetic observables, and computational techniques

The triaxial projected shell model has also been used for neutrinoless double-\(\beta\) decay. In that application, mother and daughter states were built from triaxially deformed quasiparticle vacua with angular-momentum projection and quasiparticle configuration mixing, and odd-odd intermediate nuclei were treated explicitly as two-quasiparticle proton-neutron configurations. The decay matrix element was decomposed as
$$
M^{0\nu}= M_F^{0\nu}-M_{GT}^{0\nu}+M_T^{0\nu},
$$
and the numerical analysis found that the total NME is dominated by the Gamow–Teller channel at about \(85\%\), with the Fermi term contributing about \(14\%\) and the tensor term about \(1\%\). Configuration mixing reduced the NME by \(2\%\)–\(7\%\) depending on nucleus, explicit odd-odd intermediate states enhanced it by \(4\%\)–\(20\%\), and varying the triaxial deformation \(\gamma\) from \(0^\circ\) to \(60^\circ\) changed the NME by \(14\%\) to \(511\%\), depending on the decay candidate [2105.02649].

Electromagnetic observables are evaluated directly from projected wavefunctions. Standard reduced transition probabilities,
$$
B(E2; I_i\rightarrow I_f)=\frac{1}{2I_i+1}\left|\langle \Psi^{I_f}\|\hat Q_2\|\Psi^{I_i}\rangle\right|^2,
$$
and
$$
B(M1; I_i\rightarrow I_f)=\frac{\mu_N^2}{2I_i+1}\left|\langle \sigma_f, I_f \| \hat{\mathcal M}_1 \| \sigma_i, I_i\rangle\right|^2,
$$
are used repeatedly to diagnose collectivity, alignment, and band mixing. In \(^{229}\)Th, a modern PSM treatment with quadrupole, octupole, and hexadecapole deformation reproduced the low-energy spectrum and obtained
$$
B(M1; ^{229m}\mathrm{Th}\to \text{g.s.}) = 0.0240\ \text{W.u.},
$$
which was stated to agree well with the recently measured radiative lifetime of \(^{229m}\)Th. The same calculation found that the \(5/2^+\) ground state is dominated by the \(n\,5/2[633]\) one-quasiparticle component at about \(77\%\), while the \(3/2^+\) isomer is dominated by the \(n\,3/2[631]\) component at about \(83\%\), thereby explaining large in-band \(E2\) strengths and weak cross-band \(E2\) transitions in terms of small wave-function overlap [2508.19686].

A significant computational ingredient in modern implementations is the Pfaffian algorithm. It is used to compute rotational overlaps in the TPSM treatment of neutrinoless double-\(\beta\) decay and to evaluate reduced one-body transition densities in the \(^{229}\)Th calculation [2105.02649; 2508.19686]. This suggests that efficient overlap technology is integral to extending the projected-shell-model program to odd-mass systems, large configuration spaces, and transition operators beyond standard spectroscopy.

The literature also records explicit limitations. In \(^{103}\)Nb, the remaining discrepancy in the two-\(\gamma\) band was attributed to \(\gamma\)-softness and to the use of a single fixed \(\epsilon'\), with projection over multiple \(\epsilon'\) values and generator-coordinate mixing proposed as a better treatment [1004.1873]. In \(^{130}\)Cs, one observed band was argued to lie beyond the present model space because it likely has larger deformation [1806.06548]. Such cases indicate that PSM accuracy depends not only on projection and configuration mixing, but also on the adequacy of the chosen intrinsic manifold.

Source: https://www.emergentmind.com/topics/projected-shell-model