Nuclear Shell Model
It has been found that the nuclei with proton number or neutron number equal to
certain numbers 2,8,20,28,50,82 and 126 behave in a different manner when compared to
other nuclei having neighboring values of Z or N. Hence these numbers are known as magic
numbers. This is found to be in accordance with the observed nature of elements with filled
shells. Thus Physicists looked at such a possibility in case of filling of nucleons in the nucleus.
Thus a new model of nucleus has emerged. This model is known as the Shell model.
Experimental evidences for the existence of magic numbers;
1. The binding energy of magic numbered nuclei is much larger than the neighboring
nuclei. Thus larger energy is required to separate a single nucleon from such nuclei.
2. Number of stable nuclei with a given value of Z and N corresponding to the magic
number are much larger than the number of stable nuclei with neighboring values of Z
and N. For example, Sn with Z=50 has 10 stable isotopes, Ca with Z=20 has six stable
isotopes.
3. Naturally occurring isotopes whose nuclei contain magic numbered Z or N have greater
relative abundances. For example, Sr-88 with N=50, Ba-138 with N=82 and Ce-140 with
N=82 have relative abundances of 82.56%, 71.66% and 88.48% respectively.
4. Three naturally occurring radioactive series decay to the stable end product Pb with
Z=82 in three isotopic forms having N=126 for one of them.
5. Neutron absorbing cross section is very low for the nuclei having magic numbered
neutron number.
6. Nuclei with the value of N just one more than the magic number spontaneously emit a
neutron (when excited by preceding beta-decay) E.g., O-17, K-87 and Xe-137.
7. Nuclei with magic numbers of neutrons or protons have their first excited states at
higher energies than in cases of the neighboring nuclei.
8. Electric quadrupole moment of magic numbered nuclei is zero indicating the spherical
symmetry of nucleus for closed shells.
9. Energy of alpha or beta particles emitted by magic numbered radioactive nuclei is larger
than that from other nuclei.
Lecture points of Dr. H R Sreepad, Associate Professor, Govt College (Autonomous), Mandya - Page 1
Independent particle model
Consider a nucleon moving independently in the harmonic oscillator potential which is
spherically symmetric. The Schrodinger equation given below can be solved in the Cartesian
coordinate system as well as in the spherical coordinate system.
with -------(1)
In the Cartesian coordinate system, it can be separated into three linear harmonic oscillators and
consequently, the energy of the three-dimensional harmonic oscillator is the sum of the energies
of the three linear harmonic oscillators.
-----(2)
where N is the principal quantum number which will assume positive integral values including
zero. From Eq. (2), we observe that. the energy depends only on the quantum number N and not
on nx, ny and nz• So, there is a degeneracy of energy levels. Since the proton (or neutron) has spin
,
1/2ħ there are two possible spin orientations for each state specified by a set of quantum numbers
nx, ny and nz. Applying the Pauliexclusion principle, we obtain the number of protons or neutrons
that can have a particular energy as shown in Table I.
Table I: Number (NN) of protons or neutrons in each state with the principal quantum Number N.
The multiplicative factor 2 in the last column is due to two spin states.
Lecture points of Dr. H R Sreepad, Associate Professor, Govt College (Autonomous), Mandya - Page 2
In general
----(3)
However, it is found not convenient to work in Cartesian coordinate system. Since the potential
is spherically symmetric, one can attempt to solve the Schrodinger equation in spherical
coordinates.
Let us write 2 in spherical coordinates.
-------(4)
where L2 is the square of the angular momentum operator
----------------------------------(5)
which has the spherical harmonics Ylm(θ,φ) as eigen functions with eigen values l(l+1)ħ2.
-------------------------------(6)
Substituting (4) in (1) we get,
--------------(7)
Lecture points of Dr. H R Sreepad, Associate Professor, Govt College (Autonomous), Mandya - Page 3
Writing the solution of ψ(r) of eq. (7) as a product of radial and angular functions
-------------------------------(8)
And using the eigen value equation (6), we obtain the radial equation
-----------(9)
Applying the bound state boundary condition that the radial function should vanish at infinity,
we obtain the discrete energy levels of the three dimensional harmonic oscillator.
------------------------------------(10)
Comparing (10) with the energy levels (2) obtained in the Cartesian coordinate system,
--------------------(11)
One can find
--------------------(12)
The eigenvalues and eigenfunctions in spherical coordinate system depend on two quantum
numbers: the radial quantum number n and the orbital quantum number I. It follows that for a
given total quantum number N, the orbital quantum number I can take only even values if N is
even or only odd values if N is odd.
The additional quantum number m occurring in Eq. (6.8) is the magnetic quantum number. It
may be remarked that the accidental degeneracy in the states with a given I but with different m
having the same energy, occurs in any spherical potential. However, in the case of the oscillator
potential, there is a special degeneracy as well, i.e., states with different I values but with the
same N have the same energy. Accordingly, the number of protons or neutrons with a given
value of N is
------------------(13)
where the factor 2 is due to the two spin states. Equation (6.12) gives the different I values
allowed for a given N. The permitted values of I are all even or odd and consequently, the parity
(given by (-1)') of all the nucleons with a given quantum number N is the same.
The three dimensional harmonic oscillator energy levels are labelled by a pair of quantum
numbers (n, I). Also the spectroscopic notation s, p, d, f, g, h, .. , is used to denote states with I =
0,1,2,3,4,5, .. '.
Lecture points of Dr. H R Sreepad, Associate Professor, Govt College (Autonomous), Mandya - Page 4
(n, l) : Os; Op; (Od, Is); (Of, Ip); (Og, Id, 2s);
(Oh, If, 2p); (Oi, Ig, 2d, 3s); (OJ, Ih, 2f, 3p). ---------------(14)
The bracketed levels on the right hand side are degenerate. We also give below the normalized
radial functions for a few (n, l) values.
where
and
The radial quantum number nr denotes that it is the nth time that the I value occurs in the scheme.
Alternatively, nr can be interpreted as the number of nodes that occurs in the radial wave
Lecture points of Dr. H R Sreepad, Associate Professor, Govt College (Autonomous), Mandya - Page 5
function including the one at infinity. The bracketed levels on the right hand side of (6.20) are
degenerate. We shall use the notation (n, I) or (nr, I) as the situation demands.
In Table II, we present the number of protons or neutrons that can occupy each oscillator state
with total quantum number N.
Table II: Number of protons or neutrons occupying each state with total quantum number N and
the orbital substates with quantum numbers nr, I. The energy E is in units of ħω.
The numbers in the last column should correspond to the closed shell nuclei and hence to the
magic numbers but the discrepancy arises after the first three shells. Further, there is a
degeneracy of states with different I values. Such a degeneracy does not exist in the case of
square well potential. Empirically, let us assume a level split of different I orbitals with the
energy given by
assuming an additional force which lowers the states of larger I. A value D = -0.0225/ ħω is
found to be satisfactory.
However, the introduction of this additional force only removes the degeneracy with respect to I
but will not yield the correct magic numbers above 20. So, Mayer and Jenson independently in
1949 introduced a spin-orbit coupling l . 8. Originally, the spin-orbit force was introduced
phenomenologically to reproduce the magic numbers but now it is believed to arise from the
tensor force of the nucleon-nucleon interaction.
Lecture points of Dr. H R Sreepad, Associate Professor, Govt College (Autonomous), Mandya - Page 6
The single-particle energy level scheme with spin-orbit coupling is shown in Fig. I for the
nuclear shell model yielding correctly the magic numbers.
Lecture points of Dr. H R Sreepad, Associate Professor, Govt College (Autonomous), Mandya - Page 7
Fig 1:
Lecture points of Dr. H R Sreepad, Associate Professor, Govt College (Autonomous), Mandya - Page 8