Abstract

Laser-induced breakdown spectroscopy (LIBS) has emerged as a promising technique for both quantitative and qualitative analysis of elements in a wide variety of samples. However, conventional LIBS suffers from a high limit of detection (LoD) compared with other analytical techniques. This review briefly discusses several methods that demonstrate the applicability and prospects for trace element detection while lowering the LoD when coupled with LIBS. This review compares the enhancement mechanisms, advantages, and limitations of these techniques. Finally, the recent development and application of LIBS coupled techniques for trace element detection are also discussed for various samples such as metal alloys, biomaterials, rare earth elements, explosives, drinking water, and water bodies.

1. Introduction

Any element with an average concentration below 100 parts per million (<100 ppm) or less than 100 micrograms per gram (<100 μg/gm) is considered a trace element. Trace elements, even present in small concentrations, can be either beneficial or hazardous, and they can considerably modify materials’ characteristics. For example, heavy metals, such as arsenic and cadmium, can have an adverse effect on health even at trace levels. On the other hand, trace elements play a significant role in vital steel properties such as strength, corrosion resistance, hardness, and toughness. Therefore, reliable and robust detection of elements at a trace level is critical in various scientific, industrial, and technological sectors.

Laser-induced breakdown spectroscopy (LIBS) is a novel atomic emission spectroscopic technique for analyzing materials. LIBS uses a high-power pulsed laser to ablate a tiny part of the sample under investigation to generate a weakly ionized plasma on the surface. As this plasma expands and subsequently cools, the characteristic atomic/ionic spectral lines (popularly referred to as LIBS spectral lines) are emitted by the excited atomic/ionic species present in the microplasma. Since the optical emission from the microplasma contains the spectral signature of all the elements present in the sample material, one can quickly determine the elemental composition of the sample under study by analyzing its LIBS spectra. LIBS has the exceptional ability to perform a multielemental real-time analysis, which is not possible by any other conventional methods [15]. However, due to its relatively modest sensitivity, trace element detection by LIBS is still challenging, which is considered a significant limitation of this technique.

To date, several studies have been conducted to improve the sensitivity of the LIBS technique by including either additional experimental techniques or analytical approaches. The first part of this comprehensive review focuses on a detailed discussion of various techniques coupled with LIBS that enable trace element detection at the ppm level or lower. One of the most frequent approaches for increasing the sensitivity of LIBS is utilizing a double-pulse laser. Another popular enhancing strategy is to employ an additional energy source/tunable laser to enhance the plasma to lower the detection limit of a trace element. Examples of this enhancing strategy are: LIBS coupled with laser-induced fluorescence (LIBS-LIF) [6], spark-discharge-assisted laser-induced breakdown spectroscopy (SD-LIBS), LIBS coupled with magnetic/electric field, and surface-enhanced LIBS (SENLIBS). Another effective technique to boost the signal and lower the LoD is combining LIBS with liquid-to-solid phase conversion (LIBS-PT). Examples of LIBS-PT are adsorption coupled LIBS, SENLIBS, and NELIBS, where nanoparticles (NP) and analyte solutions are dried on a substrate [7]. Lastly, LIBS coupled with multivariate analysis is necessary to discriminate compounds such as explosives, alloys, and biomaterials from other materials where identification of multiple trace elements is essential.

In this review, we present different variants of LIBS that were proposed to overcome the challenge faced for trace element detection, with basic principles and succinct descriptions, along with their compelling advantages and inevitable disadvantages.

The second part of the review elaborates on the applications and recent developments of these LIBS coupled approaches in different fields. These techniques have become a prominent analytical approach for the real-time monitoring of metallurgical products and metal alloys. Moreover, these techniques have excellent potential to analyze micronutrients and heavy metals in food and plants and detect rare earth elements (REEs) in geological minerals or ores. In addition, LIBS-PT has been found to have a great potential for monitoring toxic elements such as heavy metals in freshwater and wastewater, which is indispensable to help prevent direct or indirect exposure to these toxic trace elements. Furthermore, LoDs achieved under various experimental conditions in different fields have also been compiled in this review.

2. Different Types of LIBS Systems for Trace Element Detection

Several approaches for enhancing LIBS sensitivity and lowering LoDs have been proposed for trace element detection having an average concentration at parts per million (ppm) or sub-ppm level. These approaches have given unique features to the LIBS technique either by adding a second energy source to traditional LIBS or combining LIBS with other techniques.

2.1. Single-Pulse LIBS

Figure 1 depicts a typical schematic diagram of the experimental setup for single-pulse LIBS. The excitation laser is a Q-switched Nd:YAG laser with a pulse duration of 8 ns, a repetition rate of 10 Hz, and maximum pulse energy of 850 mJ, and an output at the fundamental wavelength of 1064 nm [8, 9]. A harmonic generator employing KDP crystal can generate up to fourth harmonics of the fundamental [10, 11]. The development of sensitive optical detectors such as intensified CCD (ICCD) has tremendously improved LIBS’s ability to determine trace elements. After a software-controlled adjustable time delay, the Nd:YAG laser pulse electrically triggers the intensified CCD (ICCD) camera [12, 13]. This strategy greatly diminishes the strong background formed by the high-temperature plasma and hence offers the LIBS technique to detect atomic/ionic emission lines of trace elements with ease and precession [14].

2.2. Modification of LIBS System for Sensitive Trace Element Detection
2.2.1. Double-Pulse LIBS (DP-LIBS)

Double-pulse LIBS (DP-LIBS) employs two laser pulses from two different laser sources or from the same source to ablate the sample, resulting in enhanced atomic emission intensity and thereby lower LoDs [16, 17].

In collinear double-pulse mode (Figure 2(a)), the first pulse ablates the sample; plasma is induced; and then the second pulse will reheat the plasma, or a portion will be used to create another plasma. Plasma production, plasma temperature, and electron density can be maximized by choosing a suitable interpulse delay resulting in enhanced atomic emission [1820]. Other possible geometrical configurations include orthogonal preablation, orthogonal reheating, and dual pulse crossed beam modes used in DP LIBS to increase the intensities of atomic emission lines [21].

