A comprehensive study for mass attenuation coefficients of different parts of the human body through Monte Carlo methods

A comprehensive study for mass attenuation coefficients of different parts of the human body through Monte Carlo methods 1 corresp first-author Corresponding Author: elermis@hotmail.com Corresponding Author: elermis@hotmail.com elermis@hotmail.com 2 2 1 Ege University, Faculty of Science, Physics Department, 35100, Bornova, Izmir, TURKEY Ege University, Faculty of Science, Physics Department, 35100, Bornova, Izmir, TURKEY Uludag University, Faculty of Arts and Science, Physics Department, 16059, Gorukle, Bursa, TURKEY Uludag University, Faculty of Arts and Science, Physics Department, 16059, Gorukle, Bursa, TURKEY The gamma-ray mass attenuation coefficients of blood, bone, lung, eye lenses, adipose, tissue, muscle, brain, and skin were calculated at different energies (60, 80, 150, 400, 500, 600, 1000, 1250, 1500, and 2000 keV) by various theoretical methods,such as FLUKA, GEANT4 Monte Carlo (MC), and the XCOM program in this work. Calculated coefficients were also compared with the National Institute of Standards and Technology (NIST) values. Obtained results were highly in accordance with each other and NIST values.Our results showed that FLUKA was quite convenient in comparison to GEANT4 in the calculation of the mass attenuation coefficients of the used human body samples for low-energy photons (60, 80, and 150 keV) when compared to the NIST values. Gamma-ray mass attenuation coefficient FLUKA GEANT4 XCOM. 1 Introduction 1 When a beam of X or gamma-rays pass through the matter, the removal of photons from the beam is called attenuation. Absorption and scattering of the primary photons cause to attenuation. The linear attenuation coefficient ( µ ) is defined as the fraction of photons removed from a monoenergetic beam of X- or gamma-rays per unit thickness of material. It is typically expressed in units of inverse centimeters (cm -1 ). For a monoenergetic beam of photons,the relationship between the number of incident photons ( N 0 ) and those that are transmitted ( N ) through a thickness x without interaction are exponential. <math display="block" id="M1"><mrow><mi>N</mi><mo>=</mo><msub><mi>N</mi><mn>0</mn></msub><msup><mi>e</mi><mrow><mo>−</mo><mo>µ</mo><mi>x</mi></mrow></msup></mrow></math> The numbers of atoms per volume affect the probability of interaction for a given thickness. The linear attenuation coefficient is normalized by dividing the density of the material to overcome the dependency of the material. This is calledthe mass attenuation coefficient. The units of the mass attenuation coefficient are cm 2 /g [ 1 1 ] . XCOM is a program which generates mass attenuation coefficients for desired energies from 1 keV up to 100 GeV and also elements, compounds, and mixtures. It does not only generate attenuation coefficients but also partial cross-sections for incoherent, coherent scattering, photoelectric absorption, and pair production from the nucleus of atom or atomic electrons [ 2 2 ] . The simulation to estimate the mass attenuation coefficient can also be done with the use of well-known simulation programs such as GEANT4 and FLUKA [ 3 3 - 6 6 3-6 ] .These tools are based on the Monte Carlo (MC) methods to simulate the interaction of particles with their traversing medium. Their application areas vary from a wide range of topics including space, accelerators, medical, high energy, and particle physics [ 7 7 - 10 10 7-10 ] . Tomal et al. [ 11 11 ] experimentally determined the linear attenuation coefficient of the breast tissue. Linear attenuation coefficients of tissue were theoretically calculated by Böke [ 12 12 ] .Akar et al. [ 13 13 ] investigated mass attenuation coefficients of bone, muscle, fat, and water by an experimental method. In medical imaging, such as PET, an attenuation correction for gamma photons is fulfilled to enhance the spatial resolution i.e. image quality. This correction is performed by taking into consideration the mass attenuation coefficients of the related parts of the human body. In this regard, mass attenuation coefficients of the related biological materials have great importance in this process. For this reason, another MC method (FLUKA) was utilized for the first time to determine the coefficients for blood, bone, lung, eye lenses, adipose, tissue, muscle, brain, and skin materials at different energiesin this study as an alternative method,and the obtained results were compared with GEANT4, XCOM, and NIST values. 2 Methods 1 GEANT4 is a C++ based MC simulation code. There are three mandatory classes for the geometry description (DetectorConstruction), physics (PhysicsList), and generated particles (PrimaryGeneratorAction), while there are additional user classes to get information from each step, event, and run (SteppingAction.cc, EventAction.cc and RunAction.cc).FLUKA is another well-known code which is based on the FORTRAN language. Even if there are built-in scoring cards to evaluate requested quantities, it also has various routines to retrieve information from different processes. In this work,userini.f, userin.f, and mgdraw.f routines and theirentries (bxdraw, endraw)were used. For both MC codes, we utilized the advantage of built-in physics lists, PENELOPE (PENetration and Energy LOss of Positrons and Electrons) physics for GEANT4 (version tag 4.10.p02) and PRECISIOn physicsfor FLUKA (version tag 2011.2c). It takes into account detailed photoelectric edge treatments, Compton scattering with inelastic factors and computed without fully taking into account binding and orbital electron motion, Rayleigh scattering, and fluorescence.PENELOPE is essential for low energies, and it includes Compton scattering, photoelectric effects, Rayleigh scattering, gamma conversion, bremsstrahlung, ionization, and positron annihilation [ 14 14 ] . More details about this physics lists can be found in the literature [ 15 15 , 16 16 ] .While performing simulations for each photon energy impinging on different materials, we also took advantage of parallel job executions by using the power of PYTHON scripting language, as well as in the analysis phase of the results. The primary photons impinging on the materials are monoenergetic photons, which are defined as point like particles without any divergence. When the photon transverses in the material, it loses its energy bywell-known processessuch as Compton scattering, photoelectric effects, and pair production.Blood, bone, lung, eye lenses, adipose, tissue, brain, muscle, and skin materials,which have the dimensions of 10 cm (width) x 10 cm (height) x 5 cm (thickness), as sketched in Fig. 1 Fig. 1 ,were selected to investigate the photon attenuation. All material definitions,such as stoichiometry and density,were kept the same in order to compare results amongst performed simulations. The elemental concentrations (% weight) of the materials are presented in Table 1 Table 1 .These concentrations were determined according tothe ICRP report [ 17 17 ] .Photons with different energies were trackedfrom the surface of the material in the simulation. Afterwards, we determined the number of absorbed photons by subtracting the transferred photons from the incident onesin the material to evaluate mass attenuation coefficients from the known Eq.(1). 1 1 Elemental concentration (% weight) 1 1 Materials 1 1   1 1   1 1   1 1   1 1   1 1   1 1   1 1   1 1   1 1 Adipose tissue 1 1 Blood 1 1 Eye lens 1 1 Lung 1 1 Muscle 1 1 Tissue 1 1 Bone 1 1 Brain 1 1 Skin 1 1 H 1 1 0.11947 1 1 0.10186 1 1 0.09926 1 1 0.10127 1 1 0.10063 1 1 0.10447 1 1 0.06398 1 1 0.11066 1 1 0.10058 1 1 C 1 1 0.63724 1 1 0.10002 1 1 0.19371 1 1 0.10231 1 1 0.10783 1 1 0.23219 1 1 0.27800 1 1 0.12542 1 1 0.22825 1 1 N 1 1 0.00797 1 1 0.02964 1 1 0.05327 1 1 0.02865 1 1 0.02768 1 1 0.02488 1 1 0.02700 1 1 0.01328 1 1 0.04642 1 1 O 1 1 0.23233 1 1 0.75941 1 1 0.65375 1 1 0.75707 1 1 0.75477 1 1 0.63023 1 1 0.41001 1 1 0.73772 1 1 0.61900 1 1 Na 1 1 0.00050 1 1 0.00185 1 1 - 1 1 0.00184 1 1 0.00075 1 1 0.00113 1 1 - 1 1 0.00184 1 1 0.00007 1 1 Mg 1 1 0.00002 1 1 0.00004 1 1 - 1 1 0.00073 1 1 0.00019 1 1 0.00013 1 1 0.00200 1 1 0.00015 1 1 0.00006 1 1 Si 1 1 - 1 1 0.00003 1 1 - 1 1 - 1 1 - 1 1 - 1 1 - 1 1 - 1 1 - 1 1 P 1 1 - 1 1 0.00035 1 1 - 1 1 0.00080 1 1 0.00180 1 1 0.00133 1 1 0.07000 1 1 0.00354 1 1 0.00033 1 1 S 1 1 0.00016 1 1 0.00185 1 1 - 1 1 0.00225 1 1 0.00241 1 1 0.00199 1 1 0.00200 1 1 0.00177 1 1 0.00159 1 1 Cl 1 1 0.00119 1 1 0.00278 1 1 - 1 1 0.00266 1 1 - 1 1 0.00134 1 1 - 1 1 0.00236 1 1 0.00267 1 1 K 1 1 - 1 1 0.00163 1 1 - 1 1 0.00194 1 1 - 1 1 0.00199 1 1 - 1 1 0.00310 1 1 0.00085 1 1 Ca 1 1 - 1 1 0.00006 1 1 - 1 1 0.00009 1 1 - 1 1 0.00023 1 1 0.14700 1 1 0.00009 1 1 0.00015 1 1 Fe 1 1 - 1 1 0.00045 1 1 - 1 1 0.00037 1 1 - 1 1 0.00005 1 1 - 1 1 0.00005 1 1 0.00001 1 1 Zn 1 1  - 1 1  0.00001 1 1  - 1 1  0.00001 1 1  - 1 1  0.00003 1 1  - 1 1  0.00001 1 1  0.00001 XCOM program (ver. 3.1) was also used to calculate the gamma-ray mass attenuation coefficients of the selected materials. This program used ICRU Report 44 [ 18 18 ] for material concentration. In the program, material types were first defined by their elemental fractions, which are exactly the same as in FLUKA and GEANT4, and then the gamma-ray energies were specified. The coefficients of the selected materials were finally calculated by the XCOM program. In the MC calculation process, 10 6 photons were sent to the used samples. They were tracked and the transmitted photons through the samples were determined. The linear attenuation coefficients were calculated according to Eq. (1) for each sample. Their mass attenuation coefficients were determined by dividing the obtained linear attenuation coefficients by the sample densities. This process was repeated for 10 3 cycles for each sample. The averages of the results were found and the uncertainties were obtained by calculating their standard deviations. Allcalculated coefficients from FLUKA, GEANT4, and XCOM programs were compared to NIST values [ 19 19 ] . The results are given in the next section. 3 Results 1 The calculated mass attenuation coefficients (µ ρ ) through FLUKA, GEANT4, XCOM, and NIST values(except for skin) were presented in Tables 2 Tables 2 - 10 10 . In these tables,the calculated coefficients viathe FLUKA and GEANT4 programs were given with their standard deviations. 1 1 µ ρ values according to photon energies (keV) 1 1 Method 1 1   1 1   1 1   1 1   1 1  FLUKA 1 1  GEANT4 1 1  XCOM 1 1  NIST 1 1 60 1 1 0.19806±0.00019 1 1 0.18615±0.00039 1 1 0.19740 1 1 0.19740 1 1 80 1 1 0.17985±0.00001 1 1 0.17707±0.00002 1 1 0.18050 1 1 0.18000 1 1 150 1 1 0.15058±0.00010 1 1 0.14824±0.00001 1 1 0.15060 1 1 0.15000 1 1 400 1 1 0.10653±0.00001 1 1 0.10641±0.00008 1 1 0.10670 1 1 0.10620 1 1 500 1 1 0.09881±0.00007 1 1 0.09742±0.00013 1 1 0.09740 1 1 0.09696 1 1 600 1 1 0.08951±0.00004 1 1 0.08951±0.00001 1 1 0.09009 1 1 0.08965 1 1 1,000 1 1 0.07060±0.00013 1 1 0.07279±0.00005 1 1 0.07113 1 1 0.07078 1 1 1,250 1 1 0.06332±0.00006 1 1 0.06329±0.00001 1 1 0.06361 1 1 0.06330 1 1 1,500 1 1 0.05778±0.00008 1 1 0.05763±0.00012 1 1 0.05789 1 1 0.05760 1 1 2,000 1 1 0.05016±0.00001 1 1 0.04979±0.00014 1 1 0.04964 1 1 0.04940 1 1 µ ρ values according to photon energies (keV) 4 1 Methods 1 1   1 1 FLUKA 1 1 GEANT4 1 1 XCOM 1 1 NIST 1 1 60 1 1 0.20108±0.00003 1 1 0.18962±0.00011 1 1 0.20190 1 1 - 1 1 80 1 1 0.18023±0.00012 1 1 0.17506±0.00002 1 1 0.18090 1 1 - 1 1 150 1 1 0.14921±0.00017 1 1 0.14590±0.00019 1 1 0.14880 1 1 - 1 1 400 1 1 0.10489±0.00003 1 1 0.10445±0.00023 1 1 0.10500 1 1 - 1 1 500 1 1 0.09693±0.00010 1 1 0.09548±0.00002 1 1 0.09586 1 1 - 1 1 600 1 1 0.08669±0.00007 1 1 0.08815±0.00004 1 1 0.08863 1 1 - 1 1 1,000 1 1 0.07069±0.00005 1 1 0.06908±0.00001 1 1 0.06998 1 1 - 1 1 1,250 1 1 0.06252±0.00004 1 1 0.06361±0.00007 1 1 0.06257 1 1 - 1 1 1,500 1 1 0.05720±0.00001 1 1 0.05682±0.00002 1 1 0.05694 1 1 - 1 1 2,000 1 1 0.04868±0.00007 1 1 0.04870±0.00005 1 1 0.04888 1 1 - In Figs. 2 Figs. 2 , 3 3 , 4 4 , 5 5 , 6 6 , 7 7 , 8 8 , 9 9 and 10 10 , calculated mass attenuation coefficients versus the photon energies of each absorber material are shown.The comparison of the results from the programs with the NIST values can be seen in these figures and tables as follows. 4 Discussion and Conclusion 1 The gamma-ray mass attenuation coefficients of blood, bone, lung, eye lenses, adipose, tissue, muscle, brain, and skin were calculated at different energies from low energy to high energy (60, 80, 150, 400, 500, 600, 1000, 1250, 1500 and 2000 keV) through FLUKA, GEANT4 MC, and XCOM programs. Calculated mass attenuation coefficients were also compared with NIST values. In the low energy region, the calculated mass attenuation coefficients via GEANT4 were less than those of FLUKA when they were compared to NIST values, as seen in Figs. 2 Figs. 2 , 3 3 , 4 4 , 5 5 , 6 6 , 7 7 , 8 8 , 9 9 and 10 10 and Tables 2 Tables 2 , 3 3 , 4 4 , 5 5 , 6 6 , 7 7 , 8 8 , 9 9 and 10 10 . The results show that the FLUKA and GEANT4 values have differences at low energies up to 8% under the same conditions. But the difference decreases about 1%, which is reasonable for the higher energies. The difference at low energies comes from the variety of physics models used in each MC codes. Fundamentally, both MC codes use the EPDL (Evaluated Photon Data Library) library, which is mostly related to the NIST standard reference data products [ 20 20 ] . FLUKA uses the EPDL97 library for photon cross-section values, with the exception of Compton scattering, without fully considering binding and orbital electron motion [ 21 21 ] and performs detailed photoelectric edge treatments. GEANT4 also uses EPDL97 libraries for photon cross sections for various processes like the photoelectric effect, Compton scattering, the Rayleigh effect, and Bremsstrahlung; sub-shell integrated cross sections for the photoelectric effect and ionization; energy spectra of the secondary particles for electron processes; scattering functions for the Compton effect; form factors for the Rayleigh effect [ 22 22 ] . One reason for the difference between the FLUKA and GEANT4 MC codes at low energies comes from the treatment in this energy regime at which FLUKA handles a detailed photoelectric edge evaluation. On the other hand, the comparison of the MC results with the NIST values shows that calculated mass attenuation coefficients through FLUKA are closer to the NIST values than GEANT4 in the low energy region. It can be concluded that both MC codes are generally good candidates within small relative errors for medical applications. 1 1 µ ρ values according to photon energies (keV) 4 1 Method 1 1   1 1 FLUKA 1 1 GEANT4 1 1 XCOM 1 1 NIST 1 1 60 1 1 0.20330±0.00020 1 1 0.19302±0.00013 1 1 0.20500 1 1 0.20570 1 1 80 1 1 0.18345±0.00009 1 1 0.17421±0.00012 1 1 0.18240 1 1 0.18270 1 1 150 1 1 0.14866±0.00001 1 1 0.14603±0.00020 1 1 0.14920 1 1 0.14920 1 1 400 1 1 0.10593±0.00018 1 1 0.10402±0.00015 1 1 0.10520 1 1 0.10520 1 1 500 1 1 0.09453±0.00002 1 1 0.09592±0.00006 1 1 0.09598 1 1 0.09598 1 1 600 1 1 0.08956±0.00001 1 1 0.08845±0.00016 1 1 0.08873 1 1 0.08874 1 1 1,000 1 1 0.07088±0.00006 1 1 0.07050±0.00003 1 1 0.07006 1 1 0.07007 1 1 1,250 1 1 0.06252±0.00001 1 1 0.06226±0.00003 1 1 0.06265 1 1 0.06265 1 1 1,500 1 1 0.05751±0.00008 1 1 0.05678±0.00001 1 1 0.05701 1 1 0.05701 1 1 2,000 1 1 0.04870±0.00002 1 1 0.04876±0.00003 1 1 0.04896 1 1 0.04896 1 1 µ ρ values according to photon energies (keV) 4 1 Methods 1 1   1 1 FLUKA 1 1 GEANT4 1 1 XCOM 1 1 NIST 1 1 60 1 1 0.20135±0.00008 1 1 0.18781±0.00003 1 1 0.20070 1 1 0.20130 1 1 80 1 1 0.18078±0.00001 1 1 0.17071±0.00013 1 1 0.18030 1 1 0.18030 1 1 150 1 1 0.14807±0.00009 1 1 0.14558±0.00021 1 1 0.14860 1 1 0.14820 1 1 400 1 1 0.10482±0.00005 1 1 0.10367±0.00008 1 1 0.10490 1 1 0.10460 1 1 500 1 1 0.09683±0.00002 1 1 0.09603±0.00006 1 1 0.09576 1 1 0.09547 1 1 600 1 1 0.08854±0.00028 1 1 0.09041±0.00003 1 1 0.08853 1 1 0.08827 1 1 1,000 1 1 0.06961±0.00006 1 1 0.07035±0.00013 1 1 0.06991 1 1 0.06969 1 1 1,250 1 1 0.06138±0.00013 1 1 0.06274±0.00001 1 1 0.06251 1 1 0.06232 1 1 1,500 1 1 0.05615±0.00002 1 1 0.05699±0.00006 1 1 0.05688 1 1 0.05671 1 1 2,000 1 1 0.04869±0.00001 1 1 0.04795±0.00007 1 1 0.04884 1 1 0.04869 1 1 µ ρ values according to photon energies (keV) 4 1 Methods 1 1   1 1 FLUKA 1 1 GEANT4 1 1 XCOM 1 1 NIST 1 1 60 1 1 0.20323±0.00011 1 1 0.18919±0.00006 1 1 0.20520 1 1 0.20530 1 1 80 1 1 0.18380±0.00004 1 1 0.17419±0.00008 1 1 0.18240 1 1 0.18260 1 1 150 1 1 0.14833±0.00008 1 1 0.14702±0.00003 1 1 0.14910 1 1 0.14930 1 1 400 1 1 0.10461±0.00003 1 1 0.10382±0.00021 1 1 0.10510 1 1 0.10530 1 1 500 1 1 0.09585±0.00006 1 1 0.09596±0.00008 1 1 0.09592 1 1 0.09607 1 1 600 1 1 0.08911±0.00004 1 1 0.08915±0.00003 1 1 0.08869 1 1 0.08882 1 1 1,000 1 1 0.07071±0.00015 1 1 0.07006±0.00011 1 1 0.07002 1 1 0.07013 1 1 1,250 1 1 0.06182±0.00001 1 1 0.06199±0.00015 1 1 0.06262 1 1 0.06271 1 1 1,500 1 1 0.05684±0.00002 1 1 0.05649±0.00003 1 1 0.05698 1 1 0.05706 1 1 2,000 1 1 0.04909±0.00003 1 1 0.04831±0.00010 1 1 0.04893 1 1 0.04900 1 1 µ ρ values according to photon energies (keV) 4 1 Methods 1 1   1 1 FLUKA 1 1 GEANT4 1 1 XCOM 1 1 NIST 1 1 60 1 1 0.20207±0.00007 1 1 0.19035±0.00017 1 1 0.20330 1 1 0.20480 1 1 80 1 1 0,18323±0.00014 1 1 0.17421±0.00020 1 1 0.18160 1 1 0.18230 1 1 150 1 1 0.14823±0.00009 1 1 0.14527±0.00010 1 1 0.14900 1 1 0.14920 1 1 400 1 1 0.10550±0.00002 1 1 0.10434±0.00003 1 1 0.10510 1 1 0.10520 1 1 500 1 1 0.09622±0.00011 1 1 0.09588±0.00011 1 1 0.09592 1 1 0.09598 1 1 600 1 1 0.08842±0.00003 1 1 0.08939±0.00008 1 1 0.08868 1 1 0.08874 1 1 1,000 1 1 0.07003±0.00007 1 1 0.06998±0.00001 1 1 0.07002 1 1 0.07007 1 1 1,250 1 1 0.06351±0.00002 1 1 0.06180±0.00005 1 1 0.06261 1 1 0.06265 1 1 1,500 1 1 0.05635±0.00013 1 1 0.05668±0.00003 1 1 0.05698 1 1 0.05701 1 1 2,000 1 1 0.04780±0.00003 1 1 0.04710±0.00003 1 1 0.04893 1 1 0.04896 1 1 µ ρ values according to photon energies (keV) 4 1 Methods 1 1   1 1 FLUKA 1 1 GEANT4 1 1 XCOM 1 1 NIST 1 1 60 1 1 0.20167±0.00007 1 1 0.19044±0.00017 1 1 0.20330 1 1 0.20480 1 1 80 1 1 0.18317±0.00014 1 1 0.17369±0.00020 1 1 0.18170 1 1 0.18230 1 1 150 1 1 0.14863±0.00009 1 1 0.14492±0.00010 1 1 0.14930 1 1 0.14920 1 1 400 1 1 0.10446±0.00002 1 1 0.10591±0.00003 1 1 0.10540 1 1 0.10520 1 1 500 1 1 0.09524±0.00011 1 1 0.09592±0.00011 1 1 0.09619 1 1 0.09598 1 1 600 1 1 0.08882±0.00003 1 1 0.08923±0.00008 1 1 0.08893 1 1 0.08873 1 1 1,000 1 1 0.06922±0.00007 1 1 0.07041±0.00001 1 1 0.07022 1 1 0.07006 1 1 1,250 1 1 0.06294±0.00002 1 1 0.06226±0.00005 1 1 0.06279 1 1 0.06265 1 1 1,500 1 1 0.05801±0.00013 1 1 0.05674±0.00003 1 1 0.05714 1 1 0.05701 1 1 2,000 1 1 0.04866±0.00003 1 1 0.04915±0.00003 1 1 0.04905 1 1 0.04895 1 1 µ ρ values according to photon energies (keV) 4 1 Methods 1 1   1 1 FLUKA 1 1 GEANT4 1 1 XCOM 1 1 NIST 1 1 60 1 1 0.27581±0.00005 1 1 0.25330±0.00002 1 1 0.27520 1 1 0.31480 1 1 80 1 1 0.20855±0.00004 1 1 0.19483±0.00001 1 1 0.20870 1 1 0.22290 1 1 150 1 1 0.14942±0.00002 1 1 0.14325±0.00001 1 1 0.14910 1 1 0.14800 1 1 400 1 1 0.10136±0.00004 1 1 0.09948±0.00004 1 1 0.10180 1 1 0.09910 1 1 500 1 1 0.09263±0.00010 1 1 0.09333±0.00003 1 1 0.09275 1 1 0.09022 1 1 600 1 1 0.08547±0.00002 1 1 0.08548±0.00011 1 1 0.0856 1 1 0.08332 1 1 1,000 1 1 0.06738±0.00002 1 1 0.06759±0.00004 1 1 0.06758 1 1 0.06566 1 1 1,250 1 1 0.06061±0.00003 1 1 0.06013±0.00007 1 1 0.06043 1 1 0.05871 1 1 1,500 1 1 0.05478±0.00013 1 1 0.05556±0.00003 1 1 0.05501 1 1 0.05346 1 1 2,000 1 1 0.04710±0.00004 1 1 0.04737±0.00006 1 1 0.04733 1 1 0.04607 1 1 µ ρ values according to photon energies (keV) 4 1 Methods 1 1   1 1 FLUKA 1 1 GEANT4 1 1 XCOM 1 1 NIST 1 1 60 1 1 0.20274±0.00005 1 1 0.19398±0.00007 1 1 0.20650 1 1 0.20580 1 1 80 1 1 0.18351±0.00002 1 1 0.17661±0.00033 1 1 0.18380 1 1 0.18310 1 1 150 1 1 0.14752±0.00003 1 1 0.14690±0.00023 1 1 0.15030 1 1 0.14980 1 1 400 1 1 0.10443±0.00002 1 1 0.10563±0.00009 1 1 0.10600 1 1 0.10560 1 1 500 1 1 0.09456±0.00003 1 1 0.09716±0.00009 1 1 0.09672 1 1 0.09640 1 1 600 1 1 0.08729±0.00001 1 1 0.09047±0.00001 1 1 0.08942 1 1 0.08913 1 1 1,000 1 1 0.06932±0.00005 1 1 0.07082±0.00004 1 1 0.07061 1 1 0.07037 1 1 1,250 1 1 0.06163±0.00007 1 1 0.06321±0.00008 1 1 0.06314 1 1 0.06293 1 1 1,500 1 1 0.05671±0.00001 1 1 0.05790±0.00004 1 1 0.05745 1 1 0.05726 1 1 2,000 1 1 0.04799±0.00010 1 1 0.04893±0.00005 1 1 0.04933 1 1 0.04917 Our results showed thatFLUKA was somewhat successful in comparison to GEANT4 in the calculation of the mass attenuation coefficients of the human body samples used here for low-energy photons (60, 80, and 150 keV) when compared to the NIST values.Additionally, Robert et al. [ 23 23 ] expressed that GEANT4 was originally developed for high energy physics applications, whereas FLUKA was developed and successfully applied both in the high and the low energy ranges. Our results have supported this expression in terms of gamma photon attenuations in the human body samples. Medhat [ 24 24 ] has calculated mass attenuation coefficients of three biological samples (blood, bone, and muscle) by GEANT4. When we compared our GEANT4 results with his results, it was seen that our results were in agreement with them, even if the used energies and material concentrations were different from each other. Further more, comparison showed that Medhat’s coefficient for muscle at 81 keV increased (0.220 cm 2 g -1 ) unexpectedly while the coefficient was 0.200 cm 2 g -1 for 59.5 keV gamma energy. However, our coefficient for muscle at 80 keV has normally decreased as the energy was increased.This showed that our result for muscle at 80 keVwas more acceptable than his result. As the incident photon energy increases, a decrease in the attenuation coefficient of a material is expected. It is suggested from our results that similar calculations for different human body or biological samples should be repeated and the obtained results should be supported by experimental data for an accurate and satisfactory conclusion. Our analysis serves as a starting point for better understanding the discrepancies between the mass attenuation coefficient results of the GEANT4 and FLUKA programs for the human body samples, especially at low energies. The XCOM method was also used to calculate the gamma-ray mass attenuation coefficients of the materials. The mass attenuation coefficients of brain from XCOM were foundto be higher than those of other programs,especially at low energies ( Table 9 Table 9 , Fig. 9 Fig. 9 ), although its elemental concentration (% weight) was the same as FLUKA and GEANT4.Since the NIST value and an experimental result were not available for the brain, it was impossible to compareour resultswith reference data. But this conclusion needs to be checked and repeated with further studies. Different theoretical methods, such as Monte Carlo N-Particle Transport Code (MCNP), GEANT4,and XCOM have been used to calculate the gamma-ray mass attenuation coefficients of various elementsand compounds at various energies [ 25 25 - 27 27 25-27 ] . The FLUKA program was additionally used to determine the coefficients here, unlike these studies, as a novelty. The calculations through FLUKA, GEANT4, and XCOM, in this work, were performed for more materials and energy values than these studies also. Consequently, the study was carried out to determine individual absorption radiation dose and interpret the criteria for radiation damage of the body in the future studies.It has been noticed that FLUKA results agreed betterin the selected energy range of gammas with NIST values than the others.It can be proposed from this conclusion that FLUKA MC program can be efficientlyused in the determination of gamma-ray mass attenuation coefficients of the samples of the human body, as well as GEANT4 and XCOM. In addition, it is believed that the FLUKA MC program can be a useful alternative tool for medical physics applications. References Bushberg J T , Seibert J A , Leidholdt E M , Boone J M . The essential physics of medical imaging . New York (USA) : William-Wilkins Press , 2001 . The essential physics of medical imaging Berger MJ , Hubbell JH . Photon cross section on a personal computer (XCOM). NBSIR87-3597 , 1987 . Photon cross section on a personal computer (XCOM). NBSIR87-3597 Böhlen T T , Cerutti F , Chin M P W , et al . The FLUKA code: developments and challenges for high energy and medical applications . Nucl Data Sheets , 2014 , 120 : 211 - 214 . DOI: 10.1016/j.nds.2014.07.049 http://doi.org/10.1016/j.nds.2014.07.049 The FLUKA code: developments and challenges for high energy and medical applications Ferrari A , Sala PR , Fasso A , et al . Fluka: A multi-particle transport code. CERN, INFN/TC_05/11, SLAC-R-773 , 2005 . Fluka: A multi-particle transport code. CERN, INFN/TC_05/11, SLAC-R-773 Agostinelli S , Allison J , Amako K , et al . Geant4—a simulation toolkit . NuclInstrum Meth A , 2003 , 506 : 250 - 303 . DOI: 10.1016/S0168-9002(03)01368-8 http://doi.org/10.1016/S0168-9002(03)01368-8 Geant4—a simulation toolkit Allison J , Amako K , Apostolakis J , et al . Geant4 developments and applications . IEEE Trans NuclSci , 2006 , 53 : 270 - 278 . DOI: 10.1109/TNS.2006.869826 http://doi.org/10.1109/TNS.2006.869826 Geant4 developments and applications Ermis E E , Celiktas C , Pilicer E . A method to enhance coincidence time resolution with applications for medical imaging systems (TOF/PET) . RadiatMeas , 2014 , 62 : 52 - 59 . DOI: 10.1016/j.radmeas.2014.01.013 http://doi.org/10.1016/j.radmeas.2014.01.013 A method to enhance coincidence time resolution with applications for medical imaging systems (TOF/PET) Alpat B , Pilicer E , Blasko S , et al . Total and Partial Fragmentation cross-section of 500 MeV/nucleon carbon ions on different target materials . IEEE Trans NuclSci , 2013 , 60 : 4673 - 4682 .  DOI: 10.1109/TNS.2013.2284855  http://doi.org/10.1109/TNS.2013.2284855  Total and Partial Fragmentation cross-section of 500 MeV/nucleon carbon ions on different target materials Tapan I , Pilicer F B . Silicon tracker simulation for the Turkish Accelerator Center particle factory . NuclInstrum Meth A , 2014 , 765 : 240 - 243 . DOI: 10.1016/j.nima.2014.05.100  http://doi.org/10.1016/j.nima.2014.05.100  Silicon tracker simulation for the Turkish Accelerator Center particle factory Alpat B , Pilicer E , Servoli L , et al . Full Geant4 and FLUKA simulations of an e-LINAC for its use in particle detectors performance tests . JINST , 2012 , 7 : P03013 . DOI: 10.1088/1748-0221/7/03/P03013 http://doi.org/10.1088/1748-0221/7/03/P03013 Full Geant4 and FLUKA simulations of an e-LINAC for its use in particle detectors performance tests Tomal A , Mazarro I , Kakuno E M , et al . Experimental detemination of linear attenuation coefficient of normal, benign and malignant breast tissues . RadiatMeas , 2010 , 45 : 1055 - 1059 . DOI: 10.1016/j.radmeas.2010.08.008 http://doi.org/10.1016/j.radmeas.2010.08.008 Experimental detemination of linear attenuation coefficient of normal, benign and malignant breast tissues Böke A . Linear attenuation coefficients of tissues from 1 keV to 150 keV . RadiatPhysChem , 2014 , 102 : 49 - 59 . DOI: 10.1016/j.radphyschem.2014.04.006 http://doi.org/10.1016/j.radphyschem.2014.04.006 Linear attenuation coefficients of tissues from 1 keV to 150 keV Akar A , Baltas H , Cevik U , et al . Measurement of attenuation coefficients for bone, muscle, fat and water at 140, 364 and 662 keV γ-ray energies . JQSRT , 2006 , 102 : 203 - 211 . DOI: 10.1016/j.jqsrt.2006.02.007 http://doi.org/10.1016/j.jqsrt.2006.02.007 Measurement of attenuation coefficients for bone, muscle, fat and water at 140, 364 and 662 keV γ-ray energies CERN Geant4 collaboration physics reference manual . http://geant4.web.cern.ch/geant4/UserDocumentation/UsersGuides/PhysicsReferenceManual/fo/PhysicsReferenceManual.pdf http://geant4.web.cern.ch/geant4/UserDocumentation/UsersGuides/PhysicsReferenceManual/fo/PhysicsReferenceManual.pdf Accessed 22 May 2015 . CERN Geant4 collaboration physics reference manual Battistoni G , Cerutti F , Fasso A , et al . The FLUKA code: description and benchmarking . AIP Conference Proceeding , 2007 , 86 : 31 - 49 . DOI: 10.1063/1.2720455 http://doi.org/10.1063/1.2720455 The FLUKA code: description and benchmarking Moses W W , Derenzo S E . Prospects for time-of-flight PET using LSO scintillator . IEEE Trans NuclSci , 1999 , 46 : 474 - 478 . DOI: 10.1109/23.775565  http://doi.org/10.1109/23.775565  Prospects for time-of-flight PET using LSO scintillator ICRP Report . http://www.icrp.org/docs/ICRP%20Annual%20Report%202012.pdf http://www.icrp.org/docs/ICRP%20Annual%20Report%202012.pdf Accessed 22 May 2015 . ICRP Report ICRU , Tissue Substitutes in Radiation Dosimetry and Measurement, Report 44 of the International Commission on Radiation Units and Measurements (Bethesda, MD) , 1989 . Tissue Substitutes in Radiation Dosimetry and Measurement, Report 44 of the International Commission on Radiation Units and Measurements National Institute of Standards and Technology (NIST) X-ray mass attenuation coefficients . http://physics.nist.gov/PhysRefData/XrayMassCoef/tab4.html http://physics.nist.gov/PhysRefData/XrayMassCoef/tab4.html Accessed 22 May 2015 . National Institute of Standards and Technology (NIST) X-ray mass attenuation coefficients Cullen D E , Hubbell J H , Kissel L . The Evaluated Photon Data Library 97 version (EPDL97) . http://www.ge.infn.it/geant4/temp/saracco/cor/EPDL97.pdf http://www.ge.infn.it/geant4/temp/saracco/cor/EPDL97.pdf Accessed 22 May 2015 . The Evaluated Photon Data Library 97 version (EPDL97) Roach G , Tickner J , Haarlem Y V . Discrepancies in atomic shell and fluorescent X-ray energies in the Evaluated Photon Data Library EPDL97 . X-ray Spectrom , 2012 , 41 : 279 - 283 . DOI: 10.1002/xrs.2392 http://doi.org/10.1002/xrs.2392 Discrepancies in atomic shell and fluorescent X-ray energies in the Evaluated Photon Data Library EPDL97 Botta F , Mairani A , Battistoni G , et al . Calculation of electron and isotopes dose point kernels with fluka Monte Carlo code for dosimetry in nuclear medicine therapy . Med Phys , 2011 , 38 : 3944 - 3954 . DOI: 10.1118/1.3586038 http://doi.org/10.1118/1.3586038 Calculation of electron and isotopes dose point kernels with fluka Monte Carlo code for dosimetry in nuclear medicine therapy Robert C , Dedes G , Battistoni G , et al . Distributions of secondary particles in proton and carbon-ion theraphy: a comparison between GATE/Geant4 and FLUKA Monte Carlo codes . Phys Med Biol , 2013 , 58 : 2879 - 2899 . DOI: 10.1088/0031-9155/58/9/2879 http://doi.org/10.1088/0031-9155/58/9/2879 Distributions of secondary particles in proton and carbon-ion theraphy: a comparison between GATE/Geant4 and FLUKA Monte Carlo codes Medhat M E . Comprehensive study of photon through different construction matters by Monte Carlo simulation . RadiatPhysChem , 2015 , 107 : 65 - 74 . DOI: 10.1016/j.radphyschem.2014.09.005 http://doi.org/10.1016/j.radphyschem.2014.09.005 Comprehensive study of photon through different construction matters by Monte Carlo simulation Esfandiari M , Shirmardi S P , Medhat M E . Element analysis and calculation of the attenuation ceofficients for gold, bronze and water matrixes using MCNP, WinXCOM and experimental data . RadiatPhysChem , 2014 , 99 : 30 - 36 . DOI: 10.1016/j.radphyschem.2014.02.011 http://doi.org/10.1016/j.radphyschem.2014.02.011 Element analysis and calculation of the attenuation ceofficients for gold, bronze and water matrixes using MCNP, WinXCOM and experimental data Singh V P , Medhat M E , Badiger N M . Photon attenuation coefficients of thermoluminescentdosimetric materials by Geant4 toolkit, XCOM program and experimental data: a comparison study . Ann Nucl Energy , 2014 , 68 : 96 - 100 . DOI: 10.1016/j.anucene.2014.01.011 http://doi.org/10.1016/j.anucene.2014.01.011 Photon attenuation coefficients of thermoluminescentdosimetric materials by Geant4 toolkit, XCOM program and experimental data: a comparison study Singh V P , Medhat M E , Badiger N M . Photon energy absorption coefficients for nuclear track detectors using Geant4 Monte Carlo simulation . RadiatPhysChem 2015 , 106 : 83 - 87 . DOI: 10.1016/j.radphyschem.2014.07.001 http://doi.org/10.1016/j.radphyschem.2014.07.001 Photon energy absorption coefficients for nuclear track detectors using Geant4 Monte Carlo simulation 10.1007/s41365-016-0053-2 Supported by Scientific Research Project of Ege University under Project No. 2014 FEN 026 and Uludag University under Project No.OUAP (F) -2012/26. 2015-08-01 2015-08-01 2015-09-04 2016-05-07 2016-06-20 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. 2016 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-016-0053-2 http://dx.doi.org/10.1007/s41365-016-0053-2 . 2016 Nuclear Science and Techniques 03 27