Abstract

We describe our results obtained from stoichiometric ratio studies of three different energetic, inorganic samples (ammonium perchlorate (AP), boron potassium nitrate (BPN), and ammonium nitrate (AN)) using the technique of laser-induced breakdown spectroscopy (LIBS) with nanosecond pulses. Signal collection was independently executed using both gated and nongated spectrometers. The oxygen peak at 777.31 nm (O) and nitrogen peaks at 742.50 nm (N1), 744.34 nm (N2), and 746.91 nm (N3) were used for evaluating the O/N ratios. Temporal analysis of plasma parameters and ratios was carried out for the gated data. O/N1, O/N2, and O/N3 ratios retrieved from the gated AP data were in excellent agreement with the actual stoichiometry. In the case of gated BPN data, O/N2 and O/N3 ratios were in good agreement. The stoichiometry results obtained with nongated spectrometer, although less accurate than that obtained with gated spectrometer, suggest that it can be used in applications where fair accuracy is sufficient. Our results strongly indicate that non-gated LIBS technique is worthwhile in the kind of applications where precision classification is not required.

1. Introduction

Nanosecond (ns) laser-induced breakdown spectroscopy (LIBS) has been extensively employed in elemental analysis, impurity detection [1, 2], and identification/classification of materials [3, 4]. In recent years this technique has been applied to study a variety of materials such as alloys, biological samples, pharmaceuticals, and explosives amongst various others [58]. The strength of this technique emanates from prerequisites that an ideal analytical tool should possess—minimal sample preparation, simultaneous multielement detection, ability to analyze all three states of matter, the necessity of short-time scales involved for measurement and analysis, and most importantly the capability of remote detection. In most of the cases requirement for sample preparation is negligible. This technique involves the interaction of an intense laser pulse with a sample leading to plasma generation [912]. The material absorbs light through multiphoton/inverse Bremsstrahlung process that heats the sample surface further to which electrons acquire energy, and as the energy absorbed grows, it leads to collisions and breakdown. Breakdown is followed by the plasma expansion and consequently the creation of shock waves. The plasma expansion rate is observed to be highest towards the focusing lens [13]. During the cooling stage constituent ions/electrons of expanding plasma recombine to form neutral atoms, and some may combine to form molecules. The energy difference appears as light and heat. The light is collected through an optical fiber connected to spectrometer, and is known as the LIBS signal. The LIBS signal evolves rapidly in time and can be divided into two regimes. In the first few hundreds of nanoseconds, the signal is dominated by a broad white light Continuum with sharp peaks overriding it. The origin of the continuum is attributed to Bremsstrahlung and the recombination radiation. The sharp peaks seen in LIBS spectra are signatures of atomic/molecular transitions. After about few hundreds of ns, the Continuum dies down and signal consists of only sharp peaks, which are the signatures of species present in the sample. In order to achieve an excellent signal to background ratio, the data is generally collected with a few 100 ns delay with respect to the exciting laser pulse. However, there have been few reports where a nongated detection scheme was also employed [2, 1416] for recording the spectra. The use of a nongated spectrometer is enticing since it can provide rapid means of studying the stoichiometry of various compounds and thereby eliminating unwanted compounds during the detection of hazardous materials. Though the signal to background ratio is less in the case of nongated detection, compared to that of gated scheme, its use in instrumentation can help cut down the operating expenses drastically. Nongated technique is very easy to operate for a nontechnical person since the light collection part of the experimental scheme is very simple.