Two angled laser beams overlap at the sample surface in the crossed beam configuration (Figure 2(b)) [22]. In orthogonal reheating DP LIBS configuration (Figure 2(c)), plasma is produced by the first laser pulse propagating orthogonal to the target surface, whereas the resulting plasma is reheated by a successive delayed pulse focused parallel to the target surface [23]. In orthogonal preablation configuration (Figure 2(d)), the first laser pulse is focused parallel to and at some distance above the sample surface, resulting in the formation of air plasma [24]. Then a delayed second laser pulse, orthogonal to the first pulse, is focused on the target surface at normal incidence for ablation [25].

2.2.2. Femtosecond LIBS (fs-LIBS)

After the chirped pulse amplification technique (CPA) was developed, the promise of ultrashort pulses for LIBS analysis has increased manifolds. Typically, a mode-locked Ti:sapphire laser is employed as an fs laser source (Figure 3).

Because of the ultrashort duration of the fs pulse, it deposits all the energy onto the target surface without interacting with the resulting plasma [27], which contrasts with ns laser ablation [28]. Unlike ns laser ablation, fs laser pulses are too short for producing ionization, sample heating, and vaporization until the ending of the laser pulse. The fs laser pulses are much shorter than the heat conduction time of electron and energy transfer time from electron to ion, both of which are in the picosecond range. Therefore, femtosecond LIBS offers significantly reduced heat-affected region and thermal damage due to insignificant hydrodynamic motion and heat conduction throughout the laser pulse [29]. As a consequence, controlled material modification and removal are possible [27]. Other advantages of fs-LIBS over ns LIBS include higher spatial resolution, decreased continuum emission, and reduced plasma mixing from the atmosphere [30]. Employing filaments produced by UV femtosecond laser radiation enables beam propagation over long distances (5–100 m) [31] and has potential for remote LIBS helping rapid elemental analysis of distant targets [32]. On the other hand, the disadvantages of fs-LIBS are the much higher complexity of the fs laser system and significantly higher cost.

2.3. LIBS Coupled with an Additional Energy Source
2.3.1. LIBS Coupled with Laser-Induced Fluorescence (LIBS-LIF)

The LIBS-LIF approach, first proposed by Measures and Kwong, utilizes two lasers. The first laser forms plasma, and the second laser re-excites atoms of the element of interest and eventually induces fluorescence emission from the atoms [33]. The hyphenated LIBS-LIF approach preserves the advantages of LIBS, such as online applications, while offering a much lower LoD, well below the ppm range. The LoD can be minimized by optimizing the ablation and excitation laser pulses’ flux density and optimizing the delay between them [6].

2.3.2. Magnetic-/Electric-Field-Assisted LIBS

The high-intensity magnetic fields help enhance the emission from laser-induced breakdown plasma [34, 35]. Therefore, the energy imparted to the plasma can be well utilized by applying the magnetic field as it influences ionization and the lifetime of plasma and limits plasma expansion [35]. A homogeneous electric field also enhances the intensities of the atomic/ionic emission line and the S/N ratio. Therefore, highly congested emission lines can be identified easily [36, 37].

2.3.3. Spark-Discharge-Assisted LIBS (SD-LIBS)

In spark-discharge-assisted LIBS (SD-LIBS), the strong electric field of the electrical spark causes an electron avalanche, which causes the plasma plume to expand inside an electrode pair [21]. SD-LIBS demonstrates significant enhancement in line intensities as well as considerable improvement of signal-to-noise (S/N) ratios [38] because SD-LIBS increases the electron density of the plasma instantly after the initial ablation of sample material [39].

2.3.4. Microwave-Assisted LIBS (MA-LIBS)

In this technique, the microwave produced from an adjacent microwave cavity is resonantly coupled to the plasma (Figure 4). As a result, unbound electrons are excited, resulting in “hot electrons” that eventually excite more atoms and ions by collision [40]. Hence, plasma lifetime increases (from μs to ms), resulting in enhanced sensitivity [41]. However, remote interrogation of a target becomes complicated as the cavity-generated microwaves must be transmitted over short distances by antennas or waveguides.

2.3.5. Resonance-Enhanced LIBS (RELIBS)

Resonance-enhanced LIBS (RELIBS) is a variant of the two-pulse scheme in which plasma is generated and heated by photo-resonant excitation of the host species in the plume. The excitation wavelength is set to a strong absorption line of one of the main species by an optical parametric oscillator (OPO) laser. Through multistep collision, energy is transferred to the analyte [43]. When the fluence of the initial laser pulse is near the ablation threshold, the LoD improves by the most significant amount compared to ordinary LIBS [44]. RELIBS is exceptionally well suited for minimally destructive analyses. RELIBS has an advantage over LIBS-LIF because it can determine many species simultaneously [45].

2.4. LIBS Coupled with Other Techniques
2.4.1. LIBS-Raman

In a hyphenated LIBS-Raman system, the same laser is used for LIBS and Raman spectroscopy in sequence in a single set-up, simultaneously providing information about the material’s elemental composition and molecular composition [46]. Unlike typical CzernyTurner-based spectrographs, the ICCD is coupled with a broad band-pass (250–900 nm) echelle spectrograph with no moving dispersive elements. In addition, the echelle spectrograph provides high resolution with a fixed and very narrow width (10 × 50 μm) of the slit [46]. The LIBS-Raman approach can be a useful tool for sorting and identification of plastics [47], classification of inks and pigments [48], space exploration [49], analysis of kidney stones [50], and mineralogical and environmental applications [51].

2.4.2. Nanoparticle-Enhanced LIBS (NELIBS)

