Calculation of photon shielding properties for some neutron shielding materials

Calculation of photon shielding properties for some neutron shielding materials 1 2 corresp first-author Corresponding author. E-mail addresses: Ahmed_el_khayatt@yahoo.com, Corresponding author. E-mail addresses: Ahmed_el_khayatt@yahoo.com, Ahmed_el_khayatt@yahoo.com Physics Department, College of Science, Ibn Saud University (ISU), Saudi Arabia Physics Department, College of Science, Ibn Saud University (ISU), Saudi Arabia Reactor Physics Department, NRC, Atomic Energy Authority, 13759 Cairo, Egypt Reactor Physics Department, NRC, Atomic Energy Authority, 13759 Cairo, Egypt The objective of the present study is to calculate photon shielding parameters for seven polyethylene-based neutron shielding materials. The parameters include the effective atomic number ( Z eff ), the effective electron density ( N eff ) for photon interaction and photon energy absorption, and gamma-ray kerma coefficient ( k γ ). The calculations of Z eff are presented as a single-valued and are energy-dependent. While Z eff values were calculated via simplistic power-law method, the energy-dependent Z eff for photon interaction ( Z PI-eff ) and photon energy absorption ( Z PEA-eff ) are obtained via the direct method for energy ranges of 1 keV–100 GeV and 1 keV–20 MeV, respectively. The k γ coefficients are calculated by summing the contributions of the major partial photon interactions for energy range of 1 keV–100 MeV. In most cases, data are presented relative to pure polyethylene to allow direct comparison over a range of energy. The results show that combination of polyethylene with other elements such as lithium and aluminum leads to neutron shielding material with more ability to absorb neutron and γ-rays. Also, the kerma coefficient first increases with Z of the additive element at low photon energies, and then converges with pure polyethylene at energies greater than 100 keV. Neutron shielding materials Effective atomic number Kerma coefficient γ-rays. 1 Introduction 1 Energetic radiations are of serious concern in nuclear powers and accelerator facilities, medical or industrial X-ray machines, and radioisotope production projects. So the radiation shielding is still an attractive topic for research, it aims to preserving both human safety and structural material which may be compromised from radiation exposure. Pure polyethylene (PE), with its high hydrogen concentration, is often used as a moderator to slow fast neutrons to thermal region. The combination of PE and other materials such as boron, lithium or silicon, etc ., makes it an effective neutron shielding material for different purposes. For example borated PE is a useful for neutron shielding in areas of low and intermediate neutron fluxes. Adding lithium to PE can be a shielding for both γ-rays and neutrons. Mixing the PE and silicon leads to neutron shielding with high resistance to fire. Finally, γ-ray shielding parameters are important for the neutron shielding materials that are used in the presence of γ-rays or may be exposed to γ-rays from the (n, γ) reactions. Effective atomic number ( Z eff ) and effective electron density ( N eff ) are useful parameters for describing radiation interaction with composite matter, in terms of equivalent elements. Many authors reported on γ-ray attenuation parameters for ‎concretes [ 1 1 - 3 3 1-3 ], glasses [ 4 4 - 7 7 4-7 ], lunar soil samples [ 8 8 ], building materials [ 9 9 ], alloys [ 10 10 ], low Z -materials [ 11 11 , 12 12 ] and other materials [ 13 13 - 14 14 13-14 ]. However, plenty of them are usually restricted to a candidate γ-ray shielding material or to low- Z materials, which are more important in medical physics. Moreover, the γ-ray shielding properties seem to be limited to photon attenuation properties only, such as linear ‎attenuation coefficient ( μ ), mass attenuation coefficients ( μ / ρ ), and effective atomic numbers. In any case, kerma, K , subject had scarcely treated [ 15 15 ]. The kerma , an acronym for “kinetic energy released in materials”, has replaced the traditional exposure as the shielding design parameter [ 16 16 ]. In fact, gamma heating is the local energy deposition from γ-ray interactions. Often this can be done by estimating the kerma. To relate the radiation passing through a unit volume of a material of interest (fluence, Φ ) to the energy release, K , in the material; the phrase “fluence-to-kerma factors” has been introduced by International Commission on Radiological Units and Measurements, ICRU [ 17 17 ] and widely used [ 18 18 - 20 20 18-20 ]. However, instead of using kerma factor, the kerma coefficient ( k γ , kerma per fluence) will be used through this paper, as the word coefficient implies a physical dimension whereas the word factor does not [ 21 21 ]. Kerma coefficient is of interest for biomedical applications, because it is used to convert photon fluence to kerma (absorbed dose). Also kerma coefficient, which is the key response function for nuclear heating, is important in many nuclear applications, particularly for fission and fusion power reactors [ 18 18 , 20 20 , 22 22 ]. In this work, we aimed at calculating photon shielding energy dependent parameters for seven polyethylene-based neutron shielding materials. The parameters include Z eff and N eff for both photon interaction and photon energy absorption for the energy regions of 1 keV–100 GeV and 1 keV–20 MeV, respectively. The k γ was calculated for the energy region of 1 keV–100 MeV. For many cases, the results are presented relative to pure PE to allow direct comparison over a range of energy. Finally, we are also discussing the single-valued Z eff and N eff generated by the simple power law method, whereby the elemental constituents of a material are summed (weighted according to their atomic percentage or to their relative electron fraction) and raised to a power. 2. Computational method and theoretical basis 1 2.1. The single-valued effective atomic number 2 A simple way to evaluate Z eff is the use of a simple power law: <math display="block" id="M1"><mrow><msub><mi>Z</mi><mrow><mtext>eff</mtext></mrow></msub><mo>=</mo><msup><mrow><mrow><mo>(</mo><mrow><mstyle displaystyle="true"><munderover><mo>∑</mo><mi>i</mi><mrow/></munderover><mrow><mi>a</mi><msup><mi>i</mi><mtext>e</mtext></msup><msubsup><mi>Z</mi><mi>i</mi><mi>m</mi></msubsup></mrow></mstyle></mrow><mo>)</mo></mrow></mrow><mrow><mn>1</mn><mo>/</mo><mi>m</mi></mrow></msup><mo>,</mo></mrow></math> where Z i is atomic number of the i th element and a i e its fractional electronic content. In fact, the researchers consider different values for the exponent m , such as m was 2.94, 3.1, 3.4 and 3.5 by Mayneord [ 23 23 ], Hine [ 24 24 ], Tsaï and Cho [ 25 25 ], and Sellakumar et al . [ 26 26 ]. In addition, the corresponding electron density N eff is given by <math display="block" id="M2"><mrow><msub><mi>N</mi><mrow><mtext>eff</mtext></mrow></msub><mo>=</mo><mrow><mo>(</mo><mrow><mstyle displaystyle="true"><munderover><mo>∑</mo><mi>i</mi><mrow/></munderover><mrow><mi>a</mi><msup><mi>i</mi><mtext>e</mtext></msup><mo stretchy="false">(</mo><msub><mi>Z</mi><mi>i</mi></msub><mo>/</mo><msub><mi>A</mi><mi>i</mi></msub><mo stretchy="false">)</mo></mrow></mstyle></mrow><mo>)</mo></mrow><mo>,</mo></mrow></math> where A i is the atomic mass and N A is the Avogadro’s constant. Also, the Z eff can be given by atomic percentage of each element, α i at , which is defined by the mass percentage w i and A i : <math display="block" id="M3"><mrow><msubsup><mi>α</mi><mi>i</mi><mrow><mtext>at</mtext></mrow></msubsup><mo>=</mo><mfrac><mrow><mo stretchy="false">(</mo><msub><mi>w</mi><mi>i</mi></msub><mo>/</mo><msub><mi>A</mi><mi>i</mi></msub><mo stretchy="false">)</mo></mrow><mrow><mstyle displaystyle="true"><mo>∑</mo><mrow><msub><mrow/><mi>i</mi></msub><mo stretchy="false">(</mo><msub><mi>w</mi><mi>i</mi></msub><mo>/</mo><msub><mi>A</mi><mi>i</mi></msub><mo stretchy="false">)</mo></mrow></mstyle></mrow></mfrac><mo>.</mo></mrow></math> Literature review showed two more definitions: Eq.(4) by Puumalainen et al . [ 27 27 ] and Eq. (5) by Manninen et al . [ 28 28 ]. <math display="block" id="M4"><mrow><msub><mi>Z</mi><mrow><mtext>eff</mtext></mrow></msub><mo>=</mo><msup><mrow><mrow><mo>(</mo><mrow><mfrac><mrow><mstyle displaystyle="true"><msub><mo>∑</mo><mi>i</mi></msub><mrow><msubsup><mi>a</mi><mi>i</mi><mrow><mtext>at</mtext></mrow></msubsup><msubsup><mi>Z</mi><mi>i</mi><mn>3</mn></msubsup></mrow></mstyle></mrow><mrow><mstyle displaystyle="true"><msub><mo>∑</mo><mi>i</mi></msub><mrow><msubsup><mi>a</mi><mi>i</mi><mrow><mtext>at</mtext></mrow></msubsup><msub><mi>Z</mi><mi>i</mi></msub></mrow></mstyle></mrow></mfrac></mrow><mo>)</mo></mrow></mrow><mrow><mn>1</mn><mo>/</mo><mn>2</mn></mrow></msup><mo>,</mo></mrow></math> <math display="block" id="M5"><mrow><msub><mi>Z</mi><mrow><mtext>eff</mtext></mrow></msub><mo>=</mo><mstyle displaystyle="true"><munder><mo>∑</mo><mi>i</mi></munder><mrow><msubsup><mi>a</mi><mi>i</mi><mrow><mtext>at</mtext></mrow></msubsup><msub><mi>Z</mi><mi>i</mi></msub><mo>.</mo></mrow></mstyle></mrow></math> It can easily be shown that a i e and α i at are very nearly equal. Finally, Murty suggested a more simple expression for calculating the Z eff [ 29 29 ], <math display="block" id="M6"><mrow><msub><mi>Z</mi><mrow><mtext>eff</mtext></mrow></msub><msup><mrow><mrow><mo>(</mo><mrow><mstyle displaystyle="true"><msub><mo>∑</mo><mi>i</mi></msub><mrow><msub><mi>w</mi><mi>i</mi></msub><msubsup><mi>Z</mi><mi>i</mi><mrow><mn>3.1</mn></mrow></msubsup></mrow></mstyle></mrow><mo>)</mo></mrow></mrow><mrow><mn>1</mn><mo>/</mo><mn>3.1</mn></mrow></msup></mrow></math> On the other hand, the average electron density is given by <math display="block" id="M7"><mrow><mo><</mo><mi>N</mi><mo>></mo><mo>=</mo><msub><mi>N</mi><mtext>A</mtext></msub><mfrac><mrow><mo><</mo><mi>A</mi><mo>></mo></mrow><mrow><mo><</mo><mi>Z</mi><mo>></mo></mrow></mfrac><mo>,</mo></mrow></math> where < A > and < Z > is the mean atomic mass and the mean atomic number, respectively: <math display="block" id="M8"><mrow><mo><</mo><mi>A</mi><mo>></mo><mo>=</mo><mstyle displaystyle="true"><msub><mo>∑</mo><mi>i</mi></msub><mrow><msub><mi>f</mi><mi>i</mi></msub><msub><mi>A</mi><mi>i</mi></msub></mrow></mstyle><mo>,</mo></mrow></math> <math display="block" id="M9"><mrow><mo><</mo><mi>Z</mi><mo>></mo><mo>=</mo><mstyle displaystyle="true"><msub><mo>∑</mo><mi>i</mi></msub><mrow><msub><mi>f</mi><mi>i</mi></msub><msub><mi>Z</mi><mi>i</mi></msub></mrow></mstyle><mo>,</mo></mrow></math> where A is the atomic mass, Z is the atomic number, and f i is the molar fraction of the i th element (normalized so that ∑ f i =1). This simple approach, however, is approximately valid at low energies where photoelectric absorption is dominating and is overly simplistic for many applications. 2.3. Calculation of energy-dependent Z eff and N eff 2 Energy-dependent Z eff has been described elsewhere [ 30 30 ] in a more rigorous fashion, here only the final equations are given. <math display="block" id="M10"><mrow><msub><mi>Z</mi><mrow><mtext>PI-eff</mtext></mrow></msub><mo>=</mo><mfrac><mrow><msub><mrow><mstyle displaystyle="true"><munder><mo>∑</mo><mi>i</mi></munder><mrow><msub><mi>f</mi><mi>i</mi></msub><msub><mi>A</mi><mi>i</mi></msub><mo stretchy="false">(</mo><mfrac><mi>μ</mi><mi>ρ</mi></mfrac><mo stretchy="false">)</mo></mrow></mstyle></mrow><mi>i</mi></msub></mrow><mrow><msub><mrow><mstyle displaystyle="true"><munder><mo>∑</mo><mi>i</mi></munder><mrow><msub><mi>f</mi><mi>i</mi></msub><mfrac><mrow><msub><mi>A</mi><mi>i</mi></msub></mrow><mrow><msub><mi>Z</mi><mi>i</mi></msub></mrow></mfrac><mo stretchy="false">(</mo><mfrac><mi>μ</mi><mi>ρ</mi></mfrac><mo stretchy="false">)</mo></mrow></mstyle></mrow><mi>i</mi></msub></mrow></mfrac><mo>,</mo></mrow></math> where ( μ en / ρ ) is the mass attenuation coefficient. NXcom program has been used to calculate the mass attenuation coefficients for the energy range of 1 keV–100 GeV [ 31 31 ]. The effective electron density, N PI-eff , is given by (in electron per gram): <math display="block" id="M11"><mrow><msub><mi>N</mi><mrow><mtext>PI-eff</mtext></mrow></msub><mo>=</mo><msub><mi>N</mi><mtext>A</mtext></msub><mfrac><mrow><msub><mi>Z</mi><mrow><mtext>PI-eff</mtext></mrow></msub></mrow><mrow><mo><</mo><mi>A</mi><mo>></mo></mrow></mfrac><mo>.</mo></mrow></math> The effective atomic number ( Z PEA-eff ) and the effective electron density ( N PEA-eff ) for photon energy absorption can be obtained from Eqs. (10) and (11), respectively. Mass energy-absorption coefficients are calculated for the studied object at photon energies using the photon data from Hubbell and Seltzer [ 32 32 ]. 2.2. Computation of photon kerma coefficients. 2 The kerma , K, is the quotient of d E tr by d m , where d E tr is the mean sum of the initial energies of all charged particles liberated by indirectly ionizing radiations such as neutrons and photons in a mass d m of a material [ 33 33 ], thus <math display="block" id="M12"><mrow><mi>K</mi><mo>=</mo><mtext>d</mtext><msub><mi>E</mi><mrow><mi>t</mi><mi>r</mi></mrow></msub><mo>/</mo><mtext>d</mtext><mi>m</mi></mrow></math> For a fluence, Φ , of uncharged particles of energy E , K in a specified material is given by <math display="block" id="M13"><mrow><mi>K</mi><mo>=</mo><mi>Φ</mi><mi>E</mi><mo stretchy="false">(</mo><msub><mi>μ</mi><mrow><mi>t</mi><mi>r</mi></mrow></msub><mo>/</mo><mi>ρ</mi><mo stretchy="false">)</mo><mtext> </mtext><mo>=</mo><mi>Ψ</mi><mo stretchy="false">(</mo><msub><mi>μ</mi><mrow><mi>t</mi><mi>r</mi></mrow></msub><mo>/</mo><mi>ρ</mi><mo stretchy="false">)</mo></mrow></math> where ( μ tr / ρ ) is the mass energy-transfer coefficient of the material for these particles, and Ψ is the energy fluence. The kerma per fluence, K / Φ , is termed the kerma coefficient, k , for uncharged particles of energy E in a specified material, thus <math display="block" id="M14"><mrow><mi>k</mi><mo>=</mo><mi>K</mi><mo>/</mo><mi>Φ</mi><mo>=</mo><mi>E</mi><mo stretchy="false">(</mo><msub><mi>μ</mi><mrow><mi>t</mi><mi>r</mi></mrow></msub><mo>/</mo><mi>ρ</mi><mo stretchy="false">)</mo><mo>.</mo></mrow></math> For low Z -materials (e.g., air, water and soft tissue) only a small part is expended in bremsstrahlung (collisions with atomic nuclei). Therefore, mass energy transfer coefficient ( μ tr / ρ ) and mass energy mass energy-absorption coefficient ( μ en / ρ ) are almost equal. Thereby, photon kerma coefficients (in Gy·cm 2 /photon) of such materials can be obtained by Eq.(15): <math display="block" id="M15"><mrow><msub><mi>k</mi><mi>γ</mi></msub><mo>=</mo><msub><mi>k</mi><mtext>D</mtext></msub><msub><mi>E</mi><mi>γ</mi></msub><mstyle displaystyle="true"><msub><mo>∑</mo><mi>i</mi></msub><mrow><msub><mi>w</mi><mi>i</mi></msub><msub><mrow><mo stretchy="false">[</mo><msub><mi>μ</mi><mrow><mtext>en</mtext></mrow></msub><mo stretchy="false">(</mo><msub><mi>E</mi><mi>γ</mi></msub><mo stretchy="false">)</mo><mo>/</mo><mi>ρ</mi><mo stretchy="false">]</mo></mrow><mi>i</mi></msub></mrow></mstyle></mrow></math> where k D = 1.602×10 −10 Gy·g/MeV is the energy conversion coefficient. On the other hand, when the bremsstrahlung production is not negligible, the major contributing reactions for gamma kerma coefficients (in Gy·cm 2 /photon) in the energy range of fusion systems are the photoelectric (Pe), Compton scattering ( C ), pair production reactions (pp), and coherent (no energy loss) [ 18 18 ]: <math display="block" id="M16"><mrow><msub><mi>k</mi><mi>γ</mi></msub><mo>=</mo><msub><mi>k</mi><mtext>D</mtext></msub><mstyle displaystyle="true"><msub><mo>∑</mo><mi>i</mi></msub><mrow><msub><mi>w</mi><mi>i</mi></msub><mo stretchy="false">[</mo><msubsup><mi>σ</mi><mrow><mtext>Pe</mtext></mrow><mi>i</mi></msubsup><msub><mi>E</mi><mi>γ</mi></msub><mo>+</mo><msubsup><mi>σ</mi><mrow><mtext>pp</mtext></mrow><mi>i</mi></msubsup><mo stretchy="false">(</mo><msub><mi>E</mi><mi>γ</mi></msub><mo>−</mo><mn>1.022</mn><mo stretchy="false">)</mo><mo>+</mo><msubsup><mi>σ</mi><mtext>C</mtext><mi>i</mi></msubsup><mo stretchy="false">(</mo><msub><mi>E</mi><mi>γ</mi></msub><mo>−</mo><msubsup><mi>E</mi><mi>γ</mi><mi>i</mi></msubsup><mo stretchy="false">)</mo><mo stretchy="false">]</mo></mrow></mstyle></mrow></math> where σ pe i , σ pp i and σ C i are photoelectric, pair production and Compton absorption cross-sections for element i , respectively. It is important to recognize that above 5 MeV, one in principle should include photonuclear absorption cross section σ ph.n . However, this process is not readily amenable to systematic calculation and tabulation. Hence, σ ph.n is ignored in current compilations, even though at its giant resonance peak between 5 and 40 MeV it can contribute between 2% (high- Z elements) and 6% (low- Z elements) to the total cross section σ tot [ 34 34 ]. Implicit in the use of Eq. (16) is the assumption that all the gamma photon energy in the photoelectric process is deposited locally and in pair production, 1.022 MeV (energy of electron-positron masses) of the photon energy is not available for local deposition. In the Compton scattering reaction, the photon only deposits a fraction of its energy locally because the scattered photon carries the rest away. 3. Results and Discussion 1 The chemical composition of neutron shielding materials studied in the present work is given in Table 1 Table 1 [ 35 35 ], for combinations of PE and other materials such as boron, lithium or silicon. The single-valued calculated energy dependent Z eff and N eff for all samples are listed in Table 2 Table 2 . 1 1 Sample No. 1 1 H 1 1 Li 1 1 B 1 1 C 1 1 O 1 1 Na 1 1 Mg 1 1 Al 1 1 Si 1 1 S 1 1 Cl 1 1 Ca 1 1 Ti 1 1 Mn 1 1 Fe 1 1 Zn 1 1 Sr 1 1 238 1 1 2.76 1 1 ― 1 1 25.35 1 1 20.08 1 1 24.15 1 1 ― 1 1 ― 1 1 ― 1 1 26.93 1 1 ― 1 1 ― 1 1 0.04 1 1 0.02 1 1 ― 1 1 0.41 1 1 0.26 1 1 ― 1 1 207HD 1 1 6.01 1 1 ― 1 1 0.86 1 1 18.02 1 1 48.42 1 1 0.17 1 1 0.02 1 1 24.9 1 1 0.04 1 1 0.01 1 1 ― 1 1 1.36 1 1 ― 1 1 ― 1 1 0.01 1 1 ― 1 1 0.06 1 1 210 1 1 8.54 1 1 ― 1 1 30 1 1 59.25 1 1 0.76 1 1 ― 1 1 ― 1 1 0.04 1 1 1.01 1 1 ― 1 1 ― 1 1 ― 1 1 ― 1 1 0.02 1 1 0.38 1 1 ― 1 1 ― 1 1 2015 1 1 8.6 1 1 7.5 1 1 0 1 1 57.76 1 1 26.13 1 1 0.001 1 1 ― 1 1 ― 1 1 ― 1 1 ― 1 1 0.004 1 1 0.001 1 1 ― 1 1 ― 1 1 ― 1 1 ― 1 1 ― 1 1 261 1 1 8.94 1 1 ― 1 1 10 1 1 36.65 1 1 44.41 1 1 ― 1 1 ― 1 1 ― 1 1 ― 1 1 ― 1 1 ― 1 1 ― 1 1 ― 1 1 ― 1 1 ― 1 1 ― 1 1 ― 1 1 201 1 1 11.6 1 1 ― 1 1 5 1 1 61.2 1 1 22.2 1 1 ― 1 1 ― 1 1 ― 1 1 ― 1 1 ― 1 1 ― 1 1 ― 1 1 ― 1 1 ― 1 1 ― 1 1 ― 1 1 ― 1 1 213 1 1 14.4 1 1 ― 1 1 ― 1 1 85.6 1 1 ― 1 1 ― 1 1 ― 1 1 ― 1 1 ― 1 1 ― 1 1 ― 1 1 ― 1 1 ― 1 1 ― 1 1 ― 1 1 ― 1 1 ― *Ref. [ 26 26 ] 1 1 Sample No. 1 1 < A > 1 1 < N > 6 1 Z eff 1 1 N eff * (×10 23 e·g −1 ) 1 1 < N > (×10 23 e·g −1 ) 1 1   1 1   1 1   1 1 Ref.[ 29 29 ] 1 1 Ref.[ 20 20 ] 1 1 Ref.[ 24 24 ] 1 1 Ref.[ 25 25 ] 1 1 Ref.[ 26 26 ] 1 1 Ref.[ 28 28 ] 1 1   1 1   1 1 238 1 1 10.83 1 1 5.45 1 1 10.37 1 1 10.15 1 1 10.30 1 1 10.52 1 1 10.67 1 1 9.25 1 1 3.12 1 1 3.03 1 1 207HD 1 1 8.66 1 1 4.54 1 1 10.01 1 1 9.66 1 1 9.77 1 1 9.93 1 1 10.06 1 1 9.00 1 1 3.32 1 1 3.16 1 1 213 1 1 7.23 1 1 3.87 1 1 8.82 1 1 8.70 1 1 8.62 1 1 6.95 1 1 6.96 1 1 6.66 1 1 3.41 1 1 3.22 1 1 261 1 1 6.40 1 1 3.46 1 1 6.85 1 1 6.64 1 1 6.67 1 1 6.70 1 1 6.75 1 1 6.39 1 1 3.48 1 1 3.25 1 1 215 1 1 6.23 1 1 3.34 1 1 6.41 1 1 6.22 1 1 6.26 1 1 6.29 1 1 6.34 1 1 5.97 1 1 3.46 1 1 3.23 1 1 210 1 1 6.15 1 1 3.26 1 1 6.53 1 1 6.21 1 1 6.34 1 1 6.59 1 1 6.73 1 1 5.62 1 1 3.42 1 1 3.19 1 1 201 1 1 5.42 1 1 3.01 1 1 6.33 1 1 6.09 1 1 6.12 1 1 6.15 1 1 6.20 1 1 5.82 1 1 3.61 1 1 3.34 Examples of such calculations are graphed in Fig. 1 Fig. 1 . The single-valued effective atomic numbers are only a rough approximation at low energies, and should be treated with some caution. At medium photon energies (about 0.5–3.5 MeV) where Compton scattering process is the most probable, one can notice that the Z eff and N eff are approximately equal to their mean values < Z > and < N > , respectively. Therefore the mean values can represent Z eff and N eff at medium photon energies only. 3.1. The effective atomic number 2 The calculated values of energy dependent Z PI-eff , N PI-eff , Z PEA-eff and N PEA-eff for PE and six neutron shielding materials, and the ratios of ( Z PI-eff ) PE m , ( N PI-eff ) PE m , ( Z PEA-eff ) PE m and ( N PEA-eff ) PE m , are shown in Fig.2 Fig.2 , as a function of photon energy. The energy dependence of Z PI-eff (for total photon interaction) for all the neutron shielding materials have almost the same behavior. One can see that the curves of Z eff are arranged in descending order according to their mean atomic number < Z > . Fig.2(a) Fig.2(a) ,shows three energy regions, where Z eff is almost constant “plateau” for a given material, while the transition regions are characterized by a rapid variation of Z eff . This energy behavior of Z eff reflects relative importance of the three γ-ray processes: photoelectric, Compton and pair production. Since the plateau region is formed when only one photon interaction is dominant, whereas the transition one includes more than one interaction which making significant contributions to the total interaction. Therefore, we have three plateaus corresponding to the dominance regions of the three γ-ray interactions, photoelectric absorption ( E < 0.01 MeV), Compton scattering (0.1 < E < 6 MeV), and pair production ( E > 200 MeV). In addition two transition regions (0.01 < E < 0.1 MeV), and (6 < E < 200 MeV) are observed. At a given photon energy, the interaction is proportional to Z n where n is about 4 for the photoelectric absorption, 1 for the Compton scattering, and 2 for pair production. The Z 4 dependence of the photoelectric absorption cross section leads to heavy weight to the element with the highest atomic number and the maximum value of Z eff is found in the region of photoelectric plateau. At intermediate energies, Compton scattering is the main interaction process, Z eff is close to the mean atomic number of the material, as the Compton scattering cross section of a an element is proportional to Z . The minimum value of Z eff is recorded in this intermediate plateau. Above 200 MeV, the pair production plateau is observed, where Z eff is almost constant and its mean value is smaller than that obtained for photoelectric plateau. This is due to the fact that the pair production cross section is proportional to Z 2 , giving less weight to the higher- Z elements than the photoelectric absorption cross section. The same arguments also hold for Z PEA-eff as shown in Fig. 2(c) Fig. 2(c) . To compare the Z eff values of a given PE-based shielding material with that of polyethylene, Figs.2(e–h) Figs.2(e–h) present the ratios of Z eff of the materials for photon interaction and photon energy absorption to those of PE as a function of photon energy. Figs. 2(e) and 2(g) Figs. 2(e) and 2(g) show that the materials containing higher- Z elements, such as Samples 238 and 207HD, have Z eff values systematically greater than that of PE over the full energy range studied. Such neutron shielding materials have a competitive advantage in attenuating gamma-rays. 3.2. The effective electron density 2 Figs. 2(b) and 2(d) Figs. 2(b) and 2(d) show the energy dependence of N eff for total photon interaction and photon energy absorption, respectively. Figs. 2(f) and 2(h) Figs. 2(f) and 2(h) present the N eff ratios with respect to PE. The behavior is similar to that of Z eff , because the two parameters are closely related through the Eq. (7). However, in contrast to Z eff , the obtained curves for N eff are not arranged by their average electron density < N > ( Table 1 Table 1 ). Thus, we can conclude that the average electron density of the material does not give us any correct indication of real effective electron density value. Finally, it is obvious that the different single-valued effective atomic numbers, which are estimated using simple power law, are approximately valid at low energies where photoelectric absorption is dominating as shown in Fig. 1 Fig. 1 . 3.3. Kerma coefficients 2 The energy dependence of γ-ray kerma coefficient, which depends on the photon interaction cross sections, is given in Fig. 3 Fig. 3 . Surprisingly, there is no plateau, and only rapid variation of k γ with energy is seen. The kerma coefficient curves of the studied materials have dips in the energy range of around 0.3 to 0.09 MeV, and the dip position shifts to higher energy with increasing concentration of high Z-elements. We have also seen that the K-, and L-absorption edges in graphs of 210, 207HD, and 238 samples due to the presence of high-Z elements such Zn and Sr ( Table 3 Table 3 ). Fig.3 Fig.3 also shows that the k γ curves tend to converge after departure of the photoelectric absorption dominance region ( E > 0.1 MeV), where the other photon interaction processes give less weight to the higher-Z additive elements than that given by photoelectric process. 1 1 Elements 1 1 Z 1 1 L3 1 1 L2 1 1 L1 1 1 K 1 1 Na 1 1 11 1 1 ― 1 1 ― 1 1 ― 1 1 1.07 1 1 Mg 1 1 12 1 1 ― 1 1 ― 1 1 ― 1 1 1.31 1 1 Al 1 1 13 1 1 ― 1 1 ― 1 1 ― 1 1 1.56 1 1 Si 1 1 14 1 1 ― 1 1 ― 1 1 ― 1 1 1.84 1 1 S 1 1 16 1 1 ― 1 1 ― 1 1 ― 1 1 2.47 1 1 Ca 1 1 20 1 1 ― 1 1 ― 1 1 ― 1 1 4.04 1 1 Ti 1 1 22 1 1 ― 1 1 ― 1 1 ― 1 1 4.97 1 1 Mn 1 1 25 1 1 ― 1 1 ― 1 1 ― 1 1 6.54 1 1 Fe 1 1 26 1 1 ― 1 1 ― 1 1 ― 1 1 7.11 1 1 Zn 1 1 30 1 1 1.02 1 1 1.04 1 1 1.19 1 1 9.66 1 1 Sr 1 1 38 1 1 1.94 1 1 2.01 1 1 2.22 1 1 16.10 4. Conclusions 1 1. Combination of polyethylene with other elements such as lithium and aluminum leads to neutron shielding material with more ability to absorb neutron capture γ–rays, therfore these materials are good sheilding for neutron and gamma-ray. 2. Variations of Z eff , and N eff with photon energy showed three basic energy plataeus and two transition regions which are characterized by a rapid variation of Z eff and N eff . Both plateous and transition regions were attributed to dominance of one and more photon interaction processes, respectively. 3. The three energy plateous are approximately E < 0.01 MeV, 0.1 MeV < E < 6 MeV and E > 200 MeV. 4. The single-valued effective atomic numbers are calculated by different 6-expressions and the results are approximately valid at low energies Therefore should be treated with some caution. 5. The energy dependence of γ-ray kerma coefficients showed rapid variations with energy and they have dips in the energy range of around 0.3 to 0.9 MeV. 6. All curves of γ-ray kerma coefficients converge at energies greater than 0.1 MeV. References Akkurt , A. M. El-Khayatt , Effective atomic number and electron density of marble concrete . J. Radioanal. Nucl. Chem . 259 ( 1 ), 633 - 638 ( 2013 ).   doi: 10.1016/j.anucene.2013.04.021 http://doi.org/10.1016/j.anucene.2013.04.021 Effective atomic number and electron density of marble concrete E. Yılmaz , H. Baltas , E. Kırısa , İ. Ustabas , U Cevik , A. M. El-Khayatt , Gamma-ray and neutron shielding properties of some concrete materials . Ann. Nucl. Energy 38 ( 10 ), 128 - 132 ( 2011 ). doi: 10.1016/j.anucene.2011.06.011 http://doi.org/10.1016/j.anucene.2011.06.011 Gamma-ray and neutron shielding properties of some concrete materials M. Kurudirek . Radiation shielding and effective atomic number studies in different types of shielding concretes, lead base and non-lead base glass systems for total electron interaction: A comparative study . Nucl Eng Des 280 , 440 ( 2014 ). doi: 10.1016/j.nucengdes.2014.09.020 http://doi.org/10.1016/j.nucengdes.2014.09.020 Radiation shielding and effective atomic number studies in different types of shielding concretes, lead base and non-lead base glass systems for total electron interaction: A comparative study N. Chanthima , J. Kaewkhao , P. Limsuwan . Study of photon interactions and shielding properties of silicate glasses containing Bi2O3, BaO and PbO in the energy region of 1 keV to 100 GeV . Ann Nucl Energy 41 , 119 ( 2012 ). doi: 10.1016/j.anucene.2011.10.021 http://doi.org/10.1016/j.anucene.2011.10.021 Study of photon interactions and shielding properties of silicate glasses containing Bi2O3, BaO and PbO in the energy region of 1 keV to 100 GeV S. Ruengsri , S. Insiripong , N. Sangwaranatee , J. Kaewkhao . Development of barium borosilicate glasses for radiation shielding materials using rice husk ash as a silica source . Prog Nucl Energ 83 , 99 ( 2015 ). doi: 10.1016/j.pnucene.2015.03.006 http://doi.org/10.1016/j.pnucene.2015.03.006 Development of barium borosilicate glasses for radiation shielding materials using rice husk ash as a silica source A. M. El-Khayatt , A. M. Ali , V. P. Singh , Photon attenuation coefficients of Heavy-Metal Oxide glasses by MCNP, XCOM program and experimental data: A comparison study . Nucl. Instrum. Methods Phys. Res. A   735 ,  207 - 212 ( 2014 ). doi: 10.1016/j.nima.2013.09.027 http://doi.org/10.1016/j.nima.2013.09.027 Photon attenuation coefficients of Heavy-Metal Oxide glasses by MCNP, XCOM program and experimental data: A comparison study N. Chanthima , J. Kaewkhao . Investigation on radiation shielding parameters of bismuth borosilicate glass from 1 keV to 100 GeV . Ann Nucl Energy 55 , 23 ( 2013 ). doi: 10.1016/j.anucene.2012.12.011 http://doi.org/10.1016/j.anucene.2012.12.011 Investigation on radiation shielding parameters of bismuth borosilicate glass from 1 keV to 100 GeV A. M. El-Khayatt , M. A. Al-Rajhi , Analysis of some lunar soil and rocks samples in terms of photon interaction and photon energy absorption . Adv. Space Res . 55 ( 7 ), 1816 - 1822 ( 2015 ) doi: 10.1016/j.asr.2015.01.020 http://doi.org/10.1016/j.asr.2015.01.020 . Analysis of some lunar soil and rocks samples in terms of photon interaction and photon energy absorption V. P. Singh , N. M. Badiger , A. M. El-Khayatt , Study on gamma-ray exposure buildup factors and fast neutron shielding properties of some building materials . Radiat. Eff. Defects Solids .  169 ( 6 ),  547 - 559 ( 2014 ). doi: 10.1080/10420150.2014.905942 http://doi.org/10.1080/10420150.2014.905942 . Study on gamma-ray exposure buildup factors and fast neutron shielding properties of some building materials P. Limkitjaroenporn , J. Kaewkhao , S. Asavavisithchai . Determination of mass attenuation coefficients and effective atomic numbers for Inconel 738 alloy for different energies obtained from Compton scattering . Ann Nucl Energy 53 , 64 ( 2013 ). doi: 10.1016/j.anucene.2012.08.020 http://doi.org/10.1016/j.anucene.2012.08.020 Determination of mass attenuation coefficients and effective atomic numbers for Inconel 738 alloy for different energies obtained from Compton scattering A. M. El-Khayatt , A. M. Ali , V. P. Singh , Badiger N. M. , Determination of mass attenuation coefficient of low-Z dosimetric materials .   Radiat. Eff. Defects Solids , 169 ( 12 ), 1038 - 1044 ( 2014 ). doi: 10.1080/10420150.2014.988626 http://doi.org/10.1080/10420150.2014.988626 Determination of mass attenuation coefficient of low-Z dosimetric materials M. Kurudirek . Effective atomic numbers and electron densities of some human tissues and dosimetric materials for mean energies of various radiation sources relevant to radiotherapy and medical applications . Radiat. Phys. Chem . 1021 , 39 ( 2014 ). doi: 10.1016/j.radphyschem.2014.04.033 http://doi.org/10.1016/j.radphyschem.2014.04.033 Effective atomic numbers and electron densities of some human tissues and dosimetric materials for mean energies of various radiation sources relevant to radiotherapy and medical applications M. Kurudirek . Effective atomic number, energy loss and radiation damage studies in some materials commonly used in nuclear applications for heavy charged particles such as H, C, Mg, Fe, Te, Pb and U . Radiat. Phys. Chem . 122 , 15 ( 2016 ). doi: 10.1016/j.radphyschem.2016.01.012 http://doi.org/10.1016/j.radphyschem.2016.01.012 Effective atomic number, energy loss and radiation damage studies in some materials commonly used in nuclear applications for heavy charged particles such as H, C, Mg, Fe, Te, Pb and U P. Limkitjaroenporn , J. Kaewkhao . Gamma-rays attenuation of zircons from Cambodia and South Africa at different energies: A new technique for identifying the origin of gemstone . Radiation Radiat. Phys. Chem . 103 , 67 ( 2014 ). doi: 10.1016/j.radphyschem.2014.05.035 http://doi.org/10.1016/j.radphyschem.2014.05.035 Gamma-rays attenuation of zircons from Cambodia and South Africa at different energies: A new technique for identifying the origin of gemstone El-Khayatt A. M. , Vega-Carrillo H. R . , Photon and neutron kerma coefficients for polymer gel dosimeters . Nucl. Instrum. Methods Phys. Res. A 792 , 6 - 10 ( 2015 ). doi: 10.1016/j.nima.2015.04.033 http://doi.org/10.1016/j.nima.2015.04.033 Photon and neutron kerma coefficients for polymer gel dosimeters J.E. Turner , Atoms, Radiation, and Radiation Protection , Third edition, ( WILEY- VCH Verlag GmbH & Co. KGaA , Weinheim , 2007 ). Atoms, Radiation, and Radiation Protection ICRU Report 84 (1962). Radiation quantities and units . Bethesda, MD , International Commission on Radiation Units and Measurements . ICRU Report 84 (1962). Radiation quantities and units L. Zhang , M. A. Abdou , Kerma Factor Evaluation and its Application in Nuclear Heating Experiment Analysis . Fusion Eng Des , 36 , 479 - 503 ( 1997 ). doi: 10.1016/S0920-3796(96)00704-1 http://doi.org/10.1016/S0920-3796(96)00704-1 Kerma Factor Evaluation and its Application in Nuclear Heating Experiment Analysis R.S. Caswell , J. J. Coyne , M.L. Randolph , Int J Appl Radiat Isot . 33 ( 11 ), 1227 - 1262 ( 1982 ). doi: 10.1016/0020-708X(82)90246-0 http://doi.org/10.1016/0020-708X(82)90246-0 M. A. Abdou , C. W. Maynard , Calculational Methods for Nuclear Heating-Part I: Theoretical and Computational Algorithms . Nucl Sci Eng 56 ( 4 ), 360 - 380 ( 1975 ). Calculational Methods for Nuclear Heating-Part I: Theoretical and Computational Algorithms ICRU Report 85, Fundamental quantities and units for ionizing radiation (Revised) , Journal of the ICRU 11 ( 1 ), 22 ( 2011 ). doi: 10.1093/jicru/ndr009 http://doi.org/10.1093/jicru/ndr009 . ICRU Report 85, Fundamental quantities and units for ionizing radiation (Revised) H. C. Claiborne , M. Solomito , J. J. Ritts , Heat generation by neutrons in some moderating and shielding materials . Nucl Eng Des 15 , 232 - 236 ( 1971 ). doi: 10.1016/0029-5493(71)90065-3 http://doi.org/10.1016/0029-5493(71)90065-3 Heat generation by neutrons in some moderating and shielding materials M. Faiz Khan , The Physics of Radiation Therapy . Williams and Wilkins Publishers , Maryland, USA 1994 . The Physics of Radiation Therapy G. E. Hine , Secondary electron emission and effective atomic numbers . Nucleonics 10 ( 1 ), 9 - 15 ( 1952 ) CODEN: NUCLAM; ISSN: 0096-6207 Secondary electron emission and effective atomic numbers C.M. Tsaï , Z.H. Cho , Physics of contrast mechanism and averaging effect of linear attenuation coefficients in a computerized transverse axial tomography (CTAT) transmission scanner . Physic Med. Biol ., 21 ( 4 ), 544 - 59 ( 1976 ). doi: 10.1088/0031-9155/21/4/006 http://doi.org/10.1088/0031-9155/21/4/006 Physics of contrast mechanism and averaging effect of linear attenuation coefficients in a computerized transverse axial tomography (CTAT) transmission scanner P. Sellakumar , James Jebaseelan Samuel E ., Sanjay S. Supe , Water equivalence of polymer gel dosimeters . Radiat. Phys. Chem . 76 ( 7 ), 1108 - 1115 ( 2007 ). doi: 10.1016/j.radphyschem.2007.03.003 http://doi.org/10.1016/j.radphyschem.2007.03.003 Water equivalence of polymer gel dosimeters P. Puumalainen , H. Olkkonen , P. Sikanen , Assessment of fat content of liver by a photon scattering technique . Int. J. Appl. Radiat. Isot . 28 ( 9 ), 785 - 7 ( 1977 ). doi: 10.1016/0020-708X(77)90110-7 http://doi.org/10.1016/0020-708X(77)90110-7 Assessment of fat content of liver by a photon scattering technique S. Manninen , T. Pitkanen , S. Koikkalainen , T. Paakkari , Determination of the effective atomic number using elastic and inelastic scattering of γ-rays . Int. J. Appl. Radiat. Isot . 35 ( 10 ), 965 - 968 ( 1984 ). doi: 10.1016/0020-708X(84)90212-6 http://doi.org/10.1016/0020-708X(84)90212-6 Determination of the effective atomic number using elastic and inelastic scattering of γ-rays R.C. Murty , Effective atomic numbers of heterogeneous materials . Nature 207 ,, 398 - 399 ( 1965 ). doi: 10.1038/207398a0 http://doi.org/10.1038/207398a0 Effective atomic numbers of heterogeneous materials S.R. Manohara , S.M. Hanagodimath , K.S. Thind , L. Gerward , The effectıve atomıc number revısıted ın the lıght of modern photon-ınteractıon cross-sectıon databases . Appl. Radiat. Isot . 68 ( 4-5 ), 784 - 7 ( 2010 ). doi: 10.1016/j.apradiso.2009.09.047 http://doi.org/10.1016/j.apradiso.2009.09.047 The effectıve atomıc number revısıted ın the lıght of modern photon-ınteractıon cross-sectıon databases 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 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 NISTIR 5632 ( National Institute of Standards and Technology , Gaithersburg , 1995 ). 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 NISTIR 5632 F. H. Attix . Introduction to Radiological Physics and Radiation Dosimetry . ( New York , Wiley , 1986 ). Introduction to Radiological Physics and Radiation Dosimetry J H Hubbell , Review of photon interaction cross section data in the medical and biological context . Phys. Med. Biol . 44 , R1 ( 1999 ). doi: 10.1088/0031-9155/44/1/001 http://doi.org/10.1088/0031-9155/44/1/001 Review of photon interaction cross section data in the medical and biological context Bladewerx Company , http://www.shieldwerx.com http://www.shieldwerx.com . 10.1007/s41365-017-0222-y 2016-06-15 2016-09-06 2016-09-15 2017-03-31 2017-05-01 2021-04-11 © 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. 2017 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-017-0222-y http://dx.doi.org/10.1007/s41365-017-0222-y . 2017 Nuclear Science and Techniques 05 28