High energy materials (HEM’s) are soft organic compounds with the general formula C 𝛼 H 𝛽 N 𝛾 O 𝜈 . The dominant peaks in the LIBS spectra of all such compounds comprise of carbon, hydrogen, oxygen, and nitrogen. It is, therefore, extremely difficult to distinguish between HEM’s and other simple organic materials simply based on the LIBS data. In spite of such an intricacy, Rai et al. [17, 18] in a recent study, with nitro molecules, demonstrated that organic nitrocompounds 4-nitroaniline and 4-nitrotoluene can be identified through their stoichiometry by simply measuring the O/N ratios from their corresponding LIBS spectral data. It has been established that the use of LIBS for identification of compounds on an operational level requires extensive developments in two directions (a) determination of stoichiometry of a particular compound by using intensity ratios of spectral lines of one or more pairs of elements [1820] and (b) matching the unknown spectrum with predetermined and assembled spectral library of reference materials of interest. This particular method was employed by de Lucia and Gottfried [21] along with other elemental ratios calculation as part of an algorithm to discriminate explosives and achieved good success. Several other LIBS studies involving stoichiometric approaches have also been reported recently in the literature [2228]. There have been numerous other efforts over the last decade on LIBS studies of HEM’s and other hazardous materials for identification/classification [2934]. In this paper, LIBS studies on inorganic compounds such as ammonium perchlorate (AP, NH4ClO4), boron potassium nitrate mixture (BPN, BKNO3), and ammonium nitrate (AN, NH4NO3) are reported wherein the comparison of O/N ratios was evaluated in order to evaluate the performance of nongated detection scheme vis-à-vis a gated detection method. Based on these ratios one could easily perform an initial screening of materials for differentiating nitrogen-containing (hazardous) compounds from other organic/inorganic (safe) compounds. After successful initial screening, one could further develop powerful algorithms based on other atomic/molecular line intensity ratios for the fool-proof discrimination of HEM’s.

2. Experimental Details

AP (NH4ClO4), BPN (BKNO3), and AN (NH4NO3) mixtures are principal raw materials for composite propellants used in modern rocket engines as oxidizers. They become explosive when mixed with finely divided organic materials. Rapid and simple techniques for identification of such molecules are essential to prevent untoward incidents. Circular pellets of ~20 mm diameter were prepared by filling the sample in a die and then applying a pressure of eight tones for a period of twenty minutes. The samples were mounted on a linear X-Y translation motorized stage to ensure that subsequent laser pulse hits a fresh portion of the sample. Laser pulses with energies in the range of 25–45 mJ from a second harmonic of Nd : YAG laser at 532 nm (7 ns pulse width, 10/5 Hz repetition rate) were focused on to the sample using an 80 mm convex lens. The sample surface was kept at a distance of 1 mm before the focus such that the peak intensities used were below the breakdown threshold of air thereby avoiding the breakdown of ambient air and prevents its contribution to the LIBS signal. In spite of these experimental conditions we could expect some contribution from air to the total LIBS signal thereby affecting the O/N ratios and therefore we purged the experimental setup (sample region, especially) with dry Argon gas. We have performed the experiments with nongated spectrometer in argon atmosphere and gated spectrometer in ambient air. The signal was collected using a lens system and was coupled into a fiber attached to a spectrometer. The nongated spectrometer (Ocean Optics, MAYA 2000, spectral sensitivity 200–1100 nm range) has a spectral resolution of ~1.2 nm while the gated spectrometer (ANDOR, Mechelle spectrograph ME5000, coupled to an ANDOR iSTARDH734ICCD, spectral sensitivity 200 nm–900 nm range, and resolving power of 5000) had ~0.15 nm. The spectral resolution was measured as full-width half maximum (FWHM) of a Helium-Neon laser line at 632.8 nm. Figure 1 shows the typical experimental setup for LIBS collection using a gated spectrometer. The gated spectrometer was replaced with MAYA spectrometer for the nongated detection scheme. In the gated case, spectra were recorded with different gate delays in the range of 100 ns–1900 ns. The Pockel cell’s output of the laser triggered a SRS delay generator, which in turn provided a TTL pulse to trigger the ICCD.

3. Results and Discussion

