Expand this Topic clickable element to expand a topic
Skip to content
Optica Publishing Group

Accuracy improvement of single-sample calibration laser-induced breakdown spectroscopy with self-absorption correction

Open Access Open Access

Abstract

The single sample calibration laser-induced breakdown spectroscopy (SSC-LIBS) is quite suitable for the fields where the standard sample is hard to obtain, including space exploration, geology, archaeology, and jewelry identification. But in practice, the self-absorption effect of plasma destroys the linear relationship of spectral intensity and element concentration based on the Lomakin-Scherbe formula which is the guarantee of the high accuracy of the SSC-LIBS. Thus, the self-absorption effect limits the quantitative accuracy of SSC-LIBS greatly. In this work, an improved SSC-LIBS with self-absorption correction (SSC-LIBS with SAC) is proposed for the promotion of quantitative accuracy of SSC-LIBS. The SSC-LIBS with SAC can correct the intensity ratio of spectral lines in the calculation of SSC-LIBS through relative self-absorption coefficient K without complicated preparatory information. The alloy samples and pressed ore samples were used to verify the effect of the SSC-LIBS with SAC. Compared with SSC-LIBS, for alloy samples, the average RMSEP and average ARE of SSC-LIBS with SAC decreased from 0.83 wt.% and 13.75% to 0.40 wt.% and 4.06%, respectively. For the pressed ore samples, the average RMSEP and average ARE of SSC-LIBS with SAC decreased from 4.77 wt.% and 90.48% to 2.34 wt.% and 14.60%. The experimental result indicates that SSC-LIBS with SAC has a great improvement of quantitative accuracy and better universality compared with traditional SSC-LIBS, which is a mighty promotion of the wide application of SSC-LIBS.

© 2022 Optica Publishing Group under the terms of the Optica Open Access Publishing Agreement

1. Introduction

Laser-induced breakdown spectroscopy (LIBS), as a new optical emission spectrometry, has many advantages, such as fast, real-time, in-situ, multi-element simultaneous detection, and no or simple sample preparation [1]. Thus, LIBS has been widely applied to many fields such as environmental monitoring [2], food safety [3], biomedicine [4], metallurgy [5], mining [6], and space exploration [7]. The quantitative analysis of LIBS depends on the calibration methods, including classical calibration [8], calibration-free (CF) [9], multi-energy calibration (MEC) [10], one-point and multi-line calibration (OP-MLC) proposed by our group in 2018 [11], and single-sample calibration (SSC) proposed by our group in 2019 [12]. Most of them are unsuitable for the major element quantification due to the self-absorption effect. CF-LIBS can overcome the matrix effect and quantify the major elements without standard samples. However, its accuracy is not high because the theory is quite ideal and far from reality. SSC-LIBS can achieve the total element determination using only one standard sample. It is quite suitable for the fields where standard samples are hard to obtain or prepare, such as geology [13], archaeology [14], jewelry identification [15], etc.

According to the Lomakin-Scherbe formula, SSC-LIBS must select one optically thin spectral line for each element to perform the quantitative procedure. However, the self-absorption effect always exists in the laser-induced plasma [16], which makes the spectral intensity obtained lower than the real intensity, resulting in a decrease in the accuracy of SSC-LIBS. Over the past few decades, the self-absorption effect has been investigated by many scholars. There are mainly three research directions of exploitation and compensation of self-absorption in the LIBS quantitative model. The first kind is to optimize the experimental parameters and experimental environment, such as laser parameters [17], ambient gas conditions [18], acquisition parameters [19], etc. which can make spectral lines as optical thin as possible to eliminate the effect of self-absorption on quantification. But most application situations limit the use of this method. The second kind is to add auxiliary devices, including a duplicating mirror [20], laser-stimulated absorption [21], microwave-assisted excitation [22], etc. The duplicating mirror method can be used in moderate self-absorption conditions and its sensitivity decreases if the line intensity is too low. The laser-stimulated absorption method needs an optical parametric oscillator (OPO) laser, whose cost is high. The microwave is harmful to the human body. The third kind is to correct self-absorption through the algorithm, including internal reference for self-absorption correction (IRSAC) [23], self-absorption correction algorithm based on blackbody radiation (BRR-SAC) [24], and the self-absorption coefficient (SA) method [2527], which are substantially estimation of self-absorption parameter by different reference objects, There are also several untraditional self-absorption correction methods including the curve of growth (COG) [28,29], the columnar density formulation [30,31], the C-sigma method [32,33], which obtain the element concentration information from the optically thick lines directly [34]. The above methods about self-absorption correction are mainly suitable for CF-LIBS, however, the research about the self-absorption correction method for SSC-LIBS has not been reported yet.

In this work, an improved SSC-LIBS with self-absorption correction (SSC-LIBS with SAC) is proposed. This method is based on the improvement and simplification of the self-absorption coefficient (SA) calculation formula. The intensity ratio of spectral lines in the SSC-LIBS formula was corrected directly by this method without the calculation of a specific SA coefficient. We used the average relative error (ARE) and the root mean square error of prediction (RMSEP) to evaluate the effect of the SAC of the proposed method. The results proved that the proposed method can effectively reduce the influence of the self-absorption effect on SSC-LIBS and significantly improve the quantitative accuracy, which is conducive to the promotion of SSC-LIBS.

2. Self-absorption correction method

