Introduction

Radioactive waste sentencing in support of nuclear decommissioning relies on robust characterisation of waste streams. Quantification of difficult-to-measure radionuclides typically relies on chemical separation prior to liquid scintillation analysis to ensure that the target analyte is effectively separated from any radiometric interferences. Spectral deconvolution can be applied to correct for interferences that have co-extracted with the analyte. However, for beta emitting radionuclides, spectral deconvolution is not usually effective, particularly if the identity of the co-extracted species is not known. Robust measurement therefore relies on efficient chemical separation alone. Confirmation of the effectiveness of the separation is usually achieved through visual inspection of spectra. This paper presents an alternative, less subjective, approach for the evaluation of liquid scintillation spectra in order to confirm the absence of radiometric interferences. The technique is based on an approach originally developed for identification of radioactive contaminants in emergency screening of drinking waters and uses multiple quench parameters to determine the mean energy of a beta emitting radionuclide, irrespective of the sample’s quench [1]. The application of the technique to routine analysis is demonstrated using reference standards.

Theory

A typical beta radiation spectrum acquired using a liquid scintillation counter (LSC) shows continuous energy distribution with a nuclide-characteristic maximum energy (βmax), average energy where the intensity of the recorded spectrum is the highest (βaverage) and a gradual intensity decrease towards the low-energy part of the spectrum [2]. Examples of unquenched LSC spectra for several radionuclides using 1220 Quantulus ultra-low level liquid scintillation spectrometer are presented in Fig. 1.

Fig. 1
figure 1

Examples of radionuclide spectra obtained with 1220 Quantulus liquid scintillation spectrometer. The spectra were collected for unquenched (no quenching agent added) standards

One of the most significant obstacles which each analyst faces while performing an LSC measurement is spectral quench or a gradual degradation of the obtained spectrum due to a multitude of optical, chemical or physical reasons which leads to a reduction of counting efficiency and a shift of the observed energy of a beta decay. Examples of spectra affected by quench are presented in Fig. 2. All activity calculations based on LSC measurements must contain a counting efficiency correction to overcome problems created by quench presence. Usually, quench correction is being done using an electron movement within the analysed sample which is induced by the presence of an external gamma radiation source (Compton effect) and registering the channel number below which 99% of the generated spectral counts reside (SQPE method, Fig. 3). For TriCarb instruments a tSIE term is used to describe similar quench correction technique. Another method is to observe a shift of a chosen point of a true LSC spectrum of an analysed radionuclide, for example, the point where 50% of counts are registered above and below on a channel scale (SQPI method, Fig. 4). [3] Again, for TriCarb instruments a term SIS is used to describe such approach.

Fig. 2
figure 2

Examples of spectral shift of beta spectra due to the presence of quench in analysed samples. SQPE range: 3H (642–722), 241Pu (473–756), 55Fe (673–777), 45Ca (624–746), 99Tc (570–750), 147Pm (674–760), 63Ni (661–756), 36Cl (631–683)

Fig. 3
figure 3

Example of a gamma-ray induced LSC spectrum (Compton effect) and SQPE quench level check technique

Fig. 4
figure 4

Examples of LSC spectra of 36Cl and 99Tc with SQPI(50) marked with dotted lines. 50% of all registered counts are recorded above or below the marked position. The position of the SQPI(50) on a channel/energy scale depends on the quench level in the sample

Both SQP(E) and SQPI parameters depend on a quench level (Q) present in the sample. SQPI is also a function of energy (βmax) of the radionuclide which is being measured. A combination of both parameters allows to produce a new parameter (SQPI/SQP(E)), which now only depends on the energy:

  1. (1)

    SQPE = f(Q)

  2. (2)

    SQPI = f(Q, βmax)

  3. (3)

    SQPI/SQPE = f(βmax)

Using SQPI/SQP(E) factor for single nuclide analysis enables identification of an unknown beta-emitting radionuclide by means of its βmax determination [3]

Interestingly, a ratio of SQPI factors determined at different spectral points e.g. SQPI(50) and SQPI(95) should remain constant throughout a wide quench level range. By analysing these parameters for a range of radionuclides it is possible to create a set of “peak shape factors” characteristic for each individual radionuclide. The presence of any other radionuclide in analysed purified samples will affect this “peak shape factor” by shifting the SQPI(95)/SQPI(50) away from the value determined for a clean, monoisotopic sample. Analysis of such ratio for large number of numerically recorded spectra allows a quick and easy way of verifying the purity of analysed samples and the validity of produced results.

By analysis of SQPI/SQPE ratio it is also possible to look for deviations in maximum beta radiation energy for the registered spectra and identify the samples which are outliers within larger analytical datasets.

Combination of both “peak shape factor” and SQPI/SQPE checks for each registered spectrum and building them into a simple software or calculation spreadsheet should become a useful tool in quality assurance diagnostics of the produced samples.

