Maximum Likelihood Spectrum Decomposition for Isotope Identification and Quantification

Abstract

A spectral decomposition method has been implemented to identify and quantify isotopic source terms in high-resolution gamma-ray spectroscopy in static geometry and shielding scenarios. Monte Carlo simulations were used to build the response matrix of a shielded high-purity germanium detector monitoring an effluent stream with a Marinelli configuration. The decomposition technique was applied to a series of calibration spectra taken with the detector using a multi-nuclide standard. These results are compared with decay-corrected values from the calibration certificate. For most nuclei in the standard ( <sup xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink">241</sup> Am, <sup xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink">109</sup> Cd, <sup xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink">137</sup> Cs, and <sup xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink">60</sup> Co), the deviations from the certificate values were generally no more than 6% with a few outliers as high as 10%. For <sup xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink">57</sup> Co, the radionuclide with the lowest activity, the deviations from the standard reached as high as 25%, driven by the meager statistics in the calibration spectra. In addition, a complete treatment of error propagation for the technique is presented.