Radiation Interactions: HCP

Radiation Interactions with Matter: Energy Deposition Biological effects are the end product of a long series of phenomena, set in motion by the passage of radiation through the medium. [Image removed due to copyright considerations] Interactions of Heavy Charged Particles Energy-Loss Mechanisms • The basic mechanism for the slowing down of a moving charged particle is Coulombic interactions between the particle and electrons in the medium. This is common to all charged particles • A heavy charged particle traversing matter loses energy primarily through the ionization and excitation of atoms. • The moving charged particle exerts electromagnetic forces on atomic electrons and imparts energy to them. The energy transferred may be sufficient to knock an electron out of an atom and thus ionize it, or it may leave the atom in an excited, nonionized state. • A heavy charged particle can transfer only a small fraction of its energy in a single electronic collision. Its deflection in the collision is negligible. • All heavy charged particles travel essentially straight paths in matter. [Image removed due to copyright considerations] [Tubiana, 1990] Maximum Energy Transfer in a Single Collision The maximum energy transfer occurs if the collision is head-on. [Image removed due to copyright considerations] Assumptions: • The particle moves rapidly compared with the electron. • For maximum energy transfer, the collision is head-on. • The energy transferred is large compared with the binding energy of the electron in the atom. • Under these conditions the electron is considered to be initially free and at rest, and the collision is elastic. Conservation of kinetic energy: 1 MV2 = 2 1 MV 12 + 2 1 mv 12 2 Conservation of momentum: MV = MV1 + mv1. Qmax = 4mME 1 1 MV2- MV 12 = , 2 2 ( M + m) 2 Where E = MV 2/2 is the initial kinetic energy of the incident particle. Qmax values for a range of proton energies. Except at extreme relativistic energies, the maximum fractional energy loss for a heavy charged particle is small. Maximum Possible Energy Transfer, Qmax, in Proton Collision with Electron Proton Kinetic Energy E (MeV) Qmax (MeV) 0.1 1 10 100 103 104 105 106 107 0.00022 0.0022 0.0219 0.229 3.33 136 1.06 x 104 5.38 x 105 9.21 x 106 Qmax = 4mME ( M + m) 2 Maximum Percentage Energy Transfer 100Qmax/E 0.22 0.22 0.22 0.23 0.33 1.4 10.6 53.8 92.1 Single Collision Energy Loss Spectra [Image removed due to copyright considerations] • The y axis represents the calculated probability that a given collision will result in an energy loss Q. • N.B., the maximum energy loss calculated above for the 1 MeV proton, of 21.8 keV is off the scale. • The most probable energy loss is on the order of 20 eV. • N.B., energy loss spectra for fast charged particles are very similar in the range of 10 – 70 eV. • Energy loss spectra for slow charged particles differ, the most probable energy loss is closer to the Qmax. Stopping Power • The average linear rate of energy loss of a heavy charged particle in a medium (MeV cm-1) is of fundamental importance in radiation physics, dosimetry and radiation biology. • This quantity, designated –dE/dx, is called the stopping power of the medium for the particle. • It is also referred to as the linear energy transfer (LET) of the particle, usually expressed as keV µm-1 in water. • Stopping power and LET are closely associated with the dose and with the biological effectiveness of different kinds of radiation. Calculations of Stopping Power In 1913, Niels Bohr derived an explicit formula for the stopping power of heavy charged particles. Bohr calculated the energy loss of a heavy charged particle in a collision with an electron, then averaged over all possible distances and energies. [Image removed due to copyright considerations] [Image removed due to copyright considerations] The Bethe Formula for Stopping Power. Using relativistic quantum mechanics, Bethe derived the following expression for the stopping power of a uniform medium for a heavy charged particle: - 4πk 02 z 2 e 4 n dE = dx mc 2 β 2 ⎡ 2mc 2 β 2 ⎤ − β 2⎥. ⎢ln 2 ⎣ I (1 − β ) ⎦ ko = 8.99 x 109 N m2 C-2 , (the Boltzman constant) z = atomic number of the heavy particle, e = magnitude of the electron charge, n = number of electrons per unit volume in the medium, m = electron rest mass, c = speed of light in vacuum, β = V/c = speed of the particle relative to c, I = mean excitation energy of the medium. • Only the charge ze and velocity V of the heavy charged particle enter the expression for stopping power. • For the medium, only the electron density n is important. Tables for Computation of Stopping Powers If the constants in the Bethe equation for stopping power, dE/dX, are combined, the equation reduces to the following form: dE 5.08 x 10 −31 z 2 n [ F ( β ) − ln I ev ] MeV cm-1 − = 2 dx β 1.02 x 106 β 2 −β2 where, F ( β ) = ln 2 1− β [Turner] Conveniently,….. For a given value of β, the kinetic energy of a particle is proportional to the rest mass, Table 5.2 can also be used for other heavy particles. Example: The ratio of kinetic energies of a deuteron and a proton traveling at the same speed is 1 M dV 2 M 2 = d =2 1 Mp M pV 2 2 Therefore the value of F(β) of 9.972 for a 10 MeV proton, also applies to a 20 MeV deuteron. Mean Excitation Energies Mean excitation energies, I, have been calculated using the quantum mechanical approach or measured in experiments. The following approximate empirical formulas can be used to estimate the I value in eV for an element with atomic number Z: I ≈ 19.0 eV; Z = 1 (hydrogen) I ≈ 11.2 eV + (11.7)(Z) eV; 2 ≤ Z ≤ 13 I ≈ 52.8 eV + (8.71)(Z) eV; Z > 13 For compounds or mixtures, the contributions from the individual components must be added. In this way a composite lnI value can be obtained that is weighted by the electron densities of the various elements. The following example is for water (and is probably sufficient for tissue). nlnI = ∑ NiZi lnIi, i Where n is the total number of electrons in the material (n = Σi NiZi) When the material is a pure compound, the electron densities can be replaced by the electron numbers in a single molecule. Example: Calculate the mean excitation energy of H2O Solution: I values are obtained from the empirical relations above. For H, IH = 19.0 eV, for O, IO = 11.2 + 11.7 x 8 = 105 eV. Only the ratios, NiZi/n are needed to calculate the composite I. Since H2O has 10 electrons, 2 from H and 8 from O, the equation becomes lnI = 2 x1 1x8 ln 19.0 + ln 105 = 4.312 10 10 giving I = 74.6 eV Stopping power versus distance: the Bragg Peak dE 5.08 x 10 −31 z 2 n [ F ( β ) − ln I ev ] MeV cm-1 − = 2 dx β • At low energies, the factor in front of the bracket increases as β → 0, causing a peak (called the Bragg peak) to occur. • The linear rate of energy loss is a maximum as the particle energy approaches 0. [Image removed due to copyright considerations] Rate of energy loss along an alpha particle track. • The peak in energy loss at low energies is exemplified in the Figure, above, which plots -dE/dx of an alpha particle as a function of distance in a material. • For most of the alpha particle track, the charge on the alpha is two electron charges, and the rate of energy loss increases roughly as 1/E as predicted by the equation for stopping power. • Near the end of the track, the charge is reduced through electron pickup and the curve falls off. Stopping Power of Water for Protons dE 5.08 x 10 −31 z 2 n − = [ F ( β ) − ln I ev ] MeV cm-1 2 dx β What is needed to calculate stopping power, - dE/dX? n z lnI the electron density the atomic number the mean excitation energy For protons, z = 1, The gram molecular weight of water is 18.0 g/mole and the number of electrons per molecule is 10. One m3 of water has a mass of 106 g. The density of electrons, n, is: n = 6.02 x 1023 molecules/mole x 10 6 g m −3 x 10 e-/molecule = 18.0 g / mole 3.34 x 1029 electrons/m3 As found above, for water, ln Iev = 4.312. Therefore, eq (1) gives − dE 0.170 -1 = [ F ( β ) − 4.31] MeV cm 2 dx β At 1 MeV, from Table 5.2, β2 = 0.00213 and F(β) = 7.69, therefore, − dE 0.170 -1 = [7.69 − 4.31] = 270 MeV cm dx 0.00213 The stopping power of water for a 1 MeV proton is 270 MeV cm-1 Mass Stopping Power • The mass stopping power of a material is obtained by dividing the stopping power by the density ρ. • Common units for mass stopping power, -dE/ρdx, are MeV cm2 g-1. • The mass stopping power is a useful quantity because it expresses the rate of energy loss of the charged particle per g cm-2 of the medium traversed. • In a gas, for example, -dE/dx depends on pressure, but –dE/ρdx does not, because dividing by the density exactly compensates for the pressure. • Mass stopping power does not differ greatly for materials with similar atomic composition. • Mass stopping powers for water can be scaled by density and used for tissue, plastics, hydrocarbons, and other materials that consist primarily of light elements. For Pb (Z=82), on the other hand, -dE/ρdx = 17.5 MeV cm2g-1 for 10-MeV protons. (water ~ 47 MeV cm2 g-1 for 10 MeV protons) **Generally, heavy atoms are less efficient on a g cm-2 basis for slowing down heavy charged particles, because many of their electrons are too tightly bound in the inner shells to participate effectively in the absorption of energy. Range The range of a charged particle is the distance it travels before coming to rest. The range is NOT equal to the energy divided by the stopping power. Table 5.3 [Turner] gives the mass stopping power and range of protons in water. The range is expressed in g cm-2; that is, the range in cm multiplied by the density of water (ρ = 1 g cm-3). Like mass stopping power, the range in g cm-2 applies to all materials of similar atomic composition. A useful relationship….. For two heavy charged particles at the same initial speed β, the ratio of their ranges is simply z 22 M 1 R1 ( β ) = 2 , R2 ( β ) z1 M 2 where: R1 and R2 are the ranges M1 and M2 are the rest masses and z1 and z2 are the charges If particle number 2 is a proton (M2 = 1 and z2 = 1), then the range R of the other particle is given by: R(β) = M Rp(β), z2 where Rp(β) is the proton range. Figure 5.7 shows the ranges in g cm-2 of protons, alpha particles, and electrons in water or muscle (virtually the same), bone, and lead. For a given proton energy, the range in g cm-2 is greater in Pb than in H2O, consistent with the smaller mass stopping power of Pb. [Image removed due to copyright considerations] References: This material was taken largely from J.E. Turner, Atoms, Radiation, and Radiation Protection, Wiley, New York, 1995, Chapter 5 “Interaction of heavy charged particles with matter”. Additional reading: E.L Alpen, Radiation Biophysics, Prentice Hall, Englewood Cliffs, New Jersey, 1990. M. Tubiana, J. Dutrice, A. Wambersie, Introduction to Radiobiology, Taylor and Francis, New York, 1990, Chapter 1.