Next Article in Journal
Effect of Unbalanced and Non-Linear Loads on Operation of the Turbogenerator of a Distributed Generation Unit
Previous Article in Journal
Response of Seismically Damaged Steel Reduced Beam Section Joints under Fire
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Calculating the Limits of Detection in Laser-Induced Breakdown Spectroscopy: Not as Easy as It Might Seem

Applied and Laser Spectroscopy Laboratory, Institute of Chemistry of Organometallic Compounds, National Research Council, Via G. Moruzzi 1, 56124 Pisa, Italy
*
Author to whom correspondence should be addressed.
Appl. Sci. 2023, 13(6), 3642; https://doi.org/10.3390/app13063642
Submission received: 31 January 2023 / Revised: 5 March 2023 / Accepted: 9 March 2023 / Published: 13 March 2023
(This article belongs to the Section Optics and Lasers)

Abstract

:
The objectives of this paper will be to discuss the issues related to the determination of the limits of detection (LOD) in laser-induced breakdown spectroscopy (LIBS) analytical applications. The derivation of the commonly used ‘3-sigma over slope’ rule and its evolution towards the new official definition recently adopted by the International Union of Pure and Applied Chemistry (IUPAC) will be illustrated. Methods for extending the calculation of the LOD to LIBS multivariate analysis will also be discussed, using as an example the detection of Cu traces in cast iron samples by LIBS.

1. Introduction

The limit of detection (LOD) of a given element/compound is a key parameter in analytical chemistry, representing the lower concentration of a given element that an analytical technique can detect with a reasonable certainty. The limit of detection depends on the matrix in which the analyte is present but also on the specific experimental setup used for the analysis. For this reason, the LOD is often used to compare the performances of different techniques, or variations of the same technique, in terms of their analytical capability.
In this paper, we will discuss the issues related to the calculation of the LODs in applications of the Laser-Induced Breakdown Spectroscopy (LIBS) technique [1]. LIBS is a powerful spectroscopic technique which is based on the analysis of the light emitted by a laser-induced plasma (LIP). One of its most interesting features is its ability to work on untreated samples, which makes it particularly suited for in situ analysis [2,3,4,5]. In the laboratory applications of LIBS, the possibility of working without solvents or gases is also appreciated, as well as the short time needed for the analysis and the simplicity of the experimental apparatus [6,7,8,9,10].
However, the same characteristics above described make sometimes problematic the building of LIBS calibration curves of good quality. LIBS’s unique feature of not needing to prepare a sample has as a counterpart the occurrence of strong matrix effects (the signal of the analyte depends on the composition of the sample) [11], which may negatively affect its sensitivity, as well as the trueness of the LIBS analytical results.
According to the current narrative, LIBS is characterized by high LODs for many of the elements of interest (typically in the range of parts per million or higher). The LOD of a given analyte is usually derived from a calibration curve, which reports the measured signal as a function of the analyte concentration in suitable samples of known composition (standards) [12,13]. To reduce matrix effects, it is necessary to build calibration curves using standards of known composition, having the same matrix as the samples to be analyzed.
Depending on the application, the availability of suitable matrix-matched standards may be limited. It is not uncommon to see LIBS calibration curves built with less than 10 samples, sometimes characterized by concentrations of the analyte that are very different from the concentrations of the unknown sample. In some cases, the change in the concentration of the analyte among the standards may produce non-linearities in the calibration curve [14]. Finding suitable blanks (reference samples where the analyte is absent) with the same matrix of the samples to be analyzed can be problematic, too.
Besides that, it must be stressed that the growing use of multivariate analytical tools [15] for the quantitative determination of samples’ composition calls for a reconsideration of the formulas traditionally used for the calculation of the LODs [16]. A similar effort is needed for extending the calculation of LODs in the framework of standardless/calibration-free analytical approaches [17,18].
Despite the abundant literature on the topic [14,16,19,20,21,22,23,24], LODs in Laser-Induced Breakdown Spectroscopy (LIBS) are often miscalculated [1], leading to over-optimistic estimations of the analytic capabilities of the technique. Moreover, as already mentioned, the calculation of the LOD is usually limited to the case of univariate analysis, after the construction of a suitable calibration curve [14]. At the best of our knowledge, the problem of calculating the elemental LODs in multivariate analysis has not been previously discussed within the framework of LIBS analysis.
The issues associated with the building of a calibration curve in LIBS have been extensively discussed in the past [14,25,26]; as we have seen, the number of available standards is usually small, and sometimes, it can be hard to obtain a suitable blank. Moreover, in many LIBS applications, the analyte concentrations of interest are of the order of the percent or higher, a condition that may produce non-negligible self-absorption effects [27,28]. These effects can be easily compensated, but the impact of self-absorption can be subtle, introducing non-linearities in the calibration curve that may not be easy to recognize.
In the next sections, we will discuss the derivation of the commonly used ‘3-sigma over slope’ [29] rule for the determination of the LOD of a given analyte from a univariate calibration curve and how this old definition, almost universally adopted in LIBS analysis, should be modified to consider the results of recent research in the field. Subsequently, we will discuss how the current definition of LOD in univariate calibration can be extended to the case of multivariate analysis. Finally, we will introduce a simple practical method for the estimation of the LOD in multivariate analysis, recently proposed by Oleneva et al. [23].

