Abstract

We present a study of the plasma generated by transversely excited atmospheric (TEA) CO2 laser irradiation of a basalt sample. The plasma was induced in air at atmospheric pressure. The same sample was also analyzed using a commercial LIBS system based on Nd:YAG laser and time-gated detection. The main plasma parameters, temperature, and electron number density were determined and analytical capabilities of the two systems compared. Despite differences in laser wavelength, pulse duration, applied fluence, and signal detection scheme, the two systems are comparable in terms of element detectability and limits of detection. In both cases, all elements usually present in geological samples were identified. The estimated limits of detection for most elements were below 100 ppm, while for Cu, Cr, and Sr they were around or below 10 ppm. The obtained results led to the conclusion that simple, cost-effective TEA CO2 LIBS system can find applications for geological explorations.

1. Introduction

The data on the chemical composition of geomaterials are of fundamental importance for exploitation and use of georesources, environmental studies, and space explorations. Because of that, development of rapid and accurate instrumental analytical methods for chemical analysis of geological materials is of prime importance. At the same time, the analysis of geomaterials is often a challenging analytical task, as it requires determination of a large number of elements, over a wide range of concentrations, in a variety of complex matrices.

The well-established spectrochemical techniques that are mainly used in this field are inductively coupled plasma optical emission spectroscopy (ICP-OES), inductively coupled plasma mass spectrometry (ICP-MS), and X-ray fluorescence (XRF) spectroscopy. In recent years, a considerable progress has been made in the development of another spectrochemical technique, laser-induced breakdown spectroscopy (LIBS) [1]. LIBS combines capability of providing fast multielemental analysis with no sample pretreatment, potential for in situ and remote analysis, the ability to perform surface and depth profiling, and the ability to provide isotopic ratio information additionally to elemental composition [2]. Regarding the analysis of geomaterials, LIBS has certain advantages compared to ICP and XRF. ICP-OES and ICP-MS have better analytical figures of merit such as accuracy and precision and limits of detection compared to LIBS. However, many geological samples are silicate rocks, and the main benefit of LIBS compared to ICP is that there is no need for solid sample digestion, that is, no need for complex and time-consuming dissolution procedures. On the other hand, XRF spectrometry is a rapid method for in-field elemental analysis of geological samples and requires simple and minimal sample preparation (e.g., fusion or preparation of pressed pellets) [3]. These characteristics are also common to LIBS. At the same time, compared to XRF, LIBS has higher sensitivity for elements like Mg, Si, Al, Ca, and K and especially for light elements like Li, Be, B, C, F, and Na. Also, unlike ICP and XRF, LIBS may be easily applied for depth profiling by repeated ablation of the same location. In summary, each of the aforementioned analytical techniques is characterized by its own strengths and weaknesses, and the choice of the “best” technique is dictated by the analytical problem that has to be solved. In that sense, LIBS may be considered as spectrochemical technique complemental to ICP and XRF.

Due to its unique characteristics, LIBS can be used for elemental analysis of almost any type of material. However, the vast majority of LIBS applications deal with compact solid samples and, among them, perhaps the most analyzed sample types were geological materials. LIBS has been successfully applied, both in laboratory setting and in the field, for chemical analysis of rocks, soils, sediments, and other natural materials [2, 510]. One of the main problems in quantitative LIBS analysis of complex matrices (such as geological samples) comes from the matrix effects which can reduce the accuracy and precision of the analysis [1, 2, 11]. A number of calibration strategies were proposed for compensation of matrix effects in LIBS, for example, different normalization procedures as well as calibration-free approach [2, 9].

In this work, a recently developed LIBS system, based on transversely excited atmospheric (TEA) CO2 laser and time-integrated signal detection [12], was applied to the qualitative analysis of basalt samples in air at atmospheric pressure. The basaltic rock was chosen as a representative of geological samples. Basalt is a very important rock as it is the most common rock in Earth’s crust and is also an abundant rock on other satellites and planets, for example, Moon and Mars [13]. TEA CO2 laser-based LIBS has been previously applied for the analysis of basalt under simulated Martian atmospheric conditions (CO2 gas; pressure: 9 mbar) [4]. It was shown that this LIBS system provides good analytical capabilities for geological studies under low-pressure CO2 gas. Although LIBS under different atmospheric conditions offers the possibility to improve the resolution, signal intensity, and signal-to-noise ratio or enable specific applications, it also introduces a degree of complexity in the experimental setup. Thus, the present study was undertaken to examine the possibility of using TEA CO2 laser-based LIBS for the analysis of geologic samples in air at atmospheric pressure. The same sample was also analyzed using a commercial LIBS system based on Nd:YAG laser and time-gated detection.

