Enhancement of Gamma-ray Shielding Properties in Cobalt-Doped Heavy Metal Borate Glasses: The Role of Lanthanum Oxide Reinforcement

The direct influence of La3+ ions on the gamma-ray shielding properties of cobalt-doped heavy metal borate glasses with the chemical formula 0.3CoO-(80-x)B2O3-19.7PbO-xLa2O3: x = 0, 0.5, 1, 1.5, and 2 mol% was examined herein. Several significant radiation shielding parameters were evaluated. The glass density was increased from 3.11 to 3.36 g/cm3 with increasing La3+ ion content from 0 to 2 mol%. The S5 glass sample, which contained the highest concentration of La3+ ions (2 mol%), had the maximum linear (μ) and mass (μm) attenuation coefficients for all photon energies entering, while the S1 glass sample free of La3+ ions possessed the minimum values of μ and μm. Both the half value layer (T1/2) and tenth value layer (TVL) of all investigated glasses showed a similar trend of (T1/2, TVL)S1 > (T1/2, TVL)S2 > (T1/2, TVL)S3 > (T1/2, TVL)S4 > (T1/2, TVL)S5. Our results revealed that the S5 sample had the highest effective atomic number (Zeff) values over the whole range of gamma-ray energy. S5 had the lowest exposure (EBF) and energy absorption (EABF) build-up factor values across the whole photon energy and penetration depth range. Our findings give a strong indication of the S5 sample’s superior gamma-ray shielding characteristics due to the highest contribution of lanthanum oxide.


Introduction
Gamma radiation is commonly used in medicine and in industry [1]. The biological effect of radiation is well known; consequently, massive efforts are being made in the field of radiation shielding. Lead-based shields were commonly used, but there is growing

Materials and Methods
In this study, a group of lanthanum-oxide-reinforced heavy metal borate glasses with the chemical formula 0.3CoO-(80-x)B 2 O 3 -19.7PbO-xLa 2 O 3 : x = 0, 0.5, 1, 1.5, and 2 mol% were investigated in terms of their extensive gamma-ray shielding parameters. Previously, Abul-Magd et al. [13] studied the effect of the rare-earth compound lanthanum oxide on the structural, mechanical, and optical properties of a glassy system composed of cobalt-doped heavy metals. Our objective was to conduct a theoretical follow-up analysis on those glass samples in order to provide new skills and information to the scientific community by analyzing their gamma-ray shielding properties.

