Case Study 2 — The Interacting Boson Model: Algebraic Nuclear Structure

Origins: A Bridge Between the Shell Model and Collective Models

By the early 1970s, nuclear structure theory had developed two powerful but seemingly disconnected frameworks. The shell model, rooted in the Pauli principle and individual nucleon orbits, could describe nuclei near closed shells with impressive precision. The collective model of Bohr and Mottelson, treating the nucleus as a deformable drop, captured rotational and vibrational behavior of nuclei far from closed shells. But the relationship between the two pictures — how collective deformation emerges from the single-particle degrees of freedom — remained murky.

The Interacting Boson Model (IBA), introduced by Akito Arima and Francesco Iachello in 1975 (Arima and Iachello, Physical Review Letters 35, 1069, 1975), offered an elegant bridge. The central idea was audacious in its simplicity: replace the fermion pairs of the shell model with bosons, and describe collective nuclear states using the algebra of these bosons.

The Boson Mapping: From Fermion Pairs to $s$ and $d$ Bosons

In the shell model, the low-energy structure of nuclei is dominated by paired valence nucleons. A pair of identical nucleons (two protons or two neutrons) in the same $j$-shell can couple to total angular momentum $J = 0, 2, 4, \ldots, (2j-1)$. The $J = 0$ pair is strongly favored by the short-range attractive pairing interaction, and the $J = 2$ pair represents the most important non-trivial coupling. Higher-$J$ pairs are less important because their overlap with the pairing force is weaker.

The IBA truncation retains only these two types of pairs:

• $s$ bosons ($J^\pi = 0^+$): represent $J = 0$ nucleon pairs (the pairing degree of freedom).
• $d$ bosons ($J^\pi = 2^+$): represent $J = 2$ nucleon pairs (the quadrupole degree of freedom).

Each boson carries the quantum numbers of a nucleon pair, and the total number of bosons is fixed:

$$N_B = N_s + N_d = \frac{1}{2}(N_{\text{valence protons}} + N_{\text{valence neutrons}}),$$

where "valence" means particles (or holes) counted from the nearest closed shell.

This mapping from a fermionic Hilbert space (vast, with antisymmetrization constraints) to a bosonic Hilbert space (much smaller, with only symmetrization constraints) is the source of the IBA's computational power. A shell-model calculation for a mid-shell rare-earth nucleus might require matrices of dimension $10^{10}$ or more; the corresponding IBA calculation involves matrices of dimension $\sim 10^2$.

The Algebra: U(6) and Its Subgroups

The six boson operators $\{s^\dagger, d_\mu^\dagger\}$ ($\mu = -2, -1, 0, 1, 2$) generate the unitary group $U(6)$. The 36 bilinear products $\{s^\dagger s, s^\dagger d_\mu, d_\mu^\dagger s, d_\mu^\dagger d_\nu\}$ form the generators of $U(6)$, and the IBA Hamiltonian is constructed from these generators.

The power of the algebraic approach lies in the fact that $U(6)$ has three maximal subgroup chains that each lead to the rotation group $O(3)$:

Chain I: The Vibrational Limit — $U(5)$

$$U(6) \supset U(5) \supset O(5) \supset O(3)$$

Quantum numbers: $N_B, \, n_d, \, \tau, \, n_\Delta, \, L, \, M$

The Casimir operators of $U(5)$ and its subgroups provide an analytic energy formula:

$$E(n_d, \tau, L) = \epsilon \, n_d + \alpha \, n_d(n_d + 4) + \beta \, \tau(\tau + 3) + \gamma \, L(L+1).$$

When $\alpha = \beta = \gamma = 0$, this reduces to equally spaced phonon multiplets: $E = \epsilon \, n_d$. The corrections lift the degeneracies within each $n_d$ multiplet. The U(5) spectrum closely resembles that of the geometric vibrational model (Section 8.3), with the crucial difference that the boson number $N_B$ imposes a natural cutoff — there are at most $N_B$ phonons.

