A comparative Monte Carlo simulation study on shielding features of the CaF2–CaO–B2O3–P2O5–SrO–Ta2O5 glass system against X-ray by GEANT4 and MCNPX codes
Farzad Isazadeh, Akbar Abdi Saray

TL;DR
This study compares the X-ray shielding properties of a specific glass system using two simulation codes and finds the best-performing sample.
Contribution
The study evaluates and compares the shielding performance of a novel glass system using GEANT4 and MCNPX simulations.
Findings
GEANT4 and MCNPX codes show good compatibility in calculating MAC values for the glass samples.
The Ta3 sample of the glass system is identified as the best X-ray shield among the tested samples.
At 40 keV, the MAC values from simulations differ by less than 8% compared to Phy-X/PSD results.
Abstract
Regarding to their unique physical and mechanical features, glasses and glass–ceramics are suitable materials for shielding purposes. The present study evaluates the shielding properties of the CaF2–CaO–B2O3–P2O5–SrO–Ta2O5 glass system using Monte Carlo GEANT4 and MCNPX codes for X-ray radiations with an energy range of 20 to 100 keV. MAC values of the Ta0, Ta1, Ta2, Ta2.5, and Ta3 samples of the CaF2–CaO–B2O3–P2O5–SrO–Ta2O5 glass were computed using Phy-X/PSD, GEANT4, and MCNPX codes and compared. According to the results, the programs have good compatibility with each other. For instance, in the energy of 40 keV and for the Ta2 sample, GEANT4 and MCNP codes are 1.445765406 and 1.517801204 cm2/g, respectively, indicating 7.419529525 and 2.829628418% differences with 1.562 cm2/g obtained using the Phy-X/PSD software. According to recent estimations, the Ta3 sample of the…
Genes, proteins, chemicals, diseases, species, mutations and cell lines named across the full text — each resolved to its canonical identifier and authoritative record.
Click any figure to enlarge with its caption.
Figure 10
Figure 11
Figure 12
Figure 13
Figure 14
Figure 15
Figure 1
Figure 2
Figure 3
Figure 4
Figure 5
Figure 6
Figure 7
Figure 8
Figure 9Peer Reviews
No public reviews on file for this paper yet. If you reviewed it on a platform where reviews are public (OpenReview, ICLR, NeurIPS, ICML), you can paste yours below so the community can read it here.
Videos
No videos yet. Explain this paper in a talk, walkthrough, or lecture? Add one.
Taxonomy
TopicsPlant Pathogens and Fungal Diseases · Forest Insect Ecology and Management · Mycorrhizal Fungi and Plant Interactions
Introduction
Nowadays, with the advancement of nuclear technology in various areas of science (e.g., industries, nuclear reactors, and especially medical centers), humans, plants, and animals are exposed to neutrons, heavy particles and ions, and gamma radiation. Therefore, research on protecting humans, plants, and animals from serious hazards of those radiations has been the subject of intense research^1,2^. Over the past decades, various shielding materials, e.g., wide range of stones, concretes, alloys, polymers, and glasses have been developed and used as materials with appropriate mechanical, optical, and shielding properties for protecting from radiation dangers^3–7^. Saleh et al.^8^ studied the structure, mechanical, and nuclear shielding properties of lead-free Sn–Zn–Bi alloys. Adib estimated neutron characteristics of single-crystal magnesium fluoride^9^. Also, Elnaggar^10^ investigated the effect of heavy chemical additives on the mechanical properties of cement pastes, radiation exposure the corrosion behavior of reinforced steel upon exposure to seawater. Recently, glass and glass–ceramic production has been interesting for researchers regarding their unique physical, chemical, and mechanical features (e.g., facile production, high optical transparency, low cost, low melting point, high thermal stability, chemical persistence, and easy formatting in shapes and sizes)^11,12^. The gamma and neutron radiation features of various sorts of glasses have been studied either by applying experimental methods or using simulation codes such as Geant4, MCNPX, and FLUKA^13–15^. Alzahrani et al.^16^ evaluated the radiation shielding properties of boron silicate glasses containing CrO_3_. Moreover, Alzahrani et al.^17^ studied the optical properties and radiation shielding effectiveness of Sm^3+^‑doped Zinc Boro phosphate glasses. Saleh et al.^18^ studied the effect of the PbO on the physical, structural, and protection features of TeO_2_–B_2_O_3_ glass against gamma, neutron, proton and alpha particles. Also, Saleh explored gamma and neutron shielding features of the gadolinium silico-borate glasses system^19^. According to Madhavi et al.^20^, Ta-dopped glasses have good biocompatibility, high mechanical strength, and high coefficient of friction. Also, the Ta-dopped glass system of CaF_2_–CaO–B_2_O_3_–P_2_O_5_–SrO–Ta_2_O_5_ has in-vitro bioactivity and antibacterial properties. Moreover, the shielding properties of the glass have not been studied before. Therefore, consistent with previous studies^21–24^, in this study, the shielding parameters of the CaF_2_–CaO–B_2_O_3_–P_2_O_5_–SrO–Ta_2_O_5_ glass system are evaluated using the GEANT4 toolkit and MCNPX code for the X-ray radiations. In this respect, the bioactivity and antibacterial properties of the glass have already been studied by Madhavi^20^. The calculations are considered for X-ray photons with energy range between 20 and 100 keV. Phy-X/PSD software was selected to evaluate MAC values and to compare them with the simulation results of MCNPX and GEANT4.
Materials and methods
Sample details
The CaF_2_–CaO–B_2_O_3_–P_2_O_5_–SrO–Ta_2_O_5_ glass has already been synthesized by the melt quenching technique^20^. Five glass samples were prepared with stable percentages of CaF_2_, CaO, B_2_O_3_, and SrO and different chemical compositions of P_2_O_5_ and Ta_2_O_5_^20^. The samples named as Ta_0_ (0% Ta_2_O_5_), Ta_1_ (1% Ta_2_O_5_), Ta_2_ (2% Ta_2_O_5_), Ta_2.5_ (2.5% Ta_2_O_5_), and Ta_3_ (3% Ta_2_O_5_)^20^. Table 1 represents the chemical compositions and density of each sample that is simulated using MCNPX and GEANT4 codes in this study, as reported by Madhavi^20^.Table 1. Chemical composition (mole %) and density details of the glass samples.GlassCaF_2_CaOB_2_O_3_P_2_O_5_SrOTa_2_O_5_Density (g/cm^3^)Ta_0_101010601002.57Ta_1_101010591012.752Ta_2_101010581022.83Ta_2.5_10101057.5102.52.844Ta_3_101010571032.903
GEANT4 and MCNP simulation geometry
GEANT4 toolkit
GEANT4 toolkit is a Monte Carlo package that uses random numbers during the calculation process. GEANT4 code can be used to simulate passing various sorts of particles either charged or non-charged with any range of energy within different types of materials. Moreover, an excellent feature of this code is the flexibility of choosing arbitrary physics lists and input, and output parameters^25,26^. In the present work, GEANT4-10.7, a version of GEANT4 code is employed to simulate the set-up of the X-ray radiations and calculate the shielding parameters.
MCNPX code
MCNPX, an improved Monte Carlo package by Los Alamos National Laboratory, is applied for calculating the shielding features of the CaF_2_–CaO–B_2_O_3_–P_2_O_5_–SrO–Ta_2_O_5_ glass system against photons with an energy range of 20–100 keV. Also, the code is designed to simulate a wide range of energies for different particles (alpha, photon, proton, neutron, deuteron, electron, ^3^He etc.) and various types of materials (elements, compositions, and mixtures). In MCNP code, the geometry can be illustrated in 2D and 3D dimensions, which is one of the best features of the code. The geometry part of the target, source, detector, and tally key can be defined in the input file of the code^26^. In the present study, MCNPX-2.6, a version of the MCNPX code was selected to simulate the set-up of the X-ray radiations.
The geometry details and set-up
Figure 1 displays the geometry of an X-ray machine, which was selected for evaluating of shielding properties of the CaF_2_–CaO–B_2_O_3_–P_2_O_5_–SrO–Ta_2_O_5_ glass system. As shown in Fig. 1, a glass sample prepared as a disk with 2 cm radius and 3 mm thickness is placed in the middle of the photon source with an energy range from 20 to 100 keV and a cylindrical detector with 2.54 cm radius and 2.54 cm thickness. The X-ray source is simulated as a circle plane with a radius of 1.5 mm, with the sample placed 10 cm far from the source and the detector. The physics lists of GEANT4 code that applied for the calculations of this study are G4EmStandardPhysics_option2, G4PhotoNuclearProcess, G4LowEGammaNuclearModel, and G4CascadeInterface.Figure 1. Geometry of X-ray radiation shielding simulated using MCNPX and GEANT4 codes: (a) before irradiating the samples, (b) during irradiating the sample.
Moreover, a simulation setup diagram from visual editor for GEANT4 and MCNP nuclear codes is displayed in Fig. 2. Photons exit from the top (the cubic), drop on the middle cylinder (the sample), and detect by the yellow cylinder (the detector). The flow chart of the overall simulation process of using GEANT4 and MCNP codes are illustrated in Fig. 3. The two Monte Carlo input files are including materials, source, and detector details, which with defining the tally, the cell flux of the photons can be seen in the output of the codes.Figure 2. The simulation setup diagram from Visual Editor for both GEANT4 and MCNP codes.Figure 3. Operation of the simulation process of GEANT4 and MCNP codes.
Phy-X/PSD
The Phy-X/PSD software can determine shielding values for photon energies between 1 keV and 100 GeV for any materials in single elementals or mixture composites. This program illustrates the chemical compositions and density of the materials. Besides, it calculates various shielding parameters such as MAC, LAC, MFP, Z_eff_, HVL, and TVL using theoretical MAC values^27,28^. In this study, Phy-X/PSD software was selected to calculate the shielding features of the CaF_2_–CaO–B_2_O_3_–P_2_O_5_–SrO–Ta_2_O_5_ glass system and compare with it GEANT4 and MCNPX simulated values.
Radiation parameters
Linear attenuation coefficient (LAC)
The linear attenuation coefficient (µ) is the number of particles or photons reflected from the beam per unit thickness of the sample surface. This coefficient is calculated in a unit of cm^−1^ using Eq. (1)^29^:
\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ LAC = - \frac{{\ln \left( {\frac{N}{{N_{0} }}} \right)}}{x} $$\end{document}where LAC is the linear attenuation coefficient (cm^−1^), N and N_0_ are primary and ultimate particle or photon numbers, respectively, and x is the thickness of the glass sample (cm).
Mass attenuation coefficient (MAC)
The mass attenuation coefficient is the probability of photons or particles that interact in a unit of mass of the glass sample. It can be computed using Eq. (2)^30^:
\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ MAC = \frac{LAC}{\rho } = \sum\limits_{i} {W_{i} (MAC)}_{i} $$\end{document}where LAC is the linear attenuation coefficient (cm^−1^), MAC is the mass attenuation coefficient (cm^2^/g), W_i_ is the weight fraction, and ρ denotes the glass sample density (g/cm^3^).
Half-value layer (HVL)
The thickness of the glass sample that cut the power of the particle or photon beam in half is defined as the half-value layer. Equation (3) is used to evaluate HVL^31^:
\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ HVL = \frac{\ln (2)}{{LAC}} $$\end{document}where HVL is the half-value layer (cm) and LAC is the linear attenuation coefficient (cm^−1^).
Tenth-value layer (TVL)
The tenth-value layer of the sample glass is the thickness of the glass that attenuates the intensity of the beam to one-tenth of its primary strength. The TVL is calculated by Eq. (4)^31^:
\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ TVL = \frac{\ln (10)}{{LAC}} $$\end{document}where TVL is the tenth-value layer (cm) and the LAC is the linear attenuation coefficient (cm^−1^).
Mean free path (MFP)
The distance between two successive interactions of the particles or photons through the glass sample is defined as the mean free path, which can be estimated by Eq. (5)^31^:
\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ MFP = \frac{1}{LAC} $$\end{document}where MFP is the mean free path (cm) and LAC denotes the linear attenuation coefficient (cm^−1^).
Effective atomic number
The effective atomic number (Z_eff_) of the samples is expressed by Eq. (6)^32^:
\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ Z_{eff} = \frac{{\sum\nolimits_{i} {f_{i} A_{i} \left( {MAC} \right)_{i} } }}{{\sum\limits_{j} {\frac{{A_{j} }}{{Z_{j} }}\left( {MAC} \right)_{i} } }} $$\end{document}MAC values of the GEANT4 code and Phys-x software were compared using Eq. (7)^32^:
\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ Diff = \left| {\frac{{MAC_{Phy - x} - MAC_{GEANT4orMCNP} }}{{MAC_{Phy - x} }}} \right| \times 100 $$\end{document}Results and discussion
In this research, the shielding parameters of the CaF_2_–CaO–B_2_O_3_–P_2_O_5_–SrO–Ta_2_O_5_ glass were computed for X-ray energies from 20 to 100 keV by the GEANT4 toolkit and MCNPX code. Figure 4 indicates the linear attenuation coefficient simulated results of five samples of the CaF_2_–CaO–B_2_O_3_–P_2_O_5_–SrO–Ta_2_O_5_ glass. According to Fig. 4a, LAC values at 20 keV are 22.07342 and 22.42434066 cm^−1^ for GEANT4 and MCNP codes, respectively, while decline to 0.514242 and 0.553058358 cm^−1^ in the energy of 100 keV. Moreover, by increasing the percentages of Ta_2_O_5_ in samples, the values of LAC rise for the Ta_0_, Ta_1_, Ta_2_, Ta_2.5_, and Ta_3_ glass samples. For example, in the energy of 40 keV, the GEANT4 value of LAC reported in this study for the mentioned samples are 3.348961, 3.816376, 4.091516, 4.223683, and 4.444864 cm^−1^, while they are 3.503685242, 3.963567196, 4.295377408, 4.424974352, and 4.630099796 cm^−1^ for MCNP code, respectively. The measured values for other energies are summarized in Table 2. Due to the higher atomic number of Ta than P, with increasing the amount of Ta_2_O_5_, the density of the samples rises from 2.57 to 2.903 g/cm^3^ for the Ta_0_, Ta_1_, Ta_2_, Ta_2,5_, and Ta_3_ samples. As can be seen from Table 2, increasing the density has raised the value of the LAC parameter. Ta_3_ with the highest density value compared to other samples is the best shield against incident photons. Therefore, the behavior of the LAC is depended on the density and chemical composition of the samples, and the energy of photons.Figure 4. Linear attenuation coefficient values calculated using MCNPX and GEANT4 codes for (a) Ta_0_, (b) Ta_1_, (c) Ta_2_, (d) Ta_2.5_, and (e) Ta_3_ samples of the CaF_2_–CaO–B_2_O_3_–P_2_O_5_–SrO–Ta_2_O_5_ glass.Table 2GEANT4 and MCNPX simulated values of LAC for the samples of the CaF_2_–CaO–B_2_O_3_–P_2_O_5_–SrO–Ta_2_O_5_ glass.Photon energy (keV)LAC simulated values (cm^−1^)Ta_0_Ta_1_GEANT4MCNPGEANT4MCNP2022.0734183322.4243406624.7322845925.33936631307.2831693427.4804685038.236943788.483611604403.3489614533.5036852423.8163758743.963567196501.8888055552.0109021482.1423246792.270125654601.2551220771.330625421.3964302471.494221831700.907367480.9723576861.1899627761.270451968800.7176010420.7659441490.912180740.982095674900.5967907280.638314390.7479600220.8013925821000.5142423230.5530583580.6368331960.681597894Photon energy (keV)LAC simulated values (cm^−1^)Ta_2_Ta_2.5_Ta_3_GEANT4MCNPGEANT4MCNPGEANT4MCNP2027.5944899927.5439262527.8724286828.3813834528.5092794729.61299889308.9368726289.2044379499.2183272189.4972093629.6489962329.936976289404.0915160984.2953774084.223682764.4249743524.4448642214.630099796502.3061772082.4542361132.3642212062.5267343592.4849127642.640002814601.5169477521.6089722391.563075591.6539706221.6208727151.724932673701.466098651.537171621.587629841.6626934351.7166130421.814044526801.1097654131.1759345651.1959502281.2639466711.2930850291.374044541900.8862761470.9444197820.9514153321.0112598971.0290329141.0948866031000.7437784110.7928216570.7937058850.84225380.8520668090.905962664
Furthermore, Fig. 5 presents a 3D comparison of LAC results between the five samples of the CaF_2_–CaO–B_2_O_3_–P_2_O_5_–SrO–Ta_2_O_5_ glass system.Figure 5A comparison of linear attenuation coefficient parameter of the five samples of the CaF_2_–CaO–B_2_O_3_–P_2_O_5_–SrO–Ta_2_O_5_ glass that evaluated by MCNPX and GEANT4 codes.
Figure 6. exhibits mass attenuation coefficient values estimated by MCNPX and GEANT4 codes for the Ta_0_, Ta_1_, Ta_2_, Ta_2,5_, and Ta_3_ samples of the CaF_2_–CaO–B_2_O_3_–P_2_O_5_–SrO–Ta_2_O_5_ glass. According to this figure, MAC values for the mentioned samples decline continuously from 20 to 100 keV. The values of MAC for the GEANT4 toolkit in energy of 40 keV for the Ta_0_, Ta_1_, Ta_2_, Ta_2,5_, and Ta_3_ samples are 1.303098, 1.386764, 1.445765, 1.485121, and 1.531128 cm^2^/g, respectively. Also, the simulation outputs by MAC for MCNP code for the glass samples were reported to be 1.363301651, 1.440249708, 1.517801204, 1.555898155, and 1.594936203 cm^2^/g, respectively, which are in good agreement with the GEANT4 values. Table 3 lists the values of MAC for five mentioned samples obtained using GEANT4 and MCNP codes and Phy-X/PSD program in the energy range of 20 to 100 keV. Table 3 presents the differences between GEANT4 and MCNPX codes, and the Phy-X/PSD program. As can be seen, the outputs of these three software are in good agreement for all energy range and the five samples of the glass system.Figure 6. Mass attenuation coefficient values calculated using MCNPX and GEANT4 codes for (a) Ta_0_, (b) Ta_1_, (c) Ta_2_, (d) Ta_2.5_, and (e) Ta_3_ samples of the CaF_2_–CaO–B_2_O_3_–P_2_O_5_–SrO–Ta_2_O_5_ glass.Table 3GEANT4 and MCNPX simulated values of MAC for the samples of the CaF_2_–CaO–B_2_O_3_–P_2_O_5_–SrO–Ta_2_O_5_ glass.Photon energy (keV)MAC simulated values (cm^2^/g)Ta_0_Ta_1_Ta_2_GEANT4MCNPPhy-X/PSDDiff (%)GEANT4MCNPPhy-X/PSDDiff (%)GEANT4MCNPPhy-X/PSDDiff (%)208.5888787268.7254243828.8613.0677929531.5300261628.9870220179.2076185739.3533.9107060421.5543828379.7507031769.732836139.8450.9563734721.139297821302.8339180322.9106881332.9684.5158048511.9309928082.9930755013.0827077053.1384.6248759651.762023423.1579055223.2524515723.3084.5512657741.679214872401.3030978421.3633016511.4037.1011754712.8295331051.3867644891.4402497081.4826.4368039272.8171587141.4457654061.5178012041.5627.4195295252.829628418500.7349437960.7824521980.8159.7898128253.9935953650.7784610030.824900310.8599.330335443.9696961030.8149036070.8672212410.9029.6994736833.855738202600.4883743490.5177530820.54410.284020444.8248011230.5074237810.5429585140.57111.185230184.9109432040.5360239410.5685414270.59810.408874614.926182765700.3530612760.378349294––0.4323992640.461646791––0.518056060.543170184––800.2792221950.2980327430.31912.537939576.572807860.3314610250.3568661610.38012.711533516.0878524260.3921432560.4155245810.44010.919395085.562595182900.2322142910.248371358––0.2717877990.2912037––0.3131717830.333717237––1000.2000942890.2151978050.23314.023614597.6404272230.2314074110.2476736530.26713.182188577.2383320990.2628192260.2801489950.30012.497348086.61700159Photon energy (keV)MAC simulated values (cm^2^/g)Ta_2.5_Ta_3_GEANT4MCNPPhy-X/PSDDiff (%)GEANT4MCNPPhy-X/PSDDiff (%)209.8004320279.97938939910.0912.878458711.1060410399.82062675610.2008263510.3374.9947603931.317342078303.241324623.3393844453.3944.4876431271.6091795823.3238016643.4230025113.4794.4542147061.609585778401.485120521.5558981551.6017.2589268492.8171046381.5311278751.5949362031.6416.7007464922.807056518500.8313014080.8884438670.92410.068174393.8480663130.8559809730.9094050340.9469.5446122013.868389613600.5496046380.5815649160.61210.163954714.9730528750.5583440290.5941896910.62510.70387144.929649447700.558238340.584632009––0.5913238170.624886161––800.4205169580.4444256930.47010.614319275.4413418910.4454305990.4733188220.50111.037149775.525185175900.3345342230.355576616––0.3544722410.377156942––1000.2790808320.2961511250.31712.034513366.5769321310.2935125070.3120780790.33412.166117696.563449287
Figure 7 compares MAC values for the CaF_2_–CaO–B_2_O_3_–P_2_O_5_–SrO–Ta_2_O_5_ glass samples.Figure 7A comparison of mass attenuation coefficient parameter of the five samples of the CaF_2_–CaO–B_2_O_3_–P_2_O_5_–SrO–Ta_2_O_5_ glass that evaluated by GEANT4 and MCNPX codes.
Figure 8 presents the half-value layer GEANT4 and MCNP codes calculations of the Ta_0_, Ta_1_, Ta_2_, Ta_2,5_, and Ta_3_ samples for photons with energy ranging from 20 to 100 keV. According to this figure, HVL results for each sample start at 20 keV and gently rise to 100 keV. However, the values for samples with a higher of Ta_2_O_5_ content are less than others. For instance, in incident photons with energy of 40 keV, HVL of the five samples are obtained as 0.20693, 0.181586, 0.169375, 0.164075, and 0.15591 cm for the GEANT4 program, respectively. Meanwhile, the MCNP values were obtained as 0.19779174, 0.174842501, 0.161336231, 0.156611077, and 0.1496728 cm, respectively. Table 4 presents the HVL results for other energies of photons.Figure 8. Half-value layer results calculated using GEANT4 and MCNPX codes for (a) Ta_0_, (b) Ta_1_, (c) Ta_2_, (d) Ta_2.5_, and (e) Ta_3_ samples of the CaF_2_–CaO–B_2_O_3_–P_2_O_5_–SrO–Ta_2_O_5_ glass.Table 4GEANT4 and MCNPX simulated values of HVL for the samples of the CaF_2_–CaO–B_2_O_3_–P_2_O_5_–SrO–Ta_2_O_5_ glass.Photon energy (keV)HVL simulated values (cm)Ta_0_Ta_1_GEANT4MCNPGEANT4MCNP200.0313952280.0309039190.0280200560.02734875300.0951508840.0926412560.0841331470.081686908400.2069298230.197791740.1815858870.174842501500.366898540.3446214430.3234803790.305269446600.5521375270.5208077260.4962653890.463786558700.7637478920.7127006960.5823711580.545475168800.9657176620.9047657070.759717860.705633899901.161211071.0856719060.9265201070.8647447151001.3476137021.2530323251.0881970421.016728494Photon energy (keV)HVL simulated values (cm)Ta_2_Ta_2.5_Ta_3_GEANT4MCNPGEANT4MCNPGEANT4MCNP200.0251137090.0251598120.024863280.0244174140.0243078750.023401885300.0775439050.0752897680.0751763290.0729688030.0718209420.069739524400.1693748680.1613362310.1640748230.1566110770.1559102740.1496728500.3004972890.282368920.293119780.2742670580.278883030.26249972600.456838410.4307097310.4433566770.4189917230.4275474520.401754811700.4726830630.4508279960.4364997320.4167936110.4037019310.382019289800.6244562970.5893185050.5794555520.5482826260.5359276340.504350463900.7819233340.7337838670.7283885140.6852837750.6734478460.6329422591000.9317291150.8740931760.8731193920.822792370.8133165060.764932185
Figure 9 compares the HVL values of the Ta_0_, Ta_1_, Ta_2_, Ta_2,5_, and Ta_3_ samples for a photon energy range of 20 to 100 keV.Figure 9A comparison of half-value layer parameter of the five samples of the CaF_2_–CaO–B_2_O_3_–P_2_O_5_–SrO–Ta_2_O_5_ glass that evaluated by MCNPX and GEANT4 codes.
Also, Fig. 10 presents the GEANT4 and MCNPX simulated data of the tenth-value layer for five samples of the CaF_2_–CaO–B_2_O_3_–P_2_O_5_–SrO–Ta_2_O_5_ glass irradiated by photons with an energy range of 20 to 100 keV. According to Fig. 10, in energy of 20 keV, the TVL values start and rise to 100 keV for each sample. According to GEANT4 values, in the thicknesses of 0.687552, 0.603343, 0.562771, 0.545161, and 0.518033 cm for the Ta_0_, Ta_1_, Ta_2_, Ta_2,5_, and Ta_3_ samples, the photon with an energy of 40 keV losses its power to one-tenth of its original one. Meanwhile, these values for the MCNP code are 0.657189483, 0.580937569, 0.536061183, 0.520361229, and 0.497307876 cm, respectively. The TVL values decline by increasing the amount of Ta_2_O_5_. The TVL simulated results for photon energies of 20 to 100 keV and the five samples of the glass are summarized in Table 5.Figure 10. Tenth-value layer simulated measurements using GEANT4 and MCNPX codes for (a) Ta_0_, (b) Ta_1_, (c) Ta_2_, (d) Ta_2.5_, and (e) Ta_3_ samples of the CaF_2_–CaO–B_2_O_3_–P_2_O_5_–SrO–Ta_2_O_5_ glass.Table 5GEANT4 and MCNPX simulated values of TVL for the samples of the CaF_2_–CaO–B_2_O_3_–P_2_O_5_–SrO–Ta_2_O_5_ glass.Photon energy (keV)TVL simulated values (cm)Ta_0_Ta_1_GEANT4MCNPGEANT4MCNP200.104314840.1026823990.093100380.090869877300.3161515250.3078129520.2795436210.271415666400.6875519840.6571894830.6033433730.580937569501.219069421.145050791.0748067811.014298521601.8345507061.7304532571.6489080631.540992806702.5376544152.3680432881.9350059851.812414126803.2087259633.006204952.5242640972.344562911903.858278933.6072899673.0784868512.8732298561004.4776265814.1633673213.6156800683.378216268Photon energy (keV)TVL simulated values (cm)Ta_2_Ta_2.5_Ta_3_GEANT4MCNPGEANT4MCNPGEANT4MCNP200.0834436550.0835968360.082611570.0811301220.0807661620.07775589300.2576499840.2501603150.2497833980.2424485980.2386346760.231718888400.562770630.5360611830.5451605210.5203612290.5180327180.497307876500.9984423940.9382084640.9739296340.9112889470.9266261280.872190393601.5179066581.4310906291.4731117981.392155981.4205835361.334884039701.5705526321.4979362511.4503286821.384852461.3413536061.26931013802.074839481.9580894741.9253184941.8217422821.780691171.675771799902.5980447522.4380949412.4201681612.2769469052.237620452.1030352251003.0957944732.9042913682.9010558382.7338375832.7023527612.541589387
Moreover, Fig. 11 compares the TVL parameter values for the five mentioned samples of the CaF_2_–CaO–B_2_O_3_–P_2_O_5_–SrO–Ta_2_O_5_ glass system.Figure 11A comparison of the TVL GEANT4 and MCNPX simulated values for samples of the CaF_2_–CaO–B_2_O_3_–P_2_O_5_–SrO–Ta_2_O_5_ glass.
Figure 12 exhibits the mean free path values computed by GEANT4 and MCNPX codes for five samples of the CaF_2_–CaO–B_2_O_3_–P_2_O_5_–SrO–Ta_2_O_5_ glass system in the photon energy range between 20 and 100 keV. As can be seen, the MFP values for each sample gradually decrease from 20 to 100 keV. According to Fig. 13, the MFP values decline by increasing the Ta_2_O_5_ content of the samples. The distance between to interacting of photons with energy of 40 keV estimated by GEANT4 toolkit in the Ta_0_, Ta_1_, Ta_2_, Ta_2,5_, and Ta_3_ samples are 0.2986, 0.262029, 0.244408, 0.23676, and 0.224979 cm, and for MCNP code are 0.285413766, 0.252297981, 0.232808414, 0.22599001, and 0.215978066 cm, respectively. Table 6 presents more values of the MFP for photon energy in the range of 20 to 100 keV.Figure 12. Mean free path simulated measurements using GEANT4 and MCNPX codes for (a) Ta_0_, (b) Ta_1_, (c) Ta_2_, (d) Ta_2.5_, and (e) Ta_3_ samples of the CaF_2_–CaO–B_2_O_3_–P_2_O_5_–SrO–Ta_2_O_5_ glass.Figure 13A comparison of the MFP GEANT4 and MCNPX simulated values for samples of the CaF_2_–CaO–B_2_O_3_–P_2_O_5_–SrO–Ta_2_O_5_ glass.Table 6GEANT4 and MCNPX simulated values of MFP for the samples of the CaF_2_–CaO–B_2_O_3_–P_2_O_5_–SrO–Ta_2_O_5_ glass.Photon energy (keV)MFP simulated values (cm)Ta_0_Ta_1_GEANT4MCNPGEANT4MCNP200.0453033590.0445943990.0404329810.039464286300.1373028630.1336814670.1214042520.117874326400.2986000330.2854137660.2620286980.252297981500.5294351220.4972892390.4667826540.440504251600.7967352480.7515263010.7161116730.669244672701.102089311.0284281330.8403624220.787121454801.393531981.3055782211.0962739681.018230735901.6756292491.5666261271.3369698521.2478278721001.9446085161.8081274531.5702699021.467140684Photon energy (keV)MFP simulated values (cm)Ta_2_Ta_2.5_Ta_3_GEANT4MCNPGEANT4MCNPGEANT4MCNP200.0362391190.0363056450.0358777490.0352343640.0350762990.033768954300.1118959660.1086432440.1084795510.1052940880.1036377230.100634234400.2444081790.2328084140.2367602060.225990010.2249787510.215978066500.4336180220.4074587590.4229722660.3957677610.4024286140.378787475600.6592184860.6215147630.6397643250.604605660.6169515910.579732772700.6820823420.6505454480.6298697440.6014337810.5825424690.551254385800.9010913370.8503874530.8361551980.791172620.7733443490.727778445901.12831651.0588511791.0510656770.9888654760.9717862140.9133365941001.3444864571.2613177151.2599125421.1872905771.1736168921.103798246
Figure 14 demonstrates the effective atomic number of the five samples of the CaF_2_–CaO–B_2_O_3_–P_2_O_5_–SrO–Ta_2_O_5_ glass. The values of Z_eff_ have the highest at 20 keV photon energy, which are 23.78, 24.76, 25.72, 26.18, and 26.64 for the Ta_0_, Ta_1_, Ta_2_, Ta_2.5_, and Ta_3_ samples, respectively. In this respect, the Z_eff_ values decrease with increasing photon energy and such that lowest result belongs to 100 keV energy. Also, the Ta_3_ sample with the highest density has the highest Z_eff_ value among the five samples. A comparison of Z_eff_ for the five samples of the CaF_2_–CaO–B_2_O_3_–P_2_O_5_–SrO–Ta_2_O_5_ glass system is given in Fig. 15. As can be noticed, the value of Z_eff_ rise with increasing the density of samples.Figure 14. Effective atomic number for (a) Ta_0_, (b) Ta_1_, (c) Ta_2_, (d) Ta_2.5_, and (e) Ta_3_ samples of the CaF_2_–CaO–B_2_O_3_–P_2_O_5_–SrO–Ta_2_O_5_ glass.Figure 15A comparison of effective atomic number for the five samples of the CaF_2_–CaO–B_2_O_3_–P_2_O_5_–SrO–Ta_2_O_5_ glass.
Conclusion
This research explored the shielding properties of the CaF_2_–CaO–B_2_O_3_–P_2_O_5_–SrO–Ta_2_O_5_ glass system with five different molar fractions of Ta_2_O_5_ using GEANT4 and MCNPX nuclear codes, and PhyX/PSD program. The density of the five samples was different from 2.57 to 2.903 g/cm^3^, depending on the amount of Ta_2_O_5_. The parameter of MAC for each sample was determined using the GEANT4 and MCNP codes for photons with an energy range between 20 and 100 keV and compared with the values of Phy-X/PSD software. In the energy of 40 keV and for the sample of Ta_3_, the values of Phy-X/PSD, GEANT4, and MCNP were estimated as 1.641, 1.531127875, and 1.594936203 cm^2^/g, respectively. In contrast, the Monte Carlo codes have 6.700746492 and 2.807056518% discrepancy with the online software. Besides, the value of MAC rises by increasing the Ta_2_O_5_ content. In the mentioned energy, for Ta_0_, Ta_1_, Ta_2_, Ta_2.5_, and Ta_3_ samples, the MAC values of the Phy-X/PSD program are 1.403, 1.482, 1.562, 1.601, and 1.641 cm^2^/g. Therefore, according to the results, the Ta_3_ sample was selected as the best shield against photons in comparison with the other samples.
The reference list from the paper itself. Each links out to its DOI / PubMed record.
- 1Yasmin S Kamişlioğlu M Sayyed MI Assessment of radiation shielding performance of Li 2O–Ba O–Bi 2O 3–P 2O 5 glass systems within the energy range from 0.081 Me V to 1.332 Me V via MCNP 6 code Optik 202327417052910.1016/j.ijleo.2023.170529 · doi ↗
- 2Ardiansyah A Heryanto H Armynah B Salah H Sulieman A Bradley DA Tahir D Physical, mechanical, optical, and gamma radiation shielding properties of the Ba O-based glass system prepared by the melt-quench technique: A review Radiat. Phys. Chem.202321011105910.1016/j.radphyschem.2023.111059 · doi ↗
- 3Şahin N Bozkurt M Karabul YKılıçMGüven Özdemir Z Low cost radiation shielding material for low energy radiation applications: Epoxy/Yahyali Stone composites Prog. Nucl. Energy 202113510370310.1016/j.pnucene.2021.103703 · doi ↗
- 4Shams T Eftekhar M Shirani A Investigation of gamma radiation attenuation in heavy concrete shields containing hematite and barite aggregates in multi-layered and mixed forms Constr. Build. Mater.2018182354210.1016/j.conbuildmat.2018.06.032 · doi ↗
- 5Aygün B Neutron and gamma radiation shielding Ni based new type super alloys development and production by Monte Carlo Simulation technique Radiat. Phys. Chem.202118810963010.1016/j.radphyschem.2021.109630 · doi ↗
- 6Oğul H Agar O Bulut F Kaçal MR Dilsiz K Polat H Akman FA comparative neutron and gamma-ray radiation shielding investigation of molybdenum and boron filled polymer composites Appl. Radiat. Isot.202319411073110.1016/j.apradiso.2023.11073136812766 · doi ↗ · pubmed ↗
- 7Al-Harbi N Al-Hadeethi Y Samir Bakry A Mechanical and radiation shielding features of bioactive glasses: Si O 2–Na 2O–Ca O–P 2O 5–B 2O 3 for utilization in dental applications J. Non-Cryst. Solids 202155212048910.1016/j.jnoncrysol.2020.120489 · doi ↗
- 8Saleh A Shalaby RM Abdelhakim NA Comprehensive study on structure, mechanical and nuclear shielding properties of lead-free Sn–Zn–Bi alloys as a powerful radiation and neutron shielding material Radiat. Phys. Chem.202219511006510.1016/j.radphyschem.2022.110065 · doi ↗