Fundamental Shielding Parameters
The linear attenuation coefficient term (µ) for photons of given energy in a certain material is the result of numerous physical processes that result in photon emission from the beam. A material's linear attenuation coefficient is determined by its density. By normalizing the linear attenuation coefficient (µ) with the absorber density (ρ), the dependency on the absorber density is removed. Because it is independent of the absorber density and physical state, the mass attenuation coefficient defined by µ/ρ is more fundamental than the linear coefficient value. The photons are transmitted in accordance with the Beer-Lambert law when a gamma ray interacts with a sample of thickness x (cm) with respect to the narrow beam geometry [14,15].
The I 0 and I values in Equation (1) denote the gamma-ray intensity before and after passing through a sample of thickness x (cm). Further, µ (cm −1 ) is the linear attenuation coefficient of the sample and t (cm) is the physical thickness of the shielding material.
The linear attenuation coefficient can be defined from the point of mass attenuation coefficient as follows: where (µ) s = µ ρ (cm 2 g −1 ) is the mass attenuation coefficient and ρ (g/cm 3 ) is the density of the sample.
The mass attenuation coefficient of a composite or mixture is equal to the total of the weights of each of the constituent elements [16,17].
This expression is the mixing rule, also known as Bragg's rule. In this equation, w i is the weight fraction and (µ/ρ) i is the mass attenuation coefficient of the ith element. The fraction by weight w i of a chemical compound is calculated using the equation below [18,19].
The half value layer (HVL), tenth value layer (TVL), and mean free path (MFP) are the other essential metrics for determining a material's shielding capabilities.
The HVL and TVL refer to the material thicknesses that attenuate one-half and onetenth of the photon strength, respectively. The MFP is the average distance that a photon with a certain energy can travel without any interaction [20,21].
Using the XCOM program, the µ ρ values of the samples were calculated. The direct technique was used to determine Z eff , which is a crucial parameter to be aware of when considering gamma protection properties [22,23]. The Z eff value is the parameter measured to show the meaning of the gamma ray and x-ray material absorption fractions of materials prepared for radiation shielding, and it can be calculated from the following equation [24,25].
Materials 2021, 14, 7703 4 of 17 In this equation, for the ith and jth elements, f i and f j are the respective fractional abundances according to the number of atoms, Z i and Z j denote the respective atomic numbers, and A i and A j are the respective weights [26,27].
Additionally, the effective electron density (N eff ) can be estimated using Z eff as follows: (9) where N A is Avogadro's number and M is the atomic mass of the glass. In all areas where radiation is used or maintained, reliable and highly sensitive information about the EABF and EBF parameters is needed. The EBF value represents the degree of probable air interactions between the radiation source and the detector. EABF, on the other hand, is a parameter that refers to how much energy is absorbed in the substance with which the radiation interacts [28][29][30][31]. The EABF and EBF values of the examined glasses were obtained using the GP fitting approach, as previously reported in several papers [32][33][34]. These two parameters are almost equal in terms of determining the ratio of un-collided/un-scattered photons. The GP fitting parameters obtained by interpolation of the equivalent atomic number (Z eq ) were used to calculate the EABF and EBF factors with the help of the following equations: where Z 1 and Z 2 indicate the atomic numbers of the samples for the ratios R 1 and R 2 , respectively, and where C 1 and C 2 denote the GP fitting parameters for the elements with atomic numbers Z 1 and Z 2 , respectively.

Monte Carlo Simulations and Theoretical Calculations of Gamma-ray Shielding Properties
The mass attenuation coefficients of the S1, S2, S3, S4, and S5 glasses were effectively computed using the general-purpose Monte Carlo tool MCNPX (Los Alamos National Laboratory) [35] (version 2.7.0). First, input data for MCNPX were prepared using the following fundamental components: • Card for a cell; • Card for a surface; • Source data.
Within a lead (Pb) shield block, a point source was positioned (see Figure 1). Following that, the glass specimens' elemental compositions (in weight percent) and densities (in grams per cubic centimeter) were calculated. Alternatively, as illustrated in Figure 1, the simulated point isotropic source can be considered as an extension of the overall gammaray transmission arrangement. The elemental mass fractions of the S1, S2, S3, S4, and S5 glasses under examination are listed in Table 1.
Notably, the MCNPX INPUT file [36][37][38] used to specify the elemental composition of glass specimens contained an M n variable. We were able to determine the importance of photon and electron interactions (IMP: p, e) in the cell using the results of the initial cell description technique. To tally the attenuated gamma rays produced, S1, S2, S3, S4, and S5 glasses were linked to the opposite side of the detection area (F4 Tally-Mesh). This type of tally mesh is advantageous for determining the average photon flux. After each run was repeated, a total of 10 8 particles with varying photon energies were collected for each glass sample. When all simulations were run, the MCNPX model had an uncertainty rate of less than 1%. To assess the consistency of the obtained mass attenuation coefficients from the MCNPX code, we used the well-known Phy-X/PSD [39] online calculation platform. Our findings showed that both results were in good harmony in terms of their quantitative values. Accordingly, the mass attenuation coefficients of the investigated glasses were used for the determination of other critical gamma-ray shielding parameters. of glass specimens contained an Mn variable. We were able to determine the importance of photon and electron interactions (IMP: p, e) in the cell using the results of the initial cell description technique. To tally the attenuated gamma rays produced, S1, S2, S3, S4, and S5 glasses were linked to the opposite side of the detection area (F4 Tally-Mesh). This type of tally mesh is advantageous for determining the average photon flux. After each run was repeated, a total of 10 8 particles with varying photon energies were collected for each glass sample. When all simulations were run, the MCNPX model had an uncertainty rate of less than 1%. To assess the consistency of the obtained mass attenuation coefficients from the MCNPX code, we used the well-known Phy-X/PSD [39] online calculation platform. Our findings showed that both results were in good harmony in terms of their quantitative values. Accordingly, the mass attenuation coefficients of the investigated glasses were used for the determination of other critical gamma-ray shielding parameters.

