Charged particles can be divided into two families: 1) electrons, including positrons and beta particles from radioactive decay, Auger electrons, and internal conversion electrons, and; 2) heavy charged particles such as the alpha particle, fission fragments, protons, deuterons, tritons, and mu and pi mesons.

The charged particles interact with matter primarily through the Coulomb force. Differences in the energy deposition trails between the two types of charged particles are due to the mass differences, i.e., the more easily altered trajectory of the electron.

The maximum energy that can be transferred from a charged particle to an electron in a single collision is given by 4 T m(e)/m, for about 1/500 of the particle energy per nucleon for heavy particles. Due to this small fraction and the fact that at any time the particle is interacting with more than one electron, the heavy particles appear to continuously slow down, gradually losing their kinetic energy along an unaltered linear path.

We will concern ourselves primarily with two types of interactions: collisional and radiative.

- TYPE 1, COLLISIONAL. Inelastic collision with atomic electrons. This results in excitation or ionization. These processes ultimately end with the heating of the absorber (through atomic and molecular vibrations) unless the ions and electrons can be separated using an electric field as is done in radiation detectors.
- TYPE 2, RADIATIVE. Inelastic collision with nucleus. A quantum of electromagnetic radiation is emitted (a photon). Energy loss is experienced by the particle. Important for electrons. Probability of nuclear excitation is negligible. This process is also known as
**RADIATIVE**energy loss. The acceleration of electrons near a nucleus is known as beam braking or bremsstrahlung. Bremsstrahlung will be discussed further in the next section.

Eqn 1: S= - dT/dx,

where T is the charged particle kinetic energy, -dT is the energy increment lost in infinitesimal material thickness of dx. The units of stopping power are keV/micron. The higher the stopping power, the shorter the range into the material the particle can penetrate. The quantity S is also referred to as *specific energy loss *.

The stopping power S increases as the particle velocity is decreased. The classical depiction of the charged particle interaction with an electron is depicted in the following drawing:

Classically, a passing charged particle will impart kinetic energy to the electron given by the following formula:

The faster the particle is, the less energy is given to the electron (the less time it has to spend imparting energy to the electron).

The classical (i.e, non-Quantum Mechanical) expression that describes the specific energy loss is known as the Bethe-Bloch formula and is written (Eqn 2):

The first term in the stopping number B is sufficient if v<

The Bethe-Bloch formula for electrons contains must take into account two facts: 1) that the electron has a much smaller mass than the heavy particles; and 2) that the electron is identical to the particles with which it is interacting, thus giving the possibility of the exchange of identity, i.e., the incoming particle has a probability of becoming the atomic electron with the atomic electron becoming the outgoing particle. Taking this into account, as well as relativity, gives us the Bethe-Bloch formula for electrons.

The plot of specifc energy loss (which can be related to specifc ionization) along the track of a charged particle is called a Bragg Curve. A typical Bragg curve is depicted in the following graphic for an alpha particle of several MeV of initial energy:

As the energy falls, the specific energy loss increases according to the Bethe-Bloch formula. As the energy falls below a threshold, however, an electron will attach to the alpha, dramatically lowering the specific energy loss.

**ELECTRONS**

The energy deposition of the electron increases more slowly with penetration depth due to the fact that its direction is changed so much more drastically (Johns and Cunningham, p. 213). In fact, these authors state: "There is no increase in energy deposited near the end of the tract and the Bragg peak for electrons is never observed." The explanation of Sorenson and Phelps, however, helps shed light on this topic. They state (p. 170) : "A similar increase in ionization density is seen at the end of an electron track; however, the peak occurs when the electron energy has been reduced to less than about 1 keV, and it accounts for only a small fraction of its total energy." Therefore it seems that for electrons and their tortuous paths the energy deposition is spread in the transverse direction as it progresses forward in the initial direction, until the last 1 keV which is deposited pretty much along a straight line in the forward direction.

The straggling peak is approximately Gaussian shaped, with a width that increases with the ratio Z/A. In other words, at lower atomic numbers (where the Z/A ratio is higher) we have wider widths. This is due to the fact that there is less shielding of inner electrons in the lower Z elements - there is more stopping influence per electron.

(1/rho) * (dT/dx) in units of keV/(gm/cm**2)

Eqn 4: [dT/dx] total = [dT/dx]collision + [dT/dx] radiative

(for heavy particles we only consider orbital electron interactions (an atomic process) since the probability of nuclear interaction resulting in energy loss is much smaller.)

The radiative component calculated from electromagnetic and relativistic theory is given by Eqn 5:

which is the basis for the approximation for the percentage of energy lost to radiation for monoenergetic electrons, Eqn. 6:

(dT/dx)r/ (dT/dx)total = EZ/1000 ; where E is in MeV.

, where Z is the atomic number of the absorber and E is the charged particle energy in MeV.

For example, a stream of 100 keV electrons hits a tungsten target (Z=74). Only 0.74% of the energy lost goes to emitted x-rays.

For polyenergetic beta rays, Eqn 6 becomes:

Eqn 7: (dT/dx)r/(dT/dx)total = Z E(max) / 3000

since the average beta energy is 1/3 that of the maximum beta energy.

The fraction of energy lost to radiation is inversely proportional to the charged particle's mass. Bremsstrahlung will therefore be very low for heavy charged particles (compared to electrons). Even at 100 MeV, for example, heavy particles dissipate most of their energy through the collisional process.