Experimental

A set of test samples spiked with a known amount of a monoisotopic standard solutions traceable to NPL, PTB or NIST and with varying levels of quench were counted using 1220 Quantulus liquid scintillation spectrometer. The method of a varying quench level standard solution preparation used for LSC spectrometer’s calibration was described in detail in a previous publication and a user manual for 1220 Quantulus LSC spectrometer [4, 5]. For each collected spectrum the SQPI(95) and SQPI(50) parameters were determined. Quench level in the measured sources was determined using a built-in SQPE check feature of a 1220 Quantulus spectrometer which utilises a 152Eu source [6]. The summary of radionuclides used in the experiment and quench levels are included in Table 1 together with the determined “peak shape factor” values.

Table 1 Summary of radionuclides used for clean standard SQPI(95)/SQPI(50) parameter

For all the radionuclides listed in Table 1 the SQPI(95)/SQPE parameter was calculated and the results plotted versus maximum beta energy. To check the behaviour of both “peak shape factor” and SQPI(95)/SQPE parameters in the presence of another radionuclide in the analysed sample, a series of 63Ni standards spiked with 3H (low-energy beta emitter), 99Tc (mid-energy beta emitter) and 36Cl (high-energy beta emitter) were produced. To check how a varying amount of a contaminating radionuclide affects the parameters of interest, 63Ni standards were spiked with varying activities of 3H and 99Tc to observe the effect of a lower and higher energy emitter presence. The influence of a gamma radiation emitter was tested using, a set of 63Ni standards spiked with varying amount of 85Sr (gamma emission at 514 keV). The effect of an alpha emitter presence was checked by spiking 63Ni standards with 242Pu (alpha emissions at: 4.901 MeV and 4.856 MeV, efficiency of ca. 76.5% and 23.5%, respectively) and 236U (alpha emissions at 4.824 MeV and 4.445 MeV, efficiency ca. 77.5% and 22.4%, respectively) solutions. For all the above spiked 63Ni standards “peak shape factor” and SQPI(95)/SQPE factors were calculated. The choice of 63Ni as a reference radionuclide in the above described tests was determined by standard availability and the fact that 63Ni is one of the most frequently analysed radionuclides in the samples originating from nuclear installations and decommissioning projects.

Results and discussion

The results of a “peak shape factor”, or SQPI(95)/SQPI(50) determination obtained for clean, monoisotopic standard solutions with varying quench levels are presented in Table 1. Results of a SQPI(95)/SQPE parameter determinations for monoisotopic sources are presented in Fig. 5.

Fig. 5
figure 5

SQPI(95)/SQPE parameters calculated for monoisotopic standards with varying amount of quench. Refer to Table 1 for SQPE range

The results collected in Table 1 show that the ratio between SQPI(95) and SQPI(50) remains reasonably constant for each tested radionuclide within a wide range of quench present in the analysed standards. Averaged values calculated for all the tested radionuclides are: 1.43 ± 0.04 (3H), 1.54 ± 0.02 (63Ni), 1.46 ± 0.02 (147Pm), 1.43 ± 0.03 (45Ca), 1.282 ± 0.006 (99Tc) and 1.22 ± 0.01 (36Cl), respectively. Uncertainties reported are standard deviation values (1σ).

Table 2 summarises the results of a clean, monoisotopic 63Ni standard analysis and after it was spiked with 3H, 99Tc and 36Cl. The values for mixed solutions are calculated based on numbers derived from a single spectrum. The values for clean sources are average values obtained for a series of spectra obtained for SQPE parameter between 661 and 756 which represent high and low level of quench for the spectrometer that was used during the experiment. Uncertainty values reported for clean 63Ni solutions are standard deviation values for the calculated averages (1σ). SQPI/SQPE factor is presented for two SQPI values: 95 and 50 to allow a comparison of a scale of the value change at different SQPI values. For the nuclide systems considered in this experiment, the SQPI(50)/SQPE parameter values are: 0.214, 0.380 and 0.625 for the addition of 3H, 99Tc and 36Cl, respectively. These values seem to differ more distinctively from a clean standard parameter of 0.245 ± 0.008 than the values obtained for SQPI(95)/SQPE parameter. In the latter case, the parameter values are 0.384, 0.671 and 0.864 while clean 63Ni signal produces value of 0.380 ± 0.020. Although data spread for SQP(95)/SQPE offer a wider spread the number calculated after 3H addition lies very close to a clean standard value which makes the distinction between “clean” and “dirty” spectrum difficult.

Table 2 SQPI(95)/SQPI(50) or “peak shape factor”, SQPI(95)/SQPE and SQPI(50)/SQPE factors for 63Ni sources both clean and mixed with a contaminating radionuclides: 3H, 99Tc, 36Cl