Results and Discussion
The gamma-ray attenuation characteristics of five glass samples strengthened with lanthanum oxide were explored in this research. Table 1 summarizes the molar and elemental mass fractions and densities of the specimens examined. We used two distinct approaches to determine the linear attenuation coefficients (μ) of the glass samples: Monte Carlo simulations and theoretical calculations. The densities of the S1, S2, S3, S4, and S5 glass specimens are shown in Figure 2. As shown in the graph, the glass density increased from 3.11 g/cm 3 to 3.36 g/cm 3 . The S5 sample, with the highest concentration of lanthanum oxide in its structure, had the highest glass density.

Results and Discussion
The gamma-ray attenuation characteristics of five glass samples strengthened with lanthanum oxide were explored in this research. Table 1 summarizes the molar and elemental mass fractions and densities of the specimens examined. We used two distinct approaches to determine the linear attenuation coefficients (µ) of the glass samples: Monte Carlo simulations and theoretical calculations. The densities of the S1, S2, S3, S4, and S5 glass specimens are shown in Figure 2. As shown in the graph, the glass density increased from 3.11 g/cm 3 to 3.36 g/cm 3 . The S5 sample, with the highest concentration of lanthanum oxide in its structure, had the highest glass density.
Given that the linear attenuation coefficient is a density-dependent characteristic, it is assumed that a relation exists here between density and the linear attenuation coefficient values, and, hence, the amount of lanthanum oxide. The shifting linear attenuation coefficients (cm −1 ) as a function of the incident photon energy are depicted in Figure 3. In the graph, it can be observed that as the photon energy increased, the linear attenuation coefficients quickly decreased, reaching as low as 0.06 MeV in certain cases. Most of the photon-matter interactions occur in the low-energy area, where the photoelectric effect dominates, with cross-sectional changes proportional to Z (4)(5) . The lowest-energy region exhibits the highest linear attenuation coefficients. Furthermore, the incoming photon energy, shown by the symbol E 3.5 , is well defined and inversely proportionate. However, the pair creation process dominates, and the cross section for this process is connected to Z 2 at energies greater than 1.022 MeV; as a result, it was discovered that the µ values increased somewhat in this area. Taking into consideration the variations in chemical composition of the glass specimens, Compton scattering becomes increasingly significant for medium-level energies beyond the energy level of 0.06 MeV. Since there is a linear link between the cross section of Compton scattering and atomic number Z, the µ values of the glasses declined slowly and were stable below 1 MeV. Nevertheless, as the glass density changed gradually, no significant variations in the linear attenuation coefficients were seen. We observed an interesting effect of lanthanum oxide on the photon resistance of glass samples at various energies. Our results reveal that the S5 sample, which contained the highest concentration of lanthanum oxide, had the maximum linear attenuation coefficients for all entering photon energies. This is explained by the S5 sample's glass density, containing the largest amount of lanthanum oxide in the glass structure. Meanwhile, another critical statistic for gamma-ray shielding, namely, the mass attenuation coefficient (µ m ), was calculated. Given that the linear attenuation coefficient is a density-dependent characteristic, it is assumed that a relation exists here between density and the linear attenuation coefficient values, and, hence, the amount of lanthanum oxide. The shifting linear attenuation coefficients (cm −1 ) as a function of the incident photon energy are depicted in Figure 3. In the graph, it can be observed that as the photon energy increased, the linear attenuation coefficients quickly decreased, reaching as low as 0.06 MeV in certain cases. Most of the photon-matter interactions occur in the low-energy area, where the photoelectric effect dominates, with cross-sectional changes proportional to Z (4)(5) . The lowest-energy region exhibits the highest linear attenuation coefficients. Furthermore, the incoming photon energy, shown by the symbol E 3.5 , is well defined and inversely proportionate. However, the pair creation process dominates, and the cross section for this process is connected to Z 2 at energies greater than 1.022 MeV; as a result, it was discovered that the μ values increased somewhat in this area. Taking into consideration the variations in chemical composition of the glass specimens, Compton scattering becomes increasingly significant for medium-level energies beyond the energy level of 0.06 MeV. Since there is a linear link between the cross section of Compton scattering and atomic number Z, the μ values of the glasses declined slowly and were stable below 1 MeV. Nevertheless, as the glass density changed gradually, no significant variations in the linear attenuation coefficients were seen. We observed an interesting effect of lanthanum oxide on the photon resistance of glass samples at various energies. Our results reveal that the S5 sample, which contained the highest concentration of lanthanum oxide, had the maximum linear attenuation coefficients for all entering photon energies. This is explained by the S5 sample's glass density, containing the largest amount of lanthanum oxide in the glass structure. Meanwhile, another critical statistic for gamma-ray shielding, namely, the mass attenuation coefficient (μm), was calculated.   Figure 4 depicts the shifting trend in mass attenuation coefficients as a function of incoming photon energy. When the energy was increased from 0.015 MeV to 15 MeV, the mass attenuation coefficients dropped. This might be a result of difficulties encountered during the absorption of high-energy gamma rays with a large penetration factor. On the other hand, S5 exhibits the greatest m values over the whole gamma-ray energy range studied. This might be explained by the fact that S5 s glass structure had the largest amount of La (see Table 1). Apart from increasing the glass density of S5, our results indicated that around 5.12 weight percent La reinforcement enhanced the elemental characteristics of S5 in terms of density-independent gamma-ray attenuation, namely, the mass attenuation coefficient (µ m ). The half value layer term (also known as the HVL) is significant in radiation shielding research since it allows for the quantification of the material thickness required to halve the initial gamma-ray intensity. This is because radiation studies necessitate that shielding requirements be determined in advance based on the type and energy of the radiation used. As a result, the amount of the half value layer required for each type of prospective shielding material should be determined on the basis of a more complete understanding of gamma-ray attenuation capabilities during the incident gamma ray's contact with the attenuator specimen. The fluctuation trend of the half value layer (cm) values of the examined glasses as a function of incident photon energy is depicted in Figure 5. As expected, the required half value layer grows with increasing gamma-ray energy. This is a frequently seen effect of increased gamma-ray energy and, hence, penetrating dominance of accompanying gamma-ray photons. In another sense, larger shields may be capable of deflecting powerfully penetrating gamma rays. Our findings indicate that the S5 sample meets the absolute minimum requirements for glass thickness. This is yet another strong indication of the S5 sample's superior gamma-ray shielding characteristics due to the highest contribution of lanthanum oxide. Another significant parameter in the radiation shielding field is the tenth value layer (TVL), which is the material thickness needed to attenuate 10% of the initial incident gamma-ray intensity.  The half value layer term (also known as the HVL) is significant in radiation shielding research since it allows for the quantification of the material thickness required to halve the initial gamma-ray intensity. This is because radiation studies necessitate that shielding requirements be determined in advance based on the type and energy of the radiation used. As a result, the amount of the half value layer required for each type of prospective shielding material should be determined on the basis of a more complete understanding of gamma-ray attenuation capabilities during the incident gamma ray's contact with the attenuator specimen. The fluctuation trend of the half value layer (cm) values of the examined glasses as a function of incident photon energy is depicted in Figure 5. As expected, the required half value layer grows with increasing gamma-ray energy. This is a frequently seen effect of increased gamma-ray energy and, hence, penetrating dominance of accompanying gamma-ray photons. In another sense, larger shields may be capable of deflecting powerfully penetrating gamma rays. Our findings indicate that the S5 sample meets the absolute minimum requirements for glass thickness. This is yet another strong indication of the S5 sample's superior gamma-ray shielding characteristics due to the highest contribution of lanthanum oxide. Another significant parameter in the radiation shielding field is the tenth value layer (TVL), which is the material thickness needed to attenuate 10% of the initial incident gamma-ray intensity. As shown in Table 2, the half value layer (T1/2) values of the investigated S5 samples were compared with those of Glass1 (Cr2O3-doped BS glass) [40], Glass2 (obsidian glass doped with CeO2) [41], Glass3 (BaO-TiO2-SiO2-B2O3 glass) [42], Glass4 (Li2O-K2O-B2O3-PbO glass) [43], Glass5 (Li2O-K2O-B2O3-HMO (HMO = SrO/TeO2/PbO/Bi2O3)) [44], Glass6 (xBaO-(0.30-x) MgO-0.10Na2O-0.10Al2O3-0.50B2O3 glass) [45], and standard shielding materials (ordinary concrete: OC [46], and hematite-serpentine concrete: HSC [47]); the T1/2 values of the S5 samples were lower than those of all other samples, even OC and HSC.  As shown in Table 2 [45], and standard shielding materials (ordinary concrete: OC [46], and hematite-serpentine concrete: HSC [47]); the T 1/2 values of the S5 samples were lower than those of all other samples, even OC and HSC. Figure 6 depicts the variation of the TVL in centimeters of the investigated glasses at selected energy levels. From Figure 6, it is clear that the change in TVL for all studied glasses (S1-S5) with low energy was small and their values tended to be close to zero in accordance with the photoelectric effect (PE) cross section dominance. With increasing energy, the values of TVL were enhanced due to the dominance of both processes of Compton scattering (CS) and pair production (PP) interactions. The S1 glass sample with ρ = 3.11 g/cm 3 had the maximum TVL values, while the S5 glass sample with ρ = 3.36 g/cm 3 had the minimum values at all selected energies. Therefore, the TVL values of the investigated glasses exhibited a reverse trend to the µ m . Thus, (TVL) S1 > (TVL) S2 > (TVL) S3 > (TVL) S4 > (TVL) S5 . This result confirms that the S5 sample can be considered superior for gamma-ray shielding among the investigated samples. Gamma rays, upon reaching the atomic structure, are projected to engage progressively with the electrons of attenuator materials. This is a natural way for incident energy to be dispersed. However, the mean distance traveled by the gamma rays between two sequential interactions must be known. Fortunately, the mean free path (MFP) is a convenient parameter that concisely defines the mean distance between two sequential encounters. In other words, a smaller mean free path indicates a better shielding substance.  Figure 6 depicts the variation of the TVL in centimeters of the investigated glasses at selected energy levels. From Figure 6, it is clear that the change in TVL for all studied glasses (S1-S5) with low energy was small and their values tended to be close to zero in accordance with the photoelectric effect (PE) cross section dominance. With increasing energy, the values of TVL were enhanced due to the dominance of both processes of Compton scattering (CS) and pair production (PP) interactions. The S1 glass sample with ρ = 3.11 g/cm 3 had the maximum TVL values, while the S5 glass sample with ρ = 3.36 g/cm 3 had the minimum values at all selected energies. Therefore, the TVL values of the investigated glasses exhibited a reverse trend to the μm. Thus, (TVL)S1 > (TVL)S2 > (TVL)S3 > (TVL)S4 > (TVL)S5. This result confirms that the S5 sample can be considered superior for gamma-ray shielding among the investigated samples. Gamma rays, upon reaching the atomic structure, are projected to engage progressively with the electrons of attenuator materials. This is a natural way for incident energy to be dispersed. However, the mean distance traveled by the gamma rays between two sequential interactions must be known. Fortunately, the mean free path (MFP) is a convenient parameter that concisely defines the mean distance between two sequential encounters. In other words, a smaller mean free path indicates a better shielding substance.  The difference in MFP (cm) values for the tested glasses as a function of incident photon energy is depicted in Figure 7. The MFP is a significant statistic in the field of radiation research, especially in studies of radiation shielding. This is because the findings provide unique information in terms of a more accurate estimate of the mean distance for an adjacent incident gamma ray's interaction with the material environment. As a consequence, one may argue that decreasing the value results in a more attenuating environment for energetic gamma rays. We determined the MFP values of all the glass samples examined in this research. The difference in MFP (cm) values for the tested glasses as a function of incident photon energy is depicted in Figure 7. The MFP is a significant statistic in the field of radiation research, especially in studies of radiation shielding. This is because the findings provide unique information in terms of a more accurate estimate of the mean distance for an adjacent incident gamma ray's interaction with the material environment. As a consequence, one may argue that decreasing the value results in a more attenuating environment for energetic gamma rays. We determined the MFP values of all the glass samples examined in this research. In addition to the superior gamma-ray attenuation parameters obtained for the S5 sample, a similar trend in the MFP values was observed. For the S5 sample, minimal MFP values were found at all photon energies investigated. The fluctuation in the effective atomic number (Zeff) values of the examined glasses as a function of incident photon energy is depicted in Figure 8. In addition to the superior gamma-ray attenuation parameters obtained for the S5 sample, a similar trend in the MFP values was observed. For the S5 sample, minimal MFP values were found at all photon energies investigated. The fluctuation in the effective atomic number (Z eff ) values of the examined glasses as a function of incident photon energy is depicted in Figure 8. The effective atomic number is a critical quantity that offers precise information about the effective atomic number for attenuation at specified energy values. S5 was observed to have the maximum effective atomic number values at all energies examined, The effective atomic number is a critical quantity that offers precise information about the effective atomic number for attenuation at specified energy values. S5 was observed to have the maximum effective atomic number values at all energies examined, as illustrated in Figure 8. This can be explained by the increased amount of lanthanum oxide reinforcement, which increased the glasses' overall atomic number from S1 to S5. As a result, the S5 sample's total atomic number changed significantly due to differences in the glass structure between the reduced (B 2 O 3 ) and increased (L 2 O 3 ) replacements. Additionally, the average atomic number of the two substitute materials can be determined. The variation in the net atomic numbers of the glasses tested is the primary explanation for the differences in their effective numbers. Our results reveal that the S5 sample had the highest Z eff values over the whole range of gamma-ray energy. Additionally, the effective electron density (N eff ), which presents the number of electrons per unit mass, was evaluated for the S1-S5 glasses. The change in the effective electron number (N eff ) values of the S1-S5 glasses is shown in Figure 9 as a function of incident photon energy (MeV). From Figure 9, in the low region of photon energy, Neff changed in a non-monotonic trend with photon energy until reaching a sudden jump near the absorption edge of Pb (0.0880 MeV). This behavior can be attributed to the PE process. In the photon energy zone from 0.1 to 1 MeV, a quick decrease in the Neff values was observed in all investigated samples; this trend is related to the CS process, which dominates in this region. In the energy zone greater than 2 MeV, an increase in the Neff values was observed and attributed to the PP process, which dominates in this region. Our data indicate definitively that the S5 sample, which contained the highest concentration of lanthanum oxide, had a significant advantage in terms of gamma-ray attenuation. A total of 2 mole percent replacement of B2O3 for L2O3 resulted in a 0.23 g/cm 3 density shift. Furthermore, this modification altered the gamma-ray shielding properties of the glass samples evaluated. The term "Build-up factor" is required for an effective evaluation of gamma attenuation and may affect the measurement's quality. Gamma ray measurement is required for nuclear technology since it is used in industry, medicine, agriculture, education, research, and military applications. Additionally, it is necessary for the building of radiation protective structures that protect human health. When gamma radiation travels via shielding material, two types of radiation are produced: un-collided photons and colliding photons. As a result, the accumulation factor is an essential statistic for gamma ray measurement. It is defined as the ratio of the total number of particles at a given point to the total number of particles that have not collided at that location. Figures 10 and 11 0.01 0.

