Therapeutic Applications of Monte Carlo Calculations in Nuclear Medicine Series in Medical Physics and Biomedical Engineering (IoP, 2003)
Therapeutic Applications of Monte Carlo Calculations in Nuclear Medicine Series in Medical Physics and Biomedical Engineering (IoP, 2003)
Therapeutic Applications of Monte Carlo Calculations in Nuclear Medicine Series in Medical Physics and Biomedical Engineering (IoP, 2003)
THERAPEUTIC APPLICATIONS
OF MONTE CARLO CALCULATIONS
IN NUCLEAR MEDICINE
Edited by
Habib Zaidi, PhD
Division of Nuclear Medicine,
Geneva University Hospital,
Switzerland
Series Editors:
C G Orton, Karmanos Cancer Institute and Wayne State University, Detroit,
USA
J A E Spaan, University of Amsterdam, The Netherlands
J G Webster, University of Wisconsin-Madison, USA
IFMBE
The IFMBE was established in 1959 to provide medical and biological
engineering with an international presence. The Federation has a long history
of encouraging and promoting international cooperation and collaboration in
the use of technology for improving the health and life quality of man.
The IFMBE is an organization that is mostly an affiliation of national
societies. Transnational organizations can also obtain membership. At present
there are 42 national members and one transnational member with a total
membership in excess of 15 000. An observer category is provided to give
personal status to groups or organizations considering formal affiliation.
Objectives
. To reflect the interests of the affiliated organizations.
. To generate and disseminate information of interest to the medical and
biological engineering community and international organizations.
. To provide an international forum for the exchange of ideas and concepts.
. To encourage and foster research and application of medical and bio-
logical engineering knowledge and techniques in support of life quality
and cost-effective health care.
. To stimulate international cooperation and collaboration on medical
and biological engineering matters.
. To encourage educational programmes which develop scientific and
technical expertise in medical and biological engineering.
Activities
The IFMBE has published the journal Medical and Biological Engineering
and Computing for over 35 years. A new journal Cellular Engineering was
established in 1996 in order to stimulate this emerging field in biomedical
engineering. In IFMBE News members are kept informed of the develop-
ments in the Federation. Clinical Engineering Update is a publication of
our division of Clinical Engineering. The Federation also has a division for
Technology Assessment in Health Care.
Every three years the IFMBE holds a World Congress on Medical Physics
and Biomedical Engineering, organized in cooperation with the IOMP and
the IUPESM. In addition, annual, milestone and regional conferences are
organized in different regions of the world, such as Asia Pacific, Baltic,
Mediterranean, African and South American regions.
The administrative council of the IFMBE meets once or twice a year and is
the steering body for the IFMBE. The council is subject to the rulings of the
General Assembly which meets every three years.
For further information on the activities of the IFMBE, please contact Jos
A E Spaan, Professor of Medical Physics, Academic Medical Centre,
University of Amsterdam, PO Box 22660, Meibergdreef 9, 1105 AZ, Amster-
dam, The Netherlands. Tel: 31 (0) 20 566 5200. Fax: 31 (0) 20 691 7233. E-
mail: [email protected]. WWW: http://www.ifmbe.org/.
IOMP
The IOMP was founded in 1963. The membership includes 64 national
societies, two international organizations and 12 000 individuals. Member-
ship of IOMP consists of individual members of the Adhering National
Organizations. Two other forms of membership are available, namely
Affiliated Regional Organization and Corporate members. The IOMP is
administered by a Council, which consists of delegates from each of the
Adhering National Organizations; regular meetings of Council are held
every three years at the International Conference on Medical Physics
(ICMP). The Officers of the Council are the President, the Vice-President
and the Secretary-General. IOMP committees include: developing countries,
education and training; nominating; and publications.
Objectives
. To organize international cooperation of medical physics in all its
aspects, especially in developing countries.
. To encourage and advise on the formation of national organizations of
medical physics in those countries which lack such organizations.
Activities
Official publications of the IOMP are Physiological Measurement, Physics
in Medicine and Biology and the Series in Medical Physics and Biomedical
Engineering, all published by the Institute of Physics Publishing. The
IOMP publishes a bulletin Medical Physics World twice a year.
Two council meetings and one General Assembly are held every three years at
the ICMP. The most recent ICMPs were held in Kyoto, Japan (1991), Rio de
Janeiro, Brazil (1994), Nice, France (1997) and Chicago, USA (2000). These
conferences are normally held in collaboration with the IFMBE to form the
World Congress on Medical Physics and Biomedical Engineering. The IOMP
also sponsors occasional international conferences, workshops and courses.
For further information contact: Hans Svensson, PhD, DSc, Professor,
Radiation Physics Department, University Hospital, 90185 Umeå, Sweden.
Tel: (46) 90 785 3891. Fax: (46) 90 785 1588. E-mail: Hans.Svensson@radfys.
umu.se. WWW: http://www.iomp.org.
Contents
LIST OF CONTRIBUTORS xv
PREFACE xviii
vii
viii Contents
16 SUMMARY 348
BIOSKETCHES 350
INDEX 361
List of contributors
xv
xvi List of contributors
Katherine S Kolbert, MS
Memorial Sloan-Kettering Cancer Center, Department of Medical Physics,
1275 York Avenue, New York 10021, USA
This book provides a review of the Monte Carlo method as it is applied in the
field of therapeutic nuclear medicine. Driven in part by the remarkable
increase in computing power and its ready and inexpensive availability,
this is a relatively new yet rapidly expanding field. Likewise, although the
use of radionuclides for therapy has origins dating back almost to the
discovery of natural radioactivity itself, radionuclide therapy and, in
particular, targeted radionuclide therapy has only recently emerged as a
promising approach for therapy of cancer and, to a lesser extent, other
diseases. An effort has, therefore, been made to place the reviews provided
in this book in a broader context. The effort to do this is reflected by the
inclusion of chapters that do not directly address Monte Carlo techniques,
but rather provide an overview of issues that are closely related to therapeutic
nuclear medicine and to the potential role of Monte Carlo in this field. A
brief overview of each chapter is provided below.
Chapters 1, 2 and 6 review the fundamental theory of Monte Carlo,
techniques that have been used in nuclear medicine dosimetry and Monte
Carlo codes currently available for carrying out absorbed dose calculations.
Unlike external radiotherapy or brachytherapy, the kinetics and biodistribu-
tion of internally administered radionuclides must be measured. Imaging is,
therefore, an integral part of nuclear medicine dosimetry. As such, imaging
considerations are incorporated in almost every chapter of the book. The
contributions of Monte Carlo in nuclear medicine imaging are specifically
reviewed in chapter 3. The anatomical input required for Monte Carlo-
based dose calculations may be provided either as idealized geometric
representations of anatomy or as discrete, voxel-based representations that
reflect the actual shape and composition of a human. Again, these
considerations are present throughout the book, but chapters 5 and 13
focus on this specifically for the whole body and for the bone marrow,
respectively, the latter being a particularly important and challenging portion
of the anatomy. The role of Monte Carlo in providing a means for validating
the overall scheme of imaging-based, patient-specific absorbed dose
calculations is reviewed in chapter 12. Chapters 9, 14 and 15 examine specific
xviii
Preface xix
of this book is rapidly evolving and the editors felt it important to minimize
the time required to get this book into press. That this problem was so very
well overcome speaks of the commitment and dedication of the contributors.
We found compilation of this book to be a rewarding and educational
experience and hope that the reader is left with the same experience.
Chapter 1
1.1. INTRODUCTION
1
2 The Monte Carlo method: theory and computational issues
dice-throwing method (the Urim and Thummim) to divine the will of God in
a given situation [3]. Dr de Buffon showed that by randomly tossing needles
on to a table on which parallel lines were drawn, one could determine experi-
mentally the value of pi [4] (it is said that this experiment started by de
Buffon throwing French stick loaves over his shoulder on to a tiled floor
and counting the number of times the loaves fell across the lines between
tiles on the floor). In 1899, Lord Rayleigh showed how the use of a one-
dimensional random walk could provide an approximate solution to a
parabolic differential equation [5]. In 1931 Kolmogorov showed the
relationship between Markov stochastic processes and certain integro-
differential equations [6].
Ulam in his autobiography [7] writes that shortly after his return to Los
Alamos, probably in 1946–1947, he gave two seminars that ‘had good or
lucky ideas and led to successful further development’. One was on what
became known as the Monte Carlo method. He describes an earlier discus-
sion of the basic idea with von Neumann as they drove from Los Alamos
to Lamy. This part of the book is vague on chronology—the seminar may
have been as late as 1948. The abbreviated bibliography lists one paper by
Ulam on Monte Carlo [8]. About 1948 Fermi, Metropolis and Ulam
obtained Monte Carlo estimates for the eigenvalues of the Schrödinger
equation [1] (the reference does not include Fermi, but he was involved in
the derivation). Ulam and von Neumann coined the name ‘Monte Carlo’
in about 1944 in the Manhattan Project during World War II, taken from
the name of the well-known city in Monaco famous for its games of
chance. Both of them led the modern Monte Carlo age and were pioneers
in the development of Monte Carlo techniques and their realizations on
digital computers.
The earliest application of the Monte Carlo method to radiation trans-
port problems appears to be a study by Spencer [9] of the effects of polariza-
tion on multiple (successive) Compton scatterings, a problem further
examined later on using the Boltzmann equation and diffusion theory [10].
This is further emphasized in the pioneering work of Kahn [11] for predicting
neutron fluxes associated with nuclear reactors designed and brought on line
in the 1940s. Mayer used one of the first automatic computers, the ENIAC,
as early as 1949 for neutron transport Monte Carlo calculations as reported
at a 1949 symposium on the Monte Carlo method [12].
the reported average behaviour of the system improves, meaning that the
statistical uncertainty decreases. Almost any complex system in principle
can be modelled; perhaps we wish to model the number of cars passing a
particular intersection during certain times of the day, to optimize traffic
management, or the number of people that will make transactions in a
bank, to evaluate the advantages of different queuing systems. If we know
from experience the distribution of events that occur in our system, we can
generate a pdf and sample it randomly to simulate the real system. Detailed
descriptions of the general principles of the Monte Carlo method are given in
a number of publications [13, 14] and will not be repeated here. The major
components of a Monte Carlo method are briefly described below. These
components comprise the foundation of most Monte Carlo applications.
The following sections will explore them in more detail. An understanding
of these major components will provide a sound foundation for the developer
to construct his own Monte Carlo method. The primary components of a
Monte Carlo simulation method include the following:
(i) Probability density functions (pdfs): the physical system must be
described by a set of pdfs.
(ii) Random number generator: a source of random numbers uniformly
distributed on the unit interval must be available.
(iii) Sampling rule: a prescription for sampling from the specified pdfs.
(iv) Scoring: the outcomes must be accumulated into overall tallies or scores
for the quantities of interest.
(v) Error estimation: an estimate of the statistical error (variance) as a func-
tion of the number of trials and other quantities must be determined.
(vi) Variance reduction techniques: methods for reducing the variance in the
estimated solution to reduce the computational time for Monte Carlo
simulation.
(vii) Parallelization and vectorization algorithms to allow Monte
Carlo methods to be implemented efficiently on advanced computer
architectures.
types of interest must be supplied for each material present. The model also
consists of algorithms used to compute the result of interactions (changes in
particle energy, direction, etc.) based on the physical principles that describe
the interaction of radiation with matter and the cross section data provided.
Therefore, it is extremely important to use an accurate transport model, as
the Monte Carlo result is only as valid as the data supplied.
where is the density of the material, wi the fraction by weight of the ith atomic
constituent as specified in ICRU report 44 [43], and ð=Þi the mass attenuation
coefficients. Different photon cross section libraries show quite large variations
as compared to the most recent EPDL97 data files (figure 1.2). The EPDL97
10 The Monte Carlo method: theory and computational issues
library is the most up-to-date, complete and consistent library available and is
already a standard in the nuclear reactor industry. It is recommended that
Monte Carlo developers only use the most recent version of this library. Its
use as a standard in the simulation of medical imaging systems will help to elim-
inate potential differences between the results obtained with different codes.
Further evaluation of the effect of the photon cross section library on
actual simulation of PET data was performed by generating data sets for the
three-dimensional Hoffman brain phantom contained in a cylindrical homo-
geneous water phantom (20 cm diameter, 17.5 cm height) with apparent relative
concentrations of 4, 1 and 0 for grey matter, white matter and the ventricles,
respectively, simulating the activity distributions found in normal human
brain for cerebral blood flow and metabolism studies currently employed in
PET [44]. Comparisons between the profiles through a simulated two-dimen-
sional projection of the three-dimensional brain phantom estimated using
different libraries are illustrated in figure 1.3. The statistical analysis of the
data sets revealed the existence of a statistically significant difference between
the resulting projections generated using different libraries ( p < 0:025) [42].
electrons are generally assumed to deposit all their energy at the point of
interaction and their bremsstrahlung production is negligible. Thus, electron
transport has not received particular attention in nuclear imaging appli-
cations of the Monte Carlo method. However, a number of investigators
in the field have considered the simulation of electron transport, mainly
for dosimetry calculations [49–53], and for this reason a brief description
of electron transport techniques is given in this section.
Strictly, all the different electron interaction types should be considered
in detail in a Monte Carlo simulation, much in the same way as it is done with
photon interactions. This is a technique called microscopic simulation [54],
which is used at low electron energies, for example in electron microscopy.
The complexity of the procedures used in microscopic techniques varies
considerably, although a common approach is to neglect bremsstrahlung
interactions due to the low energies involved. The most common simple
models used are based on the simulation of all the single-scattering events,
calculating the electron step-length between consecutive collisions with the
elastic mean-free-path. Energy losses are determined from the Bethe theory
of stopping power [55] and energy-loss straggling is accounted for approxi-
mately. These models have been improved by taking inelastic collisions
into account.
When the electron energy increases above a few tens of keV, the large
number of interactions that occur during electron slowing-down makes it
unrealistic to simulate all the single physical interactions. Instead, multiple
collisions are taken into account in what resembles a series of ‘snapshots’ of
an electron history taken at certain time or energy intervals; single inter-
actions of different kinds are considered between two consecutive
snapshots. This grouping scheme constitutes the macroscopic [54] or
condensed-history techniques [56]. In a seminal work, Berger [56] estab-
lished the basis of these techniques for charged-particle transport, where
algorithms are based on multiple scattering theory (Moliere, Goudsmit–
Saunderson) and stopping power calculations to describe the grouped
angular deflections and energy losses, respectively. He divided algorithms
into two broad classes, class I and class II, which are distinguished
mainly by the way in which energy losses are grouped and, more specifi-
cally, by how the individual interactions leading to a large energy transfer,
namely the production of energetic knock-on electrons and bremsstrahlung
photons, are described.
The Monte Carlo system ITS3 [57] with its physics based on the original
ETRAN code by Berger and Seltzer [58], and EGS4 [33] are, respectively, the
most widely used examples of class I and class II algorithms in medical
physics. Unbiased comparisons of the two systems for a variety of situations
and energies have been made by Andreo [59]. In recent years, the most
advanced electron transport algorithms, with physics that has overcome
that of older Monte Carlo systems, combine macroscopic and microscopic
Electron transport 15
techniques, where the latter are used in situations that require sophisticated
boundary-crossing algorithms [60, 61].
data, at the level of 1% or so, associated with the computer revolution that
has improved so dramatically many other aspects of medical radiation
physics.
A direct Monte Carlo simulation using true probability functions may require
an unacceptably long time to produce statistically relevant results. Emission of
photons is isotropic, so directional parameters may be sampled uniformly
within their individual ranges. Nuclear imaging systems have a low geometri-
cal efficiency because of the small solid angle defined by the collimator and/or
the small axial aperture. Therefore, the calculation would be very ineffective
in terms of required computing time [72]. It is thus desirable to bias the
sampling (non-analogue sampling) by introducing different types of impor-
tance sampling and other variance reduction techniques to improve the
computational efficiency of the Monte Carlo method [73]. The results obtained
by non-analogue simulation are, however, biased by the variance reduction
technique and a correction for this is required. A particle history weight, W,
is introduced, which describes the probability of the particle following the
current path. This weight is calculated for each particle history, and used
in the calculation of the results. If an event occurs, the weight W is added
to the counter rather than incrementing the counter by one unit. Bielajew
and Rogers [71] divided variance reduction techniques into three categories:
those that concern photon transport only, those that concern electron trans-
port only, and other more general methods. The most useful techniques are
described below. It is worth noting that concerning the Poisson nature of the
activity distributions in nuclear imaging, these variance reduction approxi-
mations may result in statistically relevant deviations from an otherwise
unbiased simulation result since Monte Carlo simulated projection data
are (as a result of variance reduction) not count but probability data.
Since the maximum distance, dmax , the photon travels before interaction is
infinite and the number of mean free paths across the geometry in any
practical situation is finite, there is a large probability that photons leave
the geometry of interest without interacting. To increase the statistical
accuracy in the imparted energy calculation, we force the photons to interact
by assigning dmax a finite distance, e.g., the thickness of the detector being
simulated. A true distributed photon pathlength d within dmax can be
sampled from the equation
1
d ¼ lnð1 R½1 e dmax Þ: ð1:13Þ
The photon’s weight must be multiplied by the interaction probability
Wn þ 1 ¼ Wn ½1 e dmax : ð1:14Þ
In emission computed tomography, the photon is allowed to interact through
coherent or incoherent interactions only within the phantom since photo-
absorption does not contribute to energy imparted in the crystal. The
weight is then multiplied by the probability of the photon being scattered:
þ coh
Wn þ 1 ¼ Wn incoh ; ð1:15Þ
where incoh and coh are the cross section data for incoherent and coherent
scattering, respectively, and is the total linear attenuation coefficient.
where C is a parameter that adjusts the magnitude of the scaling and the
angle of the photon with respect to the direction of interest. The new weight-
ing factor is given by
expðd C cos Þ
Wn þ 1 ¼ W n : ð1:19Þ
ð1 C cos Þ
Note that the new weighting factor is dependent on d . If 0 < C < 1, the
particle pathlength is stretched in the forward direction, which is used for
shielding problems. For 1 < C < 0, the average distance to the next
interaction is shortened in the forward direction, which is used for surface
problems. For C ¼ 0, we recover the unbiased sampling. The optimal
choice of this parameter is dependent on the problem to be solved. The
general guideline is to avoid to use large weighting factors because they
may increase the variance.
Russian roulette and splitting are often used together with the
exponential transform, although they are still effective when used indepen-
dently. In Russian roulette, a random number is selected and compared
with a threshold, . If the random number turns out to be smaller than ,
the particle is allowed to survive but the weight should be updated accord-
ingly: Wn þ 1 ¼ Wn =. In particle splitting, a particle coming from a region
of interest can be divided into N particles, each having a new weighting:
Wn þ 1 ¼ Wn =N.
Use of geometry symmetry. The use of some of the inherent symmetry of the
geometry may realize considerable increase in efficiency. If both the source
and target configurations contain cylindrical planar or spherical–conical
simulation geometries, the use of symmetries is more obvious. Other uses
of symmetry are less obvious but the saving in computing time is worth
the extra care and coding.
1.7. SUMMARY
ACKNOWLEDGMENTS
The author would like to thank Professor Pedro Andreo and Dr Mike Stabin
for fruitful discussions. Their comments and suggestions on this chapter are
gratefully acknowledged.
REFERENCES
[1] Metropolis N and Ulam S 1949 The Monte Carlo method J. Am. Stat. Assoc. 44 335–
41
[2] McCracken D D 1955 The Monte Carlo method Sci. Am. 192 90–6
[3] Van Dam C 2000 The Urim and Thummim: A Means of Revelation in Ancient Israel ed
I Eisenbrauns (Winona Lake, IN)
[4] de Buffon G 1777 Essai d’arithmétique morale Supple´ment à l’Histoire Naturelle
vol 4
[5] Rayleigh J W S 1899 On the theory of resonance Collected Scientific Papers 1
33–75
[6] Kolmogorov A N 1933, 1956 Foundations of the Theory of Probability (2nd English
edition) transl ed N Morrison (New York: Chelsea)
[7] Ulam S M 1976 Adventures of a Mathematician (New York: Scribners)
[8] Ulam S M 1951, 1949 On the Monte Carlo Method Symposium on Large-Scale Digital
Calculating Machines (Cambridge, MA: Harvard University Press)
[9] L V Spencer 1948 ‘Polarization of multiply scattered gamma rays’ (unpublished)
(National Bureau of Standards Interim Report on Contract with US Office of
Naval Research)
[10] Spencer L V and Wolff C 1953 Penetration and diffusion of hard X-rays: Polarization
effects Phys. Rev. 90 510–14
[11] Kahn H 1949 ‘Stochastic (Monte Carlo) attenuation analysis’ Report p-88
(RAND)
[12] Mayer M 1951 Report on a Monte Carlo calculation performed with the ENIAC
Monte Carlo Method eds A S Householder, G E Forsythe and H H Germond US
National Bureau of Standards Applied Math Series No 12, pp 19–20
[13] Raeside D E 1976 Monte Carlo principles and applications Phys. Med. Biol. 21
18197
[14] Turner J E, Wright H A and Hamm R N 1985 A Monte Carlo primer for health
physicists Health Phys. 48 717–33
[15] Marsaglia G and Zaman A 1994 Some portable very-long-period random number
generators Comput. Phys. 8 117–21
24 The Monte Carlo method: theory and computational issues
[16] Ferrenberg A M, Landau D P and Wong Y J 1992 Monte Carlo simulations: hidden
errors from ‘good’ random number generators Phys. Rev. Lett. 69 3382–4
[17] Vattulainen I, Kankaala K, Saarinen J and Ala-Nissila T 1995 A comparative
study of some pseudorandom number generators Comp. Phys. Comm. 86 209–
26
[18] Booth T E 1998 The intelligent random number technique in MCNP Nucl. Sci. Eng.
100 248–54
[19] Lehmer D H 1949 Mathematical methods in large-scale computing units 2nd
Symposium on Large-Scale Digital Calculating Machinery, Cambridge (Cambridge,
MA: Harvard University Press) pp 141–6
[20] Marsaglia G and Zaman A 1993 Monkey tests for random number generators Comp.
Math. Appl. 23 1–10
[21] Knuth D E 1981 The Art of Computer Programming, Seminumerical Algorithms 2nd
edition (Reading, MA: Addison-Wesley) vol 2
[22] Press W H 1992 Numerical Recipes in C, The Art of Scientific Computing (Cambridge:
Cambridge University Press)
[23] Anderson S L 1990 Random number generators on vector supercomputers and other
advanced architectures SIAM Rev. 32 221–51
[24] Mascagni M, Cuccaro S A, Pryor D V and M L Robinson 1995 A fast, high-quality,
and reproducible lagged-Fibonacci pseudorandom number generator J. Comput.
Phys. 15 211–19
[25] Brent R P 1994 On the periods of generalized Fibonacci recurrences Math. Comput.
63 389-401
[26] Zaidi H, Labbé C and Morel C 1998 Implementation of an environment for Monte
Carlo simulation of fully 3D positron tomography on a high-performance parallel
platform Parallel Comput. 24 1523–36
[27] Zaidi H 1999 Relevance of accurate Monte Carlo modeling in nuclear medical
imaging Med. Phys. 26 574–608
[28] Hubbell J H 1999 Review of photon interaction cross section data in the medical and
biological context Phys. Med. Biol. 44 R1–22
[29] Storm E and Israel H I 1970 Photon cross sections from 1 keV to 100 MeV for
elements Z ¼ 1 to Z ¼ 100 Nucl. Data Tables A7 565–681
[30] Briesmeister J F 1997 (Los Alamos, NM: Los Alamos National Laboratory)
[31] Hubbell J H et al. 1975 Atomic form factors, incoherent scattering functions, and
photon scattering cross sections J. Phys. Chem. Ref. Data 4 471–538
[32] Cullen D E et al. 1989 ‘Tables and graphs of photon interaction cross sections from
10 eV to 100 GeV derived from the LLNL Evaluated Photon Data Library (EPDL)’
UCRL-50400 (Lawrence Livermore National Laboratory, CA)
[33] Nelson W R, Hirayama H and Rogers D W O 1985 (Stanford, CA: Stanford Linear
Accelerator Center)
[34] Trubey D K, Berger M J and Hubbell J H 1989 Photon cross sections for ENDF/B-
IV. American Nuclear Society Topical Meeting, Advances in Nuclear Engineering
Computation and Radiation Shielding, Santa Fe, NM
[35] Sakamoto Y 1993 Photon cross section data PHOTX for PEGS4 code. EGS4 users’
meeting, Japan (KEK Proceedings) pp 77–82
[36] Halbleib J A, Kensek R P, Valdez G D, Seltzer S M and Berger M J 1992 ITS: The
Integrated TIGER Series of electron/photon transport codes—version 30 IEEE
Trans. Nucl. Sci. 39 1025–30
References 25
[37] Boone J M and Chavez A E 1996 Comparison of x-ray cross sections for diagnostic
and therapeutic medical physics Med. Phys. 23 1997–2005
[38] Berger M J and Hubbell J H 1987 ‘XCOM: Photon cross sections on a personal
computer’ NBSIR 87-3597 (US National Bureau of Standards)
[39] Brun R, Bruyant F, Maire M, McPherson A C and Zanarini P 1994 ‘GEANT detec-
tor description and simulation tool’ DD/EE/84-1 (CERN)
[40] Picard Y, Thompson C J and Marrett S 1993 Improving the precision and accuracy of
Monte Carlo simulation in positron emission tomography IEEE Trans. Nucl. Sci. 39
1111–16
[41] Cullen D E, Hubbell J H and Kissel L 1997 ‘EPDL97: the Evaluated Photon Data
Library’ UCRL-50400 (Lawrence Livermore National Laboratory, CA)
[42] Zaidi H 2000 Comparative evaluation of photon cross section libraries for materials
of interest in PET Monte Carlo simulations IEEE Trans. Nucl. Sci. 47 2722–35
[43] ICRU 1989 Tissue Substitutes in Radiation Dosimetry and Measurement vol 44
(International Commission on Radiological Units and Measurements)
[44] Hoffman E J, Cutler P D, Digby W M and Mazziotta J C 1990 3-D phantom to
simulate cerebral blood flow and metabolic images for PET IEEE Trans. Nucl. Sci.
37 616–20
[45] Klein O and Nishina Y 1929 Über die Streuung von Strahlung durch freie Elektronen
nach der neuen relativistischen Quantendynamik von Dirac Z. Phys. 52 853–68
[46] Kahn H 1956 Use of Different Monte Carlo Sampling Techniques in Monte Carlo
Methods ed H A Meyer (New York: Wiley)
[47] Kaplan M S, Harrison R L and Vannoy S D 1998 Coherent scatter implementation
for SimSET IEEE Trans. Nucl. Sci. 45 3064–8
[48] Zaidi H, Scheurer A H and Morel C 1999 An object-oriented Monte Carlo simulator
for 3D cylindrical positron tomographs Comput. Methods Programs Biomed. 58
13345
[49] Stabin M and Konijnenberg M W 2000 Re-evaluation of absorbed fractions for
photons and electrons in spheres of various sizes J. Nucl. Med. 41 149–60
[50] Bouchet L G and Bolch W E 1999 A three-dimensional transport model for deter-
mining absorbed fractions of energy for electrons within cortical bone J. Nucl.
Med. 40 2115–24
[51] Bouchet L G, Jokisch D W and Bolch W E 1999 A three-dimensional transport model
for determining absorbed fractions of energy for electrons within trabecular bone
J. Nucl. Med. 40 1947–66
[52] Furhang E E, Chui C S, Kolbert K S, Larson S M and Sgouros G 1997 Implementa-
tion of a Monte Carlo dosimetry method for patient-specific internal emitter therapy
Med. Phys. 24 1163–72
[53] Yoriyaz H, Stabin M and dos Santos A 2001 Monte Carlo MCNP-4B-based
absorbed dose distribution estimates for patient-specific dosimetry J. Nucl. Med. 42
662–9
[54] Andreo A 1991 Monte Carlo techniques in medical radiation physics Phys. Med. Biol.
36 861–920
[55] Bethe H A 1930 Zur Theorie des Durchgangs schneller Korpuskularstrahlen durch
Materie Ann. Phys., Lpz. 5 325
[56] Berger M J 1963 Monte Carlo calculation of the penetration of diffusion of fast
charged particles Methods in Computational Physics ed B Alder, S Feinbach and
M Rotenberg (New York: Academic) pp 135–215
26 The Monte Carlo method: theory and computational issues
[57] J A Halbleib et al. 1992 ‘ITS version 3.0: The integrated TIGER Series of coupled
electron/photon Monte Carlo transport codes’ Report SAND91-1634 (Sandia
National Laboratories)
[58] Berger M J and Seltzer S M 1968 ‘Electron and photon transport programs: I. Intro-
duction and notes on program DATAPAC 4; II. Notes on program ETRAN-15’
Reports NBS 9836 and 9837. Report ORNL-RSIC CCC-107 (Oak Ridge: Oak
Ridge National Laboratory; National Bureau of Standards)
[59] Andreo P 1995 A Comparison Between the ITS3 and EGS4 Monte Carlo Codes in
Radiation Physics, Advanced Monte Carlo Computer Programs for Radiation Trans-
port (Saclay: Publications Service OECD)
[60] Salvat F, Fernandez-Varea J M, Sempau J and Mazurier J 1999 Practical aspects of
Monte Carlo simulation of charged particle transport: mixed algorithms and variance
reduction techniques Radiat. Environ. Biophys. 38 15–22
[61] Kawrakow I 2000 Accurate condensed history Monte Carlo simulation of electron
transport. I. EGSnrc, the new EGS4 version Med. Phys. 27 485–98
[62] ICRU 1984 ‘Stopping powers for electrons and positrons’ Report 37 (International
Commission on Radiological Units and Measurements)
[63] ICRU 1993 ‘Stopping powers and ranges for protons and alpha particles’ Report 49
(International Commission on Radiological Units and Measurements)
[64] Berger M J 1993 ‘ESTAR, PSTAR, and ASTAR: computer programs for calculating
stopping-power and range tables for electrons, protons, and helium ions’ NIST
Report NISTIR 4999 (National Institute of Standards and Technology)
[65] Berger M J, Seltzer S M, Domen S R and Lamperti P J 1975 Stopping-power ratios
for electron dosimetry with ionization chambers. Biomedical Dosimetry (Vienna:
IAEA) pp 589–609
[66] Nahum A E 1978 Water/air mass stopping power ratios for megavoltage photon and
electron beams Phys. Med. Biol. 23 24–38
[67] Andreo P and Brahme A 1986 Stopping power data for high-energy photon beams
Phys. Med. Biol. 31 839–58
[68] Ding G X, Rogers D W and Mackie T R 1995 Calculation of stopping-power ratios
using realistic clinical electron beams Med. Phys. 22 489–501
[69] Malamut C, Rogers D W and Bielajew A F 1991 Calculation of water/air stopping-
power ratios using EGS4 with explicit treatment of electron–positron differences
Med. Phys. 18 1222–8
[70] ICRU 1984 ‘Radiation dosimetry: electron beams with energies between 1 and 50 MeV’
Report 35 (International Commission on Radiological Units and Measurements)
[71] Bielajew A F and Rogers D W O 1989 Variance-reduction techniques Monte Carlo
Transport of Electrons and Photons ed T Jenkins, A Rindi, A Nahum and D W O
Rogers (New York: Plenum) pp 407–19
[72] Haynor D R, Harrison R L and Lewellen T K 1990 Improving the efficiency of
emission tomography using variance reduction techniques IEEE Trans. Nucl. Sci.
37 749–53
[73] Zubal I G and Harrell C R 1992 Voxel-based Monte Carlo calculations of nuclear
medicine images and applied variance reduction techniques Image Vision Computing
10 342–8
[74] Haynor D R, Harrison R L and Lewellen T K 1991 The use of importance sampling
techniques to improve the efficiency of photon tracking in emission tomography
simulations Med. Phys. 18 990–1001
References 27
[75] Bielajew A F and Rogers D W O 1987 PRESTA: The parameter reduced electron-step
transport algorithm for electron Monte Carlo transport Nucl. Instrum. Methods Phys.
Res. B18 165–81
[76] Ma C M and Nahum A E 1993 Calculation of absorbed dose ratios using correlated
Monte Carlo sampling Med. Phys. 20 1189–99
[77] Holmes M A et al. 1993 The application of correlated sampling to the computation of
electron beam dose distributions in heterogeneous phantoms using the Monte Carlo
method Phys. Med. Biol. 38 675–88
Chapter 2
2.1. INTRODUCTION
The use of the Monte Carlo method to simulate radiation transport has
become the most accurate means of predicting absorbed dose distributions
and other quantities of interest in the radiation treatment of cancer patients
using either external or radionuclide radiotherapy. The same trend has
occurred for the estimation of the absorbed dose in diagnostic procedures
using radionuclides. There is broad consensus in accepting that the earliest
Monte Carlo calculations in medical radiation physics have been made in
the area of nuclear medicine, where the technique has been used for
dosimetry modelling and computations. Formalism and data based on
Monte Carlo calculations, developed by the Medical Internal Radiation
Dose (MIRD) Committee of the Society of Nuclear Medicine, have been
published in a series of supplements to the Journal of Nuclear Medicine,
the first one being released in 1968 [1]. Some of these pamphlets made
extensive use of Monte Carlo calculations to derive specific absorbed frac-
tions for electron and photon sources uniformly distributed in organs of
mathematical phantoms. Interest in Monte Carlo-based dose calculations
with and emitters has been revived with the application of radiolabelled
monoclonal antibodies to radioimmunotherapy (RIT). As a consequence of
this generalized use, many questions are being raised, primarily about the
need and potential of Monte Carlo techniques, but also about how accurate
it really is, and what it would take to apply it clinically and make it available
widely to the nuclear medicine community at large.
In this chapter, we try to answer many of these questions, sometimes
without explicitly stating them. Many other questions will be answered
when Monte Carlo techniques are implemented and used for more routine
calculations and for in-depth investigations. A survey is presented of the
28
Monte Carlo techniques in medical radiation physics 29
Monte Carlo techniques have become one of the most popular tools in
different areas of medical physics following the development and subsequent
implementation of powerful computing systems for clinical use. In particu-
lar, they have been extensively applied to simulate processes involving
random behaviour and to quantify physical parameters that are difficult or
even impossible to calculate analytically or to determine by experimental
measurements.
Following the review by Andreo [2], the applications of the Monte Carlo
method in medical radiation physics cover almost all topics including radia-
tion protection, nuclear medicine, diagnostic radiology and radiotherapy,
with an increasing enthusiastic interest in exotic and exciting new applica-
tions such as intravascular radiation therapy and boron neutron capture
therapy and synovectomy. The use of Monte Carlo techniques in the last
two areas falls outside the scope of this survey whereas applications in
radiation protection and nuclear medicine imaging and dosimetry are
covered in the following sections. The present authors are not aware of a
comprehensive review on applications of the Monte Carlo method in diag-
nostic radiology and would recommend a recent paper by Boone et al. [6]
for an up-to-date validation study providing a general overview of the
field. Interested readers in radiotherapy physics are encouraged to check
the reviews by Ma and Jiang [7] for electron beams, and by Ahnesjo and
Aspradakis [8] for photon beams. With the rapid development of computer
technology, Monte Carlo-based treatment planning for radiation therapy is
becoming practical.
An area which has received considerable attention during the past
few years is intravascular brachytherapy (see chapter 14). The irradiation
of the coronary arteries is being studied by a number of researchers world-
wide, with the hope of reducing the rate of restenosis. Currently, there is
30 Monte Carlo techniques in nuclear medicine dosimetry
capture synovectomy (GNCS), has also been conducted using Monte Carlo
simulations where 10 B and 157 Gd have been compared as isotopes for
accelerator-based NCS [14].
Figure 2.1. Scientific and technical strategy for recording accurate functional
images. In bold, the parts where Monte Carlo simulation plays an important
role. Reprinted with permission from AAPM [3].
SPECT [17, 18], PET [19, 20] and MET [21]. However, due to computer
limitations, the method has not yet fully lived up to its potential. With the
advent of high-speed supercomputers, the field has received increased attention,
particularly with parallel algorithms, which have much higher execution rates.
Applications of Monte Carlo techniques 33
Smith et al. [38] used Monte Carlo modelling to study photon detection
kernels, which characterize the probabilities that photons emitted by radio-
isotopes in different parts of the source region will be detected at particular
projection pixels of the projection images for the case of parallel-hole
collimators. The authors also proposed a reconstruction method using the
three-dimensional kernels in which projection measurements in three
adjacent planes are used simultaneously to estimate the source activity of
the centre plane. The search for unified reconstruction algorithms led to
the development of inverse Monte Carlo (IMC) reconstruction techniques.
The principal merits of IMC are that, like direct Monte Carlo, the method
can be applied to complex and multivariable problems, and variance reduc-
tion procedures can be applied [39]. Floyd et al. [40] used IMC to perform
tomographic reconstruction for SPECT with simultaneous compensation
for attenuation, scatter, and distance-dependent collimator resolution.
The Monte Carlo code Eidolon [20], for example, has been developed
using modern software development tools. This simulator is an efficient
tool that can be used to generate data sets in a controllable manner in
order to assess different reconstruction algorithms and scatter correction
methods [37]. As the ‘best’ algorithm can only be selected with respect to a
certain task, different ‘basic’ performance measures can be used. Image
degrading effects are illustrated using simulated projections of the digitized
three-dimensional Hoffman brain phantom [41]. A slice of this phantom is
shown in figure 2.4(A). The ratio between the activity in white, grey matter
and ventricles has been chosen as 1 :4 :0, respectively. The projections of
this phantom at different levels of fidelity are generated. The strengths of
the image degrading factors are characteristic of an [18 F]-FDG-brain
study. Figure 2.4(F–I) shows the effects of different aspects of image degra-
dation on filtered backprojection reconstructions. The loss of resolution
caused by detector blurring (FWHM ¼ 4 mm) on projection data and FBP
reconstructions is shown in figures 2.4(C) and 2.4(G), respectively, while in
D and H, effects of detector blurring, attenuation and scatter are included
in the simulation and no corrections performed on the simulated data sets.
Finally, in E and I, effects of detector blurring, attenuation and scatter are
included and appropriate corrections for attenuation and scatter applied.
using both the EGS4 and MCNP4B Monte Carlo codes [56]. Calculation of
absorbed fractions for positron emitters relevant to neurologic studies has
also been reported [57].
The application of the Monte Carlo method to internal radiation dosi-
metry is further emphasized in two recent MIRD Pamphlets. In MIRD
Pamphlet No. 15 [58], the EGS4 Monte Carlo radiation transport code has
been used to revise substantially the dosimetric model of the adult head
and brain originally published in MIRD Pamphlet No. 5. Pamphlet No.
17 [59] crystallizes the utility of the MIRD formalism for the calculation of
the non-uniform distribution of radiation-absorbed dose in different
organs through the use of radionuclide-specific S values defined at the
voxel level. Skeletal S values and absorbed fractions estimates to both
marrow and endosteum in trabecular and cortical bone have been further
improved through the use of electron transport models [60].
Mathematical anthropomorphic phantoms are continuously being
improved. Current developments are aimed at computer phantoms that
are flexible while providing accurate modelling of patient populations. An
important contribution came from Bouchet and Bolch [61], making available
a series of five dosimetric head and brain models developed to allow more
precise dosimetry in paediatric neuroimaging procedures. A new rectal
model [62] and dynamic urinary bladder model [63] have also been proposed
more recently. To develop more patient-specific dosimetry, new mathemati-
cal models for adults of different height have been developed using anthropo-
metric data [64]. The use of dynamic anthropomorphic phantoms in Monte
Carlo simulations is becoming possible due to the increasing availability of
computing power. This includes the development of appropriate primitives
that allow accurate modelling of anatomical variations and patient
motion, like superquadrics [65] and non-uniform rational B-spline surfaces
[66]. More recently, an image-based whole-body model, called VIP-Man,
has been developed using high-resolution transversal colour photographic
images obtained from the National Library of Medicine’s Visible Human
Project [67]. The phantom is shown in figures 2.5(b) and 2.5(c) in comparison
with the conventional anthropomorphic mathematical models (figure 2.5(a)).
The EGS4 code has been also used to estimate specific absorbed fractions
from internal electron emitters with energies from 100 keV to 4 MeV [68].
Figure 2.5. Anterior view of (a) the mathematical model and the VIP-Man
model in (b) two-dimensions and (c) three-dimensions. Reprinted with
permission from the Institute of Physics [68]. (See plate 1 for colour version.)
point sources in water of relevance to RIT [78]. This study provided the
initial data required for modelling and analysing the scatter, attenuation,
and image formation processes in quantitative imaging of bremsstrahlung
for RIT dosimetry. Furhang et al. [70] generated photon point dose kernels
and absorbed fractions in water for the full photon emission spectrum of
radionuclides of interest in nuclear medicine, by simulating the transport
of particles using Monte Carlo techniques. The kernels have been then
fitted to a mathematical expression.
A unified approach for photon and particle dosimetry has been
proposed by Leichner [79] by fitting Berger’s tables for photons and electrons
to generate an empirical function that is valid for both photons and
particles. Therefore both point-kernel and Monte Carlo techniques can be
effectively employed to calculate absorbed dose to tissue from radionuclides
that emit photons or electrons. The latter are computationally much more
intensive; however, point-kernel methods are restricted to homogeneous
tissue regions that can be mathematically described by analytical geometries
whereas Monte Carlo methods have the advantage of being able to accom-
modate heterogeneous tissue regions with complex geometric shapes.
parameter standard deviations, due to data noise, and much smaller noise-
induced biases. The authors reported good agreement between regression
and traditional methods.
Welch [84] investigated and quantified the effect of typical SPECT
system resolution and photon counting statistics on the bias and precision
of dynamic cardiac SPECT parameters. Simulation of dynamic SPECT
projection data has been performed using a realistic human torso phantom
assuming both perfect system resolution and a system resolution typical of
a clinical SPECT system. The results showed that the rate constant charac-
terizing the washing of activity into the myocardium is more sensitive to
the position of the region of interest than the washout rate constant and
that the main effect of increased photon noise in the projection data is to
decrease the precision of the estimated parameters.
Computer simulations demonstrate that estimation of the kinetic
parameters directly from the projections is more accurate than the estimation
from the reconstructed images [85]. A strategy for joint estimation of
physiological parameters and myocardial boundaries has been proposed
and evaluated by simulated myocardial perfusion studies based on a simpli-
fied heart model [86]. A method allowing the estimation of kinetic parameters
directly from SPECT cone-beam projections has also been proposed and
validated with a simulated chest phantom [87]. The results showed that
myocardial uptake and washout parameters estimated by conventional
analysis of noiseless simulated cone-beam data had biases ranging between
3 and 26% and between 0 and 28%, respectively, while uncertainties of
parameter estimates with this method ranged between 0.2 and 9% for the
uptake parameters and between 0.3 and 6% for the washout parameters.
Figure 2.6. Diagram showing the essential steps required in developing a three-
dimensional internal dosimetry program for treatment planning with RIT based
on quantitative emission computed tomography where Monte Carlo simulations
play a crucial role. Adapted and reproduced with permission from AAPM [3].
46 Monte Carlo techniques in nuclear medicine dosimetry
treated with this modality is close to 30 000 [118]. From the two properties
mentioned above the physical dose distribution of protons, enabling accurate
dose conformation delivery, poses a clear advantage over conventional
radiotherapy beams; however, their radiobiological effect is only about
810% higher than that of photons or electrons. The radiobiological
superiority of ions heavier than protons, with an RBE augmentation of the
order of three to four times higher than that of conventional beams, led to
an increased interest worldwide that resulted in the construction of the
Japanese HIMAC clinical facility in Chiba, near Tokyo. HIMAC started
the treatment of patients mainly with carbon ions in 1994 [119], and more
than 1000 patients have been treated in this facility. It had become clear
that whereas heavy ions like the neon beams used in the early 1970s in
Berkeley have the largest radiobiological effects, these also appear in regions
close to the beam entrance and in the plateau region, where usually normal
tissue is situated. In addition, due to the large penetration of the fragmenta-
tion products released by the incident ions, the tail of the dose distribution
beyond the Bragg peak may be too high for sparing normal tissue beyond
the primary ion range. These two aspects suggest that the ideal ions for
radiotherapy are the light ions, and for this reason carbon has dominated
the clinical applications at HIMAC. The GSI heavy-ion physics research
facility in Darmstadt, Germany, initiated clinical treatments also with
carbon ions in 1997 [115]. GSI uses advanced three-dimensional beam
scanning techniques and PET imaging to monitor the delivery of dose by
visualizing the positrons emitted by 10 C and 11 C ions produced in the nuclear
reactions of the incident carbon ions with tissue, together with radio-
biologically optimized treatment planning. Plans for building additional
radiotherapy light-ion facilities have followed in Europe, and feasibility studies
have been conducted for such facilities in Germany, Austria, Italy, France,
Sweden and Switzerland, among others [120].
Several groups have reported the applicability of PET to in vivo dosi-
metry for proton radiotherapy using the Monte Carlo method. Different
Monte Carlo codes have been used to investigate this challenging field
including Berger’s PTRAN [121], used by Del Guerra et al. [122], FLUKA
[123], and specially designed software [114]. Parodi and Enghardt [123], in
particular, made Monte Carlo simulations of the number and spatial distri-
bution of the positron emitters created by therapeutic protons in PMMA.
Since proton therapy requires a particle fluence up to 20 times higher than
that required in carbon therapy for the same physical dose level, the authors
concluded that the total activity produced within the irradiated volume
should be twice as intense as that actually used at the carbon ion facility at
GSI. Moreover, the considered ratio between the fluence of protons and
carbon ions could increase considerably due to the higher relative biological
effectiveness of carbon ions, leading to a lower physical dose than for protons
at the same level of biological dose.
48 Monte Carlo techniques in nuclear medicine dosimetry
REFERENCES
[1] Berger M J 1968 MIRD pamphlet 2: energy deposition in water by photons from
point isotropic sources J. Nucl. Med. 9 15–25
[2] Andreo A 1991 Monte Carlo techniques in medical radiation physics Phys. Med.
Biol. 36 861–920
[3] Zaidi H 1999 Relevance of accurate Monte Carlo modeling in nuclear medical
imaging Med. Phys. 26 574–608
[4] Jenkins T M, Nelson W R and Rindi A 1998 Monte Carlo Transport of Electrons and
Photons (New York: Plenum)
[5] Ljungberg M, Strand S-E and King M A 1998 Monte Carlo Calculations in Nuclear
Medicine: Applications in Diagnostic Imaging (London: IOP Publishing)
[6] Boone J M, Buonocore M H and Cooper V N 2000 Monte Carlo validation in
diagnostic radiological imaging Med. Phys. 27 1294–304
[7] Ma C M and Jiang S B 1999 Monte Carlo modelling of electron beams from medical
accelerators Phys. Med. Biol. 44 R157–89
[8] Ahnesjo A and Aspradakis M M 1999 Dose calculations for external photon beams
in radiotherapy Phys. Med. Biol. 44 R99–155
[9] Stabin M , Konijnenberg M, Knapp F F and Spencer R H 2000 Monte Carlo
modeling of radiation dose distributions in intravascular radiation therapy Med.
Phys. 27 1086–92
[10] Reynaert N, Verhaegen F, Taeymans Y, Van Eijkeren M and Thierens H 1999
Monte Carlo calculations of dose distributions around 32P and 198Au stents for
intravascular brachytherapy Med. Phys. 26 1484–91
[11] Lee C L, Zhou X L, Kudchadker R J, Harmon F and Harker Y D 2000 A Monte
Carlo dosimetry-based evaluation of the 7Li(p,n)7Be reaction near threshold for
accelerator boron neutron capture therapy Med. Phys. 27 192–202
[12] Levin V, Andreo P and Dodd B eds 2001 Current Status of Neutron Capture Therapy.
IAEA Tec-Doc 1223 (Vienna: IAEA, International Atomic Energy Agency)
[13] Gierga D P, Yanch J C and Shefer R E 2000 Development and construction of a
neutron beam line for accelerator-based boron neutron capture synovectomy
Med. Phys. 27 203–14
[14] Gierga D P, Yanch J C and Shefer R E 2000 An investigation of the feasibility of
gadolinium for neutron capture synovectomy Med. Phys. 27 1685–92
[15] Zaidi H 1996 Quantitative SPECT: recent developments in detector response,
attenuation and scatter correction techniques Physica Medica 12 101–17
[16] Webb S, Binnie D M, Flower M A and Ott R J 1992 Monte Carlo modelling of the
performance of a rotating slit-collimator for improved planar gamma-camera
imaging Phys. Med. Biol. 37 1095–108
[17] Ljungberg M and Strand S-E 1989 A Monte Carlo program for the simulation of
scintillation camera characteristics Comput. Meth. Programs Biomed. 29 257–72
[18] de Vries D J et al. 1990 Development and validation of a Monte Carlo simulation of
photon transport in an Auger camera IEEE Trans. Med. Imaging 9 430–8
[19] Dahlbom M, Rosenquist G, Eriksson L and Bohm C 1989 A study of the possibility
of using multi-slice PET systems for 3D imaging IEEE Trans. Nucl. Sci. 36 1066–71
[20] Zaidi H, Scheurer A H and Morel C 1999 An object-oriented Monte Carlo simulator
for 3D cylindrical positron tomographs Comput. Meth. Programs Biomed. 58
13345
References 49
[21] Dahlbom M et al. 1998 A YSO/LSO phoswich array detector for single and
coincidence photon imaging IEEE Trans. Nucl. Sci. 45 1128–32
[22] Derenzo S E 1981 Monte Carlo calculations of the detection efficiency of NaI(Tl),
BGO, CsF, Ge and plastic detectors for 511 keV photons IEEE Trans. Nucl. Sci.
28 131–6
[23] Derenzo S E and Riles J K 1982 Monte Carlo calculations of the optical coupling
between bismuth germanate crystals and photomultiplier tubes IEEE Trans. Nucl.
Sci. 29 191–5
[24] Bottigli U et al. 1985 Monte Carlo simulation and experimental tests on BGO, CsF
and NaI(Tl) crystals for positron emission tomography J. Nucl. Med. Allied Sci. 29
221–7
[25] Bice A N et al. 1990 Monte Carlo simulation of BaF2 detectors used in time-of-flight
positron emission tomography IEEE Trans. Nucl. Sci. 37 696–701
[26] Lopes M I, Chepel V, Carvalho J C, Ferreira Marques R and Policarpo A J P 1995
Performance analysis based on a Monte Carlo simulation of a liquid xenon PET
detector IEEE Trans. Nucl. Sci. 42 2298–302
[27] Bradshaw J, Burnham C and Correia J 1985 Application of Monte Carlo methods to
the design of SPECT detector systems IEEE Trans. Nucl. Sci. 32 753–7
[28] Kimiaei S, Ljungberg M and Larsson S A 1997 Evaluation of optimally designed
planar-concave collimators in single-photon emission tomography Eur. J. Nucl.
Med. 24 1398–404
[29] Pavlopoulos S and Tzanakos G 1996 Design and performance evaluation of a high-
resolution small animal positron tomograph IEEE Trans. Nucl. Sci. 43 3249–55
[30] Moses W W, Virador P R G, Derenzo S E, Huesman R H and Budinger T F 1997
Design of a high-resolution, high-sensitivity PET camera for human brains and
small animals IEEE Trans. Nucl. Sci. 41 1487–91
[31] Moisan C, Tupper P, Rogers J G and de Jong J K 1996 A Monte Carlo study of the
acceptance to scattered events in a depth encoding PET camera IEEE Trans. Nucl.
Sci. 43 1974–80
[32] Ljungberg M, Strand S-E, Rajeevan N and King M A 1994 Monte Carlo simulation
of transmission studies using a planar source with a parallel collimator and a line
source with a fan-beam collimator IEEE Trans. Nucl. Sci. 41 1577–84
[33] Wegmann K et al. 1999 Investigation of the scatter contribution in single photon
transmission measurements by means of Monte Carlo simulations IEEE Trans.
Nucl. Sci. 46 1184–90
[34] Gustafsson A et al. 1998 Evaluation of attenuation corrections using Monte Carlo
simulated lung SPECT Phys. Med. Biol. 43 2325–36
[35] Zaidi H 2001 Scatter modelling and correction strategies in fully 3-D PET Nucl.
Med. Commun. 22 1181–4
[36] El Fakhri G N et al. 1999 Respective roles of scatter, attenuation, depth-dependent
collimator response and finite spatial resolution in cardiac single-photon emission
tomography quantitation: a Monte Carlo study Eur. J. Nucl. Med. 26 437–46
[37] Zaidi H 2000 Comparative evaluation of scatter correction techniques in 3D
positron emission tomography Eur. J. Nucl. Med. 27 1813–26
[38] Smith M F, Floyd C E, Jaszczak R J and Coleman R E 1992 Three-dimensional
photon detection kernels and their application to SPECT reconstruction Phys.
Med. Biol. 37 605–22
[39] Dunn W L 1981 Inverse Monte Carlo analysis J. Comput. Phys. 41 154–66
50 Monte Carlo techniques in nuclear medicine dosimetry
[40] Floyd C E, Jaszczak R J, Greer K L and Coleman R E 1986 Inverse Monte Carlo as
a unified reconstruction algorithm for ECT J. Nucl. Med. 27 1577–85
[41] Hoffman E J, Cutler P D, Digby W M and Mazziotta J C 1990 3-D phantom to
simulate cerebral blood flow and metabolic images for PET IEEE Trans. Nucl.
Sci. 37 616–20
[42] Pollard K R, Bice A N, Eary J F, Durack L D and Lewellen T K 1992 A method for
imaging therapeutic doses of iodine-131 with a clinical gamma camera J. Nucl. Med.
33 771–6
[43] Dewaraja Y K, Ljungberg M and Koral K F 2000 Accuracy of 131I tumor
quantification in radioimmunotherapy using SPECT imaging with an ultra-high-
energy collimator: Monte Carlo study J. Nucl. Med. 41 1760–7
[44] Zaidi H 1996 Assessment of thyroid volume with pinhole emission computed
tomography Physica Medica XII 97–100
[45] Huili W, Jaszczak R J and Coleman R E 1996 Monte Carlo modeling of penetration
effect for iodine-131 pinhole imaging IEEE Trans. Nucl. Sci. 43 3272–7
[46] Smith M F and Jaszczak R J 1997 The effect of gamma ray penetration on angle-
dependent sensitivity for pinhole collimation in nuclear medicine Med. Phys. 24
1701–9
[47] Dewaraja Y K, Ljungberg M and Koral K F 2001 Monte Carlo evaluation of object
shape effects in iodine-131 SPECT tumor activity quantification Eur. J. Nucl. Med.
28 900–6
[48] Dewaraja Y K, Ljungberg M and Koral K F 2000 Characterization of scatter and
penetration using Monte Carlo simulation in 131I imaging J. Nucl. Med. 41 123–30
[49] Ellett W H and Humes R M 1971 Absorbed fractions for small volumes containing
photon-emitting radioactivity J. Nucl. Med. Suppl. 5 25–32
[50] Snyder W, Ford M R and Warner G 1978 Estimates of Specific Absorbed Fractions
for Photon Sources Uniformly Distributed in Various Organs of a Heterogeneous
Phantom, NM/MIRD (Society of Nuclear Medicine Publication, New York,
Pamphlet No 5)
[51] Akabani G, Poston J W and Bolch W E 1991 Estimates of beta absorbed fractions in
small tissue volumes for selected radionuclides J. Nucl. Med. 32 835–9
[52] Siegel J A and Stabin M 1994 Absorbed fractions for electrons and beta particles in
spheres of various sizes J. Nucl. Med. 35 152–6
[53] Cristy M 1983 Applying the reciprocal dose principle to heterogeneous phantoms:
practical experience from Monte Carlo studies Phys. Med. Biol. 28 1289–303
[54] Cristy M and Eckerman K F 1987 Specific Absorbed Fractions of Energy at Various
Ages from Internal Photon Sources I Methods, II One Year Old, III Five Year Old, IV
Ten Year Old, V Fifteen Year Old Male and Adult Female, VI New-Born and VII
Adult Male ORNL/TM 8381/V1-V7 Oak Ridge National Laboratory.
[55] Poston J W, Kodimer K A and Bolch W E 1996 A revised model for the calculation
of absorbed energy in the gastrointestinal tract Health Phys. 71 307–14
[56] Stabin M and Konijnenberg M W 2000 Re-evaluation of absorbed fractions for
photons and electrons in spheres of various sizes J. Nucl. Med. 41 149–60
[57] Bice A N, Links J M, Wong D F and Wagner H N 1985 Absorbed fractions for dose
calculations of neuroreceptor PET studies Eur. J. Nucl. Med. 11 127–31
[58] Bouchet L G, Bolch W E, Weber D A, Atkins H L and Poston J W 1999 MIRD
Pamphlet No 15: Radionuclide S values in a revised dosimetric model of the adult
head and brain. Medical Internal Radiation Dose J. Nucl. Med. 40 62S–101S
References 51
[59] Bolch W E et al. 1999 MIRD pamphlet No 17: the dosimetry of nonuniform activity
distributions—radionuclide S values at the voxel level. Medical Internal Radiation
Dose Committee J. Nucl. Med. 40 11S–36S
[60] Bouchet L G, Bolch W E, Howell R W and Rao D V 2000 S values for radionuclides
localized within the skeleton J. Nucl. Med. 41 189–212
[61] Bouchet L G and Bolch W E 1999 Five pediatric head and brain mathematical
models for use in internal dosimetry J. Nucl. Med. 40 1327–36
[62] Mardirossian G et al. 1999 A new rectal model for dosimetry applications J. Nucl.
Med. 40 1524–31
[63] Thomas S R, Stabin M, Chen C T and Samaratunga R C 1999 MIRD Pamphlet No
14 revised: A dynamic urinary bladder model for radiation dose calculations Task
Group of the MIRD Committee, Society of Nuclear Medicine J. Nucl. Med. 40
102S–123S
[64] Clairand I et al. 2000 Improvement of internal dose calculations using mathematical
models of different adult heights Phys. Med. Biol. 45 2771–85
[65] Peter J, Gilland D R, Jaszczak R J and Coleman R E 1999 Four-dimensional
superquadric-based cardiac phantom for Monte Carlo simulation of radiological
imaging systems IEEE Trans. Nucl. Sci. 46 2211–17
[66] Segars W P, Lalush D S and Tsui B M W 2001 Modeling respiratory mechanics in
the MCAT and spline-based MCAT phantoms IEEE Trans. Nucl. Sci. 48 89–97
[67] Xu X G, Chao T C and Bozkurt A 2000 VIP-Man an image-based whole-body adult
male model constructed from color photographs of the Visible Human Project for
multi-particle Monte Carlo calculations Health Phys. 78 476–86
[68] Chao T C and Xu X G 2001 Specific absorbed fractions from the image-based VIP-
Man body model and EGS4-VLSI Monte Carlo code internal electron emitters
Phys. Med. Biol. 46 901–27
[69] Mackie T R, Bielajew A F, Rogers D W and Battista J J 1988 Generation of
photon energy deposition kernels using the EGS Monte Carlo code Phys. Med.
Biol. 33 1–20
[70] Furhang E E, Sgouros G and Chui C S 1996 Radionuclide photon dose kernels for
internal emitter dosimetry Med. Phys. 23 759–64
[71] Berger M J 1971 MIRD Pamphlet 7: Distribution of absorbed doses around point
sources of electrons and beta particles in water and other media J. Nucl. Med. 12
5–23
[72] Simpkin D J and Mackie T R 1990 EGS4 Monte Carlo determination of the beta
dose kernel in water Med. Phys. 17 179–186
[73] Bardies M and Myers M J 1996 Computational methods in radionuclide dosimetry
Phys. Med. Biol. 41 1941–55
[74] Halbleib J A, Kensek R P, Valdez G D, Seltzer S M and Berger M J 1992 ITS: the
Integrated TIGER Series of electron/photon transport codes—version 3.0 IEEE
Trans. Nucl. Sci. 39 1025–30
[75] Sgouros G et al. 1993 Three-dimensional dosimetry for radioimmunotherapy
treatment planning J. Nucl. Med. 34 1595–601
[76] Johnson J C, Langhorst S M, Loyalka S K, Volkert W A and Ketring A R 1992
Calculation of radiation dose at a bone-to-marrow interface using Monte Carlo
modeling techniques (EGS4) J. Nucl. Med. 33 623–8
[77] Kawrakow I 2000 Accurate condensed history Monte Carlo simulation of electron
transport. I. EGSnrc, the new EGS4 version Med. Phys. 27 485–98
52 Monte Carlo techniques in nuclear medicine dosimetry
[78] Simpkin D J, Cullom S J and Mackie T R 1992 The spatial and energy dependence of
bremsstrahlung production about beta point sources in H2O Med. Phys. 19 105–14
[79] Leichner P K 1994 A unified approach to photon and beta particle dosimetry J. Nucl.
Med. 35 1721–9
[80] Casciari J J, Graham M M and Rasey J S 1995 A modeling approach for quantifying
tumor hypoxia with [F-18]fluoromisonidazole PET time–activity data Med. Phys. 22
1127–39
[81] Millet P, Delforge J, Pappata S, Syrota A and Cinotti L 1996 Error analysis on
parameter estimates in the ligand-receptor model application to parameter
imaging using PET data Phys. Med. Biol. 41 2739–56
[82] Burger C and Buck A 1996 Tracer kinetic modelling of receptor data with
mathematical metabolite correction Eur. J. Nucl. Med. 23 539–45
[83] Thie J A, Smith G T and Hubner K F 1997 Linear least squares compartmental-
model-independent parameter identification in PET IEEE Trans. Med. Imaging 16
11–16
[84] Welch A, Smith A M and Gullberg G T 1995 An investigation of the effect of finite
system resolution and photon noise on the bias and precision of dynamic cardiac
SPECT parameters Med. Phys. 22 1829–36
[85] Zeng G L, Gullberg G T and Huesman R H 1995 Using linear time-invariant system
theory to estimate kinetic parameters directly from projection measurements IEEE
Trans. Nucl. Sci. 42 2339–46
[86] Ping-Chun C, Rogers W L, Clinthorne N H, Fessler J A and Hero A O 1994 Model-
based estimation for dynamic cardiac studies using ECT IEEE Trans. Med. Imaging
13 217–26
[87] Huesman R H, Reutter B W, Zeng G L and Gullberg G T 1998 Kinetic parameter
estimation from SPECT cone-beam projection measurements Phys. Med. Biol. 43
973–82
[88] Simpkin D J 1989 Shielding requirements for constant-potential diagnostic x-ray
beams determined by a Monte Carlo calculation Health Phys. 56 151–64
[89] Metzger R, Richardson R and Van Riper K A 1993 A Monte Carlo model for
retrospective analysis of shield design in a diagnostic x-ray room Health Phys. 65
164–71
[90] Biggs P J 1991 Calculation of shielding door thicknesses for radiation therapy
facilities using the ITS Monte Carlo program Health Phys. 61 465–72
[91] Verhaegen F, Buffa F M and Deehan C 2001 Quantifying effects of lead shielding in
electron beams: a Monte Carlo study Phys. Med. Biol. 46 757–69
[92] da Silva A X and Crispim V R 2001 Moderator–collimator–shielding design for
neutron radiography systems using 252Cf Appl. Radiat. Isot. 54 217–25
[93] Howard W B and Yanch J C 1995 Shielding design and dose assessment for
accelerator based neutron capture therapy Health Phys. 68 723–30
[94] Kramer G H, Chamberlain M J and Yiu S 1997 A study of thyroid radioiodine
monitoring by Monte Carlo simulations: implications for equipment design Phys.
Med. Biol. 42 2175–82
[95] Kramer G H and Crowley P 2000 The assessment of the effect of thyroid size and
shape on the activity estimate using Monte Carlo simulation Health Phys. 78 727–38
[96] Rodriguez Villafuerte M, Gamboa deBuen I and Brandan M E 1997 Monte Carlo
simulation of depth-dose distributions in TLD-100 under 90Sr–90Y irradiation
Health Phys. 72 574–8
References 53
[116] Amaldi U 1999 Cancer therapy with particle accelerators Nucl. Phys. A 654 375–99
[117] Chu W T 1999 Hadron therapy Biomedical Uses of Radiation—Therapeutic
Applications ed W R Hendee (Weinheim: VCH–Wiley) Part B, pp 1057–131
[118] Sisterson J 2001 Particles PTCOG Newsletter 27 January
[119] Kanai T et al. 1999 Biophysical characteristics of HIMAC clinical irradiation system
for heavy-ion radiation therapy Int. J. Radiat. Oncol. Biol. Phys. 44 201–10
[120] Rosso D 2000 A boost for hadron therapy Phys. World
[121] Berger M J 1993 Proton Monte Carlo Transport Program PTRAN Report NISTIR
5113 (Bethesda, MD: NIST)
[122] Del Guerra A, Di Domenico G and Mukhopadhayay D 1997 PET dosimetry in
proton radiotherapy: a Monte Carlo Study Appl. Radiat. Isot. 18 1617–24
[123] Parodi K and Enghardt W 2000 Potential application of PET in quality assurance of
proton therapy Phys. Med. Biol. 45 N151–6
Chapter 3
3.1. INTRODUCTION
55
56 Medical imaging techniques for radiation dosimetry
chapter. Also, monitoring radiation dose during therapy may be the purpose
of interest, but you may still choose to do it from evaluative imaging rather
than from intra-therapy imaging. This choice has the advantage of lower
radiation exposure to imaging personnel, and of having no worries about
camera deadtime.
similar to the use of case-specific rules referred to above. That is, Erdi et al.
propose an initial segmentation, after which an approximate target to back-
ground ratio is calculated. In addition, the target volume is estimated from
CT. Then, depending on the ratio and the volume, a threshold is chosen
with which to carry out the final segmentation. Another approach, which
is receiving considerable attention, is the use of fuzzy segmentation techni-
ques. They demonstrate excellent performance and produce good results as
automated, unsupervised tools for segmenting noisy images in a robust
manner. Based on the Fuzzy C-Means (FCM) algorithm, Boudraa et al.
segmented PET images into a given number of clusters to extract specific
areas of differing activity concentration (MBq cm3 ) for brain lesion quanti-
fication purposes [11].
With SPECT or PET imaging, if a registration between the radio-
pharmaceutical tomographic image set and CT or MRI is available, the
VoI outlines from CT or MRI can simply be transferred to the three-
dimensional space of the reconstructed activity distribution [12]. If a
multimodality imager is employed, a phantom calibration provides the
transformation to place the CT outline into the activity space [13]. In that
space, the outlines provide the VoI over which the total activity is found,
or, if pixel-by-pixel dose, or pixel-by-pixel dose rate, is calculated, the
volume over which to investigate the dose parameter. Image registration is
discussed further in section 3.4.7.
therapy of lymphoma, the correlation is best for treated patients who had
not previously undergone bone-marrow transplantation when the nadir in
platelet counts was used [16].
Advantages. This imaging is simple and self-calibrating.
Disadvantages. There is no resolution in space.
also requires that a CT image set, registered with the projection image, is
available.
Advantages. These related methods are simple and lead directly to
macrodose. The simplicity is an aid for making the repeated measurements
needed for time–activity curves. Newer methods can employ information
from an anatomical image.
Disadvantages. Outlining each target region in one of the projections is
subject to error. Resolution in the anterior-to-posterior direction is either
non-existent or not of the highest quality. In some cases, results are obtained
for only one large target when there are actually two or more smaller targets.
Figure 3.1. Plot of activity as a function of the radius ri of the spherical volume
of interest. OSEM is defined in the text. Two other methods of reconstruction
are also presented for comparison. FBP ¼ standard filtered backprojection.
MAP/ICD ¼ maximum a posteriori/iterative coordinate descent. Reprinted
from [29] # 2001 IEEE, with permission.
that adapt the quantification to the particular case. These other procedures
include the use of regression equations [8]. In cases where there is no uniform
surrounding activity and no nearby objects that contain activity, a volume of
interest that exactly fits the object is used but only after expansion in all direc-
tions by one pixel to account for finite SPECT resolution [33]. The regression
equations relate the true activity to the measured activity with a constant
greater than 1. The constant is smaller when VoI expansion has been
employed, as is to be expected [33]. Also for 131 I, Fleming and colleagues
measured the activity values for variously shaped objects in a large cylindri-
cal water bath experimentally. Activity values ranged from 0.8 to 29.3 MBq.
Using segmentation based on maximization of interclass variance (see
section 3.3 above), the measured values compared with the true activities
fall quite close to the line of identity [8]. With some of the procedures of
Fleming and colleagues, it appears that the segmentation of the object of
interest and the determination of its activity are intimately connected,
rather than being two distinct steps.
Tang et al. have developed a set of methods for SPECT quantification
using a research multimodality CT-SPECT imager. They employ attenuation
correction and model the collimator but do not include scatter–penetration
correction. For three neuroblastoma patients being treated with 131 I meta-
iodobenzylguanidine, each with a large and a small tumour, they list activity
concentration results from eight methods, including conjugate views [13].
Distribution of radioactivity in both space and time 65
Table 3.1. Ratio indicating the variation in activity concentration with the
evaluation method as a function of tumour size and patient number.
Herein, we exclude results from one method in a baseline mode that did not
include compensation for collimator effects. Then table 3.1 shows the ratio of
the largest estimate of activity concentration to the smallest by patient and
tumour size.
For the large tumours, the ratio ranges from 1.9 to 4.2. For the harder-
to-estimate small tumours, it ranges from 13 to 21. The authors explain that
one of these methods is likely to yield a lower bound on the activity while at
least one other is likely to yield an upper bound. Excluding the largest and
smallest result from the seven methods, one still gets a large variation as
shown in the fourth column of the table: the ratio is never less than a
factor of 1.3 for the large tumours and not less than a factor of 2.0 for the
small tumours. Clearly, for a given patient tumour, results can vary by a
considerable factor depending on the method of quantification employed.
So, the choice of method is important and, perhaps, difficult to optimize
for all targets.
Koral et al. [14] employ a conversion factor for focal activity, C, that
varies both with the radius of rotation, r, of the circular orbit of the
acquisition (which affects resolution) and also with a measured parameter,
b, that indicates the relative magnitude of count spill-in from a uniform
background concentration compared with count spill-out from the target.
Phantom measurements of the functional dependence of C on r and b are
made using a 200 cm3 sphere located inside an elliptical water-bath phantom.
In contrast to some of the methods of Tang [34], no explicit subtraction for
the effects of spill-in of counts from the background is needed. Presumably,
the noise-amplification of a subtraction is thus not present. The authors
emphasize that the parameters for the SAGE reconstruction [35] of their
calibration phantom are the same as for their target measurements, and
that the identity of the method for phantoms and patients is also true for
the attenuation compensation, the outlining of the target of interest, and
the scatter–penetration compensation. The attenuation compensation is
66 Medical imaging techniques for radiation dosimetry
Figure 3.3. These images show a slice reconstructed by iterative MLEM after
1, 5, 10 and 50 iterations. The images show a 131 I-labelled monoclonal antibody
in a patient at a slice location just above the bladder. Note the increase in noise
as the number of iterations becomes large.
As seen in figure 3.3, the number of iterations has an effect on the result
and so must be chosen carefully.
Newer algorithms use only sub-sets in a certain order before the image is
updated. This procedure has been shown to decrease the number of needed
iterations dramatically and is therefore one of a set of methods called accel-
erated methods. OSEM is such a commonly used method [32]. The iterative
methods do not have the streak artefact and are generally expected to be
more quantitatively accurate, especially in typical oncology ‘low-statistics’
studies.
Figure 3.4. Plot of linear attenuation coefficient for gamma rays from 131 I
versus the pixel CT number from an X-ray CT scanner with an effective
energy of 76 keV. It is assumed that above CT ¼ 485 HU (Hounsfield units)
the material is bonelike and below CT ¼ 25 HU tissue-like. Between
25 HU and 485 HU, it is assumed there is a transition region.
70 Medical imaging techniques for radiation dosimetry
are used to model a spatially variant scatter distribution. This so called ESSE
method utilizes an attenuation map to compensate for non-homogeneous
attenuation. A fast Monte Carlo-based scatter compensation method has
also been developed for SPECT [47–49]. This method uses a rotation-
based Monte Carlo simulation that calculates scatter using various variance
reduction methods. An application is to calculate down-scatter in dual-
isotope SPECT. With this method, only a fraction of all SPECT projections
have to be simulated in a standard manner. The other projections can be
estimated rapidly using the results of these standard calculations. For
efficiency, approximations have to be made in the method with regard to
the final scatter angle of the detected photons. Further speed-up is obtained
by combining it with convolution-based forced detection instead of standard
forced detection, which is the more common variance reduction technique.
Rotation-based Monte Carlo with convolution-based forced detection was
found to be about three orders of magnitude faster than standard Monte
Carlo with forced detection. Even though SPECT has great potential, a
major drawback is its limited spatial resolution [50], mainly caused by the
inherent design of collimators. It is rare that the reconstructed spatial
resolution using commercial collimators is less than 10 mm for commonly
used radioisotopes in radionuclide therapy. Several groups have developed
methods for compensation for distance-dependent blurring. These methods
are basically attacking the previously mentioned activity spill-in and spill-
out. The methods can be based on Fourier-based methods, such as the
energy–distance principle [51] and Wiener filtering [52], but collimator-blur
compensation is mostly incorporated in the iterative reconstruction [47,
49]. The advantage with incorporation in the iterative algorithm is that
pre-calculated response functions can be placed into lookup tables and
selected for use based on the corresponding pixel location.
In many cases in nuclear medicine, a low count-rate produces images
with disturbing noise. The common way to reduce noise is to apply two-
or three-dimensional low-pass filters. However, it is important to recognize
that low-pass filtering of projections for noise reduction also will result in
a degradation of image resolution. For small objects, blurring due to low-
pass filtering reduces image contrast and the ability to detect lesions. The
blurring also affects, for example, segmentation of volumes of interest,
which may incorrectly increase the determined volume. Furthermore, the
limited spatial resolution will produce ‘spill-in’ and ‘spill-out’ of activity
from one region to another.
Figure 3.5. Results for one slice from a three-dimensional registration between
a CT image set and a SPECT image set. The pelvic area of a non-Hodgkin’s
lymphoma patient (# 33) being treated with 131 I tositumomab is shown. The
CT slice (at the top) is from a scan that was taken 13 days before the treatment
infusion. The two tumours, shown by black outlines, had their edges manually
chosen in this image. The SPECT slice (at the bottom) is from a scan that
was taken 42 hours after the treatment infusion. The tumour outlines were
placed on the bottom image by virtue of the registration. The anterior pelvic
tumour (ant pel. tum.) could be evaluated even though its percent injected
dose per gram was only 0.68 times that of the left pelvic tumour (left pel.
tum.). 131 I activity in the urine of the bladder (not labelled) is a strong feature
of the SPECT image. (See plate 2 for colour version.)
of PET per disintegration with comparable axial fields of view is two orders of
magnitude greater than that of SPECT cameras. Another boost for sensitivity
comes from the fact that PET usually detects molecules labelled with short-
lived radioisotopes for which the decay rate is very high. Within the spectrum
of medical imaging modalities, sensitivity ranges from the detection of milli-
molar concentrations in functional MRI to pico-molar concentrations in
PET which is a difference involving a factor of 109 [74].
It is also generally held that PET offers better spatial resolution than
SPECT. Furthermore, PET has provided the possibility of quantitative
measurements of tracer concentration in vivo for a longer time than
SPECT has. However, in practice, there are several issues that must be
addressed in order to fully realize the potential of PET. That is, the
measured line integrals must be corrected for a number of background
and physical effects before reconstruction. The needed processes include
subtraction of random coincidences, detector normalization, deadtime,
attenuation and scatter corrections. Detector normalization is similar to
correction for camera–collimator-system non-uniformities based on field-
flood measurements in SPECT. Deadtime, attenuation and scatter correc-
tion are all also needed in SPECT. Since attenuation and scatter correction
are especially important in clinical PET, they are briefly discussed in section
3.4.9 below.
During the past decade, much research and development has concen-
trated on the development of improved fully three-dimensional reconstruction
algorithms and scatter compensation techniques required for quantitative
PET. Increasingly sophisticated scatter correction procedures are under
investigation: particularly those based on accurate scatter models and
iterative reconstruction-based scatter compensation approaches [76]. In
summary, PET quantification remains an open issue that requires further
R&D efforts [77].
When using imaging for dosimetry, one often assumes that the results during
therapy will be the same as those during evaluation, except for scaling by the
ratio of the administered activities. Since the administration for the first
imaging could possibly affect the pharmacokinetics, this assumption may
not hold. If it does not, the choice of administered activity might still be
made in the same way but measurements would be needed during therapy
to quantify the discrepancy. It is well known that in ablation of thyroid
remnants, measurements seem to imply a reduction in the therapy activity
compared with what one calculates by scaling. The explanation for the
observed results is still under discussion, however.
In BexxarTM RIT of lymphoma, the ratio of the therapy percent–
injected dose divided by the evaluation percent–injected dose has been
checked for variations from 1. These checks were carried out at 2 to 3 days
after the administration [81]. There is an extra biological complication
here in that the RIT consists of the administration of an anti-CD20
murine monoclonal IgG2a antibody (tositumomab) pre-dose followed by
References 79
its 131 I-radiolabelled conjugate (iodine 131 I tositumomab). The evaluation was
carried out with a reduced activity of iodine 131 I tositumomab but with the
same pre-dose. For 31 chemotherapy-relapsed patients, the hypothesis that
the ratio equals 1 was rejected using conjugate-view imaging [82]. However,
the average of the ratio was only slightly less than one (0.93). For seven
previously untreated patients, results from imaging using either conjugate
views or SPECT agreed with the hypothesis of the ratio being equal to 1
[81]. So, for both types of patient, there will not be much error on average
with assuming a ratio of 1.
ACKNOWLEDGMENTS
REFERENCES
[1] Leichner P K et al. 1993 An overview of imaging techniques and physical aspects of
treatment planning in radioimmunotherapy Med. Phys. 20 569–77
[2] Siegel J A et al. 1999 MIRD Pamphlet No 16: Techniques for quantitative
radiopharmaceutical biodistribution data acquisition and analysis for use in human
radiation dose estimates J. Nucl. Med. 40 37S–61S
[3] Snyder W S, Ford M R, Warner G G and Watson S B 1975 ‘S’, Absorbed Dose per
Unit Cumulated Activity for Selected Radionuclides and Organs. MIRD Pamphlet
No 11 (Society of Nuclear Medicine)
[4] Zaidi H and Hasegawa B H 2002 Determination of the attenuation map in emission
tomography J. Nucl. Med. 43 in press
[5] Liu A, Williams L E and Raubitschek A A 1996 A CT assisted method for absolute
quantitation of internal radioactivity Med. Phys. 23 1919–28
[6] Zaidi H 1996 Comparative methods for quantifying thyroid volume using planar
imaging and SPECT J. Nucl. Med. 37 1421–6
[7] Alaamer A S, Fleming J S and Perring S 1993 Evaluation of a technique for the
quantification of radioactivity and volume of an object using SPECT Nucl. Med.
Commun. 14 1061–70
[8] Fleming J S and Alaamer A S 1998 A rule based method for context sensitive
threshold segmentation in SPECT using simulation Phys. Med. Biol. 43 2309–23
[9] Zaidi H 1996 Organ volume estimation using SPECT IEEE Trans. Nucl. Sci. 43
2174–82
80 Medical imaging techniques for radiation dosimetry
[10] Erdi Y E et al. 1997 Segmentation of lung lesion volume by adaptive positron
emission tomography image thresholding Cancer 80 2505–9
[11] Boudraa A E et al. 1996 Delineation and quantitation of brain lesions by fuzzy
clustering in positron emission tomography Comput. Med. Imaging Graph. 20
3141
[12] Koral K F et al. 1989 Calculating radiation absorbed dose for pheochromocytoma
tumors in 131I MIBG therapy Int. J. Radiat. Oncol. Biol. Phys. 17 21118
[13] Tang H R et al. 2001 Neuroblastoma imaging using a combined CT scanner–
scintillation camera and 131I-MIBG J. Nucl. Med. 42 23747
[14] Koral K F et al. 2000 Initial results for hybrid SPECTconjugate-view tumor
dosimetry in 131I-anti-B1 antibody therapy of previously untreated patients with
lymphoma J. Nucl. Med. 41 157986
[15] Sisson J C et al. 1994 Predictors of toxicity in treating patients with neuroblastoma by
radiolabeled metaiodobenzylguanidine Eur. J. Nucl. Med. 21 4652
[16] Wahl R L et al. 1998 Iodine-131 anti-B1 antibody for B-cell lymphoma an update on
the Michigan Phase I experience J. Nucl. Med. 39 21S27S
[17] Thomas S R, Maxon H R and Kereiakes J G 1976 In vivo quantitation of lesion
radioactivity using external counting methods Med. Phys. 3 2535
[18] Hammond N D, Moldofsky P J, Beardsley M R and Mulhern C B Jr 1984 External
imaging techniques for quantitation of distribution of I-131 Fðab0 Þ2 fragments of
monoclonal antibody in humans Med. Phys. 11 77883
[19] Eary J F, Appelbaum F L, Durack L and Brown P 1989 Preliminary validation of
the opposing view method for quantitative gamma camera imaging Med. Phys. 16
3827
[20] Ogawa K, Harata H, Ichihara T, Kubo A and Hashimoto S 1991 A practical method
for position dependent Compton-scatter correction in single photon emission CT
IEEE Trans. Med. Imaging 10 40812
[21] Shulkin B L et al. 1988 Conjugate view gamma camera method for estimating tumor
uptake of iodine-131 metaiodobenzylguanidine J. Nucl. Med. 29 5428
[22] Matthay K K et al. 2001 Correlation of tumor and whole-body dosimetry with tumor
response and toxicity in refractory neuroblastoma treated with 131I-MIBG J. Nucl.
Med. 42 171321
[23] Doherty P, Schwinger R, King M and Gionet M 1985 Distribution and dosimetry of
indium-111 labeled Fðab0 Þ2 fragments in humans. Fourth International Radiopharma-
ceutical Dosimetry Symposium, Oak Ridge, TN, CONF-851113 (DE86010102)
[24] Zaidi H 1996 Quantitative SPECT: Recent developments in detector response
attenuation and scatter correction techniques Physica Medica 12 10117
[25] Rosenthal M S et al. 1995 Quantitative SPECT imaging a review and
recommendations by the Focus Committee of the Society of Nuclear Medicine
Computer and Instrumentation Council J. Nucl. Med. 36 1489513
[26] Tsui B M, Zhao X, Frey E C and McCartney W H 1994 Quantitative single-photon
emission computed tomography basics and clinical considerations Semin. Nucl. Med.
24 3865
[27] Dewaraja Y K, Ljungberg M and Koral K F 2000 Accuracy of 131I tumor
quantification in radioimmunotherapy using SPECT imaging with an ultra-high-
energy collimator: Monte Carlo study J. Nucl. Med. 41 176067
[28] Sayeram S et al. 2001 Compensation for collimator detector response in tumor
SPECT imaging using In-111 with dual energy photons J. Nucl. Med. 42 197P
References 81
[47] Beekman F J, de Jong H W and Slijpen E T 1999 Efficient SPECT scatter calculation
in non-uniform media using correlated Monte Carlo simulation Phys. Med. Biol. 44
N183N192
[48] de Jong H W and Beekman F J 2001 Rapid SPECT simulation of downscatter in non-
uniform media Phys. Med. Biol. 46 62135
[49] Beekman F J, Slijpen E T, de Jong H W and Viergever M A 1999 Estimation of the
depth-dependent component of the point spread function of SPECT Med. Phys. 26
231122
[50] King M A, Long D T and Brill B A 1991 SPECT volume quantitation: Influence of
spatial resolution source size and shape and voxel size Med. Phys. 18 101624
[51] Hawkins W G, Yang N-C and Leichner P K 1991 Validation of the circular harmonic
transform (CHT) algorithm for quantitative SPECT J. Nucl. Med. 32 14150
[52] King M A, Schwinger R B, Doherty P W and Penney B C 1984 Two-dimensional
filtering of SPECT images using the Metz and Wiener filters J. Nucl. Med. 25 123440
[53] Andersson J L R, Sundin A and Valind S 1995 A method for coregistration of PET
and MR brain images J. Nucl. Med. 36 130715
[54] Borgefors G 1988 Hierachical chamfer matching a parametrical edge matching
algorithm IEEE Trans. Pattern Anal. Mach. Intell. 10 84965
[55] Eberl S et al. 1996 Automated interstudy image registration technique for SPECT and
PET J. Nucl. Med. 37 13745
[56] Flux G 1995 Multimodality Image Registration and its Application to the Dosimetry of
Intralesional Radionuclide Therapy. Thesis. Institute of Cancer Research & Royal
Marsden Hospital, University of London
[57] Hill D L G and Hawkes D J 2000 Across-modality registration using intensity-based
cost functions. In Handbook of Medical Imaging—Processing and Analysis ed I N
Bankman (San Diego: Academic)
[58] Maes F, Collignon A, Vandermeulen D, Marchal G and Suetens P 1997
Multimodality image registration by maximization of mutual information IEEE
Trans. Med. Imaging 16 18798
[59] Maguire G Q, Noz M E, Lee E M and Schimpf J H 1986 Correlation methods for
tomographic images using two and three dimensional techniques. In Information
Processing in Medical Imaging ed S L Bacharach (Dordrecht: Martinus Nijhoff )
pp 26679
[60] Maintz J B A, van den Elsen P A and Viergever M A 2001 3D multimodality medical
image registration using morphological tools Imag. Vision Comput. 19 5362
[61] Pelizzari C A, Chen G T Y, Spelbring D R, Wiechselbaum R R and Chen C 1989
Accurate three-dimensional registration of CT, PET, and/or MR images of the
brain J. Comput. Assist. Tomogr. 13 206
[62] Rangarajan A, Chui H and Duncan J S 1999 Rigid point feature registration using
mutual information Med. Imag. Anal. 3 42540
[63] Studholme C, Hill D L G and Hawkes D J 1996 Automated 3D registration of MR
and CT images of the head Med. Imag. Anal. 1 16375
[64] Studholme C, Hill D L G and Hawkes D J 1999 An overlap invariant entropy
measure of 3D medical image alignment Pattern Recogn. 32 7186
[65] van den Elsen P A, Pol E J D and Viergever M A 1993 Medical image matching—A
review with classification IEEE Eng. Med. Biol. Mag. 12 2639
[66] Wells W M, Viola P, Atsumi H, Nakajima S and Kikinis R 1996 Multi-modal volume
registration by maximization of mutual information Med. Imag. Anal. 1 3551
References 83
[67] Woods R P, Mazziotta J C and Cherry S R 1992 Rapid automated algorithm for
aligning and reslicing PET images J. Comput. Assist. Tomogr. 16 620–33
[68] Woods R P 2000 Within-modality registration using intensity-based cost functions. In
Handbook of Medical Imaging—Processing and Analysis ed I N Bankman (San Diego:
Academic)
[69] Hutton B F, Braun M, Thurfjell L and Lau D Y H 2002 Image registration an
essential tool for nuclear medicine Eur. J. Nucl. Med. 29 559–77
[70] Meyer C R et al. 1997 Demonstration of accuracy and clinical versatility of mutual
information for automatic multimodality image fusion using affine and thin-plate
spline warped geometric deformations Med. Imag. Anal. 1 195–206
[71] Sjögreen K, Ljungberg M, Wingšrdh K, Erlandsson K and Strand S-E 2001
Registration of emission and transmission whole-body scintillation camera images
J. Nucl. Med. 42 1563–1570
[72] Sjögreen K 2001 Image Processing for Quantitative Scintillation-Camera Imaging:
Application to Radionuclide Therapy. Thesis. Lund University
[73] DeVries D J and Moore S C 2001 Registration of CT and SPECT volumes using CT
reprojections and planar transmission images Proc IEEE Nuclear Science Symposium
and Medical Imaging Conference, San Diego, CA, in press
[74] Jones T 1996 The role of positron emission tomography within the spectrum of
medical imaging Eur. J. Nucl. Med. 23 207–11
[75] Zaidi H, Diaz-Gomez M, Boudraa A E and Slosman D O 2002 Fuzzy clustering-
based segmented attenuation correction in whole-body PET imaging Phys. Med.
Biol. 47 1143–60
[76] Zaidi H 2000 Comparative evaluation of scatter correction techniques in 3D positron
emission tomography Eur. J. Nucl. Med. 27 1813–26
[77] Zaidi H 2001 Scatter modelling and correction strategies in fully 3-D PET Nucl. Med.
Commun. 22 1181–4
[78] Beyer T et al. 2000 A combined PET/CT scanner for clinical oncology J. Nucl. Med.
41 1369–79
[79] Dimitrakopoulou-Strauss A, Strauss L G and Burger C 2001 Quantitative PET studies
in pretreated melanoma patients a comparison of 6-[18F]fluoro-L-dopa with 18F-FDG
and 15O-water using compartment and noncompartment analysis J. Nucl. Med. 42 248–
56
[80] Torizuka T et al. 2000 Metabolic response of non-Hodgkin’s lymphoma to 131I-anti-
B1 radioimmunotherapy evaluation with FDG PET J. Nucl. Med. 41 999–1005
[81] Koral K F et al. 1999 I-131 anti-B1 therapy/tracer uptake ratio using a new procedure
for fusion of tracer images to computed tomography images Clin. Cancer Res. 5
3004s–3009s
[82] Zasadny K R et al. 1997 Do tracer dosimetry studies predict normal organ and
tumor uptake of I-131-anti-B1 (anti CD-20) antibody in patients receiving
radioimmunotherapy for non-Hodgkins lymphoma? Radiology 205 509–509
Chapter 4
4.1. INTRODUCTION
Radiation dosimetry deals with the determination of the amount and the
spatial and temporal distribution of energy deposited in matter by ionizing
radiation. Internal radionuclide dosimetry specifically deals with radiation
energy deposition in tissue for radionuclides within the body. Internal
dosimetry has been applied to the determination of tissue doses and related
quantities for occupational exposures in radiation protection, environmental
exposures in radiation epidemiology, and diagnostic and therapeutic
exposures in nuclear medicine. With the increasing therapeutic application
of internal radionuclides in medicine, larger administered activities are
used and higher normal-tissue doses, with associated radiation injury, may
result. Treatment planning for radionuclide therapy thus requires, in
addition to reliable dose–response relationships for target tissues and
dose–toxicity relationships for normal tissues, increasingly more accurate
and precise dose estimates for target tissue and at-risk normal tissues.
Radiation dosimetry in nuclear medicine continues to evolve—from
estimation of population- and organ-averaged doses for risk estimation for
stochastic effects to calculation of individualized patient doses that reflect
the heterogeneity of dose distributions within tumours and organs and the
probability of deterministic effects [1–4]. Reviewed in this chapter, the
standard computational methods developed for internal radiation dosimetry,
such as the ‘MIRD schema’, have generally been based on ‘average’ kinetic
and anatomic (i.e., anthropomorphic) models. These methods yield popula-
tion- and organ-averaged doses for estimation of stochastic risk associated
with diagnostic administered activities. In other chapters of this book,
strategies for individualized patient dosimetry, more appropriate for radio-
nuclide therapy, are developed.
84
Radiation quantities and units 85
dQ
X ð4:1Þ
dm
where dQ is the absolute value of the total charge of ions of one sign
produced in air when all the electrons (negatrons and positrons) liberated
by photons (X- and -rays) in air are completely stopped in air, and dm is
the mass of air in which the electrons (negatrons and positrons) were liber-
ated by photons (X- and -rays).
Exposure is essentially the ionization equivalent of air kerma (defined
below) and, as in the case of kerma, it may often be convenient to refer to
a value of exposure at a point inside a material different from air. In such
a case, the exposure specified will be that which would be determined for a
small quantity of air placed at the point of interest. Exposure is a deter-
ministic quantity and thus corresponds to a sufficiently large irradiated
volume and/or a sufficiently large amount of radiation to yield insignificant
statistical fluctuations in its measured value.
The ‘absorbed dose’, D, is probably the most widely used quantity to
characterize ‘dose’ and is defined as
dE
D ð4:2Þ
dm
where dE is the mean energy imparted by ionizing radiation to matter, and
dm is the mass of matter to which the energy is imparted.
86 Computational methods in internal radiation dosimetry
the absorbed dose (in rad) may be considered to be numerically equal to the
exposure (in R).
‘Kerma’, K, is defined as
dEtr
K ð4:7Þ
dm
where dEtr is the sum of the initial kinetic energies of all the charged ionizing
particles liberated by uncharged ionizing particles (including photons) in
matter, and dm is the mass of matter in which the charged ionizing particles
were liberated.
In the case in which the matter is air, kerma is often referred to as ‘air
kerma’ or ‘free air kerma’.
The quality as well as the quantity of radiation are important deter-
minants of the frequency and/or severity of radiogenic effects. The ‘quality’
of a radiation is related to the characteristics of the microscopic spatial
distribution of energy-deposition events, determined by the mass, charge
and energy (i.e., velocity) of the charged particles composing the radiation
or, in the case of X-rays, -rays and neutrons, the charged particles produced
by the radiation. Sparsely ionizing radiations such as X- and -rays and
intermediate-to-high-energy electrons and particles are characterized as
‘low-quality’ radiations, while densely ionizing radiations such as low-
energy electrons (e.g., Auger electrons), protons, neutrons and particles
are typically characterized as ‘high-quality’ radiations. Importantly, for the
same absorbed dose, the frequency and/or severity of radiogenic biological
effects are generally less for sparsely ionizing, low-quality radiations than
for densely ionizing, high-quality radiations.
The quality of radiation is quantitatively characterized by the ‘linear
energy transfer’, L or LET, or the ‘restricted linear energy transfer’, L or
LET of a material for charged particles:
dE
L ð4:8Þ
dl
where dE is the energy lost by a charged particle in traversing a distance in
matter due to those collisions in which the energy loss is less than the
energy cut-off , and dl is the distance traversed in matter.
The restricted linear energy transfer (L ) requires the specification of a
cut-off energy (), necessitated by primary energy-deposition events result-
ing in relatively high-energy, relatively long-range secondary electrons (i.e.,
-rays) whose energy deposition events may be considered as separate from
those along the track of the primary radiation. In radiation protection,
this feature is generally disregarded by specification of the ‘unrestricted
linear energy transfer’, L1 or LET1 (also known simply as the ‘linear
energy transfer’, L or LET), where the energy cut-off () is set equal to
infinity (1).
88 Computational methods in internal radiation dosimetry
Like the RBE, the quality factors, mean quality factors and radiation weight-
ing factors are dimensionless quantities.
For a tissue or organ, the equivalent dose is conceptually different from
the dose equivalent. The dose equivalent is based on the absorbed dose at a
point in tissue weighted by the LET-dependent distribution of quality factors
(Q(L)) at that point. The equivalent dose, in contrast, is based on the average
absorbed doses (DT;R ) in the tissue or organ weighted by the radiation
weighting factor (wR ) for the radiation(s) actually impinging on the tissue
or organ or, in the case of internal radionuclides, as it is actually emitted
by the source.
One of the most widely used and authoritative approaches for internal dose
calculations in medicine was developed by the Medical Internal Radiation
Dose (MIRD) Committee of the Society of Nuclear Medicine (SNM) and
generally is referred to as the ‘MIRD schema’ or ‘MIRD formalism’ [4,
13–16]. The International Commission on Radiological Protection (ICRP)
has developed a similar methodology and similar reference data [17]. The
MIRD schema, including notation, terminology, mathematical method-
ology, and reference data, has been disseminated in the form of the collected
MIRD Pamphlets and associated publications [2, 13, 18, 19]. With the
publication of ORNL/TM-8381/V1-7 [20], age- and gender-specific body
habiti other than the original 70 kg adult anthropomorphic model (‘Refer-
ence Man’ or ‘Standard Man’) [14, 21] are now incorporated into the
MIRD schema. In addition, several computerized versions of the MIRD
90 Computational methods in internal radiation dosimetry
Figure 4.1. ‘Step 1’ in the MIRD schema for internal radionuclide dosimetry:
compilation of the pertinent physical data (i.e., radionuclide-specific nuclear
decay data) and calculation of the pertinent physical parameters (i.e., absorbed
fractions and S factors for selected source region–target region pairs in a
mathematically formulated ‘average’ anatomic model). Adapted, in part,
from [14] and [20].
Figure 4.2. ‘Step 2’ in the MIRD schema for internal radionuclide dosimetry:
acquisition of the pertinent biological data (i.e., time–activity data for the
radiopharmaceutical under consideration in selected organs) and calculation
of the pertinent biological parameters (i.e., cumulated activities, or residence
times, in the selected organs).
region rk for radiation i emitted in source region rh , that is, the fraction of
energy of radiation i emitted in source region rh that is absorbed per unit
mass in target region rk
i ðrk rh Þ
i ðrk rh Þ ð4:16Þ
mk
and Sðrk rh Þ is the radionuclide-specific S factor for target region rk
and source region rh , that is, the absorbed dose to target region rk per unit
cumulated activity in source region rh
P
ðr rh Þ
Sðrk rh Þ i i i k ð4:17Þ
mk
An important simplification of absorbed dose calculations was intro-
duced by combining the radionuclide-specific equilibrium dose constant
i , the source region-to-target region absorbed fraction i ðrk rh Þ, and
the target region mass mk into a single quantity, the S factor, as defined by
equation (4.17) [15]. As a result, given the cumulated activity A~h in a given
source region rh , one can use the tabulated S factors to yield, by a simple
multiplication, the absorbed dose contribution D ðrk rh Þ to a target
region rk . Besides computationally simplifying the determination of
The MIRD schema 93
Figure 4.3. ‘Step 3’ in the MIRD schema for internal radionuclide dosimetry:
integration of the pertinent physical and biological parameters to calculate the
mean absorbed doses to selected target regions in the mathematically formu-
lated anatomic model. Adapted, in part, from [13].
absorbed dose, the S factor conceptually clarified this task by isolating all
non-time-dependent and, largely, non-biology-dependent dosimetric factors
into a single parameter.
The total mean absorbed dose Dðrk Þ to target region rk is then calculated
by summing of the absorbed dose contributions from all source regions rh :
X A~h Pi i i ðrk rh Þ
Dðrk Þ ¼ ð4:18aÞ
h
mk
X X
¼ A~h i i ðrk rh Þ ð4:18bÞ
h i
X
¼ A~h Sðrk rh Þ ð4:18cÞ
h
The time–activity function includes both the physical decay constant of the
radionuclide and the biological disappearance constant(s) (see the detailed
discussion below). It is usually determined by discrete serial measurements.
Despite the complexity of the multiple underlying biological processes,
the time–activity data of virtually any radiopharmaceutical in any tissue or
organ can be accurately represented by an exponential function:
X
Ah ðtÞ ¼ ðAh Þ eðh Þj t ð4:20Þ
j
where ðAh Þj is the (extrapolated) activity at time t ¼ 0 for the jth exponential
component of the time–activity function in source region rh , and ðh Þj is the
clearance constant of the jth exponential component of the time–activity
function in source region rh , that is, the fraction per unit time of activity
eliminated for the jth exponential component of the time–activity function
in source region rh .
The clearance constant, ðh Þj , is actually the ‘effective’ (e) clearance
constant because it includes the effects of both biological clearance and
physical (i.e., radioactive) decay:
ðh Þj ½ðe Þh j ð4:21aÞ
¼ p þ ½ðb Þh j ð4:21bÞ
where p is the physical decay constant of the radionuclide, and ½ðb Þh j is the
biological clearance constant of the jth exponential component of the time–
activity function in source region rh , that is, the fraction per unit time of
activity biologically eliminated for the jth exponential component of the
time–activity function in source region r.
For each exponential component, the half-life (or half-time) is related to
its corresponding clearance constant:
ln 2
T1=2 ¼ ð4:22aÞ
0:693
¼ ð4:22bÞ
and the effective half-time, Te ¼ ðT1=2 Þe , includes the effects of both the
biological and physical half-lives:
1 1
Te ¼ þ ð4:23aÞ
Tp Tb
Tp Tb
¼ ð4:23bÞ
Tp þ Tb
Dr Robert Loevinger, one of the originators of the MIRD schema, has stated
that ‘ . . . there is in principle no way of attaching a numerical uncertainty to
the profound mismatch between the patient and the model (the totality of all
assumptions that enter into the dose calculation). The extent to which the
model represents in some meaningful way a patient, or a class of patients, is
always open to question, and it is the responsibility of the clinician to make
that judgment’ [13]. In contrast to diagnostic radiopharmaceuticals, therapeutic
radiopharmaceuticals engender much smaller risk–benefit ratios and therefore
markedly reduced tolerances for such inaccuracies in dose estimation. To the
extent, then, that specific patients deviate kinetically and anatomically from
the respective kinetic and anatomic averages, tissue dose estimates will be
inaccurate and individualized (i.e., patient-specific) dosimetry may be required.
where
for water/soft tissue (see above)), np is the equilibrium dose constant for
non-penetrating radiations for 131 I (0.41 g rad/Ci h [27]).
if rh ¼ rk .
In principle, patient organ masses may be estimated by computed
tomography or magnetic resonance imaging. In practice, however, masses
of normal organs may not be available. A more practical, though certainly
less accurate, adaptation of equation (33) based on the patient and Reference
Man total body masses may then be used:
Patient- and position-specific dosimetry 99
Reference Man
patient Reference Man total body mass
ð4:34Þ
Sðrk rh Þ Sðrk rh Þ patient total
body mass
if rh ¼ rk and the patient target region ðrk Þ mass is not available, assuming,
arbitrarily, that organ mass is directly proportional to total body mass.
For a tumour-bearing organ, adaptation of the MIRD schema is some-
what more complicated. First, the self-irradiation absorbed dose and S factor
for the tumour-bearing organ can be separated into their penetrating and
non-penetrating radiation components:
Dðrh rh Þ Dp ðrh rh Þ þ Dnp ðrh rh Þ ð4:35Þ
Sðrh rh Þ Sp ðrh rh Þ þ Snp ðrh rh Þ ð4:36aÞ
P
i ðnp Þi ðnp Þi ðrh rh Þ
¼ Sp ðrh rh Þ þ ð4:36bÞ
mh
P
i ðnp Þi
¼ Sp ðrh rh Þ þ ð4:36cÞ
mh
np
¼ Sp ðrh rh Þ þ ð4:36dÞ
mh
where Sðrh rh Þ is the total self-irradiation S factor in source and target
region rh , Sp ðrh rh Þ is the self-irradiation S factor for penetrating radiations
in source and target region rh , Snp ðrh rh Þ is the self-irradiation S factor
for non-penetrating radiations in source and target region rh , ðnp Þi is the
equilibrium dose constant for non-penetrating radiation i, ðnp Þi ðrh rh Þ is
the self-irradiation absorbed fraction in non-penetrating radiations for
source and target region rh ,
ðnp Þi ðrh rh Þ ¼ 1 ð4:37Þ
P
mh is the mass of the source and target region rh , and np i ðnp Þi is the
total equilibrium dose constant for non-penetrating radiations.
Therefore,
np
Sp ðrh rh Þ ¼ Sðrh rh Þ : ð4:38Þ
mh
The patient-specific self-irradiation S factor for penetrating radiations can
then be calculated:
Reference Man
patient Reference Man target region ðrk Þ mass
: ð4:39Þ
Sðrh rh Þ Sðrh rh Þ patient target
region ðrk Þ mass
100 Computational methods in internal radiation dosimetry
use of Monte Carlo analysis has increased [39, 40, 80]. Monte Carlo-based
dosimetry can more accurately account for tissue variations in mass density
and atomic number as well as edge effects which may be important at the
periphery of the body and at soft tissue–lung and –bone interfaces [39].
For example, if the relevant distribution data were somehow available
(e.g., by autoradiography of biopsy specimens), Monte Carlo analysis
might be particularly applicable to normal lung dosimetry in radioiodine
treatment of metastatic thyroid cancer, particularly when in the form of
dosimetrically problematic miliary disease. This method remains computa-
tionally time-consuming, however [39]. Tabulations of voxel S factors,
conceptually equivalent to voxel source-kernels (the mean absorbed dose
to a target voxel per radioactive decay in a source voxel, both of which are
contained in an infinite homogeneous soft-tissue medium) [49, 59], are
becoming available [60]. In contrast to dose point-kernel- and Monte
Carlo-based techniques, the voxel S factor method does not require
specialized computer facilities and is relatively fast, and thus may emerge
as the practical method of choice for calculation of macroscopic non-
uniform dose distributions.
Once a dose distribution has been calculated, a corresponding ‘dose–
volume histogram’ can be derived [81]. Basically, a dose–volume histogram
is a graph of the fraction of the tumour or organ volume receiving a specified
dose versus the dose (differential form) or the fraction of the tumour or organ
volume receiving less than a specified dose versus the dose (integral, or
cumulative, form). It graphically presents the minimum, mean and maximum
doses and the dispersion about the mean dose. The greater this dispersion,
the more non-uniform is the dose distribution.
An important practical component of macroscopic non-uniform dosi-
metry is the ability to fuse, or register, tomographic images from multiple
modalities [46, 47, 56, 82–84]. Dose distributions, calculated from three-
dimensional activity distributions measured by scintigraphic imaging (i.e.,
SPECT or PET), may be presented as isodose contours or colour-coded
images. By image fusion, such isodose contours or colour-coded images
can be superimposed on or juxtaposed with the corresponding anatomy to
allow correlation of doses with tumour and normal organs (as imaged by
CT or MRI) [36, 56, 59, 60].
4.5. SUMMARY
Radiation dosimetry deals with the determination of the amount and the
spatial and temporal distribution of energy deposited in matter by ionizing
radiation. Internal radionuclide radiation dosimetry specifically deals with
the deposition of radiation energy in tissue due to a radionuclide within
the body. However, unlike external radiation dose (which can often be
102 Computational methods in internal radiation dosimetry
REFERENCES
[14] Snyder W, Ford M, Warner G and Fisher H Jr 1969 Estimates of absorbed fractions
for monoenergetic photon sources uniformly distributed in various organs of a
heterogeneous phantom. MIRD Pamphlet No 5 J. Nucl. Med. 10 (Suppl. 3) 5–52
[15] Snyder W, Ford M, Warner G and Watson S 1975 ‘S’, Absorbed dose per unit
cumulated activity for selected radionuclides and organs. Medical Internal
Radiation Dose (MIRD) Pamphlet No 11 (New York: Society of Nuclear Medicine)
[16] Zanzonico P, Brill A and Becker D 1995 Radiation dosimetry. In Principles of Nuclear
Medicine eds H Wagner, Z Szabo and J Buchanan (Philadelphia: Saunders) 2nd
edition pp 106–34
[17] ICRP 1988 Radiation Dose to Patients from Radiopharmaceuticals. International
Commission on Radiological Protection (ICRP) Publication 53 (New York:
Pergamon)
[18] Watson E 1999 Foreword J. Nucl. Med. 40 1S–2S
[19] Loevinger R 1989 The MIRD perspective. In Dosimetry of Administered
Radionuclides eds S Adelstein, A Kassis and R Burt (Washington, DC: American
College of Nuclear Physicians/Department of Energy) pp 29–43
[20] Cristy M and Eckerman K 1987 Specific Absorbed Fractions of Energy at Various
Ages from Internal Photon Sources (Parts I–VII). Oak Ridge National Laboratory
Report ORNL/TM-8381/V1-7 (Springfield, VA: National Technical Information
Service, Dept of Commerce)
[21] ICRP 1975 Report of the Task Group on Reference Man ICRP Publication 23 (Oxford:
Pergamon)
[22] Stabin M 1996 MIRDOSE: personal computer software for internal dose assessment
in nuclear medicine J. Nucl. Med. 37 538–46
[23] Sgouros G, Bigler R and Zanzonico P 1988 DOSCAL: a tumor-incorporating mean
absorbed dose calculation program (Abstract) J. Nucl. Med. 29 874
[24] Johnson T 1988 MABDOS: a generalized program for internal dosimetry Comput.
Meth. Programs Biomed. 37 159–67
[25] Johnson T K, McClure D and McCourt S 1999 MABDOSE. II: Validation of a
general purpose dose estimation code [In Process Citation] Med. Phys. 26 1396–403
[26] Johnson T K, McClure D and McCourt S 1999 MABDOSE. I: Characterization of a
general purpose dose estimation code [In Process Citation] Med. Phys. 26 1389–95
[27] Weber D, Eckerman K, Dillman L and Ryman J 1989 MIRD: Radionuclide Data and
Decay Schemes (New York: Society of Nuclear Medicine)
[28] Beierwaltes W H, Nishiyama R H, Thompson N W, Copp J E and Kubo A 1982
Survival time and ‘cure’ in papillary and follicular thyroid carcinoma with distant
metastases: statistics following University of Michigan therapy J. Nucl. Med. 23
561–8
[29] Benua R, Cicale N and Sonenberg M 1962 The relation of radioiodine dosimetry to
results and complications in the treatment of metastatic thyroid cancer Am. J.
Roentgenol. 87 171–82
[30] Furhang E E, Larson S M, Buranapong P and Humm J L 1999 Thyroid cancer
dosimetry using clearance fitting J. Nucl. Med. 40 131–6
[31] Marinelli L, Quimby E and Hine G 1948 Dosage determination with radioactive
isotopes. II. Practical considerations in therapy and protection Am. J. Roentgenol.
Radium Therapy 59 260–80
[32] Marinelli L, Quimby E and Hine G 1948 Dosage determination with radioactive
isotopes. I. Fundamental dosage formulae Nucleonics 2 56–66
104 Computational methods in internal radiation dosimetry
[33] Marinelli L, Quimby E and Hine G 1948 Dosage determination with radioactive
isotopes. II. Biological considerations and practical applications Nucleonics 2
449
[34] Loevinger R, Holt J G and Hine G J 1956 Internal administered radioisotopes. In
Radiation Dosimetry eds G J Hine and G L Brownell (New York: Academic)
pp 80173
[35] Zanzonico P B 2001 Radiation dosimetry of rhenium-188 Lipiodol: a workshop. In
Lecture Notes: Regional Training Workshop on Radionuclide Treatment of Liver
Cancer. NMS/11/RTC-2001-01/RCA-1 (Perth) (Perth: International Atomic Energy
Agency (IAEA)) pp 12861
[36] Erdi A K, Wessels B W, DeJager R, Erdi Y E, Atkins F B, Yorke E D, Smith L, Huang
E, Smiddy M, Murray J et al. 1994 Tumor activity confirmation and isodose curve
display for patients receiving iodine-131-labeled 16.88 human monoclonal antibody
Cancer 73 93244
[37] Erdi A K, Yorke E D, Loew M H, Erdi Y E, Sarfaraz M and Wessels B W 1998 Use of
the fast Hartley transform for three-dimensional dose calculation in radionuclide
therapy Med. Phys. 25 222633
[38] Erdi A K, Erdi Y E, Yorke E D and Wessels B W 1996 Treatment planning for radio-
immunotherapy Phys. Med. Biol. 41 200926
[39] Furhang E E, Chui C S and Sgouros G 1996 A Monte Carlo approach to patient-
specific dosimetry Med. Phys. 23 15239
[40] Furhang E E, Chui C S, Kolbert K S, Larson S M and Sgouros G 1997
Implementation of a Monte Carlo dosimetry method for patient-specific internal
emitter therapy Med. Phys. 24 116372
[41] Giap H B, Macey D J, Bayouth J E and Boyer A L 1995 Validation of a dose-
point kernel convolution technique for internal dosimetry Phys. Med. Biol. 40
36581
[42] Giap H B, Macey D J and Podoloff D A 1995 Development of a SPECT-based three-
dimensional treatment planning system for radioimmunotherapy J. Nucl. Med. 36
188594
[43] Howell R W, Rao D V and Sastry K S 1989 Macroscopic dosimetry for radioimmuno-
therapy: nonuniform activity distributions in solid tumors Med. Phys. 16 6674
[44] Humm J L 1986 Dosimetric aspects of radiolabeled antibodies for tumor therapy
J. Nucl. Med. 27 14907
[45] Humm J L and Cobb L M 1990 Nonuniformity of tumor dose in radioimmuno-
therapy [see comments] J. Nucl. Med. 31 7583
[46] Koral K F, Zasadny K R, Kessler M L, Luo J Q, Buchbinder S F, Kaminski M S,
Francis I and Wahl R L 1994 CT–SPECT fusion plus conjugate views for
determining dosimetry in iodine-131-monoclonal antibody therapy of lymphoma
patients J. Nucl. Med. 35 1714-20
[47] Koral K F, Lin S, Fessler J A, Kaminski M S and Wahl R L 1997 Preliminary results
from intensity-based CT–SPECT fusion in I-131 anti-B1 monoclonal-antibody
therapy of lymphoma Cancer 80 253844
[48] Leichner P K and Kwok C S 1993 Tumor dosimetry in radioimmunotherapy:
methods of calculation for beta particles Med. Phys. 20 52934
[49] Liu A, Williams L, Wong J and Raubitschek A 1997 A voxel source kernel (VSK)
method for rapid, patient-specific dose estimates in radioimmunotherapy (RIT)
(Abstract) J. Nucl. Med. 38(suppl) 106P
References 105
5.1. INTRODUCTION
108
Introduction 109
In 1966, Fisher and Snyder, at the Oak Ridge National Laboratory (ORNL),
reported on the development of an adult phantom for use in radiation
absorbed dose calculations [8]. The adult phantom was assumed to be
standing erect with the arms at the sides of the body. Three specific regions
were defined: the head and neck, the trunk including the arms, and the legs.
The head and neck were represented by a 14 cm 20 cm elliptical cylinder
with a height of 24 cm. The trunk and arms were modelled as a larger elliptical
cylinder, 20 cm 40 cm with a height of 70 cm. The legs below the buttocks
were modelled as a truncated elliptical cone with a height 80 cm. Regions of
little dosimetric concern were not included, e.g., the hands, feet, ears, nose,
etc. The composition of the phantom was assumed to be tissue distributed
homogeneously throughout. No attempt was made to model the lungs or
skeleton or to define the locations of specific organs in the phantom. Approxi-
mately 120 sub-regions were defined in the phantom that were used to assign
approximate values of the absorbed doses to organs located within specific
regions. In some cases, absorbed dose estimates for large organs required
the evaluation of the doses deposited in several of these regions.
Also included in this research was the development of phantoms for
humans of other ages. These phantoms represented humans with ages 0
(newborn), 1, 5, 10 and 15 years. These phantoms will be discussed in
more detail in section 5.2.6 below.
In 1967, Fisher and Snyder reported on the development of an adult
phantom with 22 internal organs and more than 100 sub-regions [9]. This
phantom represented the next step in the development of anthropomorphic
phantoms for use in dose calculations. However, this phantom also was
homogeneous and the skeleton and lung regions were not defined—in
terms of the assignment of different densities and elemental compositions
for these regions. This point was very misleading in that many of the early
drawings indicated such regions. However, a careful reading of the research
reports from the ORNL group will confirm the homogeneous composition.
Estimates of the absorbed dose to the skeleton were obtained by summing
all the energy deposited in the entire skeleton and dividing by the skeletal
volume. No special treatments were applied to obtain absorbed dose
estimates for the lungs. Calculations using the adult phantom agreed well
with those of Ellett et al. [10, 11] but were about 50% lower than those
obtained using the ICRP spherical organ methods [6].
Even though the original phantom was designed for use with internally
deposited radionuclides, Snyder saw many other applications. In addition, in
1967, he used the phantom to study the distribution of dose in the body from
external point sources of -rays [12]. He studied four photon energies (0.07,
0.15, 0.5 and 1.0 MeV) and four different source locations at distances of 1
and 2 m from the centre of the phantom.
112 Mathematical models of the human anatomy
Development of the adult human phantom by Snyder and his colleagues (see
section 5.2.3 above) was paralleled by the development of phantoms represent-
ing humans of other ages [8]. These phantoms represented children with ages
of 0 (newborn), 1, 5, 10 and 15 years. These early designs were assumed to have
outer dimensions that represented the average height, surface area and body
mass of a child of the particular age. All ‘paediatric phantoms’ were obtained
by applying a set of simple transformations to the axes of the Cartesian
coordinate system in which the adult phantom was located. These phantoms
became known as the ‘similitude phantoms’ because of their resemblance to
the adult. This approach has its limitations because children generally are
not just ‘little adults’. However, these phantoms were the first developed to
answer a real need in the nuclear medicine community [22].
Improvements in the paediatric models were closely linked with the
development of the heterogeneous adult phantom. Even though these new
phantoms were heterogeneous, the paediatric phantoms were obtained
through the same transformation method. That is, these phantoms also
were similitude phantoms. The outside dimensions were obtained by apply-
ing a series of transformations to the coordinate system and no consideration
was given to the actual organ sizes or shapes of the ‘transformed’ organs.
Although the masses of these transformed organs had to be known for the
calculation of absorbed dose, these masses were never published [23].
The limitations associated with transforming the major axes of the adult
phantom should be clear. Children are not simply small adults and their
organs are not necessarily ‘smaller adult organs’. Some organs are larger in a
child than in the adult and get smaller as the individual approaches adulthood,
e.g., the thymus. In addition, it was important that the actual size, shape,
location and mass of each organ be known for accurate dose calculations.
For these reasons, and others, a significant effort was undertaken at ORNL
during the mid-1970s to develop individual paediatric phantoms based upon
a careful review of the existing literature for each particular age. This effort
produced the next generation of mathematical phantoms that, although they
appeared to be modelled after the adult, were designed independently.
Three ‘individual phantoms’ were designed by Hwang et al. [24, 25].
These were the newborn, 1- and 5-year-old models. A separate effort was
undertaken by Jones et al. [26] for the 15-year-old and Deus and Poston
[27] undertook the design of a 10-year old after the other four designs were
complete. The development of the 10-year-old was significantly different
from those for the other four ages. In fact, this design was intended to
point the way to the next generation of more realistic phantoms. Even
though the design was completed and used for a limited number of dose
calculations, it was not popular because of the very complex geometry,
and other approaches to the development of phantoms were followed [28, 29].
The current stylized models 115
also presented new regions and improvements such as a new stylized breast
tissue region for all ages, the inclusion of the new model of the heart developed
by Coffey [42, 43], and a new model of the thyroid. While the ORNL paediatric
model series was initially published in 1981 [28], these models were not readily
utilized until 1987 with the publication of ORNL/TM-8381 [29]. In this
report, the only major change in the phantom series was that the 15-year-
old model was assumed also to represent the average adult female. For this
purpose, the breast, uterus and ovaries were modified according to published
reference average values [14]. The phantoms were used with the ETRAN
Monte Carlo photon transport code [44, 45] to calculate specific absorbed
fractions of energy in all five paediatric phantoms, as well as in the adult
male, for 12 photon energies (0.01 to 4.0 MeV). Electron transport was not
considered in these simulations and the electron energy was assumed to be
locally deposited.
models was less than that seen in the real anatomy (6–18 pyramids), it
represented a compromise between the mathematical simplicity needed for
Monte Carlo transport calculations, and the need for an improved anatomic
representation over the concentric ellipsoid-shell model of McAfee [53].
Each region was derived from dimensions and volumes given in ICRP Publi-
cation 23 [14] for both the newborn and the adult, and assumed constant
volume ratios between the different kidney sub-regions for the other ages. In
these models each medullary pyramid was modelled by half-ellipsoids (two
vertical and three horizontal) with the papillae at its tip, the pelvis by a portion
of an ellipsoid within the whole kidney, and the cortex was the remainder of
the kidney.
In both MIRD Pamphlets No 15 and 19, the EGS4 Monte Carlo trans-
port code [18, 59] was used for photon and electron transport. In these two
pamphlets, absorbed fractions of energy were tabulated for selected source
and target combinations (12 energies were simulated between 10 keV and
4 MeV). Following the MIRD method of internal dose calculation [30],
mean absorbed doses to the target regions per unit cumulated activity in
the source region (S values) were tabulated for specific radionuclides.
changing the position of the bladder and prostate gland to properly model
the relative positions of these organs. These models were developed not
only for the adult male model, but also for the other phantoms in the
ORNL series using physiological and anatomical descriptions published in
ICRP Publication 23 [14].
Because the intestinal tract and curvatures cannot be modelled with
simple geometrical shapes that can be easily coupled to a Monte Carlo trans-
port code, all models of this region have relied on thickening the wall region
to preserve anatomical wall and content mass. The critical cells for these
organs have been identified as the cells within the mucosal layer. In 1994,
Poston et al. developed a revision of the GI tract of the MIRD adult phan-
tom to better represent these sensitive targets [20, 21]. The actual wall of the
GI tract was divided in its thickness into four regions of varying radiosensi-
tivities; these layers were very difficult to model because the thickness of each
layer varied from one section to another along the different regions of the GI
tract. Poston et al. developed two methods to model this wall. The first
method divided the tissue wall into ten small, concentric layers (100 mm
thick for the adult) and the dose to each layer was recorded separately.
Then, the determination of the mucosal layer thickness for each section
would give directly the dose to the mucosa. However, since it is not possible
to determine directly the mucosal layer thickness for a specific patient, the
subdivision into 10 regions has not been used for medical dose calculation.
In a second method, Poston et al. measured the average thickness of the
mucosal layer along the GI tract from anatomic slides of cross sections of
a human GI tract. Different mucosal thicknesses for the stomach, the
small intestine and the large intestine were obtained. This layer was included
in the GI wall of the adult mathematical phantom and coupled to the EGS4
Monte Carlo transport code [18].
Stubbs et al. [62] presented calculations of the radiation-absorbed dose
to the walls of hollow organs. These authors studied all four sections of the
gastrointestinal tract but only for four radionuclides important to nuclear
medicine. They presented S values (in units of Gy/Bq s) for 90 Y, 99m Tc, 123 I
and 131 I. They concluded, as did Poston et al., that the ICRP ‘one-half’
assumption was overly conservative for non-penetrating radiation.
More recently, an improved (but very simple) geometric model for the
small intestine has been developed and used in a study of electron energy
deposition in the wall of the tract [63]. Results were obtained using this
model for monoenergetic electrons in the range 10 to 500 keV. Depth dose
profiles were developed so that dose to the more sensitive cells could be
assessed. Electron transport calculations using the computer code
MCNP4A [19] and a new model of the small intestine showed that only a
small fraction of the available energy reaches the critical cells in the crypts
of the wall. For electron energies below 330 keV, the only contribution to
absorbed dose to the stem cells came from bremsstrahlung. For higher
120 Mathematical models of the human anatomy
electron energies (i.e., >330 keV), the dose distribution curve was dominated
by monoenergetic electrons.
assumed, the dose from electrons scales as the inverse first power of the organ
mass (i.e., m1S ).
To better match a stylized model to a specific patient, models of
different sizes also can be used as proposed by Clairand et al. [68]. Using
the MIRD stylized model, they developed six new models of different
heights: three male models of heights 160, 170 and 180 cm and three
female models of heights 150, 160 and 170 cm. Using a statistical analysis
of anthropometric data gathered from 684 forensic autopsy cases for the
years 1987–1991 (all subjects selected were Caucasian adults who died of
external causes and showed no pathological changes), they derived organ
masses as a function of height [69]. Only organs weighed during autopsy
were taken into account in this study (heart, kidneys, liver, lungs, pancreas,
spleen, testes and thyroid), and other organs were obtained from the ORNL
models.
To support the dosimetry in therapeutic nuclear medicine, simple
tumour models also have been introduced in stylized phantoms [70–74].
These approaches allow for the determination of the photon dose from a
tumour to adjacent organs, but also provide the photon dose from body
organs to the tumour, which for some radionuclides can be non-negligible
[75]. Clairand et al. developed a computer program (DOSE3D) to simulate
spherical tumours inside the ORNL adult phantom [70]. With this program,
the user can select one or more tumours of different sizes, and the EGS4
Monte Carlo transport code is used for the particle transport. Similarly,
Johnson et al. developed the MABDOSE software to simulate spherical
tumours inside the ORNL phantom series [71–74]. In this code, the phantom
series is first voxelized at a user-selectable resolution (10, 5, 2, or 1 mm), and
the ETRAN Monte Carlo photon transport code used for the particle
transport.
recently reported on their efforts to merge this data set with both the MCNP
and EGS4 radiation transport codes [86]. Computational results have been
reported for internal electron sources [87], external photon sources [88],
neutrons [89, 90] and electrons [91].
At present, three tomographic computational models have been
constructed for specific use in paediatric dosimetry studies [92, 93]. Two of
these are the BABY and CHILD models developed at GSF in Germany
by Zankl and colleagues. BABY is constructed from a whole-body CT
exam of an 8-week-old female 24 h after death. CHILD was constructed
from a CT exam of a live 7-year-old female who was undergoing subsequent
whole-body irradiation prior to bone marrow transplantation. Image
segmentation was performed which assigned voxels to 54 regions in the
baby and 64 regions in the child. The models have been used to determine
paediatric organ doses in a variety of medical exposure scenarios including
nuclear medicine [94], projection radiography [95–98] and computed tomo-
graphy [99, 100]. The third paediatric model, ADELAIDE, is based upon
44 consecutive CT scans of a 14-year-old female [101]. This particular
model is restricted to the torso only, and has been used with the EGS4
transport code to assess organ doses received during CT exams.
Recent efforts at GSF have resulted in a new adult reference male
tomographic model, GOLEM. This model was segmented from whole-
body medical images of a living 38-year-old male with external dimensions
close to those of the ICRP Reference Man [102–105]. Recent German–
Japanese collaborations also have been undertaken to establish a reference
tomographic computational model for the adult Asian male [106].
5.5. SUMMARY
REFERENCES
[1] NCRP 1985 The experimental basis for absorbed-dose calculations in medical uses of
radionuclides NCRP Report No 83 (Bethesda, MD: National Council on Radiation
Protection and Measurements)
[2] Marinelli L D 1942 Dosage determination with radioactive isotopes Am. J.
Roentgenol. Rad. Therapy 47 210–16
[3] Marinelli L D, Quimby J H and Hine G J 1948 Dosage determination with
radioactive isotopes II. Practical considerations in therapy and protection Am. J.
Roentgenol. Rad. Therapy 59 260–81
[4] Loevinger R, Japha E M and Brownell G L 1965 Discrete radioisotope sources.
In Radiation Dosimetry eds Hine G J and Brownell G L (New York: Academic)
pp 693–799
[5] Loevinger R, Hiolt J G and Hine G J 1965 Internally administered radionuclides.
In Radiation Dosimetry eds Hine G J and Brownell G L (New York: Academic)
pp 801–73
[6] ICRP 1959 Report of committee II on permissible dose for internal radiation
International Commission on Radiological Protection (Oxford: Pergamon)
[7] Brownell G L, Ellett W H and Reddy A R 1968 MIRD Pamphlet No 3: Absorbed
fractions for photon dosimetry J. Nucl. Med. Suppl 29–39
[8] Fisher H L J and Snyder W S 1966 Variation of dose delivered by 137Cs as a function of
body size from infancy to adulthood ORNL-4007 (Oak Ridge, TN: Oak Ridge
National Laboratory) pp 221–28
References 127
[9] Fisher H L J and Snyder W S 1967 Distribution of dose in the body from a source of
gamma rays distributed uniformly in an organ ORNL-4168 (Oak Ridge, TN: Oak
Ridge National Laboratory) pp 24557
[10] Ellet W H, Callahan A B and Brownell G H 1964 Gamma-ray dosimetry of internal
emitters: Monte Carlo calculations of absorbed dose from point sources Brit. J.
Radiol. 37 4552
[11] Ellet W H, Callahan A B and Brownell G H 1965 Gamma-ray dosimetry of internal
emitters II: Monte Carlo calculations of absorbed dose from uniform sources Brit. J.
Radiol. 38 54144
[12] Snyder WS 1967 The variation of dose in man from exposure to a point source of
gamma rays ORNL-4168 Oak Ridge National Laboratory Oak Ridge, TN 25761
[13] Snyder W S, Ford M R, Warner G G and Fisher Jr H L 1969 MIRD Pamphlet No 5:
Estimates of absorbed fractions for monoenergetic photon sources uniformly
distributed in various organs of a heterogeneous phantom J. Nucl. Med. 10 152
[14] ICRP 1975 Report on the task group on Reference Man ICRP Publication 23
(Oxford: International Commission on Radiological Protection)
[15] ICRP 1979 Limits for intakes of radionuclides by workers Publication 30.
International Commission on Radiological Protection (Oxford: Pergamon)
[16] Snyder W S, Ford M R and Warner G G 1978 MIRD Pamphlet No 5, Revised:
Estimates of specific absorbed fractions for photon sources uniformly distributed in
various organs of a heterogeneous phantom (New York: Society of Nuclear Medicine)
[17] Warner G G and Craig A N J 1968 ALGAM: A computer program for estimating
internal dose in a man phantom ORNL-TM-2250 (Oak Ridge, TN: Oak Ridge
National Laboratory)
[18] Nelson W R, Hirayama R H and Roger D W O 1985 The EGS4 code system SLAC
Report 265 (Stanford, CA: Stanford Linear Accelerator Center)
[19] Briesmeister J F 1993 MCNP—a general Monte Carlo N-particle transport code,
version 4A LA-12625-M (Los Alamos, NM: Los Alamos National Laboratory)
[20] Poston J W Jr, Kodimer K A, Bolch W E and Poston J W Sr 1996 Calculation of
absorbed energy in the gastrointestinal tract Health Phys. 71 3006
[21] Poston J W Jr, Kodimer K A, Bolch W E and Poston J W Sr 1996 A revised model
for the calculation of absorbed energy in the gastrointestinal tract Health Phys. 71
30714
[22] Kereiakes J G, Seltzer R A, Blackburn B and Saenger E L 1965 Doses to infants and
children: a plea for a standard child Health Phys. 11 9991004
[23] Poston J W 1990 The development of early pediatric models and their application to
radiation absorbed dose calculations. In Dosimetry of Administered Radionuclides
eds S J Adelstein, A I Kassis and R.W.B) (Washington, DC: American College of
Nuclear Physicians) pp 10521
[24] Hwang J M L, Shoup R L, Warner G G and Poston J W 1976 Mathematical
description of a one- and five-year-old child for use in dosimetry calculations
ORNL/TM-5293 (Oak Ridge, TN: Oak Ridge National Laboratory)
[25] Hwang J M L, Shoup R L and Poston J W 1976 Mathematical description of a
newborn human for use in dosimetry calculations ORNL/TM-5453 (Oak Ridge,
TN: Oak Ridge National Laboratory)
[26] Jones R M, Poston J W, Hwang J L, Jones T D and Warner G G 1976 The develop-
ment and use of a fifteen-year-old equivalent mathematical phantom for internal dose
calculations ORNL/TM-5278 (Oak Ridge, TN: Oak Ridge National Laboratory)
128 Mathematical models of the human anatomy
[45] Seltzer S M 1991 Electronphoton Monte Carlo calculations: the ETRAN code
Appl. Radiat. Isot. 42 91741
[46] Bouchet L G, Bolch W E, Weber D A, Atkins H L and Poston J W Sr 1999 MIRD
Pamphlet No 15: Radionuclide S values in a revised dosimetric model of the adult
head and brain J. Nucl. Med. 40 62S101S
[47] London E D 1993 Imaging Drug Action in the Brain (Boca Raton, FL: CRC)
[48] van Heertum R L and Tikofsky R S 1995 Cerebral SPECT Imaging (New York:
Raven)
[49] Blinkov S M and Glezer I I 1968 The Human Brain in Figures and Tables (New York:
Plenum)
[50] Petegnief Y, Petiet A, Peker M C, Bonnin F, Meulemans A and Le Guludec D 1998
Quantitative autoradiography using a radioimager based on a multiwire propor-
tional chamber Phys. Med. Biol. 43 362938
[51] Willis K W, Martinez D A, Hedley-Whyte E T, Davis M A, Judy P F and Treves S
1977 Renal localization of 99mTc-stannous glucophetonate and 99mTc-stannous
dimercaptosuccinate in the rat by frozen section autoradiography. The efficiency
and resolution of technetium-99m Radiat. Res. 69 47588
[52] Snyder W S and Ford M R 1967 A dosimetric study for the administration of
neohydrin labeled with 203Hg and 197Hg ORNL-4168 (Oak Ridge, TN: Oak Ridge
National Laboratory) pp 26773
[53] McAfee J G 1969 Problems in evaluating the radiation dose for radionuclides
excreted by the kidneys. In Medical Radionuclides: Radiation Dose and Effects eds
R J Cloutier, C L Edwards and W S Snyder (Oak Ridge, TN: US Atomic Energy
Commission) pp 27194
[54] Blau M, McAfee J G, Roher R H, Snyder W S and Smith E M 1975 MIRD dose
estimate No. 6: Hg-197- and Hg-203-labeled chlormerodrin J. Nucl. Med. 16 121417
[55] Patel J S 1988 A revised model of the kidney for medical internal radiation dose
(College Station, TX: Texas A&M University)
[56] Hui T Y E 1985 A preliminary model of the circulating blood for use in radiation
dose calculations (College Station, TX: Texas A&M University)
[57] Bouchet L G, Bolch W E, Blanco P, Rajon D A, Clairand I, Konijnnenberg M W
and Wessels B W 2002 MIRD Pamphlet No 19: Absorbed fractions of energy and
radionuclide S values for six pediatric multiregion models of the kidney (in press)
[58] Boerman O C, Oyen W J and Corstens F H 2001 Between the Scylla and Charybdis
of peptide radionuclide therapy: hitting the tumor and saving the kidney Eur. J.
Nucl. Med. 28 14479
[59] Ford R L and Nelson W R 1978 The EGS code system: computer programs for Monte
Carlo simulation of electromagnetic cascade shower Report 210 (Stanford: Stanford
Linear Accelerator Center)
[60] Stabin M 1994 A model of the prostate gland for use in internal dosimetry J. Nucl.
Med. 35 51620
[61] Mardirossian G, Tagesson M, Blanco P, Bouchet LG, Stabin M, Yoriyaz H, Baza S,
Ljungberg M, Strand S E and Brill A B 1999 A new rectal model for dosimetry
applications J. Nucl. Med. 40 152431
[62] Stubbs J B, Evans J F and Stabin M 1998 Radiation absorbed doses to the walls of
hollow organs J. Nucl. Med. 39 198995
[63] Bhuiyan N U 2000 A revised dosimetric model for calculation of electron dose in the
small intestine (College Station, TX: Texas A&M University)
130 Mathematical models of the human anatomy
[64] Watson E E, Stabin M, Davis J L and Eckerman K F 1989 A model of the peritoneal
cavity for use in internal dosimetry J. Nucl. Med. 30 200211
[65] Deloar H M, Watabe H, Nakamura T, Narita Y, Yamadera A, Fujiwara T and Itoh
M 1997 Internal dose estimation including the nasal cavity and major airway for
continuous inhalation of C15O2, 15O2 and C15O using the thermoluminescent
dosimeter method J. Nucl. Med. 38 160313
[66] Clairand I, Bouchet L G and Bolch W E 2000 A new macroscopic model of the long
bones for skeletal dosimetry J. Nucl. Med. 41 1062
[67] Snyder W S, Ford M R, Warner G G and Watson S B 1975 MIRD Pamphlet No 11:
‘S,’ absorbed dose per unit cumulated activity for selected radionuclides and organs
(New York: Society of Nuclear Medicine)
[68] Clairand I, Bouchet L G, Ricard M, Durigon M, Di Paola M and Aubert B 2000
Improvement of internal dose calculations using mathematical models of different
adult heights Phys. Med. Biol 45 277185
[69] de la Grandmaison G L, Clairand I and Durigon M 2001 Organ weight in 684 adult
autopsies: new tables for a Caucasoid population Forensic Sci. Int. 119 14954
[70] Clairand I, Ricard M, Gouriou J, Di Paola M and Aubert B 1999 DOSE3D: EGS4
Monte Carlo code-based software for internal radionuclide dosimetry J. Nucl. Med.
40 151723
[71] Johnson T K, McClure D L, McCourt S L, Andl G J and Berman B D 1998
MABDOSE: A computer program for the calculation of dose. In Proceedings of the
Sixth International Radiopharmaceutical Dosimetry Symposium ed A T Schlafke-
Stelson (Oak Ridge, TN: ORAU)
[72] Johnson T K 1988 MABDOSE: A generalized program for internal radionuclide
dosimetry Comput. Meth. Progr. Biomed. 27 15967
[73] Johnson T K, McClure D and McCourt S 1999 MABDOSE. I: Characterization of a
general purpose dose estimation code Med. Phys. 26 138995
[74] Johnson T K, McClure D and McCourt S 1999 MABDOSE. II: Validation of a
general purpose dose estimation code Med. Phys. 26 1396403
[75] Johnson T K and Colby S B 1993 Photon contribution to tumor dose from
considerations of 131 I radiolabeled antibody uptake in liver, spleen, and whole
body Med. Phys. 20 166774
[76] Kawrakow I 2000 Accurate condensed history Monte Carlo simulation of electron
transport. I. EGSnrc, the new EGS4 version Med. Phys. 27 48598
[77] Zubal I G, Harrell C R, Smith E O, Rattner Z, Gindi G and Hoffer P B 1994
Computerized 3-dimensional segmented human anatomy Med. Phys. 21 299302
[78] Zubal I G and Harell C R 1992 Voxel based Monte Carlo calculations of nuclear
medicine images and applied variance reduction techniques Image Vision Comput.
10 34248
[79] Zubal I G, Harrell C R, Smith E O, Smith A L and Krischlunas P 1996 High
resolution, MRI-based, segmented, computerized head phantom. In Proceedings
of the Sixth International Radiopharmaceutical Dosimetry Symposium eds A T S-
Stelson, M Stabin and R B Sparks (Oak Ridge, TN: ORAU) pp 31924
[80] Stabin M, Yoriyaz H, Brill A B and Dawant M 1999 Monte Carlo calculation of
dose conversion factors for a new generation of dosimetry phantoms J. Nucl.
Med. 40 310P11P
[81] Yoriyaz H, dos Santos A, Stabin M and Cabezas R 2000 Absorbed fractions in a
voxel-based phantom calculated with the MCNP-4B code Med. Phys. 27 155562
References 131
[82] Stabin M and Yoriyaz H 2001 Photon specific absorbed fractions calculated in the
trunk of an adult male voxel-based phantom Health Phys. 82 2144
[83] Tagesson M, Zubal I G, Ljungberg M and Strand S E 1996 S values for subregions in
the brain. In Proceedings of the Sixth International Radiopharmaceutical Dosimetry
Symposium eds A T S-Stelson, M Stabin and R B Sparks (Oak Ridge, TN:
ORAU) pp 32530
[84] Jones D G 1997 A realistic anthropomorphic phantom for calculating organ doses
arising from external photon irradiation Radiat. Prot. Dosim. 72 2129
[85] Jones D G 1998 A realistic anthropomorphic phantom for calculating specific
absorbed fractions of energy deposited from internal gamma emitters Radiat.
Prot. Dosim. 79 41114
[86] Xu X G, Chao T C and Bozkurt A 2000 VIP-Man: an image-based whole-body adult
male model constructed from color photographs of the Visible Human Project for
multi-particle Monte Carlo calculations Health Phys. 78 47686
[87] Chao T C and Xu X G 2001 Specific absorbed fractions from the image-based VIP-
Man body model and EGS4-VLSI Monte Carlo code: internal electron emitters
Phys. Med. Biol 46 90127
[88] Chao T C, Bozkurt A and Xu X G 2001 Conversion coefficients based on the VIP-
Man anatomical model and EGS4 Health Phys. 81 16383
[89] Bozkurt A, Chao T C and Xu X G 2000 Fluence-to-dose conversion coefficients
from monoenergetic neutrons below 20 MeV based on the VIP-man anatomical
model Phys. Med. Biol. 45 305979
[90] Bozkurt A, Chao T C and Xu X G 2001 Fluence-to-dose conversion coefficients
based on the VIP-Man anatomical model and MCNPX code for monoenergetic
neutrons above 20 MeV Health Phys. 81 184202
[91] Chao T C, Bozkurt A and Xu X G 2001 Organ dose conversion coefficients for
0.110 MeV electrons calculated for the VIP-Man tomographic model Health
Phys. 81 20314
[92] Zankl M, Veit R, Petoussi N, Mannweiler E, Wittmann A and Drexler G 1994
Realistic computerized human phantoms Life Sci. Space Res. 14 42331
[93] Zankl M, Veit R, Williams G, Schneider K, Fendel H, Petoussi N and Drexler G
1988 The construction of computer tomographic phantoms and their application
in radiology and radiation protection Radiat. Environ. Biophys. 27 153–64
[94] Petoussi-Henß N and Zankl M 1998 Voxel anthropomorphic models as a tool for
internal dosimetry Radiat. Prot. Dosim. 79 41518
[95] Petoussi-Henß N, Panzer W, Zankl M and Drexler G 1995 Dose–area product and
body doses Radiat. Prot. Dosim. 57 36366
[96] Veit R and Zankl M 1992 Influence of patient size on organ doses in diagnostic
radiology Radiat. Prot. Dosim. 43 24143
[97] Veit R and Zankl M 1993 Variation of organ doses in pediatric radiology due to
patient diameter calculated with phantoms of varying voxel size Radiat. Prot.
Dosim. 49 35356
[98] Zankl M 1993 Computational models employed for dose assessment in diagnostic
radiology Radiat. Prot. Dosim. 49 33944
[99] Kalender W A, Schmidt B, Zankl M and Schmidt M 1999 A PC program for estimating
organ dose and effective dose values in computed tomography Eur. Radiol. 9 55562
[100] Zankl M, Panzer W, Petoussi-Henß N and Drexler G 1995 Organ doses for children
from computer tomographic examinations Radiat. Prot. Dosim. 57 39396
132 Mathematical models of the human anatomy
6.1. INTRODUCTION
The Monte Carlo method is widely used for solving problems involving
statistical processes and is very useful in medical physics due to the stochastic
nature of radiation emission, transport and detection processes. The method
is very useful for complex problems that cannot be modelled by computer
codes using deterministic methods or when experimental measurements
may be impractical. However, due to computer limitations, the method has
not yet fully lived up to its potential. With the advent of high-speed super-
computers, the field has received increased attention, particularly with
parallel algorithms and high performance platforms, which have much
higher execution rates. Recent advances in computer hardware and software
have lowered, if not removed, the barriers to wider acceptance of Monte
Carlo calculations by life scientists.
Many Monte Carlo programs have been in use in the field of nuclear
imaging, internal dosimetry and treatment planning with many of them
available in the public domain [1, 2]. Those Monte Carlo software
packages are excellent research tools; however, becoming familiar with
the codes can be time consuming, and most codes are not completely
free of bugs. Moreover, the simulation package must be integrated in
the local software environment. For that purpose, modifications may be
necessary, which require a detailed working knowledge of the code struc-
ture. This chapter summarizes state-of the-art developments of widely used
Monte Carlo codes in therapeutic nuclear medicine in connection with
computing facilities available for routine and research applications.
Current trends and some strategies for future development in the field
will also be discussed.
133
134 Monte Carlo codes for use in therapeutic nuclear medicine
The ‘Monte Carlo method’ describes a very broad area of science, in which
many processes, physical systems and phenomena are simulated by statistical
methods employing random numbers (see chapter 1). In the specific applica-
tion of interest here, the transport of ionizing radiation particles is simulated
by creation of particles or rays from a defined source region, generally with
random initial orientation in space, with tracking of the particles as they
travel through the system, sampling the probability density functions
(pdfs) for their interactions to evaluate their trajectories and energy deposi-
tion at different points in the system. The interactions determine the penetra-
tion and motion of particles but, more importantly, the energy deposited
during each interaction gives the radiation absorbed dose, when divided by
the appropriate values of mass. With sufficient numbers of interactions,
the mean absorbed dose at points of interest will be given with acceptable
uncertainties. The central issues include how well the real system of interest
can be simulated by a geometrical model, how many histories (i.e., how much
computer time) are needed to obtain acceptable uncertainties (usually
around 5%, no more than 10%), and how measured data can be used to
validate the theoretical calculations.
In the general sense, a problem that can be treated by a Monte Carlo
method is determined by a set of parameters that in aggregate determine
completely the behaviour of an element of the system in all situations
modelled within that system. Specifically, in applications of Monte Carlo
in therapeutic nuclear medicine, we are interested in modelling the creation
and interaction of all particles (speaking somewhat imprecisely, as photons,
which are of high interest in this situation, are best described as waves
rather than as particles) throughout their ‘lifetime’ in the body. Particles
are generally ‘created’ at random positions within a source region (which
might be an entire organ or some smaller structure), with defined initial
orientations, and then travel random distances before interacting with
atoms within the tissues of the body and depositing energy in these inter-
actions, and eventually being absorbed completely or escaping the body.
A complete characterization of the source will include its spatial distribu-
tion, the angular distribution of its emission spectrum (most often
random in 4 space, unless some unique situation exists which restricts
the orientation of the emissions—this is more common in situations
involving external sources, which may be collimated or shielded in certain
ways), and its energy distribution. Once a particle is created and begins its
journey away from the source location, the three-space location of the first
interaction must be chosen. This is done by sampling the known
probability distributions for its interactions with atoms in the body. The
probabilities of interaction for photons are very well known—photo-
electric, Compton, and pair production events occur with well defined
Historical developments 135
The first real treatment of internal dose assessment was done in the Manhattan
project, developed by K Z Morgan, W Snyder, and others [4]. Their focus was
on the protection of radiation workers and, to a degree, the general public,
from the production and use of radioactive materials, in particular in relation
to the war effort. They did not, however, employ Monte Carlo methods in
these calculations. Their model of the human body and each of its parts was
a simple homogeneous spherical construct, in which electrons were completely
absorbed and for which photon absorption could be modelled by a simple
analytical expression, knowing the absorption coefficient for the organ or
tissue of interest. More specifically, for soft tissue, bone, lung, or any mixture
thereof (such as the whole body), the photon absorbed fraction for an organ
irradiating itself is given as (1 er ), where is the absorption coefficient
for the tissue composition of interest and r is the assigned effective radius of
the object. This information was developed for the very large number of radio-
nuclides of interest to radiation safety at that time—products of the fission
process, activation products from neutron interactions in the structural and
other materials in the reactors, and others. The results were given in the
context of a system of radiation protection guidelines for workers, i.e.,
amounts of permissible intake of each radionuclide per working year (given
permissible dose limits), and permissible concentrations in air and water
(given the intake limits and assumptions about standard breathing and
drinking rates) [4].
The first breakthrough in the use of Monte Carlo methods was the
development of the Fisher–Snyder heterogeneous, hermaphrodite, anthropo-
morphic model of the human body in the 1970s [5]. This model, or
‘phantom’, consisted of spheres, ellipsoids, cones, tori, and subsections of
such objects, combined to approximate the geometry of the body and its
internal structures. For the majority of cases, electron energy was assumed
to be absorbed where created (exceptions included hollow organs with
separate wall and contents sections, and the bone and marrow, in which
‘crossfire’ between adjacent regions can occur due to the small dimensions
of the regions). For photons, Monte Carlo methods were developed using
a computer code called ALGAM [6], which created photons at random posi-
tions within any ‘source’ region (organ or tissue assumed to be contaminated
with radioactivity), gave these photons a random orientation in 4 space, and
then followed them through various Compton and photoelectric interactions
(pair productions events were quite rare, as starting energies did not exceed
4 MeV) until the photon reached a certain critical low (‘cut-off ’) energy and
was assumed to be locally absorbed, or until it escaped the surface of the
body (at which point the probability of scatter from an air molecule and
redirection towards the body was assumed to be negligibly low). With
repeated sampling of the source, which in this time generally involved only
Historical developments 137
subregions of the brain [9, 10], the eye [11], the peritoneal cavity [12], the
prostate gland [13], bone and marrow [14, 15], rectum [16], spheres of varying
size [17, 18], and others.
product of that committee), and partly to integrate the name into a new
unified system of internal and external dose assessment. This unified
system is deployed on an Internet website for rapid electronic access [25].
This site, called the RAdiation Dose Assessment Resource (RADAR),
provides decay data for over 800 radionuclides, absorbed fractions for all
available stylized phantoms and some voxel phantoms, kinetic data, dose
factors (for all phantoms and nuclides), risk information, and other data
via electronic transfer to users worldwide.
Table 6.1 lists widely used public domain Monte Carlo codes together with a
short description of their key features. Most of the packages mentioned
below run virtually on different platforms and operating systems and are
available free of charge from the authors or through the official channels
(RSIC or NEA).
Table 6.1. Key features of public domain Monte Carlo codes used in therapeutic
nuclear medicine applications.
ALGAM [6] Monte Carlo calculations using the OGRE system are employed for
photon transport. The user can specify a source arbitrarily through a
subroutine. Simulation for internal dosimetry calculations is
specifically included.
MIRDOSE [19] Dose calculations based on Monte Carlo pre-calculated specific
absorbed fractions for 10 phantom models using the ALGAMP code.
EGS4 [26] Coupled photon/electron transport in any material through
user-specified geometries. Simulation for internal dosimetry
calculations is not specifically included and requires an extensive
amount of user programming in Mortran.
ITS including Coupled photon/electron transport in any material through slabs,
TIGER cylinders or combinatorial. Simulation for internal dosimetry
CYLTRAN and calculations is not specifically included and requires an extensive
ACCEPT [27] amount of user programming in Fortran.
MCNP [28] Coupled neutron/photon/electron transport in any material through
user-generalized geometry. Simulation for internal dosimetry
calculations is not specifically included and requires an extensive
amount of user manipulation of input date files to model complex
geometries.
GEANT [29] Coupled photon/electron transport in any material through
combinatorial geometry. Simulation for internal dosimetry
calculations not specifically included and requires an extensive amount
of user programming in Cþþ (GEANT 4).
ETRAN [30] Coupled photon/electron transport in plane-parallel slab targets that
have a finite thickness in one dimension and are unbound in the other
two dimensions. Simulation for internal dosimetry calculations is not
specifically included and requires an extensive amount of user
manipulation of input date files to model complex geometries.
6.3.4. ITS
The Integrated TIGER Series (ITS) of coupled electron/photon Monte Carlo
transport codes is a powerful tool for determining state-of-the-art descriptions
of the production and transport of the electron/photon cascade in time-
independent, multi-material, multi-dimensional environments [27]. ITS is a
collection of programs sharing a common source code library that can solve
sophisticated radiation transport problems. A total of eight codes are in the
collection which can be split into six groups: the TIGER codes (for one-
dimensional slab geometries), the CYLTRAN codes (for two-dimensional
cylindrical geometries), the ACCEPT codes (for arbitrary three-dimensional
geometries), the standard codes (for normal applications), the P codes (for
applications where enhanced ionization/relaxation procedures are needed),
and the M codes (for applications which involve two- or three-dimensional
macroscopic electromagnetic fields). The user selects the appropriate code
from the library and supplies it with any special requirements and the physical
description of the problem to be solved in an input file.
6.4.1. Introduction
Currently, and for the past 20–30 years, internal dose calculations for nuclear
medicine have used the standardized hermaphroditic phantoms with stylized
geometries (spheres, ellipsoids, cylinders, etc.) described above to represent
the major organs, skeletal system, whole body, and bone marrow of ‘Refer-
ence Man’ (70 kg), ‘Reference Woman’ (57 kg), and ‘Reference’ children and
pregnant women. Doses calculated using these models do not include the
detail needed for accurate assessment of risks in therapeutic administrations
of radiopharmaceuticals and their results have been shown to be poor predic-
tors of marrow radiotoxicity. These models give only average dose to whole
organs (not dose distributions within organs or tissues with possibly nonuni-
form activity distributions) and the reported dose is applicable only to a
Limitations of current nuclear medicine dose calculations 143
99m
Table 6.2. Radiation dose estimates for the reference adult for Tc
glucoheptonate.
person whose size and weight are close to that of the reference individual
after which the model was derived. Thus the doses reported with such
models really represent the dose to a phantom, not to a patient. If the bioki-
netic data to be applied were taken from the actual patient, then these data
are patient specific. In diagnostic applications in nuclear medicine, usually a
standardized kinetic model is also applied. Typically results are shown in
table 6.2.
These dose estimates are based on a standard kinetic model for gluco-
heptonate, and AFs for the adult male phantom in the Cristy/Eckerman
phantom series. If one were to use patient-specific biokinetic data for a
nuclear medicine therapy agent with this same phantom, the result would
be the same as in table 6.2 (except that effective dose quantities may not be
used in therapy applications). One can make the dose estimates more patient
specific through mass-based adjustments to the organ doses:
. Absorbed fractions for electrons and alphas scale linearly with mass.
. Absorbed fractions for photons scale with mass to a power of 1/3.
One generally cannot:
. Account for patient-specific differences in organ geometry.
. Account for patient-specific marrow characteristics.
. Calculate dose distributions within organs.
To perform real patient-specific dose calculations, one needs a patient-
specific physical model to be used with patient-specific biokinetic data. A
‘one-dose-fits-all’ approach to radiation therapy with these internal emitter
treatments is not likely to be effective (due to the narrow range between
tumour ablation and bone marrow toxicity). Individual patients not only
have significantly different uptake and retention half-times of activity of
144 Monte Carlo codes for use in therapeutic nuclear medicine
6.5.1.3.4. Juweid et al. [49] treated patients with 131 I-labelled anti-carcino-
embryonic antigen monoclonal antibodies with administered activities based
either on the patient’s body surface area or on a prescribed red marrow
dose as determined by a pre-therapy dosimetric study.
normal organs, but do not at present use patient images. These codes
work with stylized representations of average individuals, and give average
dose to whole organs. The RTDS code [52] employs either the standard
MIRDOSE phantom set (pre-calculated dose conversion factors for stan-
dard phantoms) or its own transport algorithms in a limited body space,
based on voxel source kernels (which are useful in a single tissue type) to
produce average organ doses or dose distributions within specified organs
or tissues of the body.
Figure 6.1. Images showing the GI tract (top) and kidney–liver region
(bottom) from the Yale voxel phantom [62] (left) and the Cristy/Eckerman
geometrical phantom [7] (right).
and hence the result computed by a program can vary from one run to the
next. In the replicated approach, multiple instances of the same generator
are created (for example, one per task). Each generator uses either the
same seed or a unique seed, derived, for example, from a task identifier.
Clearly, sequences generated in this fashion are not guaranteed to be
independent and, indeed, can suffer from serious correlation problems.
However, the approach has the advantages of efficiency and ease of
implementation and should be used when appropriate. In the distributed
approach, responsibility for generating a single sequence is partitioned
among many generators, which can then be parcelled out to different
tasks. The generators are all derived from a single generator; hence, the
analysis of the statistical properties of the distributed generator is simplified.
In a review of vectorized Monte Carlo, Martin and Brown [70] described
variations of event-based algorithms together with speed-up results
published by different groups. During the past two decades, investigations
were carried out to run different Monte Carlo codes on multiple-transputer
systems [67], vector parallel supercomputers [71], parallel computers [66] and
a cluster of workstations in a local area network using PVM [72]. There are
large discrepancies in the performance ratio reported by different authors. In
particular, Miura reported a speed-up of about 8 with the vectorized ESG4
code (EGS4V) [73]. A linear decrease in computing time with the number of
processors used was also demonstrated for eight [66] and up to 32 processors
[74]. This speed-up is especially significant in SPECT simulations involving
higher energy photon emitters, where explicit modelling of the phantom
and collimator is required.
6.6. SUMMARY
It is gratifying to see in overview the progress that internal dose modelling has
made, from simple spherical geometries, through stylized models of reference
individuals, and most recently towards truly patient-specific models based on
patient images. Challenges remain, particularly in the areas of bone and
marrow dosimetry, in image quantification, and in correlating radiation dose
with effect. As these challenges are met, and experience is gained, patients can
expect to receive therapy with internal emitters, based on accurate dose calcula-
tions and sound risk/benefit evaluations made by their attending physician.
ACKNOWLEDGMENTS
One of the authors (HZ) acknowledges the support of grant SNSF 3152-
062008 from the Swiss National Science Foundation. The contents and
views in this chapter are solely the responsibility of the author and do not
necessarily represent the official views of the institution providing support.
REFERENCES
[1] Andreo A 1991 Monte Carlo techniques in medical radiation physics Phys. Med. Biol.
36 861–920
[2] Zaidi H 1999 Relevance of accurate Monte Carlo modeling in nuclear medical
imaging Med. Phys. 26 574–608
[3] Bielajew A F and Rogers D W O 1989 Monte Carlo Transport of Electrons and
Photons eds N W Jenkins, A Rindi, A Nahum and D W O Rogers (New York:
Plenum Press) pp 407–19
[4] ICRP 1959 Report of Committee II on Permissible Dose for Internal Radiation. ICRP
Publication 2 I.C.o.R. Protection (New York: Pergamon)
[5] Snyder W, Ford M R and Warner G 1978 Estimates of specific absorbed fractions for
photon sources uniformly distributed in various organs of a heterogeneous phantom.
MIRD Pamphlet No 5, revised NM/MIRD (New York: Society of Nuclear Medicine)
[6] Warner G G and Craig A M J 1968 ALGAM, a computer program for estimating
internal dose for gamma-ray sources in a man phantom ORNL-TM-2250 (Oak
Ridge National Laboratory)
[7] Cristy M and Eckerman K F 1987 Specific absorbed fractions of energy at various ages
from internal photon sources. I Methods, II one year old, III five year old, IV ten year
old, V fifteen year old male and adult female, VI new-born and VII adult male ORNL/
TM 8381/V1-V7 (Oak Ridge National Laboratory)
[8] Stabin M et al. 1995 Mathematical models and specific absorbed fractions of photon
energy in the nonpregnant adult female and at the end of each trimester of pregnancy
ORNL/TM-12907 (Oak Ridge National Laboratory)
[9] Eckerman K, Cristy M and Warner G G 1981 Dosimetric Evaluation of Brain Scanning
Agents eds E E Watson, A T Schlafke-Stelson, J L Coffey and R J Cloutier, Third
International Radiopharmaceutical Dosimetry Symposium; HHS Publication FDA
154 Monte Carlo codes for use in therapeutic nuclear medicine
81-8166, US Dept of Health and Human Services, Oak Ridge, TN (Rockville, MD:
Food and Drug Administration) pp 527–40
[10] Bouchet L G, Bolch W E, Weber D A, Atkins H L and Poston J W 1999 MIRD
Pamphlet No 15: Radionuclide S values in a revised dosimetric model of the adult
head and brain. Medical Internal Radiation Dose J. Nucl. Med. 40 62S–101S
[11] Holman B L et al. 1983 Biodistribution and dosimetry of N-isopropyl-p-[123I]iodo-
amphetamine in the primate J. Nucl. Med. 24 922–31
[12] Watson E E, Stabin M, Davis J L and Eckerman K F 1989 A model of the peritoneal
cavity for use in internal dosimetry J. Nucl. Med. 30 2002–11
[13] Stabin M 1994 A model of the prostate gland for use in internal dosimetry J. Nucl.
Med. 35 516–20
[14] Eckerman K F and Stabin M 2000 Electron absorbed fractions and dose conversion
factors for marrow and bone by skeletal regions Health Phys. 78 199–214
[15] Bouchet L G, Bolch W E, Howell R W and Rao D V 2000 S values for radionuclides
localized within the skeleton J. Nucl. Med. 41 189–212
[16] Mardirossian G et al. 1999 A new rectal model for dosimetry applications J. Nucl.
Med. 40 1524–31
[17] Siegel J A and Stabin M 1994 Absorbed fractions for electrons and beta particles in
spheres of various sizes J. Nucl. Med. 35 152–56
[18] Stabin M and Konijnenberg M W 2000 Re-evaluation of absorbed fractions for
photons and electrons in spheres of various sizes J. Nucl. Med. 41 149–60
[19] Stabin M 1996 MIRDOSE personal computer software for internal dose assessment
in nuclear medicine J. Nucl. Med. 37 538–46
[20] Cloutier R J, Watson E E, Rohrer R H and Smith E M 1973 Calculating the radiation
dosage to an organ J. Nucl. Med. 14 53–5
[21] Loevinger R, Budinger T and Watson E 1988 MIRD Primer for Absorbed Dose
Calculations S.o.N. Medicine, NM/MIRD (New York)
[22] ICRP 1979 Limits for Intakes of Radionuclides by Workers. ICRP Publication 30
I.C.o.R. Protection (New York: Pergamon)
[23] ICRP 1991 Recommendations of the International Commission on Radiological
Protection ICRP. Publication 60 I.C.o.R. Protection (New York: Pergamon)
[24] Stabin M and Sparks R B 1999 MIRDOSE4 does not exist J. Nucl. Med. Suppl. 40
309P
[25] Stabin M et al. 2001 RADAR—the radiation dose assessment resource. An online
source of dose information for nuclear medicine and occupational radiation safety
J. Nucl. Med. Suppl. 42 243P
[26] Nelson W R, Hirayama H and Rogers D W O 1985 The EGS4 Code System (Stanford
Linear Accelerator Center)
[27] Halbleib J A, Kensek R P, Valdez G D, Seltzer S M and Berger M J 1992 ITS: the
Integrated TIGER Series of electron/photon transport codes—version 3.0 IEEE
Trans. Nucl. Sci. 39 1025–30
[28] Hendricks J S and Briesmeister J F 1992 Recent MCNP developments IEEE Trans.
Nucl. Sci. 39 1035–40
[29] Brun R, Bruyant F, Maire M, McPherson A C and Zanarini P 1994 GEANT detector
description and simulation tool W5013 (CERN)
[30] Berger M J and Seltzer S M 1968 Electron and Photon Transport Programs. I. Intro-
duction and Notes on Program DATAPAC 4, II. Notes on Program ETRAN-15 Rep
NBS 9836 and 9837 (NBS)
References 155
[31] Cullen D E, Balinger C T and Perkins S T 1991 The all particle method: 1991 status
report UCRL-JC-108040 (Lawrence Livermore National Laboratory)
[32] Zaidi H 2000 Comparative evaluation of photon cross section libraries for materials
of interest in PET Monte Carlo simulations IEEE Trans. Nucl. Sci. 47 2722–35
[33] Sempau J, Acosta E, Baro J, Fernandez-Varea J M and Salvat F 1997 An algorithm
for Monte Carlo simulation of the coupled electron–photon transport Nucl. Instr.
Meth. B 132 377–90
[34] Fasso A, Ferrari A and Sala P R 2001 Proceedings of the MonteCarlo 2000 Conference
eds A Kling, F Barao, M Nakagawa, L Tavora and P Vaz (Berlin: Springer) pp 159–64
[35] Zanzonico P B, Brill A B and Becker D V 1995 Principles of Nuclear Medicine eds
H N Wagner, Z Szabo and J W Buchanan (Philadelphia, PA: Saunders) 2nd
edition, pp 106–34
[36] Benua R S, Cicale N R, Sonenberg M and Rawson R W 1962 The relation of
radioiodine dosimetry to results and complications in the treatment of metastatic
thyroid cancer Am. J. Roentgenol. 87 171–82
[37] Maxon H R et al. 1992 Radioiodine-131 therapy for well-differentiated thyroid
cancer—a quantitative radiation dosimetric approach outcome and validation in 85
patients J. Nucl. Med. 33 1132–36
[38] Samuel A M and Rajashekharrao B 1994 Radioiodine therapy for well-differentiated
thyroid cancer: a quantitative dosimetric evaluation for remnant thyroid ablation
after surgery J. Nucl. Med. 35 1944–50
[39] Schlesinger T, Flower M A and McCready V R 1989 Radiation dose assessments in
radioiodine (131I) therapy. 1. The necessity for in vivo quantitation and dosimetry in
the treatment of carcinoma of the thyroid Radiother. Oncol. 14 35–41
[40] Maxon H R 1999 Quantitative radioiodine therapy in the treatment of differentiated
thyroid cancer Q. J. Nucl. Med. 43 313–23
[41] DeNardo S J 2000 Tumor-targeted radionuclide therapy trial design driven by patient
dosimetry J. Nucl. Med. 41 104–6
[42] Snyder W, Ford M R, Warner G and Watson S 1975 ‘S’, absorbed dose per unit
cumulated activity for selected radionuclides and organs. MIRD Pamphlet No 11,
NM/MIRD (New York: Society of Nuclear Medicine)
[43] J A Siegel et al. 1999 MIRD Pamphlet No 16: Techniques for quantitative
radiopharmaceutical biodistribution data acquisition and analysis for use in human
radiation dose estimates J. Nucl. Med. 40 37S–61S
[44] M E Juweid et al. 2000 Initial experience with high-dose radioimmunotherapy of
metastatic medullary thyroid cancer using 131I-MN-14 FðabÞ2 anti-carcinoembryonic
antigen MAb and AHSCR J. Nucl. Med. 41 93–103
[45] Sgouros G 1993 Bone marrow dosimetry for radioimmunotherapy theoretical
considerations J. Nucl. Med. 34 689–94
[46] Wiseman G A, White C A, Witzig T A et al. 1999 Final dosimetry results of IDEC-
Y2B8 phase I/II 90yttrium radioimmunotherapy trial in non-Hodgkin’s lymphoma
J. Nucl. Med., Suppl. 40 64P
[47] Divgi C R et al. 1998 Phase I/II radioimmunotherapy trial with iodine-131-labeled
monoclonal antibody G250 in metastatic renal cell carcinoma Clin. Cancer Res. 4
2729–39
[48] Wahl R L, Zasadny K R, MacFarlane D et al. 1998 Iodine-131 anti-B1 antibody for
B-cell lymphoma an update on the Michigan Phase I experience J. Nucl. Med., Suppl.
39 21S–27S
156 Monte Carlo codes for use in therapeutic nuclear medicine
[68] De Matteis A and Pagnutti S 1995 Controlling correlations in parallel Monte Carlo
Parallel Comput. 21 73–84
[69] Mascagni M, Cuccaro S A, Pryor D V and Robinson M L 1995 A fast, high-quality,
and reproducible lagged-Fibonacci pseudorandom number generator J. Comput.
Phys. 15 211–19
[70] Martin W R and Brown F B 1987 Status of vectorized Monte Carlo code for particle
transport analysis Int. J. Supercomputer Appl. 1 11–32
[71] Smith M F, Floyd C E and Jaszczak R J 1993 A vectorized Monte Carlo code for
modeling photon transport in SPECT Med. Phys. 20 1121–7
[72] Kirkby D R and Delpy D T 1997 Parallel operation of Monte Carlo simulations on a
diverse network of computers Phys. Med. Biol. 42 1203–8
[73] Miura K 1987 EGS4V: vectorization of the Monte Carlo cascade shower simulation
code EGS4 Comput. Phys. Commun. 45 127–36
[74] Dewaraja Y, Ljungberg M, Majumdar A, Bose A and Koral K F 2002 A parallel
Monte Carlo code for planar and SPECT imaging: implementation, verification
and applications in I-131 SPECT Comput. Meth. Prog. Biomed. in press
Chapter 7
7.1. INTRODUCTION
158
Introduction 159
in water and other media. The definition of the variables used differentiated
the specific absorbed fraction (for sources of electrons or particles) from
‘energy-dissipation distribution’ as defined by Spencer [5] and, in particular,
allowed the introduction of ‘scaled absorbed-dose distribution’.
Thus, for a point source of electrons of energy E0 , the dose rate Rðx; E0 Þ
is given by
Rðx; E0 Þ ¼ AnkE0 ðx; E0 Þ: ð7:2Þ
This formulation (very close to the preceding one) defines ðx; E0 Þ as the ‘point
isotropic specific absorbed fraction’ expressed in g1 . Thus, the quantity
4x2 ðx; E0 Þ dx is the fraction of emitted energy absorbed in a spherical
layer of radius x and thickness dx in a medium of density (g/cm3 ). If it is
assumed that all the energy emitted is finally absorbed by the medium, then
ð1
4 x2 ðx; E0 Þ dx ¼ 1; ð7:3Þ
0
The methods used to generate DPKs have kept pace with changes in the
experimental data available and benefited from spectacular progress in the
field of numerical calculation. This section describes the different methods
used to generate point-kernels and their changes over time.
Until 1955, the only experimental data available concerned a dozen
radionuclides in air. It has been shown that these data could be represented
by the formula [1, 7]
k0 x ½1ðx=cÞ ð1xÞ
ðxÞ ¼ c 1 e þ x e x < c
ðxÞ2 c
ð7:8Þ
k0 ð1xÞ
ðxÞ ¼ e x c
x
Methods used to generate dose point-kernels 161
7.3.1. Photons
Monoenergetic. According to MIRD formalism [17], photon DPK in units
of (cGy/decay) at a distance r from a monoenergetic source of energy E
(MeV) in an unbounded homogeneous medium can be expressed as
en r
KðrÞ ¼ 1:6 108 E e Ben ðrÞ ð7:10Þ
4r2
where and en are respectively the linear attenuation and energy absorption
coefficients of the photon in the medium. The quantity Ben is the energy
absorption buildup factor, defined as the ratio between the absorbed dose
rate due to both primary and scattered photons and that due to primary
photons alone.
Using the moments method of Spencer and Fano, Berger [17] has
published Ben tables for photon energies from 15 keV to 3 MeV for distances
r from 0.05 to 20 mean free paths in water (1 mean free path ¼ 1/). At
individual photon energies within the energy range, Berger also expressed
the buildup factor as a ten-term polynomial of distance. Coefficients of the
polynomial have no obvious simple energy dependence. It has been shown
by Brownell et al. [18] that Berger’s method gives results comparable with
those obtained using a Monte Carlo code for photon transport in a large
soft-tissue volume.
Berger’s tabulated Ben values were fitted to analytic functions other
than polynomials by Kwok et al. [19] and Leichner [20]. Such fittings of
data minimize computer storage space. The parameters of the fitting function
have a simple relationship with photon energy for Kwok et al., while a
common expression applies to the absorbed fractions for photons and
particles for Leichner.
More recently, several groups of researchers have used Monte Carlo
codes to derive photon DPKs in water, namely ETRAN [13], EGS4 [21–
23] and MCNP [24, 25]. The agreement between the codes themselves and
with Berger’s calculation is within 3–5% [25–27].
Monte Carlo codes have been used directly to derive the photon DPKs of
many radionuclides including 123 I, 124 I, 125 I, 131 I, 111 In, 64 Cu, 67 Ga, 68 Ga,
99m
Tc, 131 Cs, 103 Pd and 71 Ge [21, 25, 28, 29]. Some of the DPKs have also
been parametrized to analytic functions for easy computation on personal
computers [21].
7.3.2. Electrons
For dosimetry of electrons, two major types of approach are possible,
depending on the problem involved. First, the problem to be resolved may
consist in comparing doses (or absorbed fractions) delivered to the same
target by different radionuclides. In this case, it seems preferable to calculate,
once and for all, a table of absorbed fractions for monoenergetic electrons
and then integrate these results into the emission spectra of the radionuclides
considered. Second, when various geometric configurations need to be
studied for the same radionuclide, it is often preferable to use the DPK of
the radionuclide considered.
In both cases, the application of data for monoenergetic electrons to all
electron emissions of a radionuclide should be performed with due care. First
of all, it is essential to consider the -emission spectra and not mean energy
[30]. Some spectra can be neglected if their emission percentage is too low to
contribute significantly to the delivered dose, provided that this choice is
indicated and clearly documented. Monoenergetic emissions (Auger or
conversion electrons) should be taken into account if they contribute signifi-
cantly to the delivered dose.
For example, the electron emission values given in table 7.1 for 131 I have
been published by the French Primary Laboratory of Ionizing Radiations
[31]. A comparison of these values with those cited by the MIRD Committee
[32] (in boldface in the table) indicates that the difference actually concerns
less than 1% of the energy emitted by transition. This section cites various
DPK sets used in the dosimetry of radionuclides.
0.5 keV to 10 MeV, which in practice allowed scanning of all of the necessary
energies in the domain of -emitting radionuclides. These data have been
widely used in dosimetric models or in DPK calculations for radionuclides.
Subsequently, it appeared that an incorrect sampling of energy-loss straggling
in ETRAN affected the validity of these data [13, 16]. The corrected values
now available extend from 20 keV to 20 MeV [33].
The tables for -ray dose distributions in water published by Cross et al.
[14, 34] were calculated using the ACCEPT code (derived from ETRAN).
For monoenergetic emissions, comparison with values obtained with the
EGS4 code [15] for six energies ranging from 50 keV to 3.0 MeV indicated
that differences between the two data sets were slight (a few per cent).
Fð; E0 Þ values were given for nine energies from 25 keV to 4 MeV.
an aqueous medium (123 I, 124 I, 125 I, 131 I, 111 In, 64 Cu, 67 Cu, 67 Ga, 68 Ga, 186 Re,
188
Re, 153 Sm, 117m Sn and 99m Tc). The results are expressed by providing
for each radionuclide the coefficients of a series of polynomials and
exponentials approximating the calculated values. Once again, the equations
proposed are easy to include in more complex calculation programs.
Bolch et al. [41] extended the DPK concept to cubic geometry by
presenting S value tables at the voxel level. The objective was to use these
S values in cases in which cumulated activity is determined from tomoscinti-
graphic images or quantitative autoradiographies, i.e., from numerical
images sampled in three dimensions (voxels). For each target voxel, the
mean dose is obtained by summing the doses delivered by each surrounding
voxel:
X
N
ðvoxelk Þ ¼
D Avoxelh Sðvoxelk voxelh Þ ð7:12Þ
h¼0
energy. Kwok et al. [42, 43] and Nunes et al. [44, 45] used both Monte Carlo
calculations and experiments to assess the effect of bone and air on dose
distribution of point sources of low-energy electrons and emitters near
tissue boundaries. When a point source of 147 Pm (maximum energy
225 keV) or 204 Tl (maximum energy 763 keV) was positioned at a planar
interface of cortical bone (CB) and red marrow (RM), the radiation dose
in RM within 6 mm from the interface was increased by 9 1% (SE) for
both sources, as compared with the dose in the same region inside a homo-
geneous volume of RM. The dose in a spherical shell of RM 20 mm thick
next to a spherical interface of CB/RM received from a monoenergetic
point source of electrons at the centre of a sphere of RM with a radius of
500 mm was increased by up to 22 1% (SE) compared with that for the
same scoring region in the homogeneous case. The maximum dose increase
occurred at 400 keV electron energy. The dependence of dose enhancement
on the curvature of the interface was also demonstrated. Conversely, when
a 32 P point source was positioned at a planar interface of air and soft
tissue, the radiation dose in soft tissue 5.5 mm from the boundary was reduced
by 26:44 0:02% (SE) from that in the homogeneous situation. There was
no simple modification of electron/ DPK due to tissue inhomogeneity.
An ingenious way to calculate the DPK in any homogeneous medium
from that in air (or water) was suggested by Cross and co-workers [14], and is
known as the ‘scaling factor’ method. Marcu and Prestwich [46] recently
extended the applicability of the scaling factor method to dissimilar media
with a planar interface. The investigation was done for a planar source of
32
P in water and the source was located at various distances from a planar
interface between water and any other homogeneous medium with atomic
number Z in the range 8 < Z < 50. Dose deposition discrepancies of less
than 5% were detected for the depth within which at least 95% of the emitted
energy was deposited. It will be of interest to develop procedures to apply the
scaling factor method to non-planar geometries and to more than one closely
spaced boundary.
It is rather difficult to classify the different approaches for the use of DPKs in
radionuclide dosimetry. The range of applications is rather broad (calcula-
tion of absorbed doses, absorbed fractions or S factors), the scales considered
extend from the cell to the organ, and the methods used are numerous
(numerical integration or convolution; hybrid methods associating Monte
Carlo simulation and the convolution of radioactive sources by DPKs).
The classification adopted here is based on the scale of the problem, from
cellular to organ dimensions.
A general approach for use of dose point-kernels in calculations 169
7.5. CONCLUSIONS
REFERENCES
[1] Loevinger R, Japha E M and Brownell G L 1956 Discrete Radioisotope Sources in Radia-
tion Dosimetry ed G J Hine and G L Brownell (New York: Academic) pp 693–755
[2] Berger M J 1970 Beta-ray dosimetry calculations with the use of point kernels. In
Medical Radionuclides: Radiation Dose and Effects eds R J Cloutier, C L Edwards
and W S Snyder (Washington, DC: US Atomic Energy Commission) pp 63–86
[3] Loevinger R and Berman M 1968 MIRD Pamphlet No 1: A schema for absorbed-
dose calculation for biologically distributed radionuclides J. Nucl. Med. 9 7–14
[4] Berger M J 1971 MIRD Pamphlet No 7: Distribution of absorbed doses around point
sources of electrons and beta particles in water and other media J. Nucl. Med. 12 5–23
[5] Spencer L V 1955 Theory of electron penetration Phys. Rev. 98 1597–615
[6] Spencer L V 1959 Energy Dissipation by Fast Electrons. Monograph No 1 (National
Bureau of Standards)
[7] Loevinger R and Holt J G 1956 Internally administered radioisotopes. In Radiation
dosimetry ed G J Hine and G L Brownell (New York: Academic) pp 801–73
172 Dose point-kernels for radionuclide dosimetry
8.1. INTRODUCTION
175
176 Radiobiology aspects and radionuclide selection criteria
8.2.2.1. Survival
Several in vitro assays have been described to measure the ability of cells to
proliferate. In practice, these assays measure the capacity of cells to success-
fully reproduce and, thus, to form a colony. Using a colony-forming assay, it
is possible to determine the decrease in survival, expressed as a surviving
fraction, as a function of a graded radiation dose.
Radiation survival curves are log–linear plots of surviving fraction (log)
versus dose (linear). The shape of the survival curve constructed through
such a set of survival points varies and will depend on certain biological,
physical, and chemical factors. In general, two types of dose–survival
curve have been described (figure 8.1). For the exponential survival curve,
the slope is always constant and can be expressed by
where S=S0 is the surviving fraction of the irradiated cells, D is the dose
delivered, and ¼ 1=D0 , the dose needed to reduce survival to 0.37 (for this
type of curve, D0 is also known as D37 ). Such exponential survival curves
are observed when mammalian cells are exposed to ‘densely ionizing’ radiation
(e.g., particle emissions, DNA-incorporated Auger electron emitters).
In the second type of dose–response relationship, expressed by the sig-
moidal survival curve (figure 8.1), the efficiency of cell kill is not constant: at
low doses, a slow decrease in survival is observed and the curve has a shoulder;
at higher doses, an exponential decrease in survival is seen. This type of survival
curve occurs routinely when mammalian cells are exposed to sparsely ionizing
radiation (e.g., X-rays, particles, extranuclear Auger electrons). Such curves
can be fitted using a second-order polynomial equation:
where a and b are the fit parameters and solving for D when S=S0 ¼ 0:37
yields the D37 value. Alternatively, and preferably, the curve can be fitted
Radiobiological effects 177
The delay in the progression of dividing cells through their cell cycle
following irradiation is a well known phenomenon. Usually it is reversible
and dose dependent, occurs only at specific points in the cell cycle, and is
similar for both surviving and nonsurviving cells. Cells undergoing mitosis
(M phase) continue through division basically undisturbed; those in the G1
phase of the cell cycle have very little delay, those in the S phase a moderate
delay, and those in the G2 phase maximum delay. The net result is that many
cells accumulate at the G2 =M boundary, and the ratio of mitotic cells to
nonmitotic cells (i.e., the mitotic index) is altered. The length of the delay
and the decrease in mitotic index are both functions of dose.
8.2.2.3. Redistribution
The radiosensitivity of cells is a function of their position within the cell cycle.
Cells in late S phase are most resistant (the survival curve following
irradiation has a wide shoulder), while those in the G2 and M phases are
most radiosensitive (the survival curves have a steep slope and no shoulder).
Consequently, following irradiation, the cells in the most sensitive phase
will be killed preferentially. This redistribution of cells will lead to partial
cell-cycle synchrony and a change in the overall radiosensitivity of the
cell cohort. One would expect a more radioresistant cell population, but in
reality the population rapidly becomes desynchronized and the net effect is
sensitization of the surviving population.
8.2.2.4. Repair
Mammalian cells are generally capable of repairing some of the damage
induced by radiation. This phenomenon is dose-rate dependent. As the
dose rate decreases, the capacity of cells to repair radiation damage increases
and this is manifested by a widening of the shoulder of the survival curve
accompanied by an increase in D0 (e.g., X-ray irradiation). In essence, two
types of repair have been described. The first, sublethal damage (SLD)
repair, occurs when a dose of radiation is split into two fractions and suffi-
cient time is allowed (0.5–2 h) for any/all of the radiation-induced damage
to be repaired. Naturally, if no repair is allowed to occur, for example, by
the immediate application of a second dose of radiation, the cells will die.
Sublethal damage and its repair have been shown to be important factors in
the sparing of normal tissues during fractionated radiotherapy. The second,
potentially lethal damage (PLD) repair, is observed only when mammalian
cells are grown under suboptimal conditions following a single dose of radia-
tion. Under such circumstances, an increase (rather than a decrease) in survival
is observed. This phenomenon is believed to be a result of the delayed entry of
the irradiated cells into mitosis, thereby allowing the cells to repair the PLD.
Both SLD and PLD as well as their repair have been reported only for X-ray
and irradiation (i.e., low-LET-type survival curves that have considerable
Radiobiological effects 179
shoulders) and are practically nonexistent for neutron and particle irradiation
(i.e., high-LET radiation).
Figure 8.2. Tumour growth inhibition by the bystander effect. During the 15
day period, the presence of 0:2 106 125 I-labelled cells leads to the deposition
of 2 cGy only in the growing tumours (unpublished results).
8.3.1.1. particles
Current radionuclide therapy in humans is based almost exclusively on -
particle-emitting isotopes. Typically, the electrons that are emitted from
the nuclei of the decaying atoms (1 electron/decay) have various energies
up to a maximum and thus have a distribution of ranges (table 8.2). As
each particle emitted traverses matter, it follows a contorted path, loses its
kinetic energy and eventually comes to a stop. Since the LET of these light,
charged particles is very low (0.2 keV/mm), except for the few nanometres at
the end of their range just before they stop, they are sparsely ionizing and
quite inefficient at damaging DNA and killing cells. Consequently, their
use as therapeutic agents predicates the presence of high radionuclide
concentrations within the targeted tissue and the traversal of several thou-
sand electrons per mammalian cell nucleus.
An important ramification of the long range (mm) of each emitted
electron is the production of cross-fire, a circumstance that negates the need
to target every cell within the tumour. For microscopic disease, however,
long-range and some medium-range emitters may deposit a significant fraction
of the energy of their particles outside the tumour [16]. This is particularly of
concern in the selection of a radiopharmaceutical for the palliation of bone
pain from metastatic osseous lesions [17]. Long-range emitters (e.g., 32 P) run
the risk of significantly irradiating bone marrow as well as bony lesions,
whereas short-range emitters (e.g., 33 P) are calculated to result in a significantly
lower dose to bone marrow relative to bone/bony lesions.
The matter of inhomogeneity in the distribution of a radionuclide and
its consequent dose has also been addressed. O’Donoghue [18] has derived
a method for calculating an equivalent uniform biologically effective dose
based on the absorbed-dose distribution represented by the biologically
effective dose–volume histogram. For larger tumours he postulates that a
combination of radionuclide and external beam therapy may result in
optimal dose distribution.
Many of the -particle-emitting radionuclides used for therapy also
release photons that generally do not add significantly to the dose delivered
to the target tissue. However, these photons may contribute considerably to
the whole-body dose. For example, 3.7 GBq of 131 I distributed throughout
the whole body would deposit approximately 60 cGy per day. Since the
bone marrow is usually the dose-limiting organ (200–300 cGy), the success
or failure of therapy will depend on not exceeding its MTD.
8.3.1.2. particles
An particle that is emitted by any of the radionuclides suitable/available
for therapy (table 8.3) is identical to a helium nucleus with energies ranging
from 4 to 9 MeV. The corresponding tissue ranges are 40 to 100 mm.
particles travel in straight lines, initially depositing their energy, in soft
tissue, at approximately 80 keV/mm and increasing their rate of energy
deposition towards the end of their tracks. The typical energy deposition
from an particle traversing the diameter of a cell nucleus is about
600 keV; this translates to an absorbed dose of about 0.14 Gy to the cell
nucleus per hit. These numbers depend, among other parameters, on the
size of the cell nucleus.
184 Radiobiology aspects and radionuclide selection criteria
8.3.2. Half-life
For a radiopharmaceutical with an infinite residence time in a tumour or
tissue, a radionuclide with a long physical half-life will deliver more decays
than one with a short half-life if both radiopharmaceuticals have the same
initial radioactivity. Moreover, there can be a striking difference in the
time-dependent dose rate delivered by the two. If the number of radionuclide
atoms per unit of tissue mass is n and the energy emitted (and absorbed) per
186 Radiobiology aspects and radionuclide selection criteria
agent’s being present within each of the targeted cells, i.e., it permits a certain
degree of heterogeneity. However, three factors will mainly determine
whether the dose delivered to the targeted tissue will be therapeutically
effective: (i) it is essential that the radiotherapeutic agent concentrate
within foci throughout the targeted tissue, (ii) the distances between these
hot foci must be equal to or less than twice the maximum range of the emitted
energetic particles and (iii) the concentration of the radiotherapeutic agent
within each hot focus must be sufficiently high to produce a cumulative cross-
fire dose to the surrounding targeted cells of approximately 10 000 cGy. Since
the dose is inversely proportional to the square of the distance, it is important
to note that the concentration of the therapeutic agent needed to deposit such
cytocidal doses will decrease precipitously when the distance between the hot
foci decreases.
Experimentally, these predictions have been substantiated in various
animal-tumour therapy studies. In many of these studies, radionuclides
have been targeted after conjugation to monoclonal antibodies. For example,
investigators have assessed the therapeutic efficacy of 131 I-labelled mono-
clonal antibodies in rodents bearing subcutaneous tumours. Although a
substantial proportion of cells within a tumour mass show reduced/no
expression of the targeted antigen and, therefore, are not targeted by the
radioiodinated antibody, 131 I-labelled antibodies that localize in high
concentrations in tumours are therapeutically efficacious and can lead to
total tumour eradication in some instances [49–52]. Thus, even when 131 I is
not-so-uniformly distributed within a tumour, the decay of this radionuclide
can lead to tumour sterilization as long as it is present in sufficiently high
concentrations (figure 8.3). Similar results have also been reported with
other -particle-emitting isotopes, in particular 90 Y [42, 53, 54] and 67 Cu
[55, 56].
Several ‘two- and three-step’ approaches have also been proposed to
target radioactivity to a tumour [57–61]. In general, an antibody that is
not internalized by the targeted cell is injected prior to the administration
of a small radiolabelled molecule that has a strong affinity to the antibody.
Figure 8.3. Tumour growth post-treatment with 131 I-labelled monoclonal anti-
body [51].
Experimental therapeutics 189
The most extensively studied approach utilizes the high avidity of avidin (Av)
or streptavidin (SAv) for biotin, a vitamin found in low concentration
in blood and tissues. Because the noncovalently bound Av/SAv–biotin
complex has an extremely low dissociation constant (kd ) of about 1015 M,
investigators have used these two molecules as binding pairs to bridge
molecules that have no affinities for each other and to target various radio-
nuclides. A recent report [61] that examined the efficacy of this approach
in tumour-bearing mice has demonstrated that 100% of the mice treated
with the 90 Y-labelled biotin derivative have been cured of their disease
(figure 8.4).
particle emitters are highly toxic [21, 62, 65–70], (ii) cells in monolayers
require more particle traversals than spherical cells [21, 62], (iii) only a
few (1–4) particle traversals through a mammalian cell nucleus are neces-
sary to sterilize a cell [21, 62, 70].
At the molecular level, the traversal of particles through a mammalian
cell nucleus leads to the efficient production of chromosomal aberrations and
DSB. For example, the incubation of mammalian cells with 211 At-astatide
causes a significant increase in chomosomal aberrations [62], but these
decline with the passage of time. More recently, Walicka et al. [24] have
established that more than ten DSBs are produced per decay of DNA-
incorporated 211 At, a value much higher than that obtained following the
decay of the DNA-incorporated Auger electron emitter 125 I [71, 72].
Investigators have also assessed the therapeutic potential of particle
emitters in tumour-bearing animals [31, 32, 68, 73–80]. Bloomer et al. [31,
32] have reported a dose-related prolongation in median survival when
mice bearing an intraperitoneal murine ovarian tumour are treated with
211
At-tellurium colloid administered directly into the peritoneal cavity.
While this -particle-emitting radiocolloid is curative without serious
morbidity, colloids of -particle-emitting radionuclides (32 P, 165 Dy, 90 Y)
are not (figure 8.6). In another set of in vivo studies examining the therapeutic
efficacy of a 212 Bi-labelled monoclonal antibody, the radionuclide is most
effective when used with a carrier having target specificity [75]. Finally, a
recent report by McDevitt and co-workers [80] has demonstrated that
225
Ac-labelled internalizing antibodies are therapeutically effective in mice
bearing solid prostate carcinoma or disseminated lymphoma. The targeting
ligand is radiolabelled with 225 Ac (T1=2 ¼ 10 days), a radionuclide that
decays to stable 209 Bi via a cascade of six daughters (221 Fr, 217 At, 213 Bi,
Experimental therapeutics 191
213
Po, 209 Tl, 209 Pb) with the emission of five and three particles. While the
proposed approach seems promising, the use of 225 Ac-labelled radio-
pharmaceuticals in humans is very likely to generate high normal-tissue
toxicity (as a consequence of the particle recoil-energy-induced dissociation
of the five chemically different radioactive atoms from the targeting carrier)
post decay of the 225 Ac-labelled molecules that have not been internalized by
the targeted tumour cells.
Figure 8.8. Survival fraction of V79 cells plotted as function of radiation dose
to cell nucleus. RBE of each agent is calculated at the D37 [90].
Experimental therapeutics 193
(figure 8.9). This chemo–radio combination in fact leads to 5–6-log kill and
cures approximately 30% of the animals.
REFERENCES
[1] Ward J F 1986 Mechanisms of DNA repair and their potential modification for
radiotherapy Int. J. Radiat. Oncol. Biol. Phys. 12 1027–32
[2] Goodhead D T 1985 Saturable repair models of radiation action in mammalian cells
Radiat. Res. 104 S-58–67
[3] Chapman J D 1980 Biophysical models of mammalian cell inactivation by radiation.
In Radiation Biology in Cancer Research eds R E Meyn and H R Withers (New York:
Raven) pp 21–32
[4] Thames H D and Hendry J H 1987 Fractionation in Radiotherapy (London: Taylor
and Francis)
[5] Williams M V, Denekamp J and Fowler J F 1985 A review of = ratios for
experimental tumors: implications for clinical studies of altered fractionation Int.
J. Radiat. Oncol. Biol. Phys. 11 87–96
[6] Nagasawa H and Little J B 1992 Induction of sister chromatid exchanges by extremely
low doses of -particles Cancer Res. 52 6394–6
[7] Azzam E I, de Toledo S M, Gooding T and Little J B 1998 Intercellular
communication is involved in the bystander regulation of gene expression in
human cells exposed to very low fluences of alpha particles Radiat. Res. 150 497–504
[8] Bishayee A, Rao D V and Howell R W 1999 Evidence for pronounced bystander
effects caused by nonuniform distributions of radioactivity using a novel three-
dimensional tissue culture model Radiat. Res. 152 88–97
[9] Nagasawa H and Little J B 1999 Unexpected sensitivity to the induction of
mutations by very low doses of alpha-particle radiation: evidence for a bystander
effect Radiat. Res. 152 552–7
References 195
[43] Ragni G and Szybalski W 1962 Molecular radiobiology of human cell lines II.
Effects of thymidine replacement by halogenated analogues on cell inactivation by
decay of incorporated radiophosphorus J. Mol. Biol. 4 338–46
[44] Burki H J, Roots R, Feinendegen L E and Bond V P 1973 Inactivation of
mammalian cells after disintegrations of 3H or 125I in cell DNA at 1968C Int. J.
Radiat. Biol. 24 363–75
[45] Liber H L, LeMotte P K and Little J B 1983 Toxicity and mutagenicity of X-rays and
[125I]dUrd or [3H]TdR incorporated in the DNA of human lymphoblast cells Mutat.
Res. 111 387–404
[46] Cleaver J E and Burki H J 1974 Biological damage from intranuclear carbon-14
decays: DNA single-strand breaks and repair in mammalian cells Int. J. Radiat.
Biol. 26 399–403
[47] Gutiérrez S, Carbonell E, Galofré P, Creus A and Marcos R 1998 The alkaline single-
cell gel electrophoresis (SCGE) assay applied to the analysis of radiation-induced
DNA damage in thyroid cancer patients treated with 131I Mutat. Res. 413 111–19
[48] Hengstler J G, Bockisch A, Fuchs J, Grimm W, Görges R, Oesch-Bartlomowicz B,
Zapf A-O, Lade K, Tanner B, Teichmann E, Thelen M, Gebhard S and Oesch F
2000 Induction of DNA single-strand breaks by 131I and 99mTc in human
mononuclear blood cells in vitro and extrapolation to the in vivo situation Radiat.
Res. 153 512–20
[49] Esteban J M, Schlom J, Mornex F and Colcher D 1987 Radioimmunotherapy of
athymic mice bearing human colon carcinomas with monoclonal antibody B72.3:
histological and autoradiographic study of effects on tumors and normal organs
Eur. J. Cancer Clin. Oncol. 23 643–55
[50] Riklund K E, Makiya R A, Sundström B E, Thornell L-E and Stigbrand T I 1990
Experimental radioimmunotherapy of HeLa tumours in nude mice with 131I-
labeled monoclonal antibodies Anticancer Res. 10 379–84
[51] Buchsbaum D, Khazaeli M B, Liu T, Bright S, Richardson K, Jones M and Meredith
R 1995 Fractionated radioimmunotherapy of human colon carcinoma xenografts
with 131I-labeled monoclonal antibody CC49 Cancer Res. 55 (suppl) 5881s–7s
[52] Sato N, Saga T, Sakahara H, Yao Z, Nakamoto Y, Zhang M, Kuroki M, Matsuoka
Y, Iida Y and Konishi J 1999 Intratumoral distribution of radiolabeled antibody
and radioimmunotherapy in experimental liver metastases model of nude mouse
J. Nucl. Med. 40 685–92
[53] Otte A, Mueller-Brand J, Dellas S, Nitzsche E U, Herrmann R and Maecke H R
1998 Yttrium-90-labelled somatostatin-analogue for cancer treatment Lancet 351
417–18
[54] Chinn P C, Leonard J E, Rosenberg J, Hanna N and Anderson D R 1999 Preclinical
evaluation of 90Y-labeled anti-CD20 monoclonal antibody for treatment of non-
Hodgkin’s lymphoma Int. J. Oncol. 15 1017–25
[55] DeNardo G L, Kukis D L, DeNardo S J, Shen S, Mausner L F, O’Donnell R T,
Lamborn K R, Meyers F J, Srivastava S C and Miers L A 1997 Enhancement of
67
Cu-2IT-BAT-Lym-1 therapy in mice with human Burkitt’s lymphoma (Raji)
using interleukin-2 Cancer 80 2576–82
[56] DeNardo G L, Kukis D L, Shen S, Mausner L F, Meares C F, Srivastava S C, Miers
L A and DeNardo S J 1997 Efficacy and toxicity of 67Cu-2IT-BAT-Lym-1
radioimmunoconjugate in mice implanted with human Burkitt’s lymphoma (Raji)
Clin. Cancer Res. 3 71–9
198 Radiobiology aspects and radionuclide selection criteria
[57] Hnatowich D J, Virzi F and Rusckowski M 1987 Investigations of avidin and biotin
for imaging applications J. Nucl. Med. 28 1294–302
[58] Goodwin D A, Meares C F, McCall M J, McTigue M and Chaovapong W 1988 Pre-
targeted immunoscintigraphy of murine tumors with indium-111-labeled bifunc-
tional haptens J. Nucl. Med. 29 226–34
[59] Le Doussal J-M, Gruaz-Guyon A, Martin M, Gautherot E, Delaage M and Barbet J
1990 Targeting of indium 111-labeled bivalent hapten to human melanoma mediated
by bispecific monoclonal antibody conjugates: imaging of tumors hosted in nude
mice Cancer Res. 50 3445–52
[60] Kassis A I, Jones P L, Matalka K Z and Adelstein S J 1996 Antibody-dependent
signal amplification in tumor xenografts following pretreatment with biotinylated
monoclonal antibody and avidin or streptavidin J. Nucl. Med. 37 343–52
[61] Axworthy D B, Reno J M, Hylarides M D, Mallett R W, Theodore L J, Gustavson L
M, Su F-M, Hobson L J, Beaumier P L and Fritzberg A R 2000 Cure of human
carcinoma xenografts by a single dose of pretargeted yttrium-90 with negligible
toxicity Proc. Natl Acad. Sci. USA 97 1802–7
[62] Kassis A I, Harris C R, Adelstein S J, Ruth T J, Lambrecht R and Wolf A P 1986
The in vitro radiobiology of astatine-211 decay Radiat. Res. 105 27–36
[63] Raju M R, Eisen Y, Carpenter S and Inkret W C 1991 Radiobiology of particles III.
Cell inactivation by -particle traversals of the cell nucleus Radiat. Res. 128 204–9
[64] Humm J L, Chin L M, Cobb L and Begent R 1990 Microdosimetry in
radioimmunotherapy Radiat. Prot. Dosim. 31 433–6
[65] Kurtzman S H, Russo A, Mitchell J B, DeGraff W, Sindelar W F, Brechbiel M W,
Gansow O A, Friedman A M, Hines J J, Gamson J and Atcher R W 1988 212Bismuth
linked to an antipancreatic carcinoma antibody: model for alpha-particle-emitter
radioimmunotherapy J. Natl Cancer Inst. 80 449–52
[66] Larsen R H, Bruland Ø S, Hoff P, Alstad J, Lindmo T and Rofstad E K 1994
Inactivation of human osteosarcoma cells in vitro by 211At-TP-3 monoclonal
antibody: comparison with astatine-211-labeled bovine serum albumin, free
astatine-211 and external-beam X rays Radiat. Res. 139 178–84
[67] Strickland D K, Vaidyanathan G and Zalutsky M R 1994 Cytotoxicity of -particle-
emitting m-[211At]astatobenzylguanidine on human neuroblastoma cells Cancer Res.
54 5414–9
[68] Larsen R H and Bruland Ø S 1998 Intratumour injection of immunoglobulins
labelled with the -particle emitter 211At: analyses of tumour retention,
microdistribution and growth delay Br. J. Cancer 77 1115–22
[69] Palm S, Bäck T, Claesson I, Delle U, Hultborn R, Jacobsson L, Köpf I and Lindegren S
1998 Effects of the alpha-particle emitter At-211 and low-dose-rate gamma-radiation
on the human cell line Colo-205 as studied with a growth assay Anticancer Res. 18
1671–6
[70] Nikula T K, McDevitt M R, Finn R D, Wu C, Kozak R W, Garmestani K, Brechbiel M
W, Curcio M J, Pippin C G, Tiffany-Jones L, Geerlings M W Sr, Apostolidis C, Molinet
R, Geerlings M W Jr, Gansow O A and Scheinberg D A 1999 Alpha-emitting bismuth
cyclohexylbenzyl DTPA constructs of recombinant humanized anti-CD33 antibodies:
pharmacokinetics, bioactivity, toxicity and chemistry J. Nucl. Med. 40 166–76
[71] Walicka M A, Adelstein S J and Kassis A I 1998 Indirect mechanisms contribute to
biological effects produced by decay of DNA-incorporated iodine-125 in
mammalian cells in vitro: double-strand breaks Radiat. Res. 149 134–41
References 199
9.1. INTRODUCTION
202
Introduction 203
the past few years, the remainder of this chapter will focus on a discussion of
each.
Over the past several years, researchers have taken an interest in the use of
particle emitters as a potential form of therapy for micrometastatic disease
[9–19]. The radiobiological advantages of particles are their high linear
energy transfer (LET) and independence from dose rate and oxygen effects
[20]. When combined with a suitable carrier (e.g., antibody), the short
range of particles in tissue (40–90 mm) results in a highly localized deposi-
tion of energy near the point of emission (presumably in the tumour) and
limits the volume of normal tissues irradiated.
The dosimetry of particle emitters is challenging due to the stochastic
nature of energy deposited in small, subcellular targets as shown schematically
in figure 9.1. In this diagram, some cell nuclei receive multiple particle hits,
while others receive no hits whatsoever. (In this context, we define a hit as any
particle that intersects a cell nucleus, regardless of the energy deposited.)
Additionally, an particle may traverse a variety of path lengths in passing
through a cell nucleus. For example, an particle that intersects a spherical
nucleus through the exact centre will have a path length equal to twice the
radius, whereas an particle that intersects nearly tangent to the cell nucleus
will have little or no path length. Since the energy deposited by an particle is
proportional to its path length, each particle in the above scenario produces
one hit; however, the amount of energy deposited is significantly different.
Since the distribution of hits to the cell nucleus is often characterized by
using Poisson statistics, this problem is often referred to as a compound
Poisson process—a Poisson distribution in the number of events where the
magnitude of each event varies. Other effects that influence the dosimetry of
particles include the variation in energy loss (LET) along an individual
particle’s path, the size and shape of the target, and the source configuration
(within the cell versus outside the cell, etc.). Thus, it is generally accepted
that the analysis of the cellular response to particle irradiation requires
microdosimetric methods. The importance of microdosimetry in therapeutic
applications (e.g., radioimmunotherapy) is that while two different micro-
dosimetric distributions may have the same average specific energy (absorbed
dose), they can result in significantly different levels of cell kill [21, 22].
ðx x0 Þ2 þ ðy y0 Þ2 þ ðz z0 Þ2 ¼ r2n ; ð9:4Þ
where x0 , y0 and z0 are the centre coordinates of a nucleus with radius rn .
Substituting equations (9.1)–(9.3) into equation (9.4) and solving for t
yields two solutions—t1 corresponding to the distance from the source to
where the particle enters the nucleus, and t2 corresponding to the distance
between the source and where the particle exits the nucleus. Note that when
the particle emission occurs within the nucleus (starter), t1 is chosen to be
zero, while t2 is the intercept with the nuclear surface. If the particle
terminates its track within the nucleus (stopper) the value of t2 is set to the
range of the particle. A hit occurs when the solutions for t1 and t2 are both
real, positive numbers.
Next, the energy deposited within the cell nucleus is calculated. This
energy is given by
ð R t1
dE
"1 ¼ dx: ð9:5Þ
R t2 dx
Thus, the determination of the specific energy deposited relies heavily on the
choice of stopping powers used. In a comparison of the impact of stopping
power formulation, Stinchcomb and Roeske [23] calculated specific energy
spectra for the same source–target geometry as Humm [22]. The former
case used Janni’s stopping power [33] for protons (appropriately scaled for
particles), while in the latter case Walsh’s stopping power [34] was utilized.
In general, a disparity of 1–7% was observed in the average specific energies
using these two different characterizations of the stopping power. However,
in practical situations, it is expected that the uncertainty in the activity
distribution will be much larger than that of the stopping powers used in
the Monte Carlo calculation.
The process of determining the specific energy deposited for an
individual interaction can be repeated to determine the multi-event specific
energy spectrum. In this case, all sources within the maximum range of the
particle emission are evaluated to determine whether they intersect the
nucleus and the subsequent specific energy deposited is calculated and
summed. This scenario represents one possible outcome. Next, the process
is repeated and, most likely, a different total specific energy will be calculated.
By repeating this process, the entire outcome of specific energies deposited
can be determined. A histogram of these individual results is the multi-
event distribution. An alternative is to determine the single-event specific
energy spectrum. This distribution is calculated by examining only those
particle emissions that result in exactly one hit to the cell nucleus. Each
time an individual particle intersects the nucleus, the specific energy is
determined and recorded. All of these single-event specific energies are
then histogrammed to determine the single-event spectrum. Single-event
spectra are often used in conjunction with analytical codes whereby the
multi-event spectra are calculated by multiple convolutions of the single-
event spectrum. Stinchcomb and Roeske [23] demonstrated that both
methods produce the same multi-event spectra within the uncertainty of
the calculation itself. An example of a single-event spectrum for 212 Bi
distributed on the cell surface is shown in figure 9.2. A multi-event spectrum
for 15 212 Bi sources on the cell surface is shown in figure 9.3.
Microdosimetric spectra are often combined with a model of cell
survival to estimate the survival following a particular type of irradiation.
The surviving fraction following particle irradiation is given by [22, 23, 35]
ð1
SðzÞ ¼ f ðzÞ ez=z0 dz ð9:7Þ
0
Figure 9.2. Single-event specific energy distribution for 212 Bi sources distribu-
ted on a cell surface. The cell radius is 10 mm and the nuclear radius is 5 mm. The
two peaks correspond to the two emissions from 212 Bi: 6.05 MeV (36%) and
8.78 MeV (64%).
where hni is the average number of hits to the cell nucleus, and T1 ðz0 Þ is the
Laplace transform of the single-event spectrum. Note that the same value of
z0 satisfies both equations (9.7) and (9.8).
It is important to point out that z0 is not equal to D0 , which is
determined from the slope of the cell survival curve. Rather it is a more
fundamental quantity as D0 has folded into it not only the effects of the
radiation, but also the effects of the source–target geometry. In the
applications to be discussed a recurrent theme is to use the observed
surviving fraction to determine the inherent cell sensitivity, z0 .
divided into three regions corresponding to an outer, middle and inner layer
of cells. Cell survival was assessed individually for each layer. Additionally, a
Monte Carlo code was developed to calculate specific energy spectra as a
function of depth within the spheroid assuming the activity was concentrated
on the outer 10 mm. Their results demonstrated a wide variation in the
specific energy distribution in the three different levels of the spheroid. The
average dose was a maximum near the outer spheroid surface and decreased
rapidly with depth. Conversely, the fraction of cells receiving zero hits was a
minimum near the surface and increased rapidly with depth. Additionally,
they noted a qualitative correlation that cells near the outer edges of
the spheroids had a lower surviving fraction compared with those inside.
Charlton [45] produced a calculational model for multicellular spheroids
and simulated their subsequent survival. The distribution of particle
tracks throughout cells within the spheroid was calculated and combined
with a cell survival model that takes into account the survival as a function
of LET [46]. The results of this analysis indicated that cell survival decreased
from 57 to 37% as the spheroid diameter increased from 75 to 225 mm for a
uniform source with one decay per cell, and 50% packing. Longer ranged
particle emitters increased the number of hits per cell for the larger spheroids.
Increasing the packing fraction from 40 to 70% decreased the survival from
46 to 26% in 200 mm diameter spheroids, again with one decay per cell. This
significant decrease in survival was due to the increased crossfire dose com-
ponent as cells were brought closer together. In a different simulation, a
small fraction of cells (20%) was assumed not to take up any activity due
to lack of expression of a surface antigen. In this simulation, the unlabelled
fraction did not significantly change the expected cell survival. In both of
these studies, there is evidence of a highly nonuniform specific energy
distribution that varies with distance from the spheroid surface. Thus, a
single specific energy distribution is not representative of that through the
entire tumour. Combining specific energy distributions with a model of cell
survival may provide an overall measure of the therapeutic effectiveness.
However, these cell survival models do not take into account second-order
processes such as bystander effects. These processes are more difficult to
simulate and may play a significant role for these types of geometry. The
ongoing refinement of these cell survival models is currently an area of
active research [47].
Bone marrow is often the dose-limiting organ in RIT. The dosimetry of
bone marrow is difficult due to the complex geometry as well as the presence
of tissue inhomogeneities. Thus, idealized models, as have been utilized in
previous studies, must be replaced by more realistic geometries. The work
to date on estimating specific energy spectra for bone marrow has focused
largely on using histological samples obtained from humans or animal
models. While there has been considerable interest in bone marrow micro-
dosimetry for protection purposes [48–51] this review will focus only on
Auger electron emitters 213
Auger electrons are emitted by any radionuclide that decays by either electron
capture or internal conversion. Both these decay processes result in an inner
214 Microdosimetry of targeted radionuclides
Out of this work and that of many other contributors have emerged the
following inferences:
9.3.1. Monte Carlo simulation of Auger decays and their energy deposition
There is considerable stochastic variability in the way that energy is released
from electron capture or internal conversion of decaying radionuclides. This
variability was first modelled by Monte Carlo simulation of the decay process
of 125 I [60]. Since then, Auger electron spectra have been calculated for
several other radionuclides [62]. The number and distribution of energy
between the electrons (and photons) emitted in an Auger cascade varies
between individual disintegrations as governed by the inner shell transition
probabilities. The number of permutations of possible de-excitation path-
ways through the atom can be hundreds, depending on the atomic number
of the isotope. The way in which electron capture (or internal conversion)
and the ensuing Auger cascade are computer simulated is described in
detail by Charlton and Booz [60], Humm [68] and Howell [62]. In brief, a
random number generator is used to select an initial inner shell vacancy
according to the probabilities of capture within each of the respective orbitals
K, L1, L2, L3 etc. The probabilities for all transitions to this vacancy are
normalized to one, where each allowed transition is weighted by its respective
probability according to published values for radiative [69] and non-radiative
[70] transitions, determined from quantum mechanical calculations of
oscillator strengths. A second random number is then used to select one of
these transitions, according to their respective weights. The new electron
occupancy of the orbitals is recorded and the next innermost shell vacancy
is selected. This process then repeats until all of the vacancies have moved
to the outermost shells (the lowest energy state), when further transitions
are no longer possible.
For high Z elements (e.g., tellurium) following K-capture in iodine, the
fluorescent yield is high (0.877). This means that 87.7% of the transitions to
fill a K-shell vacancy emit characteristic X-rays and only 12.3% result in a K-
shell Auger electron. This results in a variable total amount of electron
216 Microdosimetry of targeted radionuclides
energy released, contingent upon the photon transitions which occur. One
example of a possible atomic de-excitation pathway following electron
capture for 123 I spectra is schematically depicted in figure 9.4.
Simulation of 10 000 individual atomic de-excitations results in an
output of individual Auger electron and X-ray spectra, containing the
stochastics of the de-excitation process. These data can be averaged to
yield an average Auger electron emission spectrum, an example of which is
shown in figure 9.5. Note that this does not denote an emission spectrum
from an individual decay, but rather a frequency weighted average of
10 000 simulated decays. Dosimetric estimates to different compartments
of the cell have been performed using the average Auger electron spectrum,
and the results are the same as when multiple individual Auger electron
spectra are used. The stochastics of individual spectra is only pertinent
when calculating the energy deposition along the DNA.
Charlton and Humm [71] performed Monte Carlo electron track
simulations for 10 000 individual Auger electron spectra resulting from 125 I
and other radionuclides [66]. They scored the coordinates of each energy
deposition generated by the electron track structure code of Paretzke [72],
relative to the origin of the site of decay placed at the position of a base in
a computer-simulated geometric model of the DNA double helix. This
model consisted of a central cylinder representing the bases surrounded by
Auger electron emitters 217
123
Figure 9.5. An example of an average electron energy spectrum from I
following 10 000 simulated decays.
LET and the consequent radiotoxicity effects are unwarranted. However, the
extreme short range effects of Auger electron emitters do pose a new problem
in assessing the dosimetry of novel radiopharmaceuticals. Humm et al. [74]
and Goddu et al. [75] have recommended some guidelines for their evalua-
tion, which requires the initial determination of the partition of the agent
into extracellular, cell surface bound, intracellular, intranuclear and
appended to the DNA. Cellular S-factors for the former four compartments
have been calculated by Goddu et al. [76, 77] for a number of electron capture
and internal conversion decaying radionuclides. The average radiation dose
to the nucleus (or other compartment) can be determined by the piecewise
addition of each contribution, upon experimental determination of the rela-
tive compartmental partitions.
the activity distribution, cell geometry and LET response may be combined to
predict cell survival. Since cell response in vitro may not correlate with tumour
response in vivo, such an approach will require feedback to refine such models.
It is the hope that such an approach would help avoid useless therapies in
patients who may not benefit (such as those with highly resistant cells), or
aid in determining the optimal dosing within an individual patient. Such appli-
cations are potentially years away from implementation, but may provide a
method of optimizing radionuclide therapy to improve the therapeutic gain.
REFERENCES
[13] Zalutsky M R, Stabin M, Larsen R H and Bigner D D 1997 Tissue distribution and
radiation dosimetry of astatine-211-labeled chimeric 81C6, an alpha-particle-
emitting immunoconjugate Nucl. Med. Biol. 24 255–1
[14] Rotmensch J, Whitlock J L, Schwartz J L, Hines J J, Reba R C and Harper P V 1997
In vitro and in vivo studies on the development of the alpha-emitting radionuclide
bismuth 212 for intraperitoneal use against microscopic ovarian carcinoma Am. J.
Obstet. Gyne. 176 833–41
[15] Kennel S J and Mirzadeh S 1997 Vascular targeting for radioimmunotherapy with
213
Bi Radiochimica Acta 79 87–91
[16] Larsen R H, Vaidyanathan G and Zalutsky M R 1997 Cytotoxicity of alpha-particle-
emitting 5-[211At]astato-20 -deoxyuridine in human cancer cell Int. J. Radiat. Biol. 72(1)
79–90
[17] Larsen R H, Akabani G, Welsh P and Zalutsky M R 1998 The cytotoxicity and
microdosimetry of astatine-211-labeled chimeric monoclonal antibodies in human
glioma and melanoma cells in vitro Rad. Res. 149 155–62
[18] Couturier O, Faivre-Chauvet A, Filippovich I V, Thedrez P, Sai-Maurel C, Barides
M, Mishra A K, Gauvrit M, Blain G, Apostolodis C, Molinet R, Abbe J-C, Bataille
R, Wijdenes J, Chatal J-F and Cherel M 1999 Validation of 213Bi-alpha
radioimmunotherapy for multiple myeloma Clin. Cancer Res. 5 3165s–70s
[19] Sgouros G, Ballangrud A M, Jurcic J G, McDevitt M R, Humm J L, Erdi Y E,
Mehta B M, Finn R D, Larson S M and Scheinberg D A 1999 Pharmacokinetics
and dosimetry of an alpha-particle emitter labeled antibody: 213Bi-HuM195 (Anti-
CD33) in patients with leukemia J. Nucl. Med. 40(11) 1935–46
[20] Hall E J 2000 Radiobiology for the Radiologist (Philadelphia: Lippincott, Williams
and Wilkins)
[21] Humm J L 1987 A microdosimetric model of astatine-211 labeled antibodies for
radioimmunotherapy Int. J. Radiat. Oncol. Biol. Phys. 13 1767–73
[22] Humm J L and Chin L M 1993 A model of cell inactivation by alpha-particle internal
emitters Rad. Res. 134 143–50
[23] Stinchcomb T G and Roeske J C 1992 Analytic microdosimetry for radioimmuno-
therapeutic alpha emitters Med. Phys. 19(6) 1385–93
[24] Polig E 1978 The localized dosimetry of internally deposited alpha-emitters Curr.
Top. Rad. Res. 13 189–327
[25] Olko P and Booz J 1990 Energy deposition by protons and alpha particles in spherical
sites of nanometer to micrometer diameter Radiat. Environ. Biophys. 28 1–17
[26] Eckerman K F, Ryman J C, Taner A C and Kerr G D Traversal of cells by radiation
and absorbed dose fraction estimates for electrons and alpha particles. In
Proceedings of the Fourth International Radiopharmaceutical Dosimetry Symposium
eds A T Schlafke-Stelson and E E Watson (CONF-851113) (Oak Ridge, TN: Oak
Ridge Associated Universities) pp 67–86
[27] Attix F H 1986 Introduction to Radiological Physics and Radiation Dosimetry (New
York: Wiley)
[28] Ziegler J F, Biersack J P and Littmark U 1985 The Stopping and Range of Ions in
Matter (New York: Pergamon)
[29] Ziegler J F 1988 The stopping and range of ions in solids. Ion Implantation: Science
and Technology (San Diego: Academic) pp 3–61
[30] Ziegler J F 1999 The stopping of energetic light ions in matter Appl. Phys. Rev., J.
Appl. Phys. 85 1249–72
References 223
[31] Bichsel H 1992 Stopping powers of fast charged particles in heavy elements NIST
Report IR-4550 (Springfield, VA: US National Technical Information Service)
[32] ICRU 1993 Stopping powers and ranges for protons and alpha particles Report No 49
(Bethesda, MD: International Commission of Radiation Units)
[33] Janni J F 1966 Calculation of energy loss, range, path length, straggling, multiple
scattering and the probability of inelastic nuclear collisions for 0.1 to 1000 MeV
protons Technical Report No AFWL-TR-65-150
[34] Walsh P J 1970 Stopping power and range of alpha particles Health Phys. 19 312–16
[35] Charlton D E and Sephton R 1991 A relationship between microdosimetric spectra
and cell survival for high-LET irradiation Int. J. Radiat. Biol. 59(2) 447–57
[36] Fisher D R 1985 The microdosimetry of monoclonal antibodies labeled with alpha
emitters. In Proceedings of the Fourth International Radiopharmaceutical Dosimetry
Symposium eds A T Schlafke-Stelson and E E Watson (CONF-851113) (Oak Ridge,
TN: Oak Ridge Associated Universities) pp 26–36
[37] Roesch W C 1977 Microdosimetry of internal sources Radiat. Res. 70 494–510
[38] Stinchcomb T G and Roeske J C 1995 Survival of alpha particle irradiated cells as a
function of the shape and size of the sensitive target (nucleus) Radiat. Prot. Dosim.
62(3) 157–64
[39] Kvinnsland Y, Stokke T and Aurilien E 2001 Radioimmunotherapy with alpha-
particle emitters: Microdosimetry of cells with heterogeneous antigen expression
and with various diameters of cells and nuclei Radiat. Res. 155 288–96
[40] Stinchcomb T G and Roeske J C 1994 Analysis of survival of C-18 cells after irradia-
tion in suspension with chelated and ionic bismuth-212 using microdosimetry
Radiat. Res. 140 48–54
[41] Charlton D E 1999 The survival of monolayers of cells growing in clusters irradiated
by 211At appended to the cell surfaces Radiat. Res. 151 750–3
[42] Aurlien E, Larsen R H, Akabani G, Olsen D R, Zalutsky M R and Bruland O S 2000
Exposure of human osteosarcoma and bone marrow cells to tumor-targeted alpha-
particles and gamma-irradiation: Analysis of cell survival and microdosimetry Int. J.
Radiat. Biol. 76(8) 1129–41
[43] Charlton D E, Kassis A I and Adelstein S J 1994 A comparison of experimental and
calculated survival curves for V79 cells grown as monolayers or in suspension
exposed to alpha irradiation from 212Bi distributed in growth medium Radiat.
Prot. Dosim. 52 311–15
[44] Kennel S J, Stabin M, Roeske J C, Foote L J, Lankford P K, Terzaghi-Howe M,
Patterson H, Barkenbus J, Popp D M, Boll R and Mirzadeh S 1999 Radiotoxicity
of bismuth-231 bound to membranes of monolayer and spheroid cultures of
tumor cells Radiat. Res. 151 244–56
[45] Charlton D E 2000 Radiation effects in spheroids of cells exposed to alpha emitters
Int. J. Radiat. Biol. 76(11) 1555–64
[46] Charlton D E and Turner M S 1996 Use of chord lengths through the nucleus to
simulate the survival of mammalian cells exposed to high-LET alpha-radiation
Int. J. Radiat. Biol. 69(2) 213–17
[47] Iyer M and Lehnert B E 2002 Alpha-particle-induced increases in the radioresistance
of normal human bystander cells Radiat. Res. 157 3–7
[48] Humphreys E R and Humm J L 1988 A Monte-Carlo approach to the micro-
dosimetry of 224Ra in murine compact and cancellous bone Health Phys. 54(6)
607–15
224 Microdosimetry of targeted radionuclides
228
Dosimetry software design considerations 229
not. Dose estimates that differed from the ‘true’ value by factors of 2 to 10 still
resulted in numbers that fell within the diagnostic range. This allowed for the
creation of simple dose-per-unit administered activity tables such as are found
on all package inserts for modern day commercial radiopharmaceuticals.
The advent of hybridoma technology raised the possibility of attempting
to treat cancer with antibodies targeted to specific tumour antigens.
Antibodies offer a carrier medium that is potentially selective at the
molecular level. Hybridoma technology created a vehicle whereby the
antibodies could be manufactured in abundance. The possibility of tagging
these ‘magic bullets’ with radioactivity to treat cancers became suddenly a
reality. The desire to use radiolabelled antibodies in a therapeutic modality,
however, placed an entirely new demand on the dosimetry community.
Factors of 2 and 10 were no longer acceptable. Delivering an estimated dose
of 200 mCi to the bone marrow would allow a patient to survive. If this esti-
mate were off by a factor of 2, the difference could result in death.
The initial clinical trials that infused radiolabelled antibodies into
patients demonstrated that the concept of a ‘magic bullet’ was a mistaken
one. Radioactivity did not localize in the tumour exclusively: it concentrated
in other organs that were clearly visible on whole body scans. If the carrier
antibody was classified as a ‘good’ antibody, it localized preferentially in
the tumour, but NOT exclusively. Consequently, investigators were forced
to consider the following question: Given that radioactivity localizes in tissues
other than tumour, how much radioactivity can the patient tolerate? This
question had to be weighed in conjunction with its corollary: How much
radioactivity must be administered to sterilize the tumour? Having to consider
both questions simultaneously, investigators wondered if it was possible to
obtain the latter before the former. Answers required accurate knowledge
of the kinetics of a particular antibody, the physical decay parameters
associated with the radiolabel, and a thorough knowledge of the geometric
composition and anatomic makeup of an individual patient. This knowledge
did not exist in a form that was usable on a practical scale.
To begin addressing the problem of tumour dose versus normal organ dose in
the context of radiolabelled antibodies, the sources of error during the
process of forming an internal radionuclide dose estimate were analysed.
Additionally, the pace of computer technology development was evaluated.
These areas of inquiry resulted in a list of design features that an optimal
dosimetry solution should consider. The solution that was arrived at—a
practical software tool for estimating monoclonal antibody doses—was
developed in direct response to the design considerations. These design
considerations are summarized below.
230 The MABDOSE program for internal radionuclide dosimetry
Fourth, the total body remains as the single most important space in
which intravenously administered activity concentrates. This implies that
quantification of the total body should be performed. If one chooses to use
a SPECT acquisition, this becomes time prohibitive. Typically, when
SPECT has been used, investigators have decided to select for quantification
only organs that fall within the field of view of a single acquisition rotation. If
one is serious about quantifying activity at the voxel level, one should
quantify each and every voxel where activity occurs. Investigators, however,
evidently find this too great an expenditure of time since no one has ever
performed it in a clinical trial. The primary reason for not quantifying
total body activity from a series of total body SPECT acquisitions (aside
from the time commitment) is the gross errors introduced into the quantifi-
cation process, upwards of 400%!
Fifth, the most radiosensitive organ in the body is the bone marrow.
The current experience with radiolabelled antibodies suggests that binding
affinities are not so great that one gets wholesale uptake in the tumour
and none anywhere else. With intravenous infusions of radiolabelled anti-
body and the assumption that higher amounts of activity will result in
improved response rates, the bone marrow continues to be the critical
organ. This organ (more properly an organ space) is beneath the resolution
limits of modern scintillation cameras. For the foreseeable future, it must be
dealt with by means of a mathematical model that is predictive of biological
outcome. This model should be capable of incorporating measures of bone
marrow size (reserve) as well as measures of bone marrow viability
(patency) given the reality that radioimmunotherapy patients have usually
undergone several rounds of chemotherapy. This treatment is only
amenable to a mathematical modelling analysis. This statement derives
from the fact that the anatomy of interest is not visualizable on a practical
basis.
Sixth, the fact that cancers are the target for infused radiation means
that one does not have the luxury of being able to use an ‘average’ organ
located at an average position within the mathematical representation of
an average person. Each patient that presents will have a unique number
of cancer foci. Each focus will have a unique position and shape within the
patient’s body. A completely general dosimetry program should have the
ability to allow for the creation and construction of these arbitrary volumes
on a patient-by-patient basis, and incorporate them into the dose calculation
scheme.
Seventh, in 1985 the speed of microprocessors was increasing regularly,
while the surface area occupied by integrated circuits was shrinking. The
likelihood for CRAY supercomputer speeds on a desktop computer system
in the near future was considered high. Calculation algorithms that
previously had been considered time prohibitive might be judged feasible if
given a second look.
232 The MABDOSE program for internal radionuclide dosimetry
10.3.1. Version 1
The development of the MABDOSE program has evolved in a series of
stages. The first version of the program was presented at the Annual Meeting
of the American Association of Physicists in Medicine in 1986 [5]. To set the
historic record straight, the original version was called MABDOS because
file names were restricted to being six letters under the DOS operating
system (The MIRDOS program that Oak Ridge produced one year later
was subjected to similar constraints). With the relaxation in subsequent
years of filename size limits to allow longer and more descriptive filenames,
the current appellation is MABDOSE. The name is an acronym that
stands for Monoclonal Antibody Dosimetry. The acronym is a bit
misleading, because it is indicative of the program’s origins rather than its
functionality. MABDOSE is a completely general internal radionuclide
dosimetry program. It was initially conceived to address the shortcomings
inherent in the standard approach when dealing with dosimetry problems
posed by monoclonal antibodies as carrier molecules for radionuclides. It
is not restricted, however, to only performing dose estimates for monoclonal
antibodies.
The first version of the MABDOSE program had defining characteristics
that included (i) selection of an arbitrary radionuclide from the spectral library
assembled by Kocher [16], (ii) creation of a graphical interface for defining
mathematical models that predates the appearance of Apple’s Macintosh
GUI, (iii) curve fitting routines as alternatives to mathematical modelling,
(iv) incorporation of numerical integration routines to obtain cumulated
activities, (v) implementation of the ALGAM67 Monte Carlo code for
photon transport, (vi) incorporation of local deposition of non-penetrating
radiation for modelling electron and particulate radiation transport and (vii)
display of mean organ dose in a list format [7]. The program used the
234 The MABDOSE program for internal radionuclide dosimetry
equations for the Adult Reference Man phantom [17] to create a ‘logical’
target geometry, as opposed to a ‘virtual’ target geometry. The difference is
important: the former merely exists as a set of equations that are checked
for determining whether an x=y=z coordinate triplet falls within the volume’s
boundary. Thus, if the left kidney were declared a source organ, and a starting
coordinate desired, the equations for the left kidney would be consulted. A
random number generator would be used to obtain an arbitrary x=y=z co-
ordinate triplet from within the rectangular volume that just superscribed
the organ’s surface boundary. The triplet would then be checked against the
equations governing the left kidney to determine if the triplet satisfied the
defining equations.
The Monte Carlo simulation engine of the MABDOSE code was
derived from ALGAM67, a software code originally obtained from the
Oak Ridge Radiation Shielding Information Center in 1984. Historically,
ALGAM67 was the code that was used to create the MIRD11 tables. It
was deemed prudent at the time to use the same simulation code, because
ultimately the MABDOSE results would be compared back to the
MIRD11 tables as the nominal standard for internal radionuclide dosimetry
calculations.1 The simulation was structured so that each photon greater
than 10 keV was simulated for the same, user-determined, number of
histories. A random starting coordinate was selected from within a declared
source organ, and the photon tracked until it had either deposited all of its
energy or had escaped from the target phantom. Average dose to organs
identified in the Reference Man phantom were reported in a list.
The MABDOSE program was written in Fortran. The initial phase of a
dose estimate (source organ declaration, time–activity data entry, curve
fitting, curve integration, radionuclide selection) and the final phase
(combination of absorbed fractions with cumulated activities, display of
dose by organ) were completed on a Texas Instruments Professional
microcomputer. The intermediate portion of the program (Monte Carlo
simulation of radioactive decay using the declared source organs and
declared radionuclide spectrum) was performed on a Cray supercomputer.
This intermediate step necessitated the uploading of the solution’s inter-
mediate ‘state’, then downloading the results for final dose output.
10.3.2. Version 2
The second version of the program implemented a number of refinements to
Version 1. The speed of microcomputers was increasing exponentially. It was
1
In the interim, three other codes have gained prominence in simulation circles: EGS4, BCMP
and ETRAN. Their popularity does not detract from the applicability of the ALGAM67
algorithm, nor the practical extensions to it that includes non-penetrating radiation from a
nuclide’s spectrum.
Historic development of MABDOSE 235
readily seen that it was only a matter of time before microcomputers would
have calculation speeds approaching those of a CRAY. The simulation
portion of the program was ported to a microcomputer platform, resulting
in a tightly coupled program intended to execute completely from within a
single desktop hardware platform. At the same time, the entire MABDOSE
program was rewritten in the Cþþ language. This rewrite was occasioned by
the large size the original program had grown to. Changes to the code for
refinements, enhancements, and added functionality were impeded by the
use of a non-object-oriented language.
10.3.3. Version 3
The third version of the program opted for the creation of a full-fledged
virtual target geometry dataset. A three-dimensional lattice was physically
created in computer memory, with each voxel in the lattice representing
a single tissue type of the target. Each voxel in the lattice contained a
single index representative of its tissue type. The collection of lattice voxels
created a target geometry that served as the space in which the Monte
Carlo simulation was carried out. The same space served as the means for
declaring sources of activity. The source data structure encapsulated a
column vector whose elements corresponded to all voxels containing the
index representing that source. Selecting a starting coordinate was reduced
to calling a pseudo-random number generator to provide a value distributed
between 0 and 1, then multiplying by the length of the source’s column
vector.
Once created, a three-dimensional target lattice could be stored in a
database library for future recall. The rules governing the six different
Reference Man phantoms specified by Cristy and Eckerman [17] were imple-
mented in a hidden layer. If a target of the selected age and resolution existed
in the library, it was ‘checked out’ for use in the current simulation. If the
target did not exist, it was created by first allocating the lattice in memory.
The equations governing the appropriate age phantom were then consulted
to determine which organ index should be assigned each voxel.
10.3.4. Version 4
The fourth version of the program decoupled the target lattice from the
Reference Man phantom. This was accomplished by prefacing the lattice
data set with a header. The header contains a simple list of each voxel type
that will be encountered when reading in the lattice. Each voxel type is
represented by three fields: the first field is the integer index assigned to
represent that tissue type; the second is an ASCII descriptor of the tissue
type the index represents; and the third is an integer index that represents
the medium type, governing the voxel’s interaction characteristics with
236 The MABDOSE program for internal radionuclide dosimetry
radiation, and currently restricted to selection from among bone, soft tissue
and lung tissue. In a manner analogous to that of decoupling the target lattice
from the Reference Man, it is possible to generalize the assignment of media
types so that they can be incorporated after dynamic declaration (see section
10.5).
10.4.1. Overview
The majority of an internal radionuclide dose estimate deals with manipulating
data into a form that is usable by other program components. Too much
emphasis has been placed on the fact that MABDOSE uses a Monte Carlo
algorithm for photon transport based on the ALGAM67 code, and that it
does not use one of the more modern codes such as EGS4 or ECMP. The
reality is that MABDOSE accounts for all radiation: the algorithm used for
transporting penetrating (i.e., photon) radiation is based on ALGAM67; the
algorithm used for transporting non-penetrating radiation (i.e., electrons,
particles, photons <10 keV, and particles) is local absorption. It is logically
inconsistent to simulate electron transport when the resolution from nuclear
medicine scans routinely supersedes the range of particles. Instead, R&D
effort on MABDOSE has been devoted to generalizing the user interface for
ease of use, and to expanding the range of problems for which the interface
is applicable. This means that doses at boundaries between source and
target will necessarily have an error associated with them, since particles
will in reality ‘spill over’ the boundary edge and lose some of their energy in
the adjacent tissue, rather than stopping abruptly at the interface boundary.
This error is mitigated somewhat by adjacent whole body activity.
MABDOSE requires four objects be declared during a dosimetry session.
These are (i) a target lattice, (ii) a set of source volumes, (iii) cumulated
activities for each source volume in the set and (iv) a radionuclide from
which the dose is desired. Each of these four objects is discussed below in
detail. The order for making these objects known to the program is general,
with two exceptions. First, a source cannot be declared until a target lattice
has been created. Second, cumulated activities for a source cannot be
generated until the source has been defined. A declared source is not
considered defined until it has at least one time–activity point assigned to it.
Figure 10.1. The MABDOSE main window, showing descriptors for the
possible voxel ‘types’ in the bottom left most list. The source ‘Kidneys’ has
been constructed in the bottom middle list from the voxel types, ‘L Kidney’
and ‘R Kidney’. When this list is closed, the source ‘Kidneys’ will be added to
the patient source list displayed in the bottom right. Each source must have a
separate time–activity data set associated with it. Cumulated activity is
apportioned equally among voxels that comprise the source.
The target list in the left-most pane corresponds to the possibilities for
declaring and defining a source. A source is a logical construct that has
a physical space in the target lattice associated with it. The underlying
theme of a source is homogeneity of activity concentration: no matter
what sub-volume of a given source is sampled, when compared with another
sub-volume from the same source, the activity concentration will be the
same. This is realizable in a statistical sense. Thus, one could expect both
the left and right kidneys to handle a given radiopharmaceutical in exactly
the same way. One could then construct a single source, called ‘Kidneys’
from the two volumes that comprise the left and right kidney spaces.
The centre pane is essentially a temporary workspace that allows one to
construct a source from the left-most list of organelle possibilities. When a
list of all organelles that are to be treated as a single source has been
constructed, the list in the centre pane is closed, and the source added to
the patient’s current source list. The right-most pane on the lower half of
the program’s main window represents the current source list.
MABDOSE: from start to finish 239
For each declared source, an estimate must be made of the amount of activity
in that source as a function of time. This is frequently, but not exclusively,
derived from nuclear medicine scintigrams. The ability to incorporate higher
accuracy measurements from blood samples, urine collections and tissue biop-
sies motivated the development of the current MABDOSE interface. Each
source has associated with it a list of time–activity data point pairs. The inter-
face object for identifying time–activity point pairs to the MABDOSE system
is a two-column spreadsheet. Right clicking on a declared source displays the
spreadsheet associated with that source. Individual time–activity data are
input according to the abundance obtained for each source.
The time–activity data are fitted to a mathematical function, and
integrated. The integral of the time–activity curve gives the total number
of nuclear disintegrations, a quantity known as the cumulated activity.
Combined with knowledge of the radionuclide’s decay spectrum, the type
and total number of each energy in the spectrum is known. The cumulated
activity is divided by the number of voxels that compose the source being
simulated. Again, the rationale for doing this is the concept of source homo-
geneity: what defines a source is uniformity of function over the set of
discrete voxels that compose the source.
The mathematical functions that are integrated to achieve the cumu-
lated activity are arrived at using two main approaches. The first is curve
fitting. Curve fitting presupposes that the operator knows how the functional
form of the data should look, although not the function’s magnitude. There-
fore, the operator must have an appropriate function beforehand. The
functions that MABDOSE currently implements are a mono-exponential
function, a bi-exponential function, and a linear function. Additionally,
MABDOSE can perform a piecewise integration using discrete trapezoidal
or Riemann steps between data points.
The second technique for obtaining cumulated activities is mathematical
modelling. MABDOSE implements a graphical interface for composing a
mathematical model (figure 10.2). Each source has associated with it a
colour-coded rectangle. These are positioned on the screen in a format that
is convenient to the user. Directed pathways are then established by starting
a line segment within an activity source, and terminating the line segment in
an activity sink. The interrelationships specified by the created model define a
series of differential equations whose solutions are sums of exponentials.
An initial estimate of the exchange rate between linked compartments is
made. Starting from these estimates, a derivative-free version of the
Levenberg–Marquardt algorithm minimizes the sum of squares difference
between the raw data points and the differential equations by altering the
rate constants in a stepwise fashion [18]. The amplitude of each exponential
term corresponds to an initial value problem, and can be calculated given the
240 The MABDOSE program for internal radionuclide dosimetry
Figure 10.3. Dose display, accompanied by a DVH for the selected liver
volume.
of organs can be highlighted and the DVH for this composite organ list
displayed.
There are a number of enhancements that are envisioned for further develop-
ment. Each of these enhancements provides flexibility to the overall program
that increases the generality of the problem specification.
manual contouring is the only method that is currently error free, and this is
time prohibitive for implementation on a routine patient basis. Robust
segmentation algorithms that require minimal operator interaction are
necessary to surmount this hurdle. The algorithm currently being
investigated is a thin-plate spline routine [24–26]. The routine would operate
on a spiral CT data set that brackets a patient volume from head to foot.
Output from the routine would be a three-dimensional data set with all
organs surfaces segmented according to the RTOG format [27].
10.6. SUMMARY
REFERENCES
[1] Loevinger R and Berman M 1968 A schema for absorbed-dose calculations for
biologically distributed radionuclides J. Nucl. Med. 1 (suppl) 7
[2] Loevinger R and Berman M 1976 A revised schema for calculating the absorbed dose
from biologically distributed radionuclides, MIRD Pamphlet No 1, Revised. With a
foreword by R J Cloutier and E M Smith (New York: Society of Nuclear Medicine)
[3] Snyder W S, Ford M R, Warner G G and Watson S B 1975 ‘S’: absorbed dose per unit
cumulated activity for selected radionuclides and organs, MIRD Pamphlet No 11 (New
York: Society of Nuclear Medicine)
[4] Stabin M 1996 MIRDOSE: personal computer software for internal dose assessment
in nuclear medicine J. Nucl. Med. 37 538–46
References 247
[5] Johnson T K 1986 A tumor preferential dosimetry schema for generalized use in
radionuclide monoclonal antibody therapy Med. Phys. 13 586 (abst)
[6] Johnson T K and Vessella R L 1987 A generalized dosimetry schema for tumor
preferential uptake of monoclonal antibodies in radionuclide immunotherapy J.
Nucl. Med. 28(4) (suppl) 680 (abst)
[7] Johnson T K 1988 MABDOS: a generalized program for internal radionuclide
dosimetry Comput. Meth. Programs Biomed. 27 159–67
[8] Johnson T K and Vessella R L 1989 On the possibility of ‘real-time’ Monte Carlo
calculations for the estimation of absorbed dose in radioimmunotherapy Comput.
Meth. Programs Biomed. 29 205–10
[9] Johnson T K and Vessella R L 1989 On the application of parallel processing to the
computation of dose arising from the internal deposition of radionuclides Comput.
Phys. 3 69–72
[10] Johnson T K, McClure D L, McCourt S L, Andl G J, Berman B D and Newman F D
MABDOSE: a computer program for the calculation of dose. In Proceedings of the
6th International Radiopharmaceutical Dosimetry Sysmposium, Gatlinburg, TN, 7–
10 May 1996 (ORISE 99-0164, January 1999) pp 425–39
[11] Johnson T K, McClure D, McCourt S and Newman F 1997 MABDOSE:
Characterization and validation of a general purpose dosimetry code Med. Phys.
24 1025 (abst).
[12] Johnson T K 1997 Implementation of a 3-D internal dosimetry system on a PC Med.
Phys. 24 990–1 (abst)
[13] Johnson T K, McClure D L and McCourt S L 1999 MABDOSE. I. Characterization
of a general purpose dosimetry code Med. Phys. 26(7) 1389–95
[14] Johnson T K, McClure D L and McCourt S L 1999 MABDOSE II: Validation of a
general purpose dosimetry code Med. Phys. 26(7) 1396–403
[15] Zubal I G, Harrell C R, Smith E O, Rattner Z, Gindi G and Hoffer P B 1994
Computerized 3-dimensional segmented human anatomy Med. Phys. 21 299–302
[16] Kocher D C DLC-80/DRALIST data package (RSIC Data Library DLC-80, Oak
Ridge National Laboratory)
[17] Cristy M and Eckerman K F 1987 Specific absorbed fractions of energy at various ages
from internal photon sources. I. Methods (ORNL/TM-8381/V1 April)
[18] Press W H, Teukolsky S A, Vetterling W T and Flannery B P 1992 Numerical Recipes
in C: the Art of Scientific Computing 2nd edition (Cambridge: Cambridge University
Press) pp 656–706
[19] Roussin R W, Knight J R, Hubbell J H and Howerton R J 1983 Description of the
DLC-99/HUGO package of photon interaction data in ENDF/B-V format ORNL/
RSIC-46 (Oak Ridge National Laboratory)
[20] Hubbell J H 1969 Photon cross sections, attenuation coefficients and energy absorption
coefficients from 10 keV to 100 GeV NSRDS-NBS 29
[21] Seltzer S M 1993 Calculation of photon mass energy-transfer and mass energy-
absorption coefficients Radiat. Res. 136 147–70
[22] Snyder W S, Ford M R and Warner G G 1975 MIRD Pamphlet No 5, Revised:
Estimates of specific absorbed fractions for photon sources uniformly distributed in vari-
ous organs of a heterogeneous phantom (New York: Society of Nuclear Medicine)
[23] Warner G G and Craig A M Jr 1986 ALGAM: a computer program for estimating
internal dose from gamma-ray sources in a man phantom ORNL/RSIC CCC-152
(Oak Ridge National Laboratory) reviewed November 1986
248 The MABDOSE program for internal radionuclide dosimetry
11.1. INTRODUCTION
The standard formalism for patient dosimetry was developed by the Medical
Internal Radiation Dose (MIRD) Committee and is described in reference
[1], and also reviewed in chapter 4 of this book. As indicated above, to
simplify absorbed dose calculations, the MIRD Committee developed ‘S
factor’ tables and an associated procedure for their use [2]. This procedure
has been implemented in a software package, MIRDOSE3 [3]. To generate
tables of S factors for different radionuclides and source–target organ combi-
nations, a standard model of human anatomy was adopted in which organ
position, dimensions, and composition were mathematically defined. The
radioactivity was assumed to be uniformly distributed throughout each
source organ and the S factors were defined as the mean absorbed dose to
a target organ per unit cumulated activity in a source organ. Since the posi-
tion and size of tumours may vary within the body and since a standard
model of human anatomy was adopted for generating the S factor tables,
tumours are not included in the published tables. A number of approaches
have been developed for estimating the absorbed dose to tumours and the
dose contribution from tumours to normal organs. The simplest approxima-
tion is made by assuming that all electrons are deposited locally and that the
relative contribution to the tumour-absorbed dose from photons is negligi-
ble. Alternatively, the fraction of electron energy absorbed may be consid-
ered assuming the tumour can be modelled as a sphere [4, 5]. Using tables
of photon-absorbed fraction to spheres or ellipsoids, the photon self-dose
may be added by assuming that the tumour is a sphere or ellipsoid [6]. If
this assumption is made, the photon dose to and from normal organs may
also be calculated by placing the idealized tumour geometry in a defined posi-
tion relative to the standard geometry used for the S factor calculations [7, 8].
249
250 The three-dimensional internal dosimetry software package, 3D-ID
Monte Carlo-based dosimetry techniques can account for density () and
atomic number (Z) variations that are not easily accounted for using a
point-kernel approach. A Monte Carlo dosimetry algorithm was developed
for 3D-ID implementation that accounts for individual patient details by
using a SPECT or PET radionuclide activity distribution to initialize a
simulation of photons across a CT geometry. For a given radionuclide, initial
photons are generated according to the voxelized activity distribution
provided by SPECT or PET. Photon energies for the given radionuclide
are sampled from a published decay scheme and transport is performed
using the EGS4 transport code [38]. Efficient particle transport within the
discrete geometry was achieved by treating contiguous voxels having similar
Z as a single medium, and correcting for density variations by adjusting the
path length of the particle in proportion to the density, as determined from
CT. Discrimination between bone and water (which was used to represent
tissue) was carried out using a water–bone density threshold (WBT) with
WBT < 1:4 g cmÿ3 considered water or air and WBT 1:4 g cmÿ3 to be
bone. An additional gain in efficiency was obtained by scoring the absorbed
dose along photon, rather than electron, tracks. The method was verified by
creating discrete density and activity arrays for the Standard Phantom
geometry and comparing specific absorbed fraction for various organ and
energy combinations.
It is important to note that Monte Carlo-based dosimetry in which
activity, density and medium distribution are obtained directly from imaging
can lead to artefacts due to problems with the input information. For
example, images of Monte Carlo-derived dose distributions will depict
higher absorbed doses in the lungs than in tumour lesions in the lungs, even
252 The three-dimensional internal dosimetry software package, 3D-ID
Figure 11.2. Multi-image display from 3D-ID showing four transverse slices
from a set of CT-SPECT aligned images. The top row shows the CT with
ROIs drawn over three lesions identified as: right external iliac in blue, left
external iliac in green, and right common iliac in light blue. The bottom row
shows the corresponding 131 I SPECT slices with the lesions indicated. (See
plate 3 for colour version.)
the operating system for the dose calculation (currently UNIX or PC) and
the resulting batch file will be appropriately defined. An estimated number
of calculations as well as sample timings are provided as a guideline for
the user.
The dose calculation yields a three-dimensional array of dose values,
which are essentially images of dose distribution and can be viewed as
‘maps’ of the spatially varying dose in a target volume from a particular
source volume or volumes. The modular characteristic of 3D-ID provides
the flexibility to replace the current dose calculation technique with other
techniques that may be either point-kernel or Monte Carlo based.
Once the calculation has been performed, the final step in 3D-ID is the
display and analysis of the three-dimensional dose data set. This includes:
calculation of the mean absorbed dose to a target per unit cumulated activity
from a source (or sources); generation of dose–volume histograms of source
to target regions in order to summarize the dose information in a clinically
interpretable manner; and presentation of all image data in a summary
display of anatomic, functional and dose distribution images.
If absorbed dose instead of dose-rate is desired and time-sequential
SPECT or PET studies are unavailable, the absorbed dose in each target
256 The three-dimensional internal dosimetry software package, 3D-ID
Figure 11.5. The spatial absorbed dose distribution is depicted so that image
intensity corresponds to absorbed dose values. Isodose contour levels for 10%
(green), 25% (blue), 50% (red) and 75% (yellow) of the maximum dose
to the tumour mass are shown. The image was obtained from 3D-ID analysis
of 124 I PET studies in patients with thyroid carcinoma [36]. (See plate 4 for
colour version.)
Considerations in clinical implementations 257
Registration accuracy will depend upon the anatomical ROI, the agent being
imaged and the technique that is used for registration. Since it is likely
that the total number of counts associated with any given voxel will be low,
yielding a high standard deviation, the voxel-by-voxel-based integration
outlined above may be fundamentally limited by the error associated with
each voxel value. Integration on a voxel-by-voxel basis may yield unreliable
values for residence time or cumulated activity due to these errors. A
number of techniques and approximations may be adopted to overcome
some of these difficulties. It is possible to derive a residence time or cumulated
activity image, for example, by using planar imaging to obtain kinetics and
SPECT or PET imaging to obtain the spatial distribution. A residence time
or cumulated activity image may be obtained by assuming that all points
making up the spatial distribution of activity follow the same kinetics [33].
ACKNOWLEDGMENTS
This work was supported by NIH grants R01 CA62444, U01 CA58260 and
P01 CA33049 and DOE grant DE-FG02-86ER-60407.
REFERENCES
[1] Loevinger R, Budinger T F and Watson E E 1991 MIRD Primer for Absorbed Dose
Calculations revised edition (New York: Society of Nuclear Medicine)
References 259
[2] Snyder W S, Ford M R, Warner G G and Watson S B 1975 MIRD Pamphlet No 11:
‘S’ absorbed dose per unit cumulated activity for selected radionuclides and organs
(New York: Society of Nuclear Medicine)
[3] Stabin M 1996 MIRDOSE: personal computer software for internal dose assessment
in nuclear medicine J. Nucl. Med. 37 538–46
[4] Stabin M and Konijnenberg M W 2000 Re-evaluation of absorbed fractions for
photons and electrons in spheres of various sizes J. Nucl. Med. 41 149–60
[5] Siegel J A and Stabin M 1994 Absorbed fractions for electrons and beta particles in
spheres of various sizes J. Nucl. Med. 35 152–6
[6] Ellett W H and Humes R M 1971 MIRD Pamphlet No 8: Absorbed fractions for
small volumes containing photon-emitting radioactivity J. Nucl. Med. 12(5)
(suppl) 25–32
[7] Johnson T K, McClure D and McCourt S 1999 MABDOSE I: Characterization of a
general purpose dose estimation code Med. Phys. 26 1389–95
[8] Johnson T K and Colby S B 1993 Photon contribution to tumor dose from
considerations of 131 I radiolabeled antibody uptake in liver, spleen and whole
body Med. Phys. 20 1667–74
[9] Kolbert K S, Sgouros G, Scott A M et al. 1997 Implementation and evaluation of
patient-specific three-dimensional internal dosimetry J. Nucl. Med. 38 301–8
[10] Furhang E E, Sgouros G and Chui C S 1996 Radionuclide photon dose kernels for
internal emitter dosimetry Med. Phys. 23 759–64
[11] Sgouros G, Barest G, Thekkumthala J et al. 1990 Treatment planning for internal
radionuclide therapy: three-dimensional dosimetry for nonuniformly distributed
radionuclides J. Nucl. Med. 31 1884–91
[12] Liu A, Williams L E, Wong J Y and Raubitschek A A 1998 Monte Carlo-assisted
voxel source kernel method (MAVSK) for internal beta dosimetry Nucl. Med.
Biol. 25 423–33
[13] Flux G D, Webb S, Ott R J, Chittenden S J and Thomas R 1997 Three-dimensional
dosimetry for intralesional radionuclide therapy using mathematical modeling and
multimodality imaging J. Nucl. Med. 38 1059–66
[14] Giap H B, Macey D J and Podoloff D A 1995 Development of a SPECT-based three-
dimensional treatment planning system for radioimmunotherapy J. Nucl. Med. 36
1885–94
[15] Giap H B, Macey D J, Bayouth J E and Boyer A L 1995 Validation of a dose-point
kernel convolution technique for internal dosimetry Phys. Med. Biol. 40 365–81
[16] Akabani G, Hawkins W G, Eckblade M B and Leichner P K 1997 Patient-specific
dosimetry using quantitative SPECT imaging and three-dimensional discrete
Fourier transform convolution J. Nucl. Med. 38 308–14
[17] van Dieren E B, Plaizier M A, van Lingen A, Roos J C, Barendsen G W and Teule G J
1996 Absorbed dose distribution of the Auger emitters 67Ga and 125I and the beta-
emitters 67Cu, 90Y, 131I and 186Re as a function of tumor size, uptake and intracellular
distribution Int. J. Radiat. Oncol. Biol. Phys. 36 197–204
[18] van Dieren E B, van Lingen A, Roos J C, Huijgens P C, Barendsen G W and Teule G J
1994 A dosimetric model for intrathecal treatment with 131I and 67Ga Int. J. Radiat.
Oncol. Biol. Phys. 30 447–54
[19] Erdi A K, Yorke E D, Loew M H, Erdi Y E, Sarfaraz M and Wessels B W 1998 Use
of the fast Hartley transform for three-dimensional dose calculation in radionuclide
therapy Med. Phys. 25 2226–33
260 The three-dimensional internal dosimetry software package, 3D-ID
[38] Nelson W R, Hirayama H and Rogers D W O 1985 The EGS4 code-system SLAC
Report No 265 (Stanford, CA: SLAC)
[39] Erdi Y E, Mawlawi O, Larson S M et al. 1997 Segmentation of lung lesion volume by
adaptive positron emission tomography image thresholding Cancer 80 2505–9
[40] Kolbert K S, Watson T, Matei C, Xu S, Koutcher J A and Sgouros G 2002 Murine-
specific S-factors J. Nucl. Med. 43 88P
Chapter 12
12.1. INTRODUCTION
262
Monte Carlo simulation of imaging systems 263
Since SPECT and planar imaging provide the underlying information for
activity quantification and biokinetics, evaluation and validation of the dosi-
metry should also consider reconstruction methods, registration methods of
CT and SPECT/PET data, scatter and attenuation compensation methods
264 Dose calculation procedures in patient-specific radionuclide therapy
Figure 12.2. The left image shows the rendered image of the original Zubal
phantom. The middle image shows the extensions of arms and legs that Stuchly
et al made using the VH data set. The right image shows the straightening of the
arms that makes the phantom useful for patient-like whole-body simulations in
RNT.
spatial resolution [11, 12]. By fitting NURBS to patient data, the phantom
generated is more realistic than methods based on solid geometry definitions.
The NCAT phantom also includes temporal functions, such as modelling of
the heart motion and the respiratory movements. All organ and skeletal
models (except for the heart model) are based on CT information of the
VH male data set. The heart originates from a gated MRI scan of a human.
131
12.4.1. Evaluation of planar I imaging
A fundamental imaging parameter to study is the spatial resolution where the
point-spread function (PSF) describes the imaging system’s response to an
infinitely small point-source (a function). Figure 12.5 shows a comparison
between measured and simulated energy spectra and PSFs for a 131 I point
source in air for a high-energy collimator (HEGAP) and a ultra-high-
energy collimator (UHEGAP). The reason for using two collimators was
to investigate how well the UHEGAP collimator, with its thicker lead
septa, reduces septum penetration. The good agreement between measured
and simulated energy spectra was achieved only after inclusion of backscatter
photons in the simulation model [13, 14]. Validation of the energy spectrum
was also carried out for a case where there was significant object scatter using
a phantom consisting of a hot sphere centred in a cylindrical water tank [14].
The wide tails of the 131 I PSF (see figure 12.5) are due to septum penetration,
which degrades focal quantification accuracy because of the spread of
regional counts to surrounding areas. Correcting for this is especially
important in 131 I RNT where antibody uptake by surrounding organs is
significant compared with tumour uptake.
Monte Carlo is ideal for assessing penetration and scatter since, unlike
in experimental measurement, these events can be tracked separately in a
simulation. For 131 I, a quantitative comparison of scatter and penetration
has been performed for the HEGAP and the UHEGAP collimators [13].
The geometric, scatter and penetration components of 364 keV photo-peak
268 Dose calculation procedures in patient-specific radionuclide therapy
Figure 12.4. Simulated whole-body images of 131 I with the activity distribution
corresponding to a 131 I labelled monoclonal antibody distribution and four
simulated tumours. Image A represents the true distribution, obtained by analy-
tical integration. Image B represents a measured image, obtained by Monte
Carlo simulation and detector parameters representing a typical scintillation
camera system. The intensity distribution is altered, i.e., the contrast resolution
is heterogeneous, due to the variable attenuation and scattering properties of the
body.
The suppression of collimator scatter and penetration with the use of the
UHEGAP collimator is also evident in the comparison of 131 I energy spectra
and PSFs for the two collimators (figure 12.5).
131
12.4.2. Evaluation of SPECT I imaging
Validation of 131 I SPECT simulations has been carried out for a clinically
realistic situation using a computerized version of the experimental Radio-
logical Support Devices (RSD) thorax phantom. The physical phantom
was altered by inserting a water-filled plastic sphere, representing a
tumour, close to the liver. The organ and tumour relative activity concentra-
tion ratios were realistic values based on 131 I RIT patient images. Because of
the demanding computational requirements, the simulation was carried out
using the parallel implementation of SIMIND on an IBM-SP2 distributed
memory architecture [15]. There was good agreement between the measured
and simulated images, as demonstrated in figure 12.6.
Tumour activity quantification, such as the one simulated in figure 12.6,
is affected by the spread or blurring of regional counts to surrounding areas
due to the finite spatial resolution of the system. The ‘spill-out’ of counts
270
Figure 12.5. Comparison of the measured and simulated 131 I energy spectra and PSFs for Picker’s collimators HEGAP and UHEGAP.
(Reprinted by permission of the Society of Nuclear Medicine from [13].)
Dose calculation procedures in patient-specific radionuclide therapy
Evaluation of image quality degradation by Monte Carlo methods 271
from the target to the background decreases the tumour activity quantifica-
tion, as the source structure gets smaller while the ‘spill-in’ contribution from
the background to the target increases the tumour activity quantification.
The spill-out depends on the shape and size of the source structure and the
system resolution while spill-in depends on these parameters as well as the
activity level and distribution in the background. An experimental evaluation
of the effect of different tumour sizes and shapes on activity quantification
can be technically difficult to achieve, while Monte Carlo simulation is a
powerful tool to carry out such a study. The effects of object shape, size
and background on SPECT 131 I activity quantification without detector
response compensation has been reported in a recent Monte Carlo study
[16]. The results of the study are summarized in table 12.2. Here b is the
background to tumour activity concentration ratio. Note that the activity
error is 0% for one of the simulations because the calibration factor was
also determined from this simulation. The activity quantification was carried
out using a constant calibration factor and no partial volume correction. The
physical object boundary was used to define the volume of interest (VOI),
which is consistent with patient quantification procedures where CT-defined
tumour boundaries are superimposed on to the registered SPECT image.
Results of table 12.2 show that the activity error increases significantly
with decreasing tumour size because of the increased spill-out of counts.
The activity error for cylindrical tumours is consistently higher than the
error for the spherical tumours because spill-out is more significant for
non-spherical objects. Also, activity error depends significantly on b because
of the effect of spill-in.
Table 12.2. The activity error presented as the percentage difference between
the true 131 I activity and the SPECT estimated activity for different sizes,
shapes and b values.
The Monte Carlo method has been shown to be very efficient in both
calculating the absorbed dose and also as a tool for evaluating various
parts of the calculation schemes [19]. The key point is that if a simulation
of SPECT images (including all inherent limitations such as photon attenua-
tion, scatter contribution, collimator and septum penetration effects) can
accurately be made, then the output from a proposed calculation procedure
can be compared directly with the corresponding dose calculation using the
known activity and attenuation distribution on which the imaging simula-
tions were based.
Figure 12.7 shows the basic steps in such an evaluation. An advantage
of Monte Carlo methodology is that simulations can be made of realistic
distributions of, in principle, any radiopharmaceutical distribution and
that complex density distributions can be included.
Figure 12.8. Example of how the Monte Carlo method can help evaluate the
accuracy in a three-dimensional absorbed dose calculation program. The
figure shows a transversal section through the lungs (top row) and a section
through the kidneys (bottom row). From the definition of the true activity
distribution (A) ideal images without scatter and attenuation effects can be
simulated (B) and used as reference with which the particular scatter and
attenuation method should be compared (C). Note (a) the higher absorbed
dose shown as a border around the lungs which originated from in-spill counts
from the nearby located tissue and (b) the hot-spot that is a result of an air-
cavity that result in a very local high absorbed dose due to the low density of
the air.
(Scol =) and mass energy absorption coefficient (en =) [66, 67]. The system,
however, has low sensitivity and can only be used for higher absorbed doses,
as expected in therapeutic nuclear medicine. Another system that could be
suggested is the gel dosimetry system read with ordinary MR cameras.
Here one uses gels containing ferrous ions. These dosimeters are soft tissue
equivalent, however, also requiring relative high absorbed doses [68]. They
can be made in arbitrary shapes.
ACKNOWLEDGMENTS
This work was partially funded by US Public Health Service grant RO1
CA80927, the Swedish Cancer Foundation, The Gunnar Nilsson Founda-
tions and the Bertha Kamprads Foundations.
REFERENCES
[29] Strand S-E and Strandh M 1994 Parameters to consider for measurements of
absorbed doses in vivo with mini-thermoluminescent dosimeters Cancer 73(3) 985–8
[30] Strandh M and Strand S-E 1996 In vivo absorbed dose measurements with mini-
TLDs. Parameters affecting the reliability Acta Oncol. 35(6)
[31] Martin S, Lisbona A, Richard J, Morteau S, Denizot B and Bardies M 2000
Production of new thermoluminescent mini-dosimeters Phys. Med. Biol. 45(2) 479–94
[32] Aissi A and Poston J W 1984 An improved volumetric dosimeter for internal dose
verification Health Phys. 46(2) 371–6
[33] Skalka M 1971 Biochemical Indicators of Radiation Injury in Man (Vienna:
International Atomic Energy Agency)
[34] Bender M A and Gooch PC 1962 Types and rates of X-ray-induced chromosome
aberrations in human blood irradiated, in vitro Proc. Natl Acad. Sci. USA 48
52232
[35] Chambers K A, Harrington N P, Ross W M and Filion L G 1998 Relative alterations
in blood mononuclear cell populations reflect radiation injury in mice Cytometry
31(1) 45–52
[36] Pantelias G E and Maillie H D 1984 The use of peripheral blood mononuclear cell
prematurely condensed chromosomes for biological dosimetry Radiat. Res. 99(1)
140–50
[37] Lenarczyk M and Slowikowska M G 1995 The micronucleus assay using peripheral
blood reticulocytes from X-ray-exposed mice Mutat. Res. 335(3) 229–34
[38] Fenech M and Morley A A 1985 Measurement of micronuclei in lymphocytes Mutat.
Res. 147(1/2) 29–36
[39] Boyd E, Buchanan W W and Lennox B 1961 Damage to chromosomes by therapeutic
doses of radioiodine Lancet i 977–8
[40] M’Kacher R, Schlumberger M, Legal J D, Violot D, Beron-Gaillard N, Gaussen A
et al. 1998 Biologic dosimetry in thyroid cancer patients after repeated treatments
with iodine-131 J. Nucl. Med. 39(5) 825–9
[41] M’Kacher R, Legal J D, Schlumberger M, Voisin P, Aubert B, Gaillard N et al. 1996
Biological dosimetry in patients treated with iodine-131 for differentiated thyroid
carcinoma J. Nucl. Med. 37(11) 1860–4
[42] Tucker J D, Ramsey M J, Lee D A and Minkler J L 1993 Validation of chromosome
painting as a biodosimeter in human peripheral lymphocytes following acute
exposure to ionizing radiation in vitro Int. J. Radiat. Biol. 64(1) 27–37
[43] Schmid E, Zitzelsberger H, Braselmann H, Gray J W and Bauchinger M 1992
Radiation-induced chromosome aberrations analysed by fluorescence in situ
hybridization with a triple combination of composite whole chromosome-specific
DNA probes Int. J. Radiat. Biol. 62(6) 673–8
[44] Zanzonico P B 1997 Radiation dose to patients and relatives incident to 131I therapy
Thyroid 7(2) 199–204
[45] Smith T and Edmonds C J 1984 Radiation dosimetry in the treatment of thyroid
carcinoma by 131I Radiat. Prot. Dosim. 5 141–9
[46] Wuttke K, Streffer C, Muller W U, Reiners C, Biko J and Demidchik E 1996
Micronuclei in lymphocytes of children from the vicinity of Chernobyl before and
after 131I therapy for thyroid cancer Int. J. Radiat. Biol. 69(2) 259–68
[47] Watanabe N, Yokoyama K, Kinuya S, Shuke N, Shimizu M, Futatsuya R et al. 1998
Radiotoxicity after iodine-131 therapy for thyroid cancer using the micronucleus
assay J. Nucl. Med. 39(3) 436–40
284 Dose calculation procedures in patient-specific radionuclide therapy
[48] Ramirez M J, Surralles J, Galofre P, Creus A and Marcos R 1997 Radioactive iodine
induces clastogenic and age-dependent aneugenic effects in lymphocytes of thyroid
cancer patients as revealed by interphase FISH Mutagenesis 12(6) 44955
[49] Monsieurs M A, Thierens H M, van de Wiele C V, Vral A M, Meirlaen I A, de Winter
H A et al. 1999 Estimation of risk based on biological dosimetry for patients treated
with radioiodine Nucl. Med. Commun. 20(10) 911–17
[50] Livingston G K, Foster A E and Elson H R 1993 Effect of in vivo exposure to iodine-131
on the frequency and persistence of micronuclei in human lymphocytes J. Toxicol.
Environ. Health 40(2/3) 367–75
[51] Lambert V, Thierens H, Monsieurs M, Roncancio C and Laurent C 2001 Trans-
location frequencies measured in patients one year after radioactive iodine therapy
for thyrotoxicosis Int. J. Radiat. Biol. 77(6) 679–85
[52] Gutierrez S, Carbonell E, Galofre P, Creus A and Marcos R 1997 Micronuclei induc-
tion by 131I exposure: study in hyperthyroidism patients Mutat. Res. 373(1) 39–45
[53] Catena C, Villani P, Nastasi R, Conti D, Righi E, Salerno G et al. 1994 Micronuclei
and 3AB-index in patients receiving iodine-131 therapy J. Nucl. Biol. Med. 38(4)
58693
[54] McEwan A C 1977 Absorbed doses in the marrow during 131I therapy Br. J. Radiol.
50 329–31
[55] Goddu S M, Howell RW, Giuliani D C and Rao D V 1998 Biological dosimetry of
bone marrow for incorporated yttrium-90 J. Nucl. Med. 39(3) 547–52
[56] Goddu S M, Bishayee A, Bouchet L G, Bolch W E, Rao D V and Howell R W 2000
Marrow toxicity of 33P- versus 32P-orthophosphate: implications for therapy of bone
pain and bone metastases J. Nucl. Med. 41(5) 941–51
[57] Lenarczyk M, Goddu S M, Rao D V and Howell R W 2001 Biologic dosimetry of
bone marrow: induction of micronuclei in reticulocytes after exposure to 32P and
90
Y J. Nucl. Med. 42(1) 162–9
[58] Siegel J A, Wessels B, Watson E E, Stabin M, Vriesendorp H M, Bradley E W et al.
1990 Bone marrow dosimetry and toxicity for radioimmunotherapy Antibody
Immunocon. Radiopharm. 3 213–33
[59] Wong J Y, Wang J, Liu A, Odom-Maryon T, Shively J E, Raubitschek A A et al. 2000
Evaluating changes in stable chromosomal translocation frequency in patients
receiving radioimmunotherapy Int. J. Radiat. Oncol. Biol. Phys. 46(3) 599–607
[60] Gladstone D J, Lu X Q, Humm J L, Bowman H F and Chin L M 1994 A miniature
MOSFET radiation dosimeter probe Med. Phys. 21 1721–8
[61] Miller S D and Endres G W 1990 Laser-induced optically stimulated M centre
luminescence in LiF Radiat. Prot. Dosim. 33(1) 59–62
[62] Strandh M, Johansson J and Strand S E A 2001 BaF(Br,I)Eu dosimeter based on laser
induced fluorescence for in in vivo measurements (http://www.doseinfo-radar.com/1-
pos.htm)
[63] Strandh M, Strand S-E, Johansson J and Malmgren L 2002 Dosimetry based on
laser induced fluorescence in LiF for in vivo measurements (http://www.doseinfo-
radar.com/1-pos.htm)
[64] Schauer D A, Desrosiers M F, Le F G, Seltzer S M and Links J M 1994 EPR
dosimetry of cortical bone and tooth enamel irradiated with X and gamma rays:
Study of energy dependence Radiat. Res. 138 1–8
[65] Ignatiev E A, Lyubashevskii N M, Shishikina E A and Romanyukha A A 1999 EPR
dose reconstruction for bone-seeking 90Sr Appl. Radiat. Isot. 51(2) 151–9
References 285
13.1. INTRODUCTION
286
Introduction 287
Figure 13.1. Diagram of the upper half of a femur showing the trabecular bone
regions in the femur head (NMR image) and neck, and the cortical bone regions
in the diaphysis.
Trabecular bone
Bone trabeculae in parietal bone 511
Bone trabeculae in cervical vertebra 279
Bone trabeculae in femur head 232
Marrow cavity in parietal bone 389
Marrow cavity in cervical vertebra 910
Marrow cavity in femur head 1157
Cortical bone
Osteon length (between Volkmann canals) 2500
Concentric lamellae thickness 7
Osteon diameter 200–400
Haversian canal diameter 20–200
288 Dosimetry of skeleton and bone marrow
by Snyder et al. for the tabulation of dose conversion factors for various
radionuclides [18], and was subsequently used in MIRD Pamphlets Nos.
5 and 11 [14, 15]. Snyder et al. calculated the energy deposited by
photon interactions in the macroscopic skeletal model, and partitioned it
by target mass to generate the energy deposited in the sensitive tissues of
the skeletal system. This approach assumed that the range of secondary
electrons was smaller than the marrow cavity sizes and equivalence in
the mass absorption coefficients for marrow and bone. In an effort to
improve upon these assumptions, Eckerman presented a new approach in
which photon fluence calculated in a homogenous skeletal region of the
anthropomorphic phantom is coupled with dose per unit fluence response
factors derived considering the creation and transport of the secondary
electrons in a microscopic model of the skeletal region [16, 19, 20]. This
second approach was used in the tabulation of specific absorbed fractions
of energy (absorbed fraction of energy per unit mass of target region)
published in ORNL TM 8381 [16], which were subsequently used in
MIRDOSE 3 [21]. A similar approach proposed by Bouchet et al. consisted
of calculating the energy of the secondary electrons created by photon
interactions in the macroscopic model of the skeletal system and coupling
them to electron absorbed fractions of energy calculated using a
microscopic model of the skeletal system [22]. The last two methods
differed principally in the use of two different microscopic models of the
skeletal system that will be presented and compared in subsequent
sections.
For directly ionizing radiations, there are essentially two approaches to
calculate the energy deposited in the sensitive regions of the skeletal system.
The first is to focus the modelling on closely representing the anatomical
microstructures of the skeletal system and to couple the model to a Monte
Carlo particle transport code. However, because of its intricate anatomy
and its complex microstructure, it is difficult to exactly model the skeletal
system without using three-dimensional imaging, and only a few investiga-
tors have used three-dimensional data. The second approach is to model
the particle transport in such a way that the traversed distance across a
given microscopic region (bone trabeculae or marrow cavity for example)
by a particle matches exactly the real path of the particle in the skeletal
microstructure. This approach initially was proposed by the research
group of Spiers [6, 23–29] and subsequently was used in several skeletal
dosimetry models [19, 20, 30].
In the next sections of this chapter, the different microscopic models of
the cortical and trabecular bone regions that have been used with Monte
Carlo techniques to generate dosimetric data for nuclear medicine are
presented. Because a model is based upon specific assumptions, and is
derived for a reference individual, correction methods for patient-specific
applications are given as well.
290 Dosimetry of skeleton and bone marrow
representation of the arm [18]. In addition, a fixed point source of 1.0 MeV
monoenergetic electrons was studied.
This spherical model was improved by Zuzarte de Mendoca by
consideration of the endosteal layer as well as the marrow space and
trabecular dimensions [34]. She used the Monte Carlo code EGS4 [35] and
considered both surface and volume sources. These calculations produced
absorbed fractions and specific absorbed fractions of energy for six mono-
energetic electron sources (0.1, 0.2, 0.5, 1.0, 2.0 and 4.0 MeV) for source–
target combinations in the spine, the cranium and the pelvis. She found
that the absorbed fractions in bone marrow were strongly dependent on
the ratio of trabecular bone to bone marrow volume. The higher the ratio,
the lower the absorbed fraction in bone marrow.
Parry [36] continued this research and devised a modified model of the
trabecular bone for use in the EGS4 code [35]. He calculated electron
absorbed fractions for red marrow and the endosteal tissue on trabecular
bone. The model was based on measured omnidirectional chord length
distributions through trabeculae and marrow cavities in nine skeletal
locations. Absorbed fractions were calculated for only two target regions
and all other skeletal regions were estimated based on these results. In
addition, the absorbed fractions were used to calculate S values for a variety
of -emitting radionuclides at each of 15 skeletal regions containing red bone
marrow. Absorbed fractions and S values were calculated for sources in the
marrow, on the bone surface and in the bone volume. Comparisons were
made between these results and those published in ICRP Publication 30
[32] as well as those determined using MIRDOSE2 software [21].
total marrow mass that is active). This model simulated no escape of the
electron from the trabecular bone, and assumed that all energy emitted
was deposited in the trabecular bone. In 1994, results from this model
were implemented in the MIRDOSE3 program [21]. Reference masses for
the endosteum and active marrow in trabecular bone of the 15 skeletal
regions defined in the anthropomorphic phantoms [16] were derived as
described in [20] and used to convert Eckerman’s absorbed fractions of
energy to radionuclide S values. This implementation allowed for explicit
calculation of regional absorbed dose in the 15 different bone sites as well
as for the whole skeleton, for two trabecular bone target regions (red
marrow and the endosteum) and three trabecular bone sources (trabecular
bone volume, red marrow and trabecular bone surface).
was used to derive electron absorbed fractions of energy for two sources
(Haversian canal and bone surface) and three targets (Haversian canal,
endosteum and the bone matrix). From these calculations, dose conversion
factors were tabulated for these seven radionuclides as a function of the
Haversian canal radius. Although elegant in its design, this anatomical
model only considered a single size for the Haversian canal, and it did not
consider the cross-osteon irradiation.
[60]. Nonetheless, both external beam and internal emitter therapy share the
same uncertainty regarding predicted patient tolerance. Nuclear medicine
physicians in general have less experience in the therapy regime than those
who deal with external treatment approaches, and typically undertreat
patients in order to minimize the risk of even low-grade toxicity. If treatment
of cancer with radionuclides and RIT are to become a primary means of
treating cancer patients, the activity to be administered to individual patients
should be based on accurate evaluations of the dose to be delivered to the
tumour and of patient tolerance [61].
Radionuclide therapy has been most effective in the treatment of thyroid
cancer because of the high specific uptake of radioactive iodine by thyroid
tissue, but there is still no agreement on the means for choosing an optimal
therapeutic radiation dose to be delivered. The hope for RIT derives from the
notion that it can achieve selective high tumour uptake based on improved
targeting of tumour specific/associated antigens. Intravascular tumours,
including leukaemia and lymphoma, have the highest likelihood of successful
treatment with RIT, partly because tumour targeting efficiency is high, and
partly because of the high radiation sensitivity of lymphocytes, and to a
lesser extent of granulocytes. A major problem with RIT treatment of
most solid tumours has been low uptake per unit tumour mass, so that the
delivered dose is often too low for successful therapy. If one gives very
large doses intending to ablate marrow and tumour followed by marrow
transplantation, then lungs, liver or kidneys will become the dose limiting
organs. Bone marrow transplantation (BMT) is being employed routinely
in the treatment of breast cancer with high doses of chemotherapy and/or
of some solid tumours with RIT. A reasonable goal would be to develop
an improved means of defining an appropriate intermediate dose that
provides high ‘cure’/complete remission rates, without the need for BMT.
Progress is being made in the development of improved calculation
methods and tools for patient-specific treatment-planning focusing primarily
on tumour dose optimization [62–64]. Still, problems remain in the accurate
assessment of bone marrow dose and in the estimation of marrow reserve.
Better patient-specific knowledge of both dose and marrow reserve is
needed if one is to develop a good predictor of tolerance dose. Two analytical
approaches have been taken. The first is the physical method, in which one
measures the biokinetics and computes dose based on knowledge of the
amount injected, the radionuclide emission characteristics, and the measured
residence times in the different organs using accepted dose calculation
methods [12]. A complementary approach attempts to use patient-specific
characteristics to estimate residual damage from prior therapeutic treatments
to estimate the amount of added dose that would be tolerated by normal
organs, including bone marrow.
The importance of including pre-therapy patient status in toxicity
predictions has been studied by several groups [59, 65], and it is clear that
Improving the models: patient-specific dosimetry 301
prior patient therapy and disease status must be factored into treatment
decisions. Nonetheless, to date poor correlations have been observed
between administered dose and marrow toxicity, which are somewhat
improved by more carefully calculated marrow doses [54]. This suggests
the need for more patient-specific information, including better dose
estimates and better information on other patient characteristics, including
perhaps patient-specific indicators of response (i.e., biomarkers).
Marrow toxicity from RIT is manifest by haematological changes in
circulating platelets, lymphocytes, granulocytes and reticulocytes, and later
changes are seen in red blood cells. Due to the high sensitivity of the
marrow, second organ toxicity is not a problem unless doses which are
known to be lethal to the marrow are intentionally given. Attempts to
correlate haematological toxicity with marrow dose, when marrow cells are
specifically targeted, have not been particularly successful in the past, in part
due to uncertainties in the actual absorbed dose, but also due to the
difficulty in assessing marrow functional status prior to therapy [8, 13, 53–
58]. The best correlations obtained to date have had a value of r around 0.8
[53, 57, 66, 67]. Whole-body absorbed dose and red marrow absorbed dose
are usually the best indicators of haematological toxicity, as measured by
platelet toxicity grade, with red marrow dose being slightly better. However,
the poor correlation coefficients obtained from these analyses emphasize the
need to better account for prior patient therapy and bone marrow reserve.
That red marrow absorbed dose is the best (albeit not appropriate on a
patient-by-patient basis) predictor of haematological toxicity has been
confirmed using data from other antibody trials [68]. Most recently, a
survey of the literature using data from different antibodies, different patient
populations and different disease sites has confirmed this observation over a
much larger scale [69]. Blumenthal et al. showed that plasma FLT3-L (a
stromal-cell produced positive stimulatory cytokine) levels predicted excess
platelet toxicity in patients receiving RIT after chemotherapy [70].
Many other studies have found little or no correlation between calculated
marrow dose for large numbers of patients and observed toxicity, as deter-
mined by various haematological parameters (platelet nadir, percentage
platelet drop, percentage lymphocyte drop, percentage granulocyte drop,
etc.). The primary reason for this is probably that all available marrow dose
models are patterned after a normal healthy young adult with good marrow
reserve, while the treated population is generally older (with inherently
different trabecular bone structures and less marrow reserve), often with
compromised marrow due to disease, prior chemotherapy or radiotherapy.
where CFbone is the cellularity factor for a given bone, fRM;bone is the
fraction of the red marrow (RM) associated with a given bone, and
AFðMS RMÞbone is the absorbed fraction for the energy originating in
the marrow of a given bone and irradiating the marrow space of that
bone. The values of fRM;bone and CFbone may be assigned based on infor-
mation gained in individual radionuclide marrow scans; that is, the fraction
of marrow calculated to be in each of the model regions and the cellularities
may be calculated on a patient-specific basis. In the absence of patient-
specific values of cellularity derived from marrow biopsies, values from
age-dependent models (e.g., [72, 73]) may be applied or, as will be discussed
in the next section, magnetic resonance imaging methods could also be used
to determine patient-specific active marrow mass and location [74–77]. Bone
marrow uptake of RIT agents in expected marrow sites which have
diminished or absent marrow function on radionuclide scans may be
assumed to be tumour, and the computed dose calculated using a more
specific model (e.g., [41]). Time–activity integrals may be calculated and
doses estimated for the different bone marrow regions in which marrow
localization has been scored, assigning the total marrow residence time to
Improving the models: patient-specific dosimetry 303
each region using the values of fRM;bone . Whole marrow weighted average
dose and differential and integral dose–volume histograms then may be
calculated.
Determining marrow time–activity integrals for internal emitters is an
area of difficulty, as marrow activity is heterogeneous even in normal
subjects, and more so in patients with cancer, with or without a history of
previous treatment. Different imaging methods rely on either (i) the quanti-
fication of activity in a certain segment of the marrow (e.g., the lumbar spine
or sacrum) and extrapolation of total marrow activity assuming the total
amount of red marrow in the body and the fraction of the total that resides
in individual segments [13], or (ii) the assumption that the dose is derived
only from blood contributions, with the assumption of normal marrow
size [9] and/or (iii) combined blood and marrow-fixed activity methods.
The correlation between dose calculated in the different manners, biomarker
results, and clinical predictors with observed toxicity and patient cure rates
should be used to assess the potential utility of the different dose calculation
and marrow indices for improving treatment outcome. This is not currently
in standard practice, but should be considered if we are to seek more
meaningful dose and risk assessment in individual patients, as is routine in
external beam therapy.
13.5. SUMMARY
REFERENCES
[35] Nelson W R, Hirayama R H and Roger D W O 1985 The EGS4 code system SLAC
Report 265 (Stanford, CA: Stanford Linear Accelerator Center)
[36] Parry R A 1995 Revised Estimates of Electron Absorbed Fractions and Radionuclide S-
Values in Trabecular Bone (College Station, TX: Texas A&M University)
[37] Bielajew A F and Rogers D W O 1987 PRESTA, the ‘parameter reduced electron
step transport algorithm’ for electron Monte Carlo transport Nucl. Instrum. Meth. B
18 165–81
[38] Bouchet L G, Bolch W E, Howell R W and Rao D V 2000 S values for radionuclides
localized within the skeleton J. Nucl. Med. 41 189–212
[39] Patton P 2000 NMR Microscopy for Skeletal Dosimetry: an Investigation of Marrow
Cellularity on Dose Estimates (Gainesville, FL: University of Florida)
[40] Stabin M, Eckerman K F, Bolch W E, Bouchet L G and Patton P W 2002 Evolution
and status of bone and marrow dose models Cancer Biother. Radiopharmaceut. In press
[41] Samaratunga R C, Thomas S R, Hinnefeld J D, Von Kuster L C, Hyams D M,
Moulton J S, Sperling M I and Maxon H R 1995 A Monte Carlo simulation
model for radiation dose to metastatic skeletal tumor from rhenium-186(Sn)-
HEDP J. Nucl. Med. 36 336–50
[42] Jokisch D W, Patton P W, Inglis B A, Bouchet L G, Rajon D A, Rifkin J and Bolch
W E 1998 NMR microscopy of trabecular bone and its role in skeletal dosimetry
Health Phys. 75 584–96
[43] Jokisch D W, Bouchet L G, Patton P W, Rajon D A and Bolch W E 2001 Beta-
particle dosimetry of the trabecular skeleton using Monte Carlo transport within
3D digital images Med. Phys. 28 1505–18
[44] Bolch W E, Patton P W, Rajon D A, Shah A P, Jokisch D W and Inglis B A 2002
Considerations of marrow cellularity in 3-dimensional dosimetric models of the
trabecular skeleton J. Nucl. Med. 43 97–108
[45] Patton P W, Rajon D A, Shah A P, Jokisch D W, Inglis B A and Bolch W E 2002 Site-
specific variability in trabecular bone dosimetry: considerations of energy loss to
cortical bone Med. Phys. 29 6–14
[46] Patton P W, Jokisch D W, Rajon D A, Shah A P, Myers S L, Inglis B A and Bolch
W E 2002 Skeletal dosimetry via NMR microscopy: investigations of sample
reproducibility and signal source Health Phys. 82 316–26
[47] Rajon D A, Jokisch D W, Patton P W, Shah A P and Bolch W E 2000 Voxel size
effects in three-dimensional nuclear magnetic resonance microscopy performed for
trabecular bone dosimetry Med. Phys. 27 2624–35
[48] Bouchet L G and Bolch W E 1999 A three-dimensional transport model for deter-
mining absorbed fractions of energy for electrons within cortical bone J. Nucl. Med.
40 2115–24
[49] Akabani G 1993 Absorbed dose calculations in Haversian canals for several beta-
emitting radionuclides J. Nucl. Med. 34 1361–6
[50] Breems D A, Blokland E A, Neben S and Ploemacher R E 1994 Frequency analysis of
human primitive haematopoietic stem cell subsets using a cobblestone area forming
cell assay Leukemia 8 1095–104
[51] Neben S, Anklesaria P, Greenberger J and Mauch P 1993 Quantitation of murine
hematopoietic stem cells in vitro by limiting dilution analysis of cobblestone area
formation on a clonal stromal cell line Exp. Hematol. 21 438–43
[52] Bolch W E, Bouchet L G, Robertson J S, Wessels B W, Siegel J A, Howell R W,
Erdi A K, Aydogan B, Costes S, Watson E E, Brill A B, Charkes N D, Fisher D R,
308 Dosimetry of skeleton and bone marrow
[68] Sgouros G, Deland D, Loh A C, Divgi C R and Larson S M 1997 Marrow and whole-
body absorbed dose vs marrow toxicity following 131I-G250 antibody therapy in
patients with renal-cell carcinoma J. Nucl. Med. 38 252P
[69] O’Donoghue J A, Baidoo N and Sgouros G 1998 Calculated red marrow absorbed
dose predicts hematological toxicity in radioimmunotherapy J. Nucl. Med. 39 186P
[70] Blumenthal R D, Lew W, Juweid M, Alisauskas R, Ying Z and Goldenberg D M 2000
Plasma FLT3-L levels predict bone marrow recovery from myelosuppressive therapy
Cancer 88 333–43
[71] Shields A F, Grierson J R, Dohmen B M, Machulla H J, Stayanoff J C, Lawhorn-
Crews J M, Obradovich J E, Muzik O and Mangner T J 1998 Imaging
proliferation in vivo with [F-18]FLT and positron emission tomography Nat. Med.
4 1334–6
[72] Mauch P, Constine L, Greenberger J, Knospe W, Sullivan J, Liesveld J L and Deeg
H J 1995 Hematopoietic stem cell compartment: acute and late effects of radiation
therapy and chemotherapy Int. J. Radiat. Oncol. Biol. Phys. 31 1319–39
[73] Koo K H, Dussault R, Kaplan P, Kim R, Ahn I O, Christopher J, Song H R and
Wang G J 1998 Age-related marrow conversion in the proximal metaphysis of the
femur: evaluation with T1-weighted MR imaging Radiology 206 745–8
[74] Ballon D, Dyke J, Schwartz L H, Lis E, Schneider E, Lauto A and Jakubowski A
2000 Bone marrow segmentation in leukemia using diffusion and T2 weighted echo
planar magnetic resonance imaging NMR Biomed. 13 321–28
[75] Ballon D, Jakubowski A A, Tulipano P K, Graham M C, Schneider E, Aghazadeh B,
Chen Q S and Koutcher J A 1998 Quantitative assessment of bone marrow hemato-
poiesis using parametric magnetic resonance imaging Magn. Reson. Med. 39 789–800
[76] Ballon D, Jakubowski A A, Graham M C, Schneider E and Koutcher J A 1996
Spatial mapping of the percentage cellularity in human bone marrow using
magnetic resonance imaging Med. Phys. 23 243–50
[77] Ballon D, Jakubowski A, Gabrilove J, Graham M C, Zakowski M, Sheridan C and
Koutcher J A 1991 In vivo measurements of bone marrow cellularity using volume-
localized proton NMR spectroscopy Magn. Reson. Med. 19 85–95
[78] Link T M, Majumdar S, Grampp S, Guglielmi G, van Kuijk C, Imhof H, Glueer C
and Adams J E 1999 Imaging of trabecular bone structure in osteoporosis Eur.
Radiol. 9 1781–8
[79] Saito Y and Ortawara B S 2001 Development and evaluation of a real-time 3-
dimensional CT (4D-CT) scanner. In 87th Annual Meeting of the Radiological Society
of North America Chicago, IL
[80] Nonomura Y, Yasumoto M, Yoshimura R, Haraguchi K, Ito S, Akashi T and Ohashi
I 2001 Relationship between bone marrow cellularity and apparent diffusion coeffi-
cient J. Magn. Reson. Imaging 13 757–60
[81] Schick F, Einsele H, Kost R, Duda S, Jung W I, Lutz O and Claussen C D 1994
Hematopoietic reconstitution after bone marrow transplantation: assessment with
MR imaging and H-1 localized spectroscopy J. Magn. Reson. Imaging 4 71–8
[82] Schick F, Einsele H, Weiss B, Forster J, Lutz O, Kanz L and Claussen C D 1996
Assessment of the composition of bone marrow prior to and following autologous
BMT and PBSCT by magnetic resonance Ann. Hematol. 72 361–70
Chapter 14
14.1. INTRODUCTION
310
Introduction 311
formation of a fibrous mesh of platelets and red blood cells. Monocytes then
phagocytize the mesh, and smooth muscle cells migrate into the intimal
subendothelial space, leading to neointimal formation. Over weeks, the
neointima becomes less cellular and the healing site resembles a plaque. In
about 40% of patients, the neointimal hyperplasia is excessive and results
in more narrowing of the lumen than was resolved by angioplasty [1].
Various preclinical and clinical studies have shown that the delivery of
between perhaps 15 and 50 Gy of ionizing radiation to the coronary artery
wall after balloon angioplasty can markedly reduce the rate of restenosis,
with few or no complications [1]. Due to the geometry, internal sources are
thought to be more efficacious than external beams of radiation, thus a
science of ‘intravascular brachytherapy’ (IVB) or ‘intravascular radionuclide
therapy’ is under development. This form of therapy principally takes two
approaches:
1. The use of temporary implants—seeds, wires and balloons with radio-
active walls or filled with radioactive solutions to deliver between 15
and 20 Gy to a 2–3 cm length of vessels between 2 and 5 cm in diameter.
The desired dose rate is about 5 Gy/min, to deliver the dose in a short
time and reduce complications from reduced blood flow. The technique
of afterloading with high dose rate sources is common.
2. The use of permanent implants—i.e., radioactive stents, that deliver a
similar dose, but with a lower dose rate.
The effectiveness of using high dose-rate intravascular radioisotope sources
for the inhibition of arterial restenosis has been demonstrated in a variety
of animal models and clinical trials [2]. This technology has been applied
primarily to the coronary arteries for the inhibition of intimal hyperplasia
following balloon angioplasty and/or stent placement. This neointimal
hyperplasia is considered to be pathologic, as the narrowing or restenosis
of a coronary artery can be life threatening [3, 4]. The current concept of
the mechanism of action is based on the inhibition of myofibroblast
stimulation in the adventitia by radiation [3]. Vessel wall irradiation with
particles has been demonstrated to prevent the onset of restenosis after
stent implantation or angioplasty in a porcine model [5-8].
Calculation of the radiation dose delivered to the coronary artery wall, as
well as to other organs and tissues in the body, is important to fully evaluate
the possible risks and benefits of this procedure. Computational methods for
these calculations are well developed and simply need to be implemented
appropriately. Both analytical and Monte Carlo methods are available and
have been used to calculate dose distributions around these sources. Monte
Carlo methods have the advantage of explicit treatment of the radiation
transport of electrons and photons at interfaces of materials of different
density and composition (e.g., at the interface of an atherosclerotic plaque
and soft tissue). The different approaches to the delivery of this radiation
312 Monte Carlo modelling of dose distributions
dose include the use of radioactive wires, balloons with radioactive walls or
filled with radioactive solutions, and stents. There are also numerous candidate
radionuclides, including and emitters, of low and high energy, which are
currently under investigation.
the use of liquid-filled balloons has been proposed, as this inherently solves
the problem of source uniformity. In this application, however, there is the
risk of balloon rupture and release of the radioactive liquid into the blood-
stream. This would presumably act like an intravenous injection of the
material, with its subsequent distribution, retention and radiation dosimetry
being determined by the radionuclide and chemical compound (if any) to
which it is attached. The use of a balloon with radioactive walls can give
the same uniformity with a much lower risk of releasing the radioactive
materials.
The use of radioactive stents has been proposed, to provide structure
and stability to the affected artery as well as delivering the radiation dose
to the vessel wall. In these applications, to date, only emitters (principally
32
P) have been proposed. The activity may be coated on the surface of the
stent (which raises the possibility of activity leaching over time into the
bloodstream) or impregnated into the metal of the stent (by a manufacturing
process or neutron activation of extant metal components). These stents have
the advantages that low activities (kBq to MBq range) can be employed (as
the dose is delivered over many hours or days as the radionuclide decays in
situ), and that the stent conforms in shape to the lesion site. Non-radioactive
stents, however, do not inhibit the process of restenosis; in fact, there is
evidence that they may even stimulate it [1].
These are not the only candidate technologies and nuclides available or
suggested. Many strategies and nuclides have been proposed, including low
energy X- or -ray emitters, the use of inserted X-ray generation devices,
and others. As many ideas are being generated and tested, one cannot
generalize about the future. Many factors, scientific, economic and logistic,
will influence the final suitability of different technologies and will determine
those that find the greatest acceptance.
particles, etc.). Point kernel methods generally apply to only a single medium
(generally water), but can be adapted to approximate particle behaviour in
other media [9,10]. Monte Carlo methods have the advantage of explicitly
treating particle behaviour in various media. For example, in this applica-
tion, a fluid in a filled balloon, the wall of the balloon, and the vessel wall
will all have similar compositions (water or soft tissue equivalent), but a
central guide wire or source wire will be different (metallic) and the plaque
will also be different (calcified tissue, with a composition and density signifi-
cantly different from that of soft tissue). Monte Carlo methods can have the
disadvantage of only being able to treat certain predefined geometries in
combinations (e.g., cylinders, spheres).
Most early attempts to characterize the dosimetry of various intra-
vascular sources employed point-kernel approaches, while later ones
employ more Monte Carlo applications. Figure 14.1 shows isodose distribu-
tions, based on point-kernel calculations, for a linear array of 192 Ir seeds
(adapted from a study by Wiedermann et al. [8], as given in AAPM
Report No 66 [1]). Prestwich et al. [11] used analytical kernels, which were
based ultimately on measurements, to calculate dose from a 32 P coated
stent. They showed one-dimensional radial dose distributions; two-dimen-
sional isodose contours and three-dimensional dose distributions for a
3 mm 20 mm cylindrical source. Li et al. [12] used MCNP to calculate
doses from a single strut of a stent with radioisotope (48 V, 32 P) distributed
throughout, then used superposition to calculate dose distributions for an
Radiation dosimetry studies 315
entire stent containing 48 V and 32 P (figure 14.2). Janicki et al. [13] used a more
realistic geometry but still used an analytic kernel, and calculated three-
dimensional dose distributions in tissue around a 3:5 mm 20 mm wire
mesh stent. Later, Janicki et al. [9] extended this model to cases in which
the stent was surrounded by layers of materials with different compositions
and densities (figure 14.3). Duggan et al. [14] performed a calculation similar
to that of Prestwich et al. but with a Monte-Carlo generated kernel from
Simpkin et al. [15]. McLemore [16] also used MCNP to calculate dose
from a single strut with a 103 Pd coating on the outside and then used super-
position to calculate the dose from the entire stent.
Many complete sources, with realistic geometries and materials, have
now been simulated. These include entire radioactive stents [17]. Seed or
wire sources are usually encapsulated, often in high atomic number
materials. This can have a dramatic effect on the dose distribution, especially
for or low-energy photon sources. Monte Carlo simulation is the most
accurate technique for predicting the dose distribution for such sources.
Examples of simulations of encapsulated sources that are now in clinical
use include that of the Novoste BetaCath 90 Sr seed by Soares et al. [18]
and the Guidant encapsulated 32 P wire source by Mourtada et al. [19]
(figure 14.4). In addition, many other proposed sources have been simulated,
including a coil-shaped wire source [20] and a soft X-ray device [21].
As mentioned above, an important issue in the use of radioactive
wires or seeds is the centring of the source within the vessel. A study done
by Amols et al. [22] showed that centring is substantially more important
for sources than for photon sources (although it is clearly important for
both). Figure 14.5 shows that a centring error of 0.5 mm in a 5 mm diameter
316 Monte Carlo modelling of dose distributions
Figure 14.3. Comparison of calculations with all water DPK model (W) and
multilayer DPK model (ML) with radiochromic film measurements (blue)
for a 32 P stent in a block of solid water with a layer of Teflon in between the
stent and the radiochromic film. Reprinted with permission from AAPM
(Janicki et al [9]).
vessel will result in an error of 1.6 for 192 Ir (i.e., the actual dose received on
the overdosed side will be 1.6 times that calculated for a centred source) and
2.1 for 90 Y (figure 14.5).
Extensive calculations, using both the MCNP 4B [23] and EGS4 [24]
Monte Carlo codes, were made for a number of fixed artery sizes and for
Figure 14.4. Comparison of Monte Carlo and measured values for the Guidant
encapsulated 32 P wire source [19]. Reprinted with permission from AAPM.
Radiation dosimetry studies 317
wire and liquid-filled balloon sources by Stabin et al. [25]. They performed
calculations for monoenergetic photons and electrons over a broad range
of energies, then applied decay spectra for over 800 radionuclides to the
monoenergetic results, and incorporated the data into an interactive PC
computer program that provided rapid comparisons of candidate radio-
nuclides and geometries. Their blood vessel models were simple cylinders
3 cm in length, half of whose inner circumference had an atherosclerotic
plaque and half which did not. Results were given for radial dose as a
function from the vessel wall (diameters of 1.5, 3 and 4.5 mm), and compar-
isons for some candidate nuclides were given (figure 14.6).
Figure 14.6. Dose comparison for 186 Re and 188 Re balloons and an 192 Ir wire in
a 3 mm diameter vessel. In all cases, the activities were calculated so that the
vessel would receive a dose of 30 Gy in 5 min at 0.5 mm from the vessel surface
(comparisons performed on the non-plaque side of the vessel). The calculated
activities were 2.48 GBq of 186 Re, 1.02 GBq of 188 Re and 64.1 GBq of 192 Ir.
Figure 14.7. Diagram of solid water phantom used to expose radiochromic film
to the -emitting stent. Four stents were exposed. For two stents, a layer of
Teflon was placed on the stent before the radiochromic film was stacked on
top. Low and high refer to low- and high-sensitivity radiochromic film, respec-
tively. The measured thicknesses of the different layers were Teflon, 386 mm;
low-sensitivity radiochromic film, 108 mm; and high-sensitivity radiochromic
film, 257 mm. Reprinted with permission from Janicki et al [9].
Radiation dosimetry studies 319
Figure 14.9. Radiochromic film results. Reprinted with permission from Janicki
et al. [9].
320 Monte Carlo modelling of dose distributions
Figure 14.10. Comparison of the radiochromic film results and the DPK multi-
layer (ML) model. Reprinted with permission from Janicki et al. [9].
Publication 30 [26]; the model for perrhenate was the MIRD Dose Estimate
Report No 8 model (Lathrop et al. [27]) for sodium pertechnetate in nonrest-
ing subjects, which was simply adapted for use with the rhenium isotopes, as
suggested by the data of Kotzerke et al. [28]. Organ doses for the four radio-
pharmaceuticals assuming release of unit quantities into the bloodstream are
shown in table 14.1. For 32 P and 90 Y, the organ receiving the highest dose is
red marrow with a dose of 8.12 and 3.26 Gy per GBq respectively. With 186 Re
and 188 Re, the organs receiving the highest dose are the large intestine,
thyroid and urinary bladder.
Calculations show that one would require 1.04 GBq of 32 P or 0.92 GBq
90
of Y in a balloon source to give 30 Gy in 5 min at 0.5 mm in a 3 mm
diameter artery. It is clear from these data and the dose estimates in the
above table that use of 32 P or 90 Y in the forms assumed here (liquid sources
of 32 P as phosphate or ionic 90 Y) will represent a significant risk to the
patient, as the predicted absorbed doses to red marrow from the quantities
of activity assumed are near or above median lethal doses. For the rhenium
isotopes, there are significant doses predicted to the intestines, but not at
levels that would be considered life threatening. Other strategies may
improve the radiation dose picture, such as the use of perchlorate to block
thyroid uptake, or attaching the rhenium isotopes to other chemical
compounds, such as MAG3, and to obtain different biological behaviour.
In these cases, obviously different biokinetic behaviour will result in different
dose estimates, which need to be calculated in each case to evaluate the
possible radiation doses.
Radiation effects—artery walls 321
Table 14.1. Organ doses from injection of radionuclides into the bloodstream
due to balloon rupture (Gy per GBq injected).
32 90 186 188
P Y Re Re
As this is a new technology and some of the effects we might expect are
expressed over longer times, there will be a delay before we know what
complications to expect. However, there is some clinical experience with
arterial damage and complications from external beam radiation. Acute
and chronic morphological changes have been seen in femoral arteries in
dogs given 40 Gy in 10 days [29]. In one radiation therapy centre, 20 patients
had blood vessel-associated injuries associated with radiation therapy over a
20 year period, including arterial rupture and occlusion, over a 7–24 year
period after therapy [30]. To effectively employ this promising technology,
we must make use of the best information that we can obtain from
322 Monte Carlo modelling of dose distributions
REFERENCES
[1] Task Group 60 of the American Association of Physicists in Medicine 1999 AAPM
Report No 66, Intravascular brachytherapy physics Med. Phys. 26 119–52
[2] Waksman R 1999 Vascular Brachytherapy 2nd edition (Armonk, NY: Futura)
[3] Wilcox J, Waksman R, King S and Scott N 1996 The role of the adventitia in the
arterial response to angioplasty: the effect of intravascular radiation Int. J. Radiat.
Oncol. Biol. Phys. 36 789–96
[4] Shi Y, Pieniek M, Fard A, O’Brien J, Mannion J D and Zalewski A 1996 Adventitial
remodeling after coronary arterial injury Circulation 93 340–8
[5] Waksman R, Robinson K, Croker I, Gravanis M B, Palmer S J, Wang C, Cipolla G D
and King S B III 1995 Intracoronary radiation before stent implantation inhibits
neointima formation in stented porcine coronary arteries Circulation 92 1383–6
[6] Waksman R, Robinson K, Croker I, Wang C, Gravanis M B, Cipolla G D, Hillstead
R A and King S B III 1995 Intracoronary low-dose beta-irradiation inhibits
neointima formation after coronary artery balloon injury in the swine restenosis
model Circulation 92 3025–31
[7] Waksman R, Robinson K, Croker I, Gravanis M B, Cipolla G D and King S B III
1995 Endovascular low-dose irradiation inhibits neointima formation after coronary
artery balloon injury in swine Circulation 91 1533–9
[8] Wiedermann J, Marboe C, Amols H, Schwartz A and Weinberger J 1994 Intra-
coronary irradiation markedly reduces restenosis after balloon angioplasty in a
porcine model J. Am. Coll. Cardiol. 23 1491–8
[9] Janicki C, Duggan D M, Gonzalez A, Coffey C W II and Rahdert D A 1999 Dose
model for a beta-emitting stent in a realistic artery consisting of soft tissue and
plaque Med. Phys. 26 2451–60
[10] Janicki C, Duggan D M and Rahdert D A 2001 A dose-point-kernel (DPK) model
for a low energy gamma-emitting stent in an heterogeneous medium Med. Phys. 28
1397–405
[11] Prestwich W V, Kennett T J and Kus F W 1995 The dose distribution produced by a
32P-coated stent Med. Phys. 22 313–20
[12] Li A N, Eigler N L, Litvack F and Whiting J S 1998 Characterization of a positron
emitting V48 nitinol stent for intracoronary brachytherapy Med. Phys. 25 20–8
[13] Janicki C, Duggan D M, Coffey C W, Fischell D R and Fishchell T A 1997
Radiation dose from a phosphorus-32 impregnated wire mesh vascular stent Med.
Phys. 24 437–45
[14] Duggan D M, Coffey C W and Levit S 1998 Dose distribution for a 32P-impregnated
coronary stent: comparison of theoretical calculations and measurements with radio-
chromic film Int. J. Radiat. Oncol. Biol. Phys. 40 713–20
[15] Simpkin D J and Mackie T R 1990 EGS4 Monte Carlo determination of the betadose
kernel in water Med. Phys. 17 179–86
[16] McLemore L 2000 Dosimetry characterization of a palladium 103 implanted stent for
intravascular brachytherapy, Master’s Thesis, University of Texas
References 323
15.1. INTRODUCTION
324
Introduction 325
[3], and at BNL [4] and the High Flux Reactor, Petten, The Netherlands [5]
for glioblastoma. These clinical trials all use a nuclear reactor as an epither-
mal neutron source. Accelerator-based BNCT for glioblastoma is also under
development [14–17].
There are several important differences between BNCS and BNCT.
First, the target tissue in BNCS is the synovium, which is only 0.4 to
1.3 cm below the skin surface (depending on the affected joint), while brain
tumours treated by BNCT can occur at depths up to 7 cm. Second, in vivo
10
B uptake studies have indicated that boron levels in the synovium on the
order of thousands or tens of thousands of ppm are readily achievable [7].
In BNCS, the boron-labelled compound is injected directly into the fluid
adjacent to the target tissue. In BNCT, conversely, the boronated compound
is delivered systemically, so 10 B levels in the target tissue are significantly
lower, usually on the order of 30–50 ppm. Finally, the joints to which
BNCS would be most frequently applied (knee and finger joints) are not
located near any radiosensitive organs of the body, whereas this is not the
case in BNCT clinical trials in the brain. All of these factors result in signifi-
cantly different neutron beam requirements for BNCS compared with
BNCT.
Figure 15.3. Two-dimensional plot through the MCNP model of the BNCS
moderator/reflector assembly including target, cooling apparatus, phantom
and graphite side and back reflectors around the phantom.
Table 15.2. BNCS therapy parameters for human knee and human finger
joints. These results are for the optimized configuration shown in figure 15.3,
including side and back reflectors, assuming 1000 ppm 10 B uptake in the syno-
vium, and parallel-opposed irradiations.
Human knee
2.6 MeV 9 Be(d,n) 23 13 335 38
1.5 MeV 9 Be(d,n) 20 48 586 22
4.0 MeV 9 Be(p,n) 15 7.3 203 73
3.7 MeV 9 Be(p,n) 15 10 190 79
2.5 MeV 7 Li(p,n) 15 7 129 108
Human finger
2.6 MeV 9 Be(d,n) 23 6.1 227 42
1.5 MeV 9 Be(d,n) 20 24 406 24
4.0 MeV 9 Be(p,n) 15 3.6 161 62
3.7 MeV 9 Be(p,n) 15 4.8 136 69
2.5 MeV 7 Li(p,n) 15 3.3 91 108
The neutron beamline assembly shown in figure 15.4 was used in a dose
escalation study to characterize the effects of BNCS in an animal model
(see section 15.5). Prior to the initiation of these animal experiments,
however, the MCNP predictions of the neutron and photon output of the
BNCS moderator/reflector assembly were experimentally verified. Accurate
knowledge of the thermal neutron flux in the phantom is important to
quantify the boron dose and thus treatment time. Photon and fast neutron
doses are also important in terms of dose to healthy structures such as the
Experimental characterization of the BNCS beam 335
Figure 15.5. Comparison of experimental and simulation results for (a) the
thermal neutron dose rate and (b) the fast neutron dose rate for a BNCS
target position of 23 cm. Simulation and experiment have been normalized at
a depth of 4 cm.
Experimental characterization of the BNCS beam 337
Figure 15.6. Comparison of experimental and simulation results for (a) the
thermal neutron dose rate and (b) the fast neutron dose rate for a BNCS
target position of 8 cm. Simulation and experiment have been normalized at a
depth of 4 cm.
338 The Monte Carlo method as a design tool
Table 15.3. Experimentally determined neutron yields for the 1.5 MeV
9
Be(d,n) reaction, determined by scaling thermal and fast neutron doses at
two target positions.
measured and calculated results was also assessed in terms of the effective
neutron yield, as determined by scaling either the fast or thermal neutron
dose rates at both target positions. Since the neutron yield of the 1.5 MeV
9
Be(d,n) reaction is a constant, the effective neutron yield determined by
scaling either the thermal or fast neutron dose, under any initial conditions
(i.e., for any target position), should also be constant.
Table 15.3 shows the effective neutron yield, determined by scaling
thermal and fast neutron doses at two target positions. The neutron yields,
for both dose components and target position, are very consistent and differ
by no more than 13%, with an average value of 1:67 1013 n/min mA. The
uncertainty in the effective neutron yield is estimated as 1:1 1012 n/min
mA, based on the standard deviation of the four values given in table 15.3.
Determining the neutron yield in this manner should be reasonable (assuming
the neutron energy and angular distributions used in the simulation are fairly
accurate), since the level of agreement in both the thermal and fast neutron
doses is good for two target positions that lead to very different neutron spec-
tra in the phantom.
were less than the simulated values, which suggests that from the perspective
of radiation protection, the MCNP calculations are conservative. In most
cases, the measured individual dose components agreed within a factor of
2 with the simulated values. This level of agreement is considered reasonable
for shielding problems and, in particular, for shielding measurements in a
mixed radiation field. The exception to this was the thermal neutron doses
for the shielded phantom, which agreed, on average, within a factor of 4
with simulation.
There are several reasons why the simulated and measured dose values
were only expected to agree within a factor of a few. For the simulation, the
results are only as accurate as the initial source specification. As described
earlier, the neutron and photon sources may not be robust enough to
accurately predict doses far from the target. In addition, the simulation
results have been shown to be sensitive to the surrounding materials in the
radiation vault. While an attempt has been made to accurately describe the
vault in the Monte Carlo model, it is possible that further refinements in
any approximations in the model of the vault may lead to better agreement
between simulation and measurement. Furthermore, the simulation model
neglects all radiation streaming effects, i.e., it assumes that the boronated
polyethylene wall contains no cracks. Streaming effects could have a large
effect on the shielded simulations, since the dose rates have been greatly
reduced by the presence of the shield.
The determination of the photon dose rate is dependent on a correction
for the induced TLD response from thermal neutrons, which can be quite
significant. The thermal neutron response for the TLDs has been experimen-
tally characterized in a thermal neutron beam at the MIT Research Reactor.
The correction is also dependent, however, on the measured thermal neutron
dose determined by the bubble detectors. The measured thermal neutron
doses have been shown to be about a factor of 4 higher than calculated for
the shielded results, and the experimentally determined photon response is
about a factor of 3 lower than calculated. This may result from an over-
correction of the induced thermal neutron response in the TLDs.
and not the cause, which remains unknown. Hence, recurrence is expected
2ÿ5 years later. However, a repeat BNCS procedure, unlike a repeat surgical
synovectomy, is expected to be a simple and straightforward procedure.
Monte Carlo methods have been used extensively in the design and
evaluation of several aspects of BNCS. Monte Carlo calculations indicate,
for a variety of neutron-producing charged particle reactions, that BNCS
treatment times are very reasonable. Times based on a low 10 B concentration
of only 1000 ppm are generally less than 15 min for a knee joint and less than
7 min for a finger joint, assuming an accelerator current of 1 mA. These times
decrease even further when higher synovial boron concentrations are
assumed. Associated whole-body effective doses can be kept low with
appropriate shielding. Whenever possible, Monte Carlo predictions have
been verified experimentally using a combination of detector types. Given
the wide range of neutron energies encountered and the simultaneous
presence of photons, detection and dosimetry is challenging. Combinations
of ion chambers, foil activation, bubble detectors and thermoluminescent
dosimeters have been used in the experimental evaluation of Monte Carlo
predictions for BNCS. Agreement between simulation and experimental
measurement of in-phantom therapy dose-rates was within two standard
deviations.
REFERENCES
[6] Schiff M 1997 Emerging treatments for rheumatoid arthritis Am. J. Med. 102 (suppl
1A) 11S–15S
[7] Binello E 1999 Efficacy of boron neutron capture synovectomy in an animal model.
PhD Thesis. Massachusetts Institute of Technology
[8] Gumpel J M and Roles N C 1975 A controlled trial of intra-articular radiocolloids
versus surgical synovectomy in persistent synovitis Lancet 1 488–9
[9] Davis M A and Chinol M 1989 Radiopharmaceuticals for radiation synovectomy:
evaluation of two yttrium-90 particulate agents J. Nucl. Med. 30 1047–55
[10] Yanch J C, Shortkroff S, Shefer R E, Johnson S, Binello E, Gierga D, Jones A G,
Young G, Vivieros C, Davison A and Sledge C 1999 Boron neutron capture
synovectomy: treatment of rheumatoid arthritis based on the 10B(n,)7Li nuclear
reaction Med. Phys. 26 364–75
[11] Locher G L 1936 Biological effects and therapeutic possibilities of neutrons Am. J.
Roentgenol. 36 1–13
[12] Slatkin D N 1991 A history of boron neutron capture therapy of brain tumors:
postulation of a brain radiation dose tolerance limit Brain 114 1609–29
[13] Farr L, Sweet W H, Robertson J S, Foste C G, Locksley H B, Sutherland D L,
Mendelsohn M L and Stickley E E 1954 Neutron capture therapy with boron in
the treatment of glioblastoma multiforme Am. J. Roentgenol. 71 279–91
[14] Yanch J C, Zhou X-L, Shefer R E and Klinkowstein R E 1992 Accelerator-based
epithermal neutron beam design for neutron capture therapy Med. Phys. 19 709–21
[15] Allen D A, Beynon T D, Green S and James N D 1999 Toward a final design for the
Birmingham boron neutron capture therapy neutron beam Med. Phys. 26 77–86
[16] Chris Wang C-K, Blue T E and Gahbauer R 1988 A neutronic study of an accelerator-
based neutron irradiation facility for boron neutron capture therapy Nucl. Technol. 84
93–107
[17] Bleuel D L, Donahue R J, Ludewigt B A and Vujic J 1998 Designing accelerator-
based epithermal neutron beams for boron neutron capture therapy Med. Phys. 25
1725–34
[18] Klinkowstein R E, Shefer R E, Yanch J C, Howard W B, Song H, Binello E, Black-
burn B W, Daigle J L, Sears S M, Goldie C H and Ledoux R J 1997 Operation of a
high current tandem electrostatic accelerator for boron neutron capture therapy. In
Advances in Neutron Capture Therapy eds B Larson, J Crawford and R Weinreich
(Amsterdam: Elsevier) vol I pp 522–7
[19] Binello E, Shefer R E and Yanch J C 1997 Neutron beam design for boron neutron
capture synovectomy. In Advances in Neutron Capture Therapy (Amsterdam: Elsevier)
vol II pp 459–63
[20] Gierga D P, Yanch J C and Shefer R E 2000 Development and construction of a
neutron beamline for accelerator-based neutron capture synovectomy Med. Phys.
27 203–14
[21] Briesmeister J F 1997 MCNP—A General Monte Carlo N-Particle Transport Code
Version 4B 1997. Los Alamos National Laboratory Report LA-12625-M
[22] Whalen D J, Cardon D A, Uhle J L and Hendricks J S 1991 MCNP: Neutron Bench-
mark Problems Los Alamos National Laboratory Report LA-12212
[23] Whalen D J, Hollowell D E and Hendricks J S 1991 MCNP: Photon Benchmark
Problems Los Alamos National Laboratory Report LA-12196
[24] Gierga D P and Adams K J 1999 Electron/Photon Verification Calculations using
MCNP4B Los Alamos National Laboratory Report LA-13440
346 The Monte Carlo method as a design tool
Summary
Habib Zaidi and George Sgouros
Monte Carlo analysis in nuclear medicine has been used for several decades.
Over this time, its use has evolved from a tool available to only a few expert
physicists to one that may reside on a desktop computer and that is now
potentially available for routine use. This evolution has been driven by
four fundamental and critical developments:
1. The ever-increasing sophistication of imaging instrumentation and
associated image analysis and processing methodologies, an area in
which simulation studies have become essential.
2. The emergence of targeted radionuclides as a therapeutic modality,
requiring improved and patient-specific dosimetry at the macro and
micro level.
3. The availability of three- and four-dimensional computer phantoms
providing flexibility and accurate modelling of populations of patient
anatomies, and attenuation and scatter properties as well as biodistribu-
tions of radiopharmaceuticals in the patients.
4. The near exponential increase in the speed and proportional reduction
in size and cost of computing power.
These four interrelated developments have fostered the work reviewed in this
book.
The recent development of combined modality systems (e.g., SPECT/
PET/CT/MRI) has created new opportunities for Monte Carlo simulation
in imaging. These new imaging modalities in turn make it possible to
meet the new requirements on dosimetry that have been driven by thera-
peutic nuclear medicine. The increasing availability of end-user image
formation/analysis and dosimetry software packages that incorporate
Monte Carlo calculations reflects all four of the developments listed
above.
348
Summary 349
The authors in this book are drawn from the leading research groups around
the world. Habib Zaidi is senior physicist and head of the PET Instrumenta-
tion and Neuroscience Laboratory at Geneva University Hospital and his
research activities centre on modelling nuclear medical imaging systems
using the Monte Carlo method, dosimetry, image correction, reconstruction
and quantification techniques in emission tomography. Dr George Sgouros
is Associate Member, Department of Medical Physics at Memorial Sloan-
Ketting Cancer Center and his research activities are in the field of dosimetry
and systemic radiation therapy physics.
350
Biosketches 351
Katherine S Kolbert, MS
Katherine Kolbert has a master’s degree in computer science from the Courant
Institute at New York University, New York, and is currently with the Depart-
ment of Medical Physics at Memorial Sloan-Kettering Cancer Center in New
York City. She has developed a three-dimensional patient-specific dosimetry
software package as well as other image analysis tools useful to the nuclear
medicine–physics community. Copies of her software packages have been
distributed to researchers around the world. She is also an author and contri-
butor to numerous peer-reviewed papers and abstracts.
email: [email protected]
Web: http://www.ski.edu/
356
Professor in 1981 and full Professor in 1995. Since 1981 he has been a lecturer
at the University of Lund and teaching director at the Department of
Radiation Physics. Since 1997 he has been director of the Jubileum Institute.
From September 1991 until December 1992 he spent a sabbatical year at the
Nuclear Medicine Department, Brookhaven National Laboratory. He has
published about 160 regular articles, four book chapters and holds two
patents. He has served as member of various international committees, such
as the ICRU task group on internal dosimetry, the AAPM Nuclear Medicine
Task Groups on dosimetry, and the EANM task group on radionuclide
therapy dosimetry. He has served as reviewer for peer-written manuscripts
in several journals and organized international symposia and conferences.
Sven-Erik Strand’s research activities are in the field of systemic radiation
therapy physics including: registration of time-sequence three-dimensional
pre-therapy radionuclide images from SPECT/PET with CT/MRI anatomical
information for improved image quantification and to provide the basis for a
three-dimensional radiotherapy planning method, the application of radio-
activity quantification methods to determine tumour and non-target tissue
radionuclide uptake; the use of regional radionuclide pharmacokinetics from
time course quantitative SPECT/PET images, Monte Carlo methods to
produce three-dimensional radiation dose distributions. He also participated
in the EMERALD project and received the Bergman-Loevinger award in
2002 for his research in the field of internal dosimetry.
email: [email protected]
Web: http://www.radfys.lu.se/
361
362 Index
18
F-Fluorodeoxyglucose (FDG), 37 Neutrons, 2, 30, 87, 124
FORTRAN, 5, 139, 234 beam, 30, 43, 324, 327, 328, 329
Fortran-90, 264 Noise, 42, 65, 67, 68, 70, 72, 74
FWHM, 34, 36 Poisson, 18, 203, 204, 205, 219
Non-homogeneous, 70, 71, 274
Gamma rays, 56, 61, 62, 69, 111 Nonlinear, 67, 246, 298
Gaussian, 13, 63 Nonuniform, 142, 189, 209, 212, 220
GEANT, 9, 10, 141
15
Geometrical efficiency, 18 O-water, 78
Optical coupling, 34
Heart, 42, 116, 121, 179 Ordered Subsets Expectation Maximization
Human anatomy, 90, 98, 121, 249 (OSEM), 63, 64, 68
ORNL phantoms, 117, 118, 120, 121
131
I, 36, 38, 43, 45, 68, 96, 97, 98, 117, 119,
161, 163, 164, 165, 166, 186 Parallel computers, 6, 150, 151, 152, 264
111
In, 117, 161, 163, 166 Parallel processing, 31, 150
Interfile, 252 Path length, 21, 159, 167, 205, 206, 251, 252,
Intravascular radiation therapy, 29 290
Iterative reconstruction, 67, 70, 76, 77 PET scanners, 13, 35, 36
ITS, 9, 40, 43, 139, 141 design, 34
Index 363