← Part IV: Nuclear Models
Chapter 13

Deformed Nuclei

Nuclear Deformation Parameters

The shape of a deformed nucleus is described by expanding the nuclear surface radius in terms of spherical harmonics. The general parameterization is:

$$R(\theta,\phi) = R_0\left(1 + \sum_{\lambda=1}^{\infty}\sum_{\mu=-\lambda}^{\lambda} \alpha_{\lambda\mu}\, Y_\lambda^\mu(\theta,\phi)\right)$$

The most important multipoles are the quadrupole ($\lambda=2$) and hexadecapole ($\lambda=4$) deformations. The quadrupole deformation parameters $\beta_2$ and$\gamma$ are defined through:

$$\alpha_{20} = \beta_2\cos\gamma, \qquad \alpha_{2,\pm2} = \frac{1}{\sqrt{2}}\beta_2\sin\gamma$$

The parameter $\beta_2$ measures the magnitude of quadrupole deformation, while$\gamma$ specifies the shape type:

  • - $\gamma = 0^\circ$: Prolate (cigar-shaped), axially symmetric
  • - $\gamma = 60^\circ$: Oblate (disc-shaped), axially symmetric
  • - $0^\circ < \gamma < 60^\circ$: Triaxial deformation

The hexadecapole deformation $\beta_4$ adds a higher-order correction to the shape:

$$R(\theta) \approx R_0\left[1 + \beta_2 Y_2^0(\theta) + \beta_4 Y_4^0(\theta)\right]$$

Typical values for well-deformed rare-earth nuclei are $\beta_2 \approx 0.2\text{--}0.35$ and$|\beta_4| \lesssim 0.1$. The volume-conservation condition ensures that the nuclear volume remains constant under deformation: $R_0 = r_0 A^{1/3}$.

Quadrupole Deformation & Moments

The intrinsic (body-frame) electric quadrupole moment is related to the deformation parameter by:

$$Q_0 = \frac{3}{\sqrt{5\pi}}\,Z\,R_0^2\,\beta_2\left(1 + 0.36\,\beta_2 + \cdots\right)$$

The spectroscopic (lab-frame) quadrupole moment for a state with spin $J$ in the ground-state rotational band is related to $Q_0$ by the projection factor:

$$Q_{\text{spec}}(J) = -\frac{3K^2 - J(J+1)}{(J+1)(2J+3)}\,Q_0$$

where $K$ is the projection of $J$ onto the symmetry axis. For the ground-state band of an even-even nucleus ($K=0$):

$$Q_{\text{spec}}(J) = \frac{J}{(2J+3)}\,Q_0$$

Note that $Q_{\text{spec}}(0^+) = 0$ by symmetry, and for $J=2$:$Q_{\text{spec}}(2^+) = \frac{2}{7}Q_0$. A positive $Q_0$ corresponds to prolate shape, while negative indicates oblate.

The Nilsson Model (Deformed Shell Model)

The Nilsson model describes single-particle motion in a deformed potential. The Hamiltonian for an axially symmetric deformed harmonic oscillator is:

$$H = -\frac{\hbar^2}{2m}\nabla^2 + \frac{m}{2}\left[\omega_\perp^2(x^2+y^2) + \omega_z^2 z^2\right] + C\,\boldsymbol{\ell}\cdot\boldsymbol{s} + D\,\ell^2$$

The oscillator frequencies are related to the deformation parameter $\delta$:

$$\omega_\perp = \omega_0\left(1 + \frac{\delta}{3}\right), \qquad \omega_z = \omega_0\left(1 - \frac{2\delta}{3}\right)$$

Each Nilsson level is labeled by the asymptotic quantum numbers $[Nn_z\Lambda]\Omega^\pi$, where $N$ is the total oscillator quantum number, $n_z$ counts quanta along the symmetry axis, $\Lambda$ is the orbital angular momentum projection, and$\Omega = \Lambda + \Sigma$ is the total angular momentum projection ($\Sigma = \pm\frac{1}{2}$).

Key Properties of Nilsson Levels

  • - Each level has a two-fold Kramers degeneracy: $\pm\Omega$
  • - Parity: $\pi = (-1)^N$
  • - For prolate deformation, levels with large $n_z$ (aligned along z) are pushed down
  • - For oblate deformation, levels with small $n_z$ (equatorial orbits) are favored
  • - Level crossings in the Nilsson diagram create new "deformed magic numbers"

The Nilsson diagram shows single-particle energies as functions of the deformation parameter. It successfully predicts ground-state spins and parities of odd-A deformed nuclei by filling Nilsson levels up to the Fermi surface.

Rotational Bands

A deformed nucleus can rotate collectively about an axis perpendicular to the symmetry axis. The rotational energy spectrum for an axially symmetric rotor is:

$$E_{\text{rot}}(J) = \frac{\hbar^2}{2\mathscr{I}}\left[J(J+1) - K(K+1)\right]$$

For $K=0$ bands (ground-state band of even-even nuclei), only even spins appear:$J^\pi = 0^+, 2^+, 4^+, 6^+, \ldots$ The characteristic energy ratios serve as diagnostics for the rotational character:

$$\frac{E(4^+)}{E(2^+)} = \frac{20}{6} = 3.33, \quad \frac{E(6^+)}{E(2^+)} = \frac{42}{6} = 7.00, \quad \frac{E(8^+)}{E(2^+)} = \frac{72}{6} = 12.00$$

For odd-A nuclei, the rotational band is built on a Nilsson orbital with quantum number$K = \Omega$, and the band contains spins $J = K, K+1, K+2, \ldots$ (for $K \neq 1/2$). When $K = 1/2$, a decoupling parameter $a$ modifies the energies:

$$E(J) = \frac{\hbar^2}{2\mathscr{I}}\left[J(J+1) + a(-1)^{J+1/2}\left(J+\tfrac{1}{2}\right)\right]$$

Moment of Inertia

The nuclear moment of inertia lies between two extremes. The rigid-body value for a uniformly charged ellipsoid is:

$$\mathscr{I}_{\text{rigid}} = \frac{2}{5}m_N A R_0^2\left(1 + 0.31\beta_2\right)$$

The irrotational-flow moment of inertia for a superfluid nucleus is:

$$\mathscr{I}_{\text{irrot}} = \frac{9}{8\pi}m_N A R_0^2\,\beta_2^2$$

Experimentally, nuclear moments of inertia are typically 30-50% of the rigid-body value, indicating that nuclei behave as superfluids with significant pairing correlations. The Inglis cranking formula gives:

$$\mathscr{I}_{\text{crank}} = 2\hbar^2 \sum_{i,j} \frac{|\langle i|j_x|j\rangle|^2}{E_j - E_i}$$

where the sum runs over occupied (i) and unoccupied (j) single-particle states. Pairing correlations reduce the moment of inertia below the rigid-body value by creating an energy gap in the single-particle spectrum.

E2 Transition Strengths & B(E2) Values

The reduced transition probability for E2 transitions within a rotational band is directly related to the intrinsic quadrupole moment:

$$B(E2; J_i \to J_f) = \frac{5}{16\pi}\,e^2\,Q_0^2\,|\langle J_i K 2\,0 | J_f K\rangle|^2$$

For the important $0^+ \to 2^+$ ground-state transition:

$$B(E2; 0^+ \to 2^+) = \frac{5}{16\pi}\,e^2\,Q_0^2$$

Deformed nuclei have very large B(E2) values, typically 100-300 Weisskopf units (W.u.), reflecting their collective nature. The Weisskopf single-particle estimate is:

$$B_W(E2) = \frac{1}{4\pi}\left(\frac{3}{5}\right)^2 e^2 R_0^4 \approx 0.0594\,A^{4/3}\;\text{e}^2\text{fm}^4$$

The ratio $B(E2)/B_W(E2)$ measures the collectivity: values $\gg 1$ indicate coherent contributions from many nucleons.

Superdeformation & Hyperdeformation

At very high angular momentum, nuclei can be trapped in secondary minima of the potential energy surface with very large deformations:

  • Normal deformation: $\beta_2 \approx 0.2\text{--}0.3$, axis ratio ~1.3:1
  • Superdeformation: $\beta_2 \approx 0.5\text{--}0.6$, axis ratio ~2:1. First observed in $^{152}$Dy. Characterized by very regular rotational bands with nearly constant $\Delta E_\gamma$ spacing.
  • Hyperdeformation: $\beta_2 \approx 0.8\text{--}1.0$, axis ratio ~3:1. Predicted but extremely difficult to observe experimentally.

Superdeformed bands are stabilized by shell effects at large deformation, where new magic numbers appear (e.g., N or Z = 60, 80 at 2:1 axis ratio). The moment of inertia for superdeformed bands is close to the rigid-body value:

$$\mathscr{I}_{\text{SD}} \approx 0.8\text{--}0.9\,\mathscr{I}_{\text{rigid}}$$

The transition energies in a superdeformed band follow a nearly linear pattern:$E_\gamma(J \to J-2) = \frac{\hbar^2}{\mathscr{I}}(2J-1)$, giving equally spaced gamma rays, which is the experimental signature used to identify superdeformed bands.

Ground-State Correlations

The nuclear ground state contains important correlations beyond the mean-field description. Pairing correlations are described by the BCS theory adapted to nuclei:

$$|\text{BCS}\rangle = \prod_k \left(u_k + v_k\, a_k^\dagger a_{\bar{k}}^\dagger\right)|0\rangle$$

where $v_k^2$ is the occupation probability and $u_k^2 + v_k^2 = 1$. The pairing gap $\Delta$ is determined self-consistently:

$$\Delta = G\sum_k u_k v_k = G\sum_k \frac{\Delta}{2E_k}$$

where $E_k = \sqrt{(\epsilon_k - \lambda)^2 + \Delta^2}$ is the quasiparticle energy and$G$ is the pairing strength. Empirically, $\Delta \approx 12/\sqrt{A}$ MeV. Quadrupole correlations in the ground state lead to zero-point fluctuations in the deformation, which are especially important near shape-transitional nuclei.

Python Simulation: Nilsson Diagram & Rotational Bands

Simplified Nilsson diagram showing single-particle level splitting with deformation, alongside rotational band energies for different deformation parameters.

Deformed Nuclei Analysis

Python

Nilsson diagram and rotational band energies as a function of beta2

script.py104 lines

Click Run to execute the Python code

Code will be executed with Python 3 on the server

Fortran Implementation: Quadrupole Properties

Computes intrinsic quadrupole moments, B(E2) values, and moment of inertia ratios for a selection of well-deformed nuclei. Also examines superdeformation and hyperdeformation.

Deformed Nuclei Quadrupole Analysis

Fortran

Quadrupole moments, B(E2) in Weisskopf units, and superdeformation parameters

program.f9097 lines

Click Run to execute the Fortran code

Code will be compiled with gfortran and executed on the server

Rate this chapter: