# Validation of MCNPX with Experimental Results of Mass Attenuation Coefficients for Cement, Gypsum and Mixture

## Article information

## Abstract

### Background

Shielding properties of compound or mixture is presented in terms of mass attenuation coefficients using Monte Carlo simulation. Mass attenuation coefficients of cement, gypsum and the mixture of gypsum and PbCO_{3} has been investigated using monte carlo MCNPX.

### Materials and Methods

The mass attenuation coefficients of cement, gypsum and the mixture of gypsum and PbCO_{3} were calculated for photon energies 365.5, 661.6, 1,173.2, and 1,332.5 keV energies.

### Results and Discussion

The simulated values of mass attenuation coefficients were compared avaialable experimental results, theoretical values by XCOM and found good comparability of the results.

### Conclusion

Standard simulation geometry used in the present investigation would be very useful for various types of sample for shielding and dosimetry applications.

## Introduction

The use of materials for shielding applications is increasing day-by-day for medical, reactors, accelerators, industries, agriculture, space, etc. Choice of materials selection is dependent upon the requirement for exposure rate reduction, type of source, space constrains and final cost effective analysis. Cement, gypsum and lead are common and very useful material in the field of nuclear engineering for shielding applications. These materials are being used for purpose of construction of building of nuclear facility and direct shielding applications. Various researchers have theoretical and experimental investigations for radiation interaction with these materials and proposed various applications. Experimetal investigation requires not only mateirals, radiation sources and expermental setup of very high finatial budget but also radiation exposure to involved personnel. However, simulation method for investigation of radiation interaction is found radiologically safer, less time consuming, cost effective and applicable for desired energy of radiation. It is found that Monte Carlo simulation is suitable method for investigation of radiation interaction with materials in various literature elsewhere [1]. Mass attenuation coefficient is fundamental property of a material for photon interaction to represent interaction and shielding effectiveness. Mass attenuation coefficient of a element is found to be constant at particular photon energy, whereas mass attenuation coefficient of a compound or mixture depends upon composition of elements. The linear attenuation coefficient, half-value layer thickness, tenth-value layer thickness and effective atomic number are derived parameters from mass attenuation coefficient. This study aimed validation of MCNPX (version 2.4.0) code for investigation of mass attenuation coefficients of compounds and mixtures. The investigation was carried out for interaction of gamma radiation with cement, gypsum and mixture of PbCO_{3} with cement and gypsum at energies 365.5, 661.6, 1,173.2, and 1,332.5 keV. Simulated results were compared with available experimental data and found comparable in literature. The present investigation using standard simulation geometry for MCNPX (version 2.4.0) code would be very useful for studies of photon interraction with various types of samples for shielding and dosimetry applications.

## Materials and Methods

### 1. MCNPX code

In this study, MCNPX (version 2.4.0) Monte Carlo code has been used for investigations on mass attenuation coefficients of cement, gypsum and the mixture of gypsum and PbCO_{3}. MCNPX is a general purpose radiation transport code for modeling the interactions of radiation with materials and also tracks all particles at wide range energies. MCNPX is fully three-dimensional and it utilizes extended nuclear cross section libraries and uses physics models for particle types [2]. Validation of MCNPX code have been studied by considering the different applications in the field of radiation [3–11]. In present simulation, each simulation parameters have been defined in input file by considering the material properties and experimental setup. In this study, an isotropic radiation source has been defined and located at the beginning of the collimator. To avoid the backscatter into detection point, full collimation from radiation source point until the target shielding material was provided. Afterwards, energetic photons have interacted with target shielding material. Transmitted photons and their numerical quantities have been calculated in detection field. The simulation geometry of mass attenuation coefficient calculations can be seen in Figure 1. The geometrical forms and physical parameters for simulation have been defined in cell card and surface card of MCNPX input file. On the other hand, source definition has been defined with different variables such as CEL, ERG, DIR, POS, and PAR. The energy source has been defined as a point source on the central axis as well. To obtain the average cell flux in defined detection area, mesh tally (type F4) has been used. This type of tally in MCNPX code reports the sum of the contributions in considered cell. This type of tally in MCNPX scores energy deposition data in which energy deposited per unit volume from all particles is included. MCNPX calculations were completed by using Intel^{®} Core^{TM} i7 CPU 2.80 GHz computer hardware.

