PHYS 5012 Radiation Physics and Dosimetry

Transcription

1 PHYS 5012 Radiation Physics and Dosimetry Tuesday 17 March 2009

2 What are the dominant photon interactions?

3 (cont.) Compton scattering, the photoelectric effect and pair production are the three main energy transfer mechanisms in photon interactions with matter. Each of these processes can dominate under specific conditions determined chiefly by the incident photon energy hν and atomic number Z of the absorber.

4 (cont.) Klein-Nishina differential cross section for Compton scattering measures probabilty of photon re-emission into solid angle dω = dφ d cos θ as a result of a Compton interaction between an incident photon and a free electron: d e σ KN c dω = 1 ( ) ν 2 ( ν 2 r2 e ν ν + ν ) ν sin2 θ (1) can be written in terms of Thomson differential cross section, d e σ T /dω = 1 2 r2 e(1 + cos 2 θ), and a form factor: d e σ KN c dω = d eσ T dω F KN (2)

5 Klein-Nishina form factor: F KN = [1 + ε (1 cos θ)] 2 (3) { ε 2 (1 cos θ) 2 } 1 + [1 + ε(1 cos θ)](1 + cos 2 θ) The Klein-Nishina form factor plotted against scattering angle for different incident photon energies ε = hν/m ec 2. (Fig. 7.9 in Podgoršak).

6 Klein-Nishina differential cross section forward scattering (θ 0) and backward scattering (θ π) have equal probability in Thomson limit (ε 0) probability of backscattering decreases with increasing ε forward beaming of photon re-emission

7 Energy distribution of recoil electrons The differential electronic Klein-Nishina cross section can also be expressed as a function of the recoil electron kinetic energy, E K, rather than scattering angle θ since E K = E K (θ) (c.f. eqn. 22 in last lecture): d e σc KN = d eσc KN de K dω = πr2 e εhν dω dθ [ 2 where ξ K = E K /hν. dθ (4) de K 2ξ K ε(1 ξ K ) + ξk 2 ε 2 (1 ξ K ) 2 + ξk 2 ] (1 ξ K )

8 The Klein-Nishina differential electronic cross section for Compton scattering plotted as a function of recoil electron kinetic energy E K. For a given photon energy, the maximum recoil energy is indicated. Note that these photon energies are all in the limit ε 1 (Fig in Podgoršak). energy distribution peaks sharply near (E K ) max = 2ε 1 + 2ε hν c.f. eqn. (25) in last lecture

9 Total electronic Klein-Nishina cross section eσ KN c = de σc KN dω { 1 + ε = 2πr 2 e + dω = 2π ε 2 ln(1 + 2ε) 2ε +1 1 [ 2(1 + ε) 1 + 2ε 1 + 3ε (1 + 2ε) 2 d e σc KN dω d cos θ ] ln(1 + 2ε) ε } (5)

10 Limiting solutions: ε 1: eσ KN c 8π 3 r2 e ( 1 2ε ε ε ) 35 ε4... ε 0 : e σ KN c 8π 3 r2 e = e σ T = m 2 = b =Thomson limit ε 1 : e σ KN c πr 2 e(1 + 2 ln ε)/(2ε) (hν) 1 at high photon energies, Klein-Nishina cross section declines rapidly with respect to the Thomson cross section

11 Energy Transfer Mean fraction of incident photon energy hν transferred to kinetic energy E K of the recoil electron is an average of the fractional kinetic energy E K /hν weighted over the probability distribution P(θ) for Compton scattering in direction θ, integrated over all scattering angles: E K σ hν = EK hν P(θ) d cos θ (6) P(θ) d cos θ where P(θ) = 1 eσ KN c de σ KN c dω 2π dφ = eσc KN d e σ KN c dω and E K = E K (θ) is given by eqn. (22) in the last lecture. Note that P(θ)d cos θ = 1. (7)

12 We can now write where E K σ hν = 2π eσc KN = 2π eσc KN = ( eσc KN eσ KN c ) tr E K hν d e σ KN c dω d cos θ (d e σc KN ) tr dω d cos θ ( e σc KN ) tr = electronic energy transfer cross section (d e σc KN ) tr /dω = differential energy transfer cross section (8)

