---
title: Cranked Nilsson-Strutinsky Formalism
url: https://www.emergentmind.com/topics/cranked-nilsson-strutinsky-formalism
type: topic
---

# Cranked Nilsson-Strutinsky Formalism

Searching arXiv for recent and relevant CNS papers to ground the article.
arXiv_search query: "cranked Nilsson-Strutinsky formalism rotating nuclei"
The cranked Nilsson–Strutinsky (CNS) formalism is a rotating-frame description of high-spin nuclear structure in which one diagonalizes a cranked Nilsson Hamiltonian, evaluates a Strutinsky shell correction, adds a rotating-liquid-drop reference, and minimizes the resulting total energy with respect to deformation variables such as $(\varepsilon_2,\gamma,\varepsilon_4)$ at fixed rotational frequency $\omega$ or spin $I$. In the literature represented here, the term encompasses both the standard unpaired CNS approach and paired extensions such as Lipkin–Nogami and Cranked Nilsson–Strutinsky–Bogoliubov (CNSB); it is also used more loosely for schematic “CNS-like” Hamiltonians that retain fixed single-particle levels, pairing, cranking, and particle-number projection while omitting the explicit anisotropic-oscillator Nilsson potential and the Strutinsky smoothing machinery [1106.1769] [2506.18059].

## 1. Formal definition in the rotating frame

At the level of the one-body Routhian, the defining CNS step is the introduction of uniform rotation about a principal axis, usually the $x$-axis, through
\[
H'(\omega)=H_{\rm Nilsson}-\omega J_x .
\]
Equivalent formulations appear as $h'=H_{\rm Nilsson}-\omega_x J_x$ or $H^\omega_{\rm Nilsson}=H_{\rm osc}+V_{\rm Nil}-\omega\!\cdot\!J$, with $\omega$ treated as an external parameter and $J_x$ the $x$-component of the total angular-momentum operator [1106.1769] [1002.2298].

The Nilsson part is a deformed harmonic-oscillator mean field supplemented by spin–orbit and $\ell^2$ terms. Representative forms used in the cited applications are
\[
H_{\rm Nilsson}
=\frac{p^2}{2m}
+\tfrac12\,m\!\bigl(\omega_\perp^2(x^2+y^2)+\omega_z^2z^2\bigr)
-2\kappa\,\hbar\omega_0\,\boldsymbol{\ell}\!\cdot\!\boldsymbol{s}
-\kappa\mu\,\hbar\omega_0\,\ell^2 ,
\]
and
\[
H^\omega
=
h_{\rm HO}(\varepsilon_2,\gamma)
-\kappa\hbar\omega_0\bigl[2\,\ell_t\!\cdot\!s+\mu(\ell_t^2-\langle \ell_t^2\rangle_N)\bigr]
+V_4(\varepsilon_4,\gamma)
-\omega j_x .
\]
The oscillator frequencies are chosen to encode quadrupole deformation, with volume conservation imposed in the triaxial parametrizations. In the Lund convention, $\gamma\in[-60^\circ,+60^\circ]$, with $\gamma=0^\circ$ corresponding to prolate and $\gamma=60^\circ$ to oblate shapes [1205.5986] [2509.01491].

A central practical consequence of this construction is that the theory generates single-particle Routhians $e_i'(\omega,\varepsilon)$ or $e_i(\omega;\varepsilon_2,\gamma,\varepsilon_4)$ whose occupancies define the many-body configuration. In configuration-constrained implementations, those occupancies are explicitly tracked as a function of spin and deformation [2001.01922] [1205.5986].

## 2. Nilsson mean field, deformation variables, and rotational geometry

The deformation space of CNS calculations is typically three-dimensional in $(\varepsilon_2,\gamma,\varepsilon_4)$, although some studies extend the hexadecapole sector to a full five-dimensional minimization in $(\varepsilon_2,\gamma,\varepsilon_{40},\varepsilon_{42},\varepsilon_{44})$. The quadrupole variable $\varepsilon_2$ controls elongation, $\gamma$ controls triaxiality, and $\varepsilon_4$ or the $\varepsilon_{4i}$ specify hexadecapole shape degrees of freedom [1205.5986] [2209.08905].