Images presented in Fig. 6 show the mixed (63Ni + contaminant) LSC spectra. Comparison of the visual spectrum check and the shift of the multiparametral quench factors both, as expected, show deviations from the image and values obtained for clean 63Ni standard solution.

Fig. 6
figure 6

Examples of 63Ni spectra (A ≈ 7 Bq) spiked with 3H (A ≈ 9 Bq), 99Tc (A ≈ 8 Bq) and 36Cl (A ≈ 10 Bq). For comparison a pure 63Ni spectrum is also included. Corresponding multiparametral quench parameters for the presented spectra were reported in Table 2

Varying amount of contaminant present in the sample is also results in a changing value of both “peak shape factor” and SQPI/SQPE parameters. Table 3 presents the data obtained for 63Ni solution spiked with increasing amount of 3H and 99Tc. Peak shape factor changes from 1.582 to 1.845 following the additions of increasing amount of 3H. Compared to a value obtained for clean 63Ni presented in Table 2 (1.54 ± 0.02) it becomes visible that by using “peak shape factor” analysis even relatively small amount of 3H present in 63Ni sample can be detected while visual check of a slightly contaminated spectrum may not reveal the presence of a small amount of a contaminant. Similarly, the 99Tc peak shape parameter value shifts from 2.127 to 1.805. Images illustrating the changes of 63Ni spectra with the addition of 3H and 99Tc are presented in Fig. 7. Because of the maximum beta energy difference between two components of a spectrum is higher for 63Ni–99Tc system than for 63Ni–3H system, the calculated parameters shift from clean 63Ni signal is of course much higher in case of 99Tc contamination. Such dependence of the SQPI/SQPE parameters on the energy difference between two radionuclides is presented in Fig. 8.

Table 3 Peak shape factors and SQPI/SQPE parameter value change with the increasing amount of contaminant (3H or 99Tc) present in 63Ni solution
Fig. 7
figure 7

63Ni LSC spectra with varying amount of added 3H (left) and 99Tc (right). Activities of 63Ni and contaminants are presented in Table 3

Fig. 8
figure 8

SQPI(95)/SQPE parameter value change for 63Ni standard spiked with 3H, 99Tc and 36Cl

Addition of an alpha emitter to 63Ni standard shifts SQPI(95)/SQPI(50) value to 2.910 and 2.799 for 242Pu and 236U respectively. Activities of the alpha emitters used during the experiment were approximately 10 times lower than 63Ni activity. Likewise, SQPI/SQPE parameter values for both SQPI(95) and SQPI(50) show clear shift towards values characteristic for higher maximum beta energy and the calculated results were: 0.752 and 0.680 for 242Pu (SQPI(95)) and 0.258 and 0.243 for 236U (SQPI(50)).

Addition of varying amount of 85Sr to 63Ni standard allows to check the effect of gamma radiation on the resulting spectral quench parameters. In this case, visual check of the recorded spectrum may not reveal any contamination, also SQPI/SQPE parameters values remain relatively stable. SQPI(95)/SQPI(50) values change noticeably with the increasing amount of 85Sr and change from value identical to a clean 63Ni spectrum (1.541) to 1.620. Table 4 summarises the obtained results for 85Sr-63Ni system and lack of any significant visual change of the 63Ni spectrum can be observed in Fig. 9.

Table 4 Effect of 85Sr addition on spectral parameters of 63Ni spectrum
Fig. 9
figure 9

Effect of 85Sr addition (varying amounts) on spectral parameters of 63Ni spectrum. No definite visual effect can be seen

Conclusions

Multiparametral quench analysis is a promising approach to spectrum quality check during LSC analysis. It may be especially beneficial when applied to large datasets where manual checking of spectra becomes time intensive and is susceptible to errors, misinterpretation and omittance. Accurate determination of SQPI parameters rely on a good quality spectrum which is easily attainable for high activity samples but lower activity samples will require extended counting periods in order to produce reliable calculations for SQPI. Despite this, the efficiency in which an analyst can deduce potential contamination is greatly improved by not having to manually assess spectra and by reducing the likelihood of misinterpreting data and thus saves time.

Calculating “peak shape parameters” or SQPI/SQPE values is relatively easy to achieve using specially written software or from a simple calculation spreadsheet. By observing the direction of change of SQPI/SQPE parameter, it is possible to draw conclusions about the beta energy of a contaminant present in the sample. A negative shift to lower values indicates the presence of a beta emitter with the βmax energy lower than the analysed radionuclide while a positive shift indicates that contamination is caused by a higher energy beta emitter.

By careful analysis of monoiosotopic standards with varying amount of quench, it is possible to build a database for a wide range of radionuclides and optimise the quench parameters which are most useful in purity assessment of any given radionuclide.