Currently, the most common LIBS system used for the analysis of geological samples comprises Nd:YAG laser and time-gated detection [510]. Compared to it, the main advantage of TEA CO2 laser-based LIBS is simplification of the instrumentation. The characteristics of TEA CO2 laser (wavelength and pulse duration) make this laser especially suitable for time-integrated measurements, which precludes the use of delay generators and expensive gated detectors required for time-resolved measurements. Differences in laser wavelength and pulse duration for the two laser sources (10.6 μm, 100 ns for TEA CO2 compared to 1.064 nm, 7 ns for Nd:YAG) inevitably influence the ablation process, the plasma parameters, and plasma expansion dynamics. Excitation wavelength has a strong influence on the laser-target and laser-plasma interactions. Mass ablation rate per pulse () is higher for shorter wavelengths () [14]. However, when operating with infrared wavelengths, dominant mechanism for laser-plasma interaction is Inverse Bremsstrahlung (IB) absorption. The rate of IB absorption (αIB) increases as the laser wavelength increases () [14]. Thus, for Nd:YAG plasma, a higher ablation efficiency may be expected, while laser-plasma interaction is expected to be much stronger in the case of CO2 laser. Efficient IB absorption reheats the plasma, causing an additional plasma excitation and expansion. In LIBS, only a fraction of the analyte mass is excited and is capable of producing detectable optical emission. Higher excitation efficiency in CO2 plasma may compensate the lower ablation efficiency. In addition, the dimensions and lifetime of plasma increase with pulse duration. The plasma takes longer to decay and, hence, the emission lasts longer. With the long-pulse irradiation, the excitation of the plasma is rather gradual, and, hence, it is suitable to obtain highly excited plasma.

The analytical capabilities of the two LIBS systems were compared in terms of detectability and limits of detection (LODs). A simplified single-sample method was applied for estimation of LOD [15]. Elemental composition of the basalt sample was determined by ICP analysis. Also, plasma parameters, temperature, and electron number density for plasma created by CO2 and Nd:YAG laser pulses were evaluated and compared.

2. Experimental

2.1. The Sample

The rock sample used in the present work was the same as the one used in our previous study [4], that is, a Tertiary basaltic rock from the Balkan Peninsula, classified according to the TAS scheme as basaltic trachyandesite [16]. A sample was prepared by cutting the original natural basaltic rock to a piece about 35 mm long and ~5 mm thick. The sample was used as is, without any preparation. Elemental composition of the sample was determined by ICP analysis. Samples for ICP analysis were prepared by alkaline fusion procedure [17]. For the major elements, 40 mg of powdered material was fused with LiBO flux and then dissolved in dilute HNO3. For the trace elements, HF/HClO4/HNO3 dissolution was undertaken using 100 mg of sample. Five analyses were completed to determine measurement precision. Elemental composition of the sample is given in Table 1.

2.2. TEA CO2 LIBS System

The experimental setup is described in prior publication [4, 18] and is briefly summarized here. A commercial TEA CO2 laser system, developed at Vinča Institute, was used as the energy source for plasma generation on the basalt samples [19]. Pulsed TEA CO2 laser emits radiation at 10.6 μm, with an output peak power typically in the order of megawatts. The laser-optical pulse has a gain-switched spike, followed by a slowly decaying tail. The full width at half maximum (FWHM) of the spike is about 100 ns, and the tail is about 2 μs. About 35% of the total irradiated laser energy is consisted in the initial spike. Typical output pulse energy was 160 mJ. The sample surface was irradiated by laser light focused using a ZnSe lens (focal length: 130 mm). The angle of incidence of the laser beam with respect to the surface was 90°. The sample was mounted on a rotating sample holder. Single-shot craters well separated in space were achieved by rotation of the sample, while repeated ablation of the same circular path was avoided by lateral movement of the sample.