Experimental realization: $^{110}$Cd ($Z = 48$, $N = 62$; $N_B = 7$) is one of the best U(5) examples, with a clear phonon-multiplet structure in its low-energy spectrum.

Chain II: The Rotational Limit — $SU(3)$

$$U(6) \supset SU(3) \supset O(3)$$

Quantum numbers: $N_B, \, (\lambda, \mu), \, K, \, L, \, M$

The energy formula involves the Casimir operator of $SU(3)$:

$$E(\lambda, \mu, L) = -\kappa [\lambda^2 + \mu^2 + \lambda\mu + 3(\lambda + \mu)] + \kappa' L(L+1).$$

The ground-state band has $(\lambda, \mu) = (2N_B, 0)$, with $L = 0, 2, 4, \ldots, 2N_B$. Within this band, $E \propto L(L+1)$ — exactly the rigid-rotor spectrum. The beta band has $(\lambda, \mu) = (2N_B - 4, 2)$, and the gamma band also emerges from specific $SU(3)$ representations.

The SU(3) limit generates not just the energy spectrum but also the $B(E2)$ transition rates, quadrupole moments, and two-nucleon transfer amplitudes — all from the algebra, with no numerical diagonalization.

Experimental realization: $^{156}$Gd ($Z = 64$, $N = 92$; $N_B = 12$) is a textbook SU(3) nucleus, with a ground-state band closely following $I(I+1)$ spacing and transition rates matching the SU(3) predictions.

Chain III: The $\gamma$-Unstable Limit — $O(6)$

$$U(6) \supset O(6) \supset O(5) \supset O(3)$$

Quantum numbers: $N_B, \, \sigma, \, \tau, \, n_\Delta, \, L, \, M$

The energy formula is:

$$E(\sigma, \tau, L) = A \, \sigma(\sigma + 4) + B \, \tau(\tau + 3) + C \, L(L+1).$$

This limit describes a nucleus that is deformed ($\beta \neq 0$) but has no preference for axial symmetry — the potential is flat in the $\gamma$ direction. The ground-state band has a characteristic energy ratio $E(4^+)/E(2^+) = 2.5$, and the spectrum shows a distinctive pattern of closely spaced $\gamma$-band states.

Experimental realization: $^{196}$Pt ($Z = 78$, $N = 118$; $N_B = 6$) is the classic O(6) example, identified by Cizewski et al. (1978) shortly after the IBA was proposed. The measured $R_{4/2} = 2.45$ is remarkably close to the O(6) prediction, and the level scheme matches the $O(6) \supset O(5) \supset O(3)$ classification in detail.

The Casten Triangle and Shape Phase Transitions

The three dynamical symmetries sit at the vertices of a triangle in the space of IBA Hamiltonian parameters. Any real nucleus lies somewhere within this triangle, with its position determined by the competition between pairing (which favors the $s$ boson condensate of U(5)), quadrupole deformation (which drives toward SU(3)), and $\gamma$-softness (which leads to O(6)).

The Casten triangle (R. F. Casten, 1985) is parameterized by two control parameters in the simplified "consistent-Q" Hamiltonian:

$$\hat{H}_{\text{CQF}} = c \left[ (1-\zeta) \hat{n}_d - \frac{\zeta}{4N_B} \hat{Q}^\chi \cdot \hat{Q}^\chi \right],$$

where: - $\zeta = 0$: the U(5) vertex (spherical vibrator) - $\zeta = 1$, $\chi = -\sqrt{7}/2$: the SU(3) vertex (prolate rotor) - $\zeta = 1$, $\chi = 0$: the O(6) vertex ($\gamma$-unstable)

As $\zeta$ increases from 0 to 1 along the U(5)–SU(3) edge, the nucleus undergoes a first-order quantum phase transition from spherical to deformed shape. This transition is characterized by:

1. Shape coexistence near the critical point — the potential energy surface has two competing minima (spherical and deformed).
2. A discontinuous jump in the equilibrium deformation $\beta_{\text{eq}}$ at the transition.
3. An abrupt change in the level scheme from vibrational to rotational character.
4. Large $E0$ (electric monopole) transition strengths between coexisting $0^+$ states.

The transition along the U(5)–O(6) edge is second-order (continuous), with the equilibrium deformation growing smoothly from zero.

Critical-Point Symmetries: E(5) and X(5)

In 2000–2001, Iachello identified exactly solvable models at the critical points of these phase transitions:

E(5) symmetry (Iachello, Physical Review Letters 85, 3580, 2000): describes the critical point of the second-order U(5)–O(6) transition. The potential in $\beta$ is an infinite square well, and the model is solvable in the five-dimensional $\beta$-$\gamma$ space using Bessel functions. Key predictions include:

• $R_{4/2} = 2.20$
• Specific ratios for excited $0^+$ state energies and $B(E2)$ values

The nucleus $^{134}$Ba was identified as the first empirical realization of E(5) symmetry, with level energies and transition rates matching the E(5) predictions to within 5–10%.

X(5) symmetry (Iachello, Physical Review Letters 87, 052502, 2001): describes the critical point of the first-order U(5)–SU(3) transition. The potential in $\beta$ is again an infinite square well, but the $\gamma$ degree of freedom is frozen at $\gamma = 0$ (axial symmetry). Key predictions include:

• $R_{4/2} = 2.91$
• A specific ratio $E(0^+_2)/E(2^+_1) = 5.65$
• Characteristic $B(E2)$ branching ratios

The nucleus $^{152}$Sm was identified as the first X(5) candidate (Casten and Zamfir, Physical Review Letters 87, 052503, 2001), with its measured $R_{4/2} = 3.01$ falling between the vibrational and rotational limits, close to the X(5) prediction.

The IBA-2 Extension: Protons and Neutrons

The original IBA (now called IBA-1) treats proton and neutron bosons as identical. The IBA-2 extension distinguishes between proton bosons ($s_\pi, d_\pi$) and neutron bosons ($s_\nu, d_\nu$), introducing the proton-neutron interaction as an explicit degree of freedom.

IBA-2 naturally explains several features that IBA-1 cannot:

• Mixed-symmetry states (states with proton and neutron deformations oscillating out of phase), identified through their strong $M1$ transitions.
• The $F$-spin quantum number, analogous to isospin but for bosons, which classifies states by their proton-neutron symmetry.
• Systematic trends across isotope chains, as the proton and neutron boson numbers change independently.

The lowest mixed-symmetry state, the $1^+$ "scissors mode" (two nuclear halves, proton and neutron, oscillating like scissors blades), was predicted by the IBA-2 and subsequently discovered in $^{156}$Gd by Bohle et al. (1984) at 3.075 MeV using inelastic electron scattering at the Darmstadt linear accelerator. The measured $M1$ strength of $B(M1) \approx 3 \mu_N^2$ confirmed the collective proton-neutron character of this excitation. This was a striking confirmation of the proton-neutron degree of freedom in the IBA framework and opened an entirely new field of nuclear spectroscopy — the study of mixed-symmetry states.

Subsequent experiments identified the scissors mode in dozens of deformed nuclei across the rare-earth and actinide regions, with a systematic dependence of its energy and $M1$ strength on deformation: $B(M1) \propto \beta_2^2$. This scaling is a direct prediction of the IBA-2, further validating the model.

A Worked Example: $^{150}$Nd in the IBA

To illustrate the practical application of the IBA, consider $^{150}$Nd ($Z = 60$, $N = 90$). Counting from the nearest closed shells ($Z = 50$, $N = 82$):

$$N_\pi = (60 - 50)/2 = 5, \quad N_\nu = (90 - 82)/2 = 4, \quad N_B = N_\pi + N_\nu = 9.$$

