Cranked Nilsson-Strutinsky Formalism

• Cranked Nilsson-Strutinsky formalism is a rotating-frame approach that combines a deformed harmonic-oscillator potential with spin–orbit and ℓ² terms to model high-spin nuclei.
• The method involves diagonalizing a cranked Nilsson Hamiltonian, implementing Strutinsky shell correction, and minimizing total energy over deformation parameters such as ε₂, γ, and ε₄.
• This formalism is significant for accurately predicting rotational band structures, configuration constraints, and observable quantities like alignments and transition quadrupole moments in nuclear spectroscopy.

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 (Wadsworth et al., 2011, Roy, 22 Jun 2025).

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 (Wadsworth et al., 2011, Ideguchi et al., 2010).

The Nilsson part is a deformed harmonic-oscillator mean field supplemented by spin–orbit and $\omega$0 terms. Representative forms used in the cited applications are

$\omega$1

and

$\omega$2

The oscillator frequencies are chosen to encode quadrupole deformation, with volume conservation imposed in the triaxial parametrizations. In the Lund convention, $\omega$3, with $\omega$4 corresponding to prolate and $\omega$5 to oblate shapes (Kardan et al., 2012, Petrache et al., 1 Sep 2025).

A central practical consequence of this construction is that the theory generates single-particle Routhians $\omega$6 or $\omega$7 whose occupancies define the many-body configuration. In configuration-constrained implementations, those occupancies are explicitly tracked as a function of spin and deformation (Kardan et al., 2020, Kardan et al., 2012).

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

The deformation space of CNS calculations is typically three-dimensional in $\omega$8, although some studies extend the hexadecapole sector to a full five-dimensional minimization in $\omega$9. The quadrupole variable $I$0 controls elongation, $I$1 controls triaxiality, and $I$2 or the $I$3 specify hexadecapole shape degrees of freedom (Kardan et al., 2012, Pai et al., 2022).

Several equivalent parameterizations occur in the cited literature. For example, one formulation writes

$I$4

while another uses

$I$5

The formal content is the same: the cranked Hamiltonian is diagonalized in a deformed oscillator basis at each chosen shape and frequency (Wadsworth et al., 2011, Kardan et al., 2012).

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 $I$6. For $I$7Gd, the pairing-independent CNS calculation favored $I$8 over a substantial spin range, while in $I$9Rb the same configuration [3,4] supported both a near-prolate and a near-oblate minimum (Pai et al., 2022, Wadsworth et al., 2011).

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

$x$0

$x$1

and

$x$2

The shell correction is written as

$x$3

or, equivalently,

$x$4

where $x$5 is the smoothed level density obtained by folding the discrete spectrum with a smoothing kernel (Wadsworth et al., 2011, Pai et al., 2022, Kardan et al., 2012).