2. Materials and Methods

The practical case that we chose for exemplifying the different methods used for the determination of the LOD using both univariate and multivariate approaches is the determination of copper concentration in cast iron (Fe-C-Si alloy, with carbon concentration between 1.8% and 4% and silicon concentration around 1–3%) samples. The ten cast iron standards were provided by the Bundesanstalt für Materialforschung und -prüfung (BAM, Berlin, Germany) to the participants to the LIBS proficiency test held in 2016 at the LIBS Conference in Chamonix, France. The Cu concentration in the samples is reported in Table 1.
The samples were provided to the LIBS proficiency test participants in the form of powder or small chips. For the analysis, the samples were placed on double tape on a glass substrate.
The LIBS spectra were acquired using the micro-Modì instrument by Marwan Instruments, Pisa [30,31,32]. The instrument is equipped with a double-pulse laser, emitting two collinear pulses of about 20 mJ each in 20 ns Full-Width at Half Maximum (FWHM), delayed by 1 microsecond; the laser pulses are focused on the sample by a Zeiss Axio microscope. The laser spot on the sample had a diameter of about 50 μm; the irradiance on the sample was, thus, of about 50 GW/cm2 per pulse. The laser-induced plasma emission is collected by an optical fiber, placed at about 1 cm from the laser spot, at an angle of about 45 degrees with the vertical axis, and analyzed by a wideband Avantes Avaspec-2048-2 USB spectrometer covering the spectral range between 190 and 900 nm (0.1 nm resolution from 190 to 450 nm, 0.3 nm resolution from 450 to 900 nm, for a total of 3745 spectral points) with a delay of 1 microsecond after the second pulse. The acquisition time is 2 milliseconds (time-integrated measurement).
For each sample, ten spectra were acquired at 12 different points, for a total of 120 spectra per sample. The instrument was manually refocused on the sample surface at each point. The LIBS spectra were independently saved for each point of analysis and each sample. The dataset of 1200 spectra, each one corresponding to 3745 wavelengths, was analyzed using the Matlab® software (version R2022a). Although the laser was carefully refocused on each point of analysis, the LIBS signal was still affected by the irregularity of the sample surface, which produced large fluctuations in the intensity of the spectra. With the aim of removing the outlier spectra, as a preliminary treatment the spectra with integral intensities larger or smaller by a factor of two with respect to the integral of the average spectrum of the same sample were removed from the dataset (typically around 10% of the 120 spectra). Subsequently, each of the remaining spectra was normalized to its integral intensity, to compensate for the effect of the fluctuations of the ablated mass from shot to shot, due to the fluctuations of the laser energy and the changes in the laser-sample coupling.

3. Results

3.1. Calculation of LOD from Univariate Calibration Curve

The currently accepted operative definition of the LOD identifies this quantity as the minimum concentration of the analyte that can be detected, “controlling the risks” [22] of false positives (fluctuation of the blank signal mistaken for analyte signal) and false negatives (signal of the analyte mistaken for a fluctuation of the blank) (see Figure 1).
Assuming a Gaussian distribution of the signal of the blank and the sample, with the same variance σ, the signal mean intensity at the LOD corresponding to a given probability of false positives PFP and false negatives PFN can be written as:
I L O D = ( e r f c 1 ( P F P ) + e r f c 1 ( P F N ) ) σ
where e r f c 1 ( x ) is the inverse of the complementary error function (1—integral of the Gaussian function from –∞ to x). In most cases, for the calculation of the I L O D , the two probabilities of false positives and false negatives are assumed to be equal to a reasonably low value P. In that case:
I L O D = 2   e r f c 1 ( P )   σ
In LIBS, the calculation of the Limit of Detection of an analyte is almost universally performed, exploiting the traditional IUPAC (now obsolete) definition [13]:
L O D = 3 σ b
where σ is the standard deviation of the LIBS signal of a sample with zero concentration of the analyte (blank) and b is the slope of the calibration curve (sensitivity). The LOD defined in Equation (3) is interpreted as the minimum concentration of the analyte that can be safely distinguished from the blank, being characterized by a signal (ILOD = b × LOD) equal to three times the standard deviation of the signal of the blank σ.
Comparing definition (3) with (1), we see that the probability of false positives (assumed equal to the probability of false negatives, as in (2)) is P = e r f c ( 3 / 2 ) = 7.7 % .
We have already anticipated that the definition (3) has been modified by IUPAC. The new definition for the LOD of a given element in a univariate calibration is:
L O D = 3.3   σ y / x 1 + C ¯ 2 ( C i C ¯ ) 2 b
where:
σ y / x = ( I i I ^ i ) 2 N
where N represents the number of points in the calibration curve [22]. The term σ y / x is an estimate of the deviation of the signals Ii with respect to the value predicted by the calibration curve I ^ i at that concentration, under the hypothesis of homoscedasticity (equal variance) of the signals close to the LOD. C i are the concentrations of the standards and C ¯ is the average of these concentrations. The term C ¯ 2 ( C i C ¯ ) 2 represent the leverage for the blank sample. The higher the leverage, the higher would be the fluctuations of the intercept of the calibration curve, which reflects the uncertainty of the signals.
The new IUPAC definition differs from (4) in two aspects. Firstly, the probability of false positives/false negatives is set to 5%, instead of 7.7%. Therefore, the factor of 3 in Equation (3) becomes 3.3 in (4). Moreover, as already pointed out forty years ago by Long and Winefordner [24], for the calculation of the LOD, the average LIBS signal of the blank is assumed to be zero. Consequently, the uncertainty on the intercept of the calibration curve cannot be neglected in the calculation of the σ of the signals. In fact, if I = I L I B S a , for calculating the standard deviation of the signal, the uncertainty on the calibration curve intercept σa must be considered, too:
σ   σ L I B S 2 + σ a 2
The new IUPAC definition considers this effect, taking into account the leverage of the calibration curve at zero concentration. The concentration of the analyte at the LOD is assumed to be close enough to zero to consider the standard deviation of its signal equal to the one of the blank (hypothesis used for deriving Equation (2) from Equation (1)).
It must be stressed that the (reasonable) hypotheses leading to the new definition of the LOD can be easily made less strict, considering for example a Student t-distribution instead of a Gaussian, if the number of measurements is low or considering the possibility of having different standard deviations for the blank and for the sample at the LOD.

3.2. LOD of Cu in Cast Iron from Univariate Calibration

For determining the LOD of copper in cast iron using a univariate approach, we firstly built the corresponding calibration curve. We focused on the narrow spectral interval comprised between 324.4 nm and 328 nm, where the two most intense Cu I emission lines (peaks at 324.75 and 327.40 nm) are visible (see Figure 2).
For building a univariate calibration curve for Cu, we defined the LIBS signal as the sum of the areas of the two copper lines (red zones in Figure 3). The background (corresponding to zero analyte signal) was estimated from the gray area in Figure 3.
Plotting the LIBS signal versus the concentration of the samples, we obtain the calibration curve shown in Figure 4.
It should be noted that in LIBS, the signal should be ideally proportional to the number concentration of the analyte, which is different from the weight concentration [9]. However, in our case, the evident non-linearity of the calibration curve must be attributed to the effect of the self-absorption of the Cu line emission [27]. For the determination of the LOD of Cu in our experimental conditions, we will, thus, limit the range of the calibration curve to the linear region comprised between zero (background signal) and 0.55 w%, thus eliminating the two samples at concentration higher than 1 w% (see Figure 5).
In this range, the relation between Cu signal and concentration in weight is substantially linear and it is, therefore, possible to define the LOD of the analyte. Note that the calibration curve does not reach zero at zero copper concentration. This is also evident from Figure 1, where the Cu line at 324.75 nm seems to be observable even in the spectra of the S3 and S5 samples, which have very low copper concentrations (0.02% and 0.03%, respectively). This is probably due to the interference of some (weak) iron line emitted at a close wavelength, but does not represent a problem for the calculation of the LOD, as soon as the uncertainty on the value of the offset is taken into account.
By applying the formula reported in Equation (4), in our experimental conditions, we obtain an LOD for Cu in cast iron equal to about 0.35 w% (0.2 w% according to the old IUPAC definition, Equation (4)). However, some caution must be taken when interpreting this result; observing directly the intensity distribution of the blank (estimated as the background signal of sample with the lower Cu concentration in the calibration set) and of sample S9, whose Cu concentration is about one half of the estimated LOD, we observe that the Gaussian fittings of the two intensity distributions overlap well above the 5% confidence limit, but the two distributions of the blank and of the signal do not overlap at all (and thus, Cu at this concentration can be discriminated from the blank with 100% probability, see Figure 6).
This is because the distribution of the LIBS signal of the sample is strongly asymmetric towards the higher intensities, and thus, the approximation of having a Gaussian distribution does not apply (the standard deviation of the signal is also larger than the one of the background/blank). This effect has been studied a few years ago by Klus et al., which demonstrated that the statistical distribution of the LIBS line intensities in some cases is better described by a General Extreme Value Distribution (GEVD) curve, rather than by a Gaussian curve [33]. The asymmetry of the LIBS intensity distribution can be also due to the spectral selection procedure used for eliminating the outliers in the measurement.
At lower Cu concentrations (S5 sample, for example, having a Cu concentration of 0.03 w%), we observe, instead, a substantial overlapping of the intensity distributions. In this case, the distribution of the signal is more symmetric and can be fitted, although not perfectly, with a Gaussian curve (see Figure 7).

