The Nilsson model is a nuclear shell model treating the atomic nucleus as a deformed sphere. In 1953, the first experimental examples were found of rotational bands in nuclei, with their energy levels following the same J(J+1) pattern of energies as in rotating molecules. Quantum mechanically, it is impossible to have a collective rotation of a sphere, so this implied that the shape of these nuclei was nonspherical. In principle, these rotational states could have been described as coherent superpositions of particle-hole excitations in the basis consisting of single-particle states of the spherical potential. But in reality, the description of these states in this manner is intractable, due to the large number of valence particles—and this intractability was even greater in the 1950s, when computing power was extremely rudimentary. For these reasons, Aage Bohr, Ben Mottelson, and Sven Gösta Nilsson constructed models in which the potential was deformed into an ellipsoidal shape. The first successful model of this type is the one now known as the Nilsson model. It is essentially a nuclear shell model using a harmonic oscillator potential, but with anisotropy added, so that the oscillator frequencies along the three Cartesian axes are not all the same. Typically the shape is a prolate ellipsoid, with the axis of symmetry taken to be z.
For an axially symmetric shape with the axis of symmetry being the z axis, the Hamiltonian is
H= | 1 |
2 |
2 | |
m\omega | |
z |
| ||||
z |
2 | |
m\omega | |
\perp |
(x2+y
2)-c | |
1\ell ⋅ |
2-\langle\ell | |
s-c | |
2(\ell |
2\rangle | |
N). |
Here m is the mass of the nucleon, N is the total number of harmonic oscillator quanta in the spherical basis,
\ell
\ell2
\ell(\ell+1)
2\rangle | |
\langle\ell | |
N=(1/2)N(N+3) |
\ell2
The anisotropy of the potential is such that the length of an equipotential along the z is greater than the length on the transverse axes in the ratio
\omega\perp/\omegaz
\delta ≈ 0.5
The mathematical details of the deformation parameters are as follows. Considering the success of the nuclear liquid drop model, in which the nucleus is taken to be an incompressible fluid, the harmonic oscillator frequencies are constrained so that
\omegaz\omega
3 | |
0 |
\hbar\omega0 ≈ (42 MeV)A-1/3
(\omega\perp/\omega
| |||||
\delta)/(1- | |||||
z) |
4 | |
3 |
\delta)
R=\omega\perp/\omegaz
\delta=(3/2)(R2-1)/(2R2+1)
The remaining two terms in the Hamiltonian do not relate to deformation and are present in the spherical shell model as well. The spin-orbit term represents the spin-orbit dependence of the strong nuclear force; it is much larger than, and has the opposite sign compared to, the special-relativistic spin-orbit splitting. The purpose of the
\ell2
2\rangle | |
-\langle\ell | |
N |
c1=2\kappa\hbar\omega0
c2=\mu\kappa\hbar\omega0
\hbar\omega0
For ease of computation using the computational resources of the 1950s, Nilsson used a basis consisting of eigenstates of the spherical hamiltonian. The Nilsson quantum numbers are
\{N,\ell,m\ell,ms\}
2Y | |
r | |
20 |
\delta
In an odd, well-deformed nucleus, the single-particle levels are filled up to the Fermi level, and the odd particle's Ω and parity give the spin and parity of the ground state.
Because the potential is not spherically symmetric, the single-particle states are not states of good angular momentum J. However, a Lagrange multiplier
-\omega ⋅ J
\langleJx\rangle
Often one wants to calculate a total energy as a function of deformation. Minima of this function are predicted equilibrium shapes. Adding the single-particle energies does not work for this purpose, partly because kinetic and potential terms are out of proportion by a factor of two, and partly because small errors in the energies accumulate in the sum. For this reason, such sums are usually renormalized using a procedure introduced by Strutinsky.
Single-particle levels can be shown in a "spaghetti plot," as functions of the deformation. A large gap between energy levels at zero deformation indicates a particle number at which there is a shell closure: the traditional "magic numbers." Any such gap, at a zero or nonzero deformation, indicates that when the Fermi level is at that height, the nucleus will be stable relative to the liquid drop model.