There are two common methods to calculate the SA coefficient. The first one is based on the ratio between measured and theoretical intensities or broadenings of the optically thick line [26]. The other is based on multiple standard samples to establish the calibration curve by exponential fitting [27] and the SA coefficient can be calculated through the exponential part of the curve equation. Only the first method can be chosen for SSC-LIBS which can only use single standard sample. Under the local thermodynamic equilibrium (LTE) condition and the plasma is homogeneous, the SA coefficient of the goal spectral line can be expressed as [26]:

$$SA = {(\frac{{\Delta \lambda }}{{\Delta {\lambda _0}}})^{1/\alpha }}$$
where $\Delta \lambda$ is the measured line broadening, $\Delta {\lambda _0}$ is the theoretical line broadening without self-absorption. $\alpha$ is the numerical simulation constant, usually −0.54. In LIBS spectra, the collision broadening is caused by charged particles, also known as Stark broadening (Lorentz profile), which is usually dominant in the plasma [35]. If the effect of ions on Stark broadening is ignored and only the collision effect of electrons is considered, the theoretical full width at half maximum (FWHM) of the Lorentz line profile without self-absorption can be expressed as:
$$\Delta {\lambda _0} = \frac{{2{\omega _s}{n_e}}}{{{n_e}^{ref}}}$$
where ${\omega _s}$ is the Stark half-width parameter, ${n_e}$ is the electron number density of plasma, ${n_e}^{ref}$ is the reference electron number density, usually 1016 or 1017 cm-3, which depends on the condition in the calculation of ${\omega _s}$ [36]. Thus, ${n_e}$ can be expressed through a reference line, which is optically thin, by:
$${n_e} = \frac{{\Delta {\lambda _{SAF}}{n_e}^{ref}}}{{2{\omega _{SAF}}}}$$
where $\Delta {\lambda _{SAF}}$ is FWHM of the reference line, ${\omega _{SAF}}$ is the Stark half-width parameter of reference spectral line. Equation (1) can be rewritten by simultaneous Eq. (1), (2), and (3) as:

$$SA = {\left( {\frac{{\Delta \lambda }}{{\Delta {\lambda_{SAF}}}}} \right)^{{1 / \alpha }}}{\left( {\frac{{{\omega_{SAF}}}}{{{\omega_s}}}} \right)^{{1 / \alpha }}}$$

By defining the following quantities:

$$K = {\left( {\frac{{\Delta \lambda }}{{\Delta {\lambda _{SAF}}}}} \right)^{{1 / \alpha }}}$$
$$M={\left( {\frac{{{\omega_{SAF}}}}{{{\omega_s}}}} \right)^{{1 / \alpha }}}$$

Equation (4) can be rewritten as:

$$SA = K \cdot M$$

In the absence of self-absorption effect, according to the formula of SSC-LIBS [12], the concentration of the element j in the test sample is:

$${C_{jt}} = {{\frac{{{C_{js}} \cdot {I_{jt0}}}}{{{I_{js0}}}}} {\bigg /} {\sum\limits_{i = 1}^n {\frac{{{C_{is}} \cdot {I_{it0}}}}{{{I_{is0}}}}} }}$$
where ${C_{is}}$ and ${C_{js}}$ is the certified concentration of element i and j or oxide i and j in the standard sample, ${{{I_{it0}}} / {{I_{is0}}}}$ is the line intensity ratio of element i or oxide i in the test sample and the standard sample, ${{{I_{jt0}}} / {{I_{js0}}}}$ is the line intensity ratio of element j or oxide j in the test sample and the standard sample, n represents the amount of element type or oxide type of the test sample. The self-absorption coefficient can also be defined as the ratio of the real line intensity to the theoretical line intensity without self-absorption. The line intensity without self-absorption can be expressed as:

$${I_0} = \frac{I}{{SA}}$$

According to Eq. (7) and (9), in the absence of self-absorption, the line intensity ratio of element j or oxide j in the test sample and the standard sample is:

$$\frac{{{I_{jt}}_0}}{{{I_{js}}_0}} = \frac{{{{{I_{jt}}} / {({K_{jt}} \cdot {M_{jt}})}}}}{{{{{I_{js}}} / {({K_{j\textrm{s}}} \cdot {M_{j\textrm{s}}})}}}}$$

The M value of the same spectral line is equal, and the Eq. (9) can be rewritten as:

$$\frac{{{I_{jt}}_0}}{{{I_{js}}_0}} = \frac{{{{{I_{jt}}} / {{K_{jt}}}}}}{{{{{I_{js}}} / {{K_{j\textrm{s}}}}}}}$$

Identically, the line intensity ratio of element i or oxide i in the test sample and the standard sample without self-absorption is:

$$\frac{{{I_{it}}_0}}{{{I_{is}}_0}} = \frac{{{{{I_{it}}} / {{K_{it}}}}}}{{{{{I_{is}}} / {{K_{i\textrm{s}}}}}}}$$

According to Eq. (8), (11), and (12), the content of element j or oxide j in test sample after self-absorption correction is:

$${C_{jt}} = {{\frac{{{C_{js}} \cdot {{{I_{jt}}} / {{K_{jt}}}}}}{{{{{I_{js}}} / {{K_{j\textrm{s}}}}}}}} {\bigg /} {\sum\limits_{i = 1}^n {\frac{{{C_{is}} \cdot {{{I_{it}}} / {{K_{it}}}}}}{{{{{I_{is}}} / {{K_{is}}}}}}} }}$$

The K can be obtained directly by the FWHM ratio of the goal spectral line and the reference line of self-absorption free. However, the reference line without self-absorption is hard to obtain in the actual experiment. Since the hydrogen Balmer alpha (Hα) line is generated by the interaction between plasma and water vapor, it is usually regarded as the spectral line without self-absorption [37,38]. In this work, we selected the Hα line as the reference line. According to Sherbini’s method of calculating electron density based on the Hα line in 2006, formula (3) can be expressed as:

$${n_e} = 8.02 \times {10^{12}}{(\frac{{\Delta {\lambda _{{H_\alpha }}}}}{{{\alpha _{1/2}}}})^{3/2}}$$
where $\Delta {\lambda _{{H_\alpha }}}$ is FWHM of Hα line, ${\alpha _{1/2}}$ is the half width of the reduced Stark profiles of Hα line and it is a weak function of plasma electron density and temperature. The value of ${\alpha _{1/2}}$ can be obtained from Ref. [39]. The final expression of the SA coefficient can be obtained by simultaneous Eq. (1), (2), and (14):

$$SA = \frac{{\Delta {\lambda ^{{1 / \alpha }}}}}{{\Delta {\lambda _{{H_\alpha }}}^{{3 / {2\alpha }}}}}{\left( {\frac{{{n_e}^{ref}{{({{\alpha_{1/2}}} )}^{{3 / 2}}}}}{{1.604 \times {{10}^{13}}{\omega_s}}}} \right)^{{1 / \alpha }}}$$

Defining the following quantities anew:

$$K = \frac{{\Delta {\lambda ^{{1 / \alpha }}}}}{{\Delta {\lambda _{{H_\alpha }}}^{{3 / {2\alpha }}}}}$$
$$M = {\left( {\frac{{{n_e}^{ref}{{\left( {{\alpha _{1/2}}} \right)}^{{3 / 2}}}}}{{1.604 \times {{10}^{13}}{\omega _s}}}} \right)^{{1 / \alpha }}}$$
$K$ can be calculated by FWHM of the Hα line to correct the ratio of the line intensities as Eq. (8), which requires no related parameters of characteristic emission line such as electron collision coefficient. As shown in Fig. 1, it is the flow chart of the correction algorithm. Firstly, the spectra of the test sample and standard sample are collected under the same experimental conditions. Secondly, one spectral line is selected for each element or oxide in the test sample and standard sample. Thirdly, the shapes of the selected lines and the Hα line are extracted for curve fitting through Lorentz function to obtain FWHM of the lines. Then, the K value of the selected spectral line is calculated according to Eq. (5). Finally, the remain quantitative steps of SSC-LIBS are performed based on the corrected spectral line intensity ratio by Eq. (13).

3. Experimental

3.1 Samples and sample preparation

To study the self-absorption correction effect of SSC-LIBS with SAC, alloy samples and ore pressed samples were used in the experiment, in which the alloy samples were aluminum bronze standard samples (BYG1916-1, Luoyang Copper Co., Ltd., China Aluminum Testing Center). To obtain ore samples with different oxide concentration distribution, we mixed dolomite (GBW03116, National Research Center of Testing Techniques for Building Materials) and potassium feldspar (GBW070160, Jinan Zhongbiao Technology Co., Ltd.) powder in the proportion of 5:0, 4:1, 3:2, 2:3, 1:4 respectively, and pressed them under high pressure to make five ore samples. The detailed fabrication process is shown in the following steps. Table 1 lists the element content of the alloy samples and Table 2 lists the oxide content of the pressed ore samples.

  • (1) Mix dolomite and potash feldspar powders in five different proportions and mix them evenly in a high-speed oscillator for 10 min.
  • (2) Transfer the mixed ore powder into a mortar for uniform grinding.
  • (3) After mixing and grinding, take 1 g mixed ore powder on the substrate made with 9 g boric acid powder (GB/T 628-1993, Sinopharm Chemical Reagent Co., Ltd.). Press the mixed ore powder with boric acid substrate into pellets with a diameter of 40 mm under the pressure of 30Mpa.

 figure: Fig. 1.

Fig. 1. The flow diagram of SSC-LIBS with SAC

Download Full Size | PDF

Tables Icon

Table 1. The certified concentration of the elements in aluminum bronze standard samples(wt.%)

Tables Icon

Table 2. The certified concentration of the oxides in pressed ore samples(wt.%)

3.2 Experimental setup

The schematic diagram of the experimental system used in this work is shown in Fig. 2. The Q-switched Nd-YAG laser (wavelength: 532 nm, maximum power: 400 mJ, repetition rate: 1 Hz, pulse width: 8 ns, Beamtech Optronics, Nimma-400) was used to generate the laser. The horizontally emitted laser was reflected by a mirror and then focused on the sample surface through a focusing lens (f = 150 mm) to generate plasma. The spectrum emitted by the plasma was collected by the acquisition head and coupled to a six-channel fiber-optic spectrograph (Avantes B. V, AvaSpec-ULS4096CL-EVO, spectral ranges: 196–874 nm, minimum gate width: 9 μs). The photoelectric sensor in the spectrometer converted the spectral signal into the digital signal and transmitted it to the computer for subsequent processing and analysis. The time sequence between laser and spectrometer was controlled by a digital delay generator (Wuhan N & D Laser Engineering, LDG 3.0). The gate delay and width of the spectrograph were set to 2 μs and 9 μs respectively. To obtain a high spectral signal-to-noise ratio, the laser energy was optimized to 30 mJ and 70 mJ for aluminum bronze samples and ore pressed samples, respectively. In the experiment, the sample was placed on a 3D electric displacement platform (Beijing Jiangyun Juli Technology, DZY110TA-3Z) for generating a new excitation position for each pulse. To reduce the influence of laser energy on the stability of the spectral signal, 5 pulses were accumulated per spectrum and repeated 10 times.

 figure: Fig. 2.