### 2. Theoretical estimation

In present work, following formulae were used to determine the mass attenuation coefficients for compound and composite materials. When a beam of monochromatic gamma photons is attenuated on matter according to Lambert-Beer law:

Where, *I** _{0}* and

*I*are the incident and transmitted photon intensities, respectively,

*μ*(cm

^{−1}) represents linear attenuation coefficient of the material and

*t*(cm) is the thickness of the target material/sample. Rearrangement of Equation 1 yields the following equation for the linear attenuation coefficient:

The mass attenuation coefficients of the materials were calculated using the following equation.

Where, *w** _{i}* is the weight fraction, (

*μ*/

*ρ*) is the mass attenuation coefficient of the material,

*μ*

*is mass attenuation coefficient of element. The*

_{m}*w*

*can be defined as follows:*

_{i}Where, *A** _{i}* is the atomic weight of the sample,

*n*

*is a number of formula units.*

_{i}### 3. Validation

Applying the Monte Carlo method is the one of the best solution for the investigation of different complex material behaviors since experimental duplication of investigation is quite complicated. So, it is more suitable to apply some numerical methods such as Monte Carlo [12]. In this paper, a validation for input code was performed. On the other hand, WinXcom program was also used to calculate the gamma ray mass attenuation coefficients of the studied shielding materials. WinXcom [13,14] program is a user friendly calculation program and input parameter specifications are quite flexible and easy to access. In the WinXcom program, firstly, shielding material types were defined by their elemental mass fractions, which are totally the same as in MCNPX Monte Carlo code input. Secondly, the gamma ray energies have been defined. The attenuation coefficients of the selected materials were finally calculated by the program.

## Results and Discussion

MCNPX simulation input has been used for mass attenuation coefficients calculations of cement, gypsum and the mixture of gypsum and PbCO_{3} materials. The mass attenuation coefficients for attenuator samples doped by different percentages of PbCO_{3} were calculated for the four different energies 356, 662, 1,173, and 1,333 keV and shown in Figure 2. The standard XCOM data has been used for comparison with obtained MCNPX results. Figure 2 shows mass attenuation coefficinents of cement, gypsum and small doping of PbCO_{3} in both the shielding materials. It is found that the mass attenuation coefficients of pure and doped cement or gypsum are decreasing with increase in photon energy. This variation of mass attenuation coefficients can be explained using fundamental photon interaction process of photoelectric effect, compton effect and pair production for low-, intermediate- and high energy photons, respectively, which varies with atomic number of elements of compostions. From Figure 2A and 2C, it is to be noted that the mass attenuation coefficients of 100% cement and 100% gypsum estimated using MCNPX is lesser than theoretical data and comparable with experimental results. However, from Figure 2B and 2D it is observed that using MCNPX simulation, the mass attenuation coefficients of 100% cement and gypsum are lesser than 30% doping of PbCO_{3}. Doped cement and gypsum with PbCO_{3} are found with large difference of mass attenuation coefficients at low energy (photoelectric effect region) compared with high energy (pair production region). It is because of Pb element of PbCO_{3}, which dominant for interaction in photoelectric effect region (interaction cross sectioin α Z^{4–5}).

## Conclusion

Mass attenuation coefficients of cement, gypsum and the mixture of gypsum and PbCO_{3} has been investigated using monte carlo MCNPX. The simulated values of mass attenuation coefficients were compared available experimental results, theoretical XCOM values and found good comparability of the results. Standard simulation geometry for sample prepration of the present investigation would be very useful for various types of photon interaction investigations using MCNPX without experimental analysis for pure and mixture of elements. In addition, radiation safety inside medical areas, research centers, houses, nuclear shelters especially those close to nuclear power stations can be estimated by the calculation of shielding features of common cement and gypsium mixture materials.

## References

_{3}on mass attenuation coefficients of concrete by using MCNPX code. Appl Radiat Isot 2017;121:122–125.

_{2}–5TiO

_{2}–(15–x) WO

_{3}–xA

_{n}O

_{m}glasses using WinXCom and MCNP5 code. Radiat Phys Chem 2017;141:172–178.