In the cited CNS applications, standard Strutinsky choices recur. A commonly used smoothing width is approximately $x$6, sometimes stated as $x$7, and the curvature-correction polynomial order is typically $x$8 (Wadsworth et al., 2011, Ideguchi et al., 2010, Albers et al., 2016). The rotating-liquid-drop contribution is implemented either generically as $x$9 or, in some studies, through the Lublin–Strasbourg Drop with radius parameter $H'(\omega)=H_{\rm Nilsson}-\omega J_x .$0fm and diffuseness $H'(\omega)=H_{\rm Nilsson}-\omega J_x .$1fm (Kardan et al., 2020, Kardan et al., 2012).

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 $H'(\omega)=H_{\rm Nilsson}-\omega J_x .$2Ni, $H'(\omega)=H_{\rm Nilsson}-\omega J_x .$3 to $H'(\omega)=H_{\rm Nilsson}-\omega J_x .$4 in steps of $H'(\omega)=H_{\rm Nilsson}-\omega J_x .$5, $H'(\omega)=H_{\rm Nilsson}-\omega J_x .$6 to $H'(\omega)=H_{\rm Nilsson}-\omega J_x .$7 in steps of $H'(\omega)=H_{\rm Nilsson}-\omega J_x .$8, and $H'(\omega)=H_{\rm Nilsson}-\omega J_x .$9 to $h'=H_{\rm Nilsson}-\omega_x J_x$0 in steps of $h'=H_{\rm Nilsson}-\omega_x J_x$1; for $h'=H_{\rm Nilsson}-\omega_x J_x$2Gd, a general scan over $h'=H_{\rm Nilsson}-\omega_x J_x$3–0.35, $h'=H_{\rm Nilsson}-\omega_x J_x$4–$h'=H_{\rm Nilsson}-\omega_x J_x$5, and $h'=H_{\rm Nilsson}-\omega_x J_x$6–0.05 was used; for $h'=H_{\rm Nilsson}-\omega_x J_x$7I, a mesh such as $h'=H_{\rm Nilsson}-\omega_x J_x$8, $h'=H_{\rm Nilsson}-\omega_x J_x$9, and $$H^\omega_{\rm Nilsson}=H_{\rm osc}+V_{\rm Nil}-\omega\!\cdot\!J$$0 was adopted (Albers et al., 2016, Pai et al., 2022, Petrache et al., 1 Sep 2025).

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 $$H^\omega_{\rm Nilsson}=H_{\rm osc}+V_{\rm Nil}-\omega\!\cdot\!J$$1; this approximation is explicitly invoked for $$H^\omega_{\rm Nilsson}=H_{\rm osc}+V_{\rm Nil}-\omega\!\cdot\!J$$2In, $$H^\omega_{\rm Nilsson}=H_{\rm osc}+V_{\rm Nil}-\omega\!\cdot\!J$$3Gd, $$H^\omega_{\rm Nilsson}=H_{\rm osc}+V_{\rm Nil}-\omega\!\cdot\!J$$4Ni, and the high-spin applications in the $$H^\omega_{\rm Nilsson}=H_{\rm osc}+V_{\rm Nil}-\omega\!\cdot\!J$$5 region (Ideguchi et al., 2010, Pai et al., 2022, Albers et al., 2016, Kardan et al., 2012).

Other studies incorporate pairing corrections or fully paired variants. In the $$H^\omega_{\rm Nilsson}=H_{\rm osc}+V_{\rm Nil}-\omega\!\cdot\!J$$6Rb analysis, an “unpaired” CNS calculation was supplemented by a Lipkin–Nogami extension with particle-number projection, minimizing

$$H^\omega_{\rm Nilsson}=H_{\rm osc}+V_{\rm Nil}-\omega\!\cdot\!J$$7

For the odd–odd nucleus $$H^\omega_{\rm Nilsson}=H_{\rm osc}+V_{\rm Nil}-\omega\!\cdot\!J$$8Rb, the pairing corrections were found to be small ($$H^\omega_{\rm Nilsson}=H_{\rm osc}+V_{\rm Nil}-\omega\!\cdot\!J$$9 MeV at low spin) and to decrease with increasing $\omega$0 (Wadsworth et al., 2011).

The $\omega$1Lu work presents the paired CNSB formalism. There the cranked Hamiltonian is extended to

$\omega$2

with

$\omega$3

and an associated total Routhian

$\omega$4

The same study also used an empirical average pairing shift in the unpaired CNS comparison,

$\omega$5

to mimic the overall pairing drop (Kardan et al., 2020).

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,

$\omega$6

used to drive the variational solution onto the desired particle-number subspace (Roy, 22 Jun 2025).

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 $\omega$7 notation in $\omega$8Rb, where $\omega$9 denotes three $J_x$0 protons and four $J_x$1 neutrons; the $J_x$2 notation in $J_x$3Ni, where $J_x$4 and $J_x$5 are defined relative to a $J_x$6Ni core; and the more elaborate grouped-shell notation in $J_x$7Lu$J_x$8 (Wadsworth et al., 2011, Albers et al., 2016, Kardan et al., 2020).

The formalism is then confronted with experiment through energies relative to a rotating-liquid-drop reference,

$J_x$9

alignments,

$x$0

dynamic moments,

$x$1