Figure 2 illustrates the LIBS spectra of AP, BPN, and AN samples accumulated over ten acquisitions collected with gated spectrometer with 1000 ns gate width and 1000 ns gate delay. Figure 3 demonstrates the spectrum obtained with nongated spectrometer collected with an integration time of 200 msec. The spectral region (700 nm–900 nm) relevant to the stoichiometric calculations is only shown in the figure. The spectra recorded with nongated spectrometers demonstrated a continuum background as the collected signal was a result of 200 ms accumulation time. The continuum background was subtracted from the data using a simple Matlab program. The peaks have been assigned to different neutral and ionic states using the NIST database [35]. In particular, all the spectra exhibited peaks corresponding to nitrogen, oxygen and hydrogen with the both the spectrometers. Observed peaks and their wavelengths are tabulated in Table 1. Additionally, peaks corresponding to transitions within excited nitrogen (396.84 nm, 399.50 nm, 500.20 nm, and 567.16 nm) were observed with the nongated spectrometers and with the gated spectrometer for smaller delays. The peak widths were much sharper in case of gated detection due to higher resolving power. Moreover, the gated spectrometer enabled the recording of time evolution of LIBS signal. A typical time evolution data is presented in Figure 4 for the sample AP using a gate width of 200 ns and for various delays in the range of 100 ns to 1900 ns. During the first five hundred ns of detection, the signal was dominated by Continuum and sharp ionic peaks. Subsequently, the continuum subsided, and signal was dominated by atomic peaks. The time evolution of oxygen peak (at 777.31 nm) and Nitrogen peak (at 746.91 nm) for three samples is shown in Figure 5. The decay has been approximated to a single exponential. The decay constants of both these lines (O-777.31 and N-746.91 nm) were found to be ~0.36–0.44  𝜇 s and ~0.23–0.29  𝜇 s, respectively, for all the samples. As the evolution of both species followed similar dynamics, no corrections were needed when their ratios were evaluated. A drastic difference in the decay constants would lead to incorrect values of the ratios, and corrections need to be incorporated in such a scenario. We have calculated the intensity decay constants of other peaks of Oxygen, Nitrogen, and Hydrogen for all the three samples and data is summarized in Table 2. The estimation of plasma temperature was performed considering the Nitrogen peaks observed in the LIBS spectrum with gated spectrometer. Figure 6 shows a typical Saha-Boltzmann plot obtained for AP with a gate delay of 1300 ns and 200 ns gate width. The plasma temperature was estimated to be ~5173°K. The temperatures estimated for all the three samples with variation in time delay are shown in Figure 7(a). Estimation of the electron density was done with Oxygen peak at 777.31 nm. Typical values of densities obtained were ~1017 cm−3. Figure 7(b) shows the variation of electron density with delay time for all the three samples. The procedure for calculating electron density was followed using the procedure from references [10, 11]. The time evolution of plasma parameters (temperature and electron density) indicates that the cooling time and electron density decreased with time, which is clearly depicted in Figure 7.

The Oxygen peak at 777.31 nm (O) and Nitrogen peaks at 742.50 nm (N1), 744.34 nm (N2), and 746.91 nm (N3) were used for evaluating the O/N ratios. While the Oxygen lines were not resolved, the N triplet was partially resolved with the nongated spectrometers and fully resolved with the gated spectrometer. With the nongated spectrometer, we could clearly resolve only N2 and N3 peaks (N1 peak had less number of points and was discarded for further analysis) and, therefore, were used for fitting and extracting the area under the curves. Typical plots pertaining to Lorentzian fits of the peaks of Oxygen and Nitrogen, obtained with nongated spectrometer, are depicted in Figures 8(a) and 8(b), respectively. The corresponding data for gated spectrometer is presented in Figures 8(c) and 8(d). Three different O/N ratios were calculated for the gated data using three peaks of nitrogen [O/N1, O/N2 and O/N3]. A direct evaluation of the O/N intensities by considering the area under the peaks resulted in very large values. The intensity of the specific atomic peak is dependent on the several parameters such as number densities, temperature, partition function, transition probability, and statistical weight corresponding to that line. Consequently, the parameters like transition probability and statistical weight were included in ratio calculations. This is comprehensible since the density of species is represented by 𝐼 𝑖 𝑗 / A 𝑖 𝑗 g 𝑖 , where 𝐼 𝑖 𝑗 is the observed area under the peak. Since the O line was not resolved fully, we divided the observed intensity into three parts based on the counts as reported in NIST database. The ratio of 𝐼 𝑖 𝑗 / A 𝑖 𝑗 g 𝑖 corresponding to O has then been, totaled and 𝐼 𝑖 𝑗 / A 𝑖 𝑗 g 𝑖 ratio value was obtained individually for N1, N2, and N3. This totaled O value was used to obtain O/N1, O/N2, and O/N3 ratios. Table 3 summarizes the ratios obtained using the gated and nongated spectrometer.