3.3. LOD of Cu in Cast Iron from Multivariate Calibration

In ref. [31], studying the same samples here described, the authors have demonstrated the substantial improvement in the analytical performances of the LIBS technique that can be obtained using a multivariate approach, instead of limiting the analysis to a single predictor of the concentration, as is done using a univariate approach.
Adopting the same strategy, we trained a simple Artificial Neural Network [15], using as input/predictors all the spectral points in the range previously considered (32 points between 324.5 nm and 328.0 nm). We used a single neuron in the hidden layer, with a sigmoid transfer function.
In the framework of a multivariate analysis, the equivalent of a calibration curve is expressed by the regression graph, which represents the results of the ANN vs. the nominal concentrations. The curve obtained for the Cu concentration is shown in Figure 8.
As expected, the regression curve obtained using a multivariate approach is more precise than the corresponding univariate calibration curve (Figure 5). Pushing forward this analogy, the LOD corresponding to the regression curve can be calculated in the same way as for the calibration curve (Equation (4)). In this case (pseudo-univariate estimation), the LOD for Cu improves substantially, due to the reduction of σ obtained using a multivariate approach. Applying Equation (4), we obtain a LOD = 0.15 w% for the measurement of Cu in cast iron, which should be compared with the univariate LOD = 0.35 w%.
Additionally, in this case we can do a comparison between the distribution of the predicted values for a sample near the LOD, as shown in Figure 9.
Similar to the univariate case, the assumption of a Gaussian distribution is not fulfilled for the predicted Cu concentrations. Moreover, the standard deviation of the zero signal is not equal to the one of the signal at the calculated LOD.
Therefore, the actual limit of detection of Cu can be considered substantially lower than the one calculated using Equation (4).
An alternative empirical method for determining the LOD in multivariate analysis has been recently proposed by Oleneva et al. [23]. The authors considered the Mean Relative Error (MRE) of the estimated concentration of a given sample i, defined as:
M R E i = m e a n ( c i ˜ c i ) c i
where c i ˜ is the concentration estimated by the multivariate algorithm for each repetition of the measurement on the sample with certified concentration c i . We assume to have ordered the samples according to the growing concentration of the analyte, i.e., ci+1ci. Intuitively, the MRE of the samples with analyte concentration lower than the LOD would be higher than the ones for samples with concentration above the LOD. A practical definition of the LOD can, thus, be, according to the authors, the value of ci after which the MRE fluctuations.
Δ M R E = M R E i + 1 M R E i
which stabilizes below a given threshold.
Applying this method to the ANN results for Cu, we obtain the results shown in Figure 10.
The estimated LOD (about 0.12 w%) is slightly lower than the value estimated through the pseudo-univariate method (Equation (4)); this is coherent with the failure of the approximation of Gaussian distribution and homoscedasticity of the results, which leads to an overestimation of the LOD in Equation (4).
A summary of the results obtained with the different approaches for the LOD of Cu is reported in Table 2.

4. Discussion

