Kattathu Joseph
Mathew‡
Material Signatures & Isotopic Standards Group, Chemical Sciences Division, Oak Ridge National Laboratory, 1 Bethel Valley Road, Oak Ridge, TN 37831, USA. E-mail: mathewkj@ornl.gov
First published on 24th June 2025
The accuracy and precision of isotope amount ratio measurements using thermal ionization mass spectrometry (TIMS) instrumentation are described, and the measurement of Pu materials is emphasized. The mass fractionation observed for Am, Ga, Pu, and U for isotope amount ratio measurements using the total evaporation (TE) technique is compared with theoretical estimates to demonstrate the advantage of the TE methodology and to investigate systematic biases in the major isotope amount ratios of U and Pu certified reference material (CRM) standards from the U.S. provider of CRMs. The quality of the Pu isotopic data generated by TIMS instruments in an analytical laboratory is demonstrated by the application of the double ratio technique to estimate the 241Pu half-life. Analytical data on traceable Pu CRMs from the New Brunswick Laboratory (NBL), generated as part of routine measurements supporting various programs, are used for this half-life estimation. Although the 241Pu abundances in CRMs of 136, 137, 138, and 126-A are approximately 200–2000× smaller than those in the 241Pu material used in the previous Institute for Reference Materials and Measurements (IRMM) evaluation of the 241Pu half-life, the half-life value estimated in this work shows excellent agreement with the currently accepted value from the IRMM. This agreement also demonstrates the pedigree of the Pu isotopic standards from the NBL and the quality of the isotope amount ratio measurements using TIMS instrumentation. For both the major and minor Pu isotope amount ratios, this report describes the relative importance of the factors affecting the uncertainty of TIMS measurements, which are considered the gold standard in isotope ratio measurements (LA-UR-24-29199).
• Differences in the major element content and/or major isotope amount ratio.
• Differences in the elemental abundances of impurity elements.
• Small differences in the minor isotope amount ratios of the major elements.
• Chronometric ages.
• Concordance or discordance of the chronometric ages using different parent/daughter combinations.
Forensic analysis of SNMs is used to infer the intended use of the material and to understand its history, thereby enabling attribution of the source of the material. The chronometric dating of nuclear materials, whether they are intercepted (e.g., stolen fuel pellets) or collected on purpose (e.g., cotton swipes collected in operational nuclear facilities), is one of the most valuable techniques in nuclear forensic investigations.44,45 When properly implemented, chronometric methodologies can simultaneously provide the material's “age” (i.e., elapsed time since last purification of the decay products from the parent material), enrichment of the fissile isotope, actinide element concentrations, and minor isotope abundances. These characteristics are critical in establishing a material's source and processing history as well as identifying the potentially involved parties. Additionally, this information is crucial in evaluating the safeguards of SNMs and in verifying compliance with declared non-proliferation agreements by the countries involved.
Nuclear reactors produce Pu through neutron capture on 238U. Neutron irradiation times and the conditions under which the reactor is operated dictate the relative abundances of the Pu isotopes produced. The Pu isotopic composition establishes the intended use of the Pu material. Weapons-grade, fuel-grade, and reactor-grade Pu contain approximately 4–7%, 7–19%, and >19% 240Pu, respectively.8,9,20–23 In addition to offering clues about intended use, the relative abundance of the Pu isotopic distribution may offer attribution information with respect to reactor type, operational characteristics, and the feed U material. The Pu assay results are a good indicator of intended use and of the general processing capability. The Pu element content (i.e., assay) is a critical parameter for establishing an accurate age since separation of the Pu materials. The U and Pu contents in purified Pu materials change with time owing to radioactive decay of the Pu isotopes. The isotopes 234U, 235U, 236U, and 238U are produced by the α-decay of the isotopes 238Pu, 239Pu, 240Pu, and 242Pu, respectively. Quantification of the U daughter isotope and the abundances of the corresponding Pu parent isotope are essential for estimating the chronometric age of Pu materials.13,15–18,20,43 Concordant ages from different Pu/U (parent/daughter) isotope pairs—n(238Pu)/n(234U), n(239Pu)/n(235U), n(240Pu)/n(236U), and n(242Pu)/n(238U)—indicate that Pu was fully purified from U at some point in the material's history, and then U grew in without further addition of U from other sources.17 Concordant ages from multiple Pu/U isotope combinations may also exclude the scenario where multiple Pu sources have been mixed to produce the material that is being scrutinized. Mixed sources can be the result of several factors, including processing, the starting materials being mixtures of isotopically distinct materials, unintentional contamination, or the blending of materials to alter the bulk isotopic compositions. In addition to constraining the time of purification or production of the nuclear material, chronometric ages can also establish or eliminate genetic links among different nuclear material samples. The n(241Am)/n(241Pu) chronometer is another useful tool for elucidating the history of Pu materials. This chronometer often has uncertainties smaller than the corresponding U/Pu radiochronometric ages.13,15–18,20,23,31,33,43 The larger uncertainties on the U/Pu chronometers, as compared with the Am/Pu age, are primarily due to the inability to correctly perform the tailing corrections in the minor isotope amount ratios of U because of the nonavailability of standards with isotopic abundances like the isotopic distribution of the progeny U and because of the uncertainties in the tailing corrections of the Pu minor isotope amount ratios.46 The presence of 241Am in a Pu material indicates that enough time has passed for the 241Pu to have decayed to measurable quantities of 241Am. The 241Pu isotope has the shortest half-life (∼14.3 years) among all Pu isotopes. The 241Am concentration in a Pu material can provide information about the material's history. The decay relationship between 241Pu and 241Am produces an inherent chronometry pair for the determination of Pu age since separation. The n(241Am)/n(241Pu) model ages calculated for a material would correspond to the production age of the material if the same assumptions stated earlier for Pu/U chronometry pairs are valid. Chemical processing of the Pu material that may have affected U and Am separations (from the Pu material) to varying degrees of completeness results in discordance between the U/Pu and Am/Pu chronometric ages. Chronometry principles have been successfully applied to elucidate the history of U materials, and multiple parent/daughter and parent/granddaughter chronometry pairs have been developed for the interrogation of U materials.35–39,41–45 Recently, automated U/Pu separation techniques for environmental sample analysis have been developed,47 enabling the chemical processes to catch up to the advances in the mass spectrometry instrumentation.
The U.S. Department of Energy's New Brunswick Laboratory (NBL; now called the NBL Program Office, or NBL PO) and the European Commission's Institute for Reference Materials and Measurements (IRMM; now called the Joint Research Center in Geel, or JRC-Geel) in Belgium are the major providers of certified reference material (CRM) standards for analytical measurements involving SNMs. Most of the U.S. supply of SNM standards was produced by the National Bureau of Standards (NBS) using thermal ionization mass spectrometry (TIMS) instrumentation and measurement techniques available from the 1950s to 1970s. In 1981, the responsibility for the production and certification of SNM standards was transferred from NBS to the NBL. NBS transferred their stock of standard reference materials (SRMs) to the NBL in 1987. The NBL renamed these standards as CRMs to differentiate them from the nonnuclear standards certified by the National Institute of Standards and Technology (NIST). Traceable standards provided by NBL PO and JRC-Geel include the following: U assay standards (available matrices include oxide, metal, solution, UF6, and UF4); U isotopic standards (available in oxide, metal, solution, and UF6 matrices); a U trace impurity standard (only an oxide matrix is available); U ores (powder); gamma and neutron counting standards (encased U oxides); Pu assay standards (available matrices include oxide, metal, and nitrate); Pu isotopic standards (available matrices include oxide, metal, nitrate, and sulfate); a Th impurity standard (oxide); and Th ores (powder). The referenced literature48–50 provides a documented, comprehensive history of IRMM and NBL CRMs. Standards needed for nuclear forensics investigations are different from those available from NBL PO and JRC-Geel.51,52 For nuclear material control and accountability measurements a series of standards available from NBL PO and JRC-Geel are sufficient, though. The Pu isotopic standard CRMs 136, 137, and 138 are currently being recertified by NBL PO. See ESI S1† for the certified isotopic abundances and associated uncertainties for the Pu isotopic CRMs included in this investigation.
CRMs are used to establish the traceability of the analytical data to national and international measurement databases and to ensure that laboratory measurements are free from systematic biases. In nuclear material control and accountability measurements, which are used to quantify nuclear material inventories, assurance of accuracy (demonstration of the absence of biases) is of utmost importance because the material quantities have significant safeguards and safety implications. Analytical instruments underwent far-reaching improvements since the initial NBS characterization studies of SRMs. Analytical methods have also been developed to take advantage of the improved performance of the instrumentation (i.e., better detection capabilities and automatic instrument operation). The cumulative effect of the changes in the analytical methods and associated instrumentation was an order-of-magnitude improvement in the quality of isotope amount ratio measurements via mass spectrometry relative to performance levels in the initial SRM characterization measurements. The referenced literature51 provides a more detailed discussion of the need for standards with better uncertainties to support nuclear forensics investigations of U and Pu materials. In routine TIMS measurements, standards used for instrument calibration (quantifying the mass fractionation correction) often account for more than 90% of the uncertainty in the measurement results.53 The standards needed for nuclear forensics investigations, especially those for the analyte U,11 have received greater attention recently, and JRC-Geel has produced some particle standards that contain trace amounts of analytes of interest in support of nuclear forensics investigations54,55 of environmental samples.
Because of the ability to work with increasingly smaller quantities of SNMs for high-fidelity measurements that involve element contents using the isotope dilution mass spectrometry (IDMS) technique, multicollector TIMS and inductively coupled plasma mass spectrometry (ICP-MS) instruments have played a dominant role in developing and refining analytical techniques that support nuclear forensics investigations. Using quality control (QC) data associated with routine measurements of Pu in support of various programmatic customers (including nuclear forensics investigations), this study demonstrates that TIMS remains the gold standard in isotope amount ratio measurements.
The instrumental analysis followed the total evaporation (TE) methodology using the double-filament configuration, as described in previous studies.58–63 In this configuration, a filament with the sample (evaporation filament) is placed next to a filament with no sample (ionization filament) on the sample turret. In the source housing of the TIMS instrument, the ionization filament is taken to a current of approximately 5500 mA. After peak centering and the optimization of the source lens voltages using the 187Re signal, the evaporation filament current is slowly increased until an approximately 100 mV summed (239Pu + 240Pu) signal is obtained. Peak centering and source lens voltage optimization for the Pu mass range are then performed. The evaporation filament current is increased, under computer control, to maintain a summed (239Pu + 240Pu) isotopic intensity of approximately 6 V. The referenced literature63,64 provides additional details on the TE analysis of Pu. The TE method is regarded as a state-of-the-art analytical technique for Pu, U, and Am major isotope amount ratio measurements using TIMS instruments.58–65
In addition to corrections for mass fractionation effects, the U and Pu minor isotope ratios require corrections for tailing from the major isotopes. The approach of performing tailing corrections at the Pu minor isotope amount ratios using systematic biases observed in traceable standards with similar isotopic abundances to the unknown sample yields uncertainties lower than those obtained by subtracting the tailing effects using the manufacturer-specified abundance sensitivity.46 Additionally, the performance of these tailing corrections on a per turret basis captures the variability in the run conditions of the samples/standards, thereby making the minor isotope amount ratio uncertainties more realistic. Corrections based on abundance sensitivity leave an undesirable residual bias at the minor isotope ratios.46 For the data presented herein, the tailing corrections at the minor isotope ratios were performed using the approach described in detail in the referenced literature.46
Other factors that influence the realized uncertainties in TE-TIMS data are (i) whether analyte ionization occurs at low temperature or high temperature; (ii) variability of the ionization temperatures for the same filament current; (iii) filament backgrounds which influence analyses of low amounts of the analytes, and (iv) consistency in the filament loading conditions.
![]() | (1) |
The mass fractionation correction per atomic mass unit (u) from the n(235U)/n(238U) major isotope amount ratio is scaled appropriately to the other (minor) isotope amount ratios. For two isotopes evaporating from a hot surface with masses m1 and m2, the Langmuir equation (eqn (2)), which describes the vapor pressure, can be used as an approximation for the fractionation coefficient (β).
![]() | (2) |
Mass fractionation in TIMS measurements depends on the amount loaded on the filament, the chemical form of the analyte on the filament, the analysis configuration (single, double, or triple), the temperature at which ionization occurs, the presence or absence of ionization enhancers, and the rate of analyte evaporation from the filament surface. Even for an auto sequence analysis with the same preset analysis conditions before data collection, exact reproduction of all these analysis conditions from one filament to another is not possible. Therefore, mass fractionation is recognized as a systematic effect that is variable. Graham's law of effusion has also been used to describe fractionation in TIMS isotope ratio measurements.
In the following sections, the β values for the analytes Ga, U, Am, and Pu, as obtained from eqn (2), are compared with analytical data in isotope amount ratio measurements using the TE methodology. In quantifying and interpreting the significance of the systematic biases in the Pu (and U) isotope ratios of CRM standards, it is important to rule out the possibility that such biases are not caused by the variability of the mass fractionation correction factor during the instrumental analysis. Therefore, a comparison of the mass fractionation effects on the various analytes is relevant for this discussion.
Because Ga is highly volatile, the TE analysis of this analyte uses a single-filament configuration. TE analysis using the single-filament configuration with approximately 20 to 30 ng of Ga is efficient for Ga isotopic analysis using TIMS instruments. Fig. 1 shows the n(69Ga)/n(71Ga) isotope ratio profile in the TE analysis of Ga. Ratios measured in individual cycles (dotted symbols) and the sum-integrated isotope ratio (dashed line) are shown. Whereas the n(69Ga)/n(71Ga) ratio in individual cycles changed by approximately 2.23% during TE analysis, the sum-integrated ratio differed from the true (≅ certified) ratio by only approximately 0.32%. Using the atomic masses of the Ga isotopes, a mass fractionation correction of 0.71%/u is estimated from eqn (2). The mass fractionation (0.16%/u) shown in Fig. 1 is much smaller, reflecting the high volatility of Ga. The observed fractionation for Ga is significantly higher than those for actinide elements U, Am, and Pu. The high volatility of Ga makes the isotope amount ratio measurements unique in that Ga is routinely analysed using the single-filament configuration. By contrast, Am, U, and Pu are commonly analysed in the double-filament configuration. For all analytes (Ga, U, Pu, and Am) TE has been demonstrated to be the preferred analysis method independent of whether the method utilizes the single filament configuration or double filament configuration.
Among the analytes discussed herein, the fractionation correction factor for U has received the most attention in the last 50 years. In fact, the initial characterization of the Pu isotopic SRMs 946, 947, and 948 (now renamed NBL CRMs 136, 137, and 138) at NBS used gravimetric mixtures prepared from U isotopic endmembers for quantifying the fractionation. Consequently, the correctness of the certified isotope ratios of these Pu CRMs assumes that mass fractionation correction factors are similar for U and Pu isotope ratio measurements using TIMS. The IRMM has performed commendable work in the certification of a series of U materials and in improving the uncertainties realized in U isotope amount ratio measurements by taking full advantage of the capabilities of modern TIMS instrumentation. Thus, the IRMM has certified U CRMs with uncertainties that are close to the performance limits achievable with modern TIMS instrumentation (additional details on IRMM U standards are available elsewhere48). Fig. 2 shows the n(235U)/n(238U) major isotope ratio profile in the TE analysis of CRM U100. During TE analysis, although the n(235U)/n(238U) isotope amount ratio in individual cycles changed by approximately 2%, the sum-integrated ratio differed from the certified ratio by only approximately 0.042%. Notably, the deviation of the sum-integrated ratio from the certified ratio is approximately an order of magnitude smaller than that for the more volatile Ga.
As shown in Fig. 2, the heavier isotope 238U, the most abundant isotope in this isotopic CRM, is used for normalization. Fig. 3 shows that if 235U is used for normalization, then the mass fractionation correction factor changes sign, and the magnitude of the mass fractionation remains the same. This distinction is important because the isotope used for normalization can be heavier or lighter than the other isotopes being measured. Corrections for the mass fractionation effect need to be made at all isotope ratios by appropriately scaling the per- atomic mass unit fractionation.
![]() | ||
Fig. 3 The n(238U)/n(235U) isotope ratio profile during TE analysis of CRM U100 (the solid line represents the certified ratio and the dashed lines represent the 2σ uncertainty limits on the certified ratio). The same data from Fig. 2 are shown but with a different normalizing isotope. The magnitude of the mass fractionation remains the same as that in Fig. 2. However, the sign is negative, indicating that the sum-integrated ratio must be corrected upward to obtain the true ratio. |
For U, eqn (2) yields a theoretical mass fractionation correction of 0.21%/u, which is higher than the observed fractionation of 0.014%/u. Because the relative mass differences between the isotopes are similar for Am, Pu, and U, the theoretical mass fractionation corrections predicted by eqn (2) are similar for these three analytes.
Fig. 4 shows the n(241Am)/n(243Am) major isotope ratio profile in the TE analysis of Am. During TE analysis, although the n(241Am)/n(243Am) isotope ratio in individual cycles changed by approximately 1.55%, the sum-integrated ratio differed from the certified ratio by only approximately 0.028%. Eqn (2) yields a mass fractionation correction factor of 0.21%/u, which is higher than the observed fractionation of 0.014%/u (Fig. 4). The per- atomic mass unit fractionation observed at the Am isotope amount ratio is comparable to that estimated from U isotope ratio measurements.
The isotope profiles shown in Fig. 5–9 indicate that the sum-integrated isotopic ratio from TE indicates systematic biases in the n(240Pu)/n(239Pu) major isotope ratio between the different Pu CRMs. The n(240Pu)/n(239Pu) isotope ratio in individual cycles changed by 0.418%, 0.48%, and 0.529% for CRMs 137, 126-A, and 138, respectively. For CRM 128, the n(242Pu)/n(239Pu) major isotope amount ratio changed by 0.957% during TE analysis. While the changes in isotope ratios are expressed on a 1 u basis for CRMs 136, 137, 138, and 126-A, the change for CRM 128 is 3 u.
![]() | ||
Fig. 9 The n(240Pu)/n(239Pu) isotope ratio profile during TE analysis of CRM 138 (the solid line represents the certified ratio and the dashed lines represent the 2σ uncertainty limits on the certified ratio). The sum-integrated isotopic ratio (at the end of TE analysis) is above the certified ratio, contrary to what is expected for a ratio for which the lighter isotope is used as the normalizing isotope. The mass fractionation correction is significantly larger than what has been observed in any of the Pu isotopic standards. Literature data on this standard indicate that the major isotope ratio is systematically biased.64 |
The differences in the magnitude of the mass fractionation estimated from the different Pu standards represent small systematic differences (biases) between the Pu CRMs. Systematic biases of similar magnitude (between different U CRMs) were identified for the n(235U)/n(238U) major isotope ratio (as shown in Fig. 10) when the modified total evaporation (MTE) method67,68 was developed and applied to characterization efforts of U isotopic CRMs from the NBL PO. Fig. 10 shows that systematic biases of up to ±0.05% were observed in specific NBL U CRMs when all standards were normalized to CRM U500. Notably, when these NBL U (and Pu) CRMs were initially characterized, advanced analytical techniques, such as TE, were not available; TIMS instrumentation and analytical methodologies have evolved significantly in the 30 years since the initial characterization of these standards. In an excellent demonstration of the knowledge and proficiency of the NBS metrologists and the versatility of the first-generation mass spectrometry instruments used to certify these standards, the certified isotope ratios of the NBS U and Pu SRMs have stood the test of time. The improved precision and accuracy achieved by modern TIMS instrumentation enable the identification of systematic biases greater than or equal to ±0.02% in the major isotope ratios.69 This is especially true for Pu isotopic CRMs because the characterization of the Pu isotopic standards preceded most of the U CRMs characterized by NBS/NBL. Also of note is that during the initial certification of SRMs 946, 947, and 948, gravimetric mixtures produced from 235U and 238U were used to estimate the mass fractionation effects in TIMS isotope ratio measurements.
![]() | ||
Fig. 10 Relative deviations of the major isotope amount ratios (using CRM U500 as the comparator), as reported by Richter and Goldberg67 from certified ratios of U CRMs from the NBL. The solid line represents the bias-free data and the dashed lines indicate uncertainty limits typical for recently certified U standards from the IRMM and limits that can be achieved through the use of state-of-the-art analytical techniques like TE and MTE.61,68 Using CRM U005A or U050 as the comparator standard can bias the major isotope ratios outside the certified uncertainty limits of the IRMM standards. Error bars on individual CRMs indicate the within turret precision (n = 4 or more data points) measured in MTE experiments. |
Among other parameters, the mass fractionation correction factor in TIMS analyses depends on the amounts of the element loaded on the filaments, the chemical form of the element, and the current at which the element is evaporated from the filament surface. During the initial NBS characterization, the effect of these factors on the magnitude of the mass fractionation correction factor and the precision of the correction factor were investigated in great detail.50 Fractionation factors associated with both high-temperature and low-temperature evaporation from the filaments were investigated in detail during the initial NBS characterization. For the analytical data presented herein, the same amounts of the NBS/NBL Pu isotopic CRMs were analyzed following the same analytical routines for sample loading and for instrumental analyses. Therefore, similar mass fractionation corrections (per atomic mass unit) are expected from the different CRMs (136, 137, 138, 126-A, and 128), and the differences in the fractionation factors most likely represent systematic biases in the major isotope ratio of these standards.
Fig. 11 shows the relative deviations of the major isotope ratio for Pu CRMs 136, 137, and 138 from the certified ratios (published data from Romkowski et al.59 are shown). These are the first analytical data available in the literature for Pu isotopic standards analyzed using a TE methodology. A key difference to note is that whereas the isotopic data shown for the U standards in Fig. 10 have been corrected for mass fractionation using the measured n(235U)/n(238U) isotope amount ratio in CRM U500, the Pu data shown in Fig. 11 are not corrected for mass fractionation effects. Average measured isotope ratio values and standard errors for 10 independent measurements are shown in Fig. 11. Whereas CRMs 136 and 137 indicate fractionation corrections comparable to that in U010, CRM 138 indicates fractionation correction in the opposite direction because this standard is on the other side of the zero line. For the U010 data shown in Fig. 11, the major isotope ratio is expressed as n(238U)/n(235U) to make the lighter isotope the normalizing isotope comparable with the normalization for the Pu standards.
![]() | ||
Fig. 11 Relative deviations of the measured n(240Pu)/n(239Pu) isotope ratios in Pu CRMs 136, 137, and 138 (data reported by Ramkowski et al.59). This was the first demonstration that the TE technique can be used to improve the Pu major isotope ratio measurements using TIMS instruments. The solid line represents bias free data and the expanded uncertainties that correspond to the limits achieved using modern TIMS instruments (by the TE technique) for major isotope ratio measurements on Pu, U, or Am are shown for comparison. To be consistent with the Pu isotope ratio, for which the lighter isotope is used for normalization, the n(238U)/n(235U) ratio with the lighter isotope normalization is shown for U010. |
As discussed, fractionation effects have been investigated more extensively for U because of extensive certification efforts at NBS during the production of the U-series SRMs, recent certification efforts at the IRMM48 and NBL, and the development of the MTE technique.67,68 The MTE technique is superior in that it concurrently provides both major and minor U isotope ratios by simultaneously applying the TE principle and the systematic measurement of the tailing effects (close to half-masses) on either side of the minor isotope signal. This strategy allows the tailing corrections to be performed with greater precision. The MTE technique was also used in the recent certification efforts at the NBL for CRMs 115, 112-A, 125-A, and 116-A. A previous study69 provides summary data on recent U standards certified by the NBL and a discussion of the uncertainties suggested to be achievable in analytical labs using TIMS instrumentation and the MTE technique.
The n(235U)/n(238U) major isotope ratios shown in Fig. 10 are corrected for mass fractionation effects, whereas the n(240Pu)/n(239Pu) isotope ratios shown in Fig. 11 are not corrected for mass fractionation effects. If corrected for mass fractionation effects using CRM 136, then the CRM 138 data will be pushed higher to a percent relative deviation value of approximately 0.08%, which is consistent with the bias reported in the literature64 for this CRM standard.
Although control charts for the n(240Pu)/n(239Pu) major isotope ratios for all CRMs were evaluated in detail, no such evaluation was attempted for the n(238Pu)/n(239Pu), n(241Pu)/n(239Pu), and n(242Pu)/n(239Pu) minor isotope amount ratios because these ratios are affected by tailing corrections. For the Pu minor isotope ratios, an evaluation using control charts of isotopic data corrected for tailing effects using different correction methodologies for the biases resulting from tailing has been presented elsewhere.46
The n(240Pu)/n(239Pu) major isotope ratios, decay-corrected to a common date (chosen to be the same as the first measurement indicated in the dataset), were control-charted to estimate the uncertainties realized for the major isotope ratio by the TE methodology. Fig. 12 shows a plot of the major isotope amount ratio data in the high-burnup (fuel-grade) standard CRM 136. A 2σ warning limit, corresponding to a 95% confidence interval of 0.024%, was obtained for the n(240Pu)/n(239Pu) major ratio of CRM 136 using the long-term QC dataset covering approximately 5.5 years of analytical data for this standard. Two data points (representing 0.57% of the data) plotted between the 2σ and 3σ limits indicate that the uncertainty limits estimated are not overly conservative and are appropriate for the dataset.
Fig. 13 shows a plot of the major ratio data (decay-corrected to the date corresponding to the first data point in the dataset) in the high-burnup (reactor-grade) standard CRM 137. A 2σ warning limit, corresponding to a 95% confidence interval of 0.030%, was obtained for the n(240Pu)/n(239Pu) major ratio of CRM 137 using the QC dataset covering more than 5 years. The slightly larger limit for CRM 137, compared with that of CRM 136, can be attributed to the fewer number of measurements for this isotopic standard.
Fig. 14 and 15 show the major ratio data in low-burnup standards CRM 138 (fuel-grade material) and CRM 126-A (weapons-grade material). QC data from both of these low-burnup CRMs are consistent with a 2σ warning limit, corresponding to a 95% confidence interval of 0.030%. For all four CRMs, the estimated 2σ limits represent the expanded uncertainties that can be realized in a single measurement using TIMS instrumentation and a coverage factor of k = 2.
Fig. 16 shows the uncertainty budget for the n(240Pu)/n(239Pu) major isotope ratio. Precision of the isotope ratio measurement dominates the uncertainty budget: 65% of the overall uncertainty comes from this factor. The uncertainty in the major isotope ratio of the comparator standard used to quantify the fractionation correction accounts for 35% of the uncertainty budget. For the n(240Pu)/n(239Pu) isotope ratio, the uncertainty budgets for all four Pu isotopic CRMs (both high- and low-burnup Pu standards) evaluated here are similar. The uncertainty budget for the major isotope ratio obtained here is different from that estimated by Bürger et al.,53 where more than 90% of the uncertainty was estimated to come from the uncertainty in the certified major isotope ratio of the comparator standard used for estimating the mass fractionation correction factor.
The nominal uncertainty of 0.024% to 0.030% (expanded uncertainty at the 95% confidence interval) estimated for the n(240Pu)/n(239Pu) major isotope amount ratio is similar to the uncertainties assigned to the major U isotope ratios in IRMM U CRMs, demonstrating that these uncertainty limits can be achieved on a routine basis in analytical laboratories51,67,69 with excellent quality control practices. The accuracy and precision achieved in recent NBL U CRM certification measurements indicate that the uncertainty limit is also similar to the uncertainty estimate for the n(235U)/n(238U) major isotope ratio measurements.69 This uncertainty limit also implies that for Pu and U isotopic standards, systematic biases larger than 0.03% in the n(240Pu)/n(239Pu) and n(235U)/n(238U) isotope ratios are statistically significant.64 For Pu isotopic analyses, a high-intensity total evaporation (HI-TE) method72 that has recently been developed to improve the minor isotope amount ratio uncertainties yields marginally improved performance for the n(240Pu)/n(239Pu) major isotope ratio.
While the tailing corrections for the n(238Pu)/n(239Pu) minor ratio amount to approximately 1% for the high-burnup CRMs 136 and 137, corrections greater than 10% and greater than 20% are necessary for the low-burnup CRMs 126-A and 138, respectively. Because the 238Pu abundances in CRMs 126-A and 138 are approximately 10–20× smaller than those in CRMs 136 and 137, the tailing correction is a significant source of uncertainty for the n(238Pu)/n(239Pu) isotope ratio (the variability in the isotope ratio measurement remains the larger contributor to the uncertainty budget) for low-burnup Pu materials. Performing the tailing corrections on a per-turret basis captures the variability in the tailing correction factor more accurately than making the corrections using an assumed (or measured) correction factor based on manufacturer-specified abundance sensitivity.46Fig. 18 shows the uncertainty budget for the n(238Pu)/n(239Pu) minor isotope ratio using the TE methodology.
For the n(241Pu)/n(239Pu) and n(242Pu)/n(239Pu) isotope ratios, the observed biases resulting from tailing are less than 0.5%, and the uncertainty budgets for these isotope ratios are dominated by the precision of the isotope ratio measurements. The measured values of these isotope ratios, within the precision of the measurements, overlap the certified isotope ratios, suggesting that statistically significant corrections resulting from tailing are not present.46 Uncertainty budgets for n(241Pu)/n(239Pu) and n(242Pu)/n(239Pu) are similar; the precision of the isotope ratio measurement accounts for approximately 65–70% of the overall uncertainty in the isotope amount ratio. Tailing correction accounts for the remainder of the uncertainty budget in the n(241Pu)/n(239Pu) and n(242Pu)/n(239Pu) isotope amount ratios. Uncertainty in the mass fractionation correction factor is not a significant contributor to the minor Pu isotope amount ratio uncertainties.
Recently, a HI-TE method has been developed. This method further improves the uncertainty in the minor isotope amount ratios.72 In this method, approximately double the amount of Pu was loaded on the filaments (compared with the amounts mentioned in Section II for TE analysis), and a summed (239Pu + 240Pu) target intensity of 24 V was used. The minor isotope amount ratio uncertainties achieved using the HI-TE method are substantially improved (by a factor of 2 to 3) compared with those achieved using the TE technique. Fig. 19 shows the bias in the n(238Pu)/n(239Pu) minor isotope amount ratio using the HI-TE method. For HI-TE measurements, the tailing corrections (as a percent) become smaller; more importantly, the uncertainty on the tailing corrections also becomes smaller. The HI-TE method can be used to obtain state-of-the-art uncertainties for Pu minor isotope amount ratios.
In principle, the 241Pu half-life value can be estimated from changes in the n(241Pu)/n(239Pu) ratio over time.73–75 Zeigler and Ferris75 used this approach to estimate the 241Pu half-life value using data generated by six DOE laboratories (on an interlaboratory comparison standard with approximately an order-of-magnitude higher concentration of the 241Pu isotope than CRM 137, the CRM with the highest 241Pu abundances among the four isotopic CRMs included in this study) and estimated a value of 14.89 ± 0.11 years. Although the estimated half-life value from this investigation was determined using a Pu material with approximately an order of magnitude higher 241Pu abundance than the CRM with the highest 241Pu included in this study, the offset of the half-life estimate from the currently accepted value has implications for the TE method using TIMS instruments. The single isotope ratio technique that was used for the Zeigler and Ferris75 study did not yield accurate isotope ratio results because of the nature of the analytical methodology and the analytical methodology not being based on TE. For the n(241Pu)/n(239Pu) isotope ratio, eqn (2) implies a fractionation correction factor of 1.00418. Therefore, variations less than or equal to 10−3 in the mass fractionation correction factor are sufficient to increase the uncertainty in the half-life value estimated. Fig. 5 thru 9 shows that the isotope ratios changed by several percent during the course of the Pu release from the filament. In the Zeigler and Ferris75 study, these changes in the isotope amount ratios cause accuracy issues in the data. In the double ratio method, the decrease over time of an isotope is determined not by a single ratio but by using a suitably chosen double ratio. A detailed description of the double ratio technique is provided in the referenced literature.76–78
The feasibility of estimating the 241Pu half-life value by applying the double ratio technique to QC data on CRMs 136, 137, 138, and 126 has been demonstrated.79 The application of the double ratio technique to a smaller dataset than the one presented herein yielded a 241Pu half-life value of 14.349 ± 0.023 years. In addition to the single ratio and double ratio techniques, the in-growth of the 241Am daughter by IDMS has also been used to estimate the 241Pu half-life value.80
For the TIMS instrumental data presented here, the double ratio shown in eqn3 is used:
![]() | (3) |
For the double ratio defined in eqn (3), the isotopic fractionation is given by using eqn (4):
![]() | (4) |
Eqn (4) shows that for the (239Pu·241Pu)/(240Pu)2 double ratio, deviations from unity are 100× smaller than those for the single ratio, indicating negligible mass fractionation correction. Eqn (5) represents the relationship between the slope of the double ratio defined in eqn (3) and the decay constants of the three Pu isotopes:
slopedouble ratio = λ241 + λ239 − 2λ240. | (5) |
Therefore, the 241Pu decay constant can be expressed in terms of the slope of the double ratio and the decay constants of the 240Pu and 239Pu isotopes, as represented in eqn (6):
λ241 = slopedouble ratio + 2λ240 − λ239. | (6) |
The half-life τ1/2 of 241Pu can then be calculated from the decay constant using eqn (7):
![]() | (7) |
![]() | ||
Fig. 20 Plot of the (239Pu·241Pu)/(240Pu)2 ratio vs. time for CRM 136. The slope of the linear regression line through the data and the correlation coefficient of the regression line are indicated. |
For CRM 136, Fig. 21 shows a plot of the residuals from the line of best fit shown in Fig. 20. The data follow a normal distribution, as expected for a dataset with large degrees of freedom. Separations on two different days, marked by the dashed ellipsoidal, do not overlap the zero line within characteristic uncertainties of the measurements, indicating that small amounts of 241Am remain in the Pu isotopic fraction measured. The Pu isotopic separations in the last 2.5 years included an additional rinsing step to completely remove the 241Am from the Pu fraction measured for isotope ratios. The deviation from the zero line was significantly reduced (showing improvement by a factor of 2.2 in the standard deviation of the residuals) after the additional column rinse was introduced into the separation process.
![]() | ||
Fig. 21 Residuals of the individual data points from the line of best fit shown in Fig. 20 for CRM 136. A normal distribution around the “0.0” value (solid line) demonstrates the quality of the data. |
![]() | ||
Fig. 22 Plot of the (239Pu·241Pu)/(240Pu)2 ratio vs. time for CRM 137. The slope of the linear regression line through the data and the correlation coefficient of the regression line are indicated. |
Fig. 23 shows a plot of residuals from the line of best fit in Fig. 22. As was the case for CRM 136, separations on a few days do not overlap the zero line within the characteristic uncertainties of the data set, indicating that small amounts of 241Am remain in the clean Pu isotopic fraction. An additional rinsing step in the Pu isotopic purification process was introduced for the separations in the last 2.5 years to completely remove the 241Am from the Pu isotopic fraction. The deviation from the zero line was significantly reduced (showing improvement by a factor of 1.8 in the standard deviation of the residuals) after the additional column rinse was introduced into the separation process.
![]() | ||
Fig. 23 Residuals of the individual data points from the line of best fit shown in Fig. 22 for CRM 137. A normal distribution around the “0.0” value (solid line) demonstrates the quality of the data. |
![]() | ||
Fig. 24 Plot of the (239Pu·241Pu)/(240Pu)2 ratio in CRM 138 vs. time. The slope of the linear regression line through the data and the correlation coefficient of the regression line are indicated. |
![]() | ||
Fig. 25 Residuals of the individual data points from the line of best fit shown in Fig. 24 for CRM 138. A normal distribution around the “0.0” value (solid line) demonstrates the quality of the data. |
![]() | ||
Fig. 26 Plot of the (239Pu·241Pu)/(240Pu)2 ratio in CRM 126-A vs. time. The slope of the linear regression line through the data and the correlation coefficient of the regression line are indicated. |
![]() | ||
Fig. 27 Residuals of the individual data points from the line of best fit shown in Fig. 26 for CRM 126-A. A normal distribution around the “0.0” value (solid line) demonstrates the quality of the data. |
Footnotes |
† Electronic supplementary information (ESI) available. See DOI: https://doi.org/10.1039/d5ja00140d |
‡ Actinide Analytical Chemistry (C-AAC), Los Alamos National Laboratory, Los Alamos, NM 87545, USA. |
This journal is © The Royal Society of Chemistry 2025 |