NELIBS is based on the deposition of metallic NPs with a suitable interparticle distance on the solid (metal alloy and transparent samples) surface or nanoparticle-aqueous solutions (biological fluids and water) dried on a substrate. The interaction with the laser induces the oscillation of conduction electrons coherently and collectively in the metallic NPs. Resonance, termed localized surface plasmon resonance, occurs when the laser’s frequency coincides with the plasma frequency of free electrons in the NPs [52]. In the gap between the NPs and near the particle surface, an order of magnitude increase in the incident laser’s local electric field is observed [53]. This phenomenon can induce a faster breakdown than thermal vaporization as electron extraction from conductive material is much easier in this case [54]. The field enhancement depends on laser power and the concentration of NPs on the substrate’s surface [55].

NELIBS has been advantageous, particularly for studying transparent samples such as glasses and gems. Therefore, laser-sample interaction can be avoided; the NPs absorb the laser energy to initiate the breakdown quickly. NELIBS technique can perform single-shot elemental analysis with LoD below the ppb level [56].

2.4.3. LIBS Coupled with Mathematical Methods

The analytic capabilities of LIBS can be significantly improved with proper spectrum analysis, and based on the spectrum’s characteristics, appropriate mathematical methods should be chosen. Mathematical processing, such as the univariate methods, includes single linear regression [57], internal [58], or external standard [59]. The univariate analysis method establishes a correlation between the spectrum’s intensity and the element’s content. However, due to the spectral line’s self-priming and self-etching effect and interference from other peaks, the fitting with a single variable suffers from poor prediction ability. In contrast, mathematical processing such as multivariate methods is found to be superior to univariate analysis. Multivariate analysis can fully utilize the spectrum to reduce matrix interference, remove redundant information, and extract relevant information in order to develop a quantitative model for improving the sensitivity and prediction ability [60]. Partial least squares (PLS) is the most frequent multivariate analysis (chemometric method) used by LIBS for concentration measurement [61]. Some other examples of multivariate methods include principal component analysis (PCA), artificial neural network (ANN) [62], random forest [63], vector machines regression model [64], multiple linear regression (MLR) [65], soft independent modeling of class analogy (SIMCA), and least absolute shrinkage and selection operator (LASSO) [66].

2.5. LIBS Coupled with Liquid-to-Solid Phase Transformation (LIBS-PT)

In general, liquid-to-solid phase conversion is an effective technique for lowering the LoD. Examples of LIBS in combination with phase transformation (LIBS-PT) are SENLIBS, adsorption-coupled LIBS, NELIBS [7] (where NPs and aqueous samples are dried on the substrate), dispersive liquid-liquid microextraction LIBS [67], electrodeposition of heavy metal on an electrode [68], and metal precipitation followed by membrane separation [69].

2.5.1. LIBS Coupled with Adsorbate

Direct detection of dissolved trace elements by LIBS with adequate sensitivity is challenging in an aqueous solution. Detection sensitivity can be increased, and LoD can be lowered significantly if trace elements are adsorbed by a suitable adsorbent and then implementing LIBS on the composite solid matrix of trace element and adsorbent [15]. Examples of adsorbate used for liquid-to-solid conversion of samples are wood slice [70], bamboo charcoal [71], carbon planchet [72], ZnO [15], zeolite [73], porous electrospun ultrafine fibre [74], 3D nanochannel porous membrane [74], metallic substrate [75], and calcium hydroxide pellet [76]. The calibration curves usually have an initial linear part, and however, it subsequently demonstrates saturation, which suggests these curves to be Langmuir type isotherms [15, 73] (Figure 5). This characteristic saturation limits the highest concentration detected by adsorbent-coupled LIBS measurement of trace elements in liquid for a particular adsorption time and particular adsorbent [15, 73].

2.5.2. Surface-Enhanced LIBS (SENLIBS)

In the case of surface-enhanced LIBS (SENLIBS), a microdroplet of liquid is smeared and dried on the surface of metallic substrates. The use of a metallic substrate aids in the signal enhancing process, as the plasma comprises a combination of species from the metallic substrate and the sample [77]. In SENLIBS, a low boiling point substrate provides higher detection sensitivity [78]. SENLIBS is an effective enhancement technique that has been successfully used in the investigation of oil [79] and heavy metal elements in liquids [80, 81] with improved LOD and reduction of the matrix effect.

In chemical replacement SENLIBS (CR-SENLIBS), the inert and trace heavy metal ions are replaced rapidly by highly chemically active metal ions in the aqueous solution [80]. For example, Pb, Cr, Cd, and Cu in an aqueous solution can be replaced sufficiently by Mg on the magnesium alloy within only 1 min. After that, the concentration of each heavy metal can be found by SENLIBS on a magnesium alloy surface [80].

3. Comparative Study of Different LIBS Systems

Table 1 briefly summarizes the advantages, disadvantages, and key features of different types of LIBS systems discussed in Section 2.

4. Applications of LIBS for Trace Elements Detection

4.1. Analysis of Metals and Metal Alloys

Metal alloy properties, such as strength, corrosion resistance, hardness, and toughness, depend mainly on the minor and trace elements [8789]. Several approaches have been proposed to improve the analytical sensitivity of LIBS; DP-LIBS (Section 2.2.1) is one of them. Jiang et al. [90] explored the combined effect of several ambient environments in terms of gas and pressure with DP-LIBS to enhance the detection sensitivity of the technique. They found that the LoDs of C and S in steel obtained with double-pulse excitation in a vacuum were approximately three times better than those obtained with single-pulse excitation. Cui et al. [91, 92] employed a long-short DP-LIBS to determine Mn and C in the steel. Using this method, the linearity fit (R2) was improved up to 0.988 and 0.9667 for Mn and C, respectively.