13 Differential energy transfer cross section (d e σc KN ) tr dω E K = d eσc KN (9) dω hν = 1 ( ) ν 2 ( ν 2 r2 e ν ν + ν ) ( ) ν ν ν sin2 θ ν where E K = hν hν is the kinetic energy imparted to the recoil electron. From eqn. (22) in the last lecture, the θ dependence is: E K /(hν) = ε(1 cos θ)[1 + ε(1 cos θ)] 1 Note: Podgoršak has some incorrect factors of Ē K and Ē σ K in the expression for (d e σc KN ) tr /dω given by eqn. (7.42). These should be E K, except on the very last line, where Ē σ K /hν should be deleted.

14 Total energy transfer cross section: d(e σ ( e σc KN c KN ) tr ) tr = dω = 2π dω +1 1 d( e σc KN ) tr d cos θ dω see eqn. (7.51) in Podgoršak for full solution. Recall from eqn. (8) that ( e σc KN ) tr eσc KN = E K σ hν

15 Atomic cross section e σ KN c is for free electrons independent of Z at high photon energies (hν E B, where E B = electron binding energy), total Compton cross section for entire atom is aσ KN c = Z ( e σ KN c ) (10)

16 Binding energy effects assumption of free electrons breaks down for photon energies hν E B = electron binding energy overestimates effective Compton atomic cross section a σ C at low hν, especially for high Z material a σ KN c The Compton atomic cross section aσ C (solid curves) compared to the Klein- Nishina atomic cross section aσc KN = Z eσc KN (dashed curves), demonstrating electron binding effects at low hν. (Fig in Podgoršak).

17 binding energy correction to Compton atomic cross section usually has negligible effect on overall attenuation because other processes (e.g. Rayleigh scattering, photoelectric effect) are usually more important than Compton scattering at low hν and high Z The atomic cross section for Rayleigh scattering (solid curves) compared to that for Compton scattering (dotted curves) plotted against incident photon energy for varying Z atoms. (Fig in Podgoršak)

18 Compton mass attenuation coefficient σ c ρ = N A A a σ C (11) Compton mass energy transfer coefficient ( ) σc = σ σ c E K ρ ρ hν tr (12) Atomic attenuation coefficient and atomic energy transfer coefficient for lead. (Fig in Podgoršak)

19 Example: For hν = 1 MeV photons incident on lead (Z = 82, A = kg), the Compton atomic cross section is aσ C = m 2. This can be calculated directly from the expression for e σc KN given by (5) and then using a σ C = Z e σc KN (since binding energy corrections are negligible at this hν). The Compton mass attenuation coefficient is σ c ρ = N A A a σ C = kg m 2 = m 2 kg 1 The values can be checked by going to the NIST/Xcom database: physics.nist.gov/physrefdata/xcom/text/xcom.html The average fractional recoil energy is obtained by inserting ε = 1.96 into the full solution given by eqn. (7.54) in Podgoršak: E K σ hν 0.440

20 interaction between photon and tightly bound orbital electron photon completely absorbed, electron ejected momentum transfer to atom, but recoil negligible due to relatively large nuclear mass, so energy conservation is: E K = hν E B photoelectron kinetic energy (13) E B = binding energy of electron orbital approx. 80% occur with K-shell electrons resulting shell vacancy quickly filled by a higher shell electron; resulting transition energy released either as: characteristic (fluorescent) photon Auger electron probability determined by fluorescent yield

21 Atomic cross section a τ = atomic cross section for photoelectric effect function of hν, exhibits characteristic sawtooth behaviour: sharp discontinuities coinciding with E B of a particular shell absorption edges Atomic cross section for photoabsorption. indicated. (Fig in Podgoršak.) Energies of K-shell ionisation are

22 e.g. Lead has prominent edges at the following ionisation energies for respective shells: K edge 88.0 kev L 1 edge 15.9 kev L 2 edge 15.2 kev L 3 edge 13.0 kev M edge 3.9 kev

23 3 distinct regions characterise a τ: 1. in immediate vicinty of absorption edge: a τ poorly known; for K-shell electrons, a τ K ε 3 assumed 2. away from absorption edge: cross section for K-shell electrons is aτ K 32 α 4 eσ T Z n ε 7/2 (14) where α = 1/137 = fine structure constant, n is a power index for Z dependence 3. ε 1: aτ K 1.5 α 4 eσ T Z 5 ε 1 (15)