Fig. 2. The diagram of the LIBS setup

Download Full Size | PDF

3.3 Evaluation indexes

To study the improvement effect of SSC-LIBS with SAC on the quantitative accuracy of SSC-LIBS, one sample was selected as the standard sample each time, and the remaining sample was selected as the test sample. The standard sample was used for quantitative analysis of multiple test samples based on SSC-LIBS and SSC-LIBS with SAC, and the average relative error (ARE) and root mean square error of prediction (RMSEP) of the multiple quantitative analysis were calculated. The calculation formulas of ARE and RMSEP are as follows:

$$RMSEP = \sqrt {\frac{{\sum\nolimits_{i = 1}^\textrm{n} {{{({{\widehat x}_i} - {x_i})}^2}} }}{n}}$$
$$ARE(\%) = \frac{{100}}{n}\sum\limits_{i = 1}^n {\frac{{|{{{\widehat x}_i} - {x_i}} |}}{{{x_i}}}}$$
where n is the number of test samples, ${\widehat x_i}$ is the predicted concentration of element or oxide, ${x_i}$ is the certified concentration of element or oxide.

4. Results and discussion

4.1 Quantitative analysis of the alloy sample

According to formula (12), the single-sample calibration method needs to select one spectral line for each element in the sample. To avoid spectral interference and obtain high SNR (signal to noise ratio, SNR), four persistent lines (Cu II 204.38 nm, Fe I 382.04 nm, Al I 396.15 nm, and Mn I 404.14 nm) marked with red arrows in Fig. 3 were selected for the calculation of the SSC-LIBS. The four lines and Hα line (656.28 nm) were used for the fitting curve through Lorentz function. Taking aluminum bronze sample No.3 as an example, Fig. 4 shows the fitting curves of five spectral lines in the spectra of sample No.3. Table 3 lists the K values of selected spectral lines in each sample used in this experiment. According to the K value ratio of the same spectral line between the test sample and the standard sample, the spectral line intensity ratio was corrected in the quantitative process of the SSC-LIBS.

 figure: Fig. 3.

Fig. 3. Spectra of aluminum bronze samples in the spectral range of 200-410 nm.

Download Full Size | PDF

 figure: Fig. 4.

Fig. 4. The fitting curve of five spectral lines in the spectra of sample No.3: (a) – (e) for Cu II 204.38 nm, Fe I 382.04 nm, Al I 396.15 nm, Mn I 404.14 nm and H I 656.28 nm.

Download Full Size | PDF

Tables Icon

Table 3. The K value of element spectral lines of various sample

To study the effect of SSC-LIBS with SAC, we first selected the sample No.2 as the standard sample, the spectral line intensity in the calculation of SSC-LIBS was corrected through the K value obtained above. The corrected intensity was used to predict the element concentration of the remaining four samples and the prediction results were compared with the one without self-absorption correction and the certified concentration of elements. As shown in Fig. 5, the dotted line represents the certified concentrations, the blue dot and red dot represent the predicted concentrations without and with self-absorption correction, respectively (Fig. 5 (a) for Mn element, Fig. 5 (b) for Fe element). Obviously, after self-absorption correction, the element concentration predicted by the SSC-LIBS is closer to the certified concentration, which indicates that the method proposed in this paper has an obvious quantitative improvement effect for the SSC-LIBS.

 figure: Fig. 5.

Fig. 5. Comparison of the concentration calculated by SSC-LIBS with SAC (red) and by SSC-LIBS (blue) vs. the certified concentration in brass samples: (a) Mn and (b) Fe.

Download Full Size | PDF

To further study the quantitative improvement effect of this method when using different samples as standard samples, we selected samples 1–5 as the standard sample respectively and predicted the element concentration of the remaining samples before and after the self-absorption correction. Finally, we obtained the average relative error (ARE) of the four elements before and after the self-absorption correction when using different samples as standard samples.

As is shown in Fig. 6, before self-absorption correction, the predicted concentration ARE of Mn is the largest relative to other elements. This may be due to the certified concentration of the four elements in the five aluminum bronze samples are basically more than 1%, which means that the spectral lines of these elements have strong self-absorption effect, and the certified concentration of Mn element is relatively low in the four elements. For the absolute error of similar level, relative error of elements with lower certified concentration may be larger.

 figure: Fig. 6.

Fig. 6. Comparison of the ARE calculated by SSC-LIBS before SAC and after SAC for aluminum bronze samples: (a) – (e) for No.1-5 sample as the calibration sample.

Download Full Size | PDF

Besides, no matter which sample was used as the standard sample, ARE of elements in aluminum bronze overall decreased significantly compared with that before self-absorption correction, while the decrease of Cu was small. When sample 5 was used as the standard sample, the ARE even increased slightly. This may be because Cu is the matrix element in aluminum bronze, whose content is more than 80%. Furthermore, the content of Cu in five aluminum bronze samples is relatively close, and the self-absorption effect of the Cu spectrum is also relatively close, so the calculated K value of Cu is also relatively close. Thus, the fluctuation of the spectrum itself may exceed the influence of the self-absorption effect on the spectrum of Cu. After correction, the ARE of Cu cannot be reduced always. Even though, the method proposed in this paper still has a satisfactory effect on improving the accuracy of SSC-LIBS as a whole.

In addition, as shown in Fig. 7, the RMSEP of four elements when using different standard samples were calculated, whose trend is similar with the ARE. Table 4 lists the average RMSEP and average ARE when using different samples as standard samples before and after self-absorption correction. The results showed that the average RMSEP of Al, Fe, Mn and Cu decreased from 1.47, 0.52, 0.36, and 0.96 wt.% to 0.63, 0.18, 0.10, and 0.70 wt.%, respectively, and the average ARE decreased from 13.83%, 16.4%, 23.72%, and 1.03% to 5.98%, 4.96%, 4.53%, and 0.76%, respectively. Besides, four elements in five aluminum bronze samples were quantitatively analyzed by classical calibration method, the result showed that ARE of four elements were 3.40%, 2.93%, 9.34%, and 1.22%, which indicate that the method proposed in this work has a good effect on the self-absorption correction of alloy sample and achieve the quantitative accuracy close to the classical calibration method.

 figure: Fig. 7.

Fig. 7. Comparison of the RMSEP calculated by SSC-LIBS before SAC and after SAC for aluminum bronze samples: (a) – (e) for No.1-5 sample as the calibration sample.

Download Full Size | PDF

Tables Icon

Table 4. RMSEP (wt.%) and ARE (%) of elements in aluminum bronze samples through SSC-LIBS and SSC-LIBS with SAC

4.2 Quantitative analysis of the pressed ore sample

To further determine the universality of SSC-LIBS with SAC for different types of samples, we used five different oxide concentrations of ore pressed samples made with dolomite and potash feldspar. As shown in Fig. 8, To avoid spectral interference and obtain high SNR (signal to noise ratio, SNR), five persistent lines (Al I 237.31 nm, Si I 288.16 nm, Ca I 442.54 nm, Na I 589.00 nm, and K I 769.90 nm) marked by the red arrow were selected in the SSC-LIBS.

 figure: Fig. 8.

Fig. 8. Part of Spectrum of pressed ore samples in the spectral range of 235-780 nm.

Download Full Size | PDF

Similarly, we selected ore pressed samples from No. 1 to No. 5 as standard samples respectively, and the remaining four as test samples. Before and after self-absorption correction, the concentration of five oxides in the sample was quantified by SSC-LIBS, and the average RMSEP and average ARE using the different samples as standard samples were calculated. The results of Fig. 9 showed that the average RMSEP of SiO2, Al2O3, CaO, K2O, and Na2O decreased from 6.99, 2.05, 11.69, 2.11, and 0.99 wt.% to 5.37, 0.84, 4.33, 0.88, and 0.26 wt.%, respectively. The average ARE decreased from 25.62%, 34.40%, 31.14%, 212.80%, and 148.42% to 16.04%, 11.24%, 11.48%, 16.22%, and 18.04%, respectively, which indicate that SSC-LIBS with SAC proposed in this paper also has good correction effect for ore samples, and the quantitative accuracy of SSC-LIBS for ore samples is significantly improved.

 figure: Fig. 9.

Fig. 9. Comparison of the ARE and RMSEP calculated by SSC-LIBS before (black) SAC and after (red) SAC for pressed ore samples: (a) for ARE, (b) for (RMSEP)

Download Full Size | PDF

The SSC-LIBS with SAC proposed in this paper is quite suitable for the samples with relatively large element concentration gradient distribution. Besides, the method has a good quantitative accuracy improvement effect for alloy samples and ore samples. It should be noted that SSC-LIBS with the SAC method requires that the Hα line can be observed clearly in the experiment or another reference line of self-absorption free can be obtained, only in this way can the K value be calculated. Therefore, the experiment parameters need to be optimized to obtain a clear Hα line and favorable signal-to-noise ratio. Besides the K value can be used to quantify the degree of self-absorption effect of the same spectral line under different experiment parameters, which is helpful to reduce self-absorption by optimizing experimental parameters.

5. Conclusions

In this work, The SSC-LIBS with SAC has been proposed for overcoming the self-absorption problem and improving the quantitative accuracy of SSC-LIBS. The SSC-LIBS with SAC is based on the improvement and simplification of the calculation formula of the SA coefficient and directly uses the relative SA coefficient K to correct the spectral line intensity ratio in the formula of SSC-LIBS. Alloy and ore pressed samples were used to study the self-absorption correction effect of the method. For alloy samples, after self-absorption correction, the ARE of Al, Fe, Mn, and Cu decreased from 13.83%, 16.4%, 23.72%, and 1.03% to 5.98%, 4.96%, 4.53%, and 0.76%, respectively. For the pressed ore samples, the ARE of SiO2, Al2O3, CaO, K2O, and Na2O decreased from 25.62%, 34.40%, 31.14%, 212.80%, and 148.42% to 16.04%, 11.24%, 11.48%, 16.22%, and 18.04%, respectively. These results showed that the SSC-LIBS with SAC had a good self-absorption correction effect in SSC-LIBS for various types of samples and the quantitative accuracy was significantly improved. The method also has the advantages of simple steps, good universality, and does not require preparatory information such as the Stark half-width parameter. It has great significance for the popularization and application of SSC-LIBS.

Funding

National Natural Science Foundation of China (61705064, 62075069); Natural Science Foundation of Hubei Province (2018CFB773).

Disclosures

The authors declare that there are no conflicts of interest related to this article.

Data availability

Data underlying the results presented in this paper are not publicly available at this time but may be obtained from the authors upon reasonable request.

References

1. D. W. Hahn and N. Omenetto, “Laser-induced breakdown spectroscopy (libs), part ii: Review of instrumental and methodological approaches to material analysis and applications to different fields,” Appl. Spectrosc. 66(4), 347–419 (2012). [CrossRef]  

2. R. S. Harmon, R. E. Russo, and R. R. Hark, “Applications of laser-induced breakdown spectroscopy for geochemical and environmental analysis: A comprehensive review,” Spectrochim. Acta, Part B 87, 11–26 (2013). [CrossRef]  

3. R. A. Multari, D. A. Cremers, J. A. M. Dupre, and J. E. Gustafson, “Detection of biological contaminants on foods and food surfaces using laser-induced breakdown spectroscopy (libs),” J. Agric. Food Chem. 61(36), 8687–8694 (2013). [CrossRef]  

4. V. K. Singh and A. K. Rai, “Prospects for laser-induced breakdown spectroscopy for biomedical applications: A review,” Lasers Med. Sci. 26(5), 673–687 (2011). [CrossRef]  

5. Z. Qingdong, D. Fan, Z. Zhiheng, T. Yun, W. Boyun, X. Yongjun, L. Xiong, Y. Huaqing, G. Lianbo, and L. Xiangyou, “Portable fiber-optic laser-induced breakdown spectroscopy system for the quantitative analysis of minor elements in steel,” Plasma Sci. Technol. 21, 034006 (2018). [CrossRef]  

6. A. Harhira, P. Bouchard, K. Rifai, J. E. Haddad, M. Sabsabi, A. Blouin, and M. Laflamme, “Advanced laser-induced breakdown spectroscopy (libs) sensor for gold mining,” in Proceedings of the 56th Conference of Metallurgists, Vancouver, BC, Canada(2017), pp. 27-30.

7. B. Sallé, J.-L. Lacour, E. Vors, P. Fichet, S. Maurice, D. A. Cremers, and R. C. Wiens, “Laser-induced breakdown spectroscopy for mars surface analysis: Capabilities at stand-off distances and detection of chlorine and sulfur elements,” Spectrochim. Acta, Part B 59(9), 1413–1422 (2004). [CrossRef]  

8. S. Pandhija, N. Rai, A. K. Rai, and S. N. Thakur, “Contaminant concentration in environmental samples using libs and cf-libs,” Appl. Phys. B 98(1), 231–241 (2010). [CrossRef]  

9. A. Ciucci, M. Corsi, V. Palleschi, S. Rastelli, A. Salvetti, and E. Tognoni, “New procedure for quantitative elemental analysis by laser-induced plasma spectroscopy,” Appl. Spectrosc. 53(8), 960–964 (1999). [CrossRef]  

10. D. V. Babos, A. Virgilio, V. C. Costa, G. L. Donati, and E. R. Pereira-Filho, “Multi-energy calibration (mec) applied to laser-induced breakdown spectroscopy (libs),” J. Anal. At. Spectrom. 33(10), 1753–1762 (2018). [CrossRef]  

11. Z. Q. Hao, L. Liu, R. Zhou, Y. W. Ma, X. Y. Li, L. B. Guo, Y. F. Lu, and X. Y. Zeng, “One-point and multi-line calibration method in laser-induced breakdown spectroscopy,” Opt. Express 26(18), 22926–22933 (2018). [CrossRef]  

12. R. Yuan, Y. Tang, Z. Zhu, Z. Hao, J. Li, H. Yu, Y. Yu, L. Guo, X. Zeng, and Y. Lu, “Accuracy improvement of quantitative analysis for major elements in laser-induced breakdown spectroscopy using single-sample calibration,” Anal Chim Acta 1064, 11–16 (2019). [CrossRef]  

13. C. Fabre, “Advances in laser-induced breakdown spectroscopy analysis for geology: A critical review,” Spectrochim. Acta, Part B 166, 105799 (2020). [CrossRef]  

14. D. Anglos, “Laser-induced breakdown spectroscopy in art and archaeology,” Appl. Spectrosc. 55(6), 186A–205A (2001). [CrossRef]  

15. L. Garcıa-Ayuso, J. Amador-Hernández, J. Fernández-Romero, and M. L. De Castro, “Characterization of jewellery products by laser-induced breakdown spectroscopy,” Anal. Chim. Acta 457(2), 247–256 (2002). [CrossRef]  

16. Z. Hao, L. Liu, M. Shen, X. Yang, K. Li, L. Guo, X. Li, Y. Lu, and X. Zeng, “Investigation on self-absorption at reduced air pressure in quantitative analysis using laser-induced breakdown spectroscopy,” Opt. Express 24(23), 26521–26528 (2016). [CrossRef]  

17. H. Cui, Y. Tang, S. Ma, Y. Ma, D. Zhang, Z. Hu, Z. Wang, and L. Guo, “Influence of laser wavelength on self-absorption effect in laser-induced breakdown spectroscopy,” Optik 204, 164144 (2020). [CrossRef]  