The analytical benefits of using double-pulse excitation in the LIBS system to determine minor elements in aluminum alloys were evaluated by Ismail et al. [93]. Spectral lines of Mg, Ti, Cr, Mn, Fe, Ni, and Cu were selected to build the calibrations curves where the Voigt function was used for fitting. Under optimum conditions, DP-LIBS provided significantly improved LoDs within 4–100 ppm for all the examined elements in aluminum. Sun et al. [94] employed spatially resolved double-pulse excitation in the LIBS experiment for the multielement detection of aluminum alloys. Calibration curves of Fe, Cu, Mg, Mn, Zn, Sn, Pb, Ni, Ti, Cr, Ca, and Sr were obtained at six different positions (0.36 to 2.16 nm) of the plasma plume, and the lowest LoDs for each element were found at the center of the plasma (1.08 nm) in DP-LIBS. The LoDs of all the elements of interest were achieved in the range of 0.56 to 27.41 ppm. Hai et al. [95] proposed a 10 μs delay between annular and circular nanosecond laser pulse in DP-LIBS for enhancing signal intensity and reducing LoDs of trace elements Li and Mg in aluminum alloy. In comparison with SP-LIBS (Section 2.1), the proposed DP-LIBS method observed about four times signal augmentation and three times LoD reduction (Figure 6).

Santagata et al. [96] used the femtosecond DP-LIBS technique to trace Pb, Sn, and Zn in copper-based alloys. The emission intensities of the investigated trace elements were normalized with a Cu (I) emission line, and a linear regression coefficient within the range of 0.996 to 0.999 was observed.

Abbasi et al. [97] applied MF-LIBS (Section 2.3.2) to enhance the signal intensity of palladium (Pd) plasma. For MF-LIBS, a 3–4-fold enhancement of optical emission intensity was observed for Pd-I and Pd-II lines compared to LIBS only (Figure 7). It was noticed that the effect of the magnetic field on signal enhancement is element-dependent. Hao et al. [98] introduced LIBS assisted with ring-magnet confinement to improve the detection sensitivity of trace vanadium (V) and manganese (Mn) in steel. In this technique, an external magnetic field influenced the movement of electrons and ions in the plasma plume to increase plasma temperature and electron density. The LoDs of V and Mn in steel were decreased to 11 and 30 ppm by employing a ring magnet, respectively. Al Shuaili et al. [42] utilized MA-LIBS (Section 2.3.4) to detect palladium in solid samples, where 92-fold improvements in the signal intensity were noticed. The LoD for MA-LIBS was 5 ppm, which was found to be 8 times greater than the conventional LIBS.

Some spectral lines of interest such as P (178.3 nm) and S (180.7 nm) are in the vacuum UV range, and it is significantly challenging to detect these lines using an ordinary spectrometer operating in air. However, without using a vacuum chamber, the LoD of trace elements such as phosphorous in steel can still be improved if LIBS-LIF is employed (Section 2.3.1) [99]. The LoD of P in steel was 0.7 ppm, better than conventional LIBS [99]. Li et al. [100] introduced a wavelength-tunable optical parametric oscillator (OPO) laser and an Nd:YAG laser to selectively enhance boron spectral line intensities in nickel-based superalloys to improve the quantitative analysis of trace boron. B I 208.96 nm line intensity was enhanced up to 5.8 times with a 249.77 nm OPO laser. With this LIBS-LIF approach, the LoD of trace boron was 0.9 ppm. LIBS-LIF with a ground-state atom excitation (LIBS-LIFG) has a lower LoD and higher analytical performance than LIBS-LIF with excited-state atom excitation (LIBS-LIFE; Figure 8).

Goueguel et al. [44] introduced RELIBS (Section 2.3.5) approach to investigate the signal-to-noise ratio of Mg I 285.21 nm line in aluminum alloys. They investigated the effect of four optimum parameters such as excitation wavelength, ablation fluences, excitation fluences, and interpulse delay on the SNR of Mg (285.21 nm). The LoDs of Mg and Si in aluminum were 21 and 1,438 ppm, respectively, using RELIBS under optimum conditions. In a subsequent study, the same group applied RELIBS and selective excitation of Al I 309.27 nm to enhance the detection sensitivity of trace Mg and Si in aluminum alloys [45]. In this method, the excitation energy of Al I 309.27 nm was transferred to all the components of the plasma plume, including the trace elements Mg and Si, through particle collisions. Using the optimum operating conditions, the LoDs were 0.75 and 80 ppm for Mg and Si, respectively, far better than the previous experiment. In another approach, Khedr et al. [102] spread gold nanoparticles on aluminum alloy, employed NELIBS, and observed approximately a two/three times improvement of the intensities of the aluminum atomic lines, whereas, in the case of atomic lines of Mg, about eight times signal enhancement is observed.

A high repetition rate SD-LIBS (Section 2.3.3) was developed by Jiang et al. [103] to improve the detection sensitivity of trace elements in copper. A 3–9-fold improvement of the LoDs of Pb, Fe, and Al under this experimental condition was achieved, and the LoDs were 240, 113, and 30 ppm, respectively. A similar technique with lock-in signal detection was employed by Kang et al. [104] to analyze Pb and Al in Brass, and LoDs were found to be 0.112 and 0.178 ppm, respectively.

LIBS coupled with a mathematical method (Section 2.4.3) was evaluated to quantify trace elements in steel [105107]. Sarkar et al. [105] reported that with a precision of about 4–9%, the PLS model demonstrated remarkable success in identifying Mo and Co at the level of hundreds of ppm. The LoD of Cr and Ni was achieved with an accuracy of 2%. Sun et al. [106] demonstrated that the values of LoD for Si, Mn, Cr, Ni, and V in molten steel were 107, 134, 649, 430, and 88 ppm, respectively, using calibration curves obtained with PLS model, with univariate calibration coefficient (R2) greater than 0.99. In 2019, He et al. [108] proposed a genetic algorithm and backpropagation (BP) neural network model coupled with LIBS where carbon was detected in steel using C I 247.86 nm line instead of C I 193.09 nm line. To detect trace elements such as Mn, Cr, and Ni in steel, Zhang et al. [109] developed a new LIBS spectrum data treatment technique based on machine learning algorithms. The multivariate calibration curves of trace elements of interest demonstrated that the values of linearity fit (R2) were within the range of 0.99966 to 0.99997, which was very close to unity. Relative error of prediction (REPs) and relative standard deviation (RSDs) for determining Mn, Cr, and Ni in steel was obtained in the range of 1.3–7.2% and 3.96–6.68%, respectively.

