Calculation of beta absorbed fractions for iodine isotopes in

Volume 11, No 2
International Journal of Radiation Research, April 2013
Calculation of beta absorbed fractions for iodine
isotopes in ellipsoidal thyroid lobe
A.A. Mowlavi1*, M. Mirzaei2, M.R. Fornasier3, M. de Denaro3
1PhysicsDepartment,SchoolofSciences,HakimSabzevariUniversity,Sabzevar,Iran
2PhysicsDepartment,facultyofShahidChamran,Technical&VocationalUniversity,Kerman,Iran
3S.C.diFisicaSanitaria,A.O.U.“OspedaliRiuniti”diTrieste,Trieste,Italy
ABSTRACT
► Original article
*Correspondingauthor:
Dr.AliAsgharMowlavi,
Fax:+985714003170
E‐mail:[email protected]
Received: Feb. 2012
Accepted: May 2012
Int. J. Radiat. Res., April 2013;
11(2): 121-126
Background: The thyroid gland absorbs nearly all the iodine in the blood,
independently of its isotopic composi on. When a large enough ac vity of
radioac ve iodine like 131I, 132I, 133I is taken into the body, it can destroy the
healthy thyroid gland cells as well as the cancer’s cells. In fact, as it is well
known, some isotopes of iodine are used not only for acquiring thyroid
images but also for curing thyroid cancer or hyperthyroidism due to Graves’
disease. Moreover, some of them are released in nuclear accidents. The aim
of this work is to evaluate the absorbed frac on of beta rays from different
iodine radioisotopes in thyroid lobe, using Monte Carlo method. Materials and Methods: We have applied the MCNPX code to calculate the beta
absorbed frac ons for the most important iodine radioisotopes in the thyroid
lobes, supposed to have an ellipsoidal shape, with the volume varying from 1
ml to 25 ml and the material composi on suggested by ICRP. The beta rays
spectra of iodine isotopes have been taken from the LBNL website. Results: The result showed that the volume lobe varia on had a significant effect on
the absorbed frac on for beta rays in thyroid gland, up to 25% for 124I. The
absorbed frac ons of beta rays were decreasing from 123I to 131I, 130I, 126I, 135I,
133
I and 124I respec vely. Conclusions: Decreasing of the absorbed frac on
might be related to the beta rays energy spectra of the isotopes. Moreover,
for 131I results obtained for beta absorbed frac on in spheres was in
agreement with previously published results.
Keywords: Thyroid gland; Beta absorbed fraction; Iodine isotopes; MCNPX
code.
INTRODUCTION
life of 8.05 days, is readily available and has a
low cost. Therefore, it is widely used for a
Therearemanydifferentiodineisotopes(up
number of nuclear medicine procedures,
to 37) and some of them have an important
including monitoring and tracing the low of
usage formedicalapplications (1‐20).Radioactive
thyroxin from the thyroid, imaging of thyroid,
iodine isotopes have the same physical
evaluationoftheiodineuptakecurve,therapyof
properties as stable iodine‐127 and thyroid
thyroid benign diseases or cancer (1‐4). Iodine‐
gland cannot distinguish between radioactive
123hasapropergammaemission(159keV)for
and non‐radioactive iodine. Iodine‐131, with a
thyroidimagingbygammacamerawithaphysi‐
gamma emission of 365keV and a physical half‐
cal half‐life of 13.1 hours. Therefore, it is used
Mowlavi et al. / Beta absorbed fractions for iodine isotopes in thyroid
~
forstudiesofthyroidmetabolismandaffordsan
A  (1)
D


excellentimagequality.Aswellas,using 123Ifor
m
diagnostic purposes reduces the radiation risk,
Whereɸβisthebetaabsorbedfraction,Aisthe
because of the short half‐life, low intensity and
cumulated
activity, Δβ is the average energy of
low energy of Auger electron. Also, it hasn’t the
beta
spectrum
and m is the mass of target
same stunning effect on functioning of thyroid
region.
The
result
of our study can be used for
tissue as 131I, with an advantage for the
dose calculation of iodine isotopes in thyroid
subsequentradioiodinetreatment(5,6).
glandaccordingtoMIRDformula.
Iodine‐124 can be used to identify mediasti‐
nal micro metastases in thyroid carcinoma by
PET/CT imaging. Freudenberg et al. reported
MATERIALSANDMETHODS
that Iodine‐124 PET/CT imaging is a promising
techniquetoimprovetreatmentplanninginthy‐
For simulation of beta transport, we have
roid cancer.6‐8 Iodine–125 is used in permanent
applied the MCNPX code (22). It is a general
lowdoseradiationbrachytherapy,especiallyfor
purpose, continuous and discrete energy,
localized prostate cancer (9, 10). Farran had
132
generalized‐geometry, time‐dependent code to
reported the use of I in the diagnosis of
simulate particles transport, based on Monte
patients undergoing treatment for thyroid
Carlo method. It contains lexible source and
disordersin1969(15)aswellas,ithadbeenused
tally options and utilizes the last nuclear cross
by Goolden et al. many years ago in studies of
(16)
sectionlibraries.
thyroidfunction .
Iodine neutron rich isotopes which are Thecodehasbeenusedtocalculatethebeta
rays absorbed fractions of main radioactive
producedinsidethereactorsmaybereleasedin
nuclear accidents, like Chernobyl in 1986, or iodineisotopes,withmassnumberrangingfrom
123to135,inathyroidlobeofellipsoidalshape,
disasters, as Fukushima in 2011 (17‐19). There‐
withthemajoraxistwotimesoftheminoraxis.
fore,theunwantedcontaminationofradioactive
We must mention that Iodine‐123 decays by
iodine can cause serious thyroid problems, the
electron capture, therefore its Auger electron
primaryriskbeingthyroidcancer,withchildren
spectrum is used for our simulation. The tissue
moreatriskthanadults.Strandetal.published
densityofthyroidgland hasbeenconsideredof
theirstudyaboutthethyroiduptakenotonlyof
1.05g/cm3andthevolumevaryingfrom1mlto
iodine‐131 but also of iodine‐133 in the
25 ml. Considering thyroid lobes to have an
population of southern Sweden from Chernobyl
ellipsoidalshape,themeandistanceoftwolobes
accident(17).
has been taken from the study of the
Inapreviousarticle,wereportedthethyroid
131
experimental images of more than 50 patients.
volume’s in luence on I energy deposition
Figure 1 shows the scintigraphic image of a
evaluated using the Monte Carlo simulation
thyroid and the corresponding model of a
method (20). Recently, Amato et al. developed a
thyroidlobe.
Monte Carlo simulation in Geant4 to calculate
the absorbed fractions for some beta rays
emittersuchas199Au,177Lu,131I,153Sm,186Reand
90Y,uniformlydistributedinellipsoidalvolumes
of soft tissue (21). In this study, we have
calculated the beta absorbed fractions for main
of iodine radioisotopes in the thyroid lobes,
supposed to have an ellipsoidal shape and a
volumevaryingfrom1mlto25ml.
Following the MIRD approach, the average
Figure 1. A real image of thyroid and the model of one thyroid
betadosetoatargettissuecanbecalculatedas:
lobe.
Int. J. Radiat. Res., Vol. 11 No. 2, April 2013
122
Mowlavi et al. / Beta absorbed fractions for iodine isotopes in thyroid
Absorbed fraction
Absorbed Energy(MeV)
0.7
The beta rays spectra of the iodine isotopes
I-124
havebeentakenfromtheLBNLIsotopesProject
0.6
website (http://ie.lbl.gov/toi), whose last
0.5
version2.1hasbeenreleasedinJanuary2004by
0.4
I-133
FirestoneandEkströ m(23).
Tosimulatethebody,theadult70kghuman
I-126
0.3
I-135
MIRD5 phantom has been used, where the
I-130
0.2
source organ was the thyroid gland with a uni‐
I-131
(24)
form iodine radioisotope distribution . Also,
0.1
theneckhasbeensimulatedwithmoredetailed
I-123
0.0
organs including skin, adipose layer under the
0
2
4
6
8 10 12 14 16 18 20 22 24 26
skin, bone, spinal cord, thyroid lobes, and the
Thyroid lobe volume (ml)
remainingpartassofttissue.Thesofttissuehas
Figure 2. The absorbed energy of beta rays in ellipsoidal
been considered having 1.05 g/cm3 density and
thyroid lobe against the lobe volume.
theICRPcomposition.Wehaveusedadualcore
1.00
PCwith3GHzspeedforrunningMCNPXcode,60
0.98
0.96
minutespereachrunningtoreducetheerrorof
0.94
thecalculations.
0.92
0.90
0.88
RESULTS
Energy deposition of beta rays depends on
sizeandgeometryofthesourceregion,distribu‐
tion of the radioactive material, types of the ra‐
diations and energies emitted by the radionu‐
clide.Toobtaintheabsorbedfractionfβ,wehave
calculated the total absorbed energy per decay
forbetaraysofeachiodineisotopesasshownin
igure 2; then the absorbed fraction has been
evaluatedusingthefollowingequation:
I-123
I-131
I-130
I-126
I-135
I-133
I-124
0.86
0.84
0.82
0.80
0.78
0.76
0.74
0
2
4
6
8
10
Figure3showsthebetaabsorbedfractionfor
7differentiodineisotopesagainstthevolumeof
thyroidlobe.Theresultsshowthattheabsorbed
fractionofbetaraysincreasesbyvolumeforall
the isotopes, and the variation is higher for 124I
than for other isotopes. It can be seen that the
betaabsorbedfractionfor 124Irangesfrom0.75
for a small volume (1 ml) up to 0.91 for a high
volume (25 ml). For 131I the variation is from
0.95till0.98,andfor123Iit’saround1%,because
of low energy Auger electron emitted from this
radioisotope.
123
14
16
18
20
22
24
26
Figure 3. The varia on of absorbed frac on of beta rays by
thyroid lobe volume.
DISCUSSION
Beta' s energy absorbed in the lobe (2)
 
mean energy of beta ray
12
Thyroid lobe volume (ml)
We have used a suitable analytical function
whichwasintroducedbyAmatoetal. (25)forthe
relationship between the absorbed fraction fβ
and the generalized radius ρ, with 2 it
parametersρ0ands:
  
 ()  1  0s 
  
1
(3)
The volume and surface area of an ellipsoid
aregivenby:
4

V  3 abc


p
p
p 1/ p


S  4  (ab)  (bc)  (ac) 
; p  1.6075.
3



(4)
Int. J. Radiat. Res., Vol. 11 No. 2, April 2013
Mowlavi et al. / Beta absorbed fractions for iodine isotopes in thyroid
Based on the sphere radius, “Generalized
iscorrectandsuitableinordertoobtainthebeta
Radius” for ellipsoid can be de ined as three
absorbed fraction as a function of the general‐
timesofvolumetosurfaceratio:
izedradius.
Amato et al. (21) have reported the 131I beta
V
(5)
3
absorbed fractions for ellipsoidal shape with
S
different ratio of axes and various masses of
Figure 4 shows the absorbed fraction of the
thyroid.Ourresultsareingoodagreementwith
betarayvariationagainstthegeneralizedradius
the result of Amato el al. considering the beta
of ellipsoid. Also, in this igure it can be seen a
absorbedfractionsforspheres(listedintable2).
good match of the itted curves (drawn as
continuous lines) to the obtained data. The it
Table 2. Comparison of our result for beta absorbed
frac ons in spheres with the result of Amato el al. (21)
parametersoftheanalyticalfunctioninequation
3 and the reduced chi‐square (Χ2/DoF) values
ϕ2 ϕ1 foralloftheiodineisotopesarelistedintable1.
Amato σ = ϕ1 ‐ ϕ2 V (cm3) This work The obtained value of ρ0 ranges from 0.1747 to
et al. 0.0068 (a variation of about 25.7 times!), while
1
0.899
0.90
-0.001
thesvalueisnotvariedtoomuchanditsvalueis
0.5
0.940
0.94
0
around1.Theseresultsshowthattheequation3
1
0.953
0.96
-0.007
2
0.963
0.97
-0.007
1.00
10
0.978
0.98
-0.002
0.98
20
0.983
0.98
0.003
0.96
0.94
0.92
CONCLUSION
0.90
0.88

0.86
0.84
I-123
I-131
I-130
I-126
I-135
I-133
I-124
0.82
0.80
0.78
0.76
0.74
0.72
0.70
0.2
0.4
0.6
0.8
1.0
1.2
1.4
1.6
(cm)
Figure 4. Absorbed frac on of the beta rays in ellipsoids
varia on against ρ, the generalized radius (lines curves are the
fi ed func ons).
The results show that considering the lobe
volume variation has a signi icant effect on the
absorbedfractionforbetarays.Itcanbedueto
the short range of beta rays and to increasing
surfacetovolumeratioofellipsoidallobe,which
causes some beta particles near the surface to
depositafractionoftheirenergyoutsideofthe
lobe. Moreover, our results for beta absorbed
fraction in spheres are in good agreement with
theresultofAmatoelal.(21)for131I.
Table 1. Fit parameters and reduced chi‐square values for all of the iodine isotopes.
Radioisotope For 124I
For 133I
For 135I
For 126I
For 130I
For 131I
For 123I
ρ0 0.1747 ± 2.8E‐4
0.0935 ± 1.4E‐4
0.0656 ± 7E‐5
0.061 ± 7E‐5
0.0461 ± 8E‐5
0.0301 ± 1.2E‐4
0.0068 ± 2E‐5
Int. J. Radiat. Res., Vol. 11 No. 2, April 2013
s 1.1617 ± 4.72E‐3
1.0934 ± 4.37E‐3
1.0908 ± 3.12E‐3
1.0664 ± 3.3E‐3
1.0447 ± 5.1E‐3
1.0538 ± 1.13E‐2
1.0835 ± 9.04E‐5
Χ2/DoF 7.322E‐7
2.532E‐7
7.169E‐8
7.067E‐8
1.033E‐7
2.331E‐7
8.619E‐9
124
Mowlavi et al. / Beta absorbed fractions for iodine isotopes in thyroid
ACKNOWLEDGMENT
Med, 29: 760‐767.
9. Grimm PD, Blasko JC, Sylvester JE, Meier RM,
Cavanagh W (2001) 10‐year biochemical (prostate‐
The authors would like to thank Prof. D. specific an gen) control of prostate cancer with 125I
Treleani head of TRIL program at ICTP, Trieste,
brachytherapy. Int J Radia on Oncology Biol Phys,
Italy,forhissupportthiswork.
51: 31–40.
10. Stock RG, Stone NN (2010) Current topics in the
treatment of prostate cancer with low‐dose‐rate
brachytherapy. Urol Clin North Am, 37: 83‐96.
REFERENCES
11. AllahAbadia A, Daykin J, Sheppard MC, Gough CSCL,
Franklyn JA (2001) Radiological treatment of hyper‐
1. Vini L and Harmer C (2000) Radioiodine treatment
thyroidism‐progonos c factors for outcome. J Clin
for differen ated thyroid cancer. Clin Oncl (Rcoll
Endocrinol Metab, 86: 3611‐3617.
Radiol), 12: 365‐72.
12. Cooper DS, Doherty GM, Haugen BR, Kloos RT, et al.
2. Matheoud R, Canzi C, Reschini E, Zito F, Vol ni F,
(2006) Management guidelines for pa ents with
Gerundini P (2003) Tissue‐specific dosimetry for
thyroid nodules and differen ated thyroid cancer.
radioiodine therapy of the autonomous thyroid
Thyroid, 16: 109‐141.
nodule. Med Phys, 30: 791‐8.
13. Traino AC, di Mar no F, Lazzeri M. Stabin MG (2000)
3. Kim WG, Ryu JS, Kim EY, Lee JH, Baek JH, Yoon JH,
Influence of thyroid volume reduc on on calculated
Hong SH, Kim ES, Kim TY, Kim WB, Shong YK (2010)
dose in radioiodine therapy of Graves’
Empiric high‐dose 131‐Iodine therapy lacks efficacy
hyperthyroidism. Phys Med Biol, 45: 121–129.
for treated papillary thyroid cancer pa ents with
14. Wildman JC, Powsner ER, Plato PA (1980) Absorbed
detectable serum thyroglobulin, but nega ve cervi‐
frac ons for photons of 125I and 129I in the thyroid.
18
cal sonography and F‐fluorodeoxyglucose positron
Int J Appl Radiat Isot, 31: 421‐423.
emission tomography scan. J Clin Endocrinol Metab,
15. Farran HE (1969) The use of 132I in the diagnosis of
95: 1169‐1173.
pa ents undergoing treatment for thyroid disorders.
4. Sgouros G, Hobbs RF, Atkins FB, Nostrand DV, Laden‐
Br J Radiol, 42: 235.
son PW, Wahl RL (2011) Three‐dimensional radiobio‐
16. Goolden AW, Mallard JR (1958) The use of iodine
logical dosimetry (3D‐RD) with 124I PET for 131I thera‐
132 in studies of thyroid func on. Br J Radiol, 31:
py of thyroid cancer. Eur J Nucl Med Mol Imaging,
589‐95.
38 : S41–S47.
17. Strand SE, Erlandsson K, Löwenhielm P (1988)
5. Urhan M, Dadparvar S, Mavi A, Houseni M, Cham‐
Thyroid uptake of iodine‐131 and iodine‐133 from
roonrat W, Alavi A, Mandel SJ (2007) Iodine‐123 as a
Chernobyl in the popula on of southern Sweden. J
diagnos c imaging agent in differen ated thyroid
Nucl Med, 29: 1719‐2.
carcinoma: a comparison with iodine‐131 pos reat‐
18. Tagami K and Uchida S (2011) Can we remove iodine
ment scanning and serum thyroglobulin measure‐
‐131 from tap water in Japan by boiling? – Experi‐
ment. Eur J Nucl Med Mol Imaging, 34: 1012–1017.
mental tes ng in response to the Fukushima Daiichi
6. Barwick T, Murray I, Megadmi H, Drake WM, Plow‐
Nuclear Power Plant accident. Chemosphere, 84:
man PN, Akker SA, Chew SL, Grossman AB, Avril N
1282‐4.
(2010) Single photon emission computed tomogra‐
19. Jaworowski Z and Kownacka L (1988) Tropospheric
phy (SPECT)/computed tomography using Iodine‐
and stratospheric distribu ons of radioac ve iodine
123 in pa ents with differen ated thyroid cancer:
and cesium a er the Chernobyl accident. J Environ
addi onal value over whole body planar imaging and
Radioact, 6: 145‐150.
SPECT. Eur J Endocrinol, 162: 1131–1139.
20. Mowlavi AA, Fornasier MR, de Denaro M (2011) Thy‐
7. Sundram FX (2006) Clinical use of PET/CT in thyroid
roid volume’s influence on energy deposi on from
131
cancer diagnosis and management. Biomedical Im‐
I calculated by Monte Carlo (MC) simula on. Radi‐
aging and Interven on Journal, 2: 1‐5.
ol Oncol, 45: 143‐146.
8. Eschmann SM, Reischl G, Bilger K, Kupferschläger J,
21. Amato E, Lizio D, Baldari S (2009) Absorbed frac ons
et al. (2002) Evalua on of dosimetry of radioiodine
in ellipsoidal volumes for β‐ radionuclides employed
therapy in benign and malignant thyroid disorders
in internal radiotherapy. Phys Med Biol, 54: 4171‐
by means of iodine‐124 and PET. Euro J Nuclear
4180.
125
Int. J. Radiat. Res., Vol. 11 No. 2, April 2013
Mowlavi et al. / Beta absorbed fractions for iodine isotopes in thyroid
22. Waters LS (ed), (2002) “MCNPX User’s Manual, ver‐
sion 2.3.0”, LA‐UR‐02‐2607, Los Alamos, NM, USA.
23. Firestone RB and Ekström LP (1999) LBNL Isotopes
Project ‐ LUNDS University, Table of Radioac ve
Isotopes. Version 2.1. The website:
h p://
ie.lbl.gov/toi/nuclide.asp?iZA=530131
24. Snyder W, Ford M, Warner G (1969) Es mates of
Int. J. Radiat. Res., Vol. 11 No. 2, April 2013
absorbed frac ons for monoenerge c photon
sources uniformly distributed in various organs of a
heterogeneous phantom: MIRD pamphlet no. 5. J
Nucl Med, 10: 5–52.
25. Amato E, Lizio D, Baldari S (2011) Absorbed frac‐
ons for electrons in ellipsoidal volumes. Phys Med
Biol, 56: 357–365.
126