Effective electron number (Neff)
[electrons/g] Energy (MeV) S1 S2 S3 S4 S5 From Figure 9, in the low region of photon energy, N eff changed in a non-monotonic trend with photon energy until reaching a sudden jump near the absorption edge of Pb (0.0880 MeV). This behavior can be attributed to the PE process. In the photon energy zone from 0.1 to 1 MeV, a quick decrease in the N eff values was observed in all investigated samples; this trend is related to the CS process, which dominates in this region. In the energy zone greater than 2 MeV, an increase in the N eff values was observed and attributed to the PP process, which dominates in this region. Our data indicate definitively that the S5 sample, which contained the highest concentration of lanthanum oxide, had a significant advantage in terms of gamma-ray attenuation. A total of 2 mole percent replacement of B 2 O 3 for L 2 O 3 resulted in a 0.23 g/cm 3 density shift. Furthermore, this modification altered the gamma-ray shielding properties of the glass samples evaluated. The term "Build-up factor" is required for an effective evaluation of gamma attenuation and may affect the measurement's quality. Gamma ray measurement is required for nuclear technology since it is used in industry, medicine, agriculture, education, research, and military applications. Additionally, it is necessary for the building of radiation protective structures that protect human health. When gamma radiation travels via shielding material, two types of radiation are produced: un-collided photons and colliding photons. As a result, the accumulation factor is an essential statistic for gamma ray measurement. It is defined as the ratio of the total number of particles at a given point to the total number of particles that have not collided at that location. Figures 10 and 11 illustrate the variation in the EBF and EABF values for S5 glass samples. Gamma rays are absorbed mostly in the low-and high-energy bands, which account for the majority of absorption. On the other hand, at intermediate energies, Compton scattering is the dominant mechanism of photon-matter interaction. As a result, in the low-energy zone, the Compton region has the greatest EBF values. Not only was the S5 sample deficient in terms of overall EBF, it was also the only sample with regional variance in EABF. S5 had the lowest EBF and EABF values across the whole photon energy and penetration depth range, according to our findings.
A closer look at the variation in the EBF values of all glass samples at 15 MeV for 40 MFP is also presented in Figure 12. It can be seen that the EBF values decreased linearly from sample S1 to S5 at 15 MeV for 40 MFP. One may conclude that the number of colliding photons grew as the glass density and lanthanum content rose. This is also explained by the direct influence of increasing lanthanum oxide reducing EBF values and thereby increasing the gamma-ray attenuation capabilities of the examined glasses. However, as seen from Figure 12, the mean differences were not so high. This can be explained by the similar elemental compositions and densities of the glass samples. However, our results indicated that increasing La contribution had a monotonic effect on all gamma-ray shielding parameters.   A closer look at the variation in the EBF values of all glass samples at 15 MeV for 40 MFP is also presented in Figure 12. It can be seen that the EBF values decreased linearly from sample S1 to S5 at 15 MeV for 40 MFP. One may conclude that the number of colliding photons grew as the glass density and lanthanum content rose. This is also explained by the direct influence of increasing lanthanum oxide reducing EBF values and thereby increasing the gamma-ray attenuation capabilities of the examined glasses. However, as seen from Figure 12, the mean differences were not so high. This can be explained by the similar elemental compositions and densities of the glass samples. However, our results indicated that increasing La contribution had a monotonic effect on all gamma-ray shielding parameters. Finally, it is important to discuss the verification of the obtained findings, which were utilized to calculate the gamma-ray shielding parameters for the S1, S2, S3, S4, and S5 glasses. To begin, the mass attenuation coefficients of the examined glasses were computed using the MCNPX general-purpose code, as detailed in earlier sections. To ensure the integrity of our findings, we compared the mass attenuation coefficients obtained via MCNPX to those calculated using standard NIST data, Phy-X/PSD. Table 3 compares the mass attenuation coefficients determined for various energies. As can be seen, there is a high degree of numerical consistency for each particular energy value. However, we discovered some minor variations between the values from MCNPX and Phy-X/PSD. This difference may be related to the tools' conceptual differences, with MCNPX being an input-based simulation code and Phy-X/PSD being an online computation platform. In other words, a simulation environment and associated tools should be created in MCNPX code, together with their associated characteristics (see  Finally, it is important to discuss the verification of the obtained findings, which were utilized to calculate the gamma-ray shielding parameters for the S1, S2, S3, S4, and S5 glasses. To begin, the mass attenuation coefficients of the examined glasses were computed using the MCNPX general-purpose code, as detailed in earlier sections. To ensure the integrity of our findings, we compared the mass attenuation coefficients obtained via MCNPX to those calculated using standard NIST data, Phy-X/PSD. Table 3 compares the mass attenuation coefficients determined for various energies. As can be seen, there is a high degree of numerical consistency for each particular energy value. However, we discovered some minor variations between the values from MCNPX and Phy-X/PSD. This difference may be related to the tools' conceptual differences, with MCNPX being an inputbased simulation code and Phy-X/PSD being an online computation platform. In other words, a simulation environment and associated tools should be created in MCNPX code, together with their associated characteristics (see Figure 1). The dimensions chosen, the variance reduction strategies employed, and the general performance of the computer may all have an effect on the overall simulation. On the other hand, the libraries and physics lists that are utilized may have an effect on the numerical results. As shown in Table 3, however, mass attenuation coefficients were reported with a high degree of correlation. On the other hand, multiple prior investigations have shown the compatibility of the MCNPX simulation code with experimental results [48][49][50][51][52]. As a result, the MCNPX simulation data may be regarded to be verified with standard databases and to have a satisfactory degree of consistency.