The conditions and LoD obtained in several studies of ferrous and nonferrous metals are summarized in Table 2.

4.2. Analysis of Foods

LIBS is one of the most promising optical tools for food quality control [117]. With increasing concerns placed recently for food quality and safety, LIBS demonstrates excellent potential for agricultural-related sample detection due to its rapid elemental detection capability [118].

Haider et al. [119] utilized SP-LIBS (Section 2.1) to analyze the micronutrients and ecotoxic elements in rice and husk samples from Bangladesh in which micronutrients such as Cu, Fe, Na, Mn, and Zn were found as well as the toxic trace element Cd. Kim et al. [120] applied SP-LIBS to quantify nutrients in unpolished rice, and LoDs of Mg, Ca, Na, and K were 7.54, 1.76, 4.19, and 6.70 ppm, respectively, as established from standard reference materials (SRMs). Luo et al. [121] employed SP-LIBS to analyze the effect of Cu stress on other elements such as Zn, Mn, Fe, Mg, Si, and K in rice husk. The study revealed that Zn, Mn, Fe, Si, and K concentrations are strongly correlated with Cu stress in rice husk.

Pérez-Rodríguez et al. [122] applied SD-LIBS (Section 2.3.3) with the SVM method for finding the origin of botanical rice where C, Ca, Fe, Mg, N, and Na lines were used as input variables for rice variety predictive modeling. Nespeca et al. [123] utilized an SD-LIBS system to analyze impurities in honey.

Sezer et al. [124] used the LIBS system coupled with mathematical methods (Section 2.4.3) to differentiate several types of meat regarding their differences in protein fractions (actin and myosin). Using the PLS model, LoDs for adulteration of beef with chicken and pork meat are 2.84 and 3.89%, respectively. In another paper, Guo et al. [125] applied sample set portioning on the joint x-y distance (SPXY) algorithm in conjunction with the SVM model and LIBS to differentiate five meat species. Zhao et al. [126] demonstrated that LIBS with PLSR modeling has better prediction performance and lower LoD (1079 ppm) than Raman and FT-IR with PLSR. In a similar study, Wang et al. [127] applied LIBS and the chemometrics method to determine the cadmium content in rice roots. Tian et al. [128] used the LIBS system with machine learning to evaluate the phosphorous concentration in seafood, and the LoD was estimated to be around 370 ppm for normalized data. Akın et al. [129] used the LIBS system with PLS regression to determine sorghum levels in flour mixtures consisting of corn and sorghum, and the corresponding LoD was 4.36%. Table 3 compiles different LIBS conditions and corresponding LoDs of trace elements in Food.

4.3. Micro- and Macronutrients Detection in Plants

LIBS coupled with an external electric field (Section 2.3.2) was introduced to enhance the emission intensity of Cr present in the Euphorbia indica plant [141]. The spectral lines of Cr appeared in the presence of the external electric field, while those were absent in LIBS-only spectra (Figure 9). LoD of Cr was found to be 4.8, 4.6, and 3.5 ppm in roots, stem, and leaf, respectively, using CF-LIBS methods. Similar experimental conditions and methods were applied to improve the emission signals of saline metals present in Aerva javanica plants [142]. In the presence of the external electric field line, intensities of Ca, Mg, K, and Na were enhanced by a factor of 7.9, 6.7, 1.3, and 2.97, respectively. The weighted concentration of the investigated elements (Ca, Mg, K, Na, Si, and Li) was determined in the range of 0.7 to 32 wt% and 0.5 to 33 wt% in the root and shoot samples, respectively. The electric field-assisted LIBS was also used to improve the LoD of Cu in the onion (Allium cepa L.) leaves [143], and the LoD of Cu was lowered up to 0.028 ppm.

Samek et al. [144] introduced an fs-LIBS (Section 2.2.2) for spectrochemical analysis of Fe in two leaf samples: dried maize and fresh Cornus stolonifera. Experimental results revealed that Fe concentrations as low as 5 ppm could be detected using this approach.

Seven heavy metals (As, Cd, Cr, Cu, Hg, Pb, and Zn) in Sargassum fusiforme (algae) were examined by Su et al. [145] using LIBS combined with the threshold variables-PLS model. Multivariate calibration method such as partial least square (PLS) is an empirical approach, that is, any knowledge of fundamental physical properties of plasma is not required [105]. These methods have been applied for the determination of nutrient profiles as well as toxic elements in plants [118, 146149]. Nunes et al. [146] compared PLS with univariate calibration to detect nutrients in sugarcane leaves. The PLS method demonstrated better root-mean-square error prediction (RMSEP) values than univariate calibration for all the investigated elements except K and Mn. By the PLS model, LoDs of P, K, Ca, Mg, Mn, Fe, Zn, and B obtained were 0, 210, 80, 120, 6.6, 1.2, and 0.8 ppm, respectively. Likewise, Braga et al. [149] reported that for analysis of different plant materials, LoDs of all the elements of interest (B, Cu, Fe, and Mn) except Zn were improved in the PLS approach and varied from 3 to 12 ppm. Yao et al. [118] used the PLS method to enhance the prediction accuracy of Cd in polluted fresh leafy vegetables. Peng et al. [147] applied LIBS combined with multivariate calibration based on PLS to detect Cr in rice leaves. Relatively high LoDs were achieved for nutrients in plant materials when an fs-LIBS is used in conjunction with PLS [148], and the LoDs of Ca, Mg, P, Cu, Fe, Mn, and Zn were 0.007, 0.02, 0.4, 7, 12, 2, and 80 ppm, respectively. To enhance accuracy in the detection of heavy metals in mulberry leaves, Yang et al. [150] proposed a novel analysis framework consisting of a self-organizing map (SOM), successive projection algorithm (SPA), and uninformative variable elimination (UVE). LIBS was combined with chemometrics for quick and accurate quantitative analysis of heavy metal Cd in rice stems [151] and roots [127]. Table 4 summarizes the different works related to detecting nutrients and toxic elements in plants by LIBS.

4.4. Trace Element Detection in Biomaterials

Roldán et al. [155] applied the combination of LIBS and neural networks to determine the elemental composition of bones for differentiating wild deer bones from different individuals. Detected elements were P, Mg, Ca, Ba, and Sr. Moncayo et al. [156] used LIBS and neural networks (NN) to differentiate peoples by analyzing their bones and teeth fragments collected from the graveyard. In a recent study, Siozos et al. [157] applied LIBS with NN analysis to assign seventeenth- to nineteenth-century archaeological bone remains to individuals using elemental information of bone remains. Gomes et al. [158] used the one-point calibration-laser-induced breakdown spectroscopy (OPC-LIBS) for the elemental analysis of bioactive ceramic (hydroxyapatite), which has human bone-like composition and crystalline structure. In addition, an important ratio () has been studied to know about the final deposition stoichiometry of that material. Martinez and Baudelet [159] used the LIBS system to analyze human nails to measure the concentration of zinc.

Samek et al. [160] applied LIBS with the Mahalanobis distance pattern recognition algorithm to identify unknown tooth samples in real time. Using the two lines of Mg (518.36 nm) and Ca (518.89 nm), the intensity ratio can be used to identify the healthy tooth from carious ones. If , it can be considered a healthy tooth. Ahmed et al. [161] applied LIBS and found that P, Na, and Ca concentrations are higher in healthy portions than carious parts of deciduous teeth (baby teeth).

Rehse et al. [162] applied the LIBS system to study the growth of “Pseudomonas aeruginosa” bacteria in blood and bile. In their work, calcium, magnesium, and sodium were detected and utilized to distinguish the bacteria. Samuels et al. [163] employed the LIBS technique to discriminate the bacterial spores, molds, pollens, and proteins. Dell’Aglio et al. [164] investigated the possible application of NELIBS as a sensing tool where AuNP-protein (such as human serum albumin or cytochrome C) solution is deposited on a metallic titanium substrate. The maximum enhancement of Ti line intensity was observed when nanoparticle-protein corona is formed, and the corresponding number of protein units needed to construct the corona can be estimated.

Haider and Khan [165] performed SP-LIBS on the corals (the skeleton of dead microorganisms), and several elements such as Ca, C, Sr, Na, Mg, Li, Si, Cu, Ti, K, Mn, Zn, Ba, Mo, Br, and Fe were detected (Figure 10). Muhammed Shameem et al. [50] utilized LIBS-Raman (Section 2.4.1) to determine the constituent elements at renal calculi. LIBS has been used by Singh et al. [166] to classify gallbladder stones by using their constituents elements.

4.5. Analysis of Pigments and Other Archeological Samples

Cristoforetti et al. [167] applied the LIBS-Raman technique to identify pottery elements and analyze the molecular structure of the pigments in pottery samples collected from the northern part of Syria. PCA was applied by Osticioli et al. [168] to identify and differentiate between natural and synthetic ultramarine blue pigments using LIBS-Raman (2.4.1). Spectral analyses showed that Al, Na, and Si were found in all the pigment samples. Natural and artificial ultramarine blue pigments were differentiated based on calcium content since Ca was absent in synthetic pigments.

In another paper, Lazic et al. [169] applied LIBS-LIF (Section 2.3.1) to find the pigment’s elemental composition on the surface of Renaissance Umbrian pottery. LIBS has done semiquantitative analysis on all the ceramic layers, and Al, Ca, Cu, Fe, Mg, Na, Si, Sn, and Pb were detected.

Studies revealed that the LIBS technique could successfully identify organic and inorganic materials with high precision and accuracy. Another advantageous feature of LIBS is that each laser pulse can uncover a new part of the paint layer by removing a tiny amount of surface materials [170].

4.6. Detection of Explosive Residues

The handling and transportation of explosives have been shown to generate traces of explosives on surfaces [171]. The first demonstration of detecting explosives using LIBS (Figure 11) was done by De Lucia et al. [172]. Although the sensitivity of LIBS is much less than that of vapor-based techniques, which typically have LoD of the order of ppt [173], the major advantage of LIBS is its ability for real-time, stand-off detection for safety purposes (up to 130 m has been achieved) [174]. However, for stand-off detection of the explosive residue by LIBS, the high power of lasers is hazardous for the human eye. The limitation for LIBS is that the explosive detection is restricted only to the beam path.

Different multivariate chemometric techniques have been used for the discrimination of energetic material, such as PCA, PLS discriminant analysis [175], SIMCA, and neural networks. Among them, the most promising chemometric-based approach was found to be PLS-DA. In general, a “training set” is built from many LIBS spectra, which serves as a standard in chemometric-based approaches. Then the number of dimensions in the data provided in the training sets is decreased. Therefore, intricacy is also reduced, and new data can be matched accurately [176].

DP LIBS demonstrates an increase in the emission signal of explosive residues. In DP LIBS, the first pulse is efficiently absorbed, and a successive second pulse would efficiently reheat the plasma [173, 177]. Gottfried et al. [178] analyzed RDX and a mixture of 36% TNT, 63%RDX, and 1% wax (termed Composition-B) using double-pulse LIBS, and they were able to detect the simulants at 20 m using PCA as a chemometric technique. In another work, Gottfried et al. [175] investigated RDX, Composition-B, and TNT using the same experimental setup. However, in this case, they reported the identification of explosive residues even when mixed with dust.

Another way to improve the detection of explosive residue is fs-LIBS (Section 2.2.2) that minimizes background continuum and atmospheric entrainment, which are experimentally demonstrated by De Lucia et al. [26] using RDX, C-4, and Composition-B explosives residue deposited on aluminum substrates.

4.7. Detection of Rare Earth Elements (REEs)

Most of the reported REEs detection by LIBS are linked with detecting a few elements in lab-prepared or commercially available samples of lanthanide compounds. Bhatt and coworkers [179] performed univariate and multivariate analyses of six REEs (Ce, Eu, Gd, Nd, Sm, and Y) and obtained calibration curves for these REEs. LoDs are found to be 0.098, 0.052, 0.077, 0.047, 0.250, and 0.036%, for Ce, Eu, Gd, Nd, Sm, and Y, respectively. Castro et al. [180] demonstrated the direct determination of Dy, Gd, Nd, Pr, Sm, and Tb in hard disk magnets by SP-LIBS. The concentrations of some REEs such as Ho, Er, and Pr in LiCl + KCl eutectic salt matrices were measured by SP-LIBS through a quartz window atop an airtight sample cell containing argon gas [181]. Martin et al. [182] used SP-LIBS for the detection of Eu, Gd, La, Nd, Pr, Sm, and Y, where a microscope objective was used to focus the laser on the mixture of 5% oxide and 95% graphite matrix, and an innovative collection optics are used to collect the plasma’s emission. In another paper, Tampo et al. [183] used antenna-coupled MA-LIBS and observed 50-fold enhancement in the emission intensity of gadolinium ions. The observed detection limit of the MA-LIBS was found to be 40 ppm.

Whereas most of these works are associated with the detection of REEs in commercially available samples, Abedin et al. [184] reported the first multiple rare earth (Ce, La, Pr, Nd, Y, Yb, Gd, Dy, and Er) detection in a naturally occurring monazite mineral by SP-LIBS (Figure 12). Likewise, Haider and Khan [185] detected multiple rare earth (Ce, La, Pr, Nd, Yb, Gd, Dy, Er, Sm, Eu, Ho, and Y) in an enriched zircon sample by SP-LIBS. The same group also identified trace REEs (Nd, Yb, Sm, Ce, Dy, and Gd) in coal samples collected from coal mines of Bangladesh (Barapukuria) and India (Meghalaya) [186]. LIBS can measure REE contents in geological minerals and ores and provide results comparable to ICP-MS analysis. In natural geological samples, Bhatt et al. [187] employed SP-LIBS to detect and quantify REEs (Ce, La, Nd, Y, Pr, Sm, Eu, Gd, and Dy). For the quantification of Ce, La, and Nd, multivariate analysis was executed by developing PLS regression models. Using SP-LIBS as a drill core scanner, Müller et al. [188] analyzed thin sections of drill core of the geological material Storkwitz carbonatite, and the presence of Ce, La, and Nd were detected.

Bhatt et al. [21] used DP-LIBS (pulses wavelength 1,064 nm) in a collinear configuration to demonstrate enhancement of the signal intensity of Eu, Gd, Pr, and Y using Eu2O3, Gd2O3, Pr2O3, and Y2O3 pellets, respectively. In another work, employing LIBS-LIF, Shen and coworkers [189] have demonstrated the detection of U as low as 462 ppm in glass samples [16].

4.8. Detection of Trace Elements by Liquid-to-Solid Conversion

The trace element concentration can be increased significantly using a suitable adsorbent, and this liquid-to-solid conversion of samples is advantageous for more accessible and inexpensive detection of trace and toxic elements. For example, Haider et al. [15] employed the LIBS combined with adsorption for the sensitive detection of arsenic in water (Figure 13). In this study, two adsorbents, ZnO and charcoal, were employed to improve the emission line’s intensity. As an adsorbent of arsenic, ZnO provided better results than charcoal. A LoDs of 1 ppm was achieved using ZnO, while that of the charcoal was 8 ppm. It was also reported that a concentration of arsenic as low as 0.083 ppm in water was detected using preconcentration of the liquid sample and adsorption techniques.

Chen et al. [70] used wood slices as a liquid absorber in the LIBS experiment for rapid and sensitive detection of toxic heavy elements such as Cr, Cd, Pb, Mn, and Cu in water. In this study, LoD was obtained in the range of 0.029–0.59 ppm, which was 2–3 orders lower than the results obtained by direct LIBS analysis on an aqueous solution. On the other hand, Youli et al. [190] used ordinary printing paper as a liquid absorber. In this experiment, the LoD of Cr and Pd was 0.026 and 0.033 ppm, respectively, which was a slight improvement compared to the previous one where the wood slice was used as an absorber.

To enhance the detection sensitivity of LIBS for Cd content determination in drinking water, Tian et al. [191] employed a water purifier chelating resin ECS60 as the enrichment matrix. Chelating resin has a rough and porous surface, which helps it adsorb heavy metals effectively. A comparison of resin and filter paper as an adsorbent was investigated, and the results confirmed that chelating resin was more efficient than filter paper. The LoD was 0.0036 ppm, which met the sanitation standard of drinking water in China. So LIBS coupled with chelating resin ECS60 could be a safe, low-cost, and rapid method for monitoring drinking water quality. LIBS was utilized to determine boron in sub-ppm concentrations in groundwater employing graphite planchets as solid support [192]. The LoD in groundwater was found to be 0.01 ppm by using the 249.773 nm atomic line for boron.

In contrast to these adsorption methods, Zhao et al. [193] used SP-LIBS coupled with an electrical-deposition method to analyze the trace heavy metal ions in tap water. A rotating cathode system was designed to deposit the trace metals on the surface of the cathode material. An aluminum rod having high purity was selected as a cathode because it can enhance the deposition homogeneity of the metal ions and reduce the generation of gas bubbles in the electrical deposition process. Trace heavy metal ions, such as Cr3+, Mn2+, Cu2+, Zn2+, Cd2+, and Pb2+ were quantitatively analyzed in the tap water samples with an LoD of 1.6 × 10−4 − 1.35 × 10−3 ppm. These were 5–6 orders of magnitude better than the results obtained by direct tap water analysis by LIBS and 2–3 orders better than using wood slice substrates. Chen et al. [75] also used the SP-LIBS technique coupled with an electrical-deposition method to analyze the trace heavy metal ions in water. Cr3+, Mn2+, Cu2+, Zn2+, Cd2+, and Pb2+ trace heavy metal ions were quantitatively analyzed in the tap water samples with an LoD of 8.3 × 10−5 − 5.6 × 10−3 ppm. These were 5‒6 orders better than the results produced by the analysis in water by conventional LIBS.