24 Energy transfer Auger electrons sometimes produced mean energy transfer to electrons as a result of photoelectric effect can be in range (no Auger electron) hν E B < E τ tr < hν (Auger electron) fluorescent yield for K-shell, ω K : ω K = 1 characteristic emission only, no Auger electrons ω K = 0 no characteristic emission, Auger electrons only in general, E τ tr = hν P K ω K hν K (16) where P K = fraction of all photoelectric interactions that occur in K-shell for photons with hν > E B (K) (see Fig. 3.3 in Podgoršak).

25 hν K = K-shell weighted mean of all possible fluorescent transitions (L K, M K); K α usually most probable, giving hν K 0.86E B (K) K-shell binding energy, E B(K), weighted mean characteristic X-ray energy for all transitions to K-shell, hν K, and mean energy of K-shell characteristic emission, P K ω K hν K. (Fig in Podgoršak.)

26 Example: Consider hν = 0.5 MeV photons incident on lead, which has E B (K) = 88 kev. Suppose the photoelectric effect occurs and is immediately followed by a forbidden K α3 (L 1 K) transition, with ejection of an Auger electron from the L 2 shell. The total energy transferred to electrons = photoelectron energy + Auger electron energy: E tr = [ hν E B (K) ] + [ hν Kα3 E B (L 2 ) ] But hν Kα3 = E B (K) E B (L 1 ), so this simplifies to E tr = hν E B (L 1 ) E B (L 2 ) MeV 0.94hν For all K-shell transitions, the average energy transfer is E τ tr = hν P K ω K hν K 500 kev 65 kev = MeV 0.87hν

27 Mass coefficients for the photoelectric effect mass attenuation coefficient: τ ρ = N A A a τ (17) mass energy transfer coefficient: ( τk ρ ) tr = a τ K ρ E τ ( tr hν = a τ K 1 P ) Kω K hν K = a τ K ρ hν ρ f τ (18) where f τ = mean fraction of energy hν transferred to electrons

28 f τ 1 for low-z absorbers because Auger effect is more prevalent (i.e. fluorescent yield ω K 0) Mean fraction of photon energy hν transferred to electrons in a K-shell photoelectric interaction. (Fig in Podgoršak.)

29 production of electron-positron (e e + ) pair resulting from photon interaction with atomic nucleus incident photon energy must exceed threshold 2m e c 2 = 1.02 MeV triplet production (e e e + ) results when incident photon interacts with orbital electron; higher threshold energy required: 4m e c 2 = MeV

30 Relativistic kinematics Particle 4-momenta: photon: p ν = hν c (1, ˆk) where ˆk = unit vector in direction of photon 3-momentum (i.e. propagation direction). ( ) ( ) E E electron and positron: p e = c, p e, p e + = c, p e + where p e = γβm e c = p e + are the electron and positron 3-momenta (must have same magnitude, but can have different direction). Must also consider 4-momentum of atom, p a = (E a /c, p a ), which can gain recoil energy.

31 Conservation of 4-momentum: before: p ν = p e + p e + + p a after Note that the modulus of a 4-vector A = (A 0, A) is: A 2 = A µ A µ = (A 0 ) 2 + A A which implies that (p ν ) 2 = 0 always and (p e ) 2 = m 2 ec 2. If we square the conservation equation above, then (p ν ) 2 = (p e + p e + + p a ) 2 = 0 Now consider the case where p a = 0. We have (p e ) 2 + 2p e p e + + (p e +) 2 = 0 ( ) E = 2m 2 ec c 2 + p e p e + = 0

32 = 2(γ 2 1)m 2 ec 2 (1 + cos θ e ) = 0 which can only be satisfied if the electron and positron are emitted in exactly opposite directions, with separation angle θ e = π. In general, therefore, the atom must gain some recoil energy from the collision with its nucleus. Because of its relatively large mass, however, the recoil gained by the atom will be small.

