Lecture 24 - YU
Transcript of Lecture 24 - YU
Phys. 441: Nuclear Physics 1Physics Department
Yarmouk University 21163 Irbid Jordan
Chapter 8 :γ Decay
Lecture 24
http://ctaps.yu.edu.jo/physics/Courses/Phys441/Lec8-1
© Dr. Nidal Ershaidat
� ������������ ������������� ��������� ������������������ ������
٢IntroductionThe gamma decay is the emission of an energetic photon when an excited state of a nucleus decays to a lower energy state. In general a series of γγγγ decays is necessary to reach the stable ground state.These decays are encountered each time a nucleus is excited. This excitation could be the result of α or βdecays or a nuclear reaction.
* λ(F) = E/hc = E(Mev)/1240
The resulting photons are just like the atomic X-rays in nature,i.e. they are electromagnetic radiations, but are more energetic (roughly 0.1 to 10 MeV, wavelengths are between 104 F and 100 F*)
� ������������ ������������� ��������� ������������������ ������
٣Importance of γ spectroscopy
γ photons are relatively easy to detect which makes them a very popular tool for spectroscopists.
The study of the competitive process to gamma decay namely, the internal conversion, is an excellent tool to obtain the spins and parities of nuclear states!
The Study of excited states is very rich in information about the nuclei properties (See chapter 3).
٤8-1 Energetics of γ DecayConsider a nucleus, at rest, in an excited state Ei decaying to a lower energy state Ef , a γ photon (Eγ=pγ c) is emitted: A* → A + γ
8-2γγγγ++++==== ppR
rr0
Eγ is the solution of the quadratic equation 8-38-4
∆∆∆∆++++±±±±−−−−====γγγγ
21
22
cME
211cME
Equation 8-1 can be written as:∆∆∆∆E = Ei - Ef = Eγγγγ + 8-32
2
cM2
E γγγγ
Conservation of momenta:
c
E γγ ==⇒ ppRγ−=⇒ ppR
rr
Conservation of energy:Ei = Ef + Eγ + TR 8-1
� ������������ ������������� ��������� ������������������ ������
٥
Eγ ≅ ∆E∆E is typically of the order of MeV. The rest energy Mc2 is of the order of A x 103 MeV, i.e. ∆E << Mc2
Expanding the square root we have:
8-5
∆∆∆∆++++∆∆∆∆++++++++−−−−====γγγγ
2
222
cME
21
cME
11cME
From equation 8-3, if we neglect the term ∆E/Mc2, we have Eγ = ∆E.The correction to Eγ due to the recoil energy [(∆E)2/2Mc2] is negligible (10-5) and except a special case in nuclear spectroscopy*, Eγ is simply taken as equal to ∆E
* Mössbauer Spectroscopy which is dedicated to the use of this correction
( )2
2
cM2
EE
∆−∆≈
٦8-2) Generalized EM Moments
For a multipole electromagnetic moment of order L the corresponding em radiation, the moment is said to be an 2L-pole
The electromagnetic theory gives the following properties of a 2L-pole radiation
Index of the radiation1 Dipole2 Quadrupole3 Sextupole
Moment
The γ photons emitted in a nuclear desexcitation are electromagnetic radiations. This means that an em interaction isbehind this process! An em field produced by oscillations of charges and currents produces energy. In QM this energy is emitted as quanta (photons) of energy E = h ν
� ������������ ������������� ��������� ������������������ ������
٧Properties of a multipole radiationThe classical electromagnetic theory gives the following properties of a 2L-pole radiation1) The angular distribution of 2L-pole radiation is governed by the
Legendre Polynomials P2L=(cosθ).For a dipole (L=1) P2 ∝ (3cos2θ – 1) andFor a quadrupole (L=2) P4 ∝ (35 cos4θ – 30 cos2θ +3)2) The parity of the radiation field is:
(-1)L+1 for an electric multipole and (-1)L for a magnetic multipole.
For an electric dipole parity is +1 and for a magnetic dipole parity is -1
3) The radiation power is
( ) ( )( )[ ]
( )[ ]22L2
20
Lmc!!1L2L
c1L2LP σ
ω
+ε+=σ
+
γ 8-6
� ������������ ������������� ��������� ������������������ ������
٨Generalized multipole moment
In the previous expression, m(σL) is the amplitude of the time varying electric or magnetic moment.
This moment for L=1 differs from the electric dipole moment p and the magnetic dipole moment µ through some relatively unimportant numerical factors of order unity.
( ) ( )( )[ ]
( )[ ]22L2
20
Lmc!!1L2L
c1L2LP σ
ω
+ε+=σ
+
γ *
* (2L+1)!! Is a double factorial and given by (2L-1)!!= (2L+1)x(2L-1) x(2L-3)…3x1
σ = E or M to represent electric or magnetic radiation respectively
� ������������ ������������� ��������� ������������������ ������
٩Multipole moment in Quantum Mechanics
Where we replaced the classical m(σL) by the matrix element mfi(σL) associated to the (quantum) operator m(σσσσL), defined by:
The classical properties we have seen are the same in quantum mechanics. We discuss here how the expression of the radiated power (8-6) is written in QM
(((( )))) (((( )))) (((( ))))(((( ))))[[[[ ]]]] (((( ))))[[[[ ]]]]2
1L2
20
Lmc!!1L2L
1L2LPL σσσσ
ωωωω
++++εεεε++++====
ωωωωσσσσ====σσσσλλλλ
++++
fihh
The decay probability per unit time (the decay constant) for the emission of photon of energy ω is given by:h
(((( )))) (((( )))) vdLmLmoperator
*∫∫∫∫ ψψψψσσσσψψψψ====σσσσ iffi 321 8-8
8-7
� ������������ ������������� ��������� ������������������ ������
١٠Estimation of λ(σL) – Weisskopf estimates
- The Electric Transitions Case
In order to compute λ(σL) we need to know the initial and final wave functions, i.e. ψi and ψf
Let’s consider that the transition is due to a single proton which changes from one shell-model state to another
For L = 1 (dipole) radiationThe operator m(EL) becomes ez
The multipole operator m(EL) contains a term of the form e rL Ylm(θ,φ)
For L = 2 (quadrupole) radiationThe operator m(EL) becomes e(3z2 – r2)
� ������������ ������������� ��������� ������������������ ������
١١
Approximation of ψi and ψfRadial part of the (electric) transition probabilityIf we take the radial parts of ψi and ψf to be constant for r<R (the nucleus radius) and to be = 0 for r > R then the Radial part of the transition probability is of the form:
L
3
3L
R
0
2
R
0
L2
R3L
3
R31
R3L
1
drr
drrr
++++====++++====
++++
∫∫∫∫
∫∫∫∫8-9
Angular part of the (electric) transition probabilityWe can reasonably take the angular integrals to be equal to 1 and the EL transition probability is estimated to be:
(((( )))) (((( ))))(((( ))))[[[[ ]]]]
L2
21L2
0
2
2 Rc3L
3c
Ec4
e!!1L2L1L8
LE
++++
πεπεπεπε++++
++++ππππ≅≅≅≅λλλλ++++
hh8-10
١٢λ(E1), λ(E2), λ(E3), …Taking R = R0 A1/3 with R0 = 1.2 F, the lower multipole orders are given in the table below.
3
0
2
1031.7FMeV197FMeV44.1
c4e −−−−××××========
πεπεπεπε h
(((( )))) (((( ))))[[[[ ]]]] (((( )))) 33222323
32 EA2.1103
43
1971
1031.7!!3
161E ××××
××××
ππππ≅≅≅≅λλλλ −−−−
L λ(EL)1 λ(E1) = 1.0x1014 A2/3 E3
2 λ(E2) = 7.3x107 A4/3 E5
3 λ(E3) = 34 A2 E7
4 λ(E4) = 1.1x10-5 A2 E7
λ is in s-1 when E is expressed in MeV
8-11
� ������������ ������������� ��������� ������������������ ������
١٣
λ(M1), λ(M2), λ(M3), …A similar treatment for the magnetic transitions gives for λ(ML) the following expression:
(((( )))) (((( ))))(((( ))))[[[[ ]]]]
2L2
21L2
0
22
2p
2
p2
Rc2L
3c
E
c4e
cmc
1L1
!!1L2L1L8
LM
−−−−++++
++++
××××
πεπεπεπε
++++
−−−−µµµµ++++
++++ππππ≅≅≅≅λλλλ
h
h
h
8-12
We replace the factor by 10 where µp is the nuclear
magnetic moment of the proton on which depends the moment
operator and this gives the following values for s λ in s-1 when E is
expressed in MeV
2
p 1L1
++++
−−−−µµµµ
١٤
Estimation of λ(M1), λ(M2), λ(M3), …L λ(ML)
1 λ(M1) = 5.6x1013 E3
2 λ(M2) = 3.5x107 A2/3 E5
3 λ(M3) = 16 A4/3 E7
4 λ(M4) = 4.5x10-6 A2 E9
See discussion in Krane p. 332
8-13
λ is in s-1 when E is expressed in MeV
� ������������ ������������� ��������� ������������������ ������
١٥λ(EL) and λ(ML)L λ(EL)1 λ(E1) = 1.0x1014 A2/3 E3
2 λ(E2) = 7.3x107 A4/3 E5
3 λ(E3) = 34 A2 E7
4 λ(E4) = 1.1x10-5 A2 E7
λ(ML)
λ(M1) = 5.6x1013 E3
λ(M2) = 3.5x107 A2/3 E5
λ(M3) = 16 A4/3 E7
λ(M4) = 4.5x10-6 A2 E9
Two major conclusions:1) The lower multipolarities are dominant. For L > 3 the decay rates of a
γ transition become 10-5 smaller and indicate these transitions due to an L-pole > 4 have infinitesimal probabilities.
2) For a given multipole order, electric transitions are twice as big as the magnetic transition in medium and heavy nuclei (A > 70)
١٦
Weisskopf estimates vs. ExperimentAlthough Weisskopf estimates are not true theoretical calculations, they provide us with a good tool to compare transition probabilitiesIf the observed transition rates are many order of magnitude smaller than the Weisskopf estimate then we “suspect” that a poor match between the initial and final wave function is slowing the transition. On the other hand, if the observed transition rates were much greater than the Weisskopf estimate then we might “guess”that more than one nucleon is involved in the transition.
See next lecture for more details
Phys. 441: Nuclear Physics 1Physics Department
Yarmouk University 21163 Irbid Jordan8-2 : Lifetimes for γ Emission,8-3 : Selection Rules
Lecture 25
http://ctaps.yu.edu.jo/physics/Courses/Phys441/Lec8-1
© Dr. Nidal Ershaidat
8-2 : Lifetimes for γ Emission
١٩Weisskopf Estimates vs. ExperimentExample: Fig. 8-1 shows the energy levels of the (even-even) 72Se (Z=34) isotope
1) Evaluation of the partial decay rate for γγγγ emission
All details are shown in an energy level scheme: spin-parity, energies and/or γ transition energies and half-lives of excited states.
Energies and γ transition energies are given in keV
� ������������ ������������� ��������� ������������������ ������
٢٠Example: 72Se 3rd excited stateWe’ll have a closer look on the 3rd excited state of 72Se.Spin-parity = 2+
The decay constant λt = ln2/t1/2= 8.0x1010 s-1
Energy = 1317 keVt1/2 = 8.7x10-12 s
Transition 3rd-1st = 455 keVTransition 3rd-2nd = 380 keV
Neglecting the (internal) conversion factors λt is simply the sum of the decay rates of the three transition that depopulate this excited state!, i.e.: λt = λγ,1317 + λγ,445 + λγ,380
The measured relative intensities(or branching ratios) are λγ,1317: λγ,445:λγ,380 = 51:39:10
� ������������ ������������� ��������� ������������������ ������
٢١Comparison with Weisskopf Estimates
λγ,1317 = 0.51 * λt = 4.1x1010 s-1
λγ,455 = 0.39 * λt = 3.1x1010 s-1
λγ,380 = 0.10 * λt = 0.8x1010 s-1
The partial decay rates of the 3 transitions (1317,455,380) are:
The following tables give the calculations of λ(EL) and λ(ML) for the energies involved in this example (A=72)
٢٢
λ(EL) (s-1)L
1 λ(E1) = 3.95x1015
2
3 λ(E3) = 1.20x106
λ(ML) (s-1)
λ(M1) = 1.28x1014
λ(M2) = 2.40x109
λ(M3) = 3.30x104
λγ,1317 = 4.1x1010 s-1
Energy Level 1.317 MeV
λ(E2) = 8.70x1010
Weisskopf Estimates for E = 455 keV
Weisskopf Estimates for E = 1317 keV
λ(EL) (s-1)L
1 λ(E1) = 1.63x1014
2
3 λ(E3) = 7.11x102
λ(ML) (s-1)
λ(M1) = 5.28x1012
λ(M2) = 1.18x107
λ(M3) = 1.94x101
λγ,455 = 3.1x1010 s-1
Energy Level 0.455 MeV
λ(E2) = 4.26x108
X
� ������������ ������������� ��������� ������������������ ������
٢٣Weisskopf Estimates for E = 380 keVλ(EL) (s-1)L
1 λ(E1) = 9.50x1013
2
3 λ(E3) = 2.02x102
λ(ML) (s-1)
λ(M1) = 3.07x1012
λ(M2) = 4.80x106
λ(M3) = 5.48
λγ,380 = 0.8x1010 s-1
Energy Level 0.380 MeV
λ(E2) = 1.73x108
The previous calculations indicate that the favored transitions are the E2 ones. But they also show that the measured values are one order of magnitude greater than Weisskopf Estimate. There is a strong evidence for the collective structure of the nucleus!, since Weisskopf used the shell (a single individual particle) model tomake his estimations.
� ������������ ������������� ��������� ������������������ ������
٢٤The M4 Transitions (Systematics) Case
This figure shows in particular the good agreement with the expected E-9 dependence.
Fig. 8-2 represents the experimental data for different nuclei.
The straight line represents the Weisskopf estimate (See previous lecture, Relations 8-13,slide # 14)
τ(M4) = 1.54x105 A2 E-9
Fig. 8-2 : log(τ A2) vs. E(in KeV)
8-3 : Selection Rules
� ������������ ������������� ��������� ������������������ ������
٢٦Multipoles and Angular MomentaAn em field produced by oscillations of charges and currents produces also angular momentum. In QM this angular momentum carried by the quanta of energy (photons of energy E = h ν) is quantized.
The rate at which this angular momentum is radiated, is proportional to the rate at which energy is radiated.
The proportionality is preserved if each emitted photon carries a definite angular momentum.
Conclusion: a multipole of order L transfers an angular a multipole of order L transfers an angular momentum of L per photon.momentum of L per photon.
A multipole operator of order L includes a spherical harmonic Ylm(θ,φ), which is associated with an angular momentum L.
٢٧Angular Momentum and Parity Selection Rules
Example: For Ii = 3/2 and If = 5/2, the possible values for L are:1, 2, 3 and 4 and the radiated field would be a mixture of dipole, quadrupole, octupole and hexadecapole radiation!
Thanks to the rules of addition of angular momenta we know that L could only have restricted values.
8-15| Ii – If | ≤ L ≤ Ii + If
The relative parity of the initial and final levels determine the type of the emitted radiation (electric or magnetic). The following table resumes the parities related to em radiations
Consider a γ transition from an initial excited state of angular momentum Ii and parity πi to a final state (If, πf). Assume Ii ≠ If.(Spin-parity for these states are and respectively)iππππ
iI fππππfI
Conservation of angular momentum is expressed by:LIIrrr
++++==== fi 8-14
٢٨Parity and EM TransitionsThe following table resumes what we know about the parity associated to electric and magnetic transition.
8-16L Electric Transition Magnetic TransitionEven + -Odd - +
Parity
Ii+
If L
+
∆π
∆π = no- -
The following notations are used when studying parity changes –see Table 8-17
++ Even
-- Odd+ ∆π = yes-
- +
8-17
The two tables are used to determine the type of the emitted radiation (electric or magnetic).
٢٩ExampleLet’s Take again the previous example Ii = 3/2 and If = 5/2.
8-18
L Transition∆π = no
2 E23 M3
E44
1 M1L Transition
∆π = yes
2 M23 E3
M44
1 E1
For L = 1 and ∆π = no, the transition cannot be electric because the associated parity is (-)1 negative and this would give a final parity different from the initial one. In this case the transition is the magnetic dipole M1 transition.In the case ∆π = yes, the transition cannot be magnetic; the associated parity is (-)1+1 positive. In this case the transition is the electric dipole E1 transition. See Table 8-18 for the other cases.
٣٠
The Case Ii = If.
Selection RulesThe precedent example suggests restrictions on the possible transitions. These restrictions are called selection rulesselection rules
∆π = no even electric, odd magnetic
| Ii – If | ≤ L ≤ Ii + If ; Ii ≠ If.8-19
∆π = yes odd electric, even magnetic
This exception to the previous selection rules occurs because itwould mean that L = 0 is a possible value, and there are no monopole transitions in which a single photon is emitted!*
* The magnetic monopole does not exist. For an electric monopole, a spherical distribution of charge (like a single point charge) the Coulomb field is not affected by radial oscillations and thus no corresponding radiation is produced.
The lowest multipole order allowed in the case where Ii ≠ If. is L = 1
٣١Pure Multipole TransitionAnother interesting case is when (Ii≠0, If.=0) or (Ii=0, If. ≠ 0)L is equal to Ii for the first type of transition (or If for the second one).For an even Z – even N nucleus (like 72Se) the first excited state 2+
decays to the ground state 0+ through the emission of a pure E2 (quadrupole) transition.
� ������������ ������������� ��������� ������������������ �����
٣٢Pure Multipole TransitionAnother interesting case is when (Ii≠0, If.=0) or (Ii=0, If. ≠ 0)L is equal to Ii for the first type of transition (or If for the second one).For an even Z – even N nucleus (like 72Se) the first excited state 2+
decays to the ground state 0+ through the emission of a pure E2 (quadrupole) transition.
Exercise: Find the type of transitions for the decay from the fourth state to the ground state in the case of 72Se
� ������������ ������������� ��������� ������������������ �����
٣٣
Internal Conversion
The Ii = If = 0 Case – Internal Conversion
In this process, the available energy is transmitted to an orbital electron which in turn is ejected. Orbital electrons with wave functions penetrating the nucleus field are concerned.
Here the only possible value for L is zero. This case is not permitted for radiative transitions.The transition between 2 (excited) states with spin = 0* occurs through the competitive process we mentioned before, namely the internal conversion.
* A few even Z – even N nuclei have 0+ first excited states. Those are forbidden to decay to the 0+ ground state by γ emission.
٣٤
λ(σL, A=125,E=1 MeV)
λ(E2) = 1.825x109 s-1
λ(M3) = 1.0x104
λ(E4) = 1.72 s-1
λ(σL)/λ(M1)
11.4x10-3
≅2.1x10-10
1.3x10-13
Selection Rules and Weisskopf EstimatesIn table 8-18 we showed which transitions are possible between the 2 states Ii = 3/2 and If = 5/2. Several multipoles are permitted. For example in the ∆π = no case, M1, E2, M3 and E4 are allowed.To decide which are the ones we observe, we shall make a simple calculation using Weisskopf estimates.Let’s assume a medium-weight nucleus (A=125)* and E = 1 MeV. According to Weisskopf, the transition probabilities
λ(M1) = 5.6x1013 s-1 L Transition∆π = no
2 E23 M3
E44
1 M1
* A2/3 = 25, A4/3 = 625
٣٥
λ(σL, A=125,E=1 MeV)
λ(M2) = 4.375x109 s-1
λ(E3) = 5.31x104
λ(M4) = 0.70 s-1
λ(σL)/λ(E1)
12.3x10-7
≅2.1x10-10
2.1x10-17
Selection Rules and Weisskopf EstimatesThe previous results show that the lower orders are dominant. They also indicate that this transition could be composed of M1 radiation with possibly a small mixture of E2.For the ∆π = yes case, the calculation show that the E1 is dominant and the other modes are most likely not to occur!
λ(E1) = 1.25x1016 s-1 L Transition∆π = yes
2 M23 E3
M44
1 E1
٣٦Expectations based on the single-particle Estimates1) The lowest permitted multipoles usually dominate
3) Emission of multipole L+1 is less probable than emission of multipole L by a factor of the order of 10-5
2) Electric multipole emission is 2 orders of magnitude more probable than the same magnetic multipole emission.
4) Points 2 and 3 combined give the following relations (L’ = L +1)(((( ))))(((( ))))
(((( ))))(((( ))))
(((( ))))(((( ))))
325 101010MLEL
EL'EL
ML'EL −−−−−−−− ====××××====
λλλλλλλλ××××
λλλλλλλλ====
λλλλλλλλ
8-20(((( ))))(((( ))))
(((( ))))(((( ))))
(((( ))))(((( ))))
725 101010ELML
ML'ML
EL'ML −−−−−−−−−−−− ====××××====
λλλλλλλλ××××
λλλλλλλλ====
λλλλλλλλ
Important RemarkImportant Remark : The precedent calculations are based on a single particle model with simple approximations. We observe in the lab. transitions in which λ(E1) > λ(M1) especially in transitions between vibrational and rotational collective states