Several equivalent parameterizations occur in the cited literature. For example, one formulation writes
\[
\omega_x=\omega_0(\varepsilon_2,\gamma)\Bigl[1-\tfrac23\,\varepsilon_2\cos\bigl(\gamma+\tfrac{2\pi}3\bigr)\Bigr],\quad
\omega_y=\omega_0(\varepsilon_2,\gamma)\Bigl[1-\tfrac23\,\varepsilon_2\cos\bigl(\gamma-\tfrac{2\pi}3\bigr)\Bigr],\quad
\omega_z=\omega_0(\varepsilon_2,\gamma)\Bigl[1-\tfrac23\,\varepsilon_2\cos\gamma\Bigr],
\]
while another uses
\[
\omega_\perp=\omega_0\Bigl(1+\tfrac{2}{3}\,\varepsilon_2\cos\gamma\Bigr),\qquad
\omega_z=\omega_0\Bigl(1-\tfrac{4}{3}\,\varepsilon_2\cos\gamma\Bigr).
\]
The formal content is the same: the cranked Hamiltonian is diagonalized in a deformed oscillator basis at each chosen shape and frequency [1106.1769] [1205.5986].

This deformation dependence is not merely kinematic. In the applications summarized here, equilibrium shapes are identified through local or global minima on Total Routhian Surfaces (TRS) or potential-energy surfaces. The resulting minima can correspond to near-prolate, near-oblate, triaxial short-axis, intermediate-axis, or long-axis rotation, depending on the sign and magnitude of $\gamma$. For $^{140}$Gd, the pairing-independent CNS calculation favored $\gamma=-30^\circ$ over a substantial spin range, while in $^{76}$Rb the same configuration [3,4] supported both a near-prolate and a near-oblate minimum [2209.08905] [1106.1769].

## 3. Strutinsky shell correction and total Routhian surfaces

The Strutinsky step separates the total energy into a smooth macroscopic part and an oscillatory shell contribution. Representative expressions are
\[
E_{\rm tot}(\varepsilon_2,\gamma;\omega)
=
E_{\rm LD}(\varepsilon_2,\gamma)
+\delta E_{\rm shell}(\varepsilon_2,\gamma;\omega)
+\delta E_{\rm pair}(\omega),
\]
\[
E_{\rm tot}(I;\varepsilon_2,\gamma,\varepsilon_4)
=
E_{\rm LD}(\varepsilon_2,\gamma,\varepsilon_4)
+\delta E_{\rm shell},
\]
and
\[
E_{\rm tot}(I)=E_{\rm sh}(I)+E_{\rm rld}(I).
\]
The shell correction is written as
\[
\delta E_{\rm shell}=\sum_i e_i-\widetilde E
\]
or, equivalently,
\[
\delta E_{\rm shell}
=
\sum_{i=1}^{N_{\rm occ}} e_i
-\int_{-\infty}^{e_F} e\,\tilde g(e)\,de,
\]
where $\tilde g(e)$ is the smoothed level density obtained by folding the discrete spectrum with a smoothing kernel [1106.1769] [2209.08905] [1205.5986].

In the cited CNS applications, standard Strutinsky choices recur. A commonly used smoothing width is approximately $1.2\,\hbar\omega_0$, sometimes stated as $(1.2\text{–}1.4)\,\hbar\omega_0$, and the curvature-correction polynomial order is typically $p=6$ [1106.1769] [1002.2298] [1609.00294]. The rotating-liquid-drop contribution is implemented either generically as $E_{\rm LD}$ or, in some studies, through the Lublin–Strasbourg Drop with radius parameter $r_0=1.16\,$fm and diffuseness $a=0.6\,$fm [2001.01922] [1205.5986].

