Photon attenuation parameters for some tissues from Geant4 simulation, theoretical calculations and experimental data: A comparative study

Photon attenuation parameters for some tissues from Geant4 simulation, theoretical calculations and experimental data: A comparative study 1 corresp first-author Corresponding author, harslan@sakarya.edu.tr Corresponding author, harslan@sakarya.edu.tr harslan@sakarya.edu.tr Electrical and Electronics Engineering, Sakarya University of Applied Sciences, Sakarya, Turkey Electrical and Electronics Engineering, Sakarya University of Applied Sciences, Sakarya, Turkey Mass attenuation coefficients, effective atomic numbers, effective electron densities and Kerma relative to air for adipose, muscle and bone tissues have been investigated in the photon energy region from 20 keV up to 50 MeV with Geant4 simulation package and theoretical calculations. Based on Geant4 results of the mass attenuation coefficients, the effective atomic numbers for the tissue models have been calculated. The calculation results have been compared with the values of the Auto- Z eff program and with other studies available in the literature. Moreover, Kerma of studied tissues relative to air has been determined and found to be dependent on the absorption edges of the tissue constituent elements. Attenuation coefficient Effective atomic number Geant4 simulation Tissue model 1 Introduction 1 During more than three decades, due to the digital computing developments and wide use of radioactive sources, the attenuation measurements of X-rays and gamma rays have become an important part in numerous diverse fields of physics. Their interactions with the human tissues and organs have also occupied a special significance in medical physics by virtue of their applications in diagnostic and radiation therapy. Human organs are made of tissues that are formed by specialized cells. Each cell is composed of elements with low atomic numbers [ 1 1 - 3 3 1-3 ]. Therefore, the doses received by tissues imputed to ionization, scattering and absorption processes occured with such elements [ 4 4 ]. The mass attenuation coefficient ( μ / ρ ) and the mass energy-absorption coefficient ( μ en / ρ ) need to be known for dosimetry studies. Such basic quantities are discussed in Section 3 and detailed information about them can be found in [ 5 5 ]. The first tabulations of the mass energy-absorption coefficient were the efforts of Berger [ 6 6 ] and Allison [ 7 7 ]. Berger’s tabulation included 15 elements, air, and water at the energy region between 3 keV and 10 MeV,whereas Allison’s tabulation for seven elements, air, and water for energies from 10 keV to 100 MeV. Later on, Hubbel published data [ 8 8 ] on mass attenuation and mass energy-absorption coefficients for H, C, N, O, Ar, and seven mixtures from 0.1 keV to 20 MeV. That work has been extended two times [ 9 9 , 10 10 ] to cover all elements (for Z =1 to 92) and 48 additional substances over a wide energy range from 1 keV to 20 MeV. As a technological alternative, the XCOM program was developed to calculate the attenuation coefficients of photon for elements, compounds, and mixtures over a wide range of energies [ 11 11 ]. This program has been converted to WinXcom software [ 12 12 ] based on the mixture rule for all elements and any mixture at selected energies. In addition, Chantler [ 13 13 ] provided new tabulation that was improved as a programme FFAST [ 14 14 ]. The Monte Carlo method is commonly used to simulate radiation interaction with matter and to obtain the mass attenuation coefficient. There are some strong toolkits like Geant4 and MCNP that are currently used in the field of radiation physics. The validation of Geant4 [ 15 15 ] simulation package has been frequently examined for determining the attenuation coefficients by using different types of materials [ 16 16 - 19 19 16-19 ] and tissue-equivalent models [ 20 20 , 21 21 ]. Moreover, this work is also a validation of Geant4 as a powerful alternative tool to calculate attenuation parameters. In addition to Geant4, MCNP is also preferred to simulate the attenuation of photon in various media and gives very satisfactory results (see, for example, [ 22 22 - 26 26 22-26 ]). Effective atomic number ( Z eff ), effective electron density ( N eff ) and kinetic energy released per unit mass (Kerma) are also serviceable parameters to characterize photon interaction with human tissues. The variation of Z eff with photon energies was suggested by Hine [ 27 27 ]. Thus, accurate and quick determinations of this parameter were targeted by severals experiments, methods and programs. Such experiments, in which Z eff is determined by using the attenuation coefficients, were performed for elements [ 28 28 , 29 29 ], for mixtures [ 30 30 - 32 32 30-32 ] and for tissues or tissue-equivalent materials [ 33 33 , 34 34 ]. In addition, many theoretical studies have been conducted to obtain the attenuation coefficients and effective atomic numbers for elements and mixtures. For example, Jayachandran found that the calculated value of the effective atomic number (7.23) for lithium tetraborate Li_2B_4O_7 agrees closely with the value of 7.22 for soft tissue of composition C 5 H 40 O 18 N [ 35 35 ]. Manohara, who investigated the effective atomic number for some fatty acids and cysteine, concluded that Z eff is equal to the mean atomic number over a wide energy range around 1 MeV [ 36 36 ]. Singh, who studied dosimetric materials in terms of mass energy absorption coefficients, equivalent atomic numbers and Kerma in the energy range between 0.015 and 15 MeV, concluded that Nylon has been found to be a good tissue substitute material for making tissue-phantoms of the biological tissues [ 37 37 ]. Moreover, Taylor submitted the Auto- Z eff software [ 38 38 ] as an appropriate method to evaluate the effective atomic numbers that necessitate consideration of energy-dependent behavior for composite materials such as biological tissues, alloys, polymers, compounds, and mixtures for total and partial photon interactions. Later on, Un and Caner developed a new platform (Direct- Z eff ) that computes the mass attenuation coefficients, the effective atomic numbers and the effective electron numbers per unit mass for any selected materials [ 39 39 ]. XMuDat is also a program to be used with Windows 95 or Windows NT operating systems for the evaluation of mass attenuation coefficient, mass energy-absorption coefficient and the single-value of effective atomic number for a given material for total or specific reaction channels over range of 1 keV to 50 MeV [ 40 40 ]. Similar to the XMuDat, another program called NXcom is constructed for calculating attenuation coefficients of fast neutrons and gamma-rays [ 41 41 ]. The results of the program were tested by comparing them with the ones from WinXcom and Monte Carlo simulations, and shows very good agreement. There are numerous studies in literature (see, for instance, [ 42 42 - 49 49 42-49 ]) providing a wide database for attenuation coefficients, Z eff , N eff and Kerma for biomedical materials by using various methods and programs. The present work shows the comparison between results from Geant4 simulation package, theoretical calculations and experiments for three different types of human tissues in terms of mass attenuation coefficients and effective atomic number. For this purpose, first of all, Geant4 Monte Carlo simulations have been utilized to extract the mass attenuation coefficients of photon for adipose, muscle and bone tissues. The results have been compared with the ones from theoretical calculations and experimental data available in the literature. Effective atomic numbers of studied tissue models have been computed by direct method, Auto- Z eff and XMuDat programs. In addition, mean free paths of the photons, Kerma relative to air and effective electron density for those tissues have been discussed. 2 Geant4 Simulation 1 Monte Carlo simulations of this study has been performed by using the package of Geant4.10.03.p02. Geant4 is the successor of the GEANT series of software toolkits that were first to use object-oriented programming based on C++ language. It allows users to define classes for the detector geometry, primary particle generator, and physics processes to handle the interactions of particles with matter. It provides a set of electromagnetic physics processes handling the photon-tissue interactions (Photo Electric Effect, Compton Scattering, Rayleigh Scattering and pair production) over a wide range of energy [ 50 50 ]. In this study, adipose, muscle and bone tissue models have been studied by the Geant4 code that was explained in detail in our previous work [ 51 51 ]. Briefly, they have been considered to have simple geometries, slabs of various thicknesses from 1 cm to 10 cm. The elemental composition of the studied tissues are shown in Table 1 Table 1 . The densities of tissue models have been chosen to be 0.95 g/cm 3 for adipose tissue, 1.04 g/cm 3 for muscle tissue and 1.85 g/cm 3 for bone tissue [ 1 1 , 3 3 ]. 2 1 Element 3 1 Tissue 1 1 Z 1 1 Symbol 1 1 Adipose 1 1 Muscle 1 1 Bone 1 1 1 1 1 H 1 1 0.1140 1 1 0.1020 1 1 0.0640 1 1 6 1 1 C 1 1 0.5980 1 1 0.1230 1 1 0.2780 1 1 7 1 1 N 1 1 0.0070 1 1 0.0350 1 1 0.0270 1 1 8 1 1 O 1 1 0.2780 1 1 0.7290 1 1 0.4100 1 1 11 1 1 Na 1 1 0.0010 1 1 0.0008 1 1 - 1 1 12 1 1 Mg 1 1 - 1 1 0.0002 1 1 0.0020 1 1 15 1 1 P 1 1 - 1 1 0.0020 1 1 0.0700 1 1 16 1 1 S 1 1 0.0010 1 1 0.0050 1 1 0.0020 1 1 17 1 1 Cl 1 1 0.0010 1 1 - 1 1 - 1 1 19 1 1 K 1 1 - 1 1 0.0030 1 1 - 1 1 20 1 1 Ca 1 1 - 1 1 - 1 1 0.1470 The beam source has been positioned on one of the edges of the tissue model such that the momentum directions of the photons are perpendicular to the tissue surface. Simulations have been performed for 14 different photon energies from 20 keV to 50 MeV, for each tissue type seperately. For each run, 10 6 incident photons have been injected upon the tissue models. 3 Methods of Calculation 1 3.1 Mass Attenuation Coefficient 2 A beam of photons with an incident intensity, I 0 , penetrating a layer of material with thickness, x , and density, ρ , emerges with intensity, I , given by the exponential attenuation law: <math display="block" id="M1"><mrow><mi>I</mi><mo>/</mo><msub><mi>I</mi><mn>0</mn></msub><mo>=</mo><msup><mi>e</mi><mrow><mo>−</mo><mi>μ</mi><mi>x</mi></mrow></msup></mrow></math> where μ in the unit of cm -1 is the linear attenuation coefficient of the absorber, which can be obtained from measured values of I 0 , I and x . In addition, the ratio of the linear attenuation coefficient to the density ( μ/ρ ) is called the mass attenution coefficient. μ/ρ relies heavily on theoretical values for the total cross section per atom ( σ tot ) according to the equation: <math display="block" id="M2"><mrow><mi>μ</mi><mo>/</mo><mi>ρ</mi><mo>=</mo><msub><mi>σ</mi><mrow><mtext>tot</mtext></mrow></msub><mo>/</mo><mi>u</mi><mi>A</mi></mrow></math> where u is the atomic mass unit and A is the atomic mass number of the target element. The total cross section can be written as the sum over contributions from the principal photon interactions by <math display="block" id="M3"><mrow><msub><mi>σ</mi><mrow><mtext>tot</mtext></mrow></msub><mo>=</mo><msub><mi>σ</mi><mrow><mtext>ph</mtext></mrow></msub><mo>+</mo><msub><mi>σ</mi><mrow><mtext>coh</mtext></mrow></msub><mo>+</mo><msub><mi>σ</mi><mrow><mtext>incoh</mtext></mrow></msub><mo>+</mo><msub><mi>σ</mi><mrow><mtext>pair</mtext></mrow></msub></mrow></math> where σ ph is the atomic photo-effect cross section, σ coh and σ incoh are the coherent (Rayleigh) and the incoherent (Compton) scattering cross sections, respectively, and σ pair is the cross sections for electron-positron production. The total mass attenuation coefficient, ( μ / ρ ) tot for any chemical compound or mixture of elements is given by mixture rule. <math display="block" id="M4"><mrow><msub><mrow><mo stretchy="false">(</mo><mi>μ</mi><mo>/</mo><mi>ρ</mi><mo stretchy="false">)</mo></mrow><mrow><mtext>tot</mtext></mrow></msub><mo>=</mo><mstyle displaystyle="true"><munder><mo>∑</mo><mi>i</mi></munder></mstyle><msub><mi>w</mi><mi>i</mi></msub><msub><mrow><mo stretchy="false">(</mo><mi>μ</mi><mo>/</mo><mi>ρ</mi><mo stretchy="false">)</mo></mrow><mi>i</mi></msub></mrow></math> w i and ( μ / ρ ) i are the weight fraction and mass attenuation coefficient of the i th constituent element, respectively. 3.2 Effective Atomic Number 2 Z eff has been calculated as a function of photon energy by using three different methods: direct method, interpolation method and Auto- Z eff program. In addition, XMuDat software has also been utilized to get Z eff value. According to the direct method, Z eff is equal to ratio of total atomic cross-section ( σ a ) and total electronic cross-section ( σ e ). <math display="block" id="M5"><mrow><msub><mi>Z</mi><mrow><mtext>eff</mtext></mrow></msub><mo>=</mo><msub><mi>σ</mi><mtext>a</mtext></msub><mo>/</mo><msub><mi>σ</mi><mtext>e</mtext></msub></mrow></math> The total atomic cross section, σ a , is given by: <math display="block" id="M6"><mrow><msub><mi>σ</mi><mtext>a</mtext></msub><mo>=</mo><mfrac><mrow><msub><mrow><mo stretchy="false">(</mo><mi>μ</mi><mo>/</mo><mi>ρ</mi><mo stretchy="false">)</mo></mrow><mrow><mtext>tissue</mtext></mrow></msub></mrow><mrow><msub><mi>N</mi><mtext>A</mtext></msub><mstyle displaystyle="true"><munder><mo>∑</mo><mi>i</mi></munder></mstyle><mo stretchy="false">(</mo><mfrac><mrow><msub><mi>w</mi><mi>i</mi></msub></mrow><mrow><msub><mi>A</mi><mi>i</mi></msub></mrow></mfrac><mo stretchy="false">)</mo></mrow></mfrac><mn>,</mn></mrow></math> where N A is the Avogadro’s number. For calculations, the mass attenuation coefficients have been accepted to have values estimated based on Geant4 simulation. The electronic cross section, σ e , is defined as: <math display="block" id="M7"><mrow><msub><mi>σ</mi><mtext>e</mtext></msub><mo>=</mo><mfrac><mn>1</mn><mrow><msub><mi>N</mi><mtext>A</mtext></msub></mrow></mfrac><mstyle displaystyle="true"><munder><mo>∑</mo><mi>i</mi></munder></mstyle><mo stretchy="false">[</mo><mfrac><mrow><msub><mi>f</mi><mi>i</mi></msub><msub><mi>A</mi><mi>i</mi></msub></mrow><mrow><msub><mi>Z</mi><mi>i</mi></msub></mrow></mfrac><msub><mrow><mo stretchy="false">(</mo><mi>μ</mi><mo>/</mo><mi>ρ</mi><mo stretchy="false">)</mo></mrow><mi>i</mi></msub><mo stretchy="false">]</mo><mn>,</mn></mrow></math> where f i and Z i are the fractional abundance (mass fraction) and the atomic number of the i th element in the tissue, respectively. The effective electron density ( N eff ) is the number of electrons per unit mass and expressed by the relation of: <math display="block" id="M8"><mrow><msub><mi>N</mi><mrow><mtext>eff</mtext></mrow></msub><mo>=</mo><mfrac><mrow><msub><mrow><mo stretchy="false">(</mo><mi>μ</mi><mo>/</mo><mi>ρ</mi><mo stretchy="false">)</mo></mrow><mrow><mtext>tissue</mtext></mrow></msub></mrow><mrow><msub><mi>σ</mi><mtext>e</mtext></msub></mrow></mfrac><mn>.</mn></mrow></math> The Auto- Z eff program is an user friendly interface used to compute energy dependent atomic numbers with very low uncertainties for user-specified material. The cross section matrices of photon-element interactions are constructed for energies between 0.01 and 1000 MeV. Coefficients for composite media are constructed via linear additivity of the fractional constituents and contrasted against the precalculated matrices at each energy, thereby associating an effective atomic number through interpolation of adjacent cross section data. XMuDat computer software is able to produce a single-valued effective atomic number for a given material by assuming the photoelectric effect as the main interaction process [ 40 40 ]. The method of XMuDat program to calculate the effective atomic number based on the formula: <math display="block" id="M9"><mrow><msub><mi>Z</mi><mrow><mtext>eff</mtext></mrow></msub><mo>=</mo><mstyle displaystyle="true"><munder><mo>∑</mo><mi>i</mi></munder></mstyle><msup><mrow><mo stretchy="false">(</mo><msub><mi>α</mi><mi>i</mi></msub><msub><mi>Z</mi><mi>i</mi></msub><msup><mrow/><mrow><mi>m</mi><mo>−</mo><mn>1</mn></mrow></msup><mo stretchy="false">)</mo></mrow><mrow><mn>1</mn><mo>/</mo><mo stretchy="false">(</mo><mi>m</mi><mo>−</mo><mn>1</mn><mo stretchy="false">)</mo></mrow></msup><mn>,</mn></mrow></math> where α i is the fractional number of the electrons of the i th element and m is a constant between 3 and 5. It is preferred that m is set to be 3.6 for elements with Z lt;6 and 4.1 for elements with Z gt;6 [ 52 52 ]. 3.3 Kerma 2 Kerma is the initial kinetic energy of all secondary charged particles liberated per unit mass. It is applicable to photons and neutrons and in the unit of J kg -1 = Gy, as the absorbed dose. Kerma of the compounds, mixtures, and human tissues relative to air was calculated [ 47 47 , 53 53 , 54 54 ] based on the relation: <math display="block" id="M10"><mrow><msub><mi>K</mi><mtext>r</mtext></msub><mo>=</mo><mfrac><mrow><msub><mi>K</mi><mrow><mtext>tissue</mtext></mrow></msub></mrow><mrow><msub><mi>K</mi><mrow><mtext>air</mtext></mrow></msub></mrow></mfrac><mo>=</mo><mfrac><mrow><msub><mrow><mo stretchy="false">(</mo><msub><mi>μ</mi><mrow><mtext>en</mtext></mrow></msub><mo>/</mo><mi>ρ</mi><mo stretchy="false">)</mo></mrow><mrow><mtext>tissue</mtext></mrow></msub></mrow><mrow><msub><mrow><mo stretchy="false">(</mo><msub><mi>μ</mi><mrow><mtext>en</mtext></mrow></msub><mo>/</mo><mi>ρ</mi><mo stretchy="false">)</mo></mrow><mrow><mtext>air</mtext></mrow></msub></mrow></mfrac><mn>,</mn></mrow></math> where ( μ en / ρ ) is the mass energy-absorption coefficient. In this work, the values of ( μ en / ρ ) for the elements of studied tissues have been extracted from the compilation of Hubbell [ 9 9 ], and then the values of ( μ en / ρ ) for the tissues have been calculated by mixture rule. 4 Result and Discussion 1 4.1 Mass attenuation coefficient 2 The exponential relationship between the ratio of transmitted photon intensity to that of the incident photons ( I/I 0 ), which is called trasmittance, and the tissue thickness is given by Eq. 1. Geant4 simulation results I/I 0 with energies from 0.02 MeV to 50 MeV have been plotted as functions of tissue thicknesses for adipose, muscle and bone tissue models seperately (see Figure 1 Figure 1 ). Each graph has been fitted to the function of e - μx in order to get the fitting parameter, μ (linear attenuation coefficient). The solid lines in the graphs are the fit curves. The mass attenuation coefficient ( μ/ρ ) for each tissue type and energy has been calculated by dividing μ to the density of the tissue of interest. μ/ρ values for adipose, muscle and bone tissues obtained from the simulation are given in Tables 2 Tables 2 , 3 3 and 4 4 , respectively, together with the ones from theoretical calculations and experimental studies. 1 1 Photon energy (MeV) 3 1 Mass attenuation coefficients(cm 2 =g) 1 1 1 1 Theoretical 1 1 Simulation 1 1 Experiment 1 1 0.02 1 1 1.079 1 1 0.939 1 1 1.26 [ 33 33 ] 1 1 0.03 1 1 0.306 1 1 0.264 1 1 0.32 [ 33 33 ] 1 1 0.05 1 1 0.223 1 1 0.205 1 1 0.21 [ 33 33 ] 1 1 0.06 1 1 0.198 1 1 0.187 1 1 - 1 1 0.08 1 1 0.179 1 1 0.173 1 1 0.18 [ 34 34 ] 1 1 0.15 1 1 0.150 1 1 0.148 1 1 - 1 1 0.36 1 1 0.111 1 1 0.111 1 1 - 1 1 0.66 1 1 0.086 1 1 0.088 1 1 0.08 [ 34 34 ] 1 1 0.83 1 1 0.077 1 1 0.076 1 1 0.07 [ 34 34 ] 1 1 1.17 1 1 0.065 1 1 0.066 1 1 0.06 [ 34 34 ] 1 1 1.33 1 1 0.061 1 1 0.062 1 1 0.05 [ 34 34 ] 1 1 10 1 1 0.021 1 1 0.020 1 1 - 1 1 20 1 1 0.017 1 1 0.014 1 1 - 1 1 50 1 1 0.015 1 1 0.017 1 1 - 1 1 Photon energy (MeV) 3 1 Mass attenuation coefficients(cm 2 =g) 1 1 1 1 Theoretical 1 1 Simulation 1 1 Experiment 1 1 0.02 1 1 1.702 1 1 1.517 1 1 1.75 [ 33 33 ] 1 1 0.03 1 1 0.379 1 1 0.325 1 1 - 1 1 0.05 1 1 0.244 1 1 0.221 1 1 0.30 [ 33 33 ] 1 1 0.06 1 1 0.206 1 1 0.192 1 1 0.19 [ 33 33 ] 1 1 0.08 1 1 0.182 1 1 0.174 1 1 0.17 [ 21 21 ] 1 1 0.15 1 1 0.149 1 1 0.146 1 1 - 1 1 0.36 1 1 0.110 1 1 0.109 1 1 0.09 [ 21 21 ] 1 1 0.66 1 1 0.085 1 1 0.085 1 1 0.08 [ 21 21 ] 1 1 0.83 1 1 0.076 1 1 0.078 1 1 - 1 1 1.17 1 1 0.065 1 1 0.065 1 1 0.05 [ 21 21 ] 1 1 1.33 1 1 0.061 1 1 0.062 1 1 0.05 [ 21 21 ] 1 1 10 1 1 0.022 1 1 0.025 1 1 - 1 1 20 1 1 0.018 1 1 0.017 1 1 - 1 1 50 1 1 0.017 1 1 0.012 1 1 - 1 1 Photon energy (MeV) 3 1 Mass attenuation coefficients(cm 2 =g) 1 1 1 1 Theoretical 1 1 Simulation 1 1 Experiment 1 1 0.02 1 1 6.325 1 1 - 1 1 - 1 1 0.03 1 1 0.982 1 1 0.910 1 1 - 1 1 0.05 1 1 0.419 1 1 0.379 1 1 0.49 [ 33 33 ] 1 1 0.06 1 1 0.278 1 1 0.255 1 1 0.24 [ 33 33 ] 1 1 0.08 1 1 0.207 1 1 0.193 1 1 0.19 [ 21 21 ] 1 1 0.15 1 1 0.149 1 1 0.144 1 1 - 1 1 0.36 1 1 0.107 1 1 0.106 1 1 0.09 [ 21 21 ] 1 1 0.66 1 1 0.082 1 1 0.082 1 1 0.08 [ 21 21 ] 1 1 0.83 1 1 0.074 1 1 0.074 1 1 0.07 [ 21 21 ] 1 1 1.17 1 1 0.062 1 1 0.062 1 1 0.06 [ 21 21 ] 1 1 1.33 1 1 0.058 1 1 0.058 1 1 0.05 [ 21 21 ] 1 1 10 1 1 0.023 1 1 0.025 1 1 - 1 1 20 1 1 0.20 1 1 0.021 1 1 - 1 1 50 1 1 0.019 1 1 0.019 1 1 - It can be seen from the tables that the results obtained from simulation, theoretical and experimental studies are in a very good agreement. Generally, the errors on the values from the simulations and the experiments were calculated to be less than 2% and 3%, respectively [ 21 21 , 33 33 ]. In the present study, error propogation of the fitting in the figures have been found to be less than 5%. The small discrepancies among the values can be attributed to the tissue modeling and the methods that were used to estimate the cross sections. The μ/ρ values decrease with increasing photon energies for all tissue types. This behavior depends on the total photon interaction that shows the superposition of multiple partial interactions (e.g. photoelectric absorption, Compton scattering and pair production). According to theoretical calculations, the mass attenuation coefficients experience a significant attenuation of 71.64 % for adipose tissue, 77.73 % for muscle tissue and 84.47 % for bone tissue in the energy region between 20 and 30 keV. The calculated theoretical μ/ρ values tend to be higher than the simulated values for low energies. However, these values show more agreement at a high energy region. It is clear that the higher theoretical values reflect the effect of the chemical composition of the tissue models and the mixture rule method. The photoelectric absorption is a dominant process for photon energies of less than 50 keV. In this energy region, the value of the mass attenuation coefficients rapidly goes down due to the photoelectric cross sections directly proportional to Z 4 and inversely to E 3 . At higher energies, the Compton scattering becomes the dominant mechanism that varies with Z and E -1 . Therefore, a smooth decrease of μ/ρ values was observed. For the photon energies greater than 1022 keV, the pair-production process appears and the interaction cross-section changes with Z 2 and log E . Thus, the μ/ρ values look independent of photon energy and decrease slowly. Mean free path ( λ ), defined as the average distance between two successive interactions, is simply reciprocal of the linear attenuation coefficient. λ of the photons through adipose, muscle and bone tissues are illustrated in Figure 2 Figure 2 as a function of energy. It is clear that λ values increase with energy except for the abrupt changes near the absorption edges and at the threshold energy of photon interactions such as: photoelectric effect, coherent (Rayleigh), incoherent scattering, and electron-positron production. The λ is inversely proportional with the number of photon collisions. Thus, it can be concluded that the photon collisions are relatively high in the low energy region lt; 30 keV, which makes the mass attenuation coefficients bigger. This property is also an explenation for the fact that the low energy X-rays deposit much energy (radiation dose) in the first few centimeters in the body. For higher energies, the number of collisions sharply drop down causing a decrease in the mass attenuation coefficients. 4.2 Effective atomic numbers 2 The effective atomic numbers ( Z eff ) for adipose, muscle and bone tissues in the photon energy region 0.01 to 20 MeV are shown in Fig. 3 Fig. 3 . It was found that the computed effective atomic numbers at a low energy region (lt; 30 keV) have a high uncertainty. For example, it reaches the value of ± 25% when Auto- Z eff is utilized [ 38 38 ]. On the other hand, the effective atomic numbers calculated by using the Auto- Z eff software and direct method are in a very good agreement in the energy region from 50 keV to 10 MeV where the Compton scattering is a dominant process. In this medium-energy region, the effective atomic numbers have been observed to be constant, whereas significant variations have been found in the lower (lt; 50 keV) as well as in the higher (gt; 10 MeV) regions. The effective atomic numbers calculated by using the direct method have higher values than those computed by Auto- Z eff software and interpolation method in the energy regions of both photo-absorption and pair-production. Moreover, the Z eff calculated by direct method and Auto- Z eff program have a maximum differences of 18.9% for adipose, 21.2% for muscle and 38.7% for bone in the energy region around 10 keV. The values of Z eff for adipose and muscle by using interpolation method were taken from [ 44 44 ]. Furthermore, the single Z eff values obtained from the XMuDat program are also included to the graphs shown in Figure 3 Figure 3 . For each type of tissue, it lies between the values obtained from the direct method and the one from Auto- Z eff software. The plot of effective atomic numbers versus effective electron densities for tissue models is shown in Fig. 4 Fig. 4 . As it can be expected, the relation between N eff and Z eff is smooth and linear since they are directly proportional with each other. 4.3 Kerma 2 Kerma relative to air ( K r ) as a function of photon energy for adipose, muscle and bone tissues over the energy range from 1 keV to 20 MeV, is shown in Fig. 5 Fig. 5 . Energy-dependent variation of K r represents the behavior of Z eff due to different processes such as photo-absorption, Compton effect and pair-production. Although this variation can be seen clearly at a photon energy below 100 keV, the K r values were found to be constant ( ∼ unity) for higher energies. The K r value for bone tissue, which contains the higher- Z elements like Ca, reaches up to 6.96 at 40 keV. This peak value of K r can be attributed to the photoelectric effect around the K-absorption edge of Ca. For muscle and adipose tissues, which have negligible fraction weights of high- Z elements, the values of K r are equal or less than one. 5 Conclusion 1 This article provides data, from Geant4 simulations and theoretical calculations, of mass attenuation coefficients and effective atomic numbers for adipose, muscle and bone tissues over a wide range of energies up to 50 MeV. All of these parameters are mainly dependent on the incident photon energy and the tissue constituent elements. The mass attenuation coefficient and effective atomic numbers sharply decrease where the photo-absorption is the dominant mechanism in tissue. Also, they are nearly constant at the predominance of the Compton scattering, and a smooth variation can be seen at the energy region of the pair-production process. Mass attenuation coefficients generated by using Geant4 have uncertainties less than ±1% at the low energy region and less than 8% at high energies. Effective atomic numbers calculated by direct method and Auto- Z eff program have a maximum difference of 18.9% for adipose, 21.2% for muscle and 38.7% for bone in the energy region around 10 keV. The single effective atomic number obtained by the XMuDat program exhibited values higher than those of the Auto- Z eff program and less than those of the direct method. Moreover, the elements that compose the tissue play a significant role in the calculation of Kerma. References ICRP , Basic anatomical and physiological data for use in radiological protection: reference values , ICRP Publication 89 . Ann. ICRP 32 ( 2003 ). Basic anatomical and physiological data for use in radiological protection: reference values Y.S. Kim , Human Tissues: Chemical Composition and Photon Dosimetry Data , Radiat. Res . 57 ( 1 ), 38 - 45 ( 1974 ). doi: 10.2307/3573753 http://doi.org/10.2307/3573753 Human Tissues: Chemical Composition and Photon Dosimetry Data D.R. White J. Booz , R.V. Griffith et al ., ICRU Report 44: Tissue Substitutes in Radiation Dosimetry and Measurement , Journal of the ICRU os 23 ( 1989 ). doi: 10.1093/jicru/os23.1.Report44 http://doi.org/10.1093/jicru/os23.1.Report44 ICRU Report 44: Tissue Substitutes in Radiation Dosimetry and Measurement Y. Elmahroug , B. Tellili , C. Souga et al ., Determination of total mass attenuation coefficients, effective atomic numbers and electron densities for different shielding materials , Ann. Nucl. Energy 75 , 268 - 274 ( 2015 ). doi: 10.1016/j.anucene.2014.08.015 http://doi.org/10.1016/j.anucene.2014.08.015 Determination of total mass attenuation coefficients, effective atomic numbers and electron densities for different shielding materials A. McNair , ICRU Report 33: Radiation Quantities and Units , J. Label Compd. Radiopharm . 18 ( 1981 ). doi: 10.1002/jlcr.2580180918 http://doi.org/10.1002/jlcr.2580180918 ICRU Report 33: Radiation Quantities and Units R.T. Berger , The X- or Gamma-Ray Energy Absorption or Transfer Coefficient: Tabulations and Discussion , Radiat. Res . 15 , 1 - 29 ( 1961 ). doi: 10.2307/3571063 http://doi.org/10.2307/3571063 The X- or Gamma-Ray Energy Absorption or Transfer Coefficient: Tabulations and Discussion J.W. Allison , Gamma-Radiation Absorption Coefficients of Various Materials Allowing for Bremsstrahlung and other Secondary Radiations , Austral. J. Phys 14 , 443 - 461 ( 1961 ). doi: 10.1071/PH610443 http://doi.org/10.1071/PH610443 Gamma-Radiation Absorption Coefficients of Various Materials Allowing for Bremsstrahlung and other Secondary Radiations J.H. Hubbell , Photon Mass Attenuation and Mass Energy-Absorption Coefficients for H, C, N, O, Ar, and Seven Mixtures from 0.1 keV to 20 MeV , Radiat. Res ., 70 , 58 - 81 ( 1977 ). doi: 10.2307/3574732 http://doi.org/10.2307/3574732 Photon Mass Attenuation and Mass Energy-Absorption Coefficients for H, C, N, O, Ar, and Seven Mixtures from 0.1 keV to 20 MeV J.H. Hubbell , Photon mass attenuation and energy-absorption coefficients , Int. J. Appl. Radiat. Isot . 33 ( 11 ), 1269 - 1290 ( 1982 ). doi: 10.1016/0020-708X(82)90248-4 http://doi.org/10.1016/0020-708X(82)90248-4 Photon mass attenuation and energy-absorption coefficients J.H. Hubbell , S.M. Seltzer , Tables of X-Ray mass attenuation coefficients and mass energy-absorption coefficients 1 keV to 20 MeV for Elements Z=1 to 92 and 48 Additional substances of Dosimetric Interest , NIST Standard Reference Database 126 ( 1995 ). doi: 10.18434/T4D01F http://doi.org/10.18434/T4D01F Tables of X-Ray mass attenuation coefficients and mass energy-absorption coefficients 1 keV to 20 MeV for Elements Z=1 to 92 and 48 Additional substances of Dosimetric Interest M.J. Berger , J.H. Hubbell , S.M. Seltzer et al ., XCOM: photon cross sections database , NIST Standard Reference Database 8 ( 1998 ). doi: 10.18434/T48G6X http://doi.org/10.18434/T48G6X XCOM: photon cross sections database L. Gerward , N. Guilbert , K. B. Jensen et al ., WinXCom - a program for calculating X-ray attenuation coefficients , Radiat. Phys. Chem . 71 , ( 2004 ) 653 - 654 (2004) . doi: 10.1016/j.radphyschem.2004.04.040 http://doi.org/10.1016/j.radphyschem.2004.04.040 WinXCom - a program for calculating X-ray attenuation coefficients C.T. Chantler , Theoretical form factor, attenuation, and scattering tabulation for Z= 1–92 from E= 1–10 eV to E= 0.4–1.0 MeV , J. Phys. Chem. Ref. Data , 24 ( 1 ), 71 - 643 ( 1995 ). doi: 10.1063/1.555974 http://doi.org/10.1063/1.555974 Theoretical form factor, attenuation, and scattering tabulation for Z= 1–92 from E= 1–10 eV to E= 0.4–1.0 MeV C.T. Chantler et al ., X-Ray Form Factor, Attenuation and Scattering Tables , NIST Standard Reference Database 66 ( 2005 ). doi: 10.18434/T4HS32 http://doi.org/10.18434/T4HS32 X-Ray Form Factor, Attenuation and Scattering Tables S. Agostinelli , J. Allison , K. Amako et al ., GEANT4 - A simulation toolkit , Nucl. Instrum. Methods Phys. Res. A , 506 ( 3 ), 250 - 303 ( 2003 ). doi: 10.1016/S0168-9002(03)01368-8 http://doi.org/10.1016/S0168-9002(03)01368-8 GEANT4 - A simulation toolkit M.E. Medhat , Y. Wang , Geant4 code for simulation attenuation of gamma rays through scintillation detectors , Ann. Nucl. Energy , 62 , 316 - 320 ( 2013 ). doi: 10.1016/j.anucene.2013.06.034 http://doi.org/10.1016/j.anucene.2013.06.034 Geant4 code for simulation attenuation of gamma rays through scintillation detectors V.P. Singh , M.E. Medhat N.M. Badiger , Photon attenuation coefficients of thermoluminescent dosimetric materials by Geant4 toolkit, XCOM program and experimental data: A comparison study , Ann. Nucl. Energy , 68 , 96 - 100 ( 2014 ). doi: 10.1016/j.anucene.2014.01.011 http://doi.org/10.1016/j.anucene.2014.01.011 Photon attenuation coefficients of thermoluminescent dosimetric materials by Geant4 toolkit, XCOM program and experimental data: A comparison study S.S. Obaid , M.I. Sayyed , D.K. Gaikwad et al ., Photon attenuation coefficients of different rock samples using MCNPX, Geant4 simulation codes and experimental results: a comparison study , Radiat. Eff. Defect S ., 173 ( 11-12 ), 900 - 914 ( 2018 ). doi: 10.1080/10420150.2018.1505890 http://doi.org/10.1080/10420150.2018.1505890 Photon attenuation coefficients of different rock samples using MCNPX, Geant4 simulation codes and experimental results: a comparison study A. Kumar S.P. Singh , Y. Elmahroug et al ., Gamma ray shielding studies on 26.66 B2O3–16GeO2–4Bi2O3–(53.33–x) PbO–xPbF2 glass system using MCNPX, Geant4 and XCOM , Mater. Res. Express , 5 ( 9 ), 095203 ( 2018 ). doi: 10.1088/2053-1591/aad821 http://doi.org/10.1088/2053-1591/aad821 Gamma ray shielding studies on 26.66 B2O3–16GeO2–4Bi2O3–(53.33–x) PbO–xPbF2 glass system using MCNPX, Geant4 and XCOM R.M. Lokhande , B.S. Surung , P.P. Pawar , Measurement of effective atomic number and electron density of carbohydrates by using NIST, Geant4 and NaI(Tl): A comparative study , Int. J. of Adv. Res . 5 ( 5 ), 1733 - 1740 ( 2017 ). doi: 10.21474/IJAR01/4303 http://doi.org/10.21474/IJAR01/4303 Measurement of effective atomic number and electron density of carbohydrates by using NIST, Geant4 and NaI(Tl): A comparative study M.E. Medhat , S.P. Shirmardi , V.P. Singh , Comparison of Geant 4, MCNP Simulation Codes of Studying Attenuation of Gamma rays through Biological Materials with XCOM and Experimental Data , J. Appl. Computat. Math . 3 ( 6 ), 1000179 ( 2014 ). doi: 10.4172/2168-9679.1000179 http://doi.org/10.4172/2168-9679.1000179 Comparison of Geant 4, MCNP Simulation Codes of Studying Attenuation of Gamma rays through Biological Materials with XCOM and Experimental Data M.I. Sayyed , H.O. Tekin , E.E. Altunsoy et al ., Radiation shielding study of tellurite tungsten glasses with different antimony oxide as transparent shielding materials using MCNPX code , J. Non-Cryst. Solids , 498 , 167 - 172 ( 2018 ). doi: 10.1016/j.jnoncrysol.2018.06.022 http://doi.org/10.1016/j.jnoncrysol.2018.06.022 Radiation shielding study of tellurite tungsten glasses with different antimony oxide as transparent shielding materials using MCNPX code B.O. Elbashir , M.G. Dong , M.I. Sayyed et al ., Comparison of Monte Carlo simulation of gamma ray attenuation coefficients of amino acids with XCOM program and experimental data , Results Phys ., 9 , 6 - 11 ( 2018 ). doi: 10.1016/j.rinp.2018.01.075 http://doi.org/10.1016/j.rinp.2018.01.075 Comparison of Monte Carlo simulation of gamma ray attenuation coefficients of amino acids with XCOM program and experimental data M.I. Sayyed , S.A.M. Issa , M. Büyükyıldız et al ., Determination of nuclear radiation shielding properties of some tellurite glasses using MCNP5 code , Radiat. Phys. Chem ., 150 , 1 - 8 ( 2018 ). doi: 10.1016/j.radphyschem.2018.04.014 http://doi.org/10.1016/j.radphyschem.2018.04.014 Determination of nuclear radiation shielding properties of some tellurite glasses using MCNP5 code K. Verdipoor , A. Alemi , A. Mesbahi , Photon mass attenuation coefficients of a silicon resin loaded with WO3, PbO, and Bi2O3 Micro and Nano-particles for radiation shielding , Radiat. Phys. Chem ., 147 , 85 - 90 ( 2018 ). doi: 10.1016/j.radphyschem.2018.02.017 http://doi.org/10.1016/j.radphyschem.2018.02.017 Photon mass attenuation coefficients of a silicon resin loaded with WO3, PbO, and Bi2O3 Micro and Nano-particles for radiation shielding A. Mesbahi , H. Ghiasi , Shielding properties of the ordinary concrete loaded with micro- and nano-particles against neutron and gamma radiations , Appl. Radiat. Isot ., 136 , 27 - 31 ( 2018 ). doi: 10.1016/j.apradiso.2018.02.004 http://doi.org/10.1016/j.apradiso.2018.02.004 Shielding properties of the ordinary concrete loaded with micro- and nano-particles against neutron and gamma radiations G.J. Hine , The effective atomic numbers of materials for various gamma ray interactions , Phys. Rev , 85 , 725 - 737 ( 1952 ). The effective atomic numbers of materials for various gamma ray interactions M.T. Islam , N.A. Rae , J.L. Glover et al ., Measurement of the x-ray mass attenuation coefficients of gold in the 38 – 50-keV energy range , Phys. Rev. A 81 022903 ( 2010 ). doi: 10.1103/PhysRevA.81.022903 http://doi.org/10.1103/PhysRevA.81.022903 Measurement of the x-ray mass attenuation coefficients of gold in the 38 – 50-keV energy range B. Goswami , N. Chaudhuri , Measurements of gamma-ray attenuation coefficients , Phys. Rev. A 7 1912 - 1916 ( 1973 ). doi: 10.1016/0029-554X(73)90358-3 http://doi.org/10.1016/0029-554X(73)90358-3 Measurements of gamma-ray attenuation coefficients B.S. Sidhu , A.S. Dhaliwal , K.S. Mann et al ., Study of mass attenuation coefficients, effective atomic numbers and electron densities for some low Z compounds of dosimetry interest at 59.54 keV incident photon energy , Ann. Nucl. Energy 42 , 153 - 157 ( 2012 ). doi: 10.1016/j.anucene.2011.12.015 http://doi.org/10.1016/j.anucene.2011.12.015 Study of mass attenuation coefficients, effective atomic numbers and electron densities for some low Z compounds of dosimetry interest at 59.54 keV incident photon energy H. Buhr , L. Büermann , M. Gerlach et al ., Measurement of the mass energy-absorption coefficient of air for x-rays in the range from 3 to 60keV , Phys. Med. Biol . 57 ( 24 ), 8231 - 8247 ( 2012 ). doi: 10.1088/0031-9155/57/24/8231 http://doi.org/10.1088/0031-9155/57/24/8231 Measurement of the mass energy-absorption coefficient of air for x-rays in the range from 3 to 60keV B. Akça , S.Z. Erzeneoğlu , The mass attenuation coefficients, electronic, atomic, and molecular cross sections, effective atomic numbers, and electron densities for compounds of some biomedically important elements at 59.5 keV , Sci. Technol. Nucl. Ins ., 901465 ( 2014 ). doi: 10.1155/2014/901465 http://doi.org/10.1155/2014/901465 The mass attenuation coefficients, electronic, atomic, and molecular cross sections, effective atomic numbers, and electron densities for compounds of some biomedically important elements at 59.5 keV W. Geraldelli , A. Tomal , M.E. Poletti , Characterization of tissue-equivalent materials through measurements of the linear attenuation coefficient and scattering profiles obtained with polyenergetic beams , IEEE Trans. Nucl. Sci . 60 ( 2 ), 566 - 571 ( 2013 ). doi: 10.1109/TNS.2013.2248382 http://doi.org/10.1109/TNS.2013.2248382 Characterization of tissue-equivalent materials through measurements of the linear attenuation coefficient and scattering profiles obtained with polyenergetic beams N.A.B. Amin , J. Zukhi , N.A. Kabir et al ., Determination of effective atomic numbers from mass attenuation coefficients of tissue-equivalent materials in the energy range 60 keV-1.33 MeV , J. Phys.: Conf. Ser . 851 012018 ( 2017 ). doi: 10.1088/1742-6596/851/1/012018 http://doi.org/10.1088/1742-6596/851/1/012018 Determination of effective atomic numbers from mass attenuation coefficients of tissue-equivalent materials in the energy range 60 keV-1.33 MeV C.A. Jayachandran , Calculated effective atomic number and Kerma values for tissue-equivalent and dosimetry materials , Phys. Med. Biol . 16 ( 4 ), 617 - 623 ( 1971 ). doi: 10.1088/0031-9155/16/4/005 http://doi.org/10.1088/0031-9155/16/4/005 Calculated effective atomic number and Kerma values for tissue-equivalent and dosimetry materials S.R. Manohara , S.M. Hanagodimath , K.S. Thind et al ., The effective atomic number revisited in the light of modern photon-interaction cross-section databases , Appl. Radiat. Isot . 68 ( 4–5 ), 784 - 787 ( 2010 ). doi: 10.1016/j.apradiso.2009.09.047 http://doi.org/10.1016/j.apradiso.2009.09.047 The effective atomic number revisited in the light of modern photon-interaction cross-section databases K.S. Mann , M. Kurudirek , G.S. Sidhu , Verification of dosimetric materials to be used as tissue-substitutes in radiological diagnosis , Appl. Radiat. Isot . 70 ( 4 ), 681 - 691 ( 2012 ). doi: 10.1016/j.apradiso.2011.12.008 http://doi.org/10.1016/j.apradiso.2011.12.008 Verification of dosimetric materials to be used as tissue-substitutes in radiological diagnosis M.L. Taylor , R.L. Smith , F. Dossing et al ., Robust calculation of effective atomic numbers: The Auto- Z eff software , Med. Phys . 39 ( 4 ), 1769 - 1778 ( 2012 ). doi: 10.1118/1.3689810 http://doi.org/10.1118/1.3689810 Robust calculation of effective atomic numbers: The Auto- Z eff software A. Un , T. Caner , The Direct-Zeff software for direct calculation of mass attenuation coefficient, effective atomic number and effective electron number , Ann. Nucl. Energy 65 , 158 - 165 ( 2014 ). doi: 101016/janucene201310041 http://doi.org/101016/janucene201310041 The Direct-Zeff software for direct calculation of mass attenuation coefficient, effective atomic number and effective electron number R. Nowotny , XMuDat: Photon attenuation data on PC, IAEA Report IAEA-NDS 195 ( 1998 ). XMuDat: Photon attenuation data on PC, IAEA Report IAEA-NDS 195 A.M. El-Khayatt , NXcom - A program for calculating attenuation coefficients of fast neutrons and gamma-rays , Ann. Nucl. Energy , 38 ( 1 ), 128 - 132 ( 2011 ). doi: 10.1016/j.anucene.2010.08.003 http://doi.org/10.1016/j.anucene.2010.08.003 NXcom - A program for calculating attenuation coefficients of fast neutrons and gamma-rays H.C. Manjunatha , B. Rudraswamy , Study of effective atomic number and electron density for tissues from human organs in the energy range of 1 keV–100 GeV , Health Phys . 104 ( 2 ), 158 - 162 ( 2013 ). doi: 10.1097/HP.0b013e31827132e3 http://doi.org/10.1097/HP.0b013e31827132e3 Study of effective atomic number and electron density for tissues from human organs in the energy range of 1 keV–100 GeV M. Kurudirek , Effective atomic numbers, water and tissue equivalence properties of human tissues, tissue equivalents and dosimetric materials for total electron interaction in the energy region 10keV–1GeV , Appl. Radiat. Isot ., 94 , 1 - 7 ( 2014 ). doi: 10.1016/j.apradiso.2014.07.002 http://doi.org/10.1016/j.apradiso.2014.07.002 Effective atomic numbers, water and tissue equivalence properties of human tissues, tissue equivalents and dosimetric materials for total electron interaction in the energy region 10keV–1GeV V.R. Shivaramu , Effective atomic number for photon energy absorption and photon attenuation of tissues from human organs , Med. Dosim . 27 , 1 - 9 ( 2002 ). doi: 10.1016/S0958-3947(01)00078-4 http://doi.org/10.1016/S0958-3947(01)00078-4 Effective atomic number for photon energy absorption and photon attenuation of tissues from human organs V.P. Singh , N.M. Badiger , N. Kucuk , Assessment of methods for estimation of effective atomic numbers of common human organ and tissue substitutes: waxes, plastics and polymers , Radioprotection 49 ( 2 ), 115 - 121 ( 2014 ). doi: 10.1051/radiopro/2013090 http://doi.org/10.1051/radiopro/2013090 Assessment of methods for estimation of effective atomic numbers of common human organ and tissue substitutes: waxes, plastics and polymers M. Kurudirek , T. Onaran , Calculation of effective atomic number and electron density of essential biomolecules for electron, proton, alpha particle and multi-energetic photon interactions , Radiat. Phys. Chem ., 112 , 125 - 138 ( 2015 ). doi: 10.1016/j.radphyschem.2015.03.034 http://doi.org/10.1016/j.radphyschem.2015.03.034 Calculation of effective atomic number and electron density of essential biomolecules for electron, proton, alpha particle and multi-energetic photon interactions D. Salehi , D. Sardari , M.S. Jozani , Investigation of some radiation shielding parameters in soft tissue , J. Radiat. Res. Appl. Sci . 8 ( 3 ), 439 - 445 ( 2015 ). doi: 10.1016/j.jrras.2015.03.004 http://doi.org/10.1016/j.jrras.2015.03.004 Investigation of some radiation shielding parameters in soft tissue M. Kurudirek , Effective atomic number of soft tissue, water and air for interaction of various hadrons, leptons and isotopes of hydrogen , Int. J. Radiat. Biol . 93 ( 12 ), 1299 - 1305 ( 2017 ). doi: 10.1080/09553002.2018.1388546 http://doi.org/10.1080/09553002.2018.1388546 Effective atomic number of soft tissue, water and air for interaction of various hadrons, leptons and isotopes of hydrogen D.K. Gaikwad , M.I. Sayyed , S.S. Obaid et al ., Gamma ray shielding properties of TeO2-ZnF2-As2O3-Sm2O3 glasses , J. Alloys Compd ., 765 , 451 - 458 ( 2018 ). doi: 10.1016/j.jallcom.2018.06.240 http://doi.org/10.1016/j.jallcom.2018.06.240 Gamma ray shielding properties of TeO2-ZnF2-As2O3-Sm2O3 glasses J. Apostolakis , A. Bagulya , S. Elles et al ., Validation and verification of Geant4 standard electromagnetic physics , J. Phys.: Conf. Ser . 219 , 032044 , ( 2010 ). doi: 10.1088/1742-6596/219/3/032044 http://doi.org/10.1088/1742-6596/219/3/032044 Validation and verification of Geant4 standard electromagnetic physics B.T. Tonguc , H. Arslan , M.S. Al-Buriahi , Studies on mass attenuation coefficients, effective atomic numbers and electron densities for some biomolecules , Radiat. Phys. Chem ., 153 , 86 - 91 ( 2018 ). doi: 10.1016/j.radphyschem.2018.08.025 http://doi.org/10.1016/j.radphyschem.2018.08.025 Studies on mass attenuation coefficients, effective atomic numbers and electron densities for some biomolecules D.F. Jackson , H.J. David , X-ray attenuation coefficients of elements and mixtures , Phys. Rep . 70 , 169 - 233 ( 1981 ). doi: 10.1016/0370-1573(81)90014-4 http://doi.org/10.1016/0370-1573(81)90014-4 X-ray attenuation coefficients of elements and mixtures S.R. Manohara , S.M. Hanagodimath , L. Gerward , Studies on effective atomic number, electron density and kerma for some fatty acids and carbohydrates , Phys. Med. Biol . 53 ( 20 ), N377 - 86 ( 2008 ). doi: 10.1088/0031-9155/53/20/N01 http://doi.org/10.1088/0031-9155/53/20/N01 Studies on effective atomic number, electron density and kerma for some fatty acids and carbohydrates D. Yilmaz , Y. Şahin , L. Demir , Studies on mass attenuation coefficient, mass energy absorption coefficient, and kerma for Fe alloys at photon energies of 17.44 to 51.70 keV , Turk. J. Phys . 39 ( 1 ), 81 - 90 ( 2015 ). doi: 10.3906/fiz-1408-4 http://doi.org/10.3906/fiz-1408-4 Studies on mass attenuation coefficient, mass energy absorption coefficient, and kerma for Fe alloys at photon energies of 17.44 to 51.70 keV The auhtor would like to thank Mohammed Al-Buriahi for his contributions in theoretical calculations. 10.1007/s41365-019-0617-z 2018-09-28 2018-12-10 2018-12-28 2019-05-15 2019-06-01 2021-04-10 © China Science Publishing & Media Ltd. (Science Press), Shanghai Institute of Applied Physics, the Chinese Academy of Sciences, ChineseNuclear Society and Springer Nature Singapore Pte Ltd. 2019 This is a post-peer-review, pre-copyedit version of an article published in Nuclear Science and Techniques. The final authenticated version is available online at: http://dx.doi.org/10.1007/s41365-019-0617-z http://dx.doi.org/10.1007/s41365-019-0617-z . 2019 Nuclear Science and Techniques 06 30