← Part IV: Nuclear Models
Chapter 14

Interacting Boson Model

s and d Bosons

The Interacting Boson Model (IBM), developed by Arima and Iachello, describes low-lying collective states of even-even nuclei in terms of bosons. Each boson represents a correlated pair of valence nucleons (protons or neutrons). The model uses two types of bosons:

  • s-boson: Angular momentum $\ell = 0$, representing a $J=0^+$ nucleon pair (1 component)
  • d-boson: Angular momentum $\ell = 2$, representing a $J=2^+$ nucleon pair (5 components: $m = -2, -1, 0, 1, 2$)

The total number of bosons $N$ is conserved and equals half the number of valence nucleon pairs (counting from the nearest closed shell). The boson creation and annihilation operators satisfy:

$$[s, s^\dagger] = 1, \qquad [d_m, d_{m'}^\dagger] = \delta_{mm'}, \qquad N = s^\dagger s + \sum_m d_m^\dagger d_m = n_s + n_d$$

The six operators $\{s^\dagger, d_m^\dagger\}$ span a 6-dimensional space. The bilinear products of creation and annihilation operators generate the Lie algebra of U(6), which is the starting group of the IBM.

U(6) Group Structure

The 36 generators of U(6) are the bilinear operators $\{s^\dagger s, s^\dagger \tilde{d}_m, d_m^\dagger s, d_m^\dagger \tilde{d}_{m'}\}$. The most general IBM-1 Hamiltonian conserving boson number and angular momentum can be written in terms of Casimir operators of various subgroups. The group chain starts with U(6) and must end with the rotational group:

$$\text{U}(6) \supset G \supset \text{O}(3) \supset \text{O}(2)$$

There are exactly three subgroup chains that yield analytical solutions (dynamical symmetries):

  • Chain I: $\text{U}(6) \supset \text{U}(5) \supset \text{O}(5) \supset \text{O}(3) \supset \text{O}(2)$
  • Chain II: $\text{U}(6) \supset \text{SU}(3) \supset \text{O}(3) \supset \text{O}(2)$
  • Chain III: $\text{U}(6) \supset \text{O}(6) \supset \text{O}(5) \supset \text{O}(3) \supset \text{O}(2)$

The general IBM Hamiltonian interpolates between these limits, forming a "symmetry triangle" (Casten triangle) in parameter space. Most real nuclei lie between the symmetry limits.

Casimir Operators & Energy Formulas

Each dynamical symmetry has an energy formula expressed in terms of Casimir operators of the subgroups in the chain. The Casimir operators are invariants that commute with all generators of their respective group.

U(5) Vibrational Limit

$$E(n_d, v, n_\Delta, L) = \epsilon\, n_d + \alpha\, n_d(n_d - 1) + \beta\, v(v+3) + \gamma\, L(L+1)$$

Quantum numbers: $n_d$ (d-boson number), $v$ (O(5) seniority),$L$ (angular momentum). This gives equally-spaced multiplets characteristic of a vibrational nucleus with $R_{4/2} = E(4^+_1)/E(2^+_1) \approx 2.0$.

SU(3) Rotational Limit

$$E(\lambda,\mu,K,L) = \kappa\left[\lambda^2 + \mu^2 + \lambda\mu + 3(\lambda+\mu)\right] + \kappa'\, L(L+1)$$

The quantum numbers $(\lambda,\mu)$ label the SU(3) irreducible representation. The ground-state band has $(\lambda,\mu) = (2N,0)$, giving $R_{4/2} = 10/3 \approx 3.33$, the rigid rotor value.

O(6) Gamma-unstable Limit

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

This limit describes gamma-unstable nuclei (free to oscillate in the $\gamma$ degree of freedom). The characteristic ratio is $R_{4/2} = 5/2 = 2.50$. The quantum number $\sigma$ labels O(6) representations, with $\sigma = N, N-2, \ldots$

E2 Transitions in IBM

The E2 transition operator in the IBM is:

$$T^{(E2)}_m = e_B\left[s^\dagger \tilde{d}_m + d_m^\dagger s + \chi\,(d^\dagger \times \tilde{d})^{(2)}_m\right]$$

where $e_B$ is the boson effective charge and $\chi$ is a structure parameter. The value of $\chi$ determines which dynamical symmetry limit is realized:

  • - $\chi = 0$: O(6) limit
  • - $\chi = -\sqrt{7}/2$: SU(3) limit (prolate)
  • - $\chi = +\sqrt{7}/2$: $\overline{\text{SU}(3)}$ limit (oblate)

The B(E2) values provide powerful tests of the model. Key predictions for the $2^+_1 \to 0^+_1$ transition:

$$B(E2; 2^+_1 \to 0^+_1) = e_B^2 \times \begin{cases} N & \text{U(5)} \\ \frac{N(2N+3)}{5} & \text{SU(3)} \\ \frac{N(N+4)}{5} & \text{O(6)} \end{cases}$$

These formulas show that B(E2) values scale linearly with N in the vibrational limit but quadratically in the rotational limit, reflecting enhanced collectivity in deformed nuclei.

Comparison with Experimental Data

The IBM has been remarkably successful in describing nuclear spectra. Key experimental signatures for identifying the dynamical symmetry of a nucleus:

ObservableU(5)SU(3)O(6)
$R_{4/2}$2.003.332.50
$R_{0_2/2_1}$2.00~5~3.3
B(E2) pattern$\propto N$$\propto N^2$$\propto N^2$

Classic examples of each symmetry include:

  • - U(5): $^{110}$Cd ($R_{4/2} = 2.14$), near-vibrational
  • - SU(3): $^{156}$Gd ($R_{4/2} = 3.24$), good rotor
  • - O(6): $^{196}$Pt ($R_{4/2} = 2.47$), gamma-unstable

IBM-2: Proton-Neutron Bosons

The IBM-2 distinguishes between proton ($\pi$) and neutron ($\nu$) bosons:

$$N = N_\pi + N_\nu, \qquad N_\pi = \text{proton boson number}, \quad N_\nu = \text{neutron boson number}$$

The IBM-2 Hamiltonian includes a crucial proton-neutron interaction term (the Majorana interaction) that controls the F-spin structure:

$$H_{\text{IBM-2}} = \epsilon_\pi\, \hat{n}_{d_\pi} + \epsilon_\nu\, \hat{n}_{d_\nu} + \kappa\, \hat{Q}_\pi \cdot \hat{Q}_\nu + \lambda\, \hat{M}_{\pi\nu}$$

where $\hat{Q}_\rho$ is the quadrupole operator for proton/neutron bosons and$\hat{M}_{\pi\nu}$ is the Majorana operator. IBM-2 introduces the concept of F-spin (analogous to isospin for bosons). States with maximum F-spin ($F = F_{\max} = N/2$) are fully symmetric under proton-neutron boson exchange and correspond to the low-lying collective states. Mixed-symmetry states with $F < F_{\max}$ include the scissors mode ($1^+$ state), first observed in $^{156}$Gd at ~3 MeV.

Python Simulation: IBM Energy Spectra

Energy level diagrams for the three dynamical symmetry limits of IBM-1, showing the characteristic spectral patterns and quantum numbers.

IBM-1 Dynamical Symmetries

Python

Energy spectra for U(5), SU(3), and O(6) limits

script.py139 lines

Click Run to execute the Python code

Code will be executed with Python 3 on the server

Fortran Implementation: IBM Calculator

Computes energy levels and quantum numbers for all three IBM-1 dynamical symmetry limits, including energy ratios and comparison across boson numbers.

IBM-1 Energy Level Calculator

Fortran

Analytical energy formulas for U(5), SU(3), and O(6) symmetry limits

program.f90140 lines

Click Run to execute the Fortran code

Code will be compiled with gfortran and executed on the server

Rate this chapter: