Limit of detection (LOD) estimated for different elements.
Laser-induced breakdown spectroscopy (LIBS) has gained great attention over the past two decades due to its many advantages, such as needless sample preparation, capability of remote measurement and fast multielement simultaneous analysis. However, because of its inherent uncertainty features of plasma, it is still a big challenge for LIBS community worldwide to realize high sensitivity and accurate quantitative analysis. Currently, many chemometric analytical methods have been applied to LIBS calibration analysis, including univariate regression, multivariate regression, principal component regression (PCR), partial least squares regression (PLSR) and so on. In addition, appropriate sample and spectral pretreatment can effectively improve the analytical performance (i.e., limit of detection (LOD), accuracy and repeatability) of LIBS. In this chapter, we briefly summarize the progress of these calibration methods and their applications on LIBS and provide our recommendations.
- laser-induced breakdown spectroscopy
- sample and spectral pretreatment
- calibration methods
- calibration-free laser-induced breakdown spectroscopy
Laser-induced breakdown spectroscopy (LIBS), also sometimes called laser-induced plasma spectroscopy (LIPS), has developed rapidly as an analytical technique over the past two decades. LIBS is a kind of atomic emission spectroscopy, which uses a high-energy pulsed laser as the excitation source. The laser is focused on the sample surface, thereby evaporating and atomizing the sample and generating a plasma. The light emitted by the plasma is detected by a spectrometer. One can obtain sample composition and concentration information by analyzing the plasma emission spectra. LIBS most commonly used experimental instruments such as mainly lasers, spectrometers, detectors and computers, as shown in Figure 1.
The laser which is the most widely used in LIBS is the Nd:YAG solid-state lasers operated in the Q-switch mode. Typically, it is operated at the fundamental wavelength of 1064 nm, pulse energy is 30–100 mJ, pulse width is 5~15 ns and repetition rate is 1~10 Hz. In addition, the researchers tested the effects of lasers of different types and parameters on the LIBS. Trautner et al.  investigated polyethylene (PE) and a rubber material from tire production by employing 157 nm F2 laser and 532 nm Nd:YAG laser ablation in nitrogen and argon gas background or in air. The effects of laser wavelength on depth resolution of thin-film solar cell are investigated by Choi et al.  using an ultraviolet (λ = 266 nm) and a visible (λ = 532 nm) nanosecond Nd:YAG lasers. Labutin et al.  summarize nearly two decades of studies on femtosecond laser-induced breakdown spectrometry (fs-LIBS). Picosecond pulse train and nanosecond pulse were compared for laser ablation and LIBS measurements by Lednev et al. .
Spectrometers disperse the emitted radiation of the laser-induced plasma to get a spectrum in terms of intensity as a function of the wavelength. The dominant spectrometer types used for LIBS are multichannel fiber spectrometer and echelle spectrometer coupled with an intensified CCD. The echelle spectrometer offers a wide spectral range, a high spectral resolution, and the possibility of time-resolved. The plasma parameters (plasma temperature and electron density) are constantly changing with respect to the delay time, and an echelle spectrometer with time resolution is needed when calculating these parameters. However, time-resolved broadband spectrometers are expensive and strongly dependent on external circumstances. The multichannel fiber spectrometer is robust and reliable for the use in mobile and portable LIBS instruments, providing an accredited spectral resolution, but their integration time is typically much longer than the plasma lifetime.
In recent years, with the rapid development of lasers, spectrometers and detectors and the urgent demand of in situ and online analysis, LIBS has developed rapidly. Compared with many other types of elemental analysis techniques, LIBS has obvious advantages:
Simple equipment: few instruments, low cost and easy integration.
Noncontact analysis: LIBS uses pulsed laser as the excitation source, which makes it noncontact analysis, especially in the dangerous environment or space exploration field, has broad application prospects.
No sample preparation: LIBS directly focused pulsed laser bombardment of the sample without processing the sample.
Various samples: samples can be gas, aerosols, liquids and solids.
Nondestructive analysis: The laser converges to the surface of the sample, and only a small amount of the sample is excited. It can be considered as nondestructive or near nondestructive.
Three-dimensional analysis: LIBS can collect laser at different positions on the sample surface or repeat measurements at the same location to analyze the sample surface and its different depths of the sample composition and content.
Total element analysis: The laser energy can simultaneously excite all the elements in the sample, so all elements in the sample can be analyzed simultaneously.
Remote analysis: The long-distance analysis of the LIBS can be achieved by remotely transmitting the laser energy and collecting the plasma emission spectrum through the fiber.
Online analysis: LIBS is a very fast technology that provides analytical results in seconds, making it particularly suitable for rapid analysis or online industrial monitoring.
Because of many advantages of LIBS, it has been applied to a number of analytical domains, for example, various alloys [5, 6, 7], slags [8, 9], soil , rocks [11, 12] and isotopes . We searched for all scientific papers from 1963 to 2016 on Web of Science with laser-induced breakdown spectroscopy (LIBS). The statistical results are presented in Figure 2. It can be seen that LIBS has developed rapidly since 1990.
2. Pretreatment of samples and spectra
2.1. Sample pretreatment
One of the most widely cited advantages of LIBS is that it does not require sample preparation, but this may also be the biggest limitation for improving its consistancy. In general, LIBS performance may be enhanced using two main approaches: pretreatment of samples and spectra. Many homogeneous solid samples require no sample preparation, for example, glass, alloy and plastic. For powder samples (e.g., cement , soil  and coal ) which can press the cake directly, it must be consistent with the standard sample preparation process used for calibration during the pressing process. Comparing with solid samples, the direct analysis of liquid samples by using LIBS has many disadvantages: splash, less excitation and fluctuation of liquid level. The simplest way to change a liquid sample into a solid sample is to freeze it [16, 17]. Sobral et al.  investigated the detection sensitivity of Cu, Mg, Pb, Hg, Cd, Cr and Fe traces in water and ice samples under the same experimental conditions by using LIBS. Another effective way that can be used for liquid analysis in a solid matrix configuration consists of using an absorbent substrate, for example, plant fiber spunlace nonwoven , absorbent paper , thin wood sample [21, 22] and membrane-based filter paper . Now LIBS analysis of aerosols is mainly of two categories: direct analysis and enrichment. However, the detection limit and the statistical results of the direct analysis are still relatively poor. On the other hand, the substrate-based collection does not provide as instantaneous information and does allow one to achieve lower detection limits by increasing the sample flow rates and sampling times.
2.2. Correction and removal of continuum background
The detected plasma emission spectrum at a given wavelength in a spectrum is the sum of the analyte signal and the continuum background. The analyte signal is often overwhelmed by the continuum background, which interferes the true intensity of signal and compromises spectral clarity and hence reduces the accuracy of quantitative analysis. Zou et al.  developed a modified algorithm of background removal based on wavelet transform for spectrum correction and applied to low-alloy steel samples. This method can effectively improve the quality of the signals and the accuracy of the regression model. Sun et al.  presented a method that can automatically estimate and correct varying continuum background emission. Simulations and experiments were made to successfully prove the efficiency of the method. The proposed method scarcely needs people’s intervention and can automatically and flexibly estimate varying continuum backgrounds over a very wide spectrum range. Another way to deduct a continuous background is to add a polarizer to the collected light path. Penczak et al.’s  research results show that the continuous background of the Al plasma emission spectra induced by 800 nm femtosecond pulse laser is strongly polarized. The use of a polarizer can effectively filter the continuous spectrum, thus improving the signal-to-noise ratio and the signal-to-back ratio of the characteristic spectrum.
2.3. Spectral normalization
In order to increase the stability of the signal, the analyte signal intensity can be normalized using a parameter representative of the actual plasma conditions. In general, there are three main standardized methods : (1) normalization by using the intensity of an internal standard line; (2) normalization by using a reference signal; and (3) compensation for the plasma conditions. Castro et al.  used 12 different types of data normalization to reduce the interference matrix and to improve the calibration models. Their findings show that the application of normalization modes was useful to compensate for the differences among sample matrices. Models without normalization presented two- to fivefold higher errors. Karki et al.  studied the analytical performance of six different spectrum normalization techniques, namely internal normalization, normalization with total light, normalization with background along with their three-point smoothing methods for quantification of Cr, Mn and Ni in stainless steel. The final results show the superiority of internal normalization technique over normalization with total light and normalization with background techniques irrespective of whether it is Cr, Ni or Mn analysis. Wang et al. presented three spectrum standardized methods in order to improve the reproducibility of LIBS measurements which are named the spectrum standardization approach , the sampled spectrum standardization approach  and the multivariate spectrum standardization method, respectively. In spectral standardization, a particular example is the use of acoustic signals [32, 33] or laser-induced plasma image .
2.4. Automatic identification of emission lines
LIBS can excite all the elements of the spectrum of the sample, so reliable and fast identification of emission lines in laser ablation of multicomponent samples is crucial. Labutin et al.  applied an algorithm to automatically identify emission lines in LIBS. The algorithm is implemented by three parts: simulation of the set of spectra corresponding to different temperature and electron density, searching the best correlated pair of a model spectrum and an experimental one, and attributing the peaks with certain lines. Ukwatta et al.  consider the problem of element detection as a multilabel classification problem, using support vector machines (SVMs) and artificial neural networks (ANNs) for multielement classification. The proposed algorithm is evaluated by using the LIBS image obtained from the experiment. The accuracy of the machine learning method to identify the elements correctly can reach 99%. Mateo et al.  developed the software package SALIPS, which can quickly and semiautomatically identify the spectrum peak and give the element composition of the analytical sample. The software package simulates the spectrum by using the relative intensity of the atomic line in the NIST database. In order to facilitate visual comparison, it can present both the simulated and experimental spectra on the same plot.
3. Calibration methods of LIBS
A number of calibration methods have been applied to various research fields and physical state samples of LIBS quantitative analysis. We cannot involve all published research articles, and only a few of the most commonly used calibration methods are reviewed.
3.1. Univariate analysis
The fitting area intensity corresponding to the transition between lower level and upper level of an atomic species can be expressed as:
where is the experimental parameter, is the atomic number density, is the transition probability, is the upper-level degeneracy and is the partition function at the temperature . For the same sample, if the temperature and density of each laser-induced plasma are constant, then is proportional to the elemental concentration . If there are a series of samples with different , one can establish a calibration line between spectral intensity and element concentration.
where and are model parameters, is the random error and is the number of samples. The parameter is the estimate value of , namely
In the regression analysis, the best estimate of and is obtained based on a set of and , which makes the and to the nearest degree. For example, Bhatt et al.  choose Ce II 413.38, 418.65, and 439.16 nm to establish a univariate linear calibration curve, as shown in Figure 3.
Correlation coefficient , also called the Pearson coefficient, is often used to denote the correlation between and , which is defined as:
Correlation coefficient , and it is closer to 1, indicating the better relevance. Most of the LIBS papers report which provides fast information about the correlation of the data and consequently a fast first knowledge about the prediction ability of the model since poor correlation necessarily implies poor predictive ability. However, it should be noticed that a model with a value of close to 1 may indeed have a poor accuracy for prediction .
Precision is described by the standard deviation (SD), the relative standard deviation (RSD in %) and the root-mean-square error (RMSE), which can be expressed as:
In order to describe the lower limits of a quantitative model, the limit of detection (LOD) can be calculated by the following equation:
where is the standard deviation of the background and is the reciprocal of the slope of the calibration curve. The calculated values of LOD for different elements are presented in Table 1.
3.2. Multivariate analysis
3.2.1. Multiple linear regression
For LIBS, the line of an element is not one. If there are variables and samples, then
where is the intensity of different spectral lines from the same element, is residual, is the sum of squares of the residuals and is the estimate value of . One can get , , …, when value achieves the minimum value.
For example, Chen et al.  used the multiple linear regression method to quantitatively analyze chromium in potatoes. The characteristic line of Cr can be considered that the concentration of Cr () has a relationship with the intensity of Cr and/or the other corresponding elements. They normalize the quantitative analysis of Cr by considering the influence of the Ca matrix. Four independent variables (, , and ) are used to test the performance of different linear regression methods, where is the intensity of Cr I 425.43 nm, is the sum of three Cr lines (Cr I:425.43, Cr I 427.48 and Cr I 428.97 nm), is the intensity of Ca I 431.86 nm and is the sum of five Ca lines (Ca I 422.67, Ca I 428.30, Ca I 430.25, Ca I 430.77, and Ca I 431.86 nm). Different combinations of the four independent variables were selected for unary, binary, ternary and quaternary linear regression analyses. The results of quantitative analysis of Cr element by the linear regression method with different variables are indicated in Table 2.
|Calibrate method||Input variables||Predicted value (μg/g)||Relative error (%)|
3.2.2. Principal component regression
In the LIBS quantitative analysis, the calculated concentration is affected by the lines of objective elements and other elements. In the study of empirical questions, in order to analyze the problem comprehensively and systematically, one must consider many spectral lines of many elements. Because each line reflects the information of the element concentration in varying degrees, and the lines have a certain correlation with each other, the information reflected in the calculation overlaps to some extent. In the study of multivariate problems by statistical methods, too many variables will increase the amount of calculation and increase the complexity of the problem. It is hoped that the variables involved in the process of quantitative analysis are less and the amount of information is more. Principal component analysis (PCA) is adapted to this requirement and is an ideal tool for solving such problems. Principal component analysis can transform a set of variables that may have correlation into a set of linearly uncorrelated variables by orthogonal transformation. The variable after the conversion is called the principal component. Principal component regression (PCR) is a regression analysis method for analyzing multiple regression and is based on PCA. In general, predicting the concentration by PCR can be divided into three steps: first, the PCA is performed on the data (spectral) matrix of the original independent variables, and the appropriate number of principal components is selected by finding the eigenvalue, eigenvector, variance contribution rate and cumulative contribution rate. Second, the selected principal component is analyzed by the ordinary least-squares method. Finally, the strongest possible correlations between the orthogonal PC scores and elemental composition are established. When selecting principal components by PCR, only the independent variables are taken into account, and the dependent variables are ignored. It can reduce the dimension of variables and address the problem of multiple collinearity but cannot distinguish noise when there is a lot of noise in the independent variable (signal) and lose some information of the original variables, so a better regression model will not be obtained.
Death et al.  applied PCR to determine the elemental composition of a series of run-of-mine (ROM) iron ore samples. LIBS spectral data were recorded in three separate spectral regions (250 nm, 400 nm and 750 nm) to measure major, minor and trace components of the iron ore sample pellets. Background stripping, normalization and spectral cleaning were applied to minimize RSD of the LIBS data. PCR analysis was used to produce calibration models of Fe, Al, Si, Mn, K and P. Independent LIBS measurement data are used to verify these calibration models. The model for Fe, Al, Si and K is 0.99, 0.98, 0.99 and 0.84, respectively. As an example, PCR calibration model of Fe is shown in Figure 4 .
3.2.3. Partial least squares regression
The main purpose of PCR is to extract relevant information hidden in the spectral line and then used to predict the concentration. This approach allows one to use only those independent variables, and the noise will be eliminated, so as to improve the quality of the predictive model. However, PCR still has some defects, and some useful variables whose correlation is very small are easily missed when the principal component is selected. If we choose for each component, it is too difficult. Partial least squares regression (PLSR) is a new multivariate statistical data analysis method. It mainly studies the regression modeling of multiple dependent variables to multiple independent variables. PLSR is more effective, especially when the variables are highly linearly correlated. In addition, PLSR solves the problem that the number of samples is less than the number of variables. Partial least squares (PLS) is the advantage of three analytical methods, which are PCA, canonical correlation analysis and multiple linear regression analysis. Both PLS and PCA try to extract the maximum information reflecting the data variation, but PCA only considers one independent variable matrix, while PLS has a response matrix, so it has predictive function. PLS avoids potential problems such as nonnormal distribution of data, factor indeterminacy and unidentifiable models. PLS has two types (PLS-1 and PLS-2), and PLS-1 corresponds to the case where there is only one dependent variable. PLS-2 corresponds to the case where there are several dependent variables. Although PLSR is more complex than PCR, and the tendency of overfitting is stronger, better results can be obtained by using PLS to analyze low-precision data or high-complexity systems.
The input variable of PLS can be characteristic spectral lines , partial spectral region  or full spectrum . Amador-Hernandez et al.  used PLS-1 to quantify gold and silver gold and silver in Au-Ag-Cu alloys. The influence of spectral region (266–340 nm, 266–269/326–340 nm and 269–313 nm), laser energy (3 mJ, 8 mJ), background correction and integration time on the quantitative analysis of PLS was studied, respectively (Table 3).
|Region I||Region II||Region III|
|8 mJ||3 mJ||8 mJ||3 mJ||8 mJ||3 mJ|
|Integrated spectra||WBC||3.47||3.27||3.27||2.74 (2)||2.72 (2)||2.63 (2)|
|BC||1.97||2.27||2.56||2.27||2.58 (2)||2.09 (2)|
|Time-resolved spectra||WBC||3.11||2.27||3||2.02||2.73 (2)||3.96 (2)|
|BC||3.07||2.32||3.13||2.23||2.47 (2)||3.73 (2)|
3.2.4. Artificial neural network
To overcome the poor precision of the calibration curve methods and the limitations of nonlinear problems, scholars have proposed the use of statistical methods for the quantitative analysis by LIBS. Artificial neural networks (ANNs) are computing systems inspired by the biological neural networks that constitute animal brains. Such systems learn (progressively improve performance) to do tasks by considering examples, generally without task-specific programming. The following outstanding advantages of artificial neural networks have attracted great attention in recent years: (1) it can fully approximate any complex nonlinear relationship; (2) all the quantitative or qualitative information is stored in the neurons in the network, so it has strong robustness and fault tolerance; (3) ANNS adopts the parallel distribution processing method, so that it can perform a large number of operations quickly; (4) ANNS can learn and adapt to unknown or uncertain system; and (5) it can handle both quantitative and qualitative information at the same time. An artificial neuron network (ANN) is a computational model based on the structure and functions of biological neural networks. Artificial neural networks usually contain input layer, hidden layer (competitive layer) and output layer, as shown in Figure 5.
For example, El Haddad et al.  used an artificial neuron network to analyze the heavy metals in soil and predict the concentration of element. They used average relative error of calibration REC (%) and the average relative error of prediction REP (%) to evaluate the predictive quality of the ANN models. REC and REP were preferred to RMSE because they provide percentage instead of absolute values. Before using the artificial neural network, it is necessary to optimize the parameters such as the number of neurons and training times.
where Nc is the number of samples in the calibration and Np is the number of samples in the prediction sets, respectively (Table 4).
|Output||Input element||REC (%)||REP (%)|
|Al||Al, Ca, Ba, Fe, Ti||18.7 ± 0.8||19.3 ± 2.1|
|Ca||Ca, Ba, Fe, Ti||9.4 ± 0.4||15.2 ± 0.8|
|Fe||Fe, Ba, Ca, Ti||15.5 ± 0.6||16.8 ± 0.9|
3.2.5. Support vector regression
The neural network structure design depends on the designer’s experience and prior knowledge, while support vector machine (SVM) is based on statistical theory, which has a strict theoretical and mathematical basis. The neural network learning algorithm lacks quantitative analysis and complete theoretical support, and it also needs a lot of samples to learn. SVM is often used to pattern recognition, classification and regression analysis of small samples, nonlinear and high-dimensional data and can achieve very good results. SVM is based on the principle of structural risk minimization, which can ensure that the learning machine has a good generalization ability. SVM for regression prediction is called support vector regression (SVR). SVR also can guarantee the global optimality of the algorithm and avoid the local minimum problem that the neural network cannot solve. Therefore, when there are a small number of samples, it is better to use SVR instead of neural network. It is important to note that the selection of optimized significant penalty parameter C and the kernel parameter of RBF-δ is more sensitive for the SVR model.
For example, Gu et al.  used three segmental spectra of 393–397 nm, 422–423 nm and 425–427 nm as the input variables of SVR model to predict the content of Cr in soil samples. They got better predictive results that and the absolute relative error is 2.61% and the slope of the calibration curve is closer to 1, as shown in Figure 6.
3.3. Calibration-free laser-induced breakdown spectroscopy
LIBS offers a strong potential for analysis in situ and in real time, not requiring complex sample preparation. This allows it to be applied quickly and extensively to qualitative analysis, but quantitative analysis is very difficult. Even with a given experimental configuration, the laser-induced breakdown spectrum is not only dependent on the concentration of the analyte but also dependent on the composition of the matrix and their polymerization state. Matrix effects play an important role in quantitative analysis of LIBS. In order to overcome the matrix effect, Ciucci et al.  proposed the calibration-free laser-induced breakdown spectroscopy (CF-LIBS) approach which takes the matrix into account as a part of the analytical problem. In local thermodynamic equilibrium (LTE), excited levels are populated according to the Boltzmann distribution and ionization states are populated according to the Saha-Boltzmann equilibrium equation. Each spectral line is represented as a point in a Boltzmann plane where the slope and intercept correspond to the plasma temperature and the concentration of the corresponding element, respectively (Figure 7).
Ciucci et al.  proposed CF-LIBS for the first time and used it for quantitative analysis of the composition of metallic alloys and quantitative determination of the composition of the atmosphere. CF-LIBS has been applied to many samples, such as aluminum alloys, steel and iron alloys, precious alloys for jewelry, copper alloys, archeological copper artifacts, glasses, pigments on roman frescoes and on parchments, soils and rocks, meteorites, coral skeletons and human hair. However, the accuracy of CF-LIBS is still not high.
4. The comparison of calibration methods
LIBS is an analytical technique that can inspire all the elements in the sample. Univariate analysis uses only partial spectral information and suffers from the strong effects of plasma instability. More importantly, strong matrix effects prevented to apply simple calibration curves. There is no doubt that the multivariate analysis is superior to univariate analysis. This has been proven by numerous researchers, for instance, the analysis of rocks , rare earth elements , glass , cerium oxide , alloy steel , liquid steel , soil [54, 55], soybean oil , PZT (Lead Zirconate Titanate) ceramics , Pb in navel orange , Marcellus Shale , tailing cores , geologically diverse samples , steel melt , slurry , iron ore  and pellets of plant materials .
Many multivariate analysis methods have been applied to the quantitative analysis of LIBS, especially chemometric. Generally, the most common chemometric technique applied to concentration measurement by LIBS is PLS. It has been applied to many fields of analysis, such as soil [55, 65, 66], steel [67, 68, 69], glass , rock , iron ore  and coal [71, 72]. The rest of the analysis methods are PCR [50, 73, 74], LASSO [75, 76], kNN , ANN [78, 79, 80, 81], SVM [82, 83] and so on. PLS has been implemented either to calculate the concentrations of a single element (PLS-1) or to simultaneously calculate the concentrations of more than one element (PLS-2). In addition to PLS, other linear (MLR [84, 85], PCR [50, 73] and LASSO [75, 76, 77]) and nonlinear regression methods (ANN [86, 87]) have been applied to LIBS quantitative analysis. In order to discern the most effective models for interpreting chemical abundances from LIBS spectra of geological samples, Boucher et al.  had studied nine kinds of linear and nonlinear regression methods. The advantages and disadvantages of various methods are introduced, as shown in Table 5. The final results show that nonlinear methods tend to overfit the data and predict less accurately, while the linear methods proved to be more generalizable with better predictive performance. The performance of different models for different oxides is different. At present, multivariate analysis, especially PLS, is our best choice. But it should be pointed out that multivariate quantitative analyses present a high risk of overfitting.
|PLS||Used when ||Provides a complex model in which all coefficients are linear combinations of the original channels. Involves a complex optimization problem with no simple, closed-form representation.||Linear, uses all channels (not sparse)|
|LASSO||Provides an interpretable model, selects subset of predictors with the strongest effects on the response variable. Can be used for feature selection when less data are available.||Arbitrarily chooses one covariate from a group of highly collinear covariates to use in the model and discards the rest.||Linear, sparse, eliminates noisy channels|
|Elastic net||Performs well in the p > > N case. Provides an interpretable model that is more stable than the lasso. Useful for feature selection.||Cannot be used for feature selection in situations when less data are available because it overwhelms the data with too many model variables.||Linear, sparse, eliminates noisy channels|
|PCR||Decorrelates the data and reduces its dimensionality, combating the “curse of dimensionality”||Higher-order polynomial kernels tend to overfit the training set and poorly predict the testing set in this application.||May be linear or nonlinear; both use all channels|
|SVR||Performs well with a linear kernel. Can be either linear or nonlinear depending on the kernel.||As above, polynomial kernels tend to overfit the training set and poorly predict the testing set in this application.||May be linear or nonlinear; either uses all channels|
|kNN||Requires no model training other than choosing the number of neighbors, reducing run time and making it scale well to large data sets.||Tends to overfit the training data and is only as effective as the distance metric used to compare samples.||Nonlinear, uses all channels|
CF-LIBS can overcome the influence of matrix effect, but the poor analysis accuracy has been the fatal shortcoming of CF-LIBS. This is mainly due to the fact that the laser-induced plasma is a very complex object and its realistic description is not attainable with simple mathematical models . A number of researchers have made some modifications to the CF-LIBS algorithm, such as self-absorption [89, 90]. In recent years, several research groups [91, 92, 93, 94] began using standard samples to improve the accuracy of the nonstandard analysis. Cavalcanti et al.  presented and used one-point-calibration CF-LIBS to analyze a set of copper-based samples. The results show that the new method achieves similar or even higher accuracy than the calibration curve.
This work was supported by the National Natural Science Foundation of China (Grant nos. 11075184 and 61505223) and the Knowledge Innovation Program of the Chinese Academy of Sciences (Grant no. Y03RC21124).
Trautner S, Jasik J, Parigger CG, et al. Laser-induced optical breakdown spectroscopy of polymer materials based on evaluation of molecular emission bands. Spectrochimica Acta A. 2017; 174:331-338. DOI: 10.1016/j.saa.2016.11.045
Choi J-H, Moon Y, Lee S-H, et al. Wavelength dependence of the ablation characteristics of cu (in, Ga) Se2 solar cell films and its effects on laser induced breakdown spectroscopy analysis. International Journal of Precision Engineering and Manufacturing – Green Technology. 2016; 3:167-171. DOI: 10.1007/s40684-016-0021-7
Labutin TA, Lednev VN, Ilyin AA, et al. Femtosecond laser-induced breakdown spectroscopy. Journal of Analytical Atomic Spectrometry. 2016; 31:90-118. DOI: 10.1039/c5ja00301f
Lednev VN, Pershin SM, Sdvizhenskii PA, et al. Laser induced breakdown spectroscopy with picosecond pulse train. Laser Physics Letters. 2017; 14. DOI: 10.1088/1612-202X/aa4fc0
Shakeel H, Haq SU, Aisha G, et al. Quantitative analysis of al-Si alloy using calibration free laser induced breakdown spectroscopy (CF-LIBS). Physics of Plasmas 2017; 24: 063516. DOI: 10.1063/1.4985327
Gehlen CD, Roth P, Aydin U, et al. Time-resolved investigations of laser-induced plasmas generated by nanosecond bursts in the millijoule burst energy regime. Spectrochimica Acta B. 2008; 63:1072-1076. DOI: 10.1016/j.sab.2008.08.007
Smaldone A, De Bonis A, Galasso A, et al. Inverse calibration free fs-LIBS of copper-based alloys. Zeitschrift Fur Physikalische Chemie-International Journal of Research in Physical Chemistry & Chemical Physics. 2016; 230:1201-1217. DOI: 10.1515/zpch-2015-0704
Ahamer CM, Eschlboeck-Fuchs S, Kolmhofer PJ, et al. Laser-induced breakdown spectroscopy of major and minor oxides in steel slags: Influence of detection geometry and signal normalization. Spectrochimica Acta B. 2016; 122:157-164. DOI: 10.1016/j.sab.2016.06.008
Zhang T, Wu S, Dong J, et al. Quantitative and classification analysis of slag samples by laser induced breakdown spectroscopy (LIBS) coupled with support vector machine (SVM) and partial least square (PLS) methods. Journal of Analytical Atomic Spectrometry. 2015; 30:368-374. DOI: 10.1039/c4ja00421c
Yi RX, Li JM, Yang XY, et al. Spectral interference elimination in soil analysis using laser-induced breakdown spectroscopy assisted by laser-induced fluorescence. Analytical Chemistry. 2017; 89:2334-2337. DOI: 10.1021/acs.analchem.6b03969
Yelameli M, Thornton B, Takahashi T, et al. Support Vector Machine Based Classification of Seafloor Rock Types Measured Underwater Using Laser Induced Breakdown Spectroscopy 2016. Shanghai: Oceans; 2016
Yang HX, Fu HB, Wang HD, et al. Laser-induced breakdown spectroscopy applied to the characterization of rock by support vector machine combined with principal component analysis. Chinese Physics B. 2016; 25:290-295. DOI: Artn 06520110.1088/1674-1056/25/6/065201
Hartig KC, Ghebregziabher I, Jovanovic I. Standoff detection of uranium and its isotopes by femtosecond filament laser ablation molecular isotopic spectrometry. Science Report-Uk. 2017; 7. DOI: 10.1038/srep43852
Yin H, Hou Z, Zhang L, et al. Cement raw material quality analysis using laser-induced breakdown spectroscopy. Journal of Analytical Atomic Spectrometry. 2016; 31:2384-2390. 10.1039/C6JA00323K. DOI: 10.1039/C6JA00323K
Zhang T, Yan C, Qi J, et al. Classification and discrimination of coal ash by laser-induced breakdown spectroscopy (LIBS) coupled with advanced chemometrics methods. Journal of Analytical Atomic Spectrometry. 2017. DOI: 10.1039/C7JA00218A
Caceres JO, Lopez JT, Telle HH, et al. Quantitative analysis of trace metal ions in ice using laser-induced breakdown spectroscopy. Spectrochimica Acta B. 2001; 56:831-838. DOI: 10.1016/S0584-8547(01)00173-2
Al-Adel FF, Dastageer MA, Gasmi K, et al. Optimization of a laser induced breakdown spectroscopy method for the analysis of liquid samples. Journal of Applied Spectroscopy. 2013; 80:767-770. DOI: 10.1007/s10812-013-9839-8
Sobral H, Sangines R, Trujillo-Vazquez A. Detection of trace elements in ice and water by laser-induced breakdown spectroscopy. Spectrochimica Acta B. 2012; 78:62-66. DOI: 10.1016/j.sab.2012.09.005
Chen C, Niu G, Shi Q, et al. Laser-induced breakdown spectroscopy technique for quantitative analysis of aqueous solution using matrix conversion based on plant fiber spunlaced nonwovens. Applied Optics. 2015; 54:8318-8325. DOI: 10.1364/ao.54.008318
Zhu DH, LZ W, Wang B, et al. Determination of ca and mg in aqueous solution by laser-induced breakdown spectroscopy using absorbent paper substrates. Applied Optics. 2011; 50:5695-5699. DOI: 10.1364/AO.50.005695
Chen ZJ, Godwal Y, Tsui YY, et al. Sensitive detection of metals in water using laser-induced breakdown spectroscopy on wood sample substrates. Applied Optics. 2010; 49:C87-C94. DOI: 10.1364/AO.49.000C87
Chen ZJ, Li HK, Liu M, et al. Fast and sensitive trace metal analysis in aqueous solutions by laser-induced breakdown spectroscopy using wood slice substrates. Spectrochimica Acta B. 2008; 63:64-68. DOI: 10.1016/j.sab.2007.11.010
Alamelu D, Sarkar A, Aggarwal SK. Laser-induced breakdown spectroscopy for simultaneous determination of Sm, Eu and Gd in aqueous solution. Talanta. 2008; 77:256-261. DOI: 10.1016/j.talanta.2008.06.021
Zou XH, Guo LB, Shen M, et al. Accuracy improvement of quantitative analysis in laser-induced breakdown spectroscopy using modified wavelet transform. Optics Express. 2014; 22:10233-10238. DOI: 10.1364/OE.22.010233
Sun LX, Yu HB. Automatic estimation of varying continuum background emission in laser-induced breakdown spectroscopy. Spectrochimica Acta B. 2009; 64:278-287. DOI: 10.1016/j.sab.2009.02.010
Penczak JS, Liu YM, Gordon RJ. Polarization resolved laser-induced breakdown spectroscopy of Al. The Journal of Physical Chemistry A. 2009; 113:13310-13317. DOI: 10.1021/jp904728n
Tognoni E, Cristoforetti G. Signal and noise in laser induced breakdown spectroscopy: An introductory review. Optics and Laser Technology. 2016; 79:164-172. DOI: 10.1016/j.optlastec.2015.12.010
Castro JP, Pereira-Filho ER. Twelve different types of data normalization for the proposition of classification, univariate and multivariate regression models for the direct analyses of alloys by laser-induced breakdown spectroscopy (LIBS). Journal of Analytical Atomic Spectrometry. 2016; 31:2005-2014. DOI: 10.1039/c6ja00224b
Karki V, Sarkar A, Singh M, et al. Comparison of spectrum normalization techniques for univariate analysis of stainless steel by laser-induced breakdown spectroscopy. Pramana-Journal of Physics. 2016; 86:1313-1327. DOI: 10.1007/s12043-015-1180-8
Wang Z, Li LZ, West L, et al. A spectrum standardization approach for laser-induced breakdown spectroscopy measurements. Spectrochimica Acta B. 2012; 68:58-64. DOI: 10.1016/j.sab.2012.01.005
Li LZ, Wang Z, Yuan TB, et al. A simplified spectrum standardization method for laser-induced breakdown spectroscopy measurements. Journal of Analytical Atomic Spectrometry. 2011; 26:2274-2280. DOI: 10.1039/c1ja10194c
Hrdlicka A, Zaoralkova L, Galiova M, et al. Correlation of acoustic and optical emission signals produced at 1064 and 532 nm laser-induced breakdown spectroscopy (LIBS) of glazed wall tiles. Spectrochimica Acta B. 2009; 64:74-78. DOI: 10.1016/j.sab.2008.10.043
Anabitarte F, Rodriguez-Cobo L, Lopez-Higuera JM, et al. Normalization of laser-induced breakdown spectroscopy spectra using a plastic optical fiber light collector and acoustic sensor device. Applied Optics. 2012; 51:8306-8314. DOI: 10.1364/AO.51.008306
Ni Z, Chen X, Fu H, et al. Study on quantitative analysis of slag based on spectral normalization of laser-induced plasma image. Front Phys-Beijing. 2014; 9:439-445. DOI: 10.1007/s11467-014-0433-1
Labutin TA, Zaytsev SM, Popov AM. Automatic identification of emission lines in laser-induced plasma by correlation of model and experimental spectra. Analytical Chemistry. 2013; 85:1985-1990. DOI: 10.1021/ac303270q
Ukwatta E, Samarabandu J, Hall M. Machine vision system for automated spectroscopy. Machine Vision Application. 2012; 23:111-121. journal article. DOI: 10.1007/s00138-011-0338-8
Mateo MP, Nicolas G, Pinon V, et al. Versatile software for semiautomatic analysis and processing of laser-induced plasma spectra. Spectrochimica Acta B. 2005; 60:1202-1210. DOI: 10.1016/j.sab.2005.05.010
Bhatt CR, Yueh FY, Singh JP. Univariate and multivariate analyses of rare earth elements by laser-induced breakdown spectroscopy. Applied Optics. 2017; 56:2280-2287. DOI: 10.1364/AO.56.002280
Golbraikh A, Tropsha A. Beware of q2! Journal of Molecular Graphics and Modelling. 2002; 20:269-276. DOI: 10.1016/S1093-3263(01)00123-1
Chen T, Huang L, Yao M, et al. Quantitative analysis of chromium in potatoes by laser-induced breakdown spectroscopy coupled with linear multivariate calibration. Applied Optics. 2015; 54:7807-7812. DOI: 10.1364/ao.54.007807
Death DL, Cunningham AP, Pollard LJ. Multi-element analysis of iron ore pellets by laser-induced breakdown spectroscopy and principal components regression. Spectrochimica Acta B. 2008; 63:763-769. DOI: 10.1016/j.sab.2008.04.014
Zhang TT, Shu R, Liu PX, et al. Elemental analysis of rock with remote laser induced breakdown spectroscopy. Spectroscopy and Spectral Analysis. 2017; 37:594-598. DOI: 10.3964/j.issn.1000-0593(2017)02-0594-05
Zaytsev SM, Popov AM, Chernykh EV, et al. Comparison of single- and multivariate calibration for determination of Si, Mn, Cr and Ni in high-alloyed stainless steels by laser-induced breakdown spectrometry. Journal of Analytical Atomic Spectrometry. 2014; 29:1417-1424. DOI: 10.1039/c3ja50389e
Yao SC, JD L, Li JY, et al. Multi-elemental analysis of fertilizer using laser-induced breakdown spectroscopy coupled with partial least squares regression. Journal of Analytical Atomic Spectrometry. 2010; 25:1733-1738. DOI: 10.1039/c0ja00027b
Amador-Hernandez J, Garcia-Ayuso LE, Fernandez-Romero JM, et al. Partial least squares regression for problem solving in precious metal analysis by laser induced breakdown spectrometry. Journal of Analytical Atomic Spectrometry. 2000; 15:587-593. DOI: 10.1039/B000813N
El Haddad J, Villot-Kadri M, Ismael A, et al. Artificial neural network for on-site quantitative analysis of soils using laser induced breakdown spectroscopy. Spectrochimica Acta B. 2013; 79-80:51-57. DOI: 10.1016/j.sab.2012.11.007
Gu Y-H, Zhao N-J, Ma M-J, et al. Monitoring the heavy element of Cr in agricultural soils using a mobile laser-induced breakdown spectroscopy system with support vector machine. Chinese Physics Letters. 2016; 33. DOI: 10.1088/0256-307x/33/8/085201
Ciucci A, Corsi M, Palleschi V, et al. New procedure for quantitative elemental analysis by laser-induced plasma spectroscopy. Applied Spectroscopy. 1999; 53:960-964. DOI: 10.1366/0003702991947612
Dyar MD, Fassett CI, Giguere S, et al. Comparison of univariate and multivariate models for prediction of major and minor elements from laser-induced breakdown spectra with and without masking. Spectrochimica Acta B. 2016; 123:93-104. DOI: 10.1016/j.sab.2016.07.010
Devangad P, Unnikrishnan VK, Tamboli MM, et al. Quantification of Mn in glass matrices using laser induced breakdown spectroscopy (LIBS) combined with chemometric approaches. Analytic Methods-Uk. 2016; 8:7177-7184. DOI: 10.1039/c6ay01930g
Tripathi MM, Eseller KE, Yueh FY, et al. Multivariate calibration of spectra obtained by laser induced breakdown spectroscopy of plutonium oxide surrogate residues. Spectrochimica Acta B. 2009; 64:1212-1218. DOI: 10.1016/j.sab.2009.09.003
Cong ZB, Sun LX, Xin Y, et al. Quantitative analysis of alloy steel based on laser induced breakdown spectroscopy with partial least squares method. Spectroscopy and Spectral Analysis. 2014; 34:542-547. DOI: 10.3964/3.issn.1000-0593(2014)02-0542-06
Yang Y, Wang P, Ma C. Quantitative analysis of liquid steel component by LIBS based on improved multivariate nonlinear model. Laser & Optoelectronics Progress 2016; 53: 053002_053001-053002_053007
Yang NF, Eash NS, Lee J, et al. Multivariate analysis of laser-induced breakdown spectroscopy spectra of soil samples. Soil Science. 2010; 175:447-452. DOI: 10.1097/SS.0b013e3181f516ea
Yu K, Zhao Y, Liu F, et al. Laser-induced breakdown spectroscopy for determining content of Pb and cd in soil. Transactions of the Chinese society of. Agricultural Engineering. 2016; 32:197-203
Wu Y, Sun T, Liu X, et al. Detection of chromium content in soybean oil by laser-induced breakdown spectroscopy. Laser & Optoelectronics Progress 2016; 53:043001–043001–043001–043007
Hernandez-Garcia R, Villanueva-Tagle ME, Calderon-Pinar F, et al. Quantitative analysis of lead Zirconate Titanate (PZT) ceramics by laser-induced breakdown spectroscopy (LIBS) in combination with multivariate calibration. Microchemical Journal. 2017; 130:21-26. DOI: 10.1016/j.microc.2016.07.024
Yang H, Huang L, Liu M, et al. Detection of Pb in navel orange by peel laser induced breakdown spectroscopy coupled with PLS. Chinese Journal of Analysis Laboratory. 2016; 35:760-764
Sanghapi HK, Jain J, Bol’shakov A, et al. Determination of elemental composition of shale rocks by laser induced breakdown spectroscopy. Spectrochimica Acta B. 2016; 122:9-14. DOI: 10.1016/j.sab.2016.05.011
Kuhn K, Meima JA, Rammlmair D, et al. Chemical mapping of mine waste drill cores with laser-induced breakdown spectroscopy (LIBS) and energy dispersive X-ray fluorescence (EDXRF) for mineral resource exploration. Journal of Geochemical Exploration. 2016; 161:72-84. DOI: 10.1016/j.gexplo.2015.11.005
Sun L, Yu H, Cong Z, et al. In situ analysis of steel melt by double-pulse laser-induced breakdown spectroscopy with a Cassegrain telescope. Spectrochimica Acta B. 2015; 112:40-48. DOI: 10.1016/j.sab.2015.08.008
Ayyalasomayajula KK, Dikshit V, Yueh FY, et al. Quantitative analysis of slurry sample by laser-induced breakdown spectroscopy. Analytical and Bioanalytical Chemistry. 2011; 400:3315-3322. DOI: 10.1007/s00216-011-4852-3
Hao ZQ, Li CM, Shen M, et al. Acidity measurement of iron ore powders using laser-induced breakdown spectroscopy with partial least squares regression. Optics Express. 2015; 23:7795-7801. DOI: 10.1364/oe.23.007795
Braga JWB, Trevizan LC, Nunes LC, et al. Comparison of univariate and multivariate calibration for the determination of micronutrients in pellets of plant materials by laser induced breakdown spectrometry. Spectrochimica Acta B. 2010; 65:66-74. DOI: 10.1016/j.sab.2009.11.007
Villas-Boas PR, Romano RA, de Menezes Franco MA, et al. Laser-induced breakdown spectroscopy to determine soil texture: A fast analytical technique. Geoderma. 2016; 263:195-202. DOI: 10.1016/j.geoderma.2015.09.018
Gu Y, Zhao N, Ma M, et al. Quantitative analysis of Cr in soils using LIBS with principal components regression. Journal of Optoelectronics Laser. 2016; 27:748-753
Singh M, Karki V, Sarkar A. Optimization of conditions for determination of Cr and Ni in steel by the method of laser-induced breakdown spectroscopy with the use of partial least squares regression. Journal of Applied Spectroscopy. 2016; 83:497-503. DOI: 10.1007/s10812-016-0318-x
Wu S, Zhang T, Tang H, et al. Quantitative analysis of nonmetal elements in steel using laser-induced breakdown spectroscopy combined with random forest. Analytic Methods-Uk. 2015; 7:2425-2432. DOI: 10.1039/c4ay02601b
Sarkar A, Karld V, Aggarwal SK, et al. Evaluation of the prediction precision capability of partial least squares regression approach for analysis of high alloy steel by laser induced breakdown spectroscopy. Spectrochimica Acta B. 2015; 108:8-14. DOI: 10.1016/j.sab.2015.04.002
Shi Q, Niu G, Lin Q, et al. Quantitative analysis of sedimentary rocks using laser-induced breakdown spectroscopy: Comparison of support vector regression and partial least squares regression chemometric methods. Journal of Analytical Atomic Spectrometry. 2015; 30:2384-2393. DOI: 10.1039/c5ja00255a
Dong M, Oropeza D, Chirinos J, et al. Elemental analysis of coal by tandem laser induced breakdown spectroscopy and laser ablation inductively coupled plasma time of flight mass spectrometry. Spectrochimica Acta B. 2015; 109:44-50. DOI: 10.1016/j.sab.2015.04.008
Yuan TB, Wang Z, Li Z, et al. A partial least squares and wavelet-transform hybrid model to analyze carbon content in coal using laser-induced breakdown spectroscopy. Analytica Chimica Acta. 2014; 807:29-35. DOI: 10.1016/j.aca.2013.11.027
Takahashi T, Thornton B, Sato T, et al. Temperature based segmentation for spectral data of laser-induced plasmas for quantitative compositional analysis of brass alloys submerged in water. Spectrochimica Acta B. 2016; 124:87-93. DOI: 10.1016/j.sab.2016.08.025
Yaroshchyk P, Death DL, Spencer SJ. Comparison of principal components regression, partial least squares regression, multi-block partial least squares regression, and serial partial least squares regression algorithms for the analysis of Fe in iron ore using LIBS. Journal of Analytical Atomic Spectrometry. 2012; 27:92-98. DOI: 10.1039/C1JA10164A
Bricklemyer RS, Brown DJ, Turk PJ, et al. Improved intact soil-core carbon determination applying regression shrinkage and variable selection techniques to complete spectrum laser-induced breakdown spectroscopy (LIBS). Applied Spectroscopy 2013; 67:1185-1199. 2013/09/27. DOI: 10.1366/12-06983
Dyar MD, Carmosino ML, Breves EA, et al. Comparison of partial least squares and lasso regression techniques as applied to laser-induced breakdown spectroscopy of geological samples. Spectrochimica Acta B. 2012; 70:51-67. DOI: 10.1016/j.sab.2012.04.011
Boucher TF, Ozanne MV, Carmosino ML, et al. A study of machine learning regression methods for major elemental analysis of rocks using laser-induced breakdown spectroscopy. Spectrochimica Acta B. 2015; 107:1-10. DOI: 10.1016/j.sab.2015.02.003
Li KH, Guo LB, Li JM, et al. Quantitative analysis of steel samples using laser-induced breakdown spectroscopy with an artificial neural network incorporating a genetic algorithm. Applied Optics. 2017; 56:935-941. DOI: 10.1364/AO.56.000935
Wei J, Dong J, Zhang T, et al. Quantitative analysis of the major components of coal ash using laser induced breakdown spectroscopy coupled with a wavelet neural network (WNN). Anal Methods-Uk. 2016; 8:1674-1680. DOI: 10.1039/c5ay02994e
Yang J, Yi C, Xu J, et al. Laser-induced breakdown spectroscopy quantitative analysis method via adaptive analytical line selection and relevance vector machine regression model. Spectrochimica Acta B. 2015; 107:45-55. DOI: 10.1016/j.sab.2015.02.014
Anderson RB, Morris RV, Clegg SM, et al. The influence of multivariate analysis methods and target grain size on the accuracy of remote quantitative chemical analysis of rocks using laser induced breakdown spectroscopy. Icarus. 2011; 215:608-627. DOI: 10.1016/j.icarus.2011.07.034
Zhang Y, Li Y, Y-h G, et al. LIBS quantitative analysis of Cr and Ni in iron alloys with support vector machine (SVM). Spectroscopy and Spectral Analysis. 2016; 36:2244-2248. DOI: 10.3964/j.issn.1000-0593(2016)07-2244-05
Yang H-X, H-B F, Wang H-D, et al. Laser-induced breakdown spectroscopy applied to the characterization of rock by support vector machine combined with principal component analysis. Chinese Physics B. 2016; 25. DOI: 10.1088/1674-1056/25/6/065201
da Silva RM, Milori DMBP, Ferreira EC, et al. Total carbon measurement in whole tropical soil sample. Spectrochimica Acta B. 2008; 63:1221-1224. DOI: 10.1016/j.sab.2008.09.003
Wu Y, Ouyang J, Sun T, et al. Quantitative analysis of lead in soybean oil by double pulse laser-induced breakdown spectroscopy and multivariate calibration method. Journal of Instrumental Analysis. 2016; 35:1015-1020
Li K, Guo L, Li C, et al. Analytical-performance improvement of laser-induced breakdown spectroscopy for steel using multi-spectral-line calibration with an artificial neural network. Journal of Analytical Atomic Spectrometry. 2015; 30:1623-1628. DOI: 10.1039/c5ja00089k
Rezaei F, Karimi P, Tavassoli SH. Effect of self-absorption correction on LIBS measurements by calibration curve and artificial neural network. Applied Physics B: Lasers and Optics. 2014; 114:591-600. DOI: 10.1007/s00340-013-5566-3
Tognoni E, Cristoforetti G, Legnaioli S, et al. Calibration-free laser-induced breakdown spectroscopy: State of the art. Spectrochimica Acta B. 2010; 65:1-14. DOI: 10.1016/j.sab.2009.11.006
Sun LX, Yu HB. Correction of self-absorption effect in calibration-free laser-induced breakdown spectroscopy by an internal reference method. Talanta. 2009; 79:388-395. DOI: 10.1016/j.talanta.2009.03.066
Bulajic D, Corsi M, Cristoforetti G, et al. A procedure for correcting self-absorption in calibration free-laser induced breakdown spectroscopy. Spectrochimica Acta B. 2002; 57:339-353. DOI: 10.1016/S0584-8547(01)00398-6
Gaudiuso R, Dell’Aglio M, De Pascale O, et al. Laser-induced plasma analysis of copper alloys based on local thermodynamic equilibrium: An alternative approach to plasma temperature determination and archeometric applications. Spectrochimica Acta B. 2012; 74-75:38-45. DOI: 10.1016/j.sab.2012.06.034
Cavalcanti GH, Teixeira DV, Legnaioli S, et al. One-point calibration for calibration-free laser-induced breakdown spectroscopy quantitative analysis. Spectrochimica Acta B. 2013; 87:51-56. DOI: 10.1016/j.sab.2013.05.016
Frydenvang J, Kinch KM, Husted S, et al. An optimized calibration procedure for determining elemental ratios using laser-induced breakdown spectroscopy. Analytical Chemistry. 2013; 85:1492-1500. DOI: 10.1021/ac3026652
Aragon C, Aguilera JA. CSigma graphs: A new approach for plasma characterization in laser-induced breakdown spectroscopy. Journal of Quantitative Spectroscopy and Radiative. 2014; 149:90-102. DOI: 10.1016/j.jqsrt.2014.07.026