In the previous sections, we have discussed the motives for using the new IUPAC formula for the calculation of the LOD of an analyte with a given experimental procedure, using a univariate or multivariate approach. The use of the old IUPAC formula should, thus, be deprecated, in LIBS as in other spectro-analytical applications. When the formulas for the LOD are applied to the case of LIBS, we noticed several critical points that would partially invalidate the application of the new IUPAC definition for the LOD (as well as the old one). The failure of the Gaussian approximation for the LIBS intensity distribution and the non-homoscedasticity of the LIBS signals contrast with the hypotheses at the basis of the derivation of Equation (4). The method proposed by Oleneva et al. [23], on the other hand, did not rely on the hypothesis of Gaussianity and homoscedasticity; however, it should be noted that the new IUPAC definition of LOD is based on the idea of quantifying the risk of false positives or false negatives for measurements performed near the LOD, and the Oleneva method (which can also be easily applied to a univariate calibration curve) does not provide information about these essential parameters. The definition of ‘stabilization’ of the MRE with increasing the concentration of the standards is essential arbitrary; moreover, at high concentrations of the analyte, other effects can produce errors and corresponding fluctuations of the MRE. On the other hand, a definition based on the analysis of the MRE is interesting from an operative perspective, although it seems more related to the problem of the quantification of the analyte concentration, rather than just on its detection.
The issue of determining the LOD of an analyte by LIBS is, thus, very complex, both in the univariate and in the multivariate case. The discussion and the examples shown in this paper should suggest some care in providing figures for the LOD of the analytes under study without checking the validity of the Gaussian distribution and homoscedasticity of the LIBS signals near the LOD. In particular, non-Gaussian profiles of the intensity distribution may arise from spectral selection treatments when, as in the example reported in this paper, the outlier spectra are removed. Since this kind of treatment has become customary in the analysis of large quantities of data, their effect on the calculation of the LOD should be checked.
As a general suggestion, in a multivariate as well as univariate quantitative approach, it is always appropriate to check the histogram of the relevant LIBS signals, comparing them with the histogram of the blank signal. To guarantee that false positives and false negatives would have the same probability at the LOD is on the other hand necessary to use for the calculation of only the data with a standard deviation very close to the one of the blank. Samples characterized by high concentrations of the analyte should be excluded from the calculation because of the possible non-linear effect, but also because the standard deviation of LIBS signals increases with the square root of the signal. Moreover, points at high concentration in the calibration curve contribute to the increase in the leverage, which must be considered in the calculation of the LOD.

Author Contributions

Conceptualization, V.P.; methodology, V.P.; software, V.P.; validation, V.P., F.P. and B.C. (Beatrice Campanella); formal analysis, V.P. and B.C. (Beatrice Campanella); investigation, F.P., S.L. and B.C. (Bruno Cocciaro); resources, V.P. and S.L.; data curation, F.P., B.C. (Beatrice Campanella), S.R. and V.P.; writing—original draft preparation, V.P.; writing—review and editing, F.P., S.L., G.L., B.C. (Bruno Cocciaro), B.C. (Beatrice Campanella), S.R. and V.P.; visualization, V.P.; supervision, V.P.; project administration, V.P. and S.L.; funding acquisition, V.P., S.L., B.C. (Beatrice Campanella) and S.R. All authors have read and agreed to the published version of the manuscript.

Funding

This research was partially funded by funded by the Italian Ministry of Education, University and Research, Research Projects of National Relevance PRIN 2020, Project UNDERLANDSCAPE, grant n. 2020428LS8.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The LIBS spectra are available upon reasonable request.

Acknowledgments