For mixtures of elements, one determines the "effective atomic number", Z(eff). This can be found by using the following formula, Eqn 8:

Examine the thought experiment involving alpha particle transmission depicted in the following drawing.

A source of monoenergetic alpha particles is directed at a particle detector through an absorber. We know the initially monoenergetic spectrum will demonstrate energy straggling, but the direction of the particles should not be altered and they should ALL make it to the detector, provided the thickness of the absorber is less than the thickness required to stop the slowest of the distribution. This is shown as the constant horizontal line in the graph at the bottom of the drawing. Once the absorber thickness is such that some of the slower particles are stopped, the width of the energy straggling curve presents itself in the slope of the fall-off from constant particle number detected down to zero.

The mean range is the range at which the number of particles detected is one-half the original value.

The same experiment can be performed with a beam of electrons (or a beta emitter). In this case, however, the electrons are more vulnerable to wide-angle scattering even in the thinnest of foils - so wide an angle that the scattered electrons will not be detected. The result is an immediate falloff in the number of electrons detected as a function of absorber thickness as demonstrated is this drawing:

**Energy:** Range is approximately linear with energy since the Bethe-Bloch equation for stopping power is inversely proportional to E.

**Mass:** for the same kinetic energy, the electron is much faster than the alpha due to its smaller mass, and therefore the electron has less time to spend near orbital electrons. This reduces the effect of Coulomb interactions (hence stopping power) and increases range.

**Charge:** the more charge, the more stopping power and the lower range. Range is inversely proportional to the square of the charge of the charged particle.

For example, a tritium particle with z=1 will have 1/4 the stopping power of a He-3 particle with z=2.

**Density:** The stopping power increases with increasing density. The range is inversely proportional to the density of the absorbing medium.

A formula from the Bragg-Kleeman rule can tell us the range of particles in a material for which we have no energy loss or range information if we are given range information in another material. The approximation becomes less valid as the separation in atomic weights increases. This rule is as follows (Eqn 10):

This rule follows from the fact that the range in cm (linear range) is inversely proportional to the density and we are taking a ratio of ranges - one known, one unknown. The atomic number part is the approximation and the guts of the Bragg-Kleeman rule.

The estimation of a particle's range in a given material is a two-step process:

1) determine the range in air; then

2) determine the range in the material given the Bragg-Kleeman rule .

The range of alpha particles in cm depends upon the material and the energy of the alphas. The first part takes care of the energy in air; the second part takes care of the difference in materials between air and the material of choice.
**Determination of the range in air:**

In order to compute the range of alpha particles in any material, we must use an empirical equation for the range of alpha particles in air. The following equation for range in cm is valid for alphas in air in the energy range 4 < E < 15 MeV (Eqn 11):

R = (0.005E + 0.285) E**3/2 ; R in cm

(the above equation is taken from page 202 of Lapp and Andrews textbook.

Sorenson and Phelps (p. 173) use a simpler equation for the range 4 < E < 8 (Eqn. 12):

R = 0.325 E**3/2 ; R also in cm

It is helpful to note that the density of air is 0.001293 g/cm**3 and that the effective atomic number of air is about 14.

**EXAMPLE(1):**

The alpha particles of Americium-241 have energies of 5.49 MeV and 5.44 MeV. What is the range of these particles in air?

Using the above equation, we find R(cm) = 4.04 cm.

If the material is changed, use the scaling law as follows:

If medium o is air (sqrt A = 3.82), R(1) = 3.4 x 10**-4 R(air) (sqrt(A(1)/ rho(1)) ;

since 0.001293 /3.82 = 3.4 x 10**-4

. R(air) can be calculated using the above equations. Knowing more about the material in question will tell us its density and effective atomic mass number.

**EXAMPLE(2):**

For Americium-242 what is the range in tissue if the sqrt(A) = 3 and rho = 1.0 g/cm**3?

R(1) = (3.4 x 10**-4)* (3/1)*4.04 cm = 0.0039 cm

Since the mass electron range depends only on energy, we can define it in terms of semi-empirical equations. Two semi-empirical equations are:

for the range 0.01 < E < 3.0 MeV (Eqn. 13),

R = 412 E**n ; R in mg/cm**2 ; where

n = 1.265 - 0.0954 ln E

; and

for the range 1 < E < 20 MeV (Eqn. 14),

R = 530 E - 106 ; R in mg/cm**2.

In order to demonstrate that these two empirical equations overlap and represent approximately the same electron range, we show the following graph:

**EXAMPLE (1):**

What is the maximum range of the principal Beta particle of I-131?

The maximum energy of the principal Beta particle for I-131 is 0.606 MeV. The range is then:

R = 412 (0.606)**n mg/cm**2;

where n = 1.265 - 0.0954 ln (0.606)

plugging in and calculating, we find that R = 0.345 g/cm**2. If the material was water or tissue the range would be 0.345 cm.

**EXAMPLE(2):**

The following table gives an indication of the order of magnitude of the ranges of monoenergetic alphas and electrons in both soft tissue as well as air:

The LET is related to Biologiical Damage. The severity and permenance of biological chages are directly related to the local rate of energy deposition along the particle track. The higher the LET, the higher the Q- quality factor in determining dose equivalent (Severts, where 1 Sv = 100 rem).

Dr. Justine Arthur

Dr. Joseph Kannam

Dr. David Weintraub

________________________________________________________

Douglas J. Wagenaar, Ph.D., wagenaar@nucmed.bih.harvard.edu