Developing of Lead/Polyurethane Micro/Nano Composite for Nuclear Shielding Novel Supplies: γ-Spectroscopy and FLUKA Simulation Techniques

In this work, the effect of adding Pb nano/microparticles in polyurethane foams to improve thermo-physical and mechanical properties were investigated. Moreover, an attempt has been made to modify the micron-sized lead metal powder into nanostructured Pb powder using a high-energy ball mill. Two types of fillers were used, the first is Pb in micro scale and the second is Pb in nano scale. A lead/polyurethane nanocomposite is made using the in-situ polymerization process. The different characterization techniques describe the state of the dispersion of fillers in foam. The effects of these additions in the foam were evaluated, Fourier transform infrared spectroscopy (FTIR), scanning electron microscopy (SEM), transmission electron microscopy (TEM), and X-ray diffraction (XRD) have all been used to analyze the morphology and dispersion of lead in polyurethane. The findings demonstrate that lead is uniformly distributed throughout the polyurethane matrix. The compression test demonstrates that the inclusion of lead weakens the compression strength of the nanocomposites in comparison to that of pure polyurethane. The TGA study shows that the enhanced thermal stability is a result of the inclusion of fillers, especially nanofillers. The shielding efficiency has been studied, MAC, LAC, HVL, MFP and Zeff were determined either experimentally or by Monte Carlo calculations. The nuclear radiation shielding properties were simulated by the FLUKA code for the photon energy range of 0.0001–100 MeV.


Introduction
In recent years, the use of regular or nano polymer compounds as nuclear shielding materials has attracted the attention of many researchers due to their importance [1][2][3][4][5][6][7][8][9][10].The qualities of nanopolymer compounds are comparable to those of traditional polymers, and they are ecologically benign [11].Due to polyurethane foam's (PU) properties as an effective insulating material, polyurethane foam may be modified to improve its properties [12].Which are commonly used in many industries, including building engineering, automotive, construction, and furniture [13,14] Polyurethane, often known as PU or PUR, is a polymer made of many organic units connected by urethane molecules.Polyurethanes are excellent insulators; they offer many solutions to the challenges of energy conservation and eco-design.The polyurethane sector is continually looking for methods to decrease

Apparent Density Measurement
The density of the samples was determined in accordance with EN 1602 [26].The samples were shaped like cylinders with a radius of 40 mm and a height of 80 mm.The volume of each sample was measured with a Vernier caliper, and the mass of the samples was weighed with an electronic balance (Weight Scale/Ls 1220 M, Precisa company, Dietikon, Switzerland) to calculate the apparent density.

Compressive Stress Test
To determine the compressive stress, we use a universal testing machine (DVT FU/DLC Tensile/Compression Test Device DEVOTRANS Instrument, Istanbul, Turkey).The RPUF samples had a cylindrical shape with a diameter of 80 mm and a height of approximately 80 mm.According to ASTM D1621-1 [27], the samples were measured at a compression rate of 3 mm/min as specimens.Samples were compressed to 5% of the total thickness.

Thermo-Gravimetric Analysis (TGA)
To evaluate the thermal stability of a material, a thermogravimetric analyzer (TG 209F3, NETZSCH, Selb, Germany) was used with dry nitrogen gas at a flow rate of 60 mL min −1 was used for the analysis.The samples' relative mass loss was measured from 20 • C to 700 • C at a heating rate of 10 • C min −1 .

X-ray Diffraction Analysis (XRD) Analysis
The XRD characterization has been conducted to study the crystallinity of the nanostructured RPUF.A Rigaku (Model: D2 Phaser 2nd Gen-Bruker, Billerica, MA, USA) X-ray diffractometer with a Cu-K (0.154 nm) monochromatic radiation and optimal operating parameters of 30 mA and 30 kV was used to analyze the material.All analyses were performed at room temperature, scanning from 200 to 900 in 2 s.The scanning rate was 1.80 min −1 , with a data gathering step size of 0.050.The powder was uniformly distributed on a double-sided sticky tape before being applied to the sample holder.After that, the holder was placed in the chamber to measure the angular distribution of X-ray diffraction (XRD) peaks.Using the Jade program, the obtained patterns were matched with a corresponding standard file.

Fourier-Transform Infrared Spectroscopy (FTIR) Analysis
Fourier-transform infrared (FTIR) spectra were utilized to determine organic materials.The functionalized Pb was recorded using a Shimadzu dxp 400 (Kyoto, Japan), and the FT-IR spectrum was recorded in the mid-range between 4000 and 400 cm −1 .

TGA Analysis
TGA determines the quantity and the frequency of the weight variation of the samples against temperature and time in a controlled atmosphere (e.g., The purging of nitrogen gas).TGA can be used primarily to investigate the thermal stability (the strength of the material at a given temperature), and oxidative stabilities (the oxygen absorption rate of the material).

