Open Access Article
Hunter B. Andrews
*a,
Zechariah B. Kitzhaber
a,
Chase C. Cobbleb,
Joanna McFarlane
b and
Katherine R. Johnson
b
aRadioisotope Science and Technology Division, Oak Ridge National Laboratory, Oak Ridge, Tennessee, USA. E-mail: andrewshb@ornl.gov
bNuclear Energy and Fuel Cycle Division, Oak Ridge National Laboratory, Oak Ridge, Tennessee, USA
First published on 13th January 2026
This study evaluated the potential of laser-induced breakdown spectroscopy (LIBS) for real-time monitoring of trace gas-phase iodine, which is an element of high significance in nuclear applications due to its long radioactive half-life (as iodine-129), volatility, and biological impact. In anticipation of iodine evolving into off-gas systems in molten salt reactor and nuclear fuel recycling applications, this research aimed to assess LIBS performance in flowing argon and helium matrices; optimize measurement parameters using a multichannel spectrometer; and perform calibrations to assess predictive capabilities and limits of detection (LODs). Experimental results successfully measured gas-phase iodine in flowing argon and helium; however, trace iodine was not detected in air. Optimal delay times were determined to be 10 µs for argon and 1 µs for helium, which are consistent with the expected shorter plasma lifetime in helium relative to argon. An emission line survey was provided with the 206.16, 804.37, 902.24, and 905.83 nm peaks, which were identified as the strongest emission peaks. Calibration models were successfully built in both helium and argon, achieving LODs down to 3 ppm in helium and 5 ppm in argon. The iodine emission at 905.83 nm emerged as the most robust for calibration and was subsequently applied to a time series dataset in argon. The predictive trace confirmed the feasibility of employing LIBS for continuous, online quantification of trace iodine in flowing gas systems.
Radioactive iodine (iodine-129) is of significant interest in off-gas streams due to its high volatility, high biological activity, and long half-life (1.57 × 107 years). The US Environmental Protection Agency regulations for nuclear power operations in Title 40 of the US Code of Federal Regulations (CFR), Part 190, necessitate efficient capture of krypton-85 and iodine-129 before releasing gas streams to the environment.2 This requirement has motivated a depth of research on effective abatement of iodine from gas streams using caustic scrubbing3 or solid sorbent materials.4–6 Unfortunately, options are limited for monitoring the iodine concentration in the gas phase of these experiments; for further details the readers are directed to a recent review on iodine analytical techniques from Riley et al.7 These options are limited even further when considering monitoring off-gas streams in situ at industrial MSR or reprocessing facilities.
Previous work on real-time monitoring of molten salt off-gas streams has focused on using optical spectroscopy techniques.8 Several studies from Kireev et al. have utilized laser-induced fluorescence to detect ultra-trace concentrations of gas-phase iodine and to be able to distinguish different isotopologues.9–11 Laser-induced fluorescence provided high fidelity molecular concentration information, but depending on the analyte may need specific wavelength lasers resulting in a very targeted analytical system. Raman spectroscopy also offers insight into the molecular form of the species in the off-gas stream (e.g., molecular iodine versus iodine monochloride) through measuring vibrational bands, which could be beneficial for understanding the chemistry of the off-gas. Felmy et al. demonstrated the ability to monitor gas-phase molecular iodine (I2) using Raman and fluorescence spectroscopies in a static gas cell.12 In that study, 532 and 671 nm excitation wavelengths were investigated for Raman measurements, along with chemometric modeling. These standard Raman wavelengths also resulted in iodine fluorescence, which when used with chemometric modeling resulted in limits of detection (LODs) down to 6.7 × 10−5 torr. A complementary technique to Raman is laser-induced breakdown spectroscopy (LIBS). LIBS is performed by focusing a pulsed laser into a sample stream to generate a microplasma.13 The optical emissions of this plasma are then measured to provide an elemental spectrum regardless of molecular form. Thus, tandem Raman/LIBS measurement modalities provide an understanding of molecular speciation or concentrations and total elemental composition.
Although LIBS is typically referenced as being simultaneously sensitive to nearly all elements, its application to halogens is diminished. This lower detection is caused by challenges associated with the excitation of halogen elements related to their ionization energies.14 Furthermore, the strongest emission peaks for most halogens fall in the deep ultraviolet (UV) wavelengths that are difficult to measure outside of reduced pressures and vacuum.15,16 To the authors' knowledge, only a few studies on the measurement of gas-phase iodine have been performed. Zhang et al. investigated the quantification of iodine using LIBS in a static low-pressure cell (700 Pa) by comparing nanosecond and picosecond lasers and using the deep-UV 183 nm iodine emission.17 The result was an excellent LOD for trace iodine at 60 ppb in nitrogen bulk gas. A separate study by Zhang et al. measured the signatures of volatile iodomethane. Many iodine emissions were detected, including peaks from singly ionized iodine; however, this gas stream was more concentrated in iodine than monitoring conditions would be.18 Poole and Hawari recently demonstrated the possibility of monitoring various concentrations of iodine in an argon bulk gas.19 Their work was performed by flushing argon through a Parr reactor containing volatilized iodine and directing the flow through a gas LIBS cell. One approach typically used for measuring halogens using LIBS is the formation of diatomic molecules in the plasma, which then become excited and are more easily detected. Gaft et al. successfully demonstrated measuring iodine through the formation of calcium iodide and barium iodide in the plasma plume; however, these were solid samples in which adding calcium is simple, whereas additions to flowing gas systems are complicated and may have downstream consequences associated with the added elements.20
Previous work at Oak Ridge National Laboratory established a mobile gas LIBS platform for monitoring various laboratory systems and demonstrated its use for monitoring noble gases and molten salt aerosols, which are expected in an MSR off-gas stream.21–23 These studies used a compact multichannel spectrometer to assess LODs for a robust and easily deployable system. The present study was performed to (1) investigate the feasibility of trace gas-phase iodine detection using LIBS in flowing air, argon, and He; (2) optimize LIBS settings for trace gas-phase iodine detection using LIBS; and (3) perform quantitative calibrations in various flowing bulk gases to assess LODs for online monitoring applications using the same compact spectrometer. The ability to monitor iodine in tests representative of MSRs or reprocessing systems represents an important step to meet monitoring needs to aid facilities in adhering to 40 CFR Part 190 requirements.
![]() | (1) |
LIBS univariate models were constructed by regressing peak areas against concentration. Peak areas were calculated using Simpson integration between two selected wavelength values. Calibration models were evaluated based on their coefficient of regression (R2) and root mean square error of the calibration (RMSEC). The R2 represents the variance of the peak intensities that is explained by the concentration variable in the regression model. The RMSEC, which represents the prediction capabilities, is defined as
![]() | (2) |
i is the model-predicted concentration, and n is the total number of samples. Because RMSEC metrics are in units of concentration, it is often useful to convert RMSEC values into percent RMSEC (% RMSEC) by normalizing the values to the midpoint concentration used in the calibration. In previous online monitoring studies, model efficacy was ranked using % RMSEC values, classifying them as strong (% RMSEC ≤ 5%), satisfactory (5% < % RMSEC ≤ 10%), or indicative (10% < % RMSEC ≤ 15%).25
The LODs and limits of quantification (LOQs) were determined using the 95% prediction intervals of the linear model, as described by Mermet.26 In the present study, the LOD was determined using the horizontal intersection of the fitted regression line and the upper prediction band value at a concentration of zero. The LOQ was determined similarly but using the horizontal intersection of the upper prediction band value at a concentration of zero and the lower prediction band. The benefit of this approach is that the uncertainty from the entire regression is considered rather than just the sensitivity factor and variance in the blank signal.
| Ion | λ (nm) | Relative intensity | Relative intensity | A (108 s−1) | Elow (eV) | Ehigh (eV) | glow | ghigh |
|---|---|---|---|---|---|---|---|---|
| (Argon)a | (Helium)a | |||||||
| a Relative intensities for peaks with interferences could not be determined. Relative intensities were measured on a multichannel spectrometer; thus, differing wavelength regions have differing efficiencies. Properties of the transitions (A, E, g) are reported from the NIST Atomic Spectra Database.16 | ||||||||
| I | 206.16 | 0.02 | 0.01 | 0.13 | 0.942 649 |
6.954 610 |
2 | 4 |
| I | 804.37 | 0.11 | 0.11 | 0.06 | 6.773 664 |
8.314 624 |
6 | 4 |
| I | 902.24 | 0.05 | 0.03 | 0.22 | 7.550 151 |
8.923 960 |
2 | 4 |
| I | 905.83 | 1.00 | 1.00 | 0.32 | 6.773 664 |
8.142 020 |
6 | 8 |
| I | 911.89 | — | 0.02 | 0.25 | 6.954 610 |
8.314 624 |
4 | 4 |
| I | 912.80 | — | 0.02 | 0.07 | 8.226 963 |
9.584 865 |
4 | 6 |
![]() | ||
| Fig. 2 Overlaid LIBS spectra of various concentrations of trace iodine in (a) helium and (b) argon, highlighting the strongest iodine emission peak window. | ||
Although iodine emission peaks were detected in helium and argon, none were found in air at the trace concentrations evaluated (i.e., ≤100 ppm). While further research dedicated to exploring this behavior is needed there are likely several variables involved. These may include reduced UV light transmission due to atmospheric absorbance, reduced excitation of iodine due to the rapid quenching in air plasmas, or it may be related to the excitation process itself in monoatomic inert gases versus the polyatomic molecules in air. Future studies examining the plasma properties of these different bulk gases containing iodine, as well as the spatiotemporal evolution of the plasma would provide better insight into the excitation mechanisms. Regardless of the cause, calibration in air was not performed due to a lack of iodine signal. An overlay of air, helium, and argon spectra containing nominally 70 ppm iodine at wavelength regions of interest is provided in Fig. S5. Fortunately, the detection of trace iodine in inert gases is still highly relevant to nuclear processes such as MSR off-gas systems and molten salt reprocessing of used nuclear fuel.
![]() | ||
| Fig. 3 Iodine 905.83 nm LIBS signal versus delay time in (a) helium and (c) argon. The calculated SNRs are shown in (b) helium and (d) argon. | ||
LIBS plasmas in helium and argon gas streams exhibit very different behavior because of their differences in ionization energies. Previous work explored the differences in laser-induced plasmas formed in the two bulk gases. These studies found that (1) the use of argon results in plasma temperatures nearly double those of helium, (2) electron densities in argon are typically an order of magnitude greater than those in He, and (3) argon plasma lifetimes are much greater than those of helium plasma.21 These findings are reflected in the optimization of the iodine 905.83 nm signal shown in Fig. 3. The emission intensities of iodine peaks in argon were much greater than those measured in helium. This result is reflected in the maximum SNR in argon being nearly double what was measured in helium. Furthermore, the SNR fell much faster in helium compared to argon. These differences between the two bulk gases manifest in the optimal delay times being very different: 1 µs in helium and 10 µs in argon.
The helium dataset was processed as follows: (1) normalized to the 728.2 nm helium peak to account for variation in laser energy/plasma formation; (2) smoothed using a three-point window, first-order Savitzky–Golay filter to reduce noise; and (3) summed into 500-shot accumulations. Using 500-shot accumulations was selected to balance noise reduction and time resolution (25 s). The calibration models for iodine peaks in helium are shown in Fig. 4. For the weaker emission peaks (206.16 and 902.24 nm), the measurement precision decreases as the concentration falls. These peaks also have higher LODs compared to the stronger emissions (Table 2). The 804.37 nm peak shows excellent precision and linearity. Finally, the strongest emission peak, 905.83 nm, shows excellent precision; however, the calibration model shows nonlinear behavior, tilting around 80 ppm. Thus, this peak was fit with a second-order fit to properly consider its prediction band when assessing its associated LOD.
| I 206.16 nm | I 804.37 nm | I 902.24 nm | I 905.83 nm | |
|---|---|---|---|---|
| Helium | ||||
| R2 | 0.9840 | 0.9932 | 0.9857 | 0.9952 |
| LOD | 8 | 5 | 8 | 3 |
| LOQ | 16 | 10 | 15 | 6 |
| RMSEC | 2.65 | 1.72 | 2.5 | 0.9 |
| % RMSEC | 8.4% | 5.5% | 7.9% | 2.9% |
![]() |
||||
| Argon | ||||
| R2 | 0.9893 | 0.9921 | 0.9938 | 0.9959 |
| LOD | 8 | 7 | 6 | 5 |
| LOQ | 16 | 13 | 12 | 10 |
| RMSEC | 2.48 | 2.13 | 1.89 | 1.53 |
| % RMSEC | 6.4% | 5.5% | 4.8% | 3.9% |
The argon dataset was processed in a similar fashion: (1) normalized to the 912.30 nm argon peak to account for variation in laser energy/plasma formation; (2) smoothed using a two-point window, first-order Savitzky–Golay filter to reduce noise; and (3) summed into 500-shot accumulations. The calibration models for iodine peaks in argon are shown in Fig. 5. Overall, the measurement precision in the argon bulk gas measurements was reduced compared to those in the helium bulk gas; however, the LODs for the helium calibrations were generally on par with their argon counterparts. The slight increase in the LOD for the 804.37 nm iodine peak can be attributed to the interference with the argon peak shoulder (Fig. S4). The other major difference between the two bulk gases is the behavior of the 905.83 nm iodine peak. In argon, the nonlinear behavior was not seen in the concentration regime that was studied. This different behavior is likely associated with the different plasma temperatures, electron density, and plasma self-absorption in the two gases.
The figures of merit for iodine calibrations in helium and argon are summarized in Table 2. All models resulted in strong linear fits with R2 > 0.95. For all emission peaks in both bulk gases, the LODs, determined from their prediction bands, fall below 10 ppm. Beyond detection, the LOQs fall within 10–16 ppm; this range indicates the potential for LIBS to monitor trace gas-phase iodine without overly specialized spectrometers or low-pressure gas cells. The % RMSEC values were compared to determine which peak would best be used for prediction. For helium, the % RMSEC values of the 206.16, 804.37, and 902.24 nm peaks fall into the satisfactory category, indicating they would be useful for tracking changes but not optimal for accurate quantification. The 905.83 nm peak % RMSEC value falls below 5% in helium, signifying it would be the best emission for analytical measurements, although the second-order behavior of this peak would prevent the model from being extrapolated beyond the range of the calibration sample concentrations. In argon, the 206.16 and 902.24 nm models show slightly improved % RMSEC values due to their increased intensities (Table 1). Again, the 905.83 nm peak exhibits strong predictive capabilities with a % RMSEC below 5%; however, in argon, this model is linear, which would allow for cautious extrapolation.
As further demonstration that LIBS can be used for real-time monitoring of iodine in gas streams, the iodine 905.83 nm model was applied to the entire data series collected during the argon calibration. The concentration predictions for this test are shown in Fig. 6. The highlighted green regions indicate where the spectral signals stabilized and were used for the calibration model. This data series encompasses >1 hour of nearly continuous monitoring and transient periods when the iodine concentration was varied.
In the absence of a secondary quantification technique, a gas mixture model was used to provide comparative concentration transient profiles. Using the concentration change between 25–35 min, the time constant for the gas system was estimated to be nominally 1 min. This time constant was then used to model the mixing gas concentration using the known concentration steps. As seen in Fig. 6, the LIBS calibration model clearly resolves the changes in flow rate ratios caused by concentration changes. The prediction precision also appears stable, corresponding to the strong predictive performance, as indicated by the % RMSEC value shown in Table 2. However, the LIBS model does begin to underpredict iodine concentrations at low concentrations (<20 ppm) as the LOD is approached. Finally, because this demonstration was performed only using the near-infrared channel in the multichannel spectrometer, which also covers the wavelengths needed for noble gas monitoring, the potential exists for a significantly reduced form factor for a LIBS monitoring system tailored for noble gases and iodine.21
Gas-phase iodine was successfully measured in argon and helium matrices, although it was not detected at trace concentrations in air. This lack of detection in air requires further investigation as it likely relates to the excitation mechanisms of the iodine in air versus inert gases. Additionally, an emission line survey was reported to serve as a reference for future gas-phase iodine studies. LIBS plasmas vary greatly depending on the bulk gas used, requiring different collection settings in each gas. The optimal delay times were determined to be 1 and 10 µs in helium and argon, respectively. This result corresponds to the expected plasma lifetime of helium being short and argon being long. Following this optimization, calibration models were developed in each gas with calculated LODs as low as 3 ppm in helium and 5 ppm in argon. Finally, the strongest model (iodine 905.83 nm) was applied to a time series dataset collected in argon. The resulting prediction trace demonstrated the feasibility of real-time quantification of trace iodine in flowing inert gas systems, particularly relevant to molten salt applications.
Although this study successfully demonstrated the potential for real-time monitoring of trace iodine in inert flowing gases, several future improvements will need to be explored and further considered before field deployment. Specifically, the spectral behavior of iodine in sample streams containing a more complex composition (e.g., noble gases, aerosols) should be explored to understand matrix effects and potential suppression of iodine emission peaks from other species. Additionally, alternate equipment may result in enhanced sensitivity. For example, inert gas-purged UV spectrometers may enhance the intensity of the 206.16 nm peak, which is theoretically the strongest iodine emission peak between 200 and 1000 nm. Furthermore, spectrometers with intensified charge-coupled device detectors have been demonstrated to enhance the sensitivity of LIBS measurements (particularly halides), although this benefit comes at the sacrifice of spectral coverage because these spectrometers can typically only monitor reduced wavelength ranges. For cases where very low iodine concentrations need be monitored, laser-induced fluorescence may be performed directly on the LIBS plasma to further pump transitions of interest.
Supplementary information (SI) is available. See DOI: https://doi.org/10.1039/d5ja00481k.
Footnote |
| † This manuscript has been authored by UT-Battelle, LLC, under contract DE-AC05-00OR22725 with the US Department of Energy (DOE). The US government retains and the publisher, by accepting the article for publication, acknowledges that the US government retains a nonexclusive, paid-up, irrevocable, worldwide license to publish or reproduce the published form of this manuscript, or allow others to do so, for US government purposes. DOE will provide public access to these results of federally sponsored research in accordance with the DOE Public Access Plan (https://energy.gov/downloads/doe-public-access-plan). |
| This journal is © The Royal Society of Chemistry 2026 |