33 Energy transfer to electrons and positrons in pair production interactions: (E K κ ) tr = hν 2m e c 2 (19) is the total kinetic energy gained by the particles (ignoring atom recoil energy). Generally, the electron and positron can be emitted with different kinetic energies, but their average energy still satisfies or, for triplet production E K pp = 1 2 (hν 2m ec 2 ) (20) E K tp = 1 2 (hν 4m ec 2 ) (21) Nuclear screening occurs for hν > 20 MeV photons that interact with the nuclear Coloumb field outside the K-shell; effective nuclear field is screened by two K-shell electrons and the interaction cross section is reduced.

34 Example: For a 2 MeV photon, the average energy of charged particles resulting from pair production in the nuclear field is pp 1 E K = 2 (hν 2m ec 2 ) = 1 (2 MeV MeV) 2 = MeV = hν and in the electron field, the average energy is E K tp = 0 because 2 MeV is less than the threshold energy 4m e c 2 = 2.04 MeV needed for triplet production.

35 Atomic cross section General form for pair production atomic cross section in field of nucleus or orbital electron is aκ = αr 2 ez 2 P(ε, Z) (22) P(ε, Z) = complicated function field energy P(ε, Z) 1. nucleus 1 ε (αz 1/3 ) ln(2ε) nucleus ε (αz 1/3 ) ln(183z 1/3 ) nucleus ε > 4 9 ln(2ε) electron ε > 4 Z 1 [ 28 9 ln(2ε) 11.3] nuclear screening only important in case 2 ε > 4 in case 3 must lie outside the limits of cases 1 and 2

36 P(ε, Z) Z 1 for electron field triplet production usually makes negligible contribution to total cross section: aκ = a κ pp + a κ tp = a κ pp [ 1 + (ηz) 1 ) (23) where η 1 as hν Atomic cross sections for pair production (solid curves) and triplet production (dashed curves) for a high-z and low-z absorber. (Fig in Podgoršak.)

37 Total atomic cross sections for pair production (including triplet production) for different Z. (Fig in Podgoršak.)

38 mass attenuation coefficient for pair production: κ ρ = N A A a κ (24) mass energy transfer coefficient: ( ) κ = κ (EK) κ tr ρ tr ρ hν = κ ρ (1 2m ec 2 ) hν (25) Mass energy transfer coefficient (solid curves) and mass attenuation coefficient (dashed curves) for pair production. (Fig in Podgoršak.)

39 of Photon Interactions Tabulation of interactions and symbols used: electronic atomic linear cross cross attenuation section section coefficient [m 2 ] [m 2 ] [m 1 ] Thomson scattering eσ T a σ T σ T Rayleigh scattering aσ R σ R Compton scattering eσ C a σ C σ C Photoelectric effect aτ τ Pair production aκ pp κ pp Triplet production eκ tp a κ tp κ tp

40 Tabulation of attenuation coefficients: Total linear and mass attenuation coefficients are a sum of the partial linear and mass attenuation coefficients for individual photon interactions: µ = τ + σ R + σ C + κ (26) µ m = τ ρ + σ R ρ + σ C ρ + κ ρ Similarly, the total atomic cross section is (27) aµ = µ m A N A = a τ + a σ R + a σ C + a κ = Z e µ (28) where e µ = total electronic cross section.

41

42 For compounds or mixtures, the mass attenuation coefficient is a summation of a weighted average of its constituents: µ i µ m = Σ i w i (29) ρ where w i = proportion by weight of i-th constituent. Example: Water, H 2 O, has w H = w O = = = For 1 MeV photons, µ H /ρ = m 2 kg 1 and µ O /ρ = m 2 kg 1 (data from the NIST/Xcom database), so µ ρ m 2 kg 1

43 Comparison of the three main photon energy transfer processes: photoelectric effect dominates at low hν and high Z pair production dominates at high hν and high Z Compton scattering dominates over a broad range in hν for low/moderate Z (including water and tissue)

44 Energy transfer coefficient is the sum of the partial energy transfer coefficients for the photoelectric effect, Compton scattering and pair production: µ tr = µ E tr hν E σ tr = τ Eτ tr hν + σ C hν + κeκ tr hν = τ f τ + σ C f σ + κ f κ (30) where each f is the average fraction of incident photon energy hν transferred to electrons by the corresponding physical process, with f τ = 1 P Kω K hν K hν f σ = Eσ tr hν f κ = 1 2m ec 2 hν (31) (32) (33)

45 What is the dominant photon interaction?