The optical emission from the plasma was viewed in the direction parallel to the target surface (Figure 1(a)). A time-integrated space-resolved laser-induced plasma spectroscopy (TISR-LIPS) was applied [12, 2023]. Spatial resolution was achieved by changing the position of the plasma along the direction of the incoming laser beam. The focusing lens and the target holder were placed on positioning stages (three-axis translation stage; travel: 25 mm; resolution: 10 μm) to maintain focus during experiments. For spatially resolved measurements, it was essential to find plasma region, where spectral continuum background intensity is low, while the discrete line signal intensity is high. The optimization was performed by monitoring the signal-to-background ratio (SBR) of several elements (Mg: 285.2 nm, Ca: 393.4 nm, Mn: 403.08 nm, and Fe: 358.12 nm lines) as a function of the viewing position. Different spatial regions of the plasma were sequentially examined by moving the target holder in steps of 100 μm. Highest SBR values were obtained at distances 2 mm away from the sample surface. Thus, in all reported experiments, the observing position was set to 2 mm from the target surface. The horizontal part of the plume was projected by an achromatic lens on the entrance slit of the monochromator (entrance slit width: 30 μm; height: 1 mm; magnification: 1 : 1). The monochromator was equipped with diffraction grating with 600 lines/mm (dispersion: 0.7 nm/mm; blaze at 330 nm in the first order). For the time-integrated measurements, Apogee Alta F1007 CCD camera was used. The CCD consists of 1024 × 122 pixels, each of 12 × 12 microns. The total active area is 12.3 × 1.46 mm.

The laser was operated at 1.3 Hz repetition rate, while the shutter of the nongated CCD detector was opened for 30 s. Measurement time was much longer than plasma lifetime; however, the main contribution to the time-integrated spectrum comes from a very limited temporal window. Averaging of many acquired spectra was a procedure used to attenuate the variation of spectral emission, that is, to compensate pulse-to-pulse variations. Each recorded LIBS spectrum is an average spectrum, obtained by accumulation of consecutive spectra from 40 different locations on the sample. This procedure was as a rule repeated in triplicate, and the resulting spectra were averaged. Thus, all reported spectra correspond to an average of 120 laser pulses. A relative standard deviation of measured LIBS signal intensity was 10–15%.

2.3. Nd:YAG LIBS System

Basalt sample was also analyzed by using commercial LIBS system based on pulsed Nd:YAG laser and time-gated detection (Figure 1(b)). This was a LIBSCAN-100 system which was purchased from Applied Photonics Ltd. (UK). As laser source, it is equipped with a Q-switched Nd:YAG laser emitting at 1064 nm (model Quantel Big Sky CFR Ultra GRM), pulse duration of 7 ns, frequency of 10 Hz, and energy up to 100 mJ. The system possesses 8 spectrometer channels covering 182–1057 nm spectral range. The acquisition is made on an integration time of 1.1 ms (the smallest) with 1.27 μs delay time from the laser pulse. The signals were averaged on 50 laser pulses. The statistical variability in the LIBS signal was approximately 10%.

A diameter of a focused laser spot was ~1 mm for TEA CO2 laser and ~1.4 mm for Nd:YAG laser, that is, in both cases well above the grain sizes of the rock sample (<1 mm) [4]. Thus, the influence of grain size on accuracy and precision of the LIBS analyses was not expected to be large [24]. The main contribution to uncertainty comes from fluctuations due to laser energy and high degree of surface roughness.

3. Results and Discussion

LIBS spectra of basalt were recorded in air at atmospheric pressure. All reported spectra are spectra averaged over many laser shots, and the measurements were taken at a number of locations on the sample surface (a fresh surface was exposed to each incident laser pulse) to avoid the problem with sample heterogeneity. Spectra obtained on the basaltic rock sample with a single TEA CO2 laser pulse using a nongated detection, covering the spectral region from 250 to 520 nm, are shown in Figure 2. The spectra were taken 2 mm from the target surface. For comparison, spectra obtained using LIBS system equipped with pulsed Nd:YAG laser and time-gated detection are also shown. Most of the elements present in the sample were identified (Table 1), including those present in concentrations around or below 100 ppm, like Co, Cr, Ni, and V. As it may be seen from Figure 2, there is no difference with regard to element detectability for the two LIBS systems used, despite the fact that the applied laser intensities differ by one order of magnitude (56 MW/cm2 and 460/740 MW/cm2 for TEA CO2 laser and Nd:YAG laser, resp.). Thus, it may be concluded that excitation conditions in CO2 plasma are suitable for detection of all elements usually present in geological samples.