After the shell correction is formed, the theory minimizes the total Routhian over deformation space. This minimization is the source of TRS and related contour plots. The deformation meshes differ by application but are explicitly tabulated in several cases: for $^{62}$Ni, $\varepsilon_2=0.00$ to $0.60$ in steps of $0.01$, $\gamma=-60^\circ$ to $+60^\circ$ in steps of $2^\circ$, and $\varepsilon_4=-0.10$ to $+0.10$ in steps of $0.01$; for $^{140}$Gd, a general scan over $\varepsilon_2=0.0$–0.35, $\gamma=-120^\circ$–$+60^\circ$, and $\varepsilon_4\approx-0.05$–0.05 was used; for $^{115}$I, a mesh such as $\varepsilon_2=0.10\!:\!0.02\!:\!0.30$, $\gamma=-60^\circ\!:\!6^\circ\!:+60^\circ$, and $\varepsilon_4=-0.02\!:\!0.01\!:+0.02$ was adopted [1609.00294] [2209.08905] [2509.01491].

## 4. Pairing, particle-number treatment, and paired extensions

A frequent misconception is that CNS is intrinsically pairing-free. The cited literature shows a more differentiated situation. In the pure high-spin CNS variant, static pairing is neglected and $\delta E_{\rm pair}=0$; this approximation is explicitly invoked for $^{107}$In, $^{140}$Gd, $^{62}$Ni, and the high-spin applications in the $A\sim160$ region [1002.2298] [2209.08905] [1609.00294] [1205.5986].

Other studies incorporate pairing corrections or fully paired variants. In the $^{76}$Rb analysis, an “unpaired” CNS calculation was supplemented by a Lipkin–Nogami extension with particle-number projection, minimizing
\[
E_{\rm pair}
=
\sum_i 2\,v_i^2\,e_i
-G\Bigl(\sum_i u_i v_i\Bigr)^2
-\lambda\Bigl(\sum_i v_i^2-N\Bigr)
+\lambda_2\Bigl(\sum_i v_i^2-N\Bigr)^2 .
\]
For the odd–odd nucleus $^{76}$Rb, the pairing corrections were found to be small ($\lesssim0.5$ MeV at low spin) and to decrease with increasing $\omega$ [1106.1769].

The $^{167}$Lu work presents the paired CNSB formalism. There the cranked Hamiltonian is extended to
\[
H'
=
H_{\rm Nilsson}-\omega j_x+\Delta(P^\dagger+P)-\lambda\hat N ,
\]
with
\[
\Delta=G\sum_k u_k v_k,\qquad N=\sum_k v_k^2 ,
\]
and an associated total Routhian
\[
E_{\rm CNSB}(\omega,\varepsilon)
=
E_{\rm LD}(\omega,\varepsilon)
+\sum_i^{\rm qp} E_{{\rm qp},i}
-\widetilde E
-\frac{\Delta^2}{G}.
\]
The same study also used an empirical average pairing shift in the unpaired CNS comparison,
\[
E_{\rm pair}^{\rm avg}(I)=-2.47\,e^{-I/29}\ {\rm MeV},
\]
to mimic the overall pairing drop [2001.01922].

A second misconception is that particle number enters only implicitly. The schematic quantum-variational study based on CNS-like Hamiltonians included an explicit particle-number-projection term,
\[
H_{\rm PNP}=\lambda(\hat N-N_0)^2,
\]
used to drive the variational solution onto the desired particle-number subspace [2506.18059].

## 5. Configuration constraints, observables, and spectroscopic interpretation

A defining practical strength of CNS is configuration tracking. The cited papers use several labeling schemes, all based on occupancies of specific shells, intruders, holes, or grouped orbitals. Examples include the $[p,n]$ notation in $^{76}$Rb, where $[3,4]$ denotes three $g_{9/2}$ protons and four $g_{9/2}$ neutrons; the $[p_1,p_2;\,n_1,n_2]$ notation in $^{62}$Ni, where $D1:[2\,0;0\,2]$ and $D2:[2\,0;0\,3]$ are defined relative to a $^{56}$Ni core; and the more elaborate grouped-shell notation in $^{167}$Lu$\,:\,[p_1(p_2p_3);\,(n_1n_2)n_3]$ [1106.1769] [1609.00294] [2001.01922].