The authors acknowledge the fundamental role of Ulrich Panne for organizing the proficiency test in Chamonix, France and providing the certified cast iron samples.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Palleschi, V. Laser-Induced Breakdown Spectroscopy: Principles of the Technique and Future Trends. ChemTexts 2020, 6, 18. [Google Scholar] [CrossRef]
  2. Colao, F.; Fantoni, R.; Lazic, V.; Paolini, A.; Fabbri, F.; Ori, G.G.; Marinangeli, L.; Baliva, A. Investigation of LIBS Feasibility for in Situ Planetary Exploration: An Analysis on Martian Rock Analogues. Planet Space Sci. 2004, 52, 117–123. [Google Scholar] [CrossRef]
  3. Ibano, K.; Nishijima, D.; Ueda, Y.; Doerner, R.P. LIBS Measurement of Trace Tantalum and Rhenium in Tungsten for In-Situ Diagnostic of Nuclear Transmutation. J. Nucl. Mater. 2019, 522, 324–328. [Google Scholar] [CrossRef]
  4. Girón, D.; Delgado, T.; Ruiz, J.; Cabalín, L.M.M.; Laserna, J.J.J. In-Situ Monitoring and Characterization of Airborne Solid Particles in the Hostile Environment of a Steel Industry Using Stand-off LIBS. Measurement 2018, 115, 1–10. [Google Scholar] [CrossRef]
  5. Bulajic, D.; Cristoforetti, G.; Corsi, M.; Hidalgo, M.; Legnaioli, S.; Palleschi, V.; Salvetti, A.; Tognoni, E.; Green, S.; Bates, D.; et al. Diagnostics of High-Temperature Steel Pipes in Industrial Environment by Laser-Induced Breakdown Spectroscopy Technique: The LIBSGRAIN Project. Spectrochim. Acta-Part B At. Spectrosc. 2002, 57, 1181–1192. [Google Scholar] [CrossRef]
  6. Markushin, Y.; Melikechi, N.; Marcano, O.A.; Rock, S.; Henderson, E.; Connolly, D. LIBS-Based Multi-Element Coded Assay for Ovarian Cancer Application; Achilefu, S., Raghavachari, R., Eds.; SPIE: Bellingham, WS, USA, 2009; p. 719015. [Google Scholar]
  7. Pagnotta, S.; Lezzerini, M.; Ripoll-Seguer, L.; Hidalgo, M.; Grifoni, E.; Legnaioli, S.; Lorenzetti, G.; Poggialini, F.; Palleschi, V. Micro-Laser-Induced Breakdown Spectroscopy (Micro-LIBS) Study on Ancient Roman Mortars. Appl. Spectrosc. 2017, 71, 721–727. [Google Scholar] [CrossRef] [PubMed]
  8. Sezer, B.; Bilge, G.; Boyaci, I.H. Capabilities and Limitations of LIBS in Food Analysis. TrAC Trends Anal. Chem. 2017, 97, 345–353. [Google Scholar] [CrossRef]
  9. Abdel-Salam, Z.; al Sharnoubi, J.; Harith, M.A. Qualitative Evaluation of Maternal Milk and Commercial Infant Formulas via LIBS. Talanta 2013, 115, 422–426. [Google Scholar] [CrossRef]
  10. Syvilay, D.; Bousquet, B.; Chapoulie, R.; Orange, M.; Bourdonnec, F.-X. le Advanced Statistical Analysis of LIBS Spectra for the Sourcing of Obsidian Samples. J. Anal. At. Spectrom. 2019, 34, 867–873. [Google Scholar] [CrossRef]
  11. Miziolek, A.W.; Palleschi, V.; Schechter, I. Laser Induced Breakdown Spectroscopy (LIBS): Fundamentals and Applications; Cambridge University Press: Cambridge, UK, 2006; Volume 9780521852, ISBN 9780511541261. [Google Scholar]
  12. McNaught, M.; Wilkinson, A. IUPAC. Compendium of Chemical Terminology (“Gold Book”); Blackwell Scientific Publications: Oxford, UK, 1997. [Google Scholar]
  13. IUPAC Nomenclature, Symbols, Units and Their Usage in Spectrochemical Analysis—II. Data Interpretation Analytical Chemistry Division. Spectrochim. Acta Part B At. Spectrosc. 1978, 33, 241–245. [CrossRef]
  14. Legnaioli, S.; Botto, A.; Campanella, B.; Poggialini, F.; Raneri, S.; Palleschi, V. Univariate Linear Methods. In Chemometrics and Numerical Methods in LIBS; John Wiley & Sons: Hoboken, NJ, USA, 2022; pp. 259–276. ISBN 978-1-119-75957-7. [Google Scholar]
  15. Palleschi, V. (Ed.) Chemometrics and Numerical Methods in LIBS; John Wiley: Hoboken, NJ, USA, 2022; ISBN 978-1-119-75957-7. [Google Scholar]
  16. Allegrini, F.; Olivieri, A.C. IUPAC-Consistent Approach to the Limit of Detection in Partial Least-Squares Calibration. Anal. Chem. 2014, 86, 7858–7866. [Google Scholar] [CrossRef]
  17. Tognoni, E.; Cristoforetti, G.; Legnaioli, S.; Palleschi, V. Calibration-Free Laser-Induced Breakdown Spectroscopy: State of the Art. Spectrochim. Acta Part B At. Spectrosc. 2010, 65, 1–14. [Google Scholar] [CrossRef]
  18. Tognoni, E.; Cristoforetti, G.; Legnaioli, S.; Palleschi, V.; Salvetti, A.; Mueller, M.; Panne, U.; Gornushkin, I. A Numerical Study of Expected Accuracy and Precision in Calibration-Free Laser-Induced Breakdown Spectroscopy in the Assumption of Ideal Analytical Plasma. Spectrochim. Acta Part B At. Spectrosc. 2007, 62, 1287–1302. [Google Scholar] [CrossRef]
  19. Boqué, R.; Larrechi, M.S.; Rius, F.X. Multivariate Detection Limits with Fixed Probabilities of Error. Chemom. Intell. Lab. Syst. 1999, 45, 397–408. [Google Scholar] [CrossRef]
  20. Ortiz, M.C.; Sarabia, L.A.; Sánchez, M.S. Tutorial on Evaluation of Type I and Type II Errors in Chemical Analyses: From the Analytical Detection to Authentication of Products and Process Control. Anal. Chim. Acta 2010, 674, 123–142. [Google Scholar] [CrossRef] [PubMed]
  21. Olivieri, A.C. Analytical Figures of Merit: From Univariate to Multiway Calibration. Chem. Rev. 2014, 114, 5358–5378. [Google Scholar] [CrossRef]
  22. Olivieri, A.C. Practical Guidelines for Reporting Results in Single- and Multi-Component Analytical Calibration: A Tutorial. Anal. Chim. Acta 2015, 868, 10–22. [Google Scholar] [CrossRef]
  23. Oleneva, E.; Khaydukova, M.; Ashina, J.; Yaroshenko, I.; Jahatspanian, I.; Legin, A.; Kirsanov, D. A Simple Procedure to Assess Limit of Detection for Multisensor Systems. Sensors 2019, 19, 1359. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  24. Long, G.L.; Winefordner, J.D. Limit of Detection. A Closer Look at the IUPAC Definition. Anal. Chem. 1983, 55, 712A–724A. [Google Scholar]
  25. Duponchel, L.; Bousquet, B.; Pelascini, F.; Motto-Ros, V. Should We Prefer Inverse Models in Quantitative LIBS Analysis? J. Anal. At. Spectrom. 2020, 35, 794–803. [Google Scholar] [CrossRef]
  26. Poggialini, F.; Campanella, B.; Legnaioli, S.; Raneri, S.; Palleschi, V. About the Use of Inverse Calibration in Laser-Induced Breakdown Spectroscopy Quantitative Analysis. Spectrochim. Acta Part B At. Spectrosc. 2020, 170, 105917. [Google Scholar] [CrossRef]
  27. Rezaei, F.; Cristoforetti, G.; Tognoni, E.; Legnaioli, S.; Palleschi, V.; Safi, A. A Review of the Current Analytical Approaches for Evaluating, Compensating and Exploiting Self-Absorption in Laser Induced Breakdown Spectroscopy. Spectrochim. Acta Part B At. Spectrosc. 2020, 169, 105878. [Google Scholar] [CrossRef]
  28. Safi, A.; Tavassoli, S.H.H.; Cristoforetti, G.; Tognoni, E.; Campanella, B.; Legnaioli, S.; Pagnotta, S.; Poggialini, F.; Palleschi, V. Exploiting Self-Absorption for Plasma Characterization in Laser-Induced Breakdown Spectroscopy Experiments: A Comparison of Two Recent Approaches. Anal. Chem. 2019, 91, 8595–8601. [Google Scholar] [CrossRef]
  29. Visentini, U.; Cristoforetti, G.; Legnaioli, S.; Salvetti, A.; Tognoni, E.; Palleschi, V. Accurate Measurement of Magnesium Content in Alpha-Olefins by Laser Induced Breakdown Spectroscopy (LIBS) Technique. Optoelectron. Lett. 2007, 3, 222–226. [Google Scholar] [CrossRef]
  30. Bertolini, A.; Carelli, G.; Francesconi, F.; Francesconi, M.; Marchesini, L.; Marsili, P.; Sorrentino, F.; Cristoforetti, G.; Legnaioli, S.; Palleschi, V.; et al. Modì: A New Mobile Instrument for in Situ Double-Pulse LIBS Analysis. Anal. Bioanal. Chem. 2006, 385, 240–247. [Google Scholar] [CrossRef]
  31. Safi, A.; Campanella, B.; Grifoni, E.; Legnaioli, S.; Lorenzetti, G.; Pagnotta, S.; Poggialini, F.; Ripoll-Seguer, L.; Hidalgo, M.; Palleschi, V. Multivariate Calibration in Laser-Induced Breakdown Spectroscopy Quantitative Analysis: The Dangers of a ‘Black Box’ Approach and How to Avoid Them. Spectrochim. Acta Part B At. Spectrosc. 2018, 144, 46–54. [Google Scholar] [CrossRef] [Green Version]
  32. Palleschi, A.; Palleschi, V. An Extended Kalman Filter Approach to Non-Linear Multivariate Analysis of Laser-Induced Breakdown Spectroscopy Spectra. Spectrochim. Acta Part B At. Spectrosc. 2018, 149, 271–275. [Google Scholar] [CrossRef]
  33. Klus, J.; Pořízka, P.; Prochazka, D.; Novotný, J.; Novotný, K.; Kaiser, J. Effect of Experimental Parameters and Resulting Analytical Signal Statistics in Laser-Induced Breakdown Spectroscopy. Spectrochim. Acta Part B At. Spectrosc. 2016, 126, 6–10. [Google Scholar] [CrossRef]