18. R. Hedwig, I. Tanra, I. Karnadi, M. Pardede, A. M. Marpaung, Z. S. Lie, K. H. Kurniawan, M. M. Suliyanti, T. J. Lie, and K. Kagawa, “Suppression of self-absorption effect in laser-induced breakdown spectroscopy by employing a penning-like energy transfer process in helium ambient gas,” Opt. Express 28(7), 9259–9268 (2020). [CrossRef]  

19. R. Yi, L. Guo, C. Li, X. Yang, J. Li, X. Li, X. Zeng, and Y. Lu, “Investigation of the self-absorption effect using spatially resolved laser-induced breakdown spectroscopy,” J. Anal. At. Spectrom. 31(4), 961–967 (2016). [CrossRef]  

20. H.-Y. Moon, K. K. Herrera, N. Omenetto, B. W. Smith, and J. D. Winefordner, “On the usefulness of a duplicating mirror to evaluate self-absorption effects in laser induced breakdown spectroscopy,” Spectrochim. Acta, Part B 64(7), 702–713 (2009). [CrossRef]  

21. J. M. Li, L. B. Guo, C. M. Li, N. Zhao, X. Y. Yang, Z. Q. Hao, X. Y. Li, X. Y. Zeng, and Y. F. Lu, “Self-absorption reduction in laser-induced breakdown spectroscopy using laser-stimulated absorption,” Opt. Lett. 40(22), 5224–5226 (2015). [CrossRef]  

22. Y. Tang, J. Li, Z. Hao, S. Tang, Z. Zhu, L. Guo, X. Li, X. Zeng, J. Duan, and Y. Lu, “Multielemental self-absorption reduction in laser-induced breakdown spectroscopy by using microwave-assisted excitation,” Opt. Express 26(9), 12121–12130 (2018). [CrossRef]  

23. L. Sun and H. Yu, “Correction of self-absorption effect in calibration-free laser-induced breakdown spectroscopy by an internal reference method,” Talanta 79(2), 388–395 (2009). [CrossRef]  

24. T. Li, Z. Hou, Y. Fu, J. Yu, W. Gu, and Z. Wang, “Correction of self-absorption effect in calibration-free laser-induced breakdown spectroscopy (cf-libs) with blackbody radiation reference,” Anal. Chim. Acta 1058, 39–47 (2019). [CrossRef]  

25. F. O. Bredice, H. O. D. Rocco, H. M. Sobral, M. Villagrán-Muniz, and V. Palleschi, “A new method for determination of self-absorption coefficients of emission lines in laser-induced breakdown spectroscopy experiments,” Appl. Spectrosc. 64(3), 320–323 (2010). [CrossRef]  

26. A. El Sherbini, T. M. El Sherbini, H. Hegazy, G. Cristoforetti, S. Legnaioli, V. Palleschi, L. Pardini, A. Salvetti, and E. Tognoni, “Evaluation of self-absorption coefficients of aluminum emission lines in laser-induced breakdown spectroscopy measurements,” Spectrochim. Acta, Part B 60(12), 1573–1579 (2005). [CrossRef]  

27. J. Li, Y. Tang, Z. Hao, N. Zhao, X. Yang, H. Yu, L. Guo, X. Li, X. Zeng, and Y. Lu, “Evaluation of the self-absorption reduction of minor elements in laser-induced breakdown spectroscopy assisted with laser-stimulated absorption,” J. Anal. At. Spectrom. 32(11), 2189–2193 (2017). [CrossRef]  

28. J. Aguilera, J. Bengoechea, and C. Aragón, “Curves of growth of spectral lines emitted by a laser-induced plasma: Influence of the temporal evolution and spatial inhomogeneity of the plasma,” Spectrochim. Acta, Part B 58(2), 221–237 (2003). [CrossRef]  

29. I. Gornushkin, J. Anzano, L. King, B. Smith, N. Omenetto, and J. Winefordner, “Curve of growth methodology applied to laser-induced plasma emission spectroscopy,” Spectrochim. Acta, Part B 54(3-4), 491–503 (1999). [CrossRef]  

30. G. Cristoforetti and E. Tognoni, “Calculation of elemental columnar density from self-absorbed lines in laser-induced breakdown spectroscopy: A resource for quantitative analysis,” Spectrochim. Acta, Part B 79-80, 63–71 (2013). [CrossRef]  

31. Z. Hu, F. Chen, D. Zhang, Y. Chu, W. Wang, Y. Tang, and L. Guo, “A method for improving the accuracy of calibration-free laser-induced breakdown spectroscopy by exploiting self-absorption,” Anal. Chim. Acta 1183, 339008 (2021). [CrossRef]  

32. C. Aragón and J. Aguilera, “Csigma graphs: A new approach for plasma characterization in laser-induced breakdown spectroscopy,” J. Quant. Spectrosc. Radiat. Transfer 149, 90–102 (2014). [CrossRef]  

33. C. Aragon and J. A. Aguilera, “Direct analysis of aluminum alloys by csigma laser-induced breakdown spectroscopy,” Analytica Chimica Acta 1009, 12–19 (2018). [CrossRef]  

34. F. Rezaei, G. Cristoforetti, E. Tognoni, S. Legnaioli, V. Palleschi, and A. Safi, “A review of the current analytical approaches for evaluating, compensating and exploiting self-absorption in laser induced breakdown spectroscopy,” Spectrochim. Acta, Part B 169, 105878 (2020). [CrossRef]  