and, when deformation information is available, transition quadrupole moments $x$2. In $x$3Gd$x$4 values extracted by the Doppler Shift Attenuation Method were compared to the CNS prediction $x$5eb for $x$6; in $x$7Ni, the measured $x$8 for band $x$9 was compared to a calculated $\omega$00–$\omega$01 (Pai et al., 2022, Albers et al., 2016).

Several recurrent spectroscopic patterns emerge in the applications.

Nucleus	CNS interpretation	Reported consequence
$\omega$02Rb	Same [3,4] configuration at different shapes	near-prolate and near-oblate coexistence
$\omega$03In	[20,2] configuration	terminating $\omega$04 band
$\omega$05Gd	[4,8] unpaired configuration	triaxial minimum at $\omega$06
$\omega$07Ni	[20,02] and [20,03]	prolate bands driven by $\omega$08
$\omega$09Lu	unpaired CNS and paired CNSB	AB and BC crossings disentangled
$\omega$10I	valence-space band 6	smooth termination at $\omega$11

For $\omega$12Rb, the [3,4] configuration exhibited two well-developed minima in the $\omega$13 plane, one near-prolate with $\omega$14 and one near-oblate with $\omega$15, separated by a barrier of order $\omega$16MeV. For $\omega$17In, the observed $\omega$18 band was described as the [20,2] configuration terminating at $\omega$19. For $\omega$20I, the smooth-termination criterion was expressed through the maximum spin

$\omega$21

and the assigned valence-space configuration led to $\omega$22 (Wadsworth et al., 2011, Ideguchi et al., 2010, Petrache et al., 1 Sep 2025).

Band crossings are another major CNS observable. In $\omega$23Lu, the AB crossing at $\omega$24MeV corresponds to alignment of a neutron $\omega$25 pair with a typical $\omega$26, whereas the BC crossing at $\omega$27MeV gives $\omega$28 in the odd-$\omega$29 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 (Kardan et al., 2020).

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 $\omega$30–168 nuclei, which tested the usual neglect of off-shell matrix elements and the standard basis truncation at $\omega$31 (Kardan et al., 2012).

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 $\omega$32 to 12 changed the total single-particle energy by up to $\omega$33keV and the smoothed energy by $\omega$34keV at high spin, so that $\omega$35 changed by $\omega$36keV even at the largest spins. These errors were characterized as negligible compared to other systematic uncertainties such as pairing and the liquid-drop model (Kardan et al., 2012).

By contrast, the treatment of hexadecapole deformation can be materially important. Restricting the minimization to a single $\omega$37 parameter changes the energy by only $\omega$38keV for near-prolate minima, but in the triaxial superdeformed region the full five-dimensional minimization can lower the energy by up to $\omega$39keV. The same study also stated that omission of pairing at $\omega$40 produces a smooth, monotonic shift $\omega$41MeV without affecting relative alignments or band crossings (Kardan et al., 2012).

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 $\omega$42 in $\omega$43In and above $\omega$44 in $\omega$45Ni (Ideguchi et al., 2010, Albers et al., 2016). 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 $\omega$46 terms, and it does not carry out a Strutinsky smoothing with curvature corrections. Instead, it assumes fixed single-particle energies $\omega$47 that mock up the Nilsson spectrum and builds Hamiltonians from

$\omega$48

together with the particle-number term

$\omega$49

The summary equation given there is

$\omega$50

The fermionic operators are mapped to qubit Pauli operators via Jordan–Wigner and solved with VQE (Roy, 22 Jun 2025).

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 $\omega$51, the pairing strength is $\omega$52, and the cranking frequency is varied over $\omega$53 in steps of $\omega$54. Model II uses a particle-number-projection strength $\omega$55 chosen large enough to fix $\omega$56 (Roy, 22 Jun 2025).

The observables used for benchmarking were the ground-state energy, the angular-momentum expectation value $\omega$57, and the entanglement entropy. The reported variational errors were typically $\omega$58, 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 (Roy, 22 Jun 2025).

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 (Roy, 22 Jun 2025).