We have calculated all the three O/N ratios with variation of time delay. Figure 9 depicts the time evolution of ratios for all the samples with gate delay in the case of gated LIBS. The ratios clearly showed an increased deviation in the initial part of delays and were closer to the actual stoichiometric ratio at gate delays beyond 1100 ns and up to 1900 ns. Line emission will be dominant and Continuum will die at these time scales. Moreover, the spectra which we obtain during this period are assumed to be in local thermal equilibrium (LTE). The ratios calculated for these compounds after 1  𝜇 s delay were clearly much closer to the stoichiometric lines of respective compounds (dashed lines in Figure 9). We have performed several trials and collected spectra with 1  𝜇 s time delay and 1  𝜇 s gate width for all the three samples. Table 3 shows the stoichiometric calculations for three compounds by considering the data from several independent spectra. The calculated mean ratio of all trails with standard deviation and %RSD values are reviewed in Table 3. The mean ratio obtained for O/N3 was found to be 3.90 in AP, 3.14 in BPN, and 2.84 in AN for the data obtained with gated spectrometer. Out of all the samples AP and BPN stoichiometric ratios matched very well with their actual stoichiometry. Ammonium nitrate is highly hygroscopic and, therefore, will absorb the moisture from atmosphere very quickly. This could, probably, be the reason for the high ratio obtained (actual ratio is 1.5). For the sample AP the three ratios O/N1, O/N2, and O/N3 (3.82, 3.70, and 3.90) were in excellent agreement with the actual stoichiometric ratio of 4, for BPN O/N2 (3.09) and O/N3 (3.14) were in good agreement with actual stoichiometric ratio of 3. That the ratio O/N1 does not match with the stoichiometry needs to investigated. All the three ratios had minimal standard deviation and %RSD. Therefore, we strongly feel that using this ratiometric method one could distinguish these samples at a preliminary level. Similar ratiometric approach was carried out with spectra collected with nongated spectrometer also. In case of nongated detection O/N3 provided a close agreement with actual stoichiometric ratio for AP and BPN. The Nitrogen triplet peaks were not resolved fully, and hence we could not calculate all the 3 ratios. Only two ratios (O/N2 and O/N3) were calculated. The calculated O/N3ratios found were 4.8 for AP, 3.44 for BPN, and 3.22 for AN. The standard deviation obtained was 0.94, 1.18, and 0.74 while %RSD was 19.63, 34.18, and 22.85, respectively, for AP, BPN, and AN. That the O/N2 ratio again did not match with actual stoichiometry needs to be investigated further and could mainly be attributed to the poorly resolved Oxygen and Nitrogen peaks. The evaluated standard deviation and %RSD values were less for gated data compared to nongated data.

Previous work [17] on nitrocompounds working on similar lines evaluated nine different ratios of O/N among which only two were matching with the actual stoichiometry. In the present case we could fairly match the O/N3 ratio even with nongated, data and at least two out of three ratios were in good agreement from the data obtained with gated technique. We feel that these results point a direction towards initial discrimination of such compounds. Combined with statistics-based correlation methods, one could achieve satisfactory and fool-proof classification schemes. As pointed out in earlier works [17], statistical techniques of linear correlation and principal component analysis for the identification of compounds could be included to improve the level of reliability provided by this data. We are in the process of recording the LIBS data of various other nitro-based molecules, using a higher-resolution spectrometer, in an attempt to strengthen our arguments for the case of using nongated spectrometer in such analyses.

4. Conclusions

LIBS studies have been performed on AP, AN, and BPN to weigh the prospects of nongated collection scheme against the gated light collection scheme. Two of these were nitro-based compounds while the third one was simple inorganic molecule. O/N ratios have been calculated from the LIBS spectra. In the case of data obtained with gated spectrometer for AP, all the three O/N ratios demonstrated a good agreement with the actual stoichiometric (O/N) ratio while for BPN two of the O/N ratios matched with the actual values. We observed that, when the spectra was recorded after a delay time of 1000 ns, the calculated O/N ratios were in good concurrence with actual stoichiometry values. In the nongated scheme the calculated O/N3 ratio values were in fine agreement with the actual stoichiometry of sample. The work successfully demonstrates the merits of using simple, yet cost effective, nongated spectrometers compared to the expensive gated spectrometer detection.

Acknowledgments

The authors acknowledge the Defense Research and Development Organization, India for providing the financial support needed for the experimental facility. They acknowledge the samples provided by HEMRL, Pune, India, and Dr. A. A. Vargeese.