Precise measurements of the absolute γ

Applied Radiation and Isotopes 102 (2015) 15–28
Contents lists available at ScienceDirect
Applied Radiation and Isotopes
journal homepage: www.elsevier.com/locate/apradiso
Precise measurements of the absolute γ-ray emission probabilities of
223
Ra and decay progeny in equilibrium
S.M. Collins a,n, A.K. Pearce a, P.H. Regan a,b, J.D. Keightley a
a
b
National Physical Laboratory, Hampton Road, Teddington, Middlesex TW11 0LW, United Kingdom
Department of Physics, University of Surrey, Guildford, Surrey GU2 7XH, United Kingdom
H I G H L I G H T S
Discrepancies found within currently published γ-ray emission probabilities.
Absolute γ-ray emission probabilities of decay series in equilibrium determined.
Significant improvement in precision of measured values.
Closer agreement between deduced and experimental α transition probabilities.
Correlation coefficients presented for γ-emissions of 223Ra, 219Rn and 211Pb.
art ic l e i nf o
a b s t r a c t
Article history:
Received 25 September 2014
Received in revised form
3 March 2015
Accepted 14 April 2015
Available online 15 April 2015
Precise measurements of the absolute γ-ray emission probabilities have been made of radiochemically
pure solutions of 223Ra in equilibrium with its decay progeny, which had been previously standardised by
4π(liquid scintillation)-γ digital coincidence counting techniques. Two high-purity germanium γ-ray
spectrometers were used which had been accurately calibrated using a suite of primary and secondary
radioactive standards. Comparison of the activity concentration determined by the primary technique
against γ-ray spectrometry measurements using the nuclear data evaluations of the Decay Data Evaluation Project exhibited a range of 18% in the most intense γ-ray emissions ( 41% probability) of the
223
Ra decay series. Absolute γ-ray emission probabilities and standard uncertainties have been determined for the decay of 223Ra, 219Rn, 215Po, 211Pb, 211Bi and 207Tl in equilibrium. The standard uncertainties of the measured γ-ray emission probabilities quoted in this work show a significant improvement over previously reported γ-ray emission probabilities. Correlation coefficients for pairs of the
measured γ-ray emission probabilities from the decays of the radionuclides 223Ra, 219Rn and 211Pb have
been determined and are presented. The α-transition probabilities of the 223Ra have been deduced from
P(γ þ ce) balance using the γ-ray emission probabilities determined in this work with some agreement
observed with the published experimental values of the α-emission probabilities.
Crown Copyright & 2015 Published by Elsevier Ltd. All rights reserved.
Keywords:
Nuclear data
223
Ra
219
Rn
215
Po
211
Pb
207
Tl
Gamma-ray emission probabilities
Gamma-ray spectrometry
Correlation coefficients
Radioactivity
1. Introduction
Radium-223 is a naturally occurring radionuclide, occupying
the later stages of the primordial decay series of 235U (see Fig. 1),
that makes up approximately 0.7200(51)% of naturally occurring
uranium (Rosman and Taylor, 1998). With a half-life of 11.4354
(17) d (Collins et al., 2015), 223Ra has undergone investigations for
use as a radiopharmaceutical, with successful clinical trials for
targeted radiotherapy of bone metastases and bone palliation that
occur from late-stage castration resistant prostate cancer
n
Corresponding author.
E-mail address: [email protected] (S.M. Collins).
http://dx.doi.org/10.1016/j.apradiso.2015.04.008
0969-8043/Crown Copyright & 2015 Published by Elsevier Ltd. All rights reserved.
(Michalski et al., 2013; Nilsson et al., 2007; Parker et al., 2013). As a
group II element, radium shares many chemical properties with
calcium and exhibits a high level of incorporation into metabolically active osteoblastic bone and tumour lesion sites (Bruland
et al., 2008; Nilsson et al., 2007). Coupled with the relatively short
dose deposition range of α-emissions this allows a highly targeted
cytotoxic dose of ionising radiation to a specific cancer site with
reduced damage to the bone marrow and other surrounding
healthy tissue, giving this treatment obvious advantages over the
relatively long energy deposition range of β-emitting bone-targeting radiopharmaceuticals e.g. 89Sr, 166Ho and 153Sm, that have
been used historically. Additionally, as a naturally occurring
radionuclide it is of interest as a potential radiotoxic hazard from
16
S.M. Collins et al. / Applied Radiation and Isotopes 102 (2015) 15–28
Fig. 1. Decay series of
naturally occurring radioactive material (NORM) and technologically enhanced NORM (TENORM) (Kathren, 1998).
The 223Ra nucleus decays by 100% α-emission to excited states
of 219Rn (Chechev, 2011c), the decay scheme of which is shown in
Fig. 2, with a small branch (4.7 10 10 relative to the α-decay
branch (Kutschera et al., 1985)) decaying via the spontaneous
emission of a 14C nuclear cluster to 209Pb (Rose and Jones, 1984).
The 219Rn further decays via a series of relatively short-lived (T1/
2 o37 min) α- and β-emitting decay progeny each decaying via
their respective excited states (Chechev, 2011a, 2011b, 2011c; Nichols, 2011; Kondev, 2013a, 2013b; Luca, 2010; Luca, 2011) with
associated γ-ray emissions, with the series terminating at the
stable nucleus of 207Pb. The 14C cluster decay mode of 223Ra has
not been investigated in the scope of this work.
In view of the growing importance of 223Ra in these applications the National Physical Laboratory (NPL) undertook a course of
work to provide an absolute standardisation of 223Ra (Keightley
et al., 2015) that would provide traceability to the SI unit of the
Becquerel. Two solutions of 223Ra were independently standardised, henceforth referred to as S1 and S2, using liquid scintillation
(LS) absolute standardisation techniques, the first performed in
2013 and the second in 2014. Initial measurements of the solutions
were made by high purity germanium (HPGe) γ-ray spectrometry
of the thirteen most intense γ-ray emissions (41%) of the 223Ra
decay series, listed in Table.1. The activity of each selected γ-ray of
the 223Ra decay series was determined using the nuclear data from
the evaluations of the Decay Data Evaluation Project (DDEP)
(Chechev, 2011c; Nichols, 2011; Kondev, 2013a; Luca, 2010). Whilst
235
U.
the weighted mean activity of the individual γ-ray emissions were
in good agreement with the activity determined by the LS standardisation techniques (see Fig. 3), the activity concentration calculated from the various γ-ray emissions had a range of 18%, as
shown in Table 1 and Fig. 4. We considered this indicated a statistically significant discrepancy in the existing published γ-ray
emission data, and that further measurements were required. This
conclusion reflected the findings of Kellett and Nichols (Kellett and
Nichols, 2013) who highlighted a disagreement between the αtransitions and values deduced from the P(γ þ ce) balance within the
219
Rn excited levels.
A review of the currently published normalised γ-ray emission
values (Blaton-Albicka et al., 1976; Briançon and Leang, 1968; Davidson and Connor, 1970a; Hesselink, 1972; Kossert et al., 2015;
Krien et al., 1970; Sheline et al., 1998) for the decay of 223Ra was
performed. In all cases, with the exception of Krien et al. (1970), no
specific information was provided regarding the full-energy peak
(FEP) efficiency calibration of the γ-ray spectrometers used. Such
FEP efficiency calibrations are critical to the accuracy of γ-ray
emission probability measurements. Comprehensive details of the
methodology used and the resulting FEP efficiency calibration
curve are therefore presented in detail in this article.
Maintaining a chemical separation of the different radionuclides present in the decay series of 223Ra is problematic due to
the short half-lives of the decay progeny and evolution of 219Rn;
therefore the measurements were made of the decay series in
equilibrium within an aqueous solution. Investigation of the γ-ray
emissions from the decay of 223Ra and its decay progeny indicated
S.M. Collins et al. / Applied Radiation and Isotopes 102 (2015) 15–28
Fig. 2. Decay scheme of
223
Ra to states in
17
219
Rn (Chechev, 2011c).
Table 1.
γ-ray nuclear data from the DDEP (Chechev, 2011a; Nichols, 2011; Kondev, 2013a;
Luca, 2010) and determined activity per unit mass of the S1 solution by HPGe γ-ray
spectrometry.
Energy/keV
Radionuclide
Intensity /%
Activity/kBq g 1
122.3
144.3
154.2
269.5
271.2
323.9
338.3
351.0
401.8
404.8
427.2
445.0
832.0
Final result
223
1.238(19)
3.36(8)
5.84(13)
14.23(32)
11.07(22)
4.06(8)
2.85(6)
13.00(19)
6.75(22)
3.83(6)
1.81(4)
1.28(4)
3.5(5)
55.4(9)
54.0(13)
53.9(12)
49.0(11)
50.9(10)
47.1(9)
48.0(10)
53.1(8)
50.8(17)
54.7(9)
54.1(12)
50.0(16)
52.1(8)
52.0(8)
Ra
Ra
223
Ra
223
Ra
219
Rn
223
Ra
223
Ra
211
Bi
219
Rn
211
Pb
211
Pb
223
Ra
211
Pb
223
Fig. 3. S1 activity per unit mass determined by various measurement techniques.
the measurements of the primary γ-ray emission of each radionuclide would be essentially unaffected by interferences from the
decay of the other radionuclides present, although some potential
convolutions would be present in the energy range of the main γray emissions of 223Ra (269.5 keV) with 219Rn (271.3 keV) and
219
Rn (401.8 keV) with 211Pb (404.8 keV).
A full uncertainty budget is detailed within this paper. All uncertainties are stated as standard uncertainties or combined
standard uncertainties as defined in the Guide to the Expression of
Uncertainty in Measurement (GUM) (BIPM, 2008).
2. Experimental method
2.1. Sample preparation
Two separate active solutions were supplied by Algeta ASA
(Norway) with a nominal activity of 50 MBq 223RaCl2 in a 10 mL
sodium citrate pH buffer solution. Following previous experience
at NPL with 223RaCl2 in this chemical format, the material was
diluted with an aqueous solution of 1 M HCl to reduce the risk of
hydrolysis and any associated loss of activity during transfer of the
active solution between vessels. Solution S1 was diluted to a
nominal activity concentration of 50 kBq g 1 in 1 M HCl, with 1 g
18
S.M. Collins et al. / Applied Radiation and Isotopes 102 (2015) 15–28
2.3. HPGe γ-ray spectrometry
Fig. 4. Activity per unit mass of the individual γ-ray emissions of 223Ra and decay
progeny in equilibrium compared to the absolute activity per unit mass determined
for S1 (solid line) with the standard uncertainty (dashed lines).
aliquots dispensed to three 2 mL ISO ampoules (ISO, 2010) for
measurement by γ-ray spectrometry and dispensed to sixteen
34 mL Wheaton type liquid scintillation vials. Solution S2 was diluted to a nominal activity concentration of 330 kBq g 1 in 1 M
HCl and dispensed to a further three 2 mL ISO ampoules (ISO,
2010) for γ-ray spectrometry and sixteen 22 ml Wheaton type liquid scintillation vials. In order to minimise the risks associated
with emanation of 219Rn, the ampoules were flame sealed and the
vial screw threads wrapped with PTFE tape.
2.2. Absolute standardisation
The S1 solution was standardised using the CIEMAT/NIST efficiency tracing method (Broda et al., 2007) and the 4π(LS)-γ digital
coincidence counting (DCC) technique (Keightley and Park, 2007;
Keightley and Watt, 2002). The S2 solution was standardised using
the 4π(LS)-γ DCC technique only. A description of the measurement methodology has been covered in-depth by Keightley et al.
(2015). The results of these measurements can be found in Table 2.
The activity per unit mass for the two solutions were determined as AS1 ¼ 52.35(17) kBq g 1 and AS2 ¼ 329.2(11) kBq g 1.
Table 2
Activity per unit mass of
techniques.
223
Two LN2 cooled HPGe γ-ray spectrometers, identified as BART
and LORAX, were used to perform the measurements. BART is an
n-type HPGe γ-spectrometer with a resolution (FWHM) of
1.79 keV at 1.3 MeV and a relative efficiency of 28%, LORAX is an
n-type semi planar HPGe γ-spectrometer with a resolution of
0.6 keV at 0.122 MeV. Both detectors were contained in identical
1.5 m 1 m 1 m Pb shields comprised of 10 cm thick aged Pb
walls covered with a 0.5 mm Cd and 0.7 mm Cu graded liner to
reduce effects from background radiation and Pb fluorescence
X-rays in the spectra. Aluminium optical bread boards were
mounted along the horizontal axis of the coffin with a kinematic
mounting plate holding a precision engineered sample holder attached to enable highly reproducible geometric source positioning
in front of the detector window.
The energy calibration of both detectors was performed using
the most intense γ-rays ( 41%) of 152Eu, covering an energy region
from 121.8 keV to 1408.0 keV. The energy calibration resulted in
peak centroids within 0.05 keV of the 152Eu γ-ray energies evaluated by the DDEP (Vanin et al., 2004).
The spectra were collected using a chain of CANBERRA analogue electronics (AFT Research amplifier 2025, Analogue-to-Digital
Converter 8715, AIM) connected to a PC running the CANBERRA
GENIE 2000 v2.1c software. The net peak area losses due to dead
time were corrected for by the electronics using the integrated
pile-up rejection (PUR)/live-time correction (LTC) circuit. An additional correction was required for pulse pile-up occurring from
random coincidence summing events which had not been captured by the integrated PUR/LTC circuit of the analogue electronics
(Rajput, 2010). This was determined empirically for each detector
using the decay corrected counts of a 99mTc source.
All spectra were analysed using CANBERRA GENIE 2000 v2.1c. The
fit to each peak was manually reviewed and adjusted where necessary using the CANBERRA GENIE Interactive Peak Fit application. The
photopeak areas were corrected for background and integrated decay
(Harms and Jerome, 2004) during the counting period.
2.3.1. Full-energy peak (FEP) efficiency calibration
The FEP detection efficiencies were determined using a suite of
traceable primary and secondary standards of γ-ray emitting
radionuclides to cover the photon energy range from 14 keV to
1836 keV. The calibration geometry was 1 g aqueous solution in a
2 mL ISO ampoule at perpendicular distances from the detector
window of 30 cm and 16.5 cm for detectors BART and LORAX
respectively.
Each calibration point was measured with multiple sources,
collecting at least 105 counts in the photopeaks of interest. Using
the evaluated nuclear data from the DDEP (BIPM, 2004) for each
radionuclide, the FEP efficiency was estimated using the equation:
ε keV =
N × Ct × Cd × Cp × Cc
t × m × A 0 × Iγ
where,
Ra solutions S1 and S2 determined by multiple
εkeV
Measurement technique
S1 activity/
kBq g 1
S2 activity/kBq g 1
CIEMAT/NIST efficiency tracing
4π(LS)-γ DCC
HPGe γ-ray spectrometry
4π APPC (α)-γ coincidence counting
4π APPC (α þβ)-γ coincidence
counting
Final result
52.4(4)
52.34(18)
52.0(8)
52.4(4)
52.5(10)
–
329.2(11)
325(5)
52.35(17)
329.2(11)
–
N
Ct
Cd
Cp
Cc
t
m
A0
Iγ
photopeak efficiency
total counts collected in photopeak
Integrated decay correction
Self-absorption correction to H20
Pulse pile-up correction
True coincidence summing correction
Live time of measurement
Mass of active material
Activity per unit mass of radioactive material
γ-ray emission probability per decay
(1)
S.M. Collins et al. / Applied Radiation and Isotopes 102 (2015) 15–28
Table 3
Full-energy peak efficiency calibration points for the HPGe γ-ray spectrometers
BART and LORAX.
Energy/keV
Radionuclide
14.41
26.24
35.49
46.54
59.54
88.03
122.06
136.47
165.86
320.08
391.70
514.00
661.66
834.84
898.04
1115.54
1173.23
1332.49
1836.05
57
Full-energy peak efficiency /%
BART
Co
Am
125
I
210
Pb
241
Am
109
Cd
57
Co
57
Co
139
Ce
51
Cr
113
Sn
85
Sr
137
Cs
54
Mn
88
Y
65
Zn
60
Co
60
Co
88
Y
241
LORAX
2
2.96(10) 10
1.29(5) 10 1
1.52(4) 10 1
1.647(23) 10 1
1.674(10) 10 1
1.67(4) 10 1
1.601(11) 10 1
1.537(24) 10 1
1.416(12) 10 1
8.54(6) 10 2
7.16(6) 10 2
5.61(5) 10 2
4.53(4) 10 2
3.716(14) 10 2
3.499(20) 10 2
2.929(22) 10 2
2.822(8) 10 2
2.534(14) 10 2
1.93(9) 10 2
9.0(3) 10 2
3.44(12) 10 1
–
4.26(6) 10 1
4.33(3) 10 1
4.29(8) 10 1
3.88(3) 10 1
3.54(6) 10 1
2.893(25) 10 1
1.107(10) 10 1
8.23(7) 10 2
5.63(5) 10 2
4.06(3) 10 2
3.038(12) 10 2
2.793(13) 10 2
2.162(16) 10 2
2.045(6) 10 2
1.770(5) 10 2
1.218(5) 10 2
As the chemical composition of the aqueous calibration standards varied, as is necessary to obtain a stable solution of the
various elements, corrections were calculated for each data point
to compensate for the difference in expected self-absorption
compared to H20. The corrections were determined using a farfield model for cylindrical samples (Parker, 1984). The linear attenuation coefficients were derived from the NIST XCOM database
(Berger et al., 1998).
The samples were set at large distances from the detector
window, such that the solid angles were low (approximately
0.026 sr and 0.075 sr for BART and LORAX respectively, assuming
the sample as a point source). Hence, no cascade summing coincidence corrections were applied. The measured FEP efficiencies
for each calibration point are listed in Table 3.
The FEP detection efficiency curve fitted to the calibration data
set was of the form:
ε(x, C) =
⎡ N
⎤
⎢
⎥
⎢ ∑ c jbj(x, λ)⎥
⎦
e⎣ j = 1
(2)
where
x = loge
E
511
(3)
with the FEP efficiency denoted ε and the basis functions bj (x)
selected were quartic b-splines (Cox, 1972). The only adjustable
parameter, for each value of j, was C. A series of six knots (where
the piecewise polynomials are joined at λ = loge(Eλ/511)) were
selected at intervals along the energy range; these knots were
positioned based on a visual inspection of the shape of the FEP
efficiency curve. The knot set was augmented by placing an additional four replicated knots both at the beginning and end of the
sequence; the total number of knots was therefore fourteen and
the number of fitted coefficients (N) was 9. The best fit was determined by a generalised least squares fit (Aitken, 1936) to the
calibration points spanning the range of energies that were covered for each detector, 14 keV(57Co)–1836 keV (88Y) for BART and
46 keV (210Pb)–1836 keV (88Y) for LORAX. The respective fits and
residuals of the fits can be seen in Fig. 5 and Fig. 6.
19
2.3.2. Measurements of 223Ra and detection of impurities
The S1 samples were measured solely on BART over a period of
approximately 20 d with a total of eight measurements; referred
to as dataset S1B. Each source was measured at least twice with a
typical measurement time of 50,000 s. A dead time of 0.50% was
recorded at the start of the measurement campaign and 0.23% at
the end of the measurement campaign.
The S2 samples were measured over a period of 29 d and 18 d
on both BART and LORAX respectively, with a total of nine measurements on each detector; the corresponding datasets are referred to as S2B and S2L. Each sample was measured three times
on each detector, with a typical counting time of 86400 seconds.
Dead times of 3.22% and 7.76% were recorded at the beginning of
the measurement campaign and 0.58% and 2.43% at the end of the
measurement campaign for BART and LORAX respectively. An example of a spectrum collected on BART can be seen in Fig. 7.
All peak areas were corrected for dead time, integrated decay,
pulse pile-up and background. Self-absorption corrections were
made to convert the FEP efficiency from H2O to 1 M HCl, this
correction was less than 0.2% at γ-ray emission energies greater
than 100 keV. The count rates of those γ-ray emissions produced
by the decay of 211Pb, 211Bi, 211Po and 207Tl have been corrected to
take into account that these radionuclides exist in a transient
equilibrium with that of 223Ra due to the half-life of 211Pb being
approximately 36.1 min (Sargent, 1939). This correction was determined using the Bateman equation (Bateman, 1910):
λp
λ p − λd
× BR
(4)
where λp and λd are the decay constants of the parent and
daughter respectively and BR is the branching ratio of the parent
to daughter. The correction of the γ ray emission probabilities for
these decay progeny was approximately 0.22%.
S1 and S2 were investigated for the presence of γ-ray emitting
impurities and the precursors 227Ac and 227Th. The presence of
227
Ac – T1/2 ¼7952.1 d (Browne, 2001) – and 227Th – T1/2 ¼18.68 d
(Browne, 2001) – in small quantities could affect the accuracy of
the determination of the γ-ray emission probabilities as they will
support the activity of the 223Ra, thus perturbing the expected
radioactive decay rate. Neither solution showed any indication of
impurities or the presence of 223Ra precursors, though additional
consideration of the precursor 227Ac was required as this radionuclide exhibits relatively weak α- and subsequent γ-ray emissions that render α- and γ-spectrometry impractical. The activity
of S1 was measured over a period of 72 d (3.6 times longer than
the measurement period for the γ-ray emission measurements)
using an ionisation chamber. There was no exhibited perturbation
in the radioactive decay or trends within the residuals of least
squares fit; Indicating that there was no significant presence of
227
Ac or 227Th within the solution. It has been assumed by comparison of the measurements of S2 over the measurement period
and the agreement between the calculated γ-ray emissions of S1
and S2 that there was no significant presence of 227Ac in this solution either.
3. Results
The absolute γ-ray emission probabilities of 223Ra and its decay
progeny were derived from the activity per unit mass determined
by absolute standardisation techniques and measurement by the
two HPGe γ-ray spectrometers. The absolute γ-ray emission intensities for each measurement were calculated using the equation:
20
S.M. Collins et al. / Applied Radiation and Isotopes 102 (2015) 15–28
Fig. 5. Full-energy peak efficiency curve and residuals of the B-spline fit for the
HPGe γ spectrometer BART.
Fig. 6. Full-energy peak efficiency curve and residuals of the B-spline fit for the
HPGe γ spectrometer LORAX.
Iγ =
N × Ct × Cd × Cp × Cc
t × m × A o × ε keV
(5)
The weighted mean of the measured absolute γ-ray emission
probability from each measurement for a dataset was deduced
using the uncertainties derived from the peak fitting software,
dead time and radioactive decay to determine the weight of contribution from each measurement. The final value was deduced as
the weighted mean of the three datasets with the FEP efficiency
uncertainty combined with the standard uncertainty of the
weighted mean of the γ-ray emission probability uncertainty from
each dataset. The reduction in the standard uncertainty of the
weighted mean was limited to the smallest of the standard uncertainty of the three datasets. The final uncertainties were determined by adding the additional systematic uncertainties in
quadrature to the standard uncertainty calculated from the
Fig. 7. An annotated spectrum of
223
Ra and decay progeny in equilibrium.
weighted mean. The uncertainty budget is described in detail in
Section 4.
As the photopeaks below 100 keV are mainly the convoluted
X-ray emissions of 223Ra and its decay progeny, only results for the
observed γ-ray emissions from the decay of 223Ra and the resulting
decay progeny that have energies greater than 100 keV are presented in Table 4. These are presented as the absolute γ-ray
emission intensities per 100 disintegrations. Additionally the
radionuclide source of the γ-ray emission is listed as deduced from
the evaluations of the DDEP (Chechev, 2011a, 2011b, 2011c; Nichols, 2011; Kondev, 2013a, 2013b; Luca, 2010, 2011).
The relative γ-ray emission probabilities of this work and those
previously published for 223Ra (Blaton-Albicka et al., 1976;
S.M. Collins et al. / Applied Radiation and Isotopes 102 (2015) 15–28
Table 4
Absolute γ-ray emission probabilities per 100 decays of
223
Ra and decay progeny in equilibrium.
Energy
(keV)
Source Iγ (%)
Energy
(keV)
Source
Iγ (%)
103.9(5)
106.7(4)
223
0.0119(6)
0.0213(11)
323.9(6)
328.4(6)
223
110.8(5)
122.3(5)
130.6(5)
144.3(5)
154.2(5)
158.7(5)
175.6(5)
177.4(5)
179.7(5)
221.4(5)
224.0(5)
249.4(5)
251.9(5)
255.1(5)
269.5(6)
271.3(6)
288.2(6)
223
Ra
Ra
219
Rn
223
Ra
223
Ra
223
Ra
223
Ra
223
Ra
223
Ra
223
Ra
219
Rn
223
Ra
223
Ra
223
Ra
223
Ra
219
Rn
223
Ra
0.0512(10)
1.312(6)
0.1478(10)
3.481(16)
6.02(3)
0.749(4)
0.01578(10)
0.0426(8)
0.1613(10)
0.0304(10)
0.0056(14)
0.0375(9)
0.0640(11)
0.0499(13)
13.37(7)
10.75(6)
0.1498(16)
333.9(6)
338.3(6)
342.9(6)
351.1(6)
355.5(6)
361.7(6)
363.0(6)
368.4(6)
371.7(6)
372.9(6)
376.2(6)
383.3(5)
386.3(5)
390.1(5)
401.8(6)
404.8(6)
427.1(6)
Ra
Ra, 211Po,
207
Tl
223
Ra
223
Ra
223
Ra, 211Pb
211
Bi
223
Ra
211
Pb
223
Ra
223
Ra
223
Ra
223
Ra
223
Ra
223
Ra
223
Ra, 219Rn
223
Ra
219
Rn
211
Pb, 215At
211
Pb
293.6(5)
313.7(6)
219
0.0688(7)
0.0276(5)
430.4(6)
432.4(6)
223
0.0206(19) 619.8(6)
0.0297(14) 623.4(5)
223
Ra
Ra
223
Rn
Pb
211
223
223
21
Ra,
Ra
211
Pb
Source
Iγ (%)
Energy
(keV)
Source
Iγ (%)
3.655(18)
438.8(6)
0.2021(16) 445.0(6)
215
0.0533(7)
1.218(6)
675.4(6)
676.9(6)
211
0.0058(6)
0.0184(5)
0.0756(6)
2.605(13)
0.1958(21)
13.17(7)
0.0124(15)
0.0341(7)
0.0192(9)
0.0134(4)
0.435(3)
0.1133(13)
0.0056(4)
0.0023(6)
0.0052(7)
0.0053(7)
6.57(3)
4.011(19)
1.890(9)
219
0.0011(5)
0.0083(3)
0.0013(5)
0.0022(4)
0.0453(5)
0.0021(6)
0.0659(8)
0.0028(9)
0.0033(6)
0.0026(6)
0.0028(6)
0.0026(7)
0.0035(4)
0.0043(5)
0.0029(13)
0.0867(12)
0.0543(7)
704.6(7)
707.8(7)
711.4(7)
727.4(7)
766.4(7)
831.9(7)
835.6(7)
865.8(6)
891.3(7)
897.8(7)
1014.7(7)
1074.5(7)
1080.1(7)
1103.3(8)
1109.5(8)
1196.2(8)
1234.3(8)
211
0.0056(12) 1270.7(8)
0.0082(8)
211
Briançon et al., 1968; Davidson and Connor, 1970a; Hesselink,
1972; Krien et al., 1970; Sheline et al., 1998), 219Rn (Blaton-Albicka
et al., 1976; Briançon and Leang, 1968; Dalmasso and Maria, 1967;
Davidson and Connor, 1970b; Krien et al., 1970; Liang et al., 1998),
215
Po (presented as absolute values for comparison) (Briançon and
Leang, 1968; Davidson and Connor, 1970b) and 211Pb (Blaton-Albicka et al., 1976; Briançon and Leang, 1968; Cothern and Connor,
1965; da Silveira et al., 1971; Dalmasso and Maria, 1967; Davidson
et al., 1967; Giannini et al., 1962a; Gorodetzky et al., 1968; Hamilton and Davies, 1968; Hindi et al., 1988; Mead and Draper,
1965; Vandenbosch et al., 1963) are presented in Tables 5–8. The
269.5 keV γ-ray emission of the 223Ra decay is convoluted with the
271.3 keV γ-ray emission of the 219Rn decay. Hence, the relative γray emission intensities of 223Ra were deduced as the γ-ray
emission intensity normalised to the 154.2 keV γ-ray emission
probabilities, as this transition has no known convolutions and
should therefore be a reliable normalisation point. The 271.3 keV
γ-ray emission has been used as the normalisation for the γ-ray
emissions of 219Rn. Though it is convoluted with the 269.5 keV of
223
Ra, the only other significant emission of 219Rn, the 401.8 keV γray emission, is also convoluted by the 404.8 keV γ-ray emission of
211
Pb. The 211Pb values have been deduced relative to the 351.1 keV
γ-ray emission.
4. Uncertainties
4.1. Uncertainty components of the γ-ray emission probabilities
A summary of the uncertainty components for the absolute γray emission probability of the 269.5 keV γ-ray emission is presented in Table 9. The assigned uncertainty values of the components were combined in quadrature to determine the final uncertainty value. The dominating uncertainty components for the
most significant γ-ray emissions of the 223Ra decay series (where
the statistical uncertainties are less than 0.1%) are due to the
standard uncertainty of the standardisation of the activity per unit
mass, the FEP efficiency calibration and efficiency stability of the
detector chain.
The FEP efficiency uncertainties were calculated from the least
Energy
(keV)
462.8(6)
487.3(5)
500.2(6)
504.1(6)
517.6(6)
522.6(6)
527.6(6)
531.4(6)
537.5(6)
542.1(6)
545.9(6)
555.9(5)
564.4(5)
569.6(7)
573.7(7)
598.6(7)
609.3(7)
Po
Ra
223
Rn
Ra
223
Ra
211
Pb
219
Rn
223
Ra
223
Ra
223
Ra
223
Ra
223
Ra
223
Ra
219
Rn
219
Rn
211
Po, 207Tl
223
Ra
223
Ra
223
Ra, 219Rn,
211
Pb
219
Rn
223
Ra
223
Pb
Rn
219
Pb
Rn
223
Ra
223
Ra
211
Pb
211
Pb
219
Rn
211
Pb
219
Rn
211
Po,
211
Pb
219
Rn
211
Pb
211
Pb
211
Pb
211
Pb
211
Pb
219
Pb
0.498(3)
0.0034(4)
0.0037(3)
0.0024(7)
0.685(4)
3.448(16)
0.00364(19)
0.00540(21)
0.00107(20)
207
Tl 0.2725(15)
0.0171(4)
0.00044(12)
0.01228(21)
0.00380(12)
0.1113(7)
0.01052(17)
0.00092(8)
0.00624(19)
squares fit of the photopeak efficiency calibration including the
input efficiency uncertainties (see Table 3) and taking into account
the correlations and co-variances due to nuclear data, standardisation, etc. The calculated uncertainty of the FEP efficiency calibration can be seen in Figs. 5 and 6 as the dashed line in the residual plots.
A systematic uncertainty of 5% was estimated for the dead-time
and pile-up correction, which when propagated using the median
recorded dead-time resulted in an overall uncertainty component
of 0.05%. An additional uncertainty component was incorporated
in the dead-time and pile-up uncertainty to take into account the
perturbation on the photopeak shape that occurs at increasing
count rates (Rajput, 2010). Within the range of dead times observed this was estimated to be insignificant.
A radioactive decay uncertainty was determined using the
evaluated uncertainty of the 223Ra half-life, in previous work by
the authors (Collins et al., 2015), of 11.4354(17) d, using the
median decay period of the measurements to the standardisation
reference time to propagate the radioactive decay uncertainty.
An estimated uncertainty of 0.10% was incorporated to account
for the small but inevitable quantity of true coincidence summing
events that may occur. As described previously in Section 3, due to
the small solid angle subtended by the detector the number of
events that would occur would be insignificant and therefore the
assigned uncertainty value would be sufficient to account for these
events.
It cannot be assumed that the detector response will remain
constant with small fluctuations that may occur in the detector
itself or in the detector electronics chain; these effects can be small
and hence obscured within the measurement series. This effect
should be accounted for especially if the measurements are made
over an extended period; the long term stability uncertainty has
been estimated using a 152Eu source which had been measured
over the measurement campaigns.
The geometric reproducibility uncertainty for the positioning of
the source in an identical geometry is 0.10% as the use of precision
engineered source holders in tandem with the kinematic optical
mounting system allows a very high level of positional
reproducibility.
22
S.M. Collins et al. / Applied Radiation and Isotopes 102 (2015) 15–28
Table 5
Normalised γ-ray emission probabilities of 223Ra observed in this work and published values. All values normalised to the 154.2 keV γ-ray emission. Energies marked with an
n
have interferences with γ-ray emissions from decay progeny in equilibrium.
Energy/keV This work
Briançon et al.
(1968)
Krien et al.
(1970)
Davidson and
Connor (1970a)
Hesselink
(1972)
Blaton-Albicka
et al. (1976)
Sheline et al.
(1998)
Kossert et al.
(2015)
χ 2 /(n − 1)
103.9(5)
106.7(4)
110.8(5)
122.3(5)
144.3(5)
154.2(5)
158.7(5)
175.6(5)
177.4(5)
179.7(5)
221.4(5)
249.4(5)
251.9(5)
255.1(5)
269.5(6)
288.2(6)
323.9(6)
328.4(6)n
333.9(6)
338.3(6)
342.9(6)n
355.5(6)
363.0(6)
368.4(6)
371.7(6)
372.9(6)
376.2(6)
383.3(5)n
386.3(5)
390.1(5)
430.4(6)n
432.4(6)
445.0(6)
487.3(5)
500.2(6)
522.6(6)
527.6(6)
531.4(6)
537.5(7)
542.1(7)
545.9(7)
573.7(7)
598.6(6)
609.3(7)n
623.4(5)
711.4(7)
727.4(7)
–
0.43(8)
1.04(17)
21(4)
57(7)
100(–)
13.0(16)
–
0.54(19)
2.8(8)
0.65(19)
0.67(19)
1.3(3)
1.1(3)
260(40)
3.0(5)
69(9)
3.7(5)
1.6(3)
50(6)
3.7(5)
–
–
–
10.1(13)
–
–
–
–
–
–
0.622(91)
22.5(20)
0.184(39)
–
–
1.30(16)
–
–
–
–
–
1.48(19)
0.93(13)
0.15(8)
0.065(19)
–
–
0.34(7)
0.81(7)
21.3(6)
57.9(17)
100(–)
12.2(4)
0.24(7)
0.83(8)
2.6(7)
0.54(10)
–
0.66(17)
–
243(7)
2.82(13)
71.5(20)
3.70(18)
1.85(15)
51(12)
4.14(23)
–
–
–
8.66(24)
–
–
–
–
–
–
0.63(10)
26.8(20)
0.24(7)
–
–
1.31(12)
–
–
–
–
–
1.85(17)
1.31(17)
–
–
–
–
0.42(11)
1.08(15)
26(4)
61(9)
100(–)
14.7(22)
–
–
2.9(5)
–
–
1.1(4)
1.0(4)
260(40)
2.8(5)
70(10)
3.1(7)
2.4(6)
50(8)
3.9(11)
–
–
–
10.5(19)
–
0.23(8)
0.29(11)
–
–
–
–
25(4)
–
–
–
1.3(3)
–
–
–
–
–
1.8(4)
1.21(22)
–
–
–
–
0.43(14)
0.47(20)
19.6(15)
62(6)
100(–)
11.9(11)
0.34(9)
0.79(14)
2.6(4)
0.56(10)
0.7(3)
1.1(4)
0.54(16)
225(16)
2.1(4)
60(5)
2.7(5)
1.6(3)
43(3)
1.6(4)
–
–
–
9.5(11)
–
–
0.07(5)
0.23(9)
0.05(4)
0.32(9)
0.54(14)
20.7(20)
0.18(9)
–
–
1.1(3)
–
–
–
–
–
1.5(3)
0.68(25)
–
–
–
0.35(6)
0.41(6)
1.05(15)
20(3)
57(8)
100(–)
12.1(17)
–
–
2.7(4)
–
–
1.24(23)
0.87(21)
260(30)
2.8(4)
71(8)
3.7(5)
1.42(24)
49(6)
3.9(5)
–
–
–
8.3(10)
1.9(3)
–
–
–
–
–
0.49(9)
22(3)
–
–
–
1.08(13)
–
–
–
–
–
1.65(19)
0.98(12)
–
–
–
0.34(5)
0.42(3)
1.02(8)
21.2(6)
57.3(19)
100(–)
12.2(4)
0.34(7)
0.83(8)
2.68(5)
0.63(10)
0.68(17)
0.73(24)
0.93(12)
244(8)
2.80(10)
70.0(21)
3.7(12)
1.78(11)
49.8(16)
3.9(3)
0.073(24)
0.49(12)
0.15(7)
8.5(4)
0.878(21)
0.22(7)
–
0.27(10)
0.12(5)
0.34(10)
0.61(5)
22.7(9)
0.195(25)
0.024(10)
0.024(10)
1.24(8)
0.024(10)
0.037(12)
0.024(10)
0.020(10)
0.020(10)
1.66(8)
1.00(6)
0.15(7)
0.063(17)
0.0049(24)
–
–
0.93(5)
21.63(12)
57.53(24)
100(–)
12.22(8)
0.284(16)
0.632(18)
2.54(5)
0.463(18)
0.60(5)
1.05(3)
0.73(3)
218.2(13)
2.42(3)
60.71(25)
3.28(3)
1.244(22)
43.35(17)
2.90(3)
–
0.292(10)
0.252(19)
8.3(3)
0.85(3)
0.093(5)
–
–
–
0.393(4)
0.496(21)
20.27(8)
0.100(11)
–
–
1.056(18)
–
–
–
–
–
1.425(14)
0.507(12)
–
–
–
7.8
0.8
2.3
0.9
0.2
–
0.3
0.9
4
1.2
1.3
0.1
1.2
2.1
4.2
3.6
7.3
1.6
6.6
2.8
15
15
2.3
1.7
10
440
2.0
2.7
2.8
0.9
1.2
1.3
3.0
3.7
0.1
0.6
1.9
1.5
1.4
1.9
4.0
1.6
2.3
102
0.3
0.1
8.2
0.198(10)
0.354(18)
0.850(17)
21.79(11)
57.8(3)
100(–)
12.27(6)
0.2620(17)
0.708(13)
2.678(18)
0.504(16)
0.622(14)
1.062(19)
0.829(21)
221.9(11)
2.488(27)
60.7(3)
3.36(3)
1.255(9)
43.25(21)
3.25(4)
0.206(25)
0.318(15)
0.223(5)
7.22(5)
1.882(22)
0.092(5)
0.038(9)
0.086(11)
0.088(11)
0.34(4)
0.493(23)
20.22(10)
0.137(5)
0.021(8)
0.034(9)
1.094(13)
0.047(16)
0.055(10)
0.043(9)
0.047(10)
0.049(21)
1.441(19)
0.902(10)
0.135(14)
0.061(5)
0.040(12)
Table 6
Normalised γ-ray emission probabilities of
219
Rn observed in this work and published values. All values normalised to the 271.3 keV γ-ray emission.
Energy/keV This work
Dalmasso and
Maria (1967)
Briançon et al.
(1968)
Davidson and
Connor (1970b)
Krien et al.
(1970)
Blaton-Albicka
et al. (1976)
Liang et al.
(1998)
Kossert et al.
(2015)
χ 2 /(n − 1)
130.6(5)
224.0(6)
271.3(6)
293.6(5)
401.8(6)
462.8(6)
517.6(6)
555.9(5)
564.4(5)
619.8(6)
676.9(6)
707.8(7)
835.6(7)
891.3(7)
1074.5(5)
1.40(14)
–
100(–)
0.64(6)
58(6)
–
0.44(10)
–
–
–
0.21(3)
1.18(23)
–
100(–)
0.70(14)
61(4)
–
0.44(4)
–
–
–
0.21(2)
–
–
–
–
1.05(25)
–
100(–)
0.59(15)
65.2(65)
–
0.22(5)
–
–
–
0.06(3)
–
–
–
–
1.21(10)
–
100(–)
0.51(27)
69(3)
–
–
–
–
–
–
–
–
–
–
1.10(11)
–
100(–)
0.72(5)
58.4(27)
–
0.41(3)
–
–
–
0.15(1)
–
–
–
–
1.7(2)
0.013(2)
100(–)
0.680(42)
59.0(23)
0.0015(3)
0.40(22)
0.0005(3)
0.014(3)
0.003(1)
0.160(20)
0.003(1)
0.015(3)
0.007(2)
0.003(1)
–
–
100(–)
–
60.9(5)
–
–
–
–
0.0469(78)
–
–
–
–
–
2.3
8.8
–
0.7
1.4
3.3
3.5
14
19
25
4.7
60
30
1.0
0.6
1.375(9)
0.052(13)
100(–)
0.640(7)
61.1(3)
0.010(5)
0.421(4)
0.024(6)
0.033(3)
0.052(11)
0.171(4)
0.031(4)
0.0339(17)
0.0100(19)
0.0041(10)
0.015(7)
da Silveira
et al.
(1971)
–
–
0.198(15)
30.6(4)
14.44(10)
–
–
3.776(28)
5.166(38)
25.91(17)
–
–
–
–
–
–
–
–
29.3(9)
13.9(4)
0.045(6)
–
3.6(1)
4.94(16)
26.7(8)
0.042(6)
0.129(8)
0.095(6)
0.033(4)
0.90(3)
0.072(5)
–
0.043(4)
0.20(3)
0.326(24)
30.2(14)
14.2(7)
–
0.130(8)
3.6(3)
5.1(4)
25.4(20)
0.033(4)
0.122(8)
0.090(7)
0.049(5)
0.82(6)
0.081(6)
–
0.057(5)
–
–
26(5)
12.5(25)
–
–
3.8(11)
–
24.8(25)
–
–
–
–
1.07(16)
–
–
–
Bi.
211
Briançon
et al. (1968)
Gorodetzky
et al. (1968)
Hamilton and
Davies (1968)
Dalmasso and
Maria (1967)
–
–
34(4)
22(3)
–
–
5.5(4)
6.1(4)
34.2(13)
–
0.38(19)
–
–
1.46(19)
–
–
–
A total of 83 absolute γ-ray emission probabilities have been
determined in this work. The determined probabilities show significant improvement in their precision compared to some previously published values. For example the uncertainty of the
269.5 keV γ-ray emission probability is a factor of 4 lower than
previously quoted by Sheline et al. (1998) and a factor of 2 more
precise than quoted by Kossert et al. (2015). Previously no detailed
uncertainty budgets have been published therefore it is difficult to
specify the exact reasons for the improvement. One possible reason is the use of higher activity samples of 223Ra that has allowed
313.7(6)
361.7(6)
404.8(6)
427.1(6)
504.1(6)
675.4(6)
704.6(7)
766.4(7)
831.9(7)
865.8(6)
1014.7(7)
1080.1(7)
1103.3(8)
1109.5(8)
1196.2(8)
1234.3(8)
1270.7(8)
5.1. Results
0.24(3)
–
30.0(9)
13.5(6)
0.12(2)
–
3.77(19)
5.55(28)
29.8(7)
0.050(8)
0.14(1)
0.120(12)
0.040(6)
1.15(8)
0.10(1)
0.010(2)
0.070(7)
Hindi
et al.
(1988)
5. Discussion
Davidson
et al. (1967)
Following the determination of the fitted parameters C, from Eq. (2),
and uncertainty matrix UC the covariance of the fitted parameters was
determined by the linear propagation of uncertainties for multivariate
models as described in supplement 2 of the GUM (BIPM, 2011). The
correlations between pairs of the most intense measured γ-ray
emissions (41%) associated with the decay of the radionuclides
223
Ra, 219Rn and 211Pb in the series are given in Tables 10–12.
Mead and
Draper
(1965)
(7)
Vandenbosch
et al. (1963)
Cov(xi , xj )
Var(xi )Var(xj )
Giannini et
al. (1962a)
r (x i , x j ) =
Energy/keV This work
where
Pb observed in this work and published values. All values are normalised to the 351.1 keV γ-ray emission of
(6)
211
⎡ r (P , P ) … r (P , P ) ⎤
1 N
⎢ 1 1
⎥
R=⎢
⋮
⋱
⋮
⎥
⎢⎣r (PN , P1) ⋯ r (PN , PN)⎥⎦
Table 8
Normalised γ-ray emission probabilities of
The majority of the uncertainty components are significantly
correlated between the solutions S1 and S2 (as the standardisation
methodology is identical) as well as between the two detectors.
This is due to the process of determining the efficiency curve
where the standard solutions used to calibrate the detectors were
of similar provenance i.e. ionisation chamber, and identical nuclear
data used. Failure to account for these correlations will often undermine the quality of the measurements and calculations relying
on such data. These effects are not routinely dealt with in publications of γ-ray emission probabilities.
In this work the calibration standards were all of NPL provenance and therefore the off-diagonal terms of the matrix used to
weight the generalised least-squares fit were estimated from the
ionisation chamber calibration factor uncertainties for the radionuclides concerned. Correlations due to sample preparation decay
and dead time corrections were not included as these were relatively small contributors to the final uncertainty. Correlations in
published γ-ray emission values may be more significant; however
data were not available to allow the covariance to be estimated.
The correlation matrix R as defined in GUM supplement 2 (BIPM,
2011) contains the correlation coefficients associated with pairs of
measured γ-ray emission probabilities such that:
Blaton-Albicka et al.
(1976)
4.2. Uncertainty correlations of γ-ray emission probabilities
0.19(4)
0.30(8)
30.8(15)
14.3(8)
–
0.173(15)
3.68(23)
5.04(30)
25.6(23)
0.053(15)
0.128(15)
0.083(10)
0.023(5)
0.79(8)
0.079(15)
0.005(2)
0.048(8)
9.1
0.10(5)
–
29.6(20)
13.7(10)
–
0.25(5)
3.7(3)
4.9(3)
24.1(17)
0.07(2)
0.15(1)
0.08(1)
–
0.81(6)
0.08(1)
–
0.08(1)
0.054 (4)
0.048(5)
–
–
28.6(11)
11.6(7)
–
–
2.9(1)
4.5(1)
27.4(4)
–
–
0.0025(1)
–
0.011(1)
–
–
0.0006(1)
0.064(2)
0.0533(7)
0.26(5)
–
28.0(28)
14.0(14)
–
–
3.0(3)
4.0(4)
23.0(23)
0.03(1)
0.13(2)
0.08(2)
–
0.70(15)
0.08(2)
–
0.05(1)
438.8(6)
0.21(7)
–
29.9(35)
13.9(17)
–
–
3.3(4)
5.1(6)
26.4(35)
0.0347(14)
0.125(21)
0.104(14)
–
0.87(10)
0.076(14)
–
0.042(7)
χ 2 /(n − 1)
–
–
27.4(12)
14.5(14)
–
–
3.7(2)
5.2(2)
27.4(12)
0.04(2)
0.14(2)
0.13(12)
–
1.03(10)
0.11(3)
–
0.06(2)
Kossert
et al.
(2015)
0.209(3)
0.259(4)
30.45(14)
14.34(7)
0.017(3)
0.044(3)
3.782(21)
5.20(3)
26.17(13)
0.0410(16)
0.130(3)
0.0932(16)
0.0289(9)
0.845(5)
0.0799(11)
0.0070(6)
0.0474(14)
Davidson and
Connor
(1970b)
Kossert
et al.
(2015)
Energy/keV This work Briançon
et al. (1968)
χ 2 /(n − 1)
Table 7
Absolute γ-ray emission probabilities per 100 decays of 215Po observed in this work
and published values.
23
1.2
8.0
1.1
2.3
21
59
8.7
5.9
6.8
2.2
0.9
431
5.3
2600
1.0
1.8
180
S.M. Collins et al. / Applied Radiation and Isotopes 102 (2015) 15–28
24
S.M. Collins et al. / Applied Radiation and Isotopes 102 (2015) 15–28
Table 9
Uncertainty budget for the 269.5 keV γ-ray emission of
223
Ra.
Uncertainty component
Relative uncertainty
(k ¼ 1)
Standard uncertainty of the weighted mean (includes
the full-energy peak efficiency)
Activity per unit mass
Dead time and pile-up correction
Radioactive decay correction
Geometric reproducibility
Gravimetric
Detector stability
Peak fitting
True coincidence summing
Total uncertainty
0.24%
0.33%
0.050%
0.010%
0.10%
0.10%
0.20%
0.10%
0.10%
0.50%
Table 10
Correlation coefficient matrix of the most intense γ-ray emissions (4 1%) of
223
Ra.
Energy/keV
Energy/keV
122.3
144.3
154.2
269.5
323.9
338.3
445.0
122.3
144.3
154.2
269.5
323.9
338.3
445.0
100%
85%
75%
67%
72%
73%
77%
–
100%
96%
75%
74%
75%
77%
–
–
100%
75%
72%
72%
74%
–
–
–
100%
96%
95%
85%
–
–
–
–
100%
97%
89%
–
–
–
–
–
100%
90%
–
–
–
–
–
–
100%
Table 11
Correlation coefficient matrix of the intense γ-ray emissions ( 41%) of
219
Rn.
Energy /keV
Energy/keV
271.3
401.8
271.3
401.8
100%
85%
–
100%
Table 12
Correlation coefficient matrix of the intense γ-ray emissions ( 41%) of
Fig. 8. Results of the (a) 269.5 keV and (b) 271.3 keV absolute γ-ray emissions for
the three datasets. The solid line indicates the final weighted mean result reported
and the dashed lines indicate the final uncertainty of the absolute gamma emission
probability. The error bars are composed of the statistical and efficiency uncertainties only.
211
Pb.
Energy/keV
Energy/keV
404.8
427.1
831.9
404.8
427.1
831.9
100%
95%
83%
–
100%
83%
–
–
100%
the collection of greater than 106 counts in the photopeaks,
therefore a low statistical uncertainty, of the most intense γ-ray
emissions. This has also allowed many of the relatively small γ-ray
emission probabilities to be quantified, which have only previously
been seen in γ–γ coincidence measurements (Sheline et al., 1998).
The 269.5 keV absolute γ-ray emission values determined for the
three datasets of S1B, S2B and S2L are in good agreement, shown in
Fig. 8(a), with only the statistical and FEP efficiency uncertainty included in individual dataset uncertainty. The agreement between S1B
and the results of S2B and S2L indicate a consistency in the standardisations of the solutions S1 and S2. The agreement between S2B and
S2L indicates that there is no significant difference between the FEP
efficiency calibrations of the two detectors. These two observations
Fig. 9. The peak fits of the convoluted photopeaks of the 269.5 keV and 271.3 keV
γ-ray emissions.
S.M. Collins et al. / Applied Radiation and Isotopes 102 (2015) 15–28
25
Fig. 10. The reduced χ2 of the published literature for the most intense γ-rays
normalised to the 154.2 keV (hollow squares) and the 269.5 keV (hollow triangles)
γ-ray emissions.
Fig. 12. Difference and z-score of the γ-ray emission probabilities (normalised to
the 154.2 keV γ-ray) of Krien et al. (1970) and this work before (hollow squares) and
after (hollow triangles) adjustment for the 75Se nuclear data.
It is the authors' hope that the provision of the correlation data
will help to inform any future decay data evaluation of 223Ra.
While the magnitude of the correlation will to some extent be
dependent on both the choice of radionuclides for the efficiency
calibration and on the selection of the fitting function, it is suggested that by assuming the same level of correlation is present
across other measurements of the emission probability, a more
realistic evaluation, than the usually applied assumption of no
correlation between the emission probabilities, could be achieved.
5.2. Published literature of
Fig. 11. The relative FEP efficiency curves for BART determined using the relative γray emission probabilities of 75Se; using the nuclear data of Krien et al. (1970)
(hollow triangles) and Negret and Singh (2013) (hollow squares).
provide supporting evidence in the confidence of the reliability in the
presented absolute γ-ray emission probabilities.
A potential source of error could be from the de-convolution of
the 269.5 keV and 271.3 keV γ-ray emissions by the peak fitting
software. Although the use of high resolution γ-ray spectrometry
allowed the individual peaks to be relatively well resolved (see
Fig. 9) any inaccuracy in the peak shape fitting calibration could
lead to inaccuracies in the determined peak areas. Fig. 8(b) shows
the results for the 271.3 keV absolute γ-ray emission from the
219
Rn decay, showing a consistency in the results between the
datasets. It can therefore be inferred that no significant errors
occurred from the de-convolution of these two photopeaks by the
software.
The effect of the diffusion of 219Rn from the aqueous phase to
gas phase in the air space at the top of the ampoules was considered and then discounted, as this effect has been previously
investigated by Cessna and Zimmerman (2010) and Bayer (2012).
223
Ra γ-ray emission probabilities
A comparison of the values reported for the 223Ra γ-ray emissions in this work to those in the currently published literature,
which had been normalised to the 154.2 keV γ-ray emission, are
shown in Table 5. A review of the reduced χ2 for each of the main
γ-ray emissions (41%) shows that for γ-ray emissions of the
122.3 keV, 144.3 keV and 158.7 keV all the reported values are in
good agreement; with reduced χ2 values less than one. Alternatively, the reduced χ2 values for the γ-ray emissions of the
269.5 keV, 323.9 keV, 338.3 keV and 445.0 keV indicate significant
differences within the datasets. The reported values of this work
and Kossert et al. (2015) are in good agreement for these γ-ray
emissions. In the majority of cases the reported values of the remaining published literature show a significantly positive bias;
10% and 15% for the 269.5 keV and 323.9 keV γ-ray emissions
respectively.
Normalisation of the γ-ray emission to the 269.5 keV γ-ray indicates that the discrepancies in the datasets reverse; whereby the
datasets of the 122.3 keV, 144.3 keV and 158.7 keV γ-rays become
discrepant and vice-versa. This is shown in Fig. 10.
The switch in the discrepancies would suggest that the relative
ratio of the efficiency in the region of the 269.5 keV to the region
of 154.2 keV differs between that of the FEP efficiency calibration
used in this work and Kossert et al. (2015) to the earlier publications. In the majority of the previous literature it is impossible to
test the veracity of this supposition as only Krien et al. (1970) has
reported any details about the methodology of the FEP efficiency
calibration. The details provided by Krien et al. (1970) have
26
S.M. Collins et al. / Applied Radiation and Isotopes 102 (2015) 15–28
Table 13
Comparison of published experimental probabilities of the α transition to those deduced from the P(γ þ ce) balance determined from the evaluations of the DDEP (Chechev,
2011a) and the measured γ-ray emission probabilities determined in this work.
α-particle energy/MeV
5.747
5.716
5.607
5.540
5.502
5.434
Feeding level/keV
126.77
158.64
269.48
338.27
376.26
445.03 þ 446.82
This work
Chechev (2011c)
Pilger (1957)
Walen et al. (1962)
Davidson and Connor (1972)
Iα
Giannini et al.
(1962b)
Iα
Iα
Iα
Iα
Iα
10.94(24)
52.0(8)
24.3(3)
9.70(11)
0.67(2)
2.19(6)
10.0(3)
49.6(9)
25.8(6)
10.60(17)
0.74(3)
2.10(9)
10.5(–)
50.4(–)
23.6(–)
10.3(–)
0.86(–)
2.4(–)
10.2(–)
48.0(–)
25.7(–)
10.2(–)
1.3(–)
2.5(–)
8.85(18)
52.2(11)
25.3(5)
8.85(18)
0.78(–)
2.24(–)
9.50(58)
52.8(8)
24.2(4)
9.16(30)
1.00(15)
2.27(20)
allowed investigation of the difference in the FEP efficiency calibrations to be made.
5.3. Krien et al. (1970)
Krien et al. (1970) describe that the FEP efficiency calibration
was performed using a 75Se source to determine a relative efficiency over the photon energy range 66–400 keV. While 75Se is
commonly used for performing efficiency calibrations of HPGe γray spectrometers, there are significant problems with the use of
75
Se in this case, where true coincidence summing corrections can
play an important role; this can be a source of significant error.
The relative γ-ray emission probabilities of 75Se used by Krien
et al. (Krien et al., 1970) differ significantly from the evaluated
relative γ-ray emission probabilities determined by Negret and
Singh (2013). The differences in the relative γ-ray emission probabilities for the most intense γ-rays range from 0.4% to 14%.
To estimate the effect of this difference in the nuclear data, a
75
Se source was measured on BART. The relative FEP efficiencies
were determined using both sets of emission probabilities (assuming that the 125.99 keV γ-ray reported by Krien et al. (1970)
was a typographical error and should represent the 136.0 keV γray). The relative FEP efficiency calibration points were fitted with
a 3rd order and 4th order polynomial for the Krien et al. (1970) and
Negret and Singh (2013) data respectively. The respective fits of
the data points are shown in Fig. 11. It was observed that the FEP
efficiency calibration, relative to the 269.5 keV, using the Krien
et al. (1970) 75Se nuclear data was significantly different to that
using the nuclear data evaluated by Negret and Singh (2013).
Using the two FEP efficiency curves, an estimate of the effect on
the γ-ray emission probabilities reported by Krien et al. (1970) was
made, and the results are shown in Fig.12. While the estimated
corrections do not completely resolve the differences in the reported values, the z-scores (Devore, 2011) show an improvement
in the agreement between the γ-ray emission probabilities of Krien
et al. (1970) and this work.
The analysis shows that the inaccuracy of the nuclear data used
and hence the accuracy of the FEP efficiency calibration can lead to
a significant error on the accuracy of the final values. No additional
information is available at this time to investigate further the
differences in the values e.g. true-coincidence summing effects.
Hence, these corrections should only be considered as qualitative
and not as a ‘true’ correction.
5.4.
provisional rebalancing of the decay scheme was attempted to see if
the observed differences in absolute γ-ray emission probabilities
could explain this disagreement. The energy, positioning and multipolarity of the γ-ray transitions were taken from the Evaluated
Nuclear Structure Data File (ENSDF) adopted levels and γ-rays
(Browne, 2001). The total conversion coefficients were recalculated
using the BrIcc code v2.3S (Kibédi et al., 2008) based on the adopted
energies and multipolarities. For γ-rays not observed in this work,
relative intensities were also adopted from ENSDF. The compiled
data were assembled into an ENSDF file and processed with the
GABS code (Browne and Baglin, 2004). The results are presented in
Table 13. The differences between the α-emission probabilities thus
calculated and the experimental data are generally lower than
previously reported (Chechev, 2011c). It should be noted these values are purely indicative and a full decay-scheme evaluation incorporating data from this work as well as that from PTB (Kossert
et al., 2015) is merited to determine to what extent the discrepancies have been truly resolved.
6. Conclusion
A total of 83 absolute γ-ray emission probability values have
been determined experimentally for the 223Ra decay series, derived from absolute standardisations of two 223Ra solutions. The
use of two high accuracy calibrated HPGe γ-spectrometers has
allowed the precise measurement of the γ-ray emission probabilities presented, with reduced quoted uncertainties of the most
significant γ-ray emissions of the decay series. In this work significant discrepancies have been shown in the current absolute γray emission values of 223Ra decay series.
Correlation coefficient matrices have been presented for the
first time showing the correlations within the uncertainties between pairs of the most intense γ-ray emission probabilities of
223
Ra and decay progeny, to enable improved decay data evaluations in the future.
The differences in the experimental α-transition probabilities
and those deduced by the P(γ þ ce) balance using the γ-ray emission
probabilities are reported here and show some improvement in
the agreement to the published experimental values of Davidson
and Connor (1970a). There remains a paucity of precise measurements of the α-transition probabilities and further investigation of the 223Ra nuclear data is recommended to finalise the 223Ra
decay scheme.
α-transitions
As previously reported by Kellett and Nichols (2013) there is
disagreement between the P(γ þ ce) deduced α transition probabilities and those determined experimentally (Davidson and Connor, 1970a; Giannini et al., 1962b; Pilger, 1957; Walen et al., 1962). A
Disclaimer
Identification of commercial services or products does not
imply recommendation or endorsement by the National Physical
S.M. Collins et al. / Applied Radiation and Isotopes 102 (2015) 15–28
Laboratory, nor does it imply that the products or services identified are necessarily the best available for the purpose.
Acknowledgements
The authors would like to thank Andrew Fenwick, Kelley Ferreira and Lynsey Keightley for the preparation of the samples,
Dr. Stefaan Pommé for providing the figure of the decay scheme,
Mily Brewer for translation of texts and Algeta ASA (Norway) for
provision of the 223Ra solutions. This work was supported by the
National Measurement System Programmes Unit of the UK's Department for Business, Innovation, and Skills and the European
Metrology Research Programme (EMRP). The EMRP is jointly
funded by the EMRP participating countries within EURAMET and
the European Union.
References
Aitken, A.C., 1936. On least squares and linear combination of observations. Proc. R.
Soc. Edinb. 55, 42–48.
Bateman, H., 1910. The solution of a system of differential equations occurrig in the
theory of radioactive transformations. Proc. Camb. Philos. Soc. 15, 423–427.
Bayer, 2012. Radium-223 Dichloride: Bayer responses to NRC questions, Bayer.
Berger, M.J., Hubbell, J.H., Seltzer, S.M., Chang, J., Coursey, J.S., Sukumar, R., Zucker, D.S. , Olsen, K., 1998. XCOM: Photon Cross Section Database, NIST Physical
Measurement Laboratory.
BIPM, 2004. Monographie BIPM-5-Table of Radionuclides, Seven Volumes, CEA/
LNE-LNHB, 91191 Gif-sur-Yvette, France and BIPM, Pavillon de Breteuil, 92312
Sèvres, France. 〈http://www.nucleide.org/DDEP.htm〉.
BIPM, 2008. Evaluation of measurement data – Guide to the expression of uncertainty in measurement, JCGM 100:2008, Bureau International des Poids et
Mesures.
BIPM, 2011. Evaluation of measurement data – Supplement 2 to the "Guide to the
expression of uncertainty in measurement" – Extension to any number of
output quantities, JCGM 102:2011, Bureau International des Poids et Mesures.
Blaton-Albicka, K., Kotlinska-Filipek, B., Matul, M., Stryczniewicz, K., Nowicki, M.,
Ruchowska-Lukasiak, E., 1976. Precision gamma-ray spectroscopy of the decay
of 223Ra and its daughter products. Nukleonika 21, 935–947.
Briançon, C., Leang, C.F., 1968. Schema des niveaux de 219Rn. J. Phys. Colloq. 29,
184–185.
Briançon, C., Leang, C.F., Walen, R., 1968. Etude du Spectre gamma Emis par le
Radium-223 et Ses Derives. C.R. Acad. Sci. (Paris) 266B, 1533–1536.
Broda, R., Cassette, P., Kossert, K., 2007. Radionuclide metrology using liquid scintillation counting. Metrologia 44, S36–S52.
Browne, E., 2001. Nuclear Data Sheets for A ¼ 215, 219, 223, 227, 231. Nucl. Data
Sheets 93, 763–1061.
Browne, E., Baglin, C., 2004. ENSDF–Evaluations: Methodology and Worked Examples On Nuclear structure and decay data: Theory and Evaluation Manual,
Part 2, 13.
Bruland, O.S., Jonasdottir, T.J., Fisher, D.R., Larsen, R.H., 2008. Radium-223: from
radiochemical development to clinical applications in targeted cancer therapy.
Curr. Radiopharm. 1, 203–208.
Cessna, J.T., Zimmerman, B.E., 2010. Standardization of radium-223 by liquid scintillation counting. Appl. Radiat. Isot. 68, 1523–1528.
Chechev, V.P., 2011a. Table of Radionuclides (Vol. 6 – A ¼ 22 to 242), Monographie
BIPM-5, Bureau International des Poids et Mesures, pp. 79–84.
Chechev, V.P., 2011b. Table of Radionuclides (Vol. 6 – A ¼ 22 to 242), Monographie
BIPM-5, Bureau International des Poids et Mesures, pp. 85-89.
Chechev, V.P., 2011c. Table of Radionuclides (Vol. 6 – A ¼ 22 to 242), Monographie
BIPM-5, Bureau International des Poids et Mesures, pp. 125-138.
Collins, S.M., Pearce, A.K., Ferreira, K.M., Fenwick, A.J., Regan, P.H., Keightley, J.D.,
2015. Direct measurement of the Half-Life of 223Ra. Appl. Radiat. Isot. 99, 46–53.
Cothern, C.R., Connor, R.D., 1965. Studies in the active deposit of actinium. Part II:
the decay of 211Pb (AcB). Can. J.Phys. 43, 383–403.
Cox, M.G., 1972. The numerical evaluation of B-Splines. IMA J. Appl. Math. 10,
134–149.
da Silveira, E.F., de Pinho, A.G., de Barros Leite, C.V., 1971. Esquema de Niveis Do
Bi211. An. Acad. Bras. Cienc. 43, 97.
Dalmasso, J., Maria, H., 1967. Gamma radiation of 219Rn, 215Po and 211Pb. C.R. Acad.
Sci (Paris) 265B, 822–825.
Davidson, W.F., Connor, R.D., 1970a. The decay of 223Ra and its daughter products.
(I): the decay of 223Ra. Nucl. Phys. A 149, 363–384.
Davidson, W.F., Connor, R.D., 1970b. The decay of 223Ra and its daughter products.
(II): the decay of 219Rn and 215Po. Nucl. Phys. A 149, 385–391.
Davidson, W.F., Cothern, C.R., Connor, R.D., 1967. Studies in the decay of the active
deposit of actinium. III: Levels in 211Bi and its daughter products. Can. J. Phys.
45, 2295–2313.
27
Devore, J.L., 2011. Probability and Statistics for Engineering and Science, eighth ed..
Brooks/Cole, California, pp. 346–348.
Giannini, M., Prosperi, D., Sciuti, S., 1962a. Decay scheme of 211Pb. Nuovo Cim. 25,
1227–1238.
Giannini, M., Prosperi, D., Sciuti, S., 1962b. Intensity measurements of alpha groups
from 211Bi, 211Po, 219Rn and 223Ra by means of solid state counter techniques.
Nuovo Cim. 25, 1314.
Gorodetzky, S., Beck, F.A., Byrski, T., Knipper, A., 1968. Détermination de moments
angulaires et de mélanges multipolaires dans le 211Bi. Nucl. Phys. A 117,
208–230.
Hamilton, W.D., Davies, K.E., 1968. The low energy levels of 211Bi. Nucl. Phys. A 114,
577–588.
Harms, A.V., Jerome, S.M., 2004. On the integrated decay and ingrowth equations
used in the measurement of radioactive decay families: the general solution.
Appl. Radiat. Isot. 61, 367–372.
Hesselink, W.H.A., 1972. Spin Assignments in the Transitional Nuclei 215Po, 219Rn,
and 223Ra from Alpha-Gamma Directional Correlation Measurements, NP19781.
Hindi, M.M., Adelberger, E.G., Kellogg, S.E., Murakami, T., 1988. Search for the lforbidden beta decay 207Tl–4 207Pbn(570 keV). Phys. Rev. C 38, 1370–1376.
ISO, 2010. Injection Equipment for Medical Use – Part 1: Ampoules for Injectables,
ISO 9187-1:2010, International Organization for Standardization.
Kathren, R.L., 1998. NORM sources and their origins. Appl. Radiat. Isot. 49, 149–168.
Keightley, J.D., Park, T.S., 2007. Digital coincidence counting for radionuclide standardization. Metrologia 44, S32–S35.
Keightley, J.D., Pearce, A.K., Fenwick, A.J., Collins, S.M., Ferreira, K.M., Johansson, L.E.,
2015. Standardisation of 223Ra by liquid scintillation counting techniques and
comparison with secondary measurements. Appl. Radiat. Isot. 95, 114–121.
Keightley, J.D., Watt, G.C., 2002. Digital coincidence counting (DCC) and its use in
the corrections for out-of-channel events in 4(π)(β)-(γ) coincidence counting.
Appl. Radiat. Isot. 56, 205–210.
Kellett, M.A., Nichols, A.L., 2013. Library of Recommended Actinide Decay Data
International Atomic Energy Agency, Vienna, Austria.
Kibédi, T., Burrows, T.W., Trzhaskovskaya, M.B., Davidson, P.M., Nestor Jr, C.W.,
2008. Evaluation of theoretical conversion coefficients using BrIcc. Nucl. Instrum. Methods Phys. Res. Sect. A: Accel. Spectrom. Detect. Assoc. Equip. 589,
202–229.
Kondev, F.G., 2013a. Table of Radionuclides (Vol. 7 – A ¼ 14 to 245), Monographie
BIPM-5, Bureau International des Poids et Mesures, pp. 135–141.
Kondev, F.G., 2013b. Table of Radionuclides (Vol. 7 – A ¼ 14 to 245), Monographie
BIPM-5, Bureau International des Poids et Mesures, pp. 113–118.
Kossert, K., Bokeloh, K., Dersch, R., Nähle, O., 2015. Activity determination of 227Ac
and 223Ra by means of liquid scintillation counting and determination of nuclear decay data. Appl. Radiat. Isot. 95, 143–152.
Krien, K., Günther, C., Bowman, J.D., Klemme, B., 1970. Experimental investigation
of the decay 223Ra - 219Rn. Nucl.Phys. A 141, 75–96.
Kutschera, W., Ahmad, I., Armato, S.G., Friedman, A.M., Gindler, J.E., Henning, W.,
Ishii, T., Paul, M., Rehm, K.E., 1985. Spontaneous 14C emission from 223Ra. Phys.
Rev. C 32, 2036–2042.
Liang, C.F., Paris, P., Sheline, R.K., 1998. Level structure of Pb-211. Phys. Rev. C 58,
3223–3227.
Luca, A., 2010. Table of Radionuclides (Vol. 5 – A ¼ 22 to 244), Monographie BIPM5, Bureau International des Poids et Mesures, pp. 41–46.
Luca, A., 2011. Table of Radionuclides (Vol. 6 – A ¼ 22 to 242), Monographie BIPM5, Bureau International des Poids et Mesures, pp. 79–84.
Mead, R.O., Draper, J.E., 1965. Energy Levels of Bi-211. Phys. Rev. C 139, B9–B16.
Michalski, J., Sartor, O., Parker, C., Shan, M., Garcia-Vargas, J., Aksnes, A., Vogelzang,
N.J., 2013. Radium-223 dichloride (Ra-223) impact on skeletal-related events,
external beam radiation therapy (EBRT), and pain in patients with castrationresistant prostate cancer (CRPC) with bone metastases: updated results from
the phase 3 ALSYMPCA trial. Int. J. Radiat. Oncol. Biol. Phys. 87, S108–S109.
Nichols, A.L., 2011. Table of Radionuclides (Vol. 6 – A ¼ 22 to 242), Monographie
BIPM-5, Bureau International des Poids et Mesures, pp. 95–103.
Negret, A., Singh, B., 2013. Nuclear Data Sheets for A ¼75. Nucl. Data Sheets 114,
841–1040.
Nilsson, S., Franzén, L., Parker, C., Tyrrell, C., Blom, R., Tennvall, J., Lennernäs, B.,
Petersson, U., Johannessen, D.C., Sokal, M., Pigott, K., Yachnin, J., Garkavij, M.,
Strang, P., Harmenberg, J., Bolstad, B., Bruland, Ø.S., 2007. Bone-targeted radium-223 in symptomatic, hormone-refractory prostate cancer: a randomised,
multicentre, placebo-controlled phase II study. Lancet Oncol. 8, 587–594.
Parker, C.C., Pascoe, S., Chodacki, A., O’Sullivan, J.M., Germá, J.R., O’Bryan-Tear, C.G.,
Haider, T., Hoskin, P., 2013. A randomized, double-blind, dose-finding, multicenter, phase 2 study of radium chloride (Ra 223) in patients with bone metastases and castration-resistant prostate cancer. Eur. Urol. 63, 189–197.
Parker, J.L., 1984. The Use of Calibration Standards and the Correction for Sample
Self-attenuation in Gamma-ray Nondestructive Assay, LA 10045, Los Alamos
National Laboratory, pp.159–194.
Pilger, R.C., 1957, Nuclear Decay Schemes in the Actinium Family (Thesis), UCRL3877, University of California .
Rajput, M., 2010. Random coincidence summing of gamma rays and dead time
count loss corrections in HPGe-detector based gamma ray spectrometry.
Radiochim. Acta 98, 99–107.
Rose, H.J., Jones, G.A., 1984. A new kind of natural radioactivity. Nature 307,
245–247.
Rosman, K.J.R., Taylor, P.D.P., 1998. Isotopic compositions of the elements. Pure
Appl. Chem. 70, 217–235.
28
S.M. Collins et al. / Applied Radiation and Isotopes 102 (2015) 15–28
Sargent, B.W., 1939. The β- and g-rays of actinium B and actinium C. Can. J. Res. 17a,
82–102.
Sheline, R.K., Liang, C.F., Paris, P., 1998. Configurations and level structure of Rn-219.
Phys. Rev. C 57, 104–112.
Vandenbosch, S.E., Baba, C.V.K., Christensen, P.R., Nielsen, O.B., Nordby, H., 1963.
The decay scheme of Pb211. Nucl. Phys. 41, 482–496.
Vanin, V.R., de Castro, R.M., Browne, E., 2004. Table of Radionuclides (Vol. 2 –
A ¼ 151 to 242), Monographie BIPM-5, Bureau International des Poids et Mesures, 1–19.
Walen, R.J., Nedovessov, V., Bastin-Scoffier, G., 1962. Spectrographie α de ∗223Ra
(AcX) et de ses dérivés. Nucl. Phys. 35, 232–252.