Gamma-ray Acquisition Setup
The investigation of shielding properties of Pb MPs, and PbNPs polyurethane composites were performed by using a calibrated p-type Hyper Pure Germanium (HPGe) from Canberra company (Atlanta, GA, USA), its resolution of 7.5% at the 662 keV.A set of standard point sources (Am-241, Ba-133, Cs-137, Eu-152, and Co-60) have energy range from 0.060 keV to 1.408 Mev, Table 1 Practically to obtain a parallel beam of photons plus to minimize the coincidence summing effects [28] and the dead time to be below 1%the source was placed at a distance of 50.86 cm from the detector surface.The investigated sample was placed very close to the detector as shown in Figure 1. standard point sources (Am-241, Ba-133, Cs-137, Eu-152, and Co-60) have energy rang from 0.060 keV to 1.408 Mev, Table 1 Practically to obtain a parallel beam of photons plu to minimize the coincidence summing effects [28] and the dead time to be below 1%th source was placed at a distance of 50.86 cm from the detector surface.The investigated sample was placed very close to the detector as shown in Figure 1.
The obtained spectra were analyzed by Genie 2000 software (Canberra company (At lanta, Georgia, United States) to give the count rate at a particular energy.The signals from the detector were collected through a measuring time with a statistical error of less than 1%.The obtained spectra were analyzed by Genie 2000 software (Canberra company (Atlanta, Georgia, United States) to give the count rate at a particular energy.The signals from the detector were collected through a measuring time with a statistical error of less than 1%.

Theoretical Background
The Radiation Shielding Parameters of interest are as follows: 1.
The linear attenuation coefficient (LAC) is the first step to evaluate the material capability for shielding, using Beer-Lambert's law Equation (1) [29]. 2.
The mass attenuation coefficient (MAC) was calculated by dividing the LAC for a given material by its density (ρ).The MAC of a composition can also be calculated using Equation (2) [25]: where (w i ) is the mass-fraction of the ith constituent element in the sample.Using the XCOM program [30][31][32], the MAC theoretically was calculated and compared with the experimental values.
Using the values of MAC, the other shielding parameters were calculated such as the half-value layer (HVL), tenth-value layer (TVL), the mean free path (MFP) and the effective atomic number Zeff were calculated using Equations ( 3)-( 6).

6.
The effective atomic number Zeff is a crucial radiation interaction metric used to describe how radiation attenuation in composite materials varies depending on the specimen's composition and photon energy.Using Equation (6) given below can be used to calculate the Z eff values [42] where fi, Ai, and Zi are weight fraction, atomic weight, and the atomic number of each constituent element in each specimen, respectively.

FLUKA Simulations
Moreover, to confirm the validity of our results we used FLUKA Simulations beside the XCOM online software (NIST Standard Reference Database 8, National Institute of Standards and Technology: Gaithersburg, MD, USA).FLUKA is a general Monte Carlo (MC) Simulation platform code.It has many requests, for the interaction and transport of photons.It uses the best mathematical models and an accurate, microscopic approach [43,44].FLUKA has a wide range of presentations, design shielding, radiation defense, dosimetry, and radiation detector simulation.The newest version of FLUKA, available on 13 January 2023 can be demanded and downloaded from the FLUKA site [45,46].FLUKA needs an auxiliary code-named Flair.Flair is a user-friendly interface for FLUKA code to simplify the entrance of FLUKA input data, organizing FLUKA and visualization of the output tables and figures.It is completely based on Python mode.The latest flair-2.3-0C code was published on 24 March 2023 and can be downloaded from the flair location site.In this activity, simulation was performed using the following FLUKA and flair factors: 1.
The beam profile is assumed to be rectangular.

2.
Abeam focused on the positive z-direction.
For BLKBODY and Void, a sphere is well-defined with R = 100,000 and R = 10,000 individually.

5.
Right Circular Cylinder code RCC declared for maple-target with HZ = 0.1 and R = 5.0.6.
The MGDRAW.F all-purpose event subroutine interface was started by the USER-DUMP card.7.
EMFCUT construction in a material is set equal to the bottom transport cut-offs in the requested area with the familiarized material.9.
The introduced samples have been described by applying the constructed in MATE-RIAL and COMPOUND cards using C, H, N, O and Pb elements.All elements are pre-defined in the FLUKA default library.10.Simulation was conducted for 10 6 primary particles and the code was run for five cycles.The results were directed to the output binary files of USRBIN and USRBDX.The average and mean values were calculated using the FLAIR RUN mode. Figure 2 shows the used simulated geometry.
Polymers 2023, 15, x FOR PEER REVIEW 4. For BLKBODY and Void, a sphere is well-defined with R = 100,000 and R individually.5. Right Circular Cylinder code RCC declared for maple-target with HZ = 0.1 5.0.6.The MGDRAW.F all-purpose event subroutine interface was started USERDUMP card.7. The VOID district: ASSIGNMA VACUUM VOID.8. EMFCUT construction in a material is set equal to the bottom transport cu the requested area with the familiarized material.9.The introduced samples have been described by applying the constructed in RIAL and COMPOUND cards using C, H, N, O and Pb elements.All elem pre-defined in the FLUKA default library.10.Simulation was conducted for 10 6 primary particles and the code was run cycles.The results were directed to the output binary files of USRBIN and U The average and mean values were calculated using the FLAIR RUN mode.
shows the used simulated geometry.
Pure Lead powders of Pb (5 µm) converted by using Photon-planetary ball 600, Zirconia Milling Balls (10 mL-19 mL), ball to powder ratio 1:1, Toluene Rec prevent collection of particles, at 600 rpm after 5 h of milling, The average final s grinding was 3 nm.

Composite Preparation
Rigid Polyurethane Foams (RPUF) are produced using the following steps (F
5. Right Circular Cylinder code RCC declared for maple-target with HZ = 0.1 5.0.6.The MGDRAW.F all-purpose event subroutine interface was started USERDUMP card.7. The VOID district: ASSIGNMA VACUUM VOID.8. EMFCUT construction in a material is set equal to the bottom transport cu the requested area with the familiarized material.9.The introduced samples have been described by applying the constructed in RIAL and COMPOUND cards using C, H, N, O and Pb elements.All elem pre-defined in the FLUKA default library.10.Simulation was conducted for 10 6 primary particles and the code was run cycles.The results were directed to the output binary files of USRBIN and U The average and mean values were calculated using the FLAIR RUN mode.
shows the used simulated geometry.
Pure Lead powders of Pb (5 µm) converted by using Photon-planetary ball 600, Zirconia Milling Balls (10 mL-19 mL), ball to powder ratio 1:1, Toluene Rec prevent collection of particles, at 600 rpm after 5 h of milling, The average final s grinding was 3 nm.

Composite Preparation
Rigid Polyurethane Foams (RPUF) are produced using the following steps (F 1.The Polyol raw material was mixed with a reactive blowing catalyst and tr tion catalyst, and mixed in a cartoon cup for 30 s with a mechanical stirrer rpm) until uniform dispersions were obtained.Pure Lead powders of Pb (5 µm) converted by using Photon-planetary ball mill PM 600, Zirconia Milling Balls (10 mL-19 mL), ball to powder ratio 1:1, Toluene Rectified to prevent collection of particles, at 600 rpm after 5 h of milling, The average final size after grinding was 3 nm.
The Polyol raw material was mixed with a reactive blowing catalyst and trimerization catalyst, and mixed in a cartoon cup for 30 s with a mechanical stirrer at (2000 rpm) until uniform dispersions were obtained.2.
Pb powder either lead microparticles (PbMPs) or lead nanoparticles (PbNPs) were mixed with the previous mixture for 30 s.

3.
Then the physical blowing agent was added to the mixture and mixed for a few seconds.4.
Finally, isocyanate was added, and mixed for a few seconds.
were cut according to the standards of the different tests performed.The compos polyurethane foams is given in Table 2.

Nano-Lead Preparation
Figure 5c shows the structure of the Pb microparticles, while Figure 5a sho sample created by the ball mill machine.The particle size data were obtained fr counting of more than 300 particles.Particle size and size distribution data can be sented graphically as shown in Figure 5b, which summarizes the particle size d tracted from the TEM images.Based on the particle size data, the particle size distr was plotted.A typical way to present particle size and its distribution is in the fo number-frequency histogram.A histogram is a bar graph that illustrates the frequ occurrence versus the size range, which shows that the average final size after g was 3 nm.The mixture was rapidly transferred into an open mold for foam rising.The RPUF was raised in the vertical direction, and then cured at room temperature for 24 h, samples were cut according to the standards of the different tests performed.The composition of polyurethane foams is given in Table 2.

Nano-Lead Preparation
Figure 5c shows the structure of the Pb microparticles, while Figure 5a shows the sample created by the ball mill machine.The particle size data were obtained from the counting of more than 300 particles.Particle size and size distribution data can be represented graphically as shown in Figure 5b, which summarizes the particle size data extracted from the TEM images.Based on the particle size data, the particle size distribution was plotted.A typical way to present particle size and its distribution is in the form of a number-frequency

Effect of Lead on Cell Morphology and Density
Figure 6a shows the pure polyurethane foam compared with the polyurethane composite with 10%wt fillers.It was shown that the polymer doped by 10% either PbNPs or PbMPs showed a slight change in volume compared to the pure sample volume while using 50% wt PbMPs and 50% PbNPs showed a serious change in volume compared with the pure sample volume, as shown in Figure 6b.Also, in the case of the polyurethane composite with Nano, it was slightly bigger than that doped with micro fillers.This may be caused due to the high surface area of nanoparticles compared with microparticles.This advantage of nanoparticles helps to ensure homogenous distribution through polyurethane foam porosity during preparation and prevents the possibility of aggregation when using microparticle size.

Effect of Lead on Cell Morphology and Density
Figure 6a shows the pure polyurethane foam compared with the polyurethane composite with 10%wt fillers.It was shown that the polymer doped by 10% either PbNPs or PbMPs showed a slight change in volume compared to the pure sample volume while using 50% wt PbMPs and 50% PbNPs showed a serious change in volume compared with the pure sample volume, as shown in Figure 6b.Also, in the case of the polyurethane composite with Nano, it was slightly bigger than that doped with micro fillers.This may be caused due to the high surface area of nanoparticles compared with microparticles.This advantage of nanoparticles helps to ensure homogenous distribution through polyurethane foam porosity during preparation and prevents the possibility of aggregation when using microparticle size.
the pure sample volume, as shown in Figure 6b.Also, in the case of the polyurethane composite with Nano, it was slightly bigger than that doped with micro fillers.This may be caused due to the high surface area of nanoparticles compared with microparticles This advantage of nanoparticles helps to ensure homogenous distribution through polyurethane foam porosity during preparation and prevents the possibility of aggregation when using microparticle size.Figure 7a shows the cell structure of a pure sample of polyurethane.The cell structures of polyurethane foam with the addition of lead nanoparticles and bulk sizes were analyzed using SEM as shown in Figure 7b-e.The use of lead provided the polyurethane with an open structure.The cells exhibited a diameter ranging between 170 and 76 µm.The microstructure of the samples with the weight fraction of fillers used shows that the cell distribution becomes more random and nonuniform compared to that in unreinforced polyurethane foam.Generally, the addition of Pb material caused a decrease in foam aperture size with a great range, which is apparent in 50% Pb, especially for micro size, as shown in Figure 7e.The porosity was high, but the cells lacked good evenness for 50% PbNPs, 10% PbMPs, and 50% PbMPs in ascending order, respectively.Compared to 50% PbNPs, other samples caused smaller cells, and the cells were compact and uniform.The difference was attributed to the surface tension and surfactant properties of the Pb.Moreover, with the addition of 10% PbNPs, the cell size pore decreased slightly, and the apparent density of the mixture correspondingly increased slightly from 34.08 to 35.64 kg/m 3 , as shown in Figure 8. Increasing the content of Pb adversely affected the cell uniformity of the mixture, except in cases of lower concentrations of PbMPs and PbNPs, where the cell size and porosity increased.The lower cell uniformity was attributed to the difference between particle size and percentage of lead, which in turn changed the cell structure of the mixture sample.However, this increase was significant because the large particle size decelerated the nucleation reaction of the foam, increased the crosslinking degree and prepolymer density, and decreased the amount of released gas. Figure 8 shows a high difference in apparent density between PbNPs and PbMPs at different percentage ratios.

Effect of Lead on Compressive Stress
Compressive Stress is an important parameter that impacts the performance characteristics of the foam, and the change in its value is presented in Figure 9a-c.The compressive strength of all materials tested was measured in a direction perpendicular to the direction of foam rise.In the current study, rigid polyurethane foam was weakened with different contents of Pb nano and microparticles (10 and 50 wt%).A compression test was performed in the range of 5% of the total sample thickness.The results showed that the addition of Pb nanoparticles impaired the compression properties of the foam.The largest increase in compressive strength was observed for the samples that contained 10 wt% fillers, and their load was about 960 kPa in the perpendicular direction.As we increase the percentage of Pb fillers from 10 wt% to 50 wt%, there is a decrease in compressive strength compared with pure foam; the compressive strength decreases from 990 to 140 kPa in the perpendicular direction.Nonetheless, similar to their Pb micro filler counterparts, the compressive strength of PU-modified foams is higher compared to that of Pb nanofiller foams.The microstructure of the samples with the weight fraction of fillers used shows that the cell distribution becomes more random and nonuniform compared to that in unreinforced polyurethane foam.The loss of uniform microstructure may result in scattered areas that are weak in compression.Once these areas begin to fail in compression, the whole cross-section progressively buckles under the compressive load.

Effect of Lead on Compressive Stress
Compressive Stress is an important parameter that impacts the performance charac teristics of the foam, and the change in its value is presented in Figure 9a-c.The compres sive strength of all materials tested was measured in a direction perpendicular to the di rection of foam rise.In the current study, rigid polyurethane foam was weakened wit different contents of Pb nano and microparticles (10 and 50 wt%).A compression test wa performed in the range of 5% of the total sample thickness.The results showed that th addition of Pb nanoparticles impaired the compression properties of the foam.The larges increase in compressive strength was observed for the samples that contained 10 wt% fillers, and their load was about 960 kPa in the perpendicular direction.As we increas the percentage of Pb fillers from 10 wt% to 50 wt%, there is a decrease in compressiv strength compared with pure foam; the compressive strength decreases from 990 to 14 kPa in the perpendicular direction.Nonetheless, similar to their Pb micro filler counter parts, the compressive strength of PU-modified foams is higher compared to that of P nanofiller foams.The microstructure of the samples with the weight fraction of fillers used shows that the cell distribution becomes more random and nonuniform compared to tha

Thermo-Gravimetric Analysis (TGA)
A thermogravimetric analysis for RPUF was studied with and without the addition of Pb filler.The heating rate was set to 10 • C/min and measured in a nitrogen atmosphere.The samples were heated to 500 • C.There was a two-stage degradation.Urethane linkages were broken down at lower temperatures (first stage) with the evolution of isocyanate and alcohol.Urea linkages and polyol bonds were cleaved when enough energy was provided at high temperature (Second stage) forming complex products.When Pb fillers were added to polyurethane foam with different content ratios, the weight loss changed as the temperature increased.The results are shown in Figure 10, and Table 3 shows the changes in the decomposition curves.The initial decomposition temperature (IDT) for RPUF was 240 • C when the weight loss was 8%.After Pb filler particles were added to RPUF, and as the concentration of Pb filler increased, IDT increased considerably to 300 • C (PbMPs 50% composite) when the weight loss of the hybrid material was 8%.As the concentration of Pb filler particles increased, the remaining mass increased after completely degrading the substance, including residual carbon, particularly at 500 • C. Figure 10 shows the large differences in weight retention between different samples, revealing that Pb filler increased the thermal stability of RPUF.The derivative thermogravimetric curves are displayed in Figure 10.As shown, the control group exhibited the lowest carbon residue amount.Moreover, RPUF containing 50% nano-and micro-pb showed the greatest amount of carbon residue.As shown in Table 3, this indicates that the addition of Pb has a significant effect on the degradation of RPUF and leads to an increase the thermal stability.In the polyurethane foam without the addition of Pb, the surface of the carbon residue was continuous.RPUF with 50% Pb shows the highest thermal stability due to its high initial decomposition temperature and has the lowest rate of degradation.The above proved that the addition of Pb filler nanoparticles effectively improved the thermal stability of hybrid materials.in unreinforced polyurethane foam.The loss of uniform microstructure may result in scattered areas that are weak in compression.Once these areas begin to fail in compression, the whole cross-section progressively buckles under the compressive load.

Thermo-Gravimetric Analysis (TGA)
A thermogravimetric analysis for RPUF was studied with and without the addition of Pb filler.The heating rate was set to 10 °C/min and measured in a nitrogen atmosphere.The samples were heated to 500 °C.There was a two-stage degradation.Urethane linkages were broken down at lower temperatures (first stage) with the evolution of isocyanate and alcohol.Urea linkages and polyol bonds were cleaved when enough energy was provided at high temperature (Second stage) forming complex products.When Pb fillers were added to polyurethane foam with different content ratios, the weight loss changed as the temperature increased.The results are shown in Figure10, and Table 3 shows the changes in the decomposition curves.The initial decomposition temperature (IDT) for RPUF was 240 °C when the weight loss was 8%.After Pb filler particles were added to RPUF, and as the concentration of Pb filler increased, IDT increased considerably to 300 °C (PbMPs 50% composite) when the weight loss of the hybrid material was 8%.As the concentration of Pb filler particles increased, the remaining mass increased after completely degrading the substance, including residual carbon, particularly at 500 °C.Figure 10 shows the large differences in weight retention between different samples, revealing that Pb filler increased the thermal stability of RPUF.The derivative thermogravimetric curves are displayed in Figure 10.As shown, the control group exhibited the lowest carbon residue amount.Moreover, RPUF containing 50% nano-and micro-pb showed the greatest amount of carbon residue.As shown in Table 3, this indicates that the addition of Pb has a significant effect on the degradation of RPUF and leads to an increase the thermal stability.In the polyurethane foam without the addition of Pb, the surface of the carbon residue was continuous.RPUF with 50% Pb shows the highest thermal stability due to its high initial decomposition temperature and has the lowest rate of degradation.The above proved that the addition of Pb filler nanoparticles effectively improved the thermal stability of hybrid materials.

XRD Analysis
Polyurethane foam with and without fillers was analyzed using X-ray diffraction (XRD), a non-destructive analytical method for identifying the atomic and molecular structure of materials.For the pure polyurethane foam, the XRD pattern in Figure 11 shows a wide diffraction peak without any discernible sharp peaks, indicating an amorphous molecular structure.This microstructural trait fits with the PUE material's translucent property.The wide diffraction peak in this material shows up at roughly 20 degrees, which is in line with what has been observed in other polyurethane materials.It was shown that the broad peak that indicates the presence of polyurethane foam slightly decreased when 10% Pb was added, while the peak passed out of sight at 50% Pb was added for both sizes.As can be seen, several remarkable peaks at 29.09, 32.30, 38.09, 49.53, and 63.63 • C are attributed to the Pb substrate [47].The resultant peaks of 50% micro Pb are sharper and have a width at half maxima lower than that of 50% Nano Pb.

023, 15, x FOR PEER REVIEW
63.63 °C are attributed to the Pb substrate [47].The resultant peaks sharper and have a width at half maxima lower than that of 50% Na

FTIR Analysis
Polyurethane rigid foam is a kind of high polymer with a basic repetitive unit urethane bond (−NHCOO−) and can be produced from the polymerization of isocyanates and polyols.The FTIR spectroscopy result of rigid PU foam in air at room temperature (RT) is shown in Figure 12. (Spectral range from (4000 to 500 cm −1 ).For isocyanates, symmetrical and asymmetric segments Figure 12. stretching vibrations of N-H correspond to the broad absorption bands near 3361.17cm −1 and 3289.80 cm −1 , while the medium-strong peak at 1508.98 cm −1 corresponds to the in-plane bending vibration of N − H.The sharp absorption peaks at around 1706.

Fluka Results
This work aimed to investigate the ionizing radiation shielding features of Lead yurethane Micro/Nano Composite.The mass attenuation coefficients (µ/ρ) were g ated by XCOM software over an extended photon energy range from 0.0595 MeV u 1.4080 MeV.The obtained data were verified via the calculated data from FLUKA M Carlo simulation code The obtained values of µ/ρ were used to evaluate several impo shielding parameters such as half value layer (HVL), mean free path (MFP), effe atomic number (Zeff).The calculated MAC of lead free and 10-50% samples coded S1 respectively, at different incident energies is plotted in Figure 13.The FLUKA simula can predict more results for the energy range 0.0001-100 MeV.The benefit of u FLUKA is that the energy interval can be distributed to any portion to give detailed re that are more precise.It can predict all the peaks (tooth) inside the examiner range.In calculation, 5000 portion intervals are preferred for the range from 0.1 keV instance ton energy higher to 100 MeV.This feature shows that the region ranging from 0.1 ke to 100 MeV has three tooth peaks with photon energies 0.0025, 0.014 and 0.1 MeV samples' results overlapped near 1-2.8MeV photon energy.It is clear that the samp Pb free and 10% Pb has the lowest MCA value, while the 50% Pb sample has the hig one.The XCOM calculations confirm the FLUKA results as shown in Table 4.

Fluka Results
This work aimed to investigate the ionizing radiation shielding features of Lead/ Polyurethane Micro/Nano Composite.The mass attenuation coefficients (µ/ρ) were generated by XCOM software over an extended photon energy range from 0.0595 MeV up to 1.4080 MeV.The obtained data were verified via the calculated data from FLUKA Monte Carlo simulation code The obtained values of µ/ρ were used to evaluate several important shielding parameters such as half value layer (HVL), mean free path (MFP), effective atomic number (Zeff).The calculated MAC of lead free and 10-50% samples coded S1-S6, respectively, at different incident energies is plotted in Figure 13.The FLUKA simulations can predict more results for the energy range 0.0001-100 MeV.The benefit of using FLUKA is that the energy interval can be distributed to any portion to give detailed results that are more precise.It can predict all the peaks (tooth) inside the examiner range.In this calculation, 5000 portion intervals are preferred for the range from 0.1 keV instance photon energy higher to 100 MeV.This feature shows that the region ranging from 0.1 keV up to 100 MeV has three tooth peaks with photon energies 0.0025, 0.014 and 0.1 MeV.All samples' results overlapped near 1-2.8MeV photon energy.It is clear that the sample of Pb free and 10% Pb has the lowest MCA value, while the 50% Pb sample has the highest one.The XCOM calculations confirm the FLUKA results as shown in Table 4.
ton energy higher to 100 MeV.This feature shows that the region ranging from 0.1 keV up to 100 MeV has three tooth peaks with photon energies 0.0025, 0.014 and 0.1 MeV.All samples' results overlapped near 1-2.8MeV photon energy.It is clear that the sample of Pb free and 10% Pb has the lowest MCA value, while the 50% Pb sample has the highest one.The XCOM calculations confirm the FLUKA results as shown in Table 4.

Experimental Results of Attenuations Coefficients
Based on the test results of the experiments, the linear attenuation coefficient (LAC), was calculated for all investigated composite samples, using the sample densities to calculate the mass attenuation coefficient (MAC), these results were tabulated in Tables 4  and 5. To confirm the validity of these experimental results.We applied XCOM and FLUKA programs to Pb-free samples Table 4.Where the mass attenuation coefficient (MAC) was simulated by FLUKA and XCOM online software and compared with the ex-

Experimental Results of Attenuations Coefficients
Based on the test results of the experiments, the linear attenuation coefficient (LAC), was calculated for all investigated composite samples, using the sample densities to calculate the mass attenuation coefficient (MAC), these results were tabulated in Tables 4  and 5. To confirm the validity of these experimental results.We applied XCOM and FLUKA programs to Pb-free samples Table 4.Where the mass attenuation coefficient (MAC) was simulated by FLUKA and XCOM online software and compared with the experimental MAC values.The deviation, or difference percentage, between the theoretical and experimental results, was given by the following Equation ( 7).

Experimental Results of Attenuations Coefficients
Based on the test results of the experiments, the linear attenuation coefficient (LAC), was calculated for all investigated composite samples, using the sample densities to calculate the mass attenuation coefficient (MAC), these results were tabulated in Tables 4 and 5.
To confirm the validity of these experimental results.We applied XCOM and FLUKA programs to Pb-free samples Table 4.Where the mass attenuation coefficient (MAC) was simulated by FLUKA and XCOM online software and compared with the experimental MAC values.The deviation, or difference percentage, between the theoretical and experimental results, was given by the following Equation (7).For the XCOM and the FLUKA, the maximum absolute relative deviation from the experimental values was 6.61% and 7.82 at 0.1218 MeV, while the minimum deviation was 0.32% and 0.46 at 0.9641 MeV, which indicated a suitable agreement between the experimental and XCOM and FLUKA results, and this is a prediction of the experimental accuracy.
The slight nonconformity between the theoretical and experimental results has been the focus and interest of many researchers [50].It can be due to the effect of the collimator size, hole, and absorber thickness.For gamma rays, the accurate measurements of MAC (cm 2 g −1 ) ideally need perfect narrow beam irradiation geometry.However, usually, the practical geometries used for the experimental investigations deviate from perfect narrowness, thereby the multiple scattered photons cause systematic differences in the measured values of MAC.Good agreement between theoretical and measured values of MAC was observed for all absorbers with a thickness less than or equal to 0.5-1 mean free paths, thus considered as optimum thickness for low-Z materials in the selected energy range.Sometimes it is very difficult to prepare samples with these specifications in order to satisfy these measurements.This leads us to use the available samples in sizes and volumes that we can suitably manufacture.
On the other hand, FLUKA uses different libraries: BROND-3.1,ENDF-VIII0, JEFF-3.3, and JENDEL-4.0.The photon interaction libraries contain records to define the interaction of photons with the elements Z = 1-100 over the energy range up to 100 MeV.These libraries have been designed to meet the needs of users to model the interaction and transport of primary photons.It should be mentioned that using the default JENDEL library through the simulation process may be one of the reasons for such deviations, especially in the low-Z mode.
In Table 5 and Figures 17 and 18, the experimental results were discussed at different energies, from 0.0595 MeV up to 1.4080 MeV, and it was clear that the highest MAC corresponded to low energy and the values decreased gradually until they reached the highest energy.For low energy photons < 150 keV where the predominant interactions are due to photon electric effect, due to the presence of Pb one will notice a remarkable change in the values of LAC in the energy region (88 to 120 keV) because Pb has K-edge at 88 keV energy.While in the energy region greater the 150Kev the Compton scattering is the predominant interaction.For greater energies, the effect of Pb on the photon attenuation properties of polyurethane decreases due to the increase in the penetration of energetic photons and the increase in mean free path.The experimental MAC results for other composites (containing micro and nano concentration) were determined and the micro and nanocomposites were compared gether as shown in Table 5.The relative difference between these results is given by following Equation (8).: From the MAC values and density, the LAC results were estimated for both m and compositions as shown in Tables 4 and 5 and the relative deviation between th results is given by the following Equation (9).: Dev2 % = lAC nano − lAC micro lAC micro ×100 It is clear that the MAC values for nano are greater than those of composites at energies with all Pb percentages (10, 20, 30, 40, and 50%), as shown in Table 5 and Fig 18 .Where the relative deviation (Dev1) was always positive, these results were due to greater distribution of Pb nanoparticles in the matrix rather than the Pb micropartic since the better distribution enhanced the attenuation ability of the composite.In ot words, the LAC values are affected and depend on the density of the composite rat than any other factor, so we can see the results of 10% and 20%.LAC values have the sa characteristics as MAC values due to the small difference in density between these composites, while in 30%, 40%, and 50% composites, the difference in density between   The experimental MAC results for other composites (containing micro and nano Pb concentration) were determined and the micro and nanocomposites were compared together as shown in Table 5.The relative difference between these results is given by the following Equation (8).: From the MAC values and density, the LAC results were estimated for both micro and compositions as shown in Tables 4 and 5 and the relative deviation between these results is given by the following Equation ( 9).: It is clear that the MAC values for nano are greater than those of composites at all energies with all Pb percentages (10, 20, 30, 40, and 50%), as shown in Table 5 and Figure 18.
Where the relative deviation (Dev1) was always positive, these results were due to the greater distribution of Pb nanoparticles in the matrix rather than the Pb microparticles, since the better distribution enhanced the attenuation ability of the composite.In other words, the LAC values are affected and depend on the density of the composite rather than any other factor, so we can see the results of 10% and 20%.LAC values have the same characteristics as MAC values due to the small difference in density between these two composites, while in 30%, 40%, and 50% composites, the difference in density between the two composites is clear, with the density of micro-composites being greater than that of nano-composites.For this reason, we can see that the LAC values for nano are smaller than those for micro composites at all energies with Pb percentages (30, 40, and 50%), as shown in Table 5   The half-value layer (HVL) and mean free path (MFP) of the pure, 10% and 50% le size percentages were calculated and tabulated in Table 6.The results indicated that HVL and MFP values increased with increasing photon energy as shown in Table 6 a Figures 20 and 21.The values of HVL and MFP decrease with increasing lead percenta in the matrix polymer composites; for example, the value of HVL at 5 microlead/polymer is the lowest value of the corresponding micro-composites at discussed energies.In other words, the half-value layer and mean free path values of composites impeded with nano-lead are lower than those of the micro-composites in energy ranges.The half-value layer (HVL) and mean free path (MFP) of the pure, 10% and 50% lead size percentages were calculated and tabulated in Table 6.The results indicated that the HVL and MFP values increased with increasing photon energy as shown in Table 6 and Figures 20 and 21.The values of HVL and MFP decrease with increasing lead percentages in the matrix polymer composites; for example, the value of HVL at 50% microlead/polymer is the lowest value of the corresponding micro-composites at all discussed energies.In other words, the half-value layer and mean free path values of the composites impeded with nano-lead are lower than those of the micro-composites in all energy ranges.Figure 22 shows the Zeff values of prepared polymer composites against photon energy.It is clear that the values of the effective atomic number (Zeff) decreased with increasing the photon energy until the energy of 0.122 MeV of the 152Eu γ-ray source affected by the k-edge of the lead element at ~0.088 MeV increased at this point and then decreased gradually with increasing the photon energy.The Zeff values also increased with increasing the Pb concentration inside the present polymer, where the Zeff of 50% was greater than 40%, etc.The maximum value was 48.43 of 50% lead/polymer at 0.0595 MeV, while the lowest Zeff value was 3.90 of pure polymer at 1.4081 MeV.The values of HVL and MFP decrease with increasing lead percentages in the matrix polymer composites.Figure 22 shows the Zeff values of prepared polymer composites against photon ergy.It is clear that the values of the effective atomic number (Zeff) decreased with creasing the photon energy until the energy of 0.122 MeV of the 152Eu γ-ray source fected by the k-edge of the lead element at ~0.088 MeV increased at this point and t decreased gradually with increasing the photon energy.The Zeff values also increa with increasing the Pb concentration inside the present polymer, where the Zeff of was greater than 40%, etc.The maximum value was 48.43 of 50% lead/polymer at 0.0 MeV, while the lowest Zeff value was 3.90 of pure polymer at 1.4081 MeV.The value HVL and MFP decrease with increasing lead percentages in the matrix polymer com sites.

Conclusions
A lead/polyurethane composite is made using the in-situ polymerization process s cessfully with different concentrations after the preparation of nano lead by a ball m machine with specific conditions, which is supported by XRD results; 50% PbMPs a 50% PbNPs density and morphology tests showed a serious change in the volume of pure sample and the pore size compared with the control sample.Generally, the addit of Pb material caused a decrease in foam aperture size with a great range, which increa the densities, especially for the micro size.The cell distribution becomes more rand and nonuniform causing impaired the compression properties of the foam.TGA indica that the addition of Pb increases the thermal stability.Results of FTIR generally show no big difference in peak position, but the resultant peaks of 50% microparticles Pb w

Conclusions
A lead/polyurethane composite is made using the in-situ polymerization process successfully with different concentrations after the preparation of nano lead by a ball mill machine with specific conditions, which is supported by XRD results; 50% PbMPs and 50% PbNPs density and morphology tests showed a serious change in the volume of the pure sample and the pore size compared with the control sample.Generally, the addition of Pb material caused a decrease in foam aperture size with a great range, which increased the densities, especially for the micro size.The cell distribution becomes more random and nonuniform causing impaired the compression properties of the foam.TGA indicates that the addition of Pb increases the thermal stability.Results of FTIR generally showed

Figure 2 .
Figure 2. The simulation setup for FLUKA code.Samples were modelled as a cylinder, 10 c eter, with various thicknesses 1-5 mm.

Figure 2 .
Figure 2. The simulation setup for FLUKA code.Samples were modelled as a cylinder, 10 cm diameter, with various thicknesses 1-5 mm.

Figure 2 .
Figure 2. The simulation setup for FLUKA code.Samples were modelled as a cylinder, 10 eter, with various thicknesses 1-5 mm.

Figure 8 .
Figure 8. Foam Density comparison with different percentages of nano and micro lead filler.

Figure 9 .
Figure 9. (a) Compressive Strength of RPUFs with Micro lead, (b) Compressive Strength of RPUFs with Nano lead, (c) Compressive Strength of RPUFs with different percentages of nano and micro lead filler.

Figure 9 .
Figure 9. (a) Compressive Strength of RPUFs with Micro lead, (b) Compressive Strength of RPUFs with Nano lead, (c) Compressive Strength of RPUFs with different percentages of nano and micro lead filler.

Figure 11 .
Figure 11.XRD-curves of RPUF with and without lead additions.

Figure 11 .
Figure 11.XRD-curves of RPUF with and without lead additions.

14 Figure 12 .
Figure 12.FTIR curves of FPUF with and without lead addition.

Figure 12 .
Figure 12.FTIR curves of FPUF with and without lead addition.

Figure 14
Figure 14 shows the dependence of the LAC (µ) with the incident photon energy (MeV).Its value ranges from 1500-~0.01 (cm −1 ) for the energy range 1 × 10 −3 up to 100 MeV.The obtained curves show the same trends as MAC results.Moreover, the mean free path values MFP of the proposed samples are given.It takes values 0.07-150 m for the energy range of 0.0001 up to 100 MeV.The HVL (T 1/2 ) and the TVL (T 1/10 ), Figure 15, of the S1-S6 samples almost take the values 4.84 × 10 −4 -6.6 × 10 −1 and 1.2 × 10 −3 -5.78 × 10 2 , respectively.The effective atomic number, Z eff of the proposed samples are given in Figure 16.It is clear that samples have a peak value near 0.1 MeV which does not appear for the Pb-free sample.Z eff takes values (6.1-434 for lead free), (40.71-5.75 for 10% Pb sample) up to (73.12-15.38 for 50% Pb sample).

Figure 14
Figure 14 shows the dependence of the LAC (µ) with the incident photon energy (MeV).Its value ranges from 1500-~0.01 (cm −1 ) for the energy range 1 × 10 −3 up to 100 MeV.The obtained curves show the same trends as MAC results.Moreover, the mean free path values MFP of the proposed samples are given.It takes values 0.07-150 m for the energy range of 0.0001 up to 100 MeV.The HVL (T1/2) and the TVL (T1/10), Figure 15, of the S1-S6 samples almost take the values 4.84 × 10 −4 -6.6 × 10 −1 and 1.2 × 10 −3 -5.78 × 10 2 , respectively.The effective atomic number, Zeff of the proposed samples are given in Figure 16.It is clear that samples have a peak value near 0.1 MeV which does not appear for the Pb-free sample.Zeff takes values (6.1-434 for lead free), (40.71-5.75 for 10% Pb sample) up to (73.12-15.38 for 50% Pb sample).

Figure 16 .
Figure 16.Variation of Z eff against photon energy for S1-S6 samples.

Figure 17
Figure 17  displays the result of the LAC (cm −1 ) of the lead-free sample, which was evaluated by interring the density factor.

Figure 19 .
Figure 19.The LAC of pure and 10% and 50% Lead/Polyurethane micro and nano Composites.

Figure 19 .
Figure 19.The LAC of pure and 10% and 50% Lead/Polyurethane micro and nano Composites.

Table 1 .
Photon energies, and half-life time for all radioactive sources for all radionuclides used in this work.

Table 1 .
Photon energies, and half-life time for all radioactive sources for all radionuclides used i this work.

Table 4 .
Experimental and theoretical values of the MAC of a lead-free sample.
* Photon energy taken from the IAEA live chart of the nuclides site[49].

Table 4 .
Experimental and theoretical values of the MAC of a lead-free sample.
* Photon energy taken from the IAEA live chart of the nuclides site[49].

Table 5 .
The attenuation results of the Lead/Polyurethane micro and nano Composites.
The Density of Micro-Pb and Nano-Pb were (0.0573 and 0.0484 g•cm −3 )
The MAC of pure and 10% and50% Lead/Polyurethane micro and nano Composites.The MAC of pure and 10% and50% Lead/Polyurethane micro and nano Composites.

Table 6 .
The half value layer and mean free path of the pure, 10% and 50% lead/polymer composites as a function of photon energy.