NCRP Report No. 108 - Conceptual Basis For Calculations of Absorbed-Dose Distributions
NCRP Report No. 108 - Conceptual Basis For Calculations of Absorbed-Dose Distributions
NCRP Report No. 108 - Conceptual Basis For Calculations of Absorbed-Dose Distributions
108
Recommendations of the
NATIONAL COUNCIL ON RADIATION
PROTECTION AND MEASUREMENTS
The idea for this report emerged, in the early 1970's, from the need
of a n NCRP Scientific Committee to characterize the beta-ray depth-
dose distribution in connection with immersion doses. It was realized,
however, that the calculation of such a distribution was only a small
part of the very much larger task concerned with the theoretical,
mathematical and computational concepts involved in the develop-
ment of absorbed-dose distributions in general. To address this issue
in an allencompassing manner, the NCRP formed Scientific Commit-
tee 52 on the Conceptual Basis of Calculations of Dose Distribution.
In either external or internal irradiation, the absorbed dose is
usually non-uniform in any structure and, in particular, in the
human body. This non-uniformity is to be distinguished from the
stochastic variations that exist even in regions where the dose is
uniform and that are the subject of microdosimetry and not this
report. Many illustrations of absorbed dose non-uniformity come to
mind: for example, the absorbed-dose distributions from hot parti-
cles, from internal emitters, from radiation therapy, from radiation
accidents and from environmental radiation. There can even be addi-
tional non-uniformity with respect to time of the non-uniform distri-
bution, for example, in the redistributions of administered radioac-
tivity in the body.
For all absorbed-dose calculations, there is a source (or sources) of
radiation and a receptor (or receptors) of some of the energy of this
radiation, with or without intervening material between the source
and receptor. The calculation of absorbed-dose distributions requires
specification of the sources and receptors, characterization of their
geometrical relationships and consideration of the physical interac-
tions of the radiations involving attenuation, scattering and the
production of secondary radiations. All these processes are consid-
ered in the basic transport equation, the general theorems and prop-
erties and the methods of solution of which are described in the
transport theory.
The report is a systematic presentation, discussion and compila-
tion of all the concepts involved. It contains some complicated mathe-
matics that will be of interest to the mathematically knowledgeable,
but that should not discourage those not mathematically inclined.
iv / PREFACE
The text of the report contains detailed explanations of all the con-
cepts and of the consequences of the equations so that, even omitting
the mathematics, a broad and comprehensive understanding can be
obtained of what is entailed in the calculation of an absorbed-dose
distribution.
The cutoff date for the report is about two years ago and, hence,
the report is lacking in the most current references. However, this
field does not evolve at a rapid pace and the current literature is,
therefore, not abundant and can be reviewed easily.
In accord with the recommendations of NCRP Report No. 82, SI
Units in Radiation Protection and Measurements, as of January 1990,
only SI units are used in the text. Readers needing factors for conver-
sion of SI to conventional units are encouraged to consult Report No.
82.
This report was prepared by NCRP Scientific Committee 52 on
Conceptual Basis of Calculations of Dose Distributions. Serving on
the Committee during the preparation of this report were:
Harald H. Rossi, Chairman
105 Larchdale Avenue
Upper Nyack, New York
R. G. Alsmiller, Jr. William C. Roesch
Engineering Physics and 1646 Butternut
Mathematics Division Richland, Washingtan
Oak Ridge National Laboratory
Oak Ridge, Bnnessee
Martin J. Berger Lewis V. Spencer
5011 Elm Street Post Office Box 87
Bethesda, Maryland Hopkinsville, Kentucky
Albrecht M. Kellerer Marco A. Zaider
GSF Institut fiir Strahlungbiologie Radiological Research
Neuherberg, Germany Laboratory
Columbia University
New York, New York
NCRP Secretariat: Thomas Fearon (1976-80)
J a m e s A. Spahn, Jr. (1981)
J a m e s T. Walker (1982-84)
Constantine J. Maletskos (1985-91)
The Council wishes to express its appreciation to the Committee
members for the time and effort devoted to the preparation of this
report.
PREFACE, 1 V
.
1 Introduction .................................... 1
1.1 The Concept of Absorbed Dose . . . . . . . . . . . . . . . . . . 1
1.2 Dose Measurement and Dose Calculation ......... 3
1.3 Elements of Dose Calculations . . . . . . . . . . . . . . . . . . 4
.
2 Transport Formalisms . . . . . . . . . . . . . . . . . . . . . . . . . . 6
2.1 Concepts in Dose Calculations . . . . . . . . . . . . . . . . . . 6
2.2 Transport Equation ........................... 8
.
3 Sources . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 14
3.1 Specification of Sources ......................... 14
3.2 Simplified Representations of Sources . . . . . . . . . . . . 15
.
4 Receptors ..................................... 17
.
5 Cross Sections . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 20
5.1 Schematization . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 20
5.2 General Aspects of Required Cross Sections . . . . . . . 22
.
6 Transport Theory-General Theorems and
Properties ..................................... 26
6.1 Integral Form of the Transport Equation . . . . . . . . . . 26
6.2 Iterative Solutions (Orders of Scattering) . . . . . . . . . 27
6.3 Density Scaling Theorem . . . . . . . . . . . . . . . . . . . . . . 28
6.4 Fano's Theorem . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 29
6.5 Energy Conservation . . . . . . . . . . . . . . . . . . . . . . . . . 29
6.6 Superposition . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 30
6.7 Adjoint Transport Equation .................... 31
6.8 Reciprocity . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 33
6.9 Transport Equations in Commonly Used Coordinate
Systems .................................... 34
.
7 Transport Theory-Methods of Solution .......... 36
7.1 Introduction ................................. 36
7.2 Radiation Equilibrium and Space-Integrated
Radiation Fields . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 38
7.3 Continuous Slowing-Down Approximation (CSDA) . . 40
7.4 Numerical Integration Over Energy ............. 43
7.5 Elementary Problems Involving Particle Direction . . 45 '
where Ei, (E,,,) is the sum of energies of all the charged and
uncharged ionizing particles that enter (leave)the volume, excluding
rest mass energies, and XQ is the algebraic sum of all changes
(decreases: positive sign; increases: negative sign) of rest-mass
energy in mass-energy transformations occumng in the volume.
A few clarifications of Equation (1.1) are necessary a t this point.
The energy imparted results from random discrete energy deposition
events by individual ionizing particles andlor their secondaries. The
quantity E is therefore stochastic in nature and governed by a (nor-
malized) probability distribution function fv(d where V is the volume
containing m. The mean value of e,
2Becauseof the discrete manner in which energy is imparted, the limiting process
V- 0 is obviouely an idealization. At all times the volume V should contain a large
number of atoms and molecules.
1.2 DOSE MEASUREMENT AND DOSE CALCULATION / 3
It follows (ICRU, Report 33) that
DP) = lim Z
v+ 0
FEV,
where the limiting process is such that ? is contained in the volume
v.
The stochastic variations of specific energy that occur even in a
medium that is uniformly irradiated must not be confused with
variations of absorbed dose in a medium that is not uniformly irradi-
ated. The latter often occurs with internally deposited radionuclides,
especially with alpha particles. For example, if 2S9pUparticulates,
having a diameter of the order of a micrometer, are lodged in tissue,
substantial variations in absorbed dose occur over comparable
dimensions.
The theoretical and experimental study of the distributions in z is
the objective of microdosimetry (see ICRU, 1983) and is beyond the
scope of the present report.
No formal definitions of other radiation quantities are given here
because they can be found in many publications and because their
meaning is usually obvious from the context. The units employed
are those of the International System (SI) except that, in accord with
common practice, the electron-volt (eV) is frequently used as a unit
of energy.
which, when integrated, yields the flux density (also called fluence
mte)
6,t) = J#bu &i,.E,2,t). (2.8)
'All the quantities used in the discussion of the transport equation refer to a single
type of particle as opposed to the more general definitions given in the previous
section. For simplicity the same notation is used.
2.2 THE TRANSPORT EQUATION 1 9
One should note that s includes here not only "true" sources, but
also contributions from other types of radiation in the field undergo-
ing transformation to the particle type to which Equation (2.14)
(3mll
refers. This is discussed later in more detail. The term - repre-
sents changes in n (per unit time) due to collisions in the medium.
Equation (2.14) is simply the balance of the particles entering or
leaving the phase-space element dEdu throughout the volume V, per
unit time.
With:
This is the classical form of the linear transport equation with the
interaction probabilities explicitly shown.
It is important to remark also that Equation (2.24)is linear because
it involves C#I only to the first power (for instance, the other well
known transport equation, the Boltzmann equation for gases,
belongs to a nonlinear transport formalism; see also Subsection 6.1).
Examples
In order to further elucidate the concepts introduced above, several
examples of practical application of the transport equation, Equation
(2.24), are given here. These examples also allow the introduction of
a number of additional concepts widely used in dose calculations.
The presentation here follows closely that of Rossi and Roesch (1962)
and Roesch (1968).
Possibly the most typical dose calculation involves the situation
where a constant radiation source is used to irradiate over a given
time interval (O,t,). One is interested in the total dose delivered.
2.2 THE TRANSPORT EQUATION 1 11
The transport equation, Equation (2.241, can be simplified by
observing that
and by defining:
where Sois the expected number of particles produced per unit vol-
ume during [O,t,] in dEdu about E and ii. With these simplifications
and definitions, Equation (2.24), on integration over time, becomes
For a field with only one type of charged particle and Q = O (no
secondary particle transport), one obtains the dose delivered during
[O,tol by inserting Equation (2.28) in the time integral of Equation
(2.13):
12 / 2. TRANSPORT FORMALISMS
pk(Eu)=
'I
- dE: p, (E,+E:) [E, -E: - Q,(E,+E:)I.
Eu
(2.35)
The first integral in Equation (2.34), with the correction for rest
mass energy changes, q,, represents, then, the sum of the initial
kinetic energies of all the charged particles created by the uncharged
particles per unit mass:
1
KG) = IdE,E u ' P u ~E,)
, - pk(Eu). (2.36)
P
K is called the kermu of the uncharged particles.
Frequently the second integral in Equation (2.34) is negligibly
small compared to K (the same holds, then, for dQJdm). The condition
when ly,(?, Ec)does not depend on F' (i.e., is constant over the field)
is called charged particle equilibrium. Then, because $'PC= 0, one
has
5.1 Schematization
The probability that a free flight will be terminated when its length
is between C and t + d t is given by
where I is the unit operator and, for reasons apparent below, the
notation So= s is introduced.
By introducing the operator (I - Tk)-I one can rearrange Equation
(6.8') as
28 1 6. TRANSPORT THEORY-THEOREMS AND PROPERTIES
Then,
etc.
The formalism described by Equations (6.9) and (6.10), known as
the Neumann expansion, has a simple interpretation in terms of
sequences of collision events. Thus, the first term, do, represents
particles arriving a t ;:directly from the source [see Equation (6.6)l
and after attenuation. The second term, dl, represents particles that
had a collision and scattered once to produce a new "source", S1=
Kdo, whose output was transported, TS,, to the point ,; and so on,
hence the name "orders of scattering." It is clear from Equations
(6.10) that the Neumann expansion, Equation (6.91, expresses an
iterative solution to the transport equation.
~ ( p=
) 1
0
P
dp1p(7-prz), (6.11)
then the transport Equations (6.8 and 6.8') will depend only on
local density in the final (source) term. But if So is also everywhere
proportional to the local density, Equations (6.8 and 6.8') can be
written entirely in terms of position-independent functions, and, as
a consequence, must be independent of position.
This theorem can be applied in plane symmetry to uniform radia-
tion sources confined to a plane. There results a plane-scaling theo-
rem which applies to configurations with differing plane-density
variations (Spencer et al., 1980).
6.6 Superposition
where f and g are arbitrary functions and the inner product may be
defined as
~t = -
-.
.v + CI.r
+
II
- d ~ '~U~~,(E-+E:Z+S~~),
for E'<E.
(6.19)
where the notation shows that 1+9is a solution of the transport equa-
tion corresponding to the point source, Equation (6.23). But $' of
Equation (6.25) is the dose at the point ;f' in the receptor due to the
point source in Equation (2.31).
6.8 RECIPROCITY 1 33
It can be shown, as might be intuitively obvious, that for a single
radiation with an extended source, s,(;f,
E,;), the dose is given by
6.8 Reciprocity
az
In expressions such as Equation (6.5), the generalized variable can
be written as follows: ?+z, &uZ, 7 -pii+z-pu,.
Almost as important is the transport equation in spherical symme-
try. The natural reference axis for this case is the radial unit vector,
- 6.9 TRANSPORT EQUATIONS 1 35
n. For spherical sources, 4 varies only along a radial line: All varia-
tions of 4 from one point to another on any of the concentric spheres
vanish. Hence 4=4(r, y,E,t), where r = 3 .Z and u, = u - n = cos
+ +
az
84 vanishes
If there is no dependence of the source on f?,the term with -
af?
and
7.1 Introduction
and (7.2)
where the ai are known and depend on the algorithm used for evalu-
ating the integral. In this way, the transport equation is replaced
with a set of algebraic equations for the unknowns di.
In the second method, one approximates the function &x) using
a finite set of known functions f,(x):
7.1 INTRODUCTION / 37
where the unknowns now are the expansion coefficients 4i. As an
example, it is customary to represent the angular dependence of the
flux density using spherical harmonics, Ye,(;):
noting that
where &(E) = 1 En
E,+1
dEfpS(E'+E),with the integration over E' start-
ing only a t E for E > En+,.
The scattering integral is replaced by a
sum, so that Equation (7.29) takes the simple and soluble form
(E, is the maximum energy of neutrons emitted):
n-1
4-4 / 7. TRANSPORT THEORY-METHODS OF SOLUTION
and
with
and
7.6 PENETRATION STUDIES 1 47
MoliBre's procedure involves writing the exponent in the scaled form
1
zAu21n(u2/b)= - x2 + -x21nx2.
B (7.39)
The two terms can be roughly identified as the Gaussian peak and
the largedeflection tail. One immediately finds that x2 = zABu2,
where B is the solution of
B - lnB = ln(bzA) (7.40)
The term exp [ + B-'x21nx21is expanded into powers of the exponent
which converge rapidly when B >> 1.There results a series of which
the first term is a Gaussian in s = e l m , and the next two terms,
which can be expressed a s functions fT1)(s),F2'(s) (tabulated in
Molihre, 1948), give additive corrections that evaluate the large-
deflection tail accurately.
7.7.2 Electrons
Curve 1 in Figure 7.3 is an equilibrium spectrum for electrons
calculated by Spencer and Fano (1954); curve 2 is based on the
Fig. 7.2. Electron energy degradation spectra in water. Scaled collision density,
nc(T,TJ (TIT,), as a function of the electron energy, T , for various initial electron
energies, To.The usual track length distribution, y(T,TJ, is equal to A(n nc(T,TJ,
where A (T) is the mean free path between inelastic collisions. The collision density
should also include a contribution from the first collisions made by previously unscat-
tered electrons; this contribution, represented by a delta function 6(T- T J , is not
shown. Also not shown is the fine structure of n,(T,TJ near T,. (After Berger, 1984.)
7.7 SPECTRAL EQUILIBRIUM AND RELATED CONCEPTS 1 53
2.04 MeV
I
DISTANCE, x, (arbitrary)
Fig. 7.7. The field due to a one-dimensional oscillatory source-strength distribu-
tion as a function of distance (Equation 7.45). (---)the resulting fluence density when
the factor of the oscillatory term of Equation (7.45) is 0.9; (. . . .) when it is 0.1.
8.1 Principles
and
j=l
o,. Imagine on a straight line, adjacent
segments of length u,lCuj, odCuj, etc., forming a total segment of
length 1. A random number, 5, uniformly distributed on (0,l)will
then be contained in one of these segments and thus select a value
xi. This may be done by using an algorithm in which
Next, let
where
64 1 8. MONTE-CARLOMETHODS
A value x can be chosen as follows:
a) Select i with probability proportional to cri [Equation (8.811;
b) Select x from ((x) [Equation (8.7)l;
C) Accept x with probability gi(x) [Equation (8.811;
d) If x is rejected go back to a).
It can be shown (Messel and Crawford, 1970) that the acceptance
probability for this method is I/&;. Very ingenious decompositions
for particular functions f(x) are available, such that the maximum
sampling efficiency is attained.
Consider, as an illustration of this procedure, the Klein-Nishina
cross section for the Compton scattering of a photon with energy k
to energy k' (the energies are in units of the electron rest mass).
From Equations (A.3) and (A.6) of Appendix A one has
where
c = k'lk
and
where
and
where x is the distance from the point to the volume element dVA.
The average dose rate in the receptor region B due to radiation
from the source region A is, therefore,
Let h(x) be a function of the distance x; then, one can reduce the
6-dimensional integral over the regions A and B to a one-dimensional
integral
where a is the total activity in A, and the right hand expression must
be used if A has zero volume.
A point-pair distance distribution usually cannot be determined if
one deals with unbounded regions. It is then practical to use a modi-
fied quantity, called the geometric reduction factor (Berger, 197017,
U&) = VBpm(x)/4~x2, VB 0. + (9.6)
If B has zero volume, U,(x) is not defined; but one can then utilize
One can show that the geometric reduction factor, UAB(x),is equal
to the probability that a random point of A ends up in B, if shifted
in a random direction by distance x; a more detailed discussion of
the point-pair distance distribution and the geometrical reduction
factor is given in Section 9.4.
With the geometrical reduction factor one can rewrite Equation
(9.5) as
(-
- -
a
vB /e(x) u-(x) b, for VB f 0,VA finite
1
-
(9.8)
Dm = a I d x ) UBA(x)&, for other situations,
but with VA > 0.
'In the case of an extended region A, the term average geometric reduction factor
has formerly been employed. The proximity function, s,(x) = VB . pAB(x),
can also be
used, and is closely related to U,(x).
9.2 RECIPROCITY THEOREM 1 77
If the receptor region B is unbounded, the mean dose rate may not
be meaningful; one may then consider the rate of energy imparted
to B as
1m ( = a 1e(x)um(x)61,
= a VB le(x) UB*(x)&,
for V,>O, VAfinite
where MAand MBare the mass of A and B, and where the dose rate,
DGA,2,), from a unit activity a t point 2, to point 2, depends on the
composition and density of the medium throughout the vicinity of 2,
and f B that is relevant to the radiation transport Gom f Ato 2,.
For the reciprocity theorem to hold, DGA, must equal DGB,
zA),i.e., an activity a t point 2, must produce the same dose at point
zB,as the same activity a t 2, produces a t 2,. For scattered radiation,
this identity will commonly not apply. However, for the unscattered
primary radiation, it can usually be assumed. It applies strictly for
the unscattered radiation if all interaction coefficients are propor-
tional to density; in this latter case DG,, 2,) is equal to a function
h(s)/x2,where h(s) depends only on the type of radiation and on the
integral, s, of the density along the straight line between 2, and zB.
9.3 Isotropic Point-Source Kernels
I I
I
10 MeV
0 20 40 60 80 100 120
DISTANCE, pm
Fig. 9.1. Integral, isotropic point-source kernels, normalized to the initial energy,
E,, for cr particles in water. The ordinate is the fraction of the emitted energy that is
deposited along the particle track, while the abscissa is the distance that the particle
has traveled. The electronic stopping power is taken from Ziegler (1980).
80 1 9. GEOMETRIC CONSLDERATIONS
Figure 9.2 gives the resulting integral, isotropic point-source ker-
nels normalized to the mean initial electron energy for beta rays and
monoenergetic electrons of various energies. The point kernels are
given in a scaled form, i.e., as a function of the variable XI%, where
x is the distance from the source and xS0 is the radius of the sphere
about the source within which 90 per cent of the emitted energy is
deposited. One notices that the scaled, integral point kernels are not
very sensitive to the energies of the beta rays and electrons.
High energy photons require a consideration of the scattered pho-
ton contribution and, depending on the nature of the problem, also
a consideration of energy transport by charged secondaries. Figure
o d.1 d.2 013 014 015 0:s d.7 0:s 019 1.0
RELATIVE DISTANCE, xlx,,
Fig. 9.2. integral, isotropic point-source kernels, normalized to the mean initial
electron energy for beta rays and monoenergetic electrons. The ordinate is the fraction
of the emitted energy that is deposited within a sphere of radius x about the point
source. The abscissa is the relative distance, d%, where % is the radius of the sphere
within which 90 per cent of the energy is deposited. The resulta are for interactions
with water, and x, distances are shown in the tabulation i n the figure. The curves
were developed from data in Berger (1971).
9.4 POINT-PAIR DISTANCE DISTRIBUTIONS 1 81
9.3 gives a numerical example for a photon, isotropic point-source
kernel. Neutrons present a similar problem where the consideration
of the scattered radiation is essential.
50-
0.01
PHOTON ENERGY, MeV
Fig. 9.3. Energy deposition around monoenergetic, point-isotropic sources of pho-
tons with energies up to 2 MeV in an unbounded water medium. The figure shows the
radii of spheres around the point source within which the indicated fractions of the
emitted photon energies are absorbed. Several different plots are possible; in this case,
the radii of the spheres are plotted as function8 of the source photon energy. The
curves were developed from data in Berger (1968).
82 9. GEOMETRIC CONSIDERATIONS
where V and S are volume and surface of the body, and 4VlS is its
mean chord length. Furthermore, for a convex body, U(x) is the
sum distribution of the lengths of rays in random directions from a
random point in A to the surface of A.
For both the autologous and the heterologous cases, the two quanti-
ties p,(x) and U,(x) are equivalent. The definition of the point-pair
distance distribution pAB(x)as the probability density of distances
between pairs of random points leads to solutions for the spheroid
and for the cylinder [Equations (C.7 and C.lO) in Appendix Cl. The
separate definition of UAB(x)in Subsection 9.1 can be formulated in
a way that is also suitable for computations: If A(x,8) is the translate
of A by a distance x in the direction 6, then U,(x) is numerically
equal to the volume of the intersection A(x,B)nB uniformly averaged
over all directions and divided by VA(for related considerations see
Enns and Ehlers, 1978). This form of the definition leads to the
solution for the sphere or for two concentric spheres A and B [see
Equations (C.3 and C.15) in Appendix Cl. The relation
US(x) = (VAUA(x)+ VB UB(x) - 2 VAUAB(x))1 (VA-VB) (9.14)
9.4 POINT-PAIR DISTANCE DISTRIBUTIONS 1 83
is then used to obtain the geometric reduction factor for the spherical
shell S = A - B [see Equation (C.4)I.
Finally, it may be noted that the geometric reduction factor can
also be defined in two-dimensional space. In this case, it can be
calculated from the point-pair distance distribution by the relation
-:1
Hf?) = D(L,,F)Q(LJdL, (10.1)
where L,, and L,, are the minimum and maximum values of L,
l o . CALCULATION OF THE DOSE EQUIVALENT / 85
for which D(L,,7) has a nonzero value. Equal values ofH are intended
to correspond to equal probabilities of deleterious effects regardless
of the type of charged particles that are producing the absorbed doses.
The recommended quality factors as a hnction of L, are given in
Table (10.1) (ICRU, 1980). In the last several years, changes in the
value of the neutron quality factor have been recommended (ICRP,
1985; ICRU, 1986; NCRP, 1987) and simplified relationships are
being considered between Q and LET (ICRP, 1991). However, the
values used here serve to illustrate the principles of the calculations
and of the examples given later.
TABLE10.1-Quality factor, Q as a fumtwn of Lma
keV/prn
3.5 or less 1
7.0 2
23. 5
53. 10
175. or greater 20
"ICRU (1970);ICRP (1977).
where
4fE,F) = the fluence per unit energy of uncharged particles a t
energy E and position F,
Em,,, Emin = the maximum and minimum energy a t which &E,?)
has a nonzero value,
c = summation over all charged particle species,
List of Symbols
- - - - - - - - - . I
(3) Nucleon-antinucleon
pair ~roduction
Fig. A.1. Locus of equal atomic cross section for Compton and photoelectric inter-
actions (u = T ) , and for Compton and pair-production interactions (u = K ) . The
incident photon energy is hv,, and Z is the atomic number of the atoms in the absorber.
Compton collisions have larger cross sections than any other mode of interaction in
the entire domain of medium-energy photons marked "Compton effect dominant"
(From Evans, 1958.)
where k is the energy of the incident photon and B the binding energy
of the atomic shell from which the electron is ejected. At a given
energy, the stronger the binding, the larger is the photoelectric cross
section. At photon energies greater than the K-shell binding energy,
approximately 80 percent of the photoelectrons are ejected from the
K shell. The photoelectric cross section (per atom) is proportional to
- -
Zmlkn,where m 4 and n 3. At low photon energies the dependence
of the photoelectric cross section on the photon energy is not mono-
tonic but exhibits sharp discontinuities a t the K -, L - . . .absorption
edges (see Figure A.2), that is, a t energies a t which the photon has
just enough energy to eject a n electron from a particular shell.
At low photon energies, the photoelectrons are ejected mainly a t
right angles with respect to the direction of the incident photon. At
high energies, the peak in the angular distribution is in the forward
direction; half of the photoelectrons are emitted within a cone whose
half-angle (with respect to the photon direction) is approximately
equal to 20" a t 1 MeV, 30" a t 0.5 MeV, and 60" a t 100 keV.
One of the most comprehensive sets of accurate photoelectric cross
sections is that calculated by Scofield (1973) based on the numerical
solution of the Dirac equation and the numerical evaluation of the
pertinent matrix elements. Excellent reviews of the theoretical and
experimental aspects of the photoelectric process have been given by
CARBON 12 ' 6 )
1
I I I 1 I I I I I
>\?< -
LEAD lZ*821
?.Ft0 -
-
-
-
cn
cn
0
-
IT:
Fig. A.2. Photon cross sections for carbon and lead. (From Hubbell et al.,
A.l PHOTON CROSS SECTIONS 1 97
Fano and Cooper (1968) and Starace (1982) with emphasis on the
low-energy region, and by Pratt, et al. (1973) covering the energy
region above 10 keV.
0 20 40 60 80 100
ATOMIC NUMBER, Z
Fig. A.3. Fluorescence yield for the K shell as a function of atomic number. ( A h
Bambynek et al., 1972.)
x {l +
y2(1 - ~ 0 ~ 0 ) ~
(1 + cos2e) [ l + y (1- c0sO)l
sine do,]
where y = klmc2 and re = e2/mc2 = 2.82 x 10-lScm is the classical
electron radius. Graphs of the Compton energy-angle relation and of
the differential and integrated Klein-Nishina cross section can be
found in a report by Nelms (1953). A thorough review of the theoreti-
cal and experimental aspects of Compton scattering was given by
Evans (1958).Figures A.4, A.5 and A.6, taken from Evans, illustrate
the number us angle distributions of scattered photons and recoil
electrons, and the number us energy distribution of the recoil elec-
trons.
Because the target electrons are neither free nor at rest, the Comp-
ton relation between energy loss and deflection angle does not hold
precisely. Photons incident with energy k and deflected through an
angle can have various energies k" after the scattering that are
distributed around the value k' given by the Compton relation. This
is illustrated in Figure A.7 with experimental and theoretical results
of Schumacher (1971) for 662-keV gamma rays scattered through
135" in lead.
The effect of the binding of the atomic electrons is to reduce the
scattering cross sections, particularly at small angles. This reduction
can be calculated, in good approximation, with the use of an incoher-
ent scattering function S(q,Z). The scattering cross section with
binding effect included, duBD,is calculated from the expression
dam = duKN S(q,Z), (A.7)
A.1 PHOTON CROSS SECTIONS 1 101
where
,- 1 Experiment of D. V. Rao
-
where T-and T+are the kinetic energies of the electron and positron
and 2mc2 is the combined rest energy of the pair. The threshold
energy for the process is k = 2mc2 = 1.022 MeV. The energy is
shared approximately randomly between the electron and positron.
Figure A.9 shows the scaled cross section as a function of the fraction
of the total kinetic energy of the pair electrons that is acquired by
the positron, calculated according to the theory of Bethe and Heitler
(1934). For high incident photon energies, the pair emerges predomi-
nantly in the forward direction. The average angle between the
direction of the photon and the directions of the electron and positron
has the order of magnitude mc2/T'.
Pair production also takes place in the electric field of the bound
atomic electrons. Recoil momentum and energy are transferred to a
single electron which is thereby ejected from the atom. In a track
visualization device (such as a photographic emulsion or a bubble
chamber), three charged particles are observed: the electron-positron
pair and the ejected atomic electron. For this reason the process is
also called triplet production. The threshold for triplet production is
k = 4 d = 2.044 MeV. The cross section for pair production in the
field of the Z atomic electrons, K,, is equal to (qJZ)K, where K , is the
cross section for pair production in the field of the atomic nuclaus.
The proportionality constant % is a quantity of order unity. As shown
in Table A.3, q,/Z actually goes to zero near the threshold and reaches
values of 1.1 to 1.2 a t very high energies.
A positron created in a pair-production event eventually recorn-
bines with an electron. The energy made available by the disappear-
A.l PHOTON CROSS SECTIONS 1 105
where
A.l PHOTON CROSS SECTIONS 1 107
is the momentum transfer. F(q,Z) represents the probability that
the recoil momentum q is transferred from the photon to the Z
electrons of the atom without any resulting excitation or ionization.
Tables of form factors derived from Hartree-Fock wave functions
have been given by Hubbell et al. (1975).
A review of the theory of Rayleigh scattering, including a treat-
ment more rigorous than the use of form factors, can be found in a
paper by Kissel et al. (1980) which also contains tabular results
for elastic scattering by lead. The elastic scattering cross section
(differential in deflection angle) for atomic numbers Z = 2 to 26 and
for photon energies from 0.5 to 200 keV, calculated with the use of
form factors, has been presented in graphical form by Brown (1975).
Angular distributions for scattering by carbon, takenfrom that refer-
ence, are shown in Figure A.lO.
The Rayleigh cross section peaks in the forward direction and is
large only a t energies where the photoelectric cross section is even
larger. Therefore, elastic scattering (involving no energy loss and
only a small deflection, and dominated by photoelectric absorption)
is often omitted from transport calculations.
Fig. A.lO. Calculated, normalized differential cross sections for coherent (Rayle-
igh) scattering of photons for various photon energies. Upper figure: energy range 0.5
to 10 keV. Lower figure:energy range 15 to 200 keV. To obtain absolute cross sectiona
in units of cmYsr, multiply values from figuresby Z2r,2, where rz = 7.94 x m2.
(After Brown, 1975.)
A . l PHOTON CROSS SECTIONS 1 109
I
r h h 7 energy S&nirg
f i r pducf i n \
\ \
~lod/&c A Annihi/ofin
+
o&olpIion
ff-ce
~1ofion
I I
t 1 1
1 - 1
Ionizafion and &itofion
Riohe/ecfhs,Auger e l e c f m ~ ~ ~
h okqh
Nuclear
Fig. A.12. Flow diagram of conversion of energy from one form to another in the
course of photon energy absorption in a medium. (From Hubbell and Berger, 1968.)
/---
/' 10-2
10-2 lo-! I 10 102
the older compilations are usually not completely obsolete, and they
also retain explanatory text of lasting value. The most up-to-date
photon cross-section data above 10 keV for dosimetry calculations
can be obtained by combining the data given in Hubbell (19771,
Hubbell et al. (1980), and Hubbell (1982). A complete compilation of
photoabsorption and coherent scattering data a t low energies (30 eV
to 10 keV), based on experimental data, has been given by Henke et
al. (1982).
N
E
0
0.8
s
N
2-
0
F
Z:
0.4
(I)
(I)
0
u
0
1.0 MoV
. -
.r
6 1.I
MeV 2 MeV
Fig. A.17. Scaled differential cross section for bremsstrahlung production in the
field of the atomic nucleus as a function of the ratio of the photon energy, k, to the
kinetic energy of the incident electron, T, in carbon and gold. The scaling factor is
described in the text. The curves for 50 keV have been plotted from the tables of Pratt
et al. (1977). The result. for 50 MeV have been calculated according to the theory of
Davies et al. (1954) and Olsen (1955).
T.MeV T,MeV
Fig. A.18. Scaled radiative energy-low cross section, q5&,, for bremsstrahlung in
the field of the atomic nucleus as a function of electron energy, T. Points below 2 MeV
are from the calculationsof Ratt et al. (1977), and pointa above 50 MeV are from the
highenergy theories of Davies et al. (1954)and Oleen (1955).Curves are from a least-
squares fit to the theoretical points. ( f i r ICRU. 1984.)
Fig. A.19. Ratio of the scaled cross sections 4 & , and that represent the
radiative energy losses resulting from the emission of bremsstrahlung in the field of
the atomic electrons and in the field of the atomic nucleus, respectively, as a function
of electron energy, T. The total radiative stopping power is proportional to p k , +
24+
,. (ARer ICRU, 1984.)
where S(T) = -dT/ds is the energy loss per unit path length (stop-
ping power). Suppose that the particle had an interaction of type j
when it had energy To.The probability pj(T)dTthat the next interac-
tion of this type will occur when the particle energy is between T
and T + dT is then given by
dT'
pj(T)dT = exp dT (A.34)
22;;
20 -
t @
4
4
GASES.
GASES,
SOLIDS.
EXPERIMENTAL
INTERPOLATED
EXPERIMENTAL
P -SOLIDS. INTERPOLATED
18-
0 10 20 30 40 50 60 70 80 90 00
ATOMIC NUMBER, Z
Fig. A.20. Dependence of the mean excitation energy, I, for elemental substances
on the atomic number 2.I is roughly proportional to 2.To indicate the residual, non-
smooth, dependence on Z in greater detail, the ratio IlZ is plotted as a function of Z.
(From ICRU, 1984.)
Fig. A.21. Electron stopping power in water. Dashed curve: mass collision stop-
ping power; dotted curve: mass radiative stopping power; solid curve: total mass
stopping power (collision plus radiative).
Figure (A.21) also gives the mass radiative stopping power ofwater
for electrons, calculated f r ~ m
the expression
PROTONS
a d r a , a'
a I Lev
ELECTRONS
Fig. A.23. Mass collision stopping power divided by z2 for electrons, protons, and
alpha particles in water vapor as functions of particle velocity expressed as the ratio
of the particle velocity, o, to the Bohr velocity, v, = ~1137.Particle-energy scales are
shown below the figure. ( f i r Berger, 1978.)
The heavier the ion, the greater is the reduction of the stopping
power due to the replacement of the nominal by the effective charge.
This is illustrated in Figure A.24 which gives the collision stopping
powers (divided by z2) of various ions in aluminum, plotted as func-
tions of the kinetic energy per atomic mass unit. The determination
or prediction of the effective charge is a difficult task. The extensive
literature on this subject has been reviewed by Betz (1972), who
suggests an empirical formula
z*
-= 1 - C exp [ - u/(u,,z~)l, (A.40)
Z
138 / APPENDIXA
which has been given by Wilson et al. (1977). This formula is a fit to
numerical results obtained with a krypton-carbon potential, and is
stated to give results in good agreement with experimental stopping-
power data for many other ion-atom combinations.
where Enis the neutron kinetic energy, and A is the ratio of the mass
of the target nucleus to the mass of the neutron. The largest possible
energy transfer is
a s ecM
cos 8 =
(1 + 2A cos OCM +A2)liz'
and is always smaller than the center-of-mass system deflection
angle Om. For hydrogen, 6 = ocM/2
The angular deflection in the center-of-mass system, at energies
away from resonances, is almost isotropic until the neutron energy
is greater than -10 A-213MeV. For isotropic scattering,
- 2
cos eCM= 0, and cos 0 = -
3A '
(A.47)
0
18Z
r
,a'
. . .
ea
. #
ia'
d .
id
. t
10'
. 1
e4 > eS ,a7
I
ENERGY, eV
Fig. A.25. Angular dependence of the cross section for the elastic scattering of
neutrons from 160.Figure A: angular distribution (normalized to unity) as a function
of the cosine of the scattering angle in the center-of-mass system for various neutron
energies. Figure B: average cosine of the scattering a n d e in the center-of-masssystem
as a function of the neutron energy. (After Cullen and Hlavac, 1972.)
-- E*
1 + A,'
where A, is the ratio of the mass of the recoiling (excited) nucleus to
the neutron mass. The de-excitation of the excited nucleus takes
place mainly via the emission of gamma rays, and the knowledge of
this gamma-ray spectrum is an essential part of the neutron cross-
section information for transport calculations.
The neutron interactions, 2) to 6), listed above are characterized
by sharp thresholds and resonances. By way of example, experimen-
tal cross sections for some of these interactions, with 160as target,
are shown in Figure A.26 as functions of the neutron energy. Figure
A.27 shows the total cross sections for scattering from hydrogen,
carbon, nitrogen and oxygen and from a 4-component tissue-equiva-
142 / APPENDIXA
I 1
," 80. -
i -
!2
ELASTIC
--
--_- - - --__
-
YI. '
L . .A
ma ID, D. 10' 1C 10' 10'
ENERGY. eV
I@' THRESMOLO
2.W MeV
10.
-.X
THF€SHOLD
8 lo' MH
@
,.
E
,,
Z Y
B8 10'
I@' - THREWOLO
I O M MeV
10' -
5 8 7 8 9 1 0
ENERGY, MeV
Fig. A26. Illustration of cross sections for various neutron interactions with '
'
0
as a function of neutron energy. Note change in energy scales between upper and
lower figures. (ARer Cullen and Hlavac, 1972.)
A.3 NEUTRON CROSS SECTIONS 1 143
a particles
a
LL
-
Y
BEHROOZ
I I I I I I
10 20 30 40 50 60 70
NEUTRON ENERGY. MeV
10
- I
16
I ! I I I --
-
8- 0 -
i
6 -
DIMBYLOW
-
-
-
I I I
10 20 30 40 50 60 70
NEUTRON ENERGY, MeV
Fig. A B . Kerma factors as a function of neutron energy for carbon and oxygen.
In comparison to Figure A.28, this figure extends to higher neutron energies and
includes the results of several investigators. ( k o m Cross, 1986.)
150 / APPENDIX A
energy loss, and also undergo occasional nuclear interactions that
can result in scattering or absorption. In the continuous slowing-
down approximation, the probability, p(T)dT, that a pion with initial
kinetic energy Towill undergo its first nuclear interaction when its
energy is between T and T + dT is given by
I I
lb 2b 3b LO Sb 60 70 80 90 100 110 120
ENERGY, MeV
Fig. A.30. Total and partial m a s sections for the interaction of negative pions
with 12C and 160as a function of the pion kinetic energy. ( f i r Biiche and F'rzybilla,
1981.)
152 / APPENDIX A
ENERGY, MeV
. 1 Cross sections for elastic scattering of pions from I2C and 160aa a
function of pion kinetic energy, evaluated from experiments (data points) and deduced
for "j0by partial wave analysis (curves). (After Biiche and Pnybilla, 1981.)
BThe relation between the kinetic energy, T, the momentum,p, and the rest mass,
MC2, is T = (p2C2 + Wc4)* - M S .
154 1 APPENDIX A
k
5
i
0
5 lo'
LU
Fig. A.32. Illustrations of hadronic total and elastic cross sections as a function
of momentum, p,,-, in the laboratory system and of energy squared, s, in the center-
of-mass system. (ARer Particle Data Group, 1986.)
TABLEA.ll--Crosssections for the interactions of hadrons with nuclei at high energies"
Nucled Nuclear'
total inelastic Nucled Nucled
Atomic Atomic c m cross collision interaction
number weight section ~ection length length
Material Z A UT UI AT A1 -
cm2 x loz4 cm2 x 10% g~ m - ~ gc ~ - ~
Hz 1 1.01 0.0387 0.033 43.3 50.8
D2 1 2.01 0.073 0.061 45.7 54.7
He 2 4.00 0.133 0.102 49.9 65.1
Li 3 6.94 0.211 0.157 54.6 73.4 $
Be 4 9.01 0.268 0.199 55.8 75.2
z
C
Nz
6
7
12.01
14.01
0.331
0.379
0.231
0.265
60.2
61.4
86.3
87.8 s
oz
Ne
A1
8
10
13
16.00
20.18
26.98
0.420
0.507
0.634
0.292
0.347
0.421
63.2
66.1
70.6
91.0
96.6
106.4
E
0
Si
Ar
14
18
28.09
39.95
0.660
0.868
0.440
0.566
70.6
76.4
106.0
117.2 g
U)
V)
Fe 26 55.85 1.120 0.703 82.8 131.9 M
Cu
Sn
29
50
63.54
118.69
1.232
1.967
0.782
1.21
85.6
100.2
134.9
163 8
Xe 54 131.3 2.120 1.29 102.8 169 U)
W 74 183.85 2.767 1.65 110.3 185
Pb 82 207.19 2.960 i .77 116.2 194 b
0
Fig. A.33. Calculated differential cross sections as a function of emitted proton energies a t a particular angle, of residual mass yields and 5
of residual excitation energies for protons incident on aluminum and tantalum. The figures show comparisons between three intranuclear-
cascade models. Solid curves and dots: JINR model (Barashenkov et al., 1969);dashed curves and squares: BNLColumbia model (Chen et d.. __
1971);dotted curves and crosses: ORNL model (Bertini, 1969).(Excerpted h m Barashenkov et d., 1972.)
160 / APPENDIXA
J+ are the longitudinal and transverse momentum of the emitted
secondary particle in the center-of-mass system, and where Em is
the total center-of-mass energy for the incident and emitted particle.
In order to simplify the sampling in Monte-Carlosimulations, Ranft
furthermore made the approximation of factorizing the production
cross section into a function of 2pJE, times a function of h.
The functional form of the approximation formulas is taken from
the Hagedorn-Ranft thermodynamic model (Hagedorn and Ranft,
1968;Ranft, 1970).In this theory ofhadronic matter at high tempera-
ture, a hadron-hadron collision gives rise to a system with a tempera-
ture corresponding to 160 MeV that emits secondary particles in a
process similar to boiling. The thermodynamic model was developed
for hadron-hadron collisions, but Ranft later adapted it to fit the
experimental proton-nucleus cross sections obtained by Allaby et al.
(1970) at 19.2 GeVlc. For secondary proton and neutron production
in nucleon-nucleus collisions, RanfYs formulas have the form
where M,,is the pion mass. Numerical values of the parameters a,,
and b, are given by Ranft (1972) and up-dated by Chirikov et al.
(1977) for interactions with H, Be, Al, Cu, and Pb targets.
Equations (A.53) and (A.54) take into account mainly very ener-
getic secondary particles that are emitted with an angular distribu-
tion strongly peaked in the forward direction. They do not take into
account the nucleons with energies up to a few hundred MeV and
the somewhat more diffuse angular distributions that are ejected
from the nucleus as the result of the intra-nuclear-cascade process.
Ranft and Routti (1972), therefore, used a particle-production spec-
trum that consists of two parts: (a)the high-energy spectrum given by
Equations (A.53) or (A.54); and (b) a low-energy spectrum calculated
from a parameterized fit to the results of the intranuclear-cascade
calculations of Bertini (1963) and Barashenkov et d.(1%9,1972).
Illustrative results pertaining to the number of particles produced
A.4 NUCLEAR CROSS SECTIONS / 161
per incident hadron and to the energy used for their production, are
shown in Table A.12 for part (a) and in Figure A.34 for part (b).
Neutral pions decay very close to the point where they are produced
because of their short mean life (0.83 x lo-"% in their rest system).
With a probability of 98.8 per cent they decay into two gamma rays,
and with probability of 1.2 per cent they d a y into a single gamma
ray and an electron-positron pair. The gamma rays, which have a
combined energy equal to the sum of the TO rest energy (134.96
MeV) and the kinetic energy, initiate a photon-electron cascade. The
combined cross sections for hadron and electromagnetic cascades and
a Monte-Carlo calculation of the production of photons, electrons,
and positrons from primary proton beams have been described by
Ranft and Nelson (1979).
The cross sections needed for the calculation of the electromagnetic
cascade have been discussed in Sections A.l and A.2, with the excep-
tion of photon-nucleon and photon-nucleus interactions a t very high
energies. Up to 25 MeV, the photoabsorption process is dominated
by the giant dipole resonance already mentioned in Section A.1. At
energies from 50 to 200 MeV, photons mainly interact with proton-
neutron clusters inside the nucleus. The interaction cross section is
proportional, then, to the cross section for the photodisintegration of
deuterium, and the cross section is denoted, therefore, as a quasi-
deuteron cross section. At energies above the threshold for pion
-
production (- 135 MeV for a0and 140 MeV for a'), the photons
interact with the constituent nucleons of the nucleus, which are
raised to a n excited state and then decay into a stable nucleon and
a pion. The reaction mechanisms in the three regimes are indicated
in Figure A.35 for the case of photon interactions with lg7Au.Photon
cross sections for interactions with 160below the pion threshold are
shown in Figure A.36, and the cross sections for aoand a +production
from protons are shown in Figure A.37.
The mean life of negative and positive pions is rather long (2.60
x lo-' s in their rest frame), so that most of them are slowed down
to rest before decaying. Depending on the initial pion energy, from
0.1 to 1per cent of the pions decay into a muon and a neutrino while in
flight. These muons are quite important for shielding considerations
because they are highly penetrating. They do not interact strongly
with nuclei and are slowed down mainly by electromagnetic interac-
tions with atoms. Joseph (1969) has calculated muon range-energy
tables according to Bethe's stopping power theory by applying a
density-effect correction, and a small correction for muon brems-
strahlung, pair production and nuclear excitation.
TABLEA.12-Multiplicities N iand inelasticities Ki of proton and pion production in proton-nucleus collisions calculated by Ranft and 1
Routti (1972) from Equations (A.53 and A.54). Ni represents the average number of emitted particles per incident proton. K, is the fraction
of the available energy (incident proton energy minus residual excitation energy of the target nucleus) that is used to produce
Momentum (lab)
the various p r t i c l e ~ . ~
10 GeVlc 20 GeV/c 100 GeV/c
5z
M
Ni K Ni K Ni Ki
H+ 3e
H2
Be
1.28
1.16
0.28
0.25
1.32
1.21
0.30 1.36 0.31 *
0.27 1.24 0.28
A1 1.27 0.26 1.35 0.30 1.40 0.32
Cu 1.21 0.22 1.37 0.26 1.46 0.28
Pb 0.69 0.10 0.69 0.09 0.67 0.08
71 +
0 1 2 3 4 5
KINETIC ENERGY, GeV
0 1 2 3 4 5 6
KINETIC ENERGY, GeV
I I I I I
10 100 1000
PHOTON ENERGY, MeV
Fig. A.35. Cross section for photonuclear interaction with '("Au ae a function of
photon energy. (After Forkman and Schroder, 1972.)
A.4 NUCLEAR CROSS SECTIONS / 165
(y.p)+ly.n)
from p-shell \ I
0
100 200 300 400 500 600 700 800 900 1000
ENERGY, MeV
Fig. A.37. Cross eections for photomeson production as a function of photon
energy. C w e a: y + p -,+,curve b: y + p -, n + n+.(After ICRU,1978.)
Appendix B
Examples of Absorbed-Dose
and Dose-Equivalent
Calculations
0 10 20 30 40 50 60 70
ENERGY, MeV
Fig. B.1. Energy distributions of neutrons from 35- and 50-MeV 2Hon Be and 67-
MeV 'H on Li. The distributions shown have been normalized to unity.
B.l ABSORBED DOSE FROM NEUTRONS 1 169
this code in Section 8.3.3, the low-energy (515 MeV) neutrons are not
transported with HETC, but are transported with codes specifically
developed for these low-energy neutrons. In the calculations reported
here, this low-energy transport was done with the Monte-Carlo code
05R (Irving et al, 1965). The details of this low-energy transport are
given by Alsmiller (1975). The code 05R was not discussed in Section
8 because now it is not used often, but the results presented here
would not be appreciably different if the codes MORSE (Emmett,
1975) or MCNP (Thompson, 1979) had been used to transport the
low-energy neutrons.
The calculations were carried out by an analog Monte-Carlo
method (without use of variance-reduction techniques). All nuclear
reaction products (neutrons, protons, heavy nuclei, and photons) are
included in the calculations. The absorbed dose in a given volume is
defined as the energy deposited in the specified volume by charged
particles.
Rogers (1973) and Theus (1973) have measured the absorbed dose
along the axis of a n extended-slab, tissue-equivalent phantom for
various circular fields of radius R a t the surface of the phantom. The
experimental absorbed dose for R values of 17,13,9,5,and 3 cm are
shown as a function of depth in the phantom by the solid curves in
Figure B.2. The absolute normalization of the experimental data is
not known, but all of the experimental data in the figure have the
same normalization in terms of absorbed dose per incident neutron.
The histograms in Figure B.2 show the calculated absorbed doses
for the various field sizes used experimentally. To estimate absorbed
dose along the center line in the Monte-Carlo calculations, it was
necessary to average the calculations over a circular area about the
center line. The radial intervals used for this averaging process are
shown in the figure for each of the field sizes considered. Because the
calculations are statistical in nature, there is a statistical error
associated with each histogram value. In a few cases, this statistical
error-one standard deviation-has been indicated by an error bar
in Figure B.2.
To compare the calculations and experimental data in Figure B.2,
it was necessary to establish a normalization. This was done by
making the calculated and experimental results agree a t a depth of
5 cm for the case of R = 9 cm.It must be emphasized that all the
experimental data have the same normalization so only a single
normalization constant was used in obtaining all of the comparisons
shown in Figure B.2. The agreement between the calculated and
experimental absorbed doses is good for all depths and field sizes.
For neutrons from 35-MeV 2Hon Be, Rogers (1973) and Theus
(1973) also have measured the absorbed dose along the central axis
8 RADIAL INTERVAL
f0
8
,-
0
9 to
X
$ 8
ui
$
0
to
0
m
a:
0
a
m to
u
8
0
0 5 40 (5 20 25
DEPTH, cm
Fig. B2. Center-line absorbed dose relative to depth for cylindrically collimated
neutrons from 35-MeV2Hon Be incident on a tissue-equivalent phantom. The absolute
calculated absorbed doses, in units of Gray per incident neutron on the phantom in
the field of radius R, may be obtained by multiplying the calculated values in the
figure by the normalization factor 1-cosrarctan(171125)1where the value 125 cm is
1-cos[arctan(R/125)1'
the distance between the neutron source and the front of the phantom.
B.l ABSORBED DOSE FROM NEUTRONS 1 171
of a tissue-equivalent phantom with a collimator that produces a 10-
x 10-cm field a t the surface of the phantom. A comparison between
the calculations and experimental data for this case is shown in
Figure B.3. To estimate the center line absorbed dose, the calculated
results have been averaged over a 6- x 6-cm square about the center
line. The calculated results in the figure have an absolute normaliza-
tion. The experimental results have an unknown normalization and,
therefore,have been made to agree with the calculated absorbed dose
at a depth of 5 cm. The error bars on the calculated results are again
statistical only and represent one standard deviation. The agreement
between the calculated and experimental center-line absorbed dose
in the case of the square field is approximately the same as that
obtained in the case of the cylindrical fields.
h
c
e
C
3
(I)
C
C)
c
(I)
E
0
c
.-
Y
DEPTH, cm
Fig. B.3. Center-line absorbed dose relative to depth for neutrons from 35-MeV
'H on Be incident on a tissue-equivalent phantom. The collimator produced a square
10- x 10-cm field at the surface of the phantom, but the calculated results have been
averaged over a 6- x 6-cmsquare about the center line.
172 1 APPENDIXB
In the case of the square field discussed in the previous paragraph,
Todd et al. (in Rogers, 1973)measured the absorbed dose as a function
of distance from the beam center line at depths of 5 and 10 cm in the
phantom. In Figure B.4, the calculated doses are given with absolute
units and the experimental results have been normalized to agree
with the calculations at the point shown in the figure. The distance
from the centerline (abscissa in Figure B.4) is measured parallel to
the small distance of the field, that is, a t the surface of the phantom,
where the field is 10 x 10 cm, the distance from the centerline to the
edge of the field is 5 cm. The calculated results at depths of 5 and 10
cm shown in the figure were obtained by averaging over the depth
intervals 4 to 6 cm and 9 to 11 cm, respectively. The agreement
between the calculated and experimental absorbed doses is quite
good even at distances from the center line that are irradiated only
by scattered radiation.
Galactic cosmic rays extend to very high energies and are conse-
quently very penetrating. This section considers the shielding of
manned space vehicles against galactic cosmic-ray protons and alpha
particles. The results presented are taken from the work of Santoro
et al. (1973).
The geometric model considered is a shield consisting ofa spherical
shell with a spherical tissue receptor a t the center. As indicated in
Figure B.5, the inside radius of the shield is 150 cm, the shield
thickness is re,and the receptor has a radius, r,, of 15 cm.
The hfferential omnidirectional flux density of protons a t solar
minimum and maximum and of alpha particles a t solar minimum
are shown in Fig. B.6. These flux density spectra are taken from
the review of McDonald (1969). The spectra applicable during solar
minimum were used in the calculations presented here. In the calcu-
lations, the radiations were assumed to be isotropically incident on
the spacecraft.
The calculations for incident protons were carried out using the
high-energy transport code HETC (Chandler and Armstrong, 1972;
Alsmiller, 1975; Alsmiller, 1976) (see also Section 8.3.3), and the
calculations for alpha particles were carried out using a slightly
modified version of HETC. The modifications in HETC to allow the
transport of incident alpha particles were those required to predict
the energy-angle distribution of particle production from alpha parti-
B.2 SHIELDING OF MANNED SPACE VEHICLES 1 173
I I I I I
0 I 2 3 4 5 6 7 8 9
DISTANCE FROM CENTER LINE, crn
Fig. B.4. Absorbed dose relative to distance from the beam center line at two
depths when neutrons from 35-MeV 2Hon Be are incident on a tissue-equivalent
phantom. The collimator produced a square 10- x 10-cm field at the surface of the
phantom.
Fig. B.6. Differential omnidirectional galactic proton flux density at solar maxi-
mum and minimum and alpha-particleflux density at solar minimum.
C)PRIMARY PROTONS
a SECONDARY PROTONS
+ HEAVY NUCLEI
X CHARGED PIONS
PHOTONS FROM NEUTRAL PIONS
2 + ELECTRONS, POSITRONS, PHOTONS
X MUONS
1 0 - ' * 1 I I I I I I I I i I I l
0 4 2 3 4 5 6 7 8 9 4 0 4 4 4 2 4 3 4 4 4 5
DEPTH IN TISSUE, g ~ r n - ~
Fig. B.7. Contributions from the various types of particles ta the absorbed-dose
rate as a function of depth in tissue when a galajic, cosmic-ray, solar-minimum
proton spectrum is isotropically incident on a 20-g-~m-~-thick
aluminum shield.
O PRIMARY PROTONS
A SECONDARY PROTONS 0 PHOTONS FROM NEUTRAL PIONS
+ ELECTRONS. POSITRONS. PHOTONS
I X CHARGED PIONS X MUONS I
10-l2 1 I I I I I I I I I I I I I I
0 2 3 , 4 5 6 7 8 9 10 11 12 13 14 15
DEPTH IN TISSUE, g ~ r n - ~
Fig. B.8. Contributions from the various types of particles to the dose-equivalent
r a t . as a function of depth in tissue when a galactic, cosmic-ray, solar-minimum
proton spectrum is isotropically incident on a 20-g-~m-~-thick aluminum shield.
@ PRIMARY IONIZATION
A SECONDAY PROTONS
5 - + HEAVY NUCLEI -
X CHARGED PIONS
Q PHOTONS F R W NEUTRAL PIONS
2 - + ELECTRONS, POSITRONS. PHOTONS -
X MUONS
l l l l l l l 1 1 1 1~ 1 ~
0 2 4 6 8 40 12 44
DEPTH IN TISSUE, g ~ r n ' ~
Fig. B.9. Contribution of the varioua types of particles to the absorbed-dose rate
a8 a function of depth in tissue when a galactic, cosmic-ray, solar-minimum, alpha-
particle spectrum is isotropically incident on a 20-g-~rn-~-thickaluminum shield
-
X CHARGED PONS
-
V)
& 10-~ %
5
k
a:
I-
2 2
W
lo-,
1
f=?I 5
% 2
8 ,,-,.
0 2 4 6 8 10 12 14
DEPTH IN TISSUE, g cm7
Fig. B.lO. Contribution of the various types of particles to the dose-equivalent
rate as a function of depth in tissue when a galactic, cosmic-ray, solar-minimum,
alpha-particle spectrum is isotropically incident on a 20-g-~m-~-thick aluminum
shield.
SHIELD MATERIAL: A1
SHIELD THICKNESS= 59 cm2
GALACTIC COSMIC-RAY SOLAR-MINIMUM
PROTON SPECTRUM INCIENT
-
I
U)
X
I
I
I
1 SHIELD MATERIAL: Al
1 SHIELD THICKNESS- 5q cmZ
t - j - 4 = GALACTIC COSMIC-RAY SOLAR-MNIMJM
I
I
--4
I I - -- 7 P ---7----
-
.-
t u n a - PARTICLE SPECTRUM INUMNT
1 I I SHIELD MATERIAL: AI
0 '-I 1
I . SHELD THICKNESS = 2 0 p cm-2
GALACTIC COSMIC-RAY SOLAR-MINIMUM
PROTON SPECTRUM INCIDENT
t I ; -1
SHIELD MATERIAL: Al
SHIELD THICKNESS ~ 2 Q 0
I GALLCTIC COSMIC-RAY SOLAR-MINIMUM
. ALPHA-PARTICLE SPECTRUM INCIDENT
I I I I I I
0 2 4 6 8 10 12 14 !6
DEPTH IN TISSUE, g ~ r n - ~
Fig. B.ll. Absorbed-dose rate relative to depth in tissue for galactic, cosmic-ray,
proton and alpha-particle spectra isotropically incident on spherical-shell shields.
Solid lines: attenuated primary particles and all secondary particles; dashed lines:
attenuated primary particles; points: unattenuated primary particles.
taken from Alsmiller and Barish (1979) and (1981). The restriction
of consideration to neutron energies of 1400 MeV arises because a t
higher energies pion production becomes important and must be
taken into account. This limit on energy range is, however, sufficient
for most applications even around very high-energy accelerators
(Lindenbaum, 1961; Distenfeld, 1966; Rindi, 1975). The neutron and
photon flux densities were converted to dose-equivalent rates using
the data given in NCRP (1971) and in Claiborne et al. (1970)
If the energy and angular distribution of the neutron source are
known, the codes DOT and GRTUNCL can be used to calculate the
dose equivalent a t the detector position. The calculations, however,
require the use of a large computing machine, and they must be
repeated completely if results for a different source are required. It
would be convenient if these results could be used to estimate the
skyshine dose equivalent from different neutron-source, energy-
B.3 SKYSHINE FOR NEUTRONS 1 181
SHIELD MATERIAL A1
I
SHIELD THICKNESS. 59 ~ r n - ~
I
GALACTIC COSMIC-RAY SOLAR-MINIMUM
PROTON SPECTRUM INCIOENT
SHIELD MATERIAL: A1
SHIELD THICKNESS- 59 ~ r n - ~
GALACTIC COSMIC-RAY SOLAR-MINIMUM
ALPHA-PARTICLE SPECTRUM INCIDENl
SHIELD MATERIAL: A1
SHIELD THICKNESS. X) g ~ r n - ~
GALACTIC COSMIC-RAY SOLAR-MINIMUM
PROTON SPECTRUM INCIDENT
SHIELD MATERIAL: A1
SHIELD THICKNESS= X) g ~ r n - ~
GALACTIC COSMIC-RAY SOLAR-MINIMUM
ALPHA PARTICLE SPECTRUM INCIDENT
I I I I
10-'O 1 1 I I
I
0 2 4 6 8 10 12 14 16
DEPTH IN TISSUE, g ~ r n - ~
Fig. B.12. Dose-equivalent rate relative to depth in tissue for galactic, cosmic-
ray, proton and alpha-particlespectra isotropically incident on spherical-shellshields.
Solid lines: attenuated primary particles and all secondary particles; dashed lines:
attenuated primary particles; points: unattenuated primary particles.
where
H(r) = The dose equivalent a t a radial distance r from the source
and on the air-ground interface (i.e., a t z = O),
S(E,cos8) = The number of source neutrons per unit solid angle
a t energy E and polar angle 8 (the source emission is assumed
to be independent of the azimuthal angle),
182 / APPENDIX B
AIR
POINT SOURCE
I
(ASSUMED TO
EMIT ONLY INTO
Fig. B.13. Schematic diagram of skyshine geometry for which importance func-
tions are calculated.
For r e R13:
d) Spheroid with two axes d and the third axis ed (Kellerer, 1984)
2, = , - d2) ; z2 = Min(x,h) .
d ~ a x ( ox2 ((2.11)
U,(x) is the geometric-reduction factor of the cross sedion ( x l d ) as
follows:
Circular cross section:
u c x [ c o s - ( ) - ] , d . (C.12)
'Tr
Integrating the expression, one obtains for the unit cube the relation
derived by Piefke (1978):
188 / APPANDIX C
where
!R-rJ<x<R+r (c.16)
R < r and x<r-R.
b -
0
0 0.2 0.4 0.6 0.8 1 .O
RELATIVE DISTANCE, x/d
Fig. C.1. Geometric-reduction factore, U(x), or the point-pair distance distribu-
tions, p(x), for the autologoua case.
a: Geometric-reduction factors for spherical shells according to Equations ((3.4) and
(C.5). The distance, x, is given relative to the outer diameter, D. The ratio of inner
diameter, d, to outer diameter of the shell is indicated on the curves. The dashed curve
corresponds to the sphere [Equation (C.3)].
b: Geometric-reduction factors for ephemids according to Equation (C.7). Thedistance,
x, is given relative to the two equal diametere, d, of the spheroid. The ratio e of the
third diameter to d is indicated on the curves. The daahed curve represents the limit
of an infinitely extended prolate spheroid, while the dotted curve represents, for
comparison, an infinitely long circular rod of diameter d [Equations (C.10) and (C.12)I.
190 / APPENDIXC
Special cases:
An outer surface of the shell (R = r,):
[ (
U,,(x) = - cos-I Max
7r
- 1, "'z-31,
- Max ( a - r , o ) r x ~ a + r . (C.27)
Let h be the distance of A from the plane C of B, and let U,(y) be the
2-dimensional geometric-reduction factor (y,<y<y,) of B relative to
the projection of A in C. Then:
U2(y) = -
2n
1 [
sin-' i in (a +; up)
194 / APPENDIXC
therefore,
X
3
0 1 2 3 4 5
RELATIVE DISTANCE, x/R
Spherical Surface
0 1 2 3 4 5
RELA-l7VE DISTANCE, x/R
Fig. C.2. Geometric reduction factors, U(x), or the point-pair distance distribu-
tions, p(x), for the heterologous case.
a: Geometric-reduction factors for a sphere, A, of radius R, and a concentric sphere,
B, of radius r according to Equation (C.15). The ratio rlR is indicated on the curves.
The distance, x, is given relative to R .
b: Geometric-reduction factors for a spherical surface, A, of radius R and a concentric
sphere, B, of radius r according to Equation (C.16). The ratio rlR is indicated on the
curves. The distance, x, is given relative to R.
196 / APPENDIX C
0 1 2 3 4 5
RELATIVE DISTANCE, x/h
-
I
0
0 1 2 3 4 5
RELATIVE DISTANCE, x/R
Honomry Members
LAUMSTON
S. TAYLOR,Honomry President
NCRP Commentaries
No. Title
1 Krypton-85 in the Atmosphere-With Specific Reference to
the Public Health Significance of the Proposed Controlled
Release at Three Mile Island (1980)
2 Preliminary Evaluation of Criteria for the Disposal of Trans -
uranic Contaminated Waste (1982)
3 Screening Techniques for Determining Compliance with
Environmental Standards-Releases of Radionuclides to
the Atmosphere (1986),Revised (1989)
4 Guidelines for the Release of Waste Water from Nuclear
Facilities with Special Reference to the Public Health Sig-
nificance of the Proposed Release of Treated Waste Waters
at Three Mile Island (1987)
5 Review of the Publication, Living Without Landfills (1989)
6 Radon Exposure of the U.S. Population-Status of the
Problem (1991)
228 1 NCRP PUBLICATIONS
Symposium Proceedings
The Control of Exposure of the Public to Ionizing Radiation
i n the Event ofAccident or Attack, Proceedings of a Sympo-
sium held April 27-29, 1981 (1982)
NCRP Statements
No. Title
1 "Blood Counts, Statement of the National Committee on
Radiation Protection," Radiology 63,428 (1954)
"Statements on Maximum Permissible Dose h m Televi-
sion Receivers and Maximum Permissible Dose to the
Skin of the Whole Body," Am. J. Roentgenol., Radium
Ther. and Nucl. Med. 84, 152 (1960) and Radiology 75,
122 (1960)
X-Ray Protection Standards for Home Television Receivers,
Interim Statement of the National Council on Radiation
Protection and Measurements (1968)
Specification of Units ofNatural Uranium and Natural Tho-
rium, Statement of the National Council on Radiation
Protection and Measurements, (1973)
NCRP PUBLICATIONS 1 231
Other Documents
The following documents of the NCRP were published outside
of the NCRP Report, Cornmen- and Statement series:
Somatic Radiation Dose for the Geneml Population, Report
of the Ad Hoc Committee of the National Council on Radi-
ation Protection and Measurements, 6 May 1959,Science,
February 19,1960, Vol. 131, No. 3399, pages 482-486
Dose Effect Modifying Factors I n Radiation Protection,
Report of Subcommittee M-4 (Relative Biological Effec-
tiveness) of the National Council on Radiation Protection
and Measurements, Report BNL 50073 (T-471) (1967)
Brookhaven National Laboratory (National Technical
Information Service Springfield, Virginia)
The following documents are now superseded andlor out of print:
NCRP Reports
No. Title
1 X-Ray Protection (1931) [Superseded by NCRP Repart No. 31
2 R d m Pm&tion (1934) [Superseded by NCRP Report No. 41
3 X-Ray Protection (1936) [Superseded by NCRP Report No. 61
4 Radium Prdection (1938) [Supersededby NCRP Repart No. 131
5 Safe H a d i n g of Radioactive Luminow Compound (1941)
[Out of Print]
Medicul X-Ray Protection Up to Two Million Volts (1949)
[Superseded by NCRP Report No. 181
Safe Handling of Radioactive Isotopes (1949) [Superseded
by NCRP Report No. 301
Recommendations for WasteDisposal of Phosphorus32 and
Iodine-131 for Medical Users (1951) [Out of Print]
Radiological Monitoring Methods and Instruments (1952)
[Superseded by NCRP Report No. 571
Maximum Permissible Amounts of Radioisotopes i n the
Human Body and Maximum Permissible Concentmtion+
in Air and Water (1953) [Superseded by NCRP Report
No. 221
Recommendations for the Disposal of Carbon-14 Wastes
(1953) [Superseded by NCRP Report No. 81.1
232 / NCRP PUBLICATIONS
NCRP Proceedings
No. Title
2 Quantitative Risk in Standards Setting, Proceedings of the
Sixteenth Annual Meeting held on April 2-3, 1980 [Out
of Print]
Index
Adjoint methods, 19, 31 Fluorescence radiation, 97
Attenuation coefficient, 111 Flux density, 7
Auger electrons, 97 flux-density distribution, 7
Bremsstrahlung, 122 Geometrical factors, 74
Charged-particle equilibrium, 13 autologous geometry, 82, 186
Compton scattering, 97 heterologous geometry, 82, 188
Continuous Slowing-Down Geometric reduction factor, 76, 81
Approximation, 23, 40, 126 compilation of, 186
Cross section, 9, 20 Green's theorem, 8
assumption for calculating, 21 Hadrons, 153
distributions of, 9 Hartree-Fock wave functions, 103
energy-transfer coefficient, 13 Importance function, 181
list of references, 25 Internal dosimetry. 68
macroscopic, 9 Kenna, 13
mean free path, 21 charged-particle equilibrium, 13
mirroecopic, 9 kerma factore. 146
scattering center. 21 Klein-Nishina equation, 103
Discrete-ordinate method, 36 Linear energy transfer (LET), 84
codes for, 49 unrestricted (LET), 84
dosimetry calculations, 60 Mean excitation energy, 131
neutron-photon transport, 49,67 Mean free path, 21
Dose, 1 Microdoeimetry, 3
absorbed-dose rate, 2 specific energy, 2
calculation of, 3, 6 Moment method, 24,47
deterministic approach, 4 Monte-Carlo technique, 6, 61
energy imparted, 1 and CSDA, 66
h m neutrons, 167 efficiency of, 63
geometric considerations, 74 eventby-event, 65
internal dosimetry, 58 transport codes, 70
ionizing radiation, 1 variance-reduction metho&, 66
kerma, 13 Neutron cmsa sections, 139
measurement of, 3 Neutron-photon transport, 49,67
Monte-Carlo method, 6.72 Nucleon-meeon tramport, 70
stochastic quantities, 1 Pair production. 103
total dose, 10 Particle current density, 7
banaport equatione, 5, 8 Particle density, 6
Dose equivalent, 84,167 Phase space,8
Doae rate, 2 Photoelectric effect, 93
Elastic .scattering, 113 Photonuclear effect, 107
electrons, 113 Pions, 147
neutrons, 140 Point-pair dietanee distribution, 76,
protons, 117 81
Electron-photon transport, 68 Point-source kernal, 76
Energy-absorption coefficient, 111 differential, 75
Energy conservation, 29 integral, 75
Energy-transfer coefficient, 13 isotropic, 78
Fano's theorem, 29 Probability density function, 61
INDEX / 235
Quality factor, 84-85 plane muree. 16
Radiation equilibrium, 51 point aourc8,16
quasi-equilibrium, 56 reciprocity, 33
traneient equilibrium, 57 simplified repre~entation,15
Random (pseudorandom) number Specific Energy, 2
generator, 62 Stochastic quantities, 1
Rayleigh scattering, 106 Stopping power, 12, 128
Receptor, 17 denaity-effect correction, 132
absorbed dose in, 74 mass stopping power, 12
anthropomorphic, 17 'Ranslation operator, 27
receptor-free field, 17 lhmsport equations, 5,8,26
reciprocity, 33 demity d i n g , 28
for two kind8 of particles, 12
Reciprocity theorem, 77 in different coordinate eyatems, 34 .
Scattering operator, 27 integral form of, 26
Shielding, 172 iterative solution of, 27
Skyahine for neutrons, 178 linear and nonlinear, 10
fhurces, 6, 14 Neumann expansion, 28
fictitious s o m a . 15 radiation equilibrium, 38
monoenergetic, 16 euperposition, 30
non-uniform, 57 transport codes, 67