KEK Internal 2000-10
November 2000
R
Lecture Note on
Photon interactions and Cross Sections
H. Hirayama
Lecture Note on
Photon Interactions and Cross Sections
Hideo Hirayama
KEK, High Energy Accelerator Research Organization,
1-1, Oho, Tsukuba, Ibaraki, 305-0801 Japan
Abstract
This is the lecture note used at the \Tutorials for Electron/Gamma Monte Carlo: building
blo cks and applications" as part of the International Conference on the Monte Carlo 2000,
Advanced Monte Carlo on Radiation Physics, Particle Transp ort Simulation and Applications
(Monte Carlo 2000) held at Lisb on Portugal, 23-26 Octob er, 2000.
1 Intro duction
This lecture contains descriptions of the individual atomic cross sections of photons for scat-
tering, absorption, and pair pro duction that are necessary for photon transp ort calculations,
together with related cross section data. Many parts on this lecture are based on the series of
works by J. H. Hubb ell[1, 2, 3, 4, 5, 6] in this eld.
2 Classi cation of Interaction
The interaction of photons with matter may b e classi ed according to:
1. The kind of target, e.g., electrons, atoms, or nuclei, with which the photons interact, and
2. The typ e of event, e.g., scattering, absorption, pair pro duction, etc:, whichtakes place.
Possible interactions are summarized in Table 1.[7 ] As an example, the cross sections of a
photon with Cu is shown in Fig. 1[8 ]. It is apparent from Fig. 1 that for the attenuation of
photons the most imp ortantinteractions are:
1. The photo electric e ect, ,
pe
2. Compton scattering, , and
C
3. Electron-p ositron pair pro duction, ( + ).
pair tr ip
Rayleigh scattering, , is usually of minor imp ortance for the broad b eam conditions typi-
R
cally found in shielding, but must b e known for the interaction of attenuation co ecient data.
The photonuclear e ect, , is mostly restricted to the region of the giant resonance around
ph:n
10 to 30 MeV where, at the resonance p eak, it may amount to as much as 10 p ercent of the
total \electronic" cross section.
Elastic nuclear scattering, inelastic nuclear scattering and Delbruc k scattering are negligible
pro cesses in photon interactions. Elastic nuclear scattering is regard as a nuclear analog to very
low energy Compton scattering by an electron. In this pro cess, a photon interacts with a nucleon
in such a manner that is a photon re-emitted with the same energy.
During inelastic nuclear scattering, the nucleus is raised to an excited level by absorbing a
photon. The excited nucleus subsequently de-excites by emitting a photon of equal or lower
energy.
The phenomenon of photon scattering by the Coulomb eldofanucleus is called Delbruc k
scattering (also called nuclear p otential scattering). It can b e thought of as virtual pair pro duc-
tion in the eld of the nucleus { that is, pair pro duction followed by annihilation of the created
pair. 1 4 10 Cu
) Total 2 2
cm 10
-24 Pair production (nucleus) (10
0 10 Compton Rayleigh scattering (γ,n) Pair production (barn/atom) (electron) σ -2 10 Photoelectric Effect
Pion Quasi-deutron production disintegration 10-4 0.01 0.1 1 10 100 1000 104 105
Photon Energy (MeV)
Figure 1: Cross section of Cu for photons b etween 10 keV and 100 GeV.[8 ]. 2
Table 1. Classi cation of elementary photon interactions.
Typ e of Scattering
interaction Absorption Elastic Inelastic
Interaction (Coherent) (Incoherent)
with:
Atomic Photo electric Rayleigh Compton
electrons e ect scattering scattering
(
2 4
Z Z (L:E :)
R
Z
pe C
5
(L:E :) Z (H:E :)
Nucleus Photonuclear Elastic Inelastic
reactions nuclear nuclear
( ,n),( ,p), scattering scattering
2 0
photo ssion, etc. ( ; ) Z ( ; )
Z
ph:n:
(h 10MeV)
Electric eld Electron-p ositron Delbruc k
surrounding pair pro duction in scattering
charged eld of nucleus,
2
particles Z
pair
(h 1:02MeV)
Electron-p ositron
pair pro duction in
electron eld,
2
Z
tr ip
(h 2:04MeV)
Nucleon-antinucleon
pair pro duction
(h 3 GeV)
Mesons Photomeson Mo di ed
pro duction, ( ; )
(h 150 MeV)
3 Photo electric E ect
Studies related to the photo electric e ect are reviewed historically by Hubb ell in NSRDS-NBS
29[1 ].
In the atomic photo e ect, a photon disapp ears and an electron is ejected from an atom.
The electron carries away all of the energy of the absorb ed photon, minus the energy binding
the electron to the atom. The K -shell electrons are the most tightly b ound, and are the most
imp ortant contributions to the atomic photo e ect cross-section in most cases. However, if the
photon energy drops b elow the binding energy of a given shell, an electron from that shell cannot
b e ejected. Hence, particularly for medium- and high-Z elements, a plot of versus the photon
pe
energy exhibits the characteristic sawto oth absorption edges, since the binding energy of each
electron subshell is attained and this pro cess is p ermitted to o ccur. These feature can b e seen
well in Fig. 2, which shows the total photo electric cross section, , for Ag together with
pe
subshell ones[9 ]. 3 6 10 L -edge Ag 3 L -edge 2 L -edge 1 Total 105 K (barn/atom) L1 L2 4 K-edge L3 10 M section
103 cross
102
1
Photoelectric 10 1 10 100
Photon Energy (keV)
Figure 2: Total and partial atomic photo e ect of Ag.[9 ].
The order of magnitude of the photo electric atomic-absorption cross section is
(
4 3
Z =(h ) low energy
(1)
pe
5
Z =h high energy :
Dramatic resonance structures of the order of 10% just ab ove the absorption edges are
well known and can be easily observed with high-resolution sp ectrometers. Owing to their
dep endence on the temp erature, chemical binding and other variable atomic environments, this
extended X-ray absorption ne structure (EXAFS) is not included in the photon transp ort
calculation. However, EXAFS can b e a ma jor analytical to ol in X-ray di raction sp ectrometry.
3.1 Relaxation pro cesses after the photo electric e ect
Avacancy created by the ejection of an electron from an inner shells is lled by an outer electron
falling into it (de-excitation); this pro cess may b e accompanied by one of the following mo des:
1. A uorescence X-ray is emitted from the atom, with a photon energy equal to the di erence
between the vacancy-site inner-shell energy level and energy level of the particular outer
shell which happ ens to supply the electron to ll the vacancy ( uorescence yield, ! ,isthe
fraction of uorescence X-ray emission).
2. The excess energy ejects an outer-shell electron from the atom. This electron is known as
an Auger electron (Auger yield, a, is the fraction of Auger emission).
3. A vacancy is lled by an electron in a higher subshell, like from the L subshell to the
2
L subshell. This pro cess is called Coster-Kronig. As a result, a new vacancy is created,
1
in which is lled by one of the mo des (Coster-Kronig yield, f , is the fraction of Coster-
Kronig). 4
The sum of uorescence yield, ! , Auger yield, a, and Coster-Kronig yield, f , is unity:
! + a + f =1 (2)
Exp erimental and theoretical uorescence-yield information have b een reviewed by Fink et
al.[10 ], Bambynek et al.[11 ], Krause[12 ] and Hubb ell[2].
The yields of uorescence, Auger electrons and Coster-Kronig after the K-, L -, L - or L -
1 2 3
photo electric e ect[13 ] are shown in Fig. 3.
1 L1-Fluorescence 1.2 L2-L1 Coster-Kronig 0.8 L3-L1 Coster-Kronig L1 Auger 1 0.6 0.8
0.6 K-Fluorescence Yields 0.4 K-Auger Yields 0.4 0.2 0.2
0 0 0 20406080100 0 20406080100 Atomic Number, Z Atomic Number, Z
1 1
0.8 0.8
0.6 0.6 L3 Fluorescence L2 Fluorescence L3 Auger
Yields 0.4 L3-L2 Coster-Kronig Yields 0.4 L2 Auger
0.2 0.2
0 0 0 20406080100 0 20406080100
Atomic Number, Z Atomic Number, Z
Figure 3: Yields of uorescence, Auger electrons and Coster-Kronig after the K-, L -, L - or
1 2
L -photo electric e ect[13 ].
3
The uorescence X-ray intensities are imp ortant for low energy photon transp ort. Sco eld
published K and L uorescence X-ray intensities calculated with the relativistic Hatree-Slater
theory for elements with Z=5 to 104[14 ]. His data are cited in a Table of Isotop es Eighth
Edition[13], but are slightly di erent with the exp erimental results, as presented by Salem et
al[15 ]. 5
4 Scattering
4.1 Compton scattering, Klein-Nishina Formula
In Compton scattering, a photon collides with an electron, loses some of its energy and is
de ected from its original direction of travel (Fig. 4). The basic theory of this e ect, assuming
the electron to b e initially free and at rest, is that of Klein and Nishina[16].
Y
hν θ hν0 O X φ
E
Figure 4: Compton scattering with a free electron at rest.
The relation b etween photon de ection and energy loss for Compton scattering is determined
from the conservation of momentum and energy b etween the photon and the recoiling electron.
This relation can b e expressed as:
h
0
h = ; (3)
h
0
1+ (1 cos )
2
m c
e
2 2
2(h ) cos
0
2
E = h h = m c ; (4)
0 e
2 2 2 2
(h + m c ) (h ) cos
0 e 0
1
; (5) cot tan =
h
0
2
1+
2
m c
e
where h is the energy of incident photon, h is the energy of scattered photon, E is the energy
0
of recoil electron, m is the rest mass of an electron, and c is the sp eed of light.
e
For unp olarized photons, the Klein-Nishina angular distribution function per steradian of
solid angle is
( )
2 KN 2 2
1 + cos d 1 h (1 cos )
c
2
( ) = r 1+
0
2 2
d 2 [1 + h (1 cos )] (1 + cos )[1 + h (1 cos )]
2
k k k 1
0
2 2 2 1 1
+ r sin (cm sr el ectr on ); (6) =
0
2 k k k
0 0
h h
0
k = ; k =
0
2 2
m c m c
e e 6 90
80 120 60 70 60 -1 50 30 40
electron 30 -1 20 θ deg. str. 2 10 cm 0 -27 0 10
, Ω Ω Ω Ω /d C σ σ σ σ d
210 330
0.01MeV 0.016MeV 240 0.2MeV 300 0.5MeV 1.2MeV270
3.0MeV
Figure 5: Di erential cross section of Compton scattering (original gure from [17 ]).
where r is the classical electron radius. This angular distribution is shown in Fig. 5.
0
Integration of Eq. 6over all angles gives the total Klein-Nishina cross section,
1+k 2(1 + k ) ln(1 + 2k ) ln(1 + 2k ) 1+3k
2 1 KN 2
+ (cm el ectr on ): =2r
C 0
2 2
k 1+2k k 2k (1 + 2k )
(7)
4.1.1 Electron binding correction in Compton scattering
The Klein-Nishina formula for Compton scattering from atomic electrons assume that the elec-
trons are free and at rest, which is a go o d approximation for photons of the order of 1 MeV or
higher, particularly for low-Z target materials.
Binding corrections have usually b een treated by the impulse approximation or by Waller-
Hatree theory[18 ], taking account not only of the K-shell, but all of the atomic electrons. This
involves applying a multiplicative correction, the so-called \incoherent scattering function",
S (q; X ), to the di erential Klein-Nishina formula,
BD KN
d d ( )
C C
( )=S (q; Z ) : (8)
d d
The momentum transfer, q , is related to the photon energies and de ection angle, , according
to
q
2
2
q = k + k 2k k cos ; (9)
0
0
2 2
where q is in m c units and k and k are the photon energies (in m c units) b efore and after
e 0 e
de ection.
Incoherent scattering functions have b een surveyed and studied by Hubb elle et al.[3,4]and
are presented for all elements Z=1 to 100, for photon energies 100 eV to 100MeV.
KN BD
for photon energies of 1, 10, and 100keV to Fig. 6 shows the ratios of
C C 7 1 100 keV
10 keV KN C σ σ σ σ / 0.1 BD C σ σ σ σ 1 keV
0.01 0 20406080100
Atomic Number, Z
BD
Figure 6: Ratio of the b ound-electron Compton scattering cross section, [9], to that for free
C
KN
electrons, ,evaluated from the Klein-Nishina formula.
C
4.1.2 Doppler broadening
In the energy region where the binding correction is necessary, the motion of the atomic electrons
around the atomic nucleus gives rise to a Doppler broadening of the apparent energy of the
incident photon, resulting in a corresp onding broadening of the Compton \mo di ed line" for
a given de ection angle of the outgoing scattered photon. The shap e of this broadened line
is called the \Compton pro le." Available data on Compton pro les have been surveyed by
Hubb ell, and are presented in Ref. [4 ].
If one is interested in the sp ectral distribution of Compton-scattered photons, and not just
the integrated cross section, Namito et al.[19 ] have made it dramatically clear by computing
the scattered sp ectrum with and without including Doppler broadening in the EGS4 co de[20 ],
b ecause such broadening must b e included in order to replicate the corresp onding measurements
using 40 keV photons from a 2.5 GeV synchrotron light source (KEK-PF) (see Fig. 7). 8 -2 10-2 10 C (a) Cu (b)
-3 -3 source Measurement
10 source 10 Measurement EGS4(CP)
per EGS4(CP)
per
EGS4(S) -1
-1 EGS4(S) 10-4 10-4 keV
keV
-1 -1 sr.
sr. -5 10 10-5 Photons
-6 Photons 10 10-6 30 32 34 36 38 40 30 32 34 36 38 40
Photon Energy, keV Photon Energy, keV
Figure 7: Comparison of the photon sp ectra scattered by C and Cu samples. The measured and
calculated values are shown in symb ols and histograms, resp ectively (k =40 keV, =90 ). The
0
EGS4 calculation including Doppler broadening is shown as EGS4(CP). The EGS4 calculation
without Doppler broadening is shown as EGS4(S). The p eak at 30 keV is the Ge K X-ray escap e
p eak of the Rayleigh-scattering p eak at 40 keV.
4.2 Rayleigh scattering
Rayleigh scattering is a pro cess by which photons are scattered by b ound atomic electrons and
in which the atom is neither ionized nor excited. The scattering from di erent parts of the
atomic charge distribution is thus \coherent," i.e., there are interference e ects.
This pro cess occurs mostly at low energies and for high-Z material, in the same region
where electron binding e ects in uence the Compton-scattering cross section. In practice, it
is necessary to consider the charge distribution of all electrons at once. This can be done
approximately through the use of an \atomic form factor", F (q; Z ), based on the Thomas-
2
Fermi, Hartree, or other mo del of the atom. The square of this form factor, [F (q; Z )] , is the
probability that the Z electrons of an atom take up the recoil momentum, q , without absorbing
any energy. In this pro cess q is well-represented by the following equation, since it is assumed
that k k =0:
0
q =2k sin : (10)
2
The di erential Rayleigh scattering cross section for unp olarized photons is
2
r
R
0
2 2 2 1 1
)[F (q; Z )] (cm sr atom ): (11) ( )= (1 + cos
d 2
The cumulative angular distributions of , based on the values of F (q; Z ) from Nelms and
R
Opp enheim[21], are displayed in Fig. 8 for C, Fe, and Hg. At high energies it is seen that Rayleigh
scattering is conned to a small angle (at 1 MeV more than half the photons are scattered bylee
than 5 ). At low energies, particularly for high-Z materials, the angular distribution is much
broader. 9 100 95% 75% C 6 50%
25% 10 (degree)
θ
1 0.01 0.1 1 Photon Energy, MeV
95% 100 Fe 75% 26 50%
25% 10 (degree)
θ
1 0.01 0.1 1 Photon Energy, MeV
95% 100 Hg 75% 80 50%
25% 10 (degree)
θ
1 0.01 0.1 1
Photon Energy, MeV
Figure 8: Op ening half-angle, , of cone containing 25%, 50%, and 95%, resp ectively, of photons
Rayleigh scattered from C, Fe, and Hg. 10
4.3 Polarization e ects
Linear p olarization is imp ortant in low-energy photon transp ort, b ecause photons are linearly
p olarized in scattering; further, the Compton and Rayleigh scattering of linearly p olarized pho-
tons are anisotropic in their azimuthal distributions in the low-energy region.
For linearly p olarized photons with a free electron, the Klein-Nishina angular distribution
function p er steradian of solid angle is[16]
2
KN
d 1 k k k
0
2 2
C
= r + 2 + 4 cos : (12)
0
d 4 k k k
0 0
Here, is the angle b etween the incident p olarization vector (~e ) and the scattered p olarization
0
vector (~e). According to Heitler[22 ], two directions are considered for ~e, i.e. ~e either in the same
plane as ~e (~e ) or p erp endicular to ~e (~e ).
0 0
?
k
The elastic scattering of a photon by one electron is called Thomson scattering. The
Thomson-scattering cross section p er electron for a linearly p olarized photon is[23]
d
T
2 2
= r cos : (13)
0
d
Linearly p olarized photon scattering for b oth Compton- and Rayleigh-scattering have been
implemented to the EGS4 co de by Namito et al[24 ]. The e ect of linear p olarization is clearly
presented in Fig. 9[25 ], whichgives a comparison b etween EGS4 calculations with and without
including linear p olarization and Doppler broadening and corresp onding measurements using 40
keV linearly p olarized photons from a KEK-PF.
10-2 10-2 Pb 40 keV Rayleigh Pb 40 keV -3 -3 Exp H 10 10 Exp H Exp V Comton Exp V EGS4 H Default EGS4 H LP+DB
EGS4 V Default -4 EGS4 V LP+DB 10-4 10 (/keV/sr./source)
(/keV/sr./source)
-5 10-5 10 Counts Counts
-6 10-6 10 32 34 36 38 40 42 32 34 36 38 40 42 Energy Deposition (keV) Energy Deposition (keV)
(a) (b)
Figure 9: Comparison of the photon sp ectra scattered by a Pb sample. Source photons are
linearly p olarized in the vertical direction. The measured and calculated values are indicated
bysymb ols and histograms, resp ectively (k =40 keV, =90 ). (a) is the comparison with the
0
EGS4 calculation without including linear p olarization and Doppler broadening and (b) is the
one with linear p olarization and Doppler broadening. 11
5 Electron-Positron Pair Pro duction
In this e ect, which is the most likely photon interaction at high energies, a photon disapp ears
in the eld of a charged particle, and an electron-p ositron pair app ears.
For the reaction M + M ! M + M + M + Q, it can b e shown from the conservation
1 2 3 4 5
of energy and momentum that the threshold energy for the reaction in lab oratory system is
Q
2 2 lab
[Q 2(M c + M c )] (14) T =
1 1
th
2
2M c
2
when M is at rest. In a pair-pro duction reaction ( + M ! M + m + m + Q),
2 e e
M = 0;
1
M = M = M;
2 3
M = M = m ;
4 5 e
so that
2 CM
Q = 2m c ( T ) (15)
e
th
and
2 2 2
2m c (m c + Mc )
e 2
lab
T = : (16)
thh
2
Mc
Thus;
1. Pair pro duction in the eld of a nucleus of mass M (M m ):
e
2
2m c
e
2 2 lab
(Mc )=2m c =1:022 )MeV ) (17) T '
e
th
2
Mc
2. Pair pro duction in the eld of an electron (M = m ):
e
2
2m c
e
lab 2 2 2
= T (m c + m c )=4m c =2:044 )MeV ) (18)
e e e
th
2
m c
e
The cross section for pair pro duction in the eld of a nucleus varies as
n
2
Z : (19)
n
For low photon energies,
N
A
) ln(h ): (20) (=
n n
uA
For high energies[26 ],
7
; (21)
n
9X
0
where X is the radiation length of the material. The quality of the approximation is shown in
0
Table 2 ( for 1000 MeV).
n 12
Table 2. Pair-pro duction cross section in the eld of a nucleus[27].
7