The formalism is then confronted with experiment through energies relative to a rotating-liquid-drop reference,
\[
E'(I)=E_{\rm tot}(I)-E_{\rm RLD}(I),
\]
alignments,
\[
i_x(\omega)=\langle J_x\rangle(\omega),
\]
dynamic moments,
\[
J^{(2)}(\omega)=\frac{d\,i_x}{d\omega},
\]
and, when deformation information is available, transition quadrupole moments $Q_t$. In $^{140}$Gd$,\,Q_t$ values extracted by the Doppler Shift Attenuation Method were compared to the CNS prediction $Q_t\approx4.8\,$eb for $\varepsilon_2=0.24,\gamma=-30^\circ$; in $^{62}$Ni, the measured $Q_t^{\rm exp}=2.2^{+1.1}_{-0.8}\ e\,{\rm b}$ for band $D1$ was compared to a calculated $Q_t^{\rm CNS}\approx1.6$–$1.8\ e\,{\rm b}$ [2209.08905] [1609.00294].

Several recurrent spectroscopic patterns emerge in the applications.

| Nucleus | CNS interpretation | Reported consequence |
|---|---|---|
| $^{76}$Rb | Same [3,4] configuration at different shapes | near-prolate and near-oblate coexistence |
| $^{107}$In | [20,2] configuration | terminating $\Delta I=2$ band |
| $^{140}$Gd | [4,8] unpaired configuration | triaxial minimum at $\gamma=-30^\circ$ |
| $^{62}$Ni | [20,02] and [20,03] | prolate bands driven by $\nu g_{9/2}$ |
| $^{167}$Lu | unpaired CNS and paired CNSB | AB and BC crossings disentangled |
| $^{115}$I | valence-space band 6 | smooth termination at $I_{\max}=67/2^+$ |

For $^{76}$Rb, the [3,4] configuration exhibited two well-developed minima in the $(\varepsilon_2,\gamma)$ plane, one near-prolate with $\gamma\approx0^\circ,\varepsilon_2\approx0.30$ and one near-oblate with $\gamma\approx60^\circ,\varepsilon_2\approx0.23$, separated by a barrier of order $1\,$MeV. For $^{107}$In, the observed $\Delta I=2$ band was described as the [20,2] configuration terminating at $I^\pi=33.5^+$. For $^{115}$I, the smooth-termination criterion was expressed through the maximum spin
\[
I_{\max}=\sum_{\rm occupied} j_i ,
\]
and the assigned valence-space configuration led to $I_{\max}=67/2$ [1106.1769] [1002.2298] [2509.01491].

Band crossings are another major CNS observable. In $^{167}$Lu, the AB crossing at $\hbar\omega\approx0.30\,$MeV corresponds to alignment of a neutron $i_{13/2}$ pair with a typical $\Delta I\approx8\hbar$, whereas the BC crossing at $\hbar\omega\approx0.45\,$MeV gives $\Delta I\approx2\hbar$ in the odd-$A$ case discussed there. The same study concluded that, except for the paired AB and BC crossings, the observed band crossings can be understood within the unpaired formalism [2001.01922].

## 6. Standard approximations, quantitative checks, and methodological boundaries

The CNS formalism is not a single immutable algorithm; it is a family of closely related implementations with identifiable approximations. The most explicit quantitative audit in the cited set is the configuration-constrained study of $A=158$–168 nuclei, which tested the usual neglect of off-shell matrix elements and the standard basis truncation at ${\cal N}_{\rm rot}\le8$ [1205.5986].

