Calculation of Mass Stopping Power and Range of Protons as Well as Important Radiation Quantities in Some Biological Human Bodyparts (Water, Muscle, Skeletal and Bone, Cortical) ()
1. Introduction
In many fields, such as radiation dosimetry, radiation biology, and many others, radiation chemistry, radiotherapy, and nuclear physics, the stopping power, energy loss, rangestraggling, and equivalent dose rate of ions in the air, tissue, and polymers are very important. The use of protons or heavier ions as an alternative to external photon beams in radiotherapy is increasing, with the reason being that they preserve the target dose, ensure a higher dose delivered to the tumor, and can transfer energy in the form of a point shot through diseased tissue due to the Bragg curve [1]. The stopping power of charged particles has been measured using a variety of ways, including direct energyloss measurements through films, backscattering from thick substrates with deposited absorbing layers, gamma resonance shift measurements, self-supporting methods, and indirect verification of the stopping power based on alpha energy losses in the air have all been reported as methods for measuring the stopping power of charged particles. [2] [3] [4].
In the present work, the electronic mass stopping power and range of proton in some biological human body parts (Water, Muscle, Skeletaland Bone, Cortical) are calculated using the Bethe-Bloch formula as reported in the references.
As it is known in any therapeutic unit with protons, it needs to calculate the absorbed dose, the equivalent dose to the tissue, and the effective dose according to the energy of the protons. Therefore, in this work, a variety of radiation quantities such as thickness, absorbed dose, equivalent dose, and effective dose of the protons in Water, Muscle, Skeletal and Bone, Cortical were also computed in proton energy range 0.04 - 200 MeV.
2. Methods
2.1. Calculations of Electronic Mass Stopping Power
Bethe was the first person to use quantum mechanical studies on stopping power. The Bethe theory of stopping power is valid when the projectile’s velocity surpasses the Bohr velocity. In Bethe’s theory, the goal is assumed to be charged particle. In Bethe’s approach to energy loss, the Born approximation is employed to represent inelastic collisions between heavy particles and atomic electrons. In this theory, the projectile heavy particle is as assumed to be structureless, whereas the target nucleus is assumed to be infinitely massive [4]. For the energy range 0.04 - 200 MeV, the Bethe mass stopping power equation [4] [5] [6] [7] was used:
(1)
where β is v/c wherev is the proton velocity and c is light velocity,I is the mean excitation energy and F(β) is given by
(2)
n is calculated using the following relation:
(3)
where Na Avogadro number, ρ is the density of substances and Ζ/A is the ratio of atomic number to the mass number of substances. The basic data for human body tissues are given in Table 1. The calculated mass stopping power of protons for Water, Muscle, Skeletal and Bone, Cortical are based on Bethe equation after substituting the constants from Table 1. In Table 2 the compositions of the human tissues are given [8].
2.2. Calculations of Range
The range of a heavy particle is the straight distance it travels within the target. Light particles like electrons and positrons scatter widely throughout the path of targets due to their low mass, making it difficult to determine their journey duration. The path length of light particles has been calculated with remarkable success using Monte Carlo methods, which are based on a broad class of computational algorithms. On the other hand, heavy particles like protons have a practically straight line path length. The range of protons can be calculated using numerical integration methods. The Continuous Slowing Down Approximation (CSDA), on the other hand, is a straightforward and extensively used method for finding a variable’s range. This study used a simple and standard method for calculating the range of heavy particles such as protons in the targets. The CSDA approach uses incident particles to constantly lose energy in the route of the targets. neglects energy loss fluctuations. The range, R for an incident proton in the CSDA method is given as [1] [4] [7]:
(4)
Table 1. Basic data for calculating mass stopping powers.
Table 2. Elemental composition of Water, Muscle, Skeletal and Bone, Cortical tissues.
where,
is the initial energy of incident charged particle in material,
is the final energy of incident charged particle in material and
is the mass stopping power.
2.3. Calculations of Thickness, Absorbed Dose, Equivalent Dose and Effective Dose
Thickness of proton in the tissue or substance is given by
(5)
where R (in g/cm2) is the range and is ρ the density of the tissue or substance [9]
Absorbed dose in (in rad): It is the transfer of an amount of energy of 100 erg per gram of the absorbent material when protons pass over it and is given by the following relation: [9]
(6)
where E is the proton energy.
Equivalent dose: is given by
(7)
where D is the absorbed dose and
is the weighting factor of radiation (proton) [9]
Effective dose: is given by
(8)
where
is the weighting factor of tissue or substance [9];
where
for Water, Muscle, Skeletal and Bone, Cortical and
for protons [10].
2.4. The Percentage Error of Difference
The percentage error of difference for the mass stopping power and range of protons in given biological human body parts is calculated by
(9)
3. Results and Discussion
The results of mass stopping power and range of protons in some biological human body parts (Water, Muscle, Skeletal and Bone, Cortical) are given in Table 3. In Figures 1-3 the mass stopping power of biological human body parts and their compositions are plotted using MathLab program and it is noted that the mass stopping power of all substances and tissues is approximately equal to t the average values of its compositions. The comparison between the present calculated electronic mass stopping powers and that of PSTAR program [8] are
Table 3. Values of electronicmass stopping power (in MeV cm2/g) of Water, Muscle, Skeletal and Bone, Cortical tissues.
Figure 1. Mass stopping power of Water and its composition.
Figure 2. Mass stopping power of Muscle, Skeletal and its composition.
Figure 3. Mass stopping power of Bone, Cortical and its composition.
shown in Figures 4-6 angood agreements between two results are observed in energy range 1 - 200 MeV.
In Table 4 the ranges of protons in some Biological human body parts (Water, Muscle, Skeletal and Bone, Cortical) are given. In Figures 7-9 a comparison between present results of ranges and that of PSTAR data are shown in energy
Figure 4. Mass stopping power of Water versus
.
Figure 5. Mass stopping power of Muscle, Skeletal versus
.
range 0.1 - 200 MeV and good agreements are observed. In Table 5 and Table 6 the thickness, LET, absorbed dose, effective and equivalent dose are given. In Table 7 and Table 8, the empirical formulae for calculating mass stopping powers and ranges for protons in some biological human body parts (Water, Muscle Skeletal and Bone, Cortical) are given with the percentage difference error.
4. Conclusion
In this work, calculations of mass stopping power and range of protons incident
Figure 6. Mass stopping power of Bone, Cortical versus
.
Figure 7. Range of Water versus
.
on the three different biological human parts (Water, Muscle, Skeletal and Bone, Cortical) have been done and the following conclusions are drawn:
1) The mass stopping power of the Water, Muscle, Skeletal and Bone, and Cortical is equal to the average value of mass stopping power of their compositions in energy range 0.04 - 200 MeV.
2) Values for mass stopping power and ranges of protons in Water, Muscle, Skeletal and Bone, Cortical are in good agreement with the data of PSTAR program. The percentage value of error difference was between 0.09% - 28.88% for
Figure 8. Range of Muscle, Skeletal versus
.
Figure 9. Range of Bone, Cortical versus
mass stopping power and was between 3.99% - 79.47% for range at proton energy ranging 0.04 - 200 MeV.
3) It was also observed that the maximum value of mass stopping power was at 0.1 MeV and after that, the mass stopping power start to decrease with increasing proton energy as expected.
Table 4. Values of range (in g/cm2) of Water, Muscle, Skeletal and Bone, Cortical tissues.
Table 5. Thickness, LET, absorbed dose, equivalent dose and effective dose of proton in Water and Muscle, Skeletal.
Table 6. Thickness, LET, absorbed dose, equivalent dose and effective dose of proton in Bone, Cortical.
Table 7. The empirical formulae for calculating mass stopping powers.
Table 8. The empirical formulae for calculating range.
4) The empirical formulae suggested for mass stopping power are simple and accurate at proton energy (1 - 200 MeV) while the empirical formulae suggested for proton range give a good result compared to calculated values at proton energy (4 - 200 MeV).
5) Also in this study, some radiation quantities were calculated, such as linear energy transfer, adsorbed dose, and effective and equivalent dose that give good information to those interested in proton therapy.
Acknowledgements
The authors thank the Department of Physics, Qassim University, for their support and encouragement of this research.