PHY401 - Nuclear and Particle Physics
Monsoon Semester 2020 Dr. Anosh Joseph, IISER Mohali
LECTURE 10
Monday, September 14, 2020 (Note: This is an online lecture due to COVID-19 interruption.)
Contents
1 Shell Model of the Nucleus I 1
1 Shell Model of the Nucleus I
The shell model of the nucleus was proposed to explain a collection of experimental data. The data showed that nuclei have nuclear shells in analogy to the atomic shell model. In 1963 Maria Mayer received the Nobel Prize in Physics for proposing the nuclear shell model of the atomic nucleus. Experiments showed that the binding energies deviated from the liquid drop model with increased binding at N or Z at the so-called magic numbers of 2, 8, 20, 28, 50, 126. The number of stable isotopes/isotones is significantly higher for nuclei with the proton/neutron number equal to the magic number. See Fig. 1. Nucleon capture cross sections are high for nuclei with one nucleon shy from the magic number (single vacancy in a closed shell), but significantly lower for nuclei with number of nucleons equal to the magic number (at the closed shell). See Fig. 2. Electric quadrupole moments should be zero for closed shell nuclei since they are spherically symmetric. Fig. 3 shows this to be so. The figure is based on measurements for odd A nuclei. The measured moments have been normalized with respect to the size and charge of each nucleus and these so called reduced quadrupole moments are plotted against the number of protons or neutrons - depending upon which is odd. It is also clear from the plot that in some cases the moments are relatively large. The quadrupole moments vanish near closed shells and reach their largest values far away from them. This points to some nuclei having shapes which are strongly non-spherical. The plot shows that prolate nuclei (Q > 0) are more common than oblately deformed ones (Q < 0). We could come up with a shell model similar to the shell model of the atom. PHY401 - Nuclear and Particle Physics Monsoon Semester 2020 ) MeV (
2 neutron 1 last
of 0 −1 energy −2 28 50 82 126 Binding 40 60 80 100 120 140 Neutron number N
Figure 1: Binding energy of the last neutron against neutron number. The measured binding energy is plotted relative to the value predicted by the semi-empirical mass formula. There are discontinuities of about 2 MeV at the shell closures.
We may assume that each nucleon occupies a well-defined energy level. A shell is described as the energy level where particles having the same energy exists. In this model, all the nucleons are paired one-to-one, neutron with a neutron, and proton with a proton. The electrons move in the atom in a central Coulomb potential emanating from the atomic nucleus. In the nucleon, on the other hand, the nucleons move inside a (mean field) potential produced by the other nucleons. That is, the motion of each nucleon is governed by the average attractive force of all the other nucleons. In both cases discrete energy levels arise, which are filled up according to the rules that obey the Pauli principle. Experimental data clearly pointed out to nuclear magic numbers at 2, 8, 20, 28, 50, 126. These are the magic numbers that show the most stable nuclei. See Fig. 4. The paired neutrons and protons in nuclear energy levels are filled when the number of neutrons or protons is equal to the magic number. The doubly magic nuclei, those with both magic proton and magic neutron numbers, are excep- tionally stable. These are the following nuclides
4 16 40 48 208 2He2, 8 O8, 20Ca20, 20Ca28, 82 Pb126
The unpaired nucleons are responsible for the properties of a nucleus. This is similar to the case in which valence electrons are responsible for different chemical properties of elements.
2 / 8 PHY401 - Nuclear and Particle Physics Monsoon Semester 2020
100 ) mb ( σ 10 section
Cross 1
28 50 82 126
20 40 60 80 100 120 Neutron number N
Figure 2: Nuclei with magic numbers of neutrons have neutron absorption cross-sections up to two orders of magnitude less than other nuclei with similar masses. This indicates that these nuclei are much less likely to absorb an additional neutron.
With the help of the shell model, we can accurately predict the properties of nuclei such as angular momentum and parity. But, for nuclei which are in a highly unstable state, the shell model needs to be modified or replaced with other models such as the collective model, liquid-drop model, and compound nucleus model. However, if we consider a nuclear shell model, with a flat bottom potential, we would get the shell gaps that explain the magic numbers at 2, 8, 20, 40, 70, 112. The first three magic numbers are in agreement with the data. But the consecutive higher ones were not. This led to a serious confusion back in the days. What might be happening? Let us look at the shell model with a flat bottom potential. The simplest of such potentials is the harmonic oscillator potential. Writing the wave function as
ψ(~r) = Rkl(r)Ylm(θ, φ), (1) where Rkl(r) is the radial wave function and Ylm are the spherical harmonics, we get the energy 1 2 2 levels, for the potential V (r) = 2 µω r
3 E = (2k + l + ) ω. (2) nl 2 ~
The energy is described by a single quantum number
n ≡ 2k + l. (3)
Since k is a nonnegative integer, for every even n, we have l = 0, 2, ··· , n − 2, n, and for every odd n we have l = 1, 3, ··· , n − 2, n.
3 / 8 PHY401 - Nuclear and Particle Physics Monsoon Semester 2020
167Er 0.30
176 0.25 Lu 10B 0.20 moment
0.15
55Mn 0.10 27Ai 9Be 115In quadrupole 2 H 79 0.05 41K Br 139La 201Hg 0.00 63Cu 209
Reduced Bi 93 −0.05 17O36Cl Nb 227Ac
8 20 28 123 82 126 −0.10 16 Sb 0 20 40 60 80 100 120 140 50 Number of odd neutrons or protons
Figure 3: Reduced electric quadrupole moments for nuclei with odd proton number Z or neutron number N plotted against this number. The solid curves are based upon the quadrupole moments of very many nuclei, of which only a few are explicitly shown here.
The magnetic quantum number m is an integer satisfying −l ≤ m ≤ +l. The degeneracy at level n is
X 1 D = (2l + 1) = (n + 1)(n + 2). (4) n 2 l
Let us use the spectroscopic notation for the orbital angular momentum l. That is, l = 0, 1, 2, ··· as s, p, d, f, g, ··· . In Fig. 4 we show a schematic representation of the shell structure of the nuclei. In the figure we have listed the number harmonic oscillator quanta in each set of shells. We have also made use of the fact that in the real nuclear potentials, the state with the same number of quanta but different l are no longer degenerate, but there are groups of shells with big energy gaps between them. This cannot predict the magic numbers beyond 20, and we need to find a different mechanism. Thus the early shell model was wrong but not completely wrong. They forgot to add an important piece.
4 / 8 PHY401 - Nuclear and Particle Physics Monsoon Semester 2020
126 1/2 3p 1/2 3/2 3p 13/2 13/2 5/2 3/2 5/2 2f 2f 7/2 9/2 9/2 7/2 1h 1h 82 3s 1/2 3s 1/2 3/2 3/2 2d 11/2 2d 11/2 5/2 7/2 7/2 5/2 1g 1g 50 9/2 9/2 1/2 1/2 2p 5/2 2p 5/2 3/2 1f 1f 3/2 26 7/2 7/2 20 3/2 3/2 1d 1/2 1d 1/2 2s 5/2 2s 5/2
8 1/2 1/2 1p 3/2 1p 3/2
1s 1/2 1s 1/2
Figure 4: Single particle energy levels calculated using Eq. (5). (See Ref. [2].) Magic numbers appear when the gaps between successive energy shells are particularly large. This diagram refers to the nucleons in the outermost shells.
This was the spin-orbit coupling term. This term captures the coupling between the orbital angular momentum ~l and the spin angular momentum ~s. The spin-orbit coupling term in nuclei reduces the energy of states with spin oriented parallel to the orbital angular momentum while increasing the energy of states with spin oriented opposite to the orbital angular momentum. There is one already known for atoms, which is the so-called spin-orbit splitting. 1 This means that the degeneracy in the total angular momentum (j = l± 2 ) is lifted by an energy term that splits the aligned from anti-aligned case. This is shown schematically in Fig. 4, where we label the states by n, l and j. The gaps between the groups of shells are in reality much larger than the spacing within one shell, making the binding-energy of a closed-shell nucleus much lower than that of its neighbors. Note that this is different from the result of spin-orbit interaction in atoms where states with spin oriented opposite to the orbital angular momentum are lower in energy. The spin-orbit interaction in atoms is understood from the electromagnetic interaction of the magnetic moment of an electron with the magnetic field resulting from orbiting a charged nucleus. The spin-orbit interaction in nuclei results from the spin-orbit part of the nuclear force.
5 / 8 PHY401 - Nuclear and Particle Physics Monsoon Semester 2020
Let us describe the potential by adding an additional term in the potential
h~l · ~si V (r) = Vcentral(r) + Vls(r) . (5) ~2 For very light nuclei A < 7, the potential can be approximated by that of a three-dimensional harmonic oscillator. That is, a flat bottom potential corresponding to a three-dimensional harmonic oscillator po- 1 2 2 tential V (r) = 2 µω r . But this potential grows significantly as the distance from the center r goes to infinity. A more realistic potential (especially for heavy nuclei) would be the so-called Woods-Saxon potential. This potential would approach a constant at this limit. This potential has the form 1 V (r) = −V . (6) 0 1 + e(r−R)/a
The central potential part Vcentral(r) can still be taken as the flat bottom well with diffuse surface (Woods-Saxon potential).
The term Vls(r) is often taken as negative derivative of V (r) with respect to r
dV (r) V (r) = − . (7) ls dr
For a flat bottom well with diffuse surface the derivative is peaked at the surface. The spin-orbit term makes the nuclear potential well wider for nucleons with spin parallel to the orbital angular momentum and less wide for nucleons with spin opposite to the orbital angular momentum. Wider well results in states of lower energies.
Coulomb repulsion adds to the proton well potential
Spin opposite to orbital
l = 0 Proton potential well Spin parallel to Neutron potential orbital well
Figure 5: Impact of the spin-orbit term.
Total angular momentum is a vector resulting from the coupling of the orbital and spin angular
6 / 8 PHY401 - Nuclear and Particle Physics Monsoon Semester 2020 momentum. We have ~j = ~l + ~s. (8)
The total angular momentum is quantized. We have the total angular momentum quantum number l + 1 if ~l parallel to ~s j = 2 (9) 1 ~ l − 2 if l antiparallel to ~s
Note that in the absence of the spin-orbit term the energy of a state does not depend on the total angular momentum ~j or its quantum number j. However, with the spin-orbit term it does. We can compute the magnitude of the ~l · ~s term in the following way. We use the fact that the square of the total, orbital, and spin angular momenta are defined by the corresponding quantum numbers
2 ~j2 = ~l + ~s = ~l2 + 2~l · ~s + ~s2 2~l · ~s = ~j2 − ~l2 − ~s2 1 =⇒ ~l · ~s = [j(j + 1) − l(l + 1) − s(s + 1)] . (10) 2
Since 1 ~l · ~s = [j(j + 1) − l(l + 1) − s(s + 1)] , (11) 2 and
l + 1 if ~l parallel to ~s j = 2 (12) 1 ~ l − 2 if l antiparallel to ~s then we must have
sl = 1 l if ~l parallel to ~s (j = l + 1 ) ~l · ~s = 2 2 (13) 1 ~ 1 −s(l + 1) = − 2 (l + 1) if l antiparallel to ~s (j = l − 2 )
This would amount to an energy splitting ∆Els which linearly increases with the angular mo- mentum as 2l + 1 ∆E = · hV (r)i. (14) ls 2 ls 1 Experimentally, it is found that Vls is negative, which means that the j = l + 2 is always below 1 the j = l − 2 level. This is in contrast to the atomic case, where the opposite occurs. The nuclear spin-orbit interaction has important consequences. The ls coupling in the atom generates the fine structure - these are small corrections of the order
7 / 8 PHY401 - Nuclear and Particle Physics Monsoon Semester 2020
α2. But the spin-orbit term in the nuclear potential leads to sizable splittings of the energy states which are indeed comparable with the gaps between the nl shells themselves. According to the spectroscopic notation we can use the orbital angular momentum quantum number to label the states. That is l = 0 for s states, l = 1 for p states, l = 2 for d states, and so on. After incorporating the spin-orbit coupling we have a modified spectroscopic notation. We add a subscript to the total momentum quantum number to indicate the coupling between the orbital and spin angular momentum.
Thus we have for l = 1, p 1 and p 3 states. 2 2 For l = 2, d 5 or d 3 states etc. 2 2 Note that for l = 0 (no orbital angular momentum) we have only s 1 since there cannot be 2 parallel or antiparallel coupling. Each state of total angular momentum j has 2j + 1 substates of the same energy, which differ by the azimuthal (magnetic) quantum number mj which runs from −j up to +j. Spin-orbit coupling causes the so-called intruder levels to drop down from the next higher shell into the structure of the previous shell. This can be seen in Fig. 4. We need to assign the parity quantum number P also to the single particle orbitals. The parity is a multiplicative quantum number. The total parity of a system is a product of the parities of all subsystems. The parity for a shell model state is defined by the product of the parities of 1. the radial wave function,
2. the spherical harmonic,
3. the spin wave function and,
4. the intrinsic parity of a nucleon. The parity of radial function, the spin function, and intrinsic parity of a nucleon are all positive. Thus, the parity of a single-particle orbital is fully defined by the parity of the spherical har- m l monics Yl and thus Pl = (−1) . The good thing is we can measure parities experimentally and they can be compared with model predictions.
References
[1] B. Povh, K. Rith, C. Scholz and F. Zetsche, Particles and Nuclei: An Introduction to the Physical Concepts, 6th edition, Springer (2008).
[2] P. F. A. Klinkenberg, “Tables of Nuclear Shell Structure,” Rev. Mod. Phys. 24, 63-73 (1952) doi:10.1103/RevModPhys.24.63
8 / 8