35. Y. Tang, L. Guo, S. Tang, Y. Chu, Q. Zeng, X. Zeng, J. Duan, and Y. Lu, “Determination of potassium in ceramic raw materials using laser-induced breakdown spectroscopy combined with profile fitting,” Appl. Opt. 57(22), 6451–6455 (2018). [CrossRef]  

36. A. Lesage, “Experimental stark widths and shifts for spectral lines of neutral and ionized atoms a critical review of selected data for the period 2001–2007,” New Astron. Rev. 52(11-12), 471–535 (2009). [CrossRef]  

37. A. El Sherbini, H. Hegazy, and T. M. El Sherbini, “Measurement of electron density utilizing the hα-line from laser produced plasma in air,” Spectrochim. Acta, Part B 61(5), 532–539 (2006). [CrossRef]  

38. J. Ashkenazy, R. Kipper, and M. Caner, “Spectroscopic measurements of electron density of capillary plasma based on stark broadening of hydrogen lines,” Phys. Rev. A 43(10), 5568–5574 (1991). [CrossRef]  

39. H. R. Griem and W. L. Barr, “Spectral line broadening by plasmas,” IEEE Trans. Plasma Sci. 3(4), 227 (1975). [CrossRef]  

Data availability

Data underlying the results presented in this paper are not publicly available at this time but may be obtained from the authors upon reasonable request.

Cited By

Optica participates in Crossref's Cited-By Linking service. Citing articles from Optica Publishing Group journals and other participating publishers are listed here.

Alert me when this article is cited.


Figures (9)

Fig. 1.
Fig. 1. The flow diagram of SSC-LIBS with SAC
Fig. 2.
Fig. 2. The diagram of the LIBS setup
Fig. 3.
Fig. 3. Spectra of aluminum bronze samples in the spectral range of 200-410 nm.
Fig. 4.
Fig. 4. The fitting curve of five spectral lines in the spectra of sample No.3: (a) – (e) for Cu II 204.38 nm, Fe I 382.04 nm, Al I 396.15 nm, Mn I 404.14 nm and H I 656.28 nm.
Fig. 5.
Fig. 5. Comparison of the concentration calculated by SSC-LIBS with SAC (red) and by SSC-LIBS (blue) vs. the certified concentration in brass samples: (a) Mn and (b) Fe.
Fig. 6.
Fig. 6. Comparison of the ARE calculated by SSC-LIBS before SAC and after SAC for aluminum bronze samples: (a) – (e) for No.1-5 sample as the calibration sample.
Fig. 7.
Fig. 7. Comparison of the RMSEP calculated by SSC-LIBS before SAC and after SAC for aluminum bronze samples: (a) – (e) for No.1-5 sample as the calibration sample.
Fig. 8.
Fig. 8. Part of Spectrum of pressed ore samples in the spectral range of 235-780 nm.
Fig. 9.
Fig. 9. Comparison of the ARE and RMSEP calculated by SSC-LIBS before (black) SAC and after (red) SAC for pressed ore samples: (a) for ARE, (b) for (RMSEP)

Tables (4)

Tables Icon

Table 1. The certified concentration of the elements in aluminum bronze standard samples(wt.%)

Tables Icon

Table 2. The certified concentration of the oxides in pressed ore samples(wt.%)

Tables Icon

Table 3. The K value of element spectral lines of various sample

Tables Icon

Table 4. RMSEP (wt.%) and ARE (%) of elements in aluminum bronze samples through SSC-LIBS and SSC-LIBS with SAC

Equations (19)

Equations on this page are rendered with MathJax. Learn more.

S A = ( Δ λ Δ λ 0 ) 1 / α
Δ λ 0 = 2 ω s n e n e r e f
n e = Δ λ S A F n e r e f 2 ω S A F
S A = ( Δ λ Δ λ S A F ) 1 / α ( ω S A F ω s ) 1 / α
K = ( Δ λ Δ λ S A F ) 1 / α
M = ( ω S A F ω s ) 1 / α
S A = K M
C j t = C j s I j t 0 I j s 0 / i = 1 n C i s I i t 0 I i s 0
I 0 = I S A
I j t 0 I j s 0 = I j t / ( K j t M j t ) I j s / ( K j s M j s )
I j t 0 I j s 0 = I j t / K j t I j s / K j s
I i t 0 I i s 0 = I i t / K i t I i s / K i s
C j t = C j s I j t / K j t I j s / K j s / i = 1 n C i s I i t / K i t I i s / K i s
n e = 8.02 × 10 12 ( Δ λ H α α 1 / 2 ) 3 / 2
S A = Δ λ 1 / α Δ λ H α 3 / 2 α ( n e r e f ( α 1 / 2 ) 3 / 2 1.604 × 10 13 ω s ) 1 / α
K = Δ λ 1 / α Δ λ H α 3 / 2 α
M = ( n e r e f ( α 1 / 2 ) 3 / 2 1.604 × 10 13 ω s ) 1 / α
R M S E P = i = 1 n ( x ^ i x i ) 2 n
A R E ( % ) = 100 n i = 1 n | x ^ i x i | x i
Select as filters


Select Topics Cancel
© Copyright 2024 | Optica Publishing Group. All rights reserved, including rights for text and data mining and training of artificial technologies or similar technologies.