Measurement of mass attenuation coefficients, effective atomic numbers, and electron densities for different parts of medicinal aromatic plants in low-energy region

Measurement of mass attenuation coefficients, effective atomic numbers, and electron densities for different parts of medicinal aromatic plants in low-energy region 1 first-author 2 corresp Corresponding author, fakman@bingol.edu.tr (F. Akman) Corresponding author, fakman@bingol.edu.tr (F. Akman) fakman@bingol.edu.tr 3 4 Department of Physics, University of Tabuk, 009662 Tabuk, Saudi Arabia. Department of Physics, University of Tabuk, 009662 Tabuk, Saudi Arabia. Department of Electronic Communication Technology, Vocational School of Technical Sciences, Bingöl University, 12000 Bingöl, Turkey. Department of Electronic Communication Technology, Vocational School of Technical Sciences, Bingöl University, 12000 Bingöl, Turkey. Department of Occuptional Health and Safety, Faculty of Health Science, Laboratory of Natural Product Research, Bingöl University, 12000 Bingöl, Turkey. Department of Occuptional Health and Safety, Faculty of Health Science, Laboratory of Natural Product Research, Bingöl University, 12000 Bingöl, Turkey. Department of Radiotherapy, Üsküdar University, 34662 İstanbul, Turkey. Department of Radiotherapy, Üsküdar University, 34662 İstanbul, Turkey. The mass attenuation coefficients (µ/ρ) for different parts (root, flower, stem, and leaf) of three medicinal aromatic plants ( Teucrium chamaedrys L. subsp. sinuatum , Rheum ribes , and Chrysophthalmum montanum ) were measured using an 241 Am photon source in a stable geometry and calculated using the Monte Carlo N-Particle Transport Code System-extended (MCNPX) code and the WinXCOM program. The experimental and theoretical MCNPX and WinXCOM values exhibited good agreement. The measured mass attenuation coefficient values were then used to compute the effective atomic number ( Z eff ) and electron density ( N E )of the samples. The results reveal that S1-S (stem of Teucrium chamaedrys L. subsp. sinuatum ) has the highest values of µ/ρand Z eff . Medicinal aromatic plant MCNPX code Mass attenuation coefficient Photon WinXCOM 1. Introduction 1 Ionizing radiation is energy emitted from radioactive atoms in the form of particles (such as alpha, beta, and neutron) or electromagnetic waves (such as X-rays and gamma rays). Gamma rays and X-rays are chargeless and massless photons; thus, they are the most penetrating forms of ionizing radiation. Humans are exposed to ionizing radiation through different processes, such as diagnostic treatment of various cancers, nuclear power generation, and medical radiation exposure from X-rays. It is well-known that exposure to ionizing radiation leads to the ionization of several biologically important macromolecules, including membranes, proteins, lipids, and nucleic acids [ 1 1 , 2 2 ].Ionizing radiation passing through living tissues results in the generation of reactive O species, which interact with DNA, causing cell damage, dysfunction, or death [ 3 3 ].In recent decades, various studies have been performed to protect humans from the hurtful effects of ionizing radiation and to develop a novel and active radio protector [ 4 4 ].In this regard, numerous compounds have been investigated for their capability as a radio-protective material, such as vitamins, amino acids, plants, glasses and herbs [ 4 4 – 9 9 4–9 ].Radio-protective materials should be highly effective, cheap, orally administered, and without toxicity implications. Natural products such as plants have all of these properties. It is well-known that medicinal plants are utilized in different traditional systems of medicine because they are considered non-toxic and are exceedingly accepted by people. The use of medici nal plants as radio-protective materials requires scientific evaluation and assessment. The precise values of photon interaction parameters, such as the mass attenuation coefficient and effective atomic number, are needed to clarify the interaction of photons with medicinal plants. The mass attenuation coefficient (µ/ρ) for a medicinal plant indicates how easily it can be penetrated by a photon. This coefficient is considered to be one of the most important parameters characterizing the interaction between ionizing radiation and a specific medium and can be used to calculate other parameters, such as the effective atomic number and electron density [ 10 10 , 11 11 ]. There are different methods for calculating the mass attenuation coefficients for a specific photon energy and material, including a theoretical approach employing the WinXCOM program [ 12 12 ], an experimental method using transmission geometry [ 13 13 – 15 15 13–15 ], and simulation codes such as Monte Carlo N-Particle Transport Code System-extended (MCNPX), Geant4, and FLUKA[ 10 10 , 16 16 , 17 17 ]. In the present study, we used the transmission geometry to measure the mass attenuation coefficients of different parts (root, flower, stem, and leaf) of three medicinal aromatic plants ( Teucrium chamaedrys L. subsp. sinuatum , Rheum ribes and Chrysophthalmum montanum ). Additionally, we used the MCNPX code and the WinXCOM program to validate the obtained values. Using the obtained mass attenuation coefficients, we calculated related parameters, such as the effective atomic number and electron density. 2. Materials and methods 1 2.1. Sample preparation 2 Details regarding the localities of the plant samples used in the study are presented in Table1 Table1 .The samples were collected from natural habitats and transferred to the laboratory. The plant samples were divided into target organs while in the wet state. The target plant organs were allowed to be dried on unprinted papers in the laboratory. A laboratory-type mill was used to pulverize the dried plant samples. The powdered plant specimens were passed through a laboratory sieve and brought to a 20-mesh particle size. The milled plant samples were pelleted using a laboratory-type hydraulic pellet press. 1 1 Plant samples 1 1 Localities 1 1 Coordinates 1 1 Herbarium numbers 1 1 Gathering times 1 1 Rheum ribes 1 1 Bingöl-Yelesen village rocky slopes 1750 m 1 1 N: 38°52′09.2″ 1 1 KOCAK-4232 1 1 20.06.2017 1 1   1 1   1 1 E: 40°19′26.5″ 1 1 Chrysophthalmum montanum 1 1 Diyarbakır-Kulp district vicinity, 1425 m 1 1 N: 38°30′23.3″ 1 1 KOCAK-4231 1 1 16.06.2016 1 1   1 1   1 1 E: 41°05′44.6″ 1 1 Teucrium chamaedrys L. subsp. sinuatum 1 1 Bingöl-Sancak district and surrounding slopes 1550 m 1 1 N:39°10′31.2″ 1 1 BIN-782 1 1 08.07.2016 1 1   1 1   1 1 E:40°22′16.8″ 1 1   1 1   The total organic content in the samples was analyzed using a spectrometer (FLASH 2000 Series CHNS / Organic Elemental Analyzer, Thermo Scientific). The total C, H, and N contents in each plant sample were determined according to the European standard EN15104:2011 [ 18 18 ], which suggests the following formula for the amount of O in the samples: <math display="block" id="M1"><mrow><mtext>O </mtext><mrow><mo>(</mo><mtext>%</mtext><mo>)</mo></mrow><mo>=</mo><mn>100</mn><mo>−</mo><msub><mrow><mrow><mo>[</mo><mrow><mtext>C% </mtext><mo>+</mo><mtext>H%</mtext><mo>+</mo><mtext>N%</mtext></mrow><mo>]</mo></mrow></mrow><mrow><mrow><mo>(</mo><mrow><mtext>Lightest Moist Basis</mtext></mrow><mo>)</mo></mrow></mrow></msub><mo>.</mo></mrow></math> The chemical compositions, densities, and codes of the present medicinal aromatic plants are listed in Table 2 Table 2 . 1 1 Sample 1 1 Used part 1 1 Sample code 1 1 Density (g/cm 3 ) 1 1 H (%) 1 1 C (%) 1 1 N (%) 1 1 O (%) 1 1 Teucrium chamaedrys L. subsp. sinuatum 1 1 root 1 1 S1-R 1 1 1.051 1 1 4.963 1 1 92.989 1 1 0.866 1 1 1.182 1 1 Teucrium chamaedrys L. subsp. sinuatum 1 1 flower 1 1 S1-F 1 1 1.152 1 1 4.560 1 1 35.080 1 1 1.821 1 1 58.539 1 1 Teucrium chamaedrys L. subsp. sinuatum 1 1 stem 1 1 S1-S 1 1 1.001 1 1 2.875 1 1 22.235 1 1 0.965 1 1 73.925 1 1 Rheum ribes 1 1 stem 1 1 S2-S 1 1 1.273 1 1 6.702 1 1 40.663 1 1 0.871 1 1 51.764 1 1 Rheum ribes 1 1 root 1 1 S2-R 1 1 1.135 1 1 4.368 1 1 65.964 1 1 0.736 1 1 28.932 1 1 Rheum ribes 1 1 leaf 1 1 S2-L 1 1 1.001 1 1 11.448 1 1 47.435 1 1 2.425 1 1 38.691 1 1 Chrysophthalmum montanum 1 1 stem 1 1 S3-S 1 1 1.367 1 1 5.308 1 1 51.544 1 1 0.674 1 1 42.474 1 1 Chrysophthalmum montanum 1 1 leaf 1 1 S3-L 1 1 1.093 1 1 5.383 1 1 56.705 1 1 1.525 1 1 36.387 1 1 Chrysophthalmum montanum 1 1 flower 1 1 S3-F 1 1 1.042 1 1 5.317 1 1 52.690 1 1 1.519 1 1 40.475 2.2. Experimental work 2 The spectrometer system comprised a Si(Li) semiconductor detector (SLP-04160P-OPT-0.3) having an energy resolution of 2.7% at 5.9 keV, a 12.5-mm 2 active area, and a 5-mm sensitive depth.A4K multichannel analyzer was used to measure the mass attenuation coefficients of the medicinal aromatic plants. An 241 Am photon source having activity of 370 kBq was used to irradiate the plant samples. This source provided the following useful photon energies: 13.92, 17.75, 20.78, 26.34, and 59.54 keV. The radioactive source enclosed in a lead cylinder with a 3-mm orifice was collimated to produce a narrow beam. The experimental narrow-beam geometry setup is shown in Fig. 1 Fig. 1 .The distance between the source and detector was 15 cm. The detector was housed in a lead shield to minimize the detection of background radiation from the surroundings. Each plant sample was placed between the source and detector. Then, the photon intensities I 0 (incident) and I (after attenuation) were measured by recording corresponding counts for a fixed preset time. For each sample, gamma-ray intensities without ( I o ) and after attenuation (I) were used to calculate the mass attenuation coefficients using the Lambert–Beer law[ 17 17 ].More extensive information regarding the experimental process is presented in our previous works [ 19 19 – 21 21 19–21 ]. A typical spectrum of 13.92 keV with and without attenuation by the S2-S, S2-L, and S2-R samples is shown in Fig. 2 Fig. 2 . 2.3. MCNPX code 2 The Monte Carlo method is one of the best ways to resolve physical problems, medical difficulties, and unknown parameters when the experimental or clinical conditions are limited or difficult to attain. MCNPX version 2.6.0 (Los Alamos National Lab, USA) was utilized for the calculations of the mass attenuation coefficients of different parts of three medicinal aromatic plants. The MCNPX code is a general-purpose Monte Carlo radiation transport code for simulating the interaction of radiation with matter in a wide energy range [ 22 22 ]. MCNPX can ensure fully three-dimensional simulation and utilizes extended nuclear cross-section libraries and physics models for particle types. The validation of the MCNPX Monte Carlo code for the investigation of radiation mass attenuation coefficients and other shielding parameters for different types of materials—such as glasses, concrete, and living biological tissues—has been reported [ 23 23 – 33 33 23–33 ]. The computational modeling of radiation interaction problems, including the radiation attenuation and shielding properties,depends on the material definitions. The input structure of the MCNPX code has three major partsdefining the significant details of the simulation. A completed MCNPX input file describes the problem geometry and defines the materials with their chemical compositions, as well as the structure of the radiation source,including details such as the source geometry andradiation distribution. The simulation geometry is constructed by defining cells. Each cell in the simulation is bounded by one or more geometric surfaces. In the present investigation, the square prism geometry was employed for modeling the parts of medicinal aromatic plants. The edge lengths of this square prism geometry were defined as 5 cm, while the axial z-length was set differentlyfor each simulation because it represents the thickness of the parts of the medicinal aromatic plants. The mass attenuation coefficients of each part of the medicinal aromatic plant sample were measured in a narrow beam transmission geometry using a point isotropic source with a collimated and monoenergetic beam. The radiation energy value of the point isotropic source was defined in the photon energy range of 13.92–59.54 keV for each calculation. The geometric center of the detection area on the central axis was considered as the location of the point isotropic source that emits gamma rays perpendicular to the front face of the medicinal aromatic plant sample in the direction of the z axis. Another important definition is the material definition, considering the atomic number, mass number, elemental mass fraction, and density of compounds or mixtures. The elemental mass fractions and densities of each medicinal aromatic plant sample are presented in Table 2 Table 2 . In the present MCNPX simulation, to obtain the absorbed dose amount in the detection field, the average flux tally (F4) was utilized. This type of tally mash gives the sum of the average flux in the cell. The initial quantity of gamma rays is set as 10 8 particles. The mass attenuation coefficient calculations were performed using Intel® Core™ i7 CPU 2.80 GHz computer hardware. The error rate in the output file was observed to be < 0.1%. Analysis of a recent investigation was performed using the D00205ALLCP03 MCNPXDATA package, whichcomprisesDLC-200/MCNPDATA cross-section libraries. This library typically includes ENDF/B-VI data in the range of 20 to 150 MeV. Fig. 3 Fig. 3 shows the total simulation geometry for the present investigation.A screenshot showing the design of the MCNPX simulation setup is displayed in Fig.4 Fig.4 . 3. Results and discussion 1 Different photon attenuation coefficients—such as the mass attenuation coefficients (µ/ρ), effective atomic number ( Z eff ), and electron density ( N e )—of different parts (root, flower, stem, and leaf) of three medicinal aromatic plants ( Teucrium chamaedrys L.subsp. sinuatum , Rheum ribes , and Chrysophthalmum montanum ) were investigated. The µ/ρ values of the present samples were obtained in three ways: (a) experimentally, using a transmission arrangement; (b) theoretically, using the WinXCOM program; and (c) using a Monte Carlo simulation (MCNPX code). The experimental and WinXCOM results for the µ/ρ of the present samples at 13.92, 17.75, 20.78, 26.34, and 59.54 keV are illustrated in Fig.5 Fig.5 , while the experimental and MCNPX results for the µ/ρ of all the samples are shown in Fig.6 Fig.6 . The experimental, WinXCOM, and MCNPX results for the µ/ρ are summarized in Table 3 Table 3 . 1 1 Code 3 1 13.92 keV 3 1 17.75 keV 3 1 20.78 keV 3 1 26.34 keV 3 1 59.54 keV 1 1   1 1 Exp 1 1 WinX 1 1 MCNPX 1 1 Exp 1 1 WinX. 1 1 MCNPX 1 1 Exp 1 1 WinX 1 1 MCNPX 1 1 Exp 1 1 WinX 1 1 MCNPX 1 1 Exp 1 1 WinX 1 1 MCNPX 1 1 S1-R 1 1 0.990±0.092 1 1 0.957 1 1 0.969 1 1 0.526±0.015 1 1 0.555 1 1 0.570 1 1 0.407±0.020 1 1 0.416 1 1 0.423 1 1 0.285±0.031 1 1 0.300 1 1 0.310 1 1 0.178±0.013 1 1 0.184 1 1 0.195 1 1 S1-F 1 1 1.641±0.052 1 1 1.713 1 1 1.722 1 1 0.926±0.062 1 1 0.912 1 1 0.922 1 1 0.654±0.012 1 1 0.636 1 1 0.659 1 1 0.400±0.016 1 1 0.405 1 1 0.413 1 1 0.184±0.008 1 1 0.192 1 1 0.200 1 1 S1-S 1 1 1.997±0.044 1 1 1.914 1 1 1.921 1 1 0.965±0.025 1 1 1.005 1 1 1.118 1 1 0.665±0.007 1 1 0.692 1 1 0.703 1 1 0.409±0.010 1 1 0.431 1 1 0.443 1 1 0.181±0.005 1 1 0.192 1 1 0.201 1 1 S2-S 1 1 1.538±0.041 1 1 1.606 1 1 1.621 1 1 0.854±0.052 1 1 0.863 1 1 0.875 1 1 0.624±0.067 1 1 0.607 1 1 0.627 1 1 0.400±0.014 1 1 0.393 1 1 0.627 1 1 0.184±0.007 1 1 0.194 1 1 0.201 1 1 S2-R 1 1 1.285±0.063 1 1 1.319 1 1 1.326 1 1 0.698±0.061 1 1 0.726 1 1 0.749 1 1 0.502±0.011 1 1 0.521 1 1 0.534 1 1 0.335±0.017 1 1 0.350 1 1 0.534 1 1 0.184±0.008 1 1 0.187 1 1 0.190 1 1 S2-L 1 1 1.351±0.082 1 1 1.390 1 1 1.401 1 1 0.739±0.014 1 1 0.766 1 1 0.776 1 1 0.559±0.017 1 1 0.551 1 1 0.576 1 1 0.354±0.026 1 1 0.371 1 1 0.576 1 1 0.193±0.012 1 1 0.200 1 1 0.210 1 1 S3-S 1 1 1.528±0.155 1 1 1.491 1 1 1.511 1 1 0.837±0.016 1 1 0.808 1 1 0.837 1 1 0.546±0.024 1 1 0.572 1 1 0.595 1 1 0.389±0.032 1 1 0.375 1 1 0.397 1 1 0.194±0.013 1 1 0.191 1 1 0.195 1 1 S3-L 1 1 1.469±0.102 1 1 1.420 1 1 1.433 1 1 0.749±0.017 1 1 0.774 1 1 0.785 1 1 0.526±0.021 1 1 0.551 1 1 0.575 1 1 0.356±0.057 1 1 0.365 1 1 0.378 1 1 0.188±0.014 1 1 0.190 1 1 0.199 1 1 S3-F 1 1 1.394±0.088 1 1 1.468 1 1 1.481 1 1 0.781±0.012 1 1 0.797 1 1 0.803 1 1 0.544±0.023 1 1 0.565 1 1 0.577 1 1 0.383±0.025 1 1 0.372 1 1 0.389 1 1 0.195±0.011 1 1 0.190 1 1 0.206 The tendency of µ/ρ is almost the same for all samples owing to the similarity of the elemental compositions among the samples (all the samples contain H, C, N, and O in different proportions). Figs.5 Figs.5 and 6 6 clearly show that the experimentally measured µ/ρ values of the medicinal aromatic plants agree well with the WinXCOM and MCNPX results. The present experimental WinXCOM and MCNPX results are similar to the observations of Singh et al. [ 25 25 ],who studied the µ/ρ of polymers in the energy range of 59.5–1,332.5 MeV using Monte Carlo simulation (MCNP-4C) and reported good agreement among theirMCNP-4C,XCOM,and experimental results. The µ/ρ values for all samples decrease with the increase of the photon energy. This means that with the increase of the energy of the incident photons, we obtain smaller attenuation; thus, more photons can penetrate the present samples. The differences among the µ/ρ values due to the medicinal aromatic plant compositions are obvious, particularly in the low-energy region. As shown in Table 3 Table 3 , S1-S (stem of Teucrium chamaedrys L. subsp. sinuatum ) has higher values of µ/ρ than the other samples because of the larger amount of O in this sample (73.925%).Additionally, Table 3 Table 3 shows that the S1-R sample has the lowest µ/ρ values, as this sample contains mainly C (92.989%), while the percent of O in this sample is low (1.182%).It is important to investigate the experimental uncertainty in the measurement of µ/ρ.The following equation was used for this: <math display="block" id="M2"><mrow><mi>Δ</mi><mrow><mo>(</mo><mrow><mfrac><mi>μ</mi><mi>ρ</mi></mfrac></mrow><mo>)</mo></mrow><mo>=</mo><mfrac><mn>1</mn><mrow><mi>ρ</mi><mi>x</mi></mrow></mfrac><msqrt><mrow><msup><mrow><mrow><mo>(</mo><mrow><mfrac><mrow><mi>Δ</mi><msub><mi>I</mi><mn>0</mn></msub></mrow><mrow><msub><mi>I</mi><mn>0</mn></msub></mrow></mfrac></mrow><mo>)</mo></mrow></mrow><mn>2</mn></msup><mo>+</mo><msup><mrow><mrow><mo>(</mo><mrow><mfrac><mrow><mi>Δ</mi><mi>I</mi></mrow><mi>I</mi></mfrac></mrow><mo>)</mo></mrow></mrow><mn>2</mn></msup><mo>+</mo><msup><mrow><mrow><mo>(</mo><mrow><mtext>In</mtext><mfrac><mrow><msub><mi>I</mi><mn>0</mn></msub></mrow><mi>I</mi></mfrac></mrow><mo>)</mo></mrow></mrow><mn>2</mn></msup><mo>.</mo><msup><mrow><mrow><mo>(</mo><mrow><mfrac><mrow><mi>Δ</mi><mi>ρ</mi><mi>x</mi></mrow><mrow><mi>ρ</mi><mi>x</mi></mrow></mfrac></mrow><mo>)</mo></mrow></mrow><mn>2</mn></msup></mrow></msqrt><mo>,</mo></mrow></math> where x and ρ represent the thickness and density of the sample, respectively, and Δ I 0 and Δ I are the uncertainties of I 0 and I , respectively. Eq. (2) was used to calculate the experimental uncertainty in the measurement of µ/ρ, and the results are summarized in Table 3 Table 3 . Clearly, the uncertainties in the experimental µ/ρ measurements are within the range of 1.09% to 10.73% (except for S1-R and S3-F at 26.34 keV). The uncertainties in the measurements are mainly due to statistical uncertainties in the determination of the intensity without ( I 0 ) and with ( I ) the absorber, as well as to the scattered photons reaching the detector. By usingthemeasured and calculateddatafor µ/ρ, the effective atomic numbers ( Z eff ) were calculatedas follows [ 34 34 ]: <math display="block" id="M3"><mrow><msub><mi>Z</mi><mrow><mtext>eff</mtext></mrow></msub><mo>=</mo><mfrac><mrow><msub><mstyle mathsize="140%" displaystyle="true"><mo>∑</mo></mstyle><mi>i</mi></msub><msub><mi>f</mi><mi>i</mi></msub><msub><mi>A</mi><mi>i</mi></msub><msub><mrow><mrow><mo>(</mo><mrow><mi>μ</mi><mo>/</mo><mi>ρ</mi></mrow><mo>)</mo></mrow></mrow><mi>i</mi></msub></mrow><mrow><msub><mstyle mathsize="140%" displaystyle="true"><mo>∑</mo></mstyle><mi>j</mi></msub><msub><mi>f</mi><mi>j</mi></msub><mfrac><mrow><msub><mi>A</mi><mi>j</mi></msub></mrow><mrow><msub><mi>Z</mi><mi>j</mi></msub></mrow></mfrac><msub><mrow><mrow><mo>(</mo><mrow><mfrac bevelled="true"><mi>μ</mi><mi>ρ</mi></mfrac></mrow><mo>)</mo></mrow></mrow><mi>j</mi></msub></mrow></mfrac><mo>,</mo></mrow></math> where f i represents the fractional abundance of the element i relative to the number of atoms, and Z i and A i are the atomic number and atomic weight, respectively. The results for the Z eff of the selected samples are plotted with respect to the photon energy in Fig. 7 Fig. 7 .In addition, the experimental Z eff values are presented in Table 4 Table 4 , together with theoretical values. The Z eff results indicate good agreement between the experiment and theory. As shown in Table 4 Table 4 , the values of Z eff for all the medicinal aromatic plants change with respect to the photon energy. In the case of Teucrium chamaedrys L. subsp. sinuatum , Z eff ranges from 5.692 to 4.047, 6.725 to 4.613, and 7.744 to 5.348 for the root, flower, and stem, respectively. In the case of Rheum ribes , Z eff ranges from 6.419 to 3.954, 6.226 to 4.499, and 5.795 to 3.107 for the stem, root, and leaf, respectively. In the case of Chrysophthalmum montanum , Z eff ranges from 6.790 to 4.547, 6.720 to 4.338, and 6.247 to 4.512 for the stem, leaf, and flower, respectively. Notably, the atomic numbers ( Z ) of the elements of the present medicinal aromatic plants change from 1 (H) to 8 (O), and a total of four elements (H, C, N, and O) are considered. The values of Z eff for the present plants lie within the range of Z for the constituent elements (1 < Z eff < 8). As indicated by the µ/ρ results, S1-S (stem of Teucrium chamaedrys L. subsp. sinuatum ) has higher values of Z eff . 1 1 Code 2 1 13.92 keV 2 1 17.75 keV 2 1 20.78 keV 2 1 26.34 keV 2 1 59.54 keV 1 1   1 1 Exp 1 1 Theo 1 1 Exp 1 1 Theo 1 1 Exp 1 1 Theo 1 1 Exp 1 1 Theo 1 1 Exp 1 1 Theo 1 1 S1-R 1 1 5.692±0.527 1 1 5.500 1 1 4.900±0.142 1 1 5.169 1 1 4.824±0.232 1 1 4.938 1 1 4.392±0.484 1 1 4.630 1 1 4.047±0.288 1 1 4.164 1 1 S1-F 1 1 6.725±0.211 1 1 7.022 1 1 6.740±0.452 1 1 6.637 1 1 6.506±0.115 1 1 6.321 1 1 5.734±0.234 1 1 5.816 1 1 4.613±0.194 1 1 4.821 1 1 S1-S 1 1 7.744±0.169 1 1 7.420 1 1 6.880±0.178 1 1 7.167 1 1 6.678±0.073 1 1 6.945 1 1 6.235±0.154 1 1 6.561 1 1 5.348±0.138 1 1 5.678 1 1 S2-S 1 1 6.419±0.172 1 1 6.703 1 1 6.130±0.376 1 1 6.195 1 1 5.961±0.640 1 1 5.801 1 1 5.294±0.186 1 1 5.212 1 1 3.954±0.146 1 1 4.176 1 1 S2-R 1 1 6.226±0.303 1 1 6.391 1 1 5.790±0.502 1 1 6.020 1 1 5.528±0.125 1 1 5.735 1 1 5.092±0.256 1 1 5.311 1 1 4.499±0.192 1 1 4.575 1 1 S2-L 1 1 5.795±0.354 1 1 5.962 1 1 5.067±0.096 1 1 5.252 1 1 4.843±0.150 1 1 4.768 1 1 3.949±0.294 1 1 4.130 1 1 3.107±0.185 1 1 3.205 1 1 S3-S 1 1 6.790±0.688 1 1 6.629 1 1 6.427±0.122 1 1 6.199 1 1 5.601±0.244 1 1 5.864 1 1 5.552±0.454 1 1 5.358 1 1 4.547±0.301 1 1 4.456 1 1 S3-L 1 1 6.720±0.466 1 1 6.499 1 1 5.868±0.130 1 1 6.066 1 1 5.463±0.219 1 1 5.733 1 1 5.112±0.818 1 1 5.239 1 1 4.338±0.329 1 1 4.382 1 1 S3-F 1 1 6.247±0.396 1 1 6.580 1 1 6.025±0.096 1 1 6.144 1 1 5.588±0.236 1 1 5.806 1 1 5.458±0.355 1 1 5.300 1 1 4.512±0.265 1 1 4.409 Z eff is related to another parameter called the effective electron density ( N E ), which represents the number of electrons per unit mass of the interacting materials and can be calculated using the following relation [ 35 35 ]: <math display="block" id="M4"><mrow><msub><mi>N</mi><mtext>E</mtext></msub><mo>=</mo><msub><mi>N</mi><mtext>A</mtext></msub><mfrac><mrow><msub><mi>n</mi><mrow><mtext>tot</mtext></mrow></msub><msub><mi>Z</mi><mrow><mtext>eff</mtext></mrow></msub></mrow><mrow><msub><mstyle mathsize="140%" displaystyle="true"><mo>∑</mo></mstyle><mi>i</mi></msub><msub><mi>n</mi><mi>i</mi></msub><msub><mi>A</mi><mi>i</mi></msub></mrow></mfrac><mo>=</mo><msub><mi>N</mi><mtext>A</mtext></msub><mfrac><mrow><msub><mi>Z</mi><mrow><mtext>eff</mtext></mrow></msub></mrow><mi>A</mi></mfrac><mo>,</mo></mrow></math> where ˂ A ˃ is the mean atomic mass, and N A is the Avogadro constant. The determined N E for the present samples are plotted in Fig.8 Fig.8 and listed in Table 5 Table 5 . As shown in Fig.8 Fig.8 , the experimental and theoretical (WinXCOM) values agree well. It is also clear from Fig.8 Fig.8 that the N E for all samples decreases with increasing energy owing to the dominance of the photoelectric effect in this energy region. 1 1 Code 2 1 13.92 keV 2 1 17.75 keV 2 1 20.78 keV 2 1 26.34 keV 2 1 59.54 keV 1 1   1 1 Exp 1 1 Theo 1 1 Exp 1 1 Theo 1 1 Exp 1 1 Theo 1 1 Exp 1 1 Theo 1 1 Exp 1 1 Theo 1 1 S1-R 1 1 0.945±0.087 1 1 0.913 1 1 0.813±0.024 1 1 0.858 1 1 0.801±0.039 1 1 0.820 1 1 0.729±0.080 1 1 0.769 1 1 0.672±0.048 1 1 0.691 1 1 S1-F 1 1 0.464±0.015 1 1 0.484 1 1 0.465±0.031 1 1 0.458 1 1 0.449±0.008 1 1 0.436 1 1 0.395±0.016 1 1 0.401 1 1 0.318±0.013 1 1 0.333 1 1 S1-S 1 1 0.845±0.018 1 1 0.809 1 1 0.750±0.019 1 1 0.782 1 1 0.728±0.008 1 1 0.758 1 1 0.680±0.017 1 1 0.716 1 1 0.583±0.015 1 1 0.619 1 1 S2-S 1 1 1.102±0.029 1 1 1.150 1 1 1.052±0.064 1 1 1.063 1 1 1.023±0.110 1 1 0.996 1 1 0.909±0.032 1 1 0.895 1 1 0.679±0.025 1 1 0.717 1 1 S2-R 1 1 1.108±0.054 1 1 1.138 1 1 1.031±0.089 1 1 1.072 1 1 0.984±0.022 1 1 1.021 1 1 0.907±0.046 1 1 0.946 1 1 0.801±0.034 1 1 0.815 1 1 S2-L 1 1 0.479±0.029 1 1 0.492 1 1 0.419±0.008 1 1 0.434 1 1 0.400±0.012 1 1 0.394 1 1 0.326±0.024 1 1 0.341 1 1 0.257±0.015 1 1 0.265 1 1 S3-S 1 1 1.383±0.140 1 1 1.350 1 1 1.309±0.025 1 1 1.263 1 1 1.141±0.050 1 1 1.194 1 1 1.131±0.093 1 1 1.091 1 1 0.926±0.061 1 1 0.908 1 1 S3-L 1 1 0.614±0.043 1 1 0.594 1 1 0.536±0.012 1 1 0.555 1 1 0.500±0.020 1 1 0.524 1 1 0.467±0.075 1 1 0.479 1 1 0.397±0.030 1 1 0.401 1 1 S3-F 1 1 0.566±0.036 1 1 0.596 1 1 0.546±0.009 1 1 0.557 1 1 0.506±0.021 1 1 0.526 1 1 0.495±0.032 1 1 0.480 1 1 0.409±0.024 1 1 0.400 4. Conclusion 1 This work aimed to measure the mass attenuation coefficient (µ/ρ) for different parts (root, flower, stem, and leaf) of three medicinal aromatic plants using an 241 Am photon source in a stable geometry. The µ/ρ for all samples under investigation wasalso calculated using the MCNPX code and the WinXCOM program. The µ/ρ values obtained using the three aforementioned methods were similar. Furthermore, the effective atomic numbers ( Z eff ) and the electron densities ( N E ) for the medicinal aromatic plants were calculated. The calculations indicate thatS1-S (stem of Teucrium chamaedrys L. subsp. sinuatum ) has the highestvalues of µ/ρ and Z eff . The numerical methods and simulation tools employed in this study can be very useful for similar future studies, particularly for cases where it is not possible to perform experimental studies. Monte Carlo simulations are a viable alternative method when experimental investigations are not possible. The recent Monte Carlo method can encourage the scientific community to undertake similar studies and to better evaluate the biological structures and their radiation interaction properties. References K.B. Kalpana , N. Devipriya , M. Srinivasan et al ., Evaluating the radioprotective effect of hesperidin in the liver of Swiss albino mice . Eur. J. Pharmacol . 658 , 206 - 212 ( 2011 ). doi: 10.1016/j.ejphar.2011.02.031 http://doi.org/10.1016/j.ejphar.2011.02.031 Evaluating the radioprotective effect of hesperidin in the liver of Swiss albino mice K.B. Kalpana , N. Devipriya , M. Srinivasan et al ., Investigation of the radioprotective efficacy of hesperidin against gamma-radiation induced cellular damage in cultured human peripheral blood lymphocytes . Mutat. Res . 676 , 54 - 61 ( 2009 ). doi: 10.1016/j.mrgentox.2009.03.005 http://doi.org/10.1016/j.mrgentox.2009.03.005 Investigation of the radioprotective efficacy of hesperidin against gamma-radiation induced cellular damage in cultured human peripheral blood lymphocytes W. Lee , W. Kang , N. Kim et al ., Radioprotective effects of a polysaccharide purified from Lactobacillus plantarum-fermented Ishigeokamurae against oxidative stress caused by gamma ray-irradiation in zebrafish in vivo model . J.Funct. Food . 28 , 83 - 89 ( 2017 ). doi: 10.1016/j.jff.2016.11.004 http://doi.org/10.1016/j.jff.2016.11.004 Radioprotective effects of a polysaccharide purified from Lactobacillus plantarum-fermented Ishigeokamurae against oxidative stress caused by gamma ray-irradiation in zebrafish in vivo model S.J. Hosseinimehr , Trends in the development of radioprotective agents . Drug Discov. Today 12 , 794 - 805 ( 2007 ). doi: 10.1016/j.drudis.2007.07.017 http://doi.org/10.1016/j.drudis.2007.07.017 Trends in the development of radioprotective agents A. Pezeshk , The effects of ionizing radiation on DNA: the role of thiols as radioprotectors . Life Sci . 74 , 2423 - 2429 ( 2004 ). doi: 10.1016/j.lfs.2003.09.068 http://doi.org/10.1016/j.lfs.2003.09.068 The effects of ionizing radiation on DNA: the role of thiols as radioprotectors R. Arora , D. Gupta , R. Chawla et al ., Radioprotection by plant products: present status and future prospects . Phytother. Res . 19 , 1 - 22 ( 2005 ). doi: 10.1002/ptr.1605 http://doi.org/10.1002/ptr.1605 Radioprotection by plant products: present status and future prospects M. Kurudirek , Effective atomic numbers, water and tissue equivalence properties of human tissues, tissue equivalents and dosimetric materials for total electron interaction in the energy region 10 keV-1 GeV . Appl. Radiat. Isotopes 94 , 1 - 7 ( 2014 ). doi: 10.1016/j.apradiso.2014.07.002 http://doi.org/10.1016/j.apradiso.2014.07.002 Effective atomic numbers, water and tissue equivalence properties of human tissues, tissue equivalents and dosimetric materials for total electron interaction in the energy region 10 keV-1 GeV M. Kurudirek , Effective atomic numbers of different types of materials for proton interaction in the energy region 1 keV-10 GeV . Nucl. Instrum. Meth. B 336 , 130 - 134 ( 2014 ). doi: 10.1016/j.nimb.2014.07.008 http://doi.org/10.1016/j.nimb.2014.07.008 Effective atomic numbers of different types of materials for proton interaction in the energy region 1 keV-10 GeV M. Kurudirek , T. Onaran , Calculation of effective atomic number and electron density of essential biomolecules for electron, proton, alpha particle and multi-energetic photon interactions . Radiat. Phys. Chem . 112 , 125 - 138 ( 2015 ). doi: 10.1016/j.radphyschem.2015.03.034 http://doi.org/10.1016/j.radphyschem.2015.03.034 Calculation of effective atomic number and electron density of essential biomolecules for electron, proton, alpha particle and multi-energetic photon interactions G. Lakshminarayana , A. Kumar , M.G. Dong et al ., Exploration of gamma radiation shielding features for titanate bismuth borotellurite glasses using relevant software program and Monte Carlo simulation code . J. Non-Cryst. Solids 481 , 65 - 73 ( 2018 ). doi: 10.1016/j.jnoncrysol.2017.10.027 http://doi.org/10.1016/j.jnoncrysol.2017.10.027 Exploration of gamma radiation shielding features for titanate bismuth borotellurite glasses using relevant software program and Monte Carlo simulation code A. Kumar , M.I. Sayyed , Mengge Dong et al ., Effect of PbO on the shielding behavior of ZnO-P2O5 glass system using Monte Carlo simulation . J. Non-Cryst. Solids 48 , 604 - 607 ( 2018 ). doi: 10.1016/j.jnoncrysol.2017.12.001 http://doi.org/10.1016/j.jnoncrysol.2017.12.001 Effect of PbO on the shielding behavior of ZnO-P2O5 glass system using Monte Carlo simulation L. Gerward , N. Guilbert , K.B. Jensen et al ., WinXCom - a program for calculating X-ray attenuation coefficients . Radiat. Phys.Chem . 71 , 653 - 654 ( 2004 ). doi: 10.1016/j.radphyschem.2004.04.040 http://doi.org/10.1016/j.radphyschem.2004.04.040 WinXCom - a program for calculating X-ray attenuation coefficients J. Kaewkhao , J. Laopaiboon , W. Chewpraditkul , Determination of effective atomic numbers and effective electron densities for Cu/Zn alloy . J. Quant. Spectrosc. Ra . 109 , 1260 - 1265 ( 2008 ). doi: 10.1016/j.jqsrt.2007.10.007 http://doi.org/10.1016/j.jqsrt.2007.10.007 Determination of effective atomic numbers and effective electron densities for Cu/Zn alloy S. Seven , I.H. Karahan , Ӧ.F. Bakkaloglu , The measurement of total mass attenuation coeffcients of CoCuNi alloys . J. Quant. Spectrosc. Ra . 83 , 237 - 242 ( 2004 ). doi: 10.1016/S0022-4073(03)00118-3 http://doi.org/10.1016/S0022-4073(03)00118-3 The measurement of total mass attenuation coeffcients of CoCuNi alloys M. Dong , X. Xue , A. Kumar et al ., A novel method of utilization of hot dip galvanizing slag using the heat waste from itself for protection from radiation . J. Hazard. Mater . 344 , 602 - 614 ( 2018 ). doi: 10.1016/j.jhazmat.2017.10.066 http://doi.org/10.1016/j.jhazmat.2017.10.066 A novel method of utilization of hot dip galvanizing slag using the heat waste from itself for protection from radiation B.O. Elbashir , M.G. Dong , M.I. Sayyed et al ., Comparison of Monte Carlo simulation of gamma ray attenuation coefficients of amino acids with XCOM program and experimental data . Results Phys . 9 , 6 - 11 ( 2018 ). doi: 10.1016/j.rinp.2018.01.075 http://doi.org/10.1016/j.rinp.2018.01.075 Comparison of Monte Carlo simulation of gamma ray attenuation coefficients of amino acids with XCOM program and experimental data M.I. Sayyed , M.Y. AlZaatreh , M.G. Dong et al ., A comprehensive study of the energy absorption and exposure buildup factors of different bricks for gamma-rays shielding . Results Phys . 7 , 2528 - 2533 ( 2017 ). doi: 10.1016/j.rinp.2017.07.028 http://doi.org/10.1016/j.rinp.2017.07.028 A comprehensive study of the energy absorption and exposure buildup factors of different bricks for gamma-rays shielding British Standards Institute-BSI (2011) : Solid Biofuels Standards n°EN 14774-1:2009, EN 14775:2009, EN 15148:2009, EN 15104:2011, EN 15289:2011 [online] . Available at http://www.standards.ie http://www.standards.ie (accessed 21 January 2011 ). Solid Biofuels Standards n°EN 14774-1:2009, EN 14775:2009, EN 15148:2009, EN 15104:2011, EN 15289:2011 [online] F. Akman , R. Durak , M.R. Kaçal et al ., Study of absorption parameters around the K edge for selected compounds of Gd .  X-ray Spectrom .  45 ( 2 ), 103 - 110 ( 2016 ). doi: 10.1002/xrs.2676 http://doi.org/10.1002/xrs.2676 Study of absorption parameters around the K edge for selected compounds of Gd F. Akman , R. Durak , M.F. Turhan et al ., Studies on effective atomic numbers, electron densities from mass attenuation coefficients near the K edge in some samarium compounds .  Appl. Radiat. Isotopes   101 , 107 - 113 ( 2015 ). doi: 10.1016/j.apradiso.2015.04.001 http://doi.org/10.1016/j.apradiso.2015.04.001 Studies on effective atomic numbers, electron densities from mass attenuation coefficients near the K edge in some samarium compounds F. Akman , M.R. Kaçal , F. Akman et al ., Determination of effective atomic numbers and electron densities from mass attenuation coefficients for some selected complexes containing lanthanides .  Can. J. Phys .  95 ( 10 ), 1005 - 1011 ( 2017 ). doi: 10.1139/cjp-2016-0811 http://doi.org/10.1139/cjp-2016-0811 Determination of effective atomic numbers and electron densities from mass attenuation coefficients for some selected complexes containing lanthanides RSICC computer code collection . MCNPX user’s manual version 2.4.0. Monte Carlo n-particle transport code system for multiple and high energy applications , 2002 . MCNPX user’s manual version 2.4.0. Monte Carlo n-particle transport code system for multiple and high energy applications H.O. Tekin , MCNP-X Monte Carlo Code application for mass attenuation coefficients of concrete at different energies by modeling 3X3inch NaI(Tl) detector and comparison with XCOM and Monte Carlo Data . Sci. Technol. Nucl. Install . 2016 , ( 2016 ). doi: 10.1155/2016/6547318 http://doi.org/10.1155/2016/6547318 MCNP-X Monte Carlo Code application for mass attenuation coefficients of concrete at different energies by modeling 3X3inch NaI(Tl) detector and comparison with XCOM and Monte Carlo Data I. Akkurt , H.O. Tekin , A. Mesbahi , Calculation of detection efficiency for the gamma detector using MCNPX . Acta Phys. Pol. A . 128 , 332 - 334 ( 2015 ). doi: 10.12693/APhysPolA.128.B-332 http://doi.org/10.12693/APhysPolA.128.B-332 Calculation of detection efficiency for the gamma detector using MCNPX V.P. Singh , S.P. Shirmardi , M.E. Medhat et al ., Determination of mass attenuation coefficient for some polymers using Monte Carlo simulation . Vacuum 119 , 284 - 288 ( 2015 ). doi: 10.1016/j.vacuum.2015.06.006 http://doi.org/10.1016/j.vacuum.2015.06.006 Determination of mass attenuation coefficient for some polymers using Monte Carlo simulation A.M. El-Khayatt , A.M. Ali , V.P. Singh , Photon attenuation coefficients of heavy-metal oxide glasses by MCNP code, XCOM program and experimental data: A comparison study . Nucl. Instrum. Meth. A 735 , 207 - 212 ( 2014 ). doi: 10.1016/j.nima.2013.09.027 http://doi.org/10.1016/j.nima.2013.09.027 Photon attenuation coefficients of heavy-metal oxide glasses by MCNP code, XCOM program and experimental data: A comparison study H.O. Tekin , V.P. Singh , T. Manici , Effects of micro-sized and nano-sized WO3 on mass attenauation coefficients of concrete by using MCNPX code . Appl. Radiat. Isotopes . 121 , 122 - 125 ( 2017 ). doi: 10.1016/j.apradiso.2016.12.040 http://doi.org/10.1016/j.apradiso.2016.12.040 Effects of micro-sized and nano-sized WO3 on mass attenauation coefficients of concrete by using MCNPX code A. Mesbahi , H. Ghiasi , Shielding properties of the ordinary concrete loaded with micro- and nanoparticles against neutron and gamma radiations . Appl. Radiat. Isotopes . 136 , 27 - 31 ( 2018 ). doi: 10.1016/j.apradiso.2018.02.004 http://doi.org/10.1016/j.apradiso.2018.02.004 Shielding properties of the ordinary concrete loaded with micro- and nanoparticles against neutron and gamma radiations H.O. Tekin , V.P. Singh , T. Manici et al ., Validation of MCNPX with experimental results of mass attenuation coefficients for cement, gypsum and mixture . J. Radiat. Prot. Res . 42 , 154 - 157 ( 2017 ). doi: 10.14407/jrpr.2017.42.3.154 http://doi.org/10.14407/jrpr.2017.42.3.154 Validation of MCNPX with experimental results of mass attenuation coefficients for cement, gypsum and mixture M.G. Dong , R. El-Mallawany , M.I. Sayyed et al ., Shielding properties of 80TeO2--5TiO2-(15-x)WO3--xAnOm glasses using WinXCom and MCNP5 code . Radiat. Phys. Chem . 141 , 172 - 178 ( 2017 ). doi: 10.1016/j.radphyschem.2017.12.012 http://doi.org/10.1016/j.radphyschem.2017.12.012 Shielding properties of 80TeO2--5TiO2-(15-x)WO3--xAnOm glasses using WinXCom and MCNP5 code G. Lakshminarayana , S.O. Baki , K.M. Kaky et al ., Investigation of structural, thermal properties and shielding parameters for multicomponent borate glasses for gamma and neutron radiation shielding applications . J. Non. Cryst. Solids . 471 , 222 - 237 ( 2017 ). doi: 10.1016/j.jnoncrysol.2017.06.001 http://doi.org/10.1016/j.jnoncrysol.2017.06.001 Investigation of structural, thermal properties and shielding parameters for multicomponent borate glasses for gamma and neutron radiation shielding applications H.O. Tekin , M.I. Sayyed , E.E. Altunsoy et al ., Shielding properties and effects of WO3 and PbO on mass attenuation coefficients by using MCNPX code . Dig. J. Nanomater. Biostruct . 12 , 861 - 867 ( 2017 ). Shielding properties and effects of WO3 and PbO on mass attenuation coefficients by using MCNPX code H.O. Tekin , M.I. Sayyed , T. Manici et al ., Photon shielding characterizations of bismuth modified borate - silicate-tellurite glasses using MCNPX Monte Carlo code . Mater. Chem. Phys . 211 , 9 - 16 ( 2018 ). doi: 10.1016/j.matchemphys.2018.02.009 http://doi.org/10.1016/j.matchemphys.2018.02.009 Photon shielding characterizations of bismuth modified borate - silicate-tellurite glasses using MCNPX Monte Carlo code M.I. Sayyed , H. Elhouichet , Variation of energy absorption and exposure buildup factors with incident photon energy and penetration depth for boro-tellurite (B2O3-TeO2) glasses . Radiat. Phys. Chem . 130 , 335 - 342 ( 2017 ). doi: 10.1016/j.radphyschem.2016.09.019 http://doi.org/10.1016/j.radphyschem.2016.09.019 Variation of energy absorption and exposure buildup factors with incident photon energy and penetration depth for boro-tellurite (B2O3-TeO2) glasses M. Kurudirek , Effective atomic numbers and electron densities of some human tissues and dosimetric materials for mean energies of various radiation sources relevant to radiotherapy and medical applications . Radiat. Phys. Chem . 102 , 139 - 146 ( 2014 ). doi: 10.1016/j.radphyschem.2014.04.033 http://doi.org/10.1016/j.radphyschem.2014.04.033 Effective atomic numbers and electron densities of some human tissues and dosimetric materials for mean energies of various radiation sources relevant to radiotherapy and medical applications 10.1007/s41365-018-0475-0 2018-03-17 2018-04-14 2018-05-07 2018-08-30 2018-10-01 2021-04-11 © China Science Publishing & Media Ltd. (Science Press), Shanghai Institute of Applied Physics, the Chinese Academy of Sciences, ChineseNuclear Society and Springer Nature Singapore Pte Ltd. 2018 This is a post-peer-review, pre-copyedit version of an article published in Nuclear Science and Techniques. The final authenticated version is available online at: http://dx.doi.org/10.1007/s41365-018-0475-0 http://dx.doi.org/10.1007/s41365-018-0475-0 . 2018 Nuclear Science and Techniques 10 29