Conclusions
In this study, we aimed to evaluate different types of gamma-ray shielding parameters using advanced simulation methods. The hypothesis of a recent simulation study was to observe an enhancement in gamma-ray shielding properties with increasing lanthanum oxide reinforcement in the glass structure. Accordingly, the influence of lanthanum oxide on the gamma-ray shielding properties of cobalt-doped heavy metal borate glasses with chemical formula 0.3CoO-(80-x)B 2 O 3 -19.7PbO-xLa 2 O 3 : x = 0, 0.5, 1, 1.5, and 2 mole% was examined using MCNPX general purpose Monte Carlo code and the Phy-X/PSD online platform. Several significant radiation shielding parameters were evaluated. The glass density was increased from 3.11 to 3.36 g/cm 3 with increasing La 3+ ion content from 0 to 2 mole%. The mass attenuation coefficients for all glasses were evaluated via the MCNPX code and Phy-X/PSD online calculation platform. Our findings showed that both results were in good harmony in terms of their quantitative values. The S5 glass sample, which contained the highest concentration of La 3+ ions (2 mole%), had the maximum linear (µ) and mass (µ m ) attenuation coefficients for all entering photon energies, while the S1 glass sample, free of La 3+ ions, possessed the minimum values of µ and µ m . The half value layer (T 1/2 ), tenth value layer (TVL), and mean free path (MFP) of all investigated glasses showed a similar trend of (T 1/ 2, TVL, MFP) S1 > (T 1/2 , TVL, MFP) S2 > (T 1/2 ,TVL, MFP) S3 > (T 1/2 ,TVL, MFP) S4 > (T 1/2 , TVL, MFP) S5 . Our results revealed that the S5 sample had the highest effective atomic number (Z eff ) values over the whole range of gammaray energy. S5 had the lowest EBF and EABF values across the whole photon energy and penetration depth range. Our simulation findings give a strong indication of the S5 sample's superior gamma-ray shielding characteristics due to its highest contribution of lanthanum oxide. On the other hand, it can be concluded that advanced simulations for radiation transport studies such as Monte Carlo simulations can be utilized for an initial assessment of candidate glass shields to determine their gamma-ray attenuation properties before manufacture.