To lower the LoD of Pb and Cr in drinking water, Ma et al. [78] used SENLIBS (2.5.2), where the effect of four substrates (Zn, Mg, Ni, and Si) was compared. Among the four substrates, Zn showed the best result. LoD was 0.0011 and 0.004 ppm for Cr and Pb, respectively, using a Zn substrate, which was lower than the sanitation standard of the drinking water in China (0.05 ppm for Cr and 0.01 ppm for Pb). The chronology of LoD of Cr and Pb using several substrates is in the sequence of Zn > Mg > Ni > Si, while LoD is correlated with the boiling point of the substrate. The detection sensitivity of heavy metals in water was better for a substrate of a lower boiling point as the amount of ablation of the target is inversely related to the boiling point. So, with the decrease of the boiling point, the content of the matrix increased, resulting in an intense collision between particles and increased electron density and plasma temperature. Consequently, spectral lines of higher intensity and improved LoD were observed. Zhang et al. [194] studied the impact of substrate temperature on the detection sensitivity of SENLIBS technology. Observed LoDs of Cr and Pb were 0.0046 and 0.0317 ppm, respectively, at 25°C. It was improved up to 0.0012 and 0.008 ppm when the substrate temperature reached 200°C.

Yang et al. [80] introduced chemical replacement coupled with the surface-enhancement LIBS (CR-SELIBS) method to overcome the poor detection sensitivity of LIBS in water. In this method, cleaned magnesium alloy was immersed into the aqueous solution containing heavy toxic metals and then dried. As a result of the chemical replacement reaction, Cu, Pb, Cd, and Cr in the liquid sample were enriched on the surface of the magnesium alloy (Figure 14), and an average thickness of 2 μm was reported. LoDs of the toxic metal ions were in the range of 0.016–0.386 ppm.

For improving the detection sensitivity of Cl and S in water solution, a novel technique termed indirect LIBS (ID-LIBS) was introduced by Ma et al. [195]. In ID-LIBS, Cl was indirectly detected by AgCl, that is, from the precipitation reaction of Ag and Cl and subsequent separation by centrifugation, whereas S was detected by BaSO4 precipitation. Tang et al. [196] proposed LIBS assisted with molecular emission for sensitive detection of F and Cl in water solution. The LoDs of F and Cl were 0.38 and 1.03 ppm, respectively.

Ma et al. [197] analyzed the effect of pH on detecting Cr and Cd in wastewater using LIBS-PT (Section 2.5.1). They observed that at low pH, the production of salt floccule occurs when an acid reacts with a metal, thereby reducing the spectral enhancement obtained through LIBS-PT on the surface of a metallic substrate. Ahmed et al. [11] allowed surface water (river water) around Dhaka to settle down in a container, and the accumulated residues at the bottom were dried to form pellets, enabling the identification of toxic trace element Cr by SP-LIBS. In a similar work by Haider et al. [198], Cr, Co, and Ni were detected in water samples collected from different water bodies around Dhaka city in Bangladesh.

Remarkably, several reports have demonstrated LoD in the tens of ppb to a few ppb levels. To achieve a detection limit down to the ppb level, Kiris et al. [199] employed a transversely excited atmospheric (TEA) CO2 laser, and a drop of analyte solutions is dried after deposition on Polytetrafluoroethylene (PTFE) substrate. The LoDs as low as 0.0053, 0.033, 0.32, and 0.12 ppm for Be, Cr, Pb, and Tl are obtained, respectively. However, they observed that the predeposition of copper oxide NP on the substrate, as typically done in NELIBS, does not improve the analytical performance. On the other hand, in the case of an Nd:YAG laser, predeposition of copper oxide NPs enhances the performance, and this NELIBS technique provides LoDs for Be, Cr, Ni, Co, and V to 0.0042, 0.1, 0.4, 0.35, and 0.20 ppm, respectively. Sensitive detection of Cu in water solution was reported by Wang et al. [200], taking advantage of the LIBS-LIF technique and using wood as an absorber. The LoD was found to be five orders lower than that obtained in the direct LIBS technique of liquid matrix, and it was as low as 0.0036 ppm.

The conditions and key findings obtained in several studies on detecting trace elements in water are summarized in Table 5.

5. Conclusion

LIBS-coupled techniques have emerged as a promising and effective tool for the spectrochemical analysis of trace elements, whereas conventional LIBS has poor performance. This review reports the state of the art of different LIBS-coupled techniques, and the recent use of these techniques in the elemental analysis is presented. The outstanding improvement of these techniques’ analytical performance for trace element detection is also highlighted. Modification of LIBS (e.g., using double pulse) or the combination of LIBS with a secondary energy source (laser, microwave, or magnetic/electric field) is an efficient approach to improve the emission signal of plasma plume and can significantly lower the LoDs. Selective excitation in the LIBS-LIF system can solve the complex spectral interferences that impede trace element detection. Integrating LIBS with Raman in a single piece of equipment can improve proficiency in detecting trace elemental content by delivering complementary data. The chemometric-based approaches have proven promising for sample classification of rock, mineral, or alloy by LIBS. However, more robust detection algorithms are still required to enable samples’ quick and accurate classification. A liquid-to-solid phase transformation can significantly boost the trace element concentration, which helps attain a lower LoD and allows for easier identification of trace and hazardous elements by LIBS. For these LIBS-linked approaches, a major challenge is maintaining reproducibility and repeatability for trace element quantification, and more in-depth research needs to be carried out in the future.

Data Availability

No data were used in this study.

Conflicts of Interest

The authors declare that they have no conflicts of interest regarding the publication of this article.