The dissociation, ionization, and excitation processes taking place in plasma are related to characteristic parameters of plasma such as temperatures () and electron number densities (). Thus, it was of interest to determine these parameters for plasma created by CO2 and Nd:YAG laser pulses. The excitation temperature was evaluated by Boltzmann’s plot method [25] using the equationwhere is the wavelength of the emitted light, is the integrated line intensity of the transition involving an upper level () and a lower level (), is statistical weight of upper level, is the transition probability, is Planck’s constant, is the velocity of light in vacuum, is the total number density, is the partition function, is excitation energy, and is the Boltzmann constant. For the determination of excitation temperature of the plasma, eleven iron ionic (Fe II) lines were selected as shown in Figure 3. Considering a relatively narrow spectral interval used, the chromatic sensitivity was not taken into account.

The line identifications and different spectroscopic parameters such as wavelength, transition strength (), energy , and accuracy of transition probability (; ; ) are listed in Table 2 based on reference data [26]. In case of TEA CO2 LIBS, the whole set of spectral lines listed in Table 2 was used to obtain Boltzmann’s plot, whereas in case of Nd:YAG LIBS, lines numbers 3 and 8 were excluded as they were not resolved.

The estimated temperature of CO2 laser-induced plasma was 14300 K, while that of Nd:YAG plasma was 19900 K and 21900 K for 50 mJ and 80 mJ pulses, respectively (Figure 4).

For both LIBS systems, under applied experimental conditions, high temperature plasma was created. Higher temperatures were obtained in Nd:YAG laser-induced plasma. However, it should be noted that estimated plasma temperatures are not directly comparable, because time-integrated measurements for CO2 laser plasma provide average values for the spatially selected plasma zone, while for the Nd:YAG LIBS the values correspond to a specific time window.

The electron number density was determined from the measured Stark width of Al I 396.15 nm line (Figure 5).

Deconvolution of Stark width from the experimental profile was done by applying the following formula [27, 28]:where , , and , , , , , , and are measured, Voigt, Lorentzian, Stark, Gaussian, instrumental, and Doppler width, respectively. For TEA CO2 LIBS, the instrumental width was measured to be  nm, and, for Nd:YAG LIBS,  nm. An approximate formula was used to determine from deconvoluted Stark widths with Stark broadening parameter taken from [29]:

The obtained values are shown in Table 3.

As seen from Table 3, the estimated electron number densities for the two laser systems are similar. However, these values are not directly comparable, as already discussed for the temperature evaluation. Nevertheless, they provide indication of the excitation conditions in plasma.

For both laser systems, spectral data were processed by calculating the integrated peak area of the emission line, the baseline continuum emission intensity, and the root mean square (RMS) noise of the continuum radiation intensity in regions adjacent to the emission line. The signal-to-noise ratio (SNR) was calculated as the integrated peak area A, divided by the width of the peak area , times the absolute value of the RMS noise: [30]. Signal-to-noise (SNR) ratio for the elements under study, obtained for TEA CO2 and Nd:YAG laser-induced basalt plasma, is shown in Table 4.

A simplified single-sample method was applied for estimation of the limit of detection (LOD) [15]. LOD was calculated using the formula , where is a known analyte concentration obtained by ICP measurements (Table 1). Comparison of estimated limits of detection obtained for the two laser systems is shown in Figure 6. As it can be seen, comparable or even better LODs were obtained for TEA CO2 LIBS system. Also, there is no obvious correlation between the energy of the Nd:YAG laser pulse and the estimated LODs.

4. Conclusions

Laser-induced breakdown spectroscopy is based on TEA CO2 laser and has been applied to elemental analysis of basalt sample. The analysis was done in air at atmospheric pressure without any sample treatment. Results of qualitative analysis were compared to those obtained using a commercial LIBS system based on pulsed Nd:YAG laser. It may be concluded that, with regard to detectability and limits of detection, less complex and cost-effective TEA CO2 LIBS system has comparable characteristics to the system based on Nd:YAG laser and time-gated detection.

Conflicts of Interest

The authors declare that there are no conflicts of interest regarding the publication of this paper.

Acknowledgments

This research was supported by (i) the Ministry of Education, Science and Technological Development of the Republic of Serbia through the project “Effects of Laser Radiation on Novel Materials in Their Synthesis, Modifications, and Analysis” (Project no. 172019) and (ii) COST- MP1208 action (“Developing the Physics and the Scientific Community for Inertial Confinement Fusion at the Time of NIF Ignition”). Andrei Stancalie acknowledges financial support of MCI through NUCLEU Programme (ctr. no. PN 1647/2016) and UEFISCDI (ctr. no. 8PM/2010).