In that analysis, neglect of off-shell couplings changed shell corrections by only a few tens of keV throughout the relevant spin range. Extending the basis from ${\cal N}_{\rm max}=8$ to 12 changed the total single-particle energy by up to $\sim120\,$keV and the smoothed energy by $\sim260\,$keV at high spin, so that $\delta E_{\rm sh}$ changed by $\lesssim140\,$keV even at the largest spins. These errors were characterized as negligible compared to other systematic uncertainties such as pairing and the liquid-drop model [1205.5986].

By contrast, the treatment of hexadecapole deformation can be materially important. Restricting the minimization to a single $\varepsilon_4$ parameter changes the energy by only $\lesssim50\,$keV for near-prolate minima, but in the triaxial superdeformed region the full five-dimensional minimization can lower the energy by up to $\sim500\,$keV. The same study also stated that omission of pairing at $I\gtrsim30$ produces a smooth, monotonic shift $\sim0.5\,$MeV without affecting relative alignments or band crossings [1205.5986].

A separate boundary concerns the spin domain of validity. Several cited applications explicitly justify neglect of pairing only in the high-spin regime, for example $I\gtrsim15\hbar$ in $^{107}$In and above $I\approx15\hbar$ in $^{62}$Ni [1002.2298] [1609.00294]. This suggests that the descriptive power of unpaired CNS is strongest where pair correlations are quenched and rotational alignment dominates.

## 7. Schematic CNS-like reductions and quantum-variational reformulations

The 2025 quantum-variational study of rotating nuclei provides a deliberately reduced, schematic version of the CNS structure rather than a full implementation. That work states explicitly that it does not present the full Nilsson Hamiltonian with an explicit anisotropic harmonic-oscillator potential, spin–orbit and $\ell^2$ terms, and it does not carry out a Strutinsky smoothing with curvature corrections. Instead, it assumes fixed single-particle energies $\epsilon_i$ that mock up the Nilsson spectrum and builds Hamiltonians from
\[
H_{\rm sp}=\sum_i \epsilon_i \hat n_i,\qquad
H_{\rm pairing}=-G\sum_i \hat P_i^\dagger \hat P_i,\qquad
H_{\rm crank}=-\omega \hat J_x,
\]
together with the particle-number term
\[
H_{\rm PNP}=\lambda(\hat N-N_0)^2 .
\]
The summary equation given there is
\[
\hat H
=
\sum_i\epsilon_i\,\hat n_i
-G\sum_i\hat P_i^\dagger\hat P_i
-\omega\,\hat J_x
+\lambda(\hat N-N_0)^2 .
\]
The fermionic operators are mapped to qubit Pauli operators via Jordan–Wigner and solved with VQE [2506.18059].

Five increasingly complex CNS-like models are described in that work, beginning from an implied non-interacting level model and progressing through pairing, explicit particle-number projection, and explicit cranking. For Models III and IV, the single-particle energies are $\{0.0,0.2,0.5,0.8\}$, the pairing strength is $G=0.6$, and the cranking frequency is varied over $\omega\in[0.0,1.2]$ in steps of $0.1$. Model II uses a particle-number-projection strength $\lambda$ chosen large enough to fix $N=N_0$ [2506.18059].

The observables used for benchmarking were the ground-state energy, the angular-momentum expectation value $J_x$, and the entanglement entropy. The reported variational errors were typically $<0.005$, with agreement between exact diagonalization and VQE in energy and angular momentum predictions. Slight variances in entanglement entropy were attributed to numerical precision and ansatz expressivity [2506.18059].

This usage clarifies an important terminological boundary. In strict form, CNS denotes the cranked Nilsson Hamiltonian plus Strutinsky shell correction and deformation minimization. In reduced or algorithmic settings, “CNS-like” may instead denote only the retained structural ingredients—single-particle level spacings, pairing correlations, cranking terms, and particle-number conservation—without the explicit Nilsson potential or the Strutinsky shell-correction machinery [2506.18059].

Source: https://www.emergentmind.com/topics/cranked-nilsson-strutinsky-formalism