|Year : 2017 | Volume
| Issue : 2 | Page : 114-121
Specific absorbed fractions of internal photon and electron emitters in a human voxel-based phantom: A monte carlo study
Ruhollah Ghahraman Asl1, Ali Asghar Parach2, Shahrokh Nasseri3, Mehdi Momennezhad3, Seyed Rasoul Zakavi4, Hamid Reza Sadoughi5
1 Bioinformatics Research Centre, Department of Nutrition and Biochemistry, Faculty of Medicine, Sabzevar University of Medical Sciences, Sabzevar, Iran
2 Department of Medical Physics, Faculty of Medicine, Shahid Sadoughi University of Medical Sciences, Yazd, Iran
3 Department of Medical Physics, Faculty of Medicine, Mashhad University of Medical Sciences, Mashhad, Iran
4 Nuclear Medicine Research Center, Mashhad University of Medical Sciences, Mashhad, Iran
5 Department of Biotechnology and Molecular Sciences, Faculty of Medicine, North Khorasan University of Medical Sciences, Bojnurd, Iran
|Date of Web Publication||27-Mar-2017|
Ruhollah Ghahraman Asl
Department of Nutrition and Biochemistry, Faculty of Medicine, Sabzevar University of Medical Sciences, Sabzevar
| Abstract|| |
The specific absorbed fraction (SAF) of energy is an essential element of internal dose assessment. Here reported a set of SAFs calculated for selected organs of a human voxel-based phantom. The Monte Carlo transport code GATE version 6.1 was used to simulate monoenergetic photons and electrons with energies ranging from 10 keV to 2 MeV. The particles were emitted from three source organs: kidneys, liver, and spleen. SAFs were calculated for three target regions in the body (kidneys, liver, and spleen) and compared with the results obtained using the MCNP4B and GATE/GEANT4 Monte Carlo codes. For most photon energies, the self-irradiation is higher, and the cross-irradiation is lower in the GATE results compared to the MCNP4B. The results show generally good agreement for photons and high-energy electrons with discrepancies within − 2% ±3%. Nevertheless, significant differences were found for cross-irradiation of photons of lower energy and electrons of higher energy due to statistical uncertainties larger than 10%. The comparisons of the SAF values for the human voxel phantom do not show significant differences, and the results also demonstrated the usefulness and applicability of GATE Monte Carlo package for voxel level dose calculations in nonuniform media. The present SAFs calculation for the Zubal voxel phantom is validated by the intercomparison of the results obtained by other Monte Carlo codes.
Keywords: GATE calculation, internal dosimetry, phantom, specific absorbed fraction value
|How to cite this article:|
Asl RG, Parach AA, Nasseri S, Momennezhad M, Zakavi SR, Sadoughi HR. Specific absorbed fractions of internal photon and electron emitters in a human voxel-based phantom: A monte carlo study. World J Nucl Med 2017;16:114-21
|How to cite this URL:|
Asl RG, Parach AA, Nasseri S, Momennezhad M, Zakavi SR, Sadoughi HR. Specific absorbed fractions of internal photon and electron emitters in a human voxel-based phantom: A monte carlo study. World J Nucl Med [serial online] 2017 [cited 2019 Dec 7];16:114-21. Available from: http://www.wjnm.org/text.asp?2017/16/2/114/203065
| Introduction|| |
The specific absorbed fraction (SAF) is the fraction of emitted energy from the source organ that is absorbed by the target organ per unit mass of target organ. The SAF of energy is an essential element of internal dose assessment. Several organs become a source of radiation after administered radiopharmaceutical to a patient during either diagnostic or therapeutic procedures in nuclear medicine. The accurate and realistic determine of the absorbed dose to the internal organs of the patient is important for radiation protection purposes. The SAF is usually useful quantity to evaluate the stochastic and deterministic effects.
In principal, a full Monte Carlo simulation together with a developed geometrical model of humans is the most accurate approach for calculation of absorbed dose fraction of energy and other dosimetric quantities. Thus, the development of more realistic anatomical models was desirable to achieve a better dose assessment. To do so, voxel phantoms were derived from computed tomography (CT) or magnetic resonance image (MRI) data of real persons and provide detailed information about the human anatomy. The Zubal  phantom is one of the earlier them that obtained from CT (thorax) and MRI (brain). Each organ in the Zubal phantom has a unique id number that identifies its voxels. The aim of this work is to the estimation of monoenergetic photon and electron SAFs for selected organs of a human voxelized phantom using GATE Monte Carlo package for application to internal dosimetry.
| Materials and Methods|| |
The adult male computational phantom, shown in [Figure 1], was used. This voxelized phantom was developed at Yale University by Zubal et al. on the basis of CT data of actual person. The phantom is consisted of a 128 × 128 × 243 array of 4 mm × 4 mm × 4 mm voxels. The Zubal phantom contained 56 segmented organs that a unique ID number (or tissue index) was given to each organ in the organ libraries. [Table 1] shows the number of voxels, volumes, and the masses of several organs in the Zubal phantom. The mass of each organ was determined as the organ volume multiplied by the organ density obtained from the International Commission on Radiation Units and Measurements 44.
|Figure 1: Coronal slice through the Zubal voxel-based anthropomorphic phantom with inserted lesions (red points)|
Click here to view
|Table 1: The number of voxels, volumes and the masses of some selected organs in Zubal voxelized phantom|
Click here to view
GATE Monte Carlo technique
All of the dosimetry simulations in this work were performed by GATE Monte Carlo package (version 6.1), which was built on Geant4.9.3p01. However, the uses of the GATE and GEANT4 have been evaluated for dosimetry applications at previous publications.,,
The physics lists for photon and electron interactions used in the GATE simulations displays in [Table 2]. The cut-off energy (i.e., the limit at which the energy is regarded to be locally absorbed) and the cut-off range were set to 1 keV and 0.1 mm for the photons and the electrons, respectively, in all simulations. This range for electron was regarded as sufficiently small compared with the size of a voxel (4 mm). We used a computer with a 3.5 GHz Intel(R) Core (TM) i7-4770K processor and 8.0 GB RAM running Linux Fedora 13 for Monte Carlo simulation. The numbers of histories for photon and electron source organs were 100 million and 50 million, respectively.
|Table 2: Physics settings for photon and electron interactions in the GATE simulationsa|
Click here to view
Calculation of the photons and electrons specific absorbed fractions
Three combinations of the source and target organs have been considered with initial photon and electron energies ranging between 10 keV and 2 MeV. The activity was assumed to be uniformly distributed in the source organ and the absorbed doses transferred to the target organs of the phantom were determined. The procedure was performed considering the activity in the kidneys, liver and the spleen of the phantom independently. With GATE, the dose actor was used to score the energy absorbed by target organs. The SAF values were then calculated for each source and target organ based on the MIRD formalism,, using below equation:
where rS is a source organ, rT is a target organ, ET is the energy absorbed in rT, ES is the energy emitted from rS, and m is the mass of rT. The same SAF values calculated using the MCNP4B Monte Carlo code,, and the results were compared for quality assurance purpose. The relative percentage difference (RD %) between two corresponding photon and electron SAF values were calculated as:
In this work, the results obtained with MCNP4B were arbitrarily considered as the reference.
| Results|| |
The SAF values for monoenergetic photons and electrons in some organs of the Zubal voxelized phantom were calculated by GATE Monte Carlo package and compared to those reported by other Monte Carlo codes. The selected organ geometry was exactly the same as the organ geometry described by Yoriyaz et al. and Parach et al. in the simulation.
Photon specific absorbed fractions
The SAF values calculated for the monoenergetic photons of 10, 20, 50, 100, 200, 500, 1000, and 2000 keV using GATE Monte Carlo package are presented in [Table 3] for three source and target combinations. The SAF values derived using MCNP4B Monte Carlo code  and the maximum relative percentage difference (RD %) between the GATE and MCNP4B data for each corresponding pair of organs (target ← source) were also included in the table. The photon SAF values derived with GATE and MCNP4B were divided into two groups, self-irradiation, and cross-irradiation. On average, photon SAFs calculated by GATE Monte Carlo was slightly higher values (+0.7%) for self-irradiation and lower values (−4.4%) for cross-irradiation compared to the MCNP4B results.
|Table 3: Specific absorbed fraction values (kg-1) for 10-2000 keV monoenergetic photons derived with GATE and reported data with MCNP4B using Zubal phantom|
Click here to view
In [Figure 2], the SAF values against the photon energy are presented for the self-irradiation data. The graphs visually reveal a good correlation between the two series of data. However, the minimum and maximum of absolute average difference in this photon energy range were 0.1% and 6.1%, respectively, for GATE compared to MCNP4B. The photon SAFs for self-irradiation decreases with increasing photon energy from 10 to 100 keV and begin to increase slightly and then begins to decrease again.
|Figure 2: Photon specific absorbed fraction values derived with the GATE and MCNP4B against the photon energy for self-irradiation of kidneys, liver, and spleen|
Click here to view
The photon SAF values against the energy range were 10–2000 keV for cross-irradiation when liver, kidneys, and spleen as source were displayed in [Figure 3]. The graphs show good correlation between the two series of data but smaller values for GATE compared to MCNP4B. The minimum and maximum of absolute average difference for cross-irradiation were 1.1% and 15.8%, respectively. In [Figure 4], the average relative differences against photon energy are plotted. This figure reveals that the maximum and the minimum differences are with the photon energy of 50 keV and 2000 keV. However, the relative differences for self- and cross-irradiation decreasing in this range of photon energy.
|Figure 3: Photon specific absorbed fraction values derived with the GATE and MCNP4B against the photon energy for cross-irradiation from (a) liver, (b) kidneys, and (c) spleen to other selected organs|
Click here to view
|Figure 4: Average of relative percentage difference between the GATE and MCNP4B for self- and cross-irradiation against the photon energy|
Click here to view
Electron specific absorbed fractions
In [Table 4], the SAF values for monoenergetic electrons of 10, 20, 50, 100, 200, 350, 500, 690, 935, 1200, and 2000 keV calculated with GATE Monte Carlo package are presented for three source and target combinations. The SAF values reported by Parach et al. and the maximum relative percentage difference (RD %) for each corresponding pair of organs were also included in the table. The comparison of electron SAFs was limited to five monoenergetic electrons since these energies were considered in the literature  to represent the average electron energy of 186 Re,24 Na,32 P,90 Y, and 19 O radionuclides, respectively. On average, electron SAFs calculated by GATE Monte Carlo was also slightly higher values (+0.7%) for self-irradiation and lower values (−4.9%) for cross-irradiation than reported data from the literature.
|Table 4: Specific absorbed fraction values (kg-1) for 10-2000 keV monoenergetic electrons derived with GATE and reported data using Zubal phantom|
Click here to view
[Figure 5] illustrates the SAF values against the electron energy for the self-absorption of liver, kidneys, and spleen of Zubal phantom. The graphs visually reveal a good correlation between the two series of data. However, the minimum and maximum of absolute average difference in this electron energy range for self-irradiation were 0% and 3%, respectively, for GATE compared to the corresponding reported data in the literature. The electron SAFs for self-irradiation are constant with increasing electron energy since short electron ranges still in the large source organs. However, slight drop-off of the SAFs at high electron energies was observed.
|Figure 5: Electron specific absorbed fraction values derived with the GATE and reported data against the electron energy for self-irradiation of kidneys, liver, and spleen. As can be seen, the comparison limited to five electron energies|
Click here to view
[Figure 6] shows the cross-absorption electron SAF values for the source in liver, kidneys, and spleen of the Zubal phantom. The graphs show good correlation between the two series of data. The minimum and maximum of absolute average difference for cross-irradiation were 3.3% and 5.7%, respectively. [Figure 6] also shows the electron irradiation of adjacent organs cannot be always neglected, even though electrons are considered as weak penetrating radiation. The values for neighboring organs such as spleen and liver cannot be negligible for electron energy above 690 keV and 935 keV, respectively.
|Figure 6: Electron specific absorbed fraction values derived with the GATE and reported data against the electron energy for cross-irradiation from (a) liver, (b) kidneys and (c) spleen to other selected organs. As can be seen, the comparison limited to five electron energies|
Click here to view
In [Figure 7], the average relative percentage differences against electron energy are plotted. In contrast with photon energy, the average relative percentage differences are increasing with increasing electron energy for both self- and cross-irradiation. In other words, the maximum relative differences between the GATE and  results related to the electron energy 1200 keV for both self- and cross-irradiation.
|Figure 7: Average of relative percentage difference between the GATE and reported data for self- and cross-irradiation against the electron energy|
Click here to view
| Discussion|| |
The estimation of monoenergetic photon and electron SAF values using GATE Monte Carlo package was performed in a humanoid voxel phantom and compared the data with the other well developed Monte Carlo codes. For better comparison, we divided the photon and electron SAF values into two groups, i.e., the self- and the cross-irradiation.
Specific absorbed fractions for photon energy
The comparison between the GATE and MCNP4B photon SAF values showed very good agreement. However, differences observed in some organs SAFs are related to difference Monte Carlo codes, which have been used for calculations. In other study, Pacilio et al. reported differences within 10% between MCNP4C and GEANT4 for self-irradiation voxel S values for both photon and electrons. Relative uncertainties <10% are generally considered reliable and those between 10% and 20% are considered questionable. By this standard, it was desired that all SAFs show uncertainties below 10%. All SAF values for selected source organs were considered reliable (≤5%). However, in this study, the absorbed dose to the organs was considered not to the voxels. Difference attenuation properties in GATE in our study compared to the MCNP4B would certainly result in higher absorption in source organs and, therefore, lower cross-irradiation to target organs. The bias between results that shown in [Figure 4] are almost independent of the photon energy and do not follow a physically explainable trend. However, the relative differences would reduce by increasing of the photon energy for both self- and cross-irradiation due to the difference in the material composition attributed to the tissues in two Monte Carlo codes. The tissue compositions were provided in GATE Monte Carlo package in the present study. It can be also concluded from [Figure 2] that organs with small mass obtain larger photon self SAFs than large mass organs (e.g., in here, spleen).
Visually assessment of SAF curves for photon self-irradiation in [Figure 2] was revealed, the source organ is also the target organ for the photon energy ranges from 10 keV to 2000 keV. Nevertheless, photon self SAF values decrease with increasing photon energy not monotonically due to the highest escape probability of the scattered photons when the energy increases.
The assessment of plotted curves in [Figure 3] also confirmed that the organs geometry (such as size and shape), density and the interorgan distance have the significant influence on the photon cross-irradiation SAF values.,
Specific absorbed fractions for electron energy
The data in [Table 4] show low relative differences (up to ± 13.7%) between two series data for monoenergetic electrons of 350, 500, 690, 935, and 1200 keV. However, the SAF values for the cross-irradiation of electrons were small in magnitude and therefore, the differences were most likely due to high statistical uncertainties in data. The average differences between data for all electron energies were low (−2.1% ±3.9%).
For electron SAFs self-irradiation [Figure 5], the trend versus energy is approximately constant due to the electron stopping power trend versus energy. But for electron energies higher than about 500 keV, the SAFs slightly decrease with energy.
The cross-irradiation SAF values for electrons increase with energy at high energies, as already reported in the literature. A close look at [Figure 6]b reveals that the SAF values for liver and spleen from kidneys as source organ are similar, but [Figure 6]a and c show the cross-absorption electron SAF values for the source in liver and spleen, respectively. In other words, [Figure 6] shows the electron irradiation of adjacent organs cannot be always neglected, even though electrons are considered as weak penetrating radiation. For example, according to the [Figure 6]a and c, the SAF values for neighboring organs such as spleen and liver cannot be negligible for electron energies above 50 keV and 350 keV, respectively. From this comparison and GATE Monte Carlo calculations of electron cross-irradiation SAF values, the assumption of MIRD  and ICRP publication No. 30 (i.e., the electrons are fully absorbed in the source organ itself) not always correct and led to underestimate the absorption of the neighboring organs around the source organs.
As already reported in the literature, the statistical uncertainties associated with the voxel S values depend on the voxel size, distance between the source and target voxel, energy, and type (electrons or photons) of starting particles.
| Conclusion|| |
There was a very good agreement between the two series of data for both photon and electron SAF values. The average relative differences decrease against the photon energy and increase against the electron energy for both self- and cross-irradiation. From results of this work, it can be concluded that the estimation of SAF values and other useful internal dosimetric quantities could be feasibly calculated for various organs of anthropomorphic voxel-based phantom using GATE Monte Carlo package with reasonably statistical uncertainty (≤5%). This comparison also has been confirmed that the SAFs for self-irradiation depended on the energy and the mass of the target and source organs, the SAFs for cross-irradiation depended on the relative position of source and target organs. The SAF values obtained using the GATE Monte Carlo package for real phantom and connected to the individual biokinetic data could make the patient 3D internal dose calculation possible with reliable uncertainty in nuclear medicine.
The work underlying this paper was financially supported by a grant from Sabzevar University of Medical Sciences.
Financial support and sponsorship
Office of the vice-president for research in Sabzevar University of Medical Sciences.
Conflicts of interest
There are no conflicts of intrest.
| References|| |
Zubal IG, Harrell CR, Smith EO, Rattner Z, Gindi G, Hoffer PB. Computerized three-dimensional segmented human anatomy. Med Phys 1994;21:299-302.
ICRU. International Commission on Radiation Units and Measurements. Tissue substitutes in radiation dosimetry and measurement. ICRU Report 44; 1989.
Visvikis D, Bardies M, Chiavassa S, Danford C, Kirov A, Lamare F, et al
. Use of the GATE Monte Carlo package for dosimetry applications. Nucl Instrum Methods A 2006;569:335-40.
Assié K, Gardin I, Véra P, Buvat I. Validation of the Monte Carlo simulator GATE for indium-111 imaging. Phys Med Biol 2005;50:3113-25.
Thiam CO, Breton V, Donnarieix D, Habib B, Maigne L. Validation of a dose deposited by low-energy photons using GATE/GEANT4. Phys Med Biol 2008;53:3039-55.
Loevinger R, Budinger TF, Watson EE. MIRD Primer for Absorbed Dose Calculations. New York: Society of Nuclear Medicine; 1988.
Bolch WE, Eckerman KF, Sgouros G, Thomas SR. MIRD pamphlet No 21: A generalized schema for radiopharmaceutical dosimetry – standardization of nomenclature. J Nucl Med 2009;50:477-84.
Yoriyaz H, dos Santos A, Stabin MG, Cabezas R. Absorbed fractions in a voxel-based phantom calculated with the MCNP-4B code. Med Phys 2000;27:1555-62.
Parach AA, Rajabi H, Askari MA. Paired organs-should they be treated jointly or separately in internal dosimetry? Med Phys 2011;38:5509-21.
Pacilio M, Lanconelli N, Lo MS, Betti M, Montani L, Torres AL, et al.
Differences among Monte Carlo codes in the calculations of voxel S values for radionuclide targeted therapy and analysis of their impact on absorbed dose evaluations. Med Phys 2009;36:1543-52.
Pelowitz DB. MCNPX User's Manual Version 2.5.0. Los Alamos, NM: Los Alamos National Laboratory; 2005.
Lanconelli N, Pacilio M, Lo Meo S, Botta F, Di Dia A, Aroche AT, et al.
A free database of radionuclide voxel S values for the dosimetry of nonuniform activity distributions. Phys Med Biol 2012;57:517-33.
Ghahraman Asl R, Nasseri S, Parach AA, Zakavi SR, Momennezhad M, Davenport D. Monte Carlo and experimental internal radionuclide dosimetry in RANDO head phantom. Australas Phys Eng Sci Med 2015;38:465-72.
Momennezhad M, Nasseri S, Zakavi SR, Parach AA, Ghorbani M, Ghahraman Asl R. A 3D Monte Carlo method for estimation of patient-specific internal organs absorbed dose for 99m
– Octreotide imaging. World J Nucl Med 2016;15:114-23.
] [Full text]
Snyder WS, Ford MR, Warner GG. Estimates of specific absorbed fractions for photon sources uniformly distributed in various organs of a heterogeneous phantom. MIRD Pamphlet No. 5 Revised. New York: Society of Nuclear Medicine; 1978.
ICRP. Limits for Intakes of Radionuclides by Workers: Part 1. ICRP Publication 30. Oxford: Pergamon; 1979.
[Figure 1], [Figure 2], [Figure 3], [Figure 4], [Figure 5], [Figure 6], [Figure 7]
[Table 1], [Table 2], [Table 3], [Table 4]