The measured $R_{4/2} = 2.93$ places this nucleus near the X(5) critical point. In the IBA, this corresponds to a Hamiltonian with $\zeta \approx 0.8$ and $\chi \approx -1.0$, intermediate between U(5) and SU(3) but closer to the deformed side.

A two-parameter IBA fit to the low-energy spectrum reproduces the energies of the ground-state band, beta band, and gamma band to within 5-10% — a remarkable achievement given that a full shell-model calculation for this nucleus would require a basis dimension of order $10^{12}$ or more. The IBA compresses this information into a model space of dimension $\sim 100$, making the physics transparent.

Modern Developments: Beyond the Standard IBA

The IBA framework continues to evolve. Recent developments include:

• The IBA with configuration mixing (IBM-CM): Extends the model to include intruder states — configurations with different numbers of active particles created by excitations across the shell gap. This is essential for understanding shape coexistence, where a nucleus simultaneously supports spherical and deformed structures at similar excitation energies. The lead isotopes around $^{186}$Pb, which show triple shape coexistence (spherical, oblate, and prolate minima within 700 keV of each other), have been successfully described within this framework.

• The proton-neutron IBA (IBA-2) with broken pairs: Includes configurations with broken nucleon pairs (quasiparticle excitations), extending the model to describe odd-mass and odd-odd nuclei through the Interacting Boson-Fermion Model (IBFM).

• Connection to density functional theory: Modern approaches use microscopic energy density functionals to compute the potential energy surface $V(\beta, \gamma)$, then map this surface onto the IBA Hamiltonian to determine the boson parameters from first principles. This eliminates the need for phenomenological parameter fitting and provides a bridge between the IBA and ab initio nuclear theory.

Assessment: Strengths and Limitations

Strengths of the IBA: - Provides analytic solutions at dynamical symmetry limits, offering deep physical insight. - Unifies vibrational, rotational, and $\gamma$-unstable behavior in a single framework. - Predicts quantum phase transitions with experimentally confirmed signatures. - Computationally tractable for any nucleus (boson space dimension $\sim 10^2$–$10^3$). - Natural connection to the shell model through the fermion-to-boson mapping.

Limitations: - The truncation to $s$ and $d$ bosons omits higher-multipole correlations ($g$ bosons with $J = 4$, needed for some observables). - The boson mapping is not unique — different prescriptions exist (Otsuka-Arima-Iachello, Ginocchio-Talmi), and the correspondence becomes less clean far from closed shells. - The model has several free parameters that must be fitted to data; it is a phenomenological model, not derived from first principles. - Intruder states (particle-hole excitations across shell gaps) require extensions beyond the standard IBA framework.

Despite these limitations, the IBA remains one of the most widely used models in nuclear structure physics, applied to hundreds of nuclei across the chart. Its algebraic elegance — the fact that the full taxonomy of collective nuclear behavior can be organized by the subgroup structure of $U(6)$ — represents one of the most beautiful applications of group theory in all of physics.

The Legacy of Arima and Iachello

The IBA earned Arima and Iachello numerous awards, including the Tom Bonner Prize in Nuclear Physics (1993) and the Wigner Medal (2003) for Iachello. The model inspired analogous algebraic approaches in molecular physics (the vibron model), hadronic physics, and even quantum chemistry.

More broadly, the IBA demonstrated that symmetry — not just as an exact invariance but as an organizational principle connecting different dynamical regimes — is the most powerful tool available for understanding quantum many-body systems. The Casten triangle, with its symmetry limits at the vertices and phase transitions along the edges, has become an icon of nuclear structure physics, reproduced in hundreds of papers and textbooks.

This case study develops the IBA framework introduced in Section 8.9 and provides context for the quantum phase transitions discussed in Section 8.9.4. The fermion-to-boson mapping connects to the pairing and deformed shell model topics of Chapter 9.