Figure 1. Schematic representation of the histogram of the signal from the blank (gray) and from a sample (red). The red region is the probability of false negatives (PFN), the gray region is the probability of false positives (PFP).
Figure 1. Schematic representation of the histogram of the signal from the blank (gray) and from a sample (red). The red region is the probability of false negatives (PFN), the gray region is the probability of false positives (PFP).
Applsci 13 03642 g001
Figure 2. Detail of the LIBS spectrum of the ten cast iron samples, showing the two most intense Cu lines (average spectra). The spectra are shifted on the vertical axis for making easier the comparison.
Figure 2. Detail of the LIBS spectrum of the ten cast iron samples, showing the two most intense Cu lines (average spectra). The spectra are shifted on the vertical axis for making easier the comparison.
Applsci 13 03642 g002
Figure 3. Spectral region considered for the determination of Cu in cast iron. The red areas represent the Cu signal used for univariate analysis; the gray one is the region used for estimating the background (zero) signal.
Figure 3. Spectral region considered for the determination of Cu in cast iron. The red areas represent the Cu signal used for univariate analysis; the gray one is the region used for estimating the background (zero) signal.
Applsci 13 03642 g003
Figure 4. Calibration curve for copper.
Figure 4. Calibration curve for copper.
Applsci 13 03642 g004
Figure 5. Linear region of the calibration curve for Cu. The red line corresponds to the best linear fit of the Cu LIBS signal (defined as in Figure 2) vs. Cu concentration.
Figure 5. Linear region of the calibration curve for Cu. The red line corresponds to the best linear fit of the Cu LIBS signal (defined as in Figure 2) vs. Cu concentration.
Applsci 13 03642 g005
Figure 6. Histograms of the blank (blue) and S9 sample Cu signal (Cu concentration = 0.15 w%, yellow). The red curves are the corresponding best fit Gaussian curves.
Figure 6. Histograms of the blank (blue) and S9 sample Cu signal (Cu concentration = 0.15 w%, yellow). The red curves are the corresponding best fit Gaussian curves.
Applsci 13 03642 g006
Figure 7. Histograms of the blank (blue) and S5 sample Cu signal (Cu concentration = 0.03 w%, yellow). The red curves are the corresponding best fit Gaussian curves.
Figure 7. Histograms of the blank (blue) and S5 sample Cu signal (Cu concentration = 0.03 w%, yellow). The red curves are the corresponding best fit Gaussian curves.
Applsci 13 03642 g007
Figure 8. Regression curve of copper in cast iron.
Figure 8. Regression curve of copper in cast iron.
Applsci 13 03642 g008
Figure 9. Distribution of the predicted values for zero concentration (black) and for a sample (yellow) at the LOD (Sample S9, Cu concentration = 0.15 w%).
Figure 9. Distribution of the predicted values for zero concentration (black) and for a sample (yellow) at the LOD (Sample S9, Cu concentration = 0.15 w%).
Applsci 13 03642 g009
Figure 10. Variation of the Mean Relative Error of Cu as a function of the Cu concentration.
Figure 10. Variation of the Mean Relative Error of Cu as a function of the Cu concentration.
Applsci 13 03642 g010
Table 1. Cu concentration in the samples.
Table 1. Cu concentration in the samples.
Sample Cu (w%)
S16.29
S20.1
S30.0256
S40.548
S50.0303
S60.2445
S71.231
S80.1276
S90.150
S100.150
Table 2. Value of the LOD as calculated according to the different definitions.
Table 2. Value of the LOD as calculated according to the different definitions.
DefinitionLOD (w%)
Univariate
Old IUPAC formula (Equation (3)) 0.2
New IUPAC formula (Equation (4)) 0.35
Multivariate
Pseudo-univariate IUPAC formula (Equation (4))0.15
Oleneva et al., method0.12
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Poggialini, F.; Legnaioli, S.; Campanella, B.; Cocciaro, B.; Lorenzetti, G.; Raneri, S.; Palleschi, V. Calculating the Limits of Detection in Laser-Induced Breakdown Spectroscopy: Not as Easy as It Might Seem. Appl. Sci. 2023, 13, 3642. https://doi.org/10.3390/app13063642

AMA Style

Poggialini F, Legnaioli S, Campanella B, Cocciaro B, Lorenzetti G, Raneri S, Palleschi V. Calculating the Limits of Detection in Laser-Induced Breakdown Spectroscopy: Not as Easy as It Might Seem. Applied Sciences. 2023; 13(6):3642. https://doi.org/10.3390/app13063642

Chicago/Turabian Style

Poggialini, Francesco, Stefano Legnaioli, Beatrice Campanella, Bruno Cocciaro, Giulia Lorenzetti, Simona Raneri, and Vincenzo Palleschi. 2023. "Calculating the Limits of Detection in Laser-Induced Breakdown Spectroscopy: Not as Easy as It Might Seem" Applied Sciences 13, no. 6: 3642. https://doi.org/10.3390/app13063642

APA Style

Poggialini, F., Legnaioli, S., Campanella, B., Cocciaro, B., Lorenzetti, G., Raneri, S., & Palleschi, V. (2023). Calculating the Limits of Detection in Laser-Induced Breakdown Spectroscopy: Not as Easy as It Might Seem. Applied Sciences, 13(6), 3642. https://doi.org/10.3390/app13063642

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop