B. Busser

Laser-induced breakdown spectroscopy (LIBS)–based imaging is becoming a promising technique in the panel of spatially resolved elemental approaches. This method has outstanding advantages, such as all-optical instrumentation, fast acquisition speed, operation at ambient atmosphere, and detection limits at the ppm scale for most elements. LIBS-based imaging has an extensive range of applications, including biology, medicine, industry, and geology. In this paper, we report on recent advances in LIBS imaging, focusing on instrumentation, performance, and applicability. Two examples are shown: first, a speleothem with application to paleoclimate studies, and second, a biological organ with implications for toxicology of new drugs based on nanoparticles. Finally, some future directions are discussed. 

	In the past few years, the application of laser-induced breakdown spectroscopy (LIBS) to microscopic elemental imaging has seen significant developments, in terms of both instrumentation and applications (1–4). In LIBS imaging, a series of laser-induced plasma is generated at different positions on the sample surface with a predefined sequence covering the region of interest. Such plasma sources allow specific optical responses resulting from the electronic relaxation of atoms and ions excited by the high plasma temperature to be elicited from the elements constituting the sample. The light emitted by the plasma is collected and analyzed using an optical spectrometer. The elemental “signal” (atomic and ionic emissions) is then extracted from the recorded spectra, and elemental maps can be obtained in a pixel-by-pixel manner (Figure 1) (5,6).	 	

	The main advantage of LIBS certainly lies in its simple implementation. This simplicity is because a single laser pulse simultaneously samples the material by laser ablation, atomizes, and then excites the vaporized mass by heating the plasma plume. This simplicity endows this technique with a series of advantages, including an all-optical design compatible with conventional optical microscopy, operation in ambient atmosphere, and fast operating speed up to kHz. In addition to its tabletop instrumentation, this technique has multielement capabilities, no restrictions in the detection of light elements, detection limits in the range of ppm for most elements, and microscopic-scale resolution (2). It also has quantitative capabilities through the use of standards with known concentrations or an external calibration using different strategies (6–9). All of these advantages make LIBS imaging highly promising, with the potential to become a reference technique in the area of spatially resolved elemental techniques as a complement to gold standard methods, such as laser ablation inductively coupled plasma mass spectrometry (LA-ICP-MS), synchrotron radiation microanalysis, or electron probe microanalysis (EPMA) (10–13).

	Recently, the LIBS community demonstrated the value of this technique for a range of applications, covering geological, industrial, and biomedical domains. For example, the capability of LIBS imaging to characterize complex samples was explored with the analysis of a multiphase hydrothermal ore sample (14). More than 30 elements were identified, with the concentration ranges extending from major to trace elements. In an application to paleoclimate studies, large-scale speleothems and coral samples (>10 cm

) were analyzed, with the aim of revealing the laminar structures related to vegetation or climatic variations (5). Several studies of industrial materials have also been conducted on crystals (7,9) or heterogeneous catalysts (8,15). The possibility of obtaining elemental images of biological tissues was also demonstrated in different biomedical situations (16–18). In this paper, we illustrate recent advances in LIBS imaging in terms of instrumentation and applicability. Two different examples are shown to illustrate the specificities and the current performance of this technique.

	The micro-LIBS system developed in our group is based on a custom microscope (Figure 2). The opto-mechanical structure has been designed with the use of cage and tube systems commercialized by Thorlabs. The device includes a laser injection line, monitoring cameras (sample surface, plasma), and illuminators (white light, laser pointer), as well as up to three collection lines connected to spectrometers with optical fibers (19). The device is also equipped with a Raman system operating at 735 nm, and allowing punctual analysis to be performed. The laser used is a nanosecond Nd:YAG laser working at the fundamental wavelength (1064 nm) with a frequency rate of 100 Hz (Centurion, Quantel). The typical pulse energy is in the range of 1 mJ for most applications. The laser energy is stabilized throughout the experiment by using a servo control loop consisting of a power meter and a computer-controlled attenuator (Figure 2a). The system is equipped with several Czerny-Turner spectrometers (Shamrock 303 and Shamrock 500, Andor technology) coupled to ICCD cameras (Istar, Andor Technology) as well as compact spectrometers (MayaPro and HR2000+ from Ocean Optics), which can be operated simultaneously. In our instrument, the laser beam is focused using a 15x magnification objective (LMM-15X-P01, Thorlabs). The typical crater size is in the range of 7 µm. For all imaging experiments, a laser-induced plasma is generated continuously while scanning the sample surface. The scanning is performed line by line, continuously, in a raster scan mode with the use of motorized x–y–z stages. The operating speed is equivalent to the laser frequency (100 Hz), and a 1-megapixel image is typically obtained in less than 3 h (5). In addition, the measurements are performed at room temperature under ambient pressure in air or in an argon atmosphere. Finally, in-house developed software written in the LabVIEW environment is used to control the entire imaging setup, thereby allowing automatic scanning sequences in any selected region of interest of a given sample with a preset laser pulse energy and spatial resolution.

	For processing the spectral dataset, our group has developed a custom LabVIEW software, named “LasMap.” This software allows the relevant intensities to be extracted from the spectra before building the corresponding elemental images. We generally use the average spectrum (computed with the whole spectral dataset) to obtain an overview of the spectral structure of the sample, allowing the selection of the lines of interest free of interference. Algorithms then extract the net intensity of the lines by taking into account the baseline level to generate a 2D matrix (relative abundance image) (20). This procedure is relatively fast for a single matrix sample, with a megapixel dataset typically processed in less than a few minutes.

	In LIBS-based imaging, the achievable performance (for example, spatial resolution and detection limits) is strongly related to the laser ablation process. To guarantee the best shot-to-shot repeatability, it is important to set up the experiment to avoid any overlap between consecutive laser shots. The lateral resolution is then ultimately governed by the laser-induced damage and the limits of detection (LOD) depend both on the mass of vaporized material (such as the crater volume) as well as the excitation capability of the laser pulse. Therefore, there is a compromise between these parameters. As an example, the evolution of the detection limit of magnesium (Mg) as a function of the crater diameter is shown in Figure 3a. This result was obtained on a reference glass by modifying the laser pulse energy. The LOD shows significant deterioration when reducing the crater size due to the reduction of the ablated mass (plasma density). With larger crater sizes, there is also an increase of the LOD due to the saturation of the detector, that was optimized for detecting weak signals. To give priority to the detection limits, we preferentially use crater sizes in the range of 6 to 10 µm. It is important to mention that laser ablation is a violent process. It is accompanied by different mechanisms, including shock wave formation and thermal diffusion through the sample, that might cause more sample deterioration than the ablation itself. The resolution will then be dependent on the material properties, particularly the hardness. However, resolution <15 µm can easily be obtained for the large majority of materials.

	Figure 3b shows the typical accessible range of LODs, expressed in weight %, obtained in this configuration for most of the elements. Sub-ppm LODs are generally achieved for alkali and alkali-earth metals. The sensitivity remains attractive for other metals, with LODs better than 20 ppm, but strongly deteriorates for organic and halogen elements. It is important to emphasize that such relative LODs are obtained in a single pulse configuration and from an ablated mass in the range of 0.1 nanogram (typical ablated mass from a single laser pulse). Expressed in absolute detection limits (for example, the smallest mass quantity that can be detected), values typically lower than the femtogram can be achieved.

	Thanks to the 100 Hz scanning speed, and because there is no constraint on the size of the sample (such as, for example, no need for an ablation cell), large surfaces can be easily analyzed. Such features may be valuable with geological samples, particularly for the characterization of speleothems. Stalagmites, stalactites, or columns precipitate from seepage waters in caves, and usually consist of calcium carbonate, with calcite being the most common mineral phase. Climate proxies, such as Mg, Sr, or Ba, can be retrieved from these materials, and may indicate changes in climate or environmental conditions that occurred during the growth periods and provide information on past temperature, precipitation, and vegetation changes over the last hundreds of thousands of years (21). For recording paleoclimate proxies, it is interesting to analyze several centimeters, or even more, of a sample while maintaining resolution on the scale of the annual speleothem growth (typical rates of 0.01 to 2  mm per year) (22).

	An example is shown in Figure 4 for a speleothem extracted from a Gallo-Roman copper ore mine situated close to the “Grotte Bleue” in Aude, France (courtesy of D. Genty and P. Cabrol of Les directions régionales de l’Environnement, de l’Aménagement et du Logement [DREAL], in Toulouse, France). At the collection site, there is a strong abundance of copper, which gives a particular blue color to this mineral. This site is also known to be rich in iron, lead, and arsenic. Amphoras discovered on the site suggested Roman use of the cave from the 1st century BC, but it is suspected that the site was also occupied in the Bronze Age (23). In the example shown in Figure 4, a mapping sequence of 2300 x 1600 pixels was performed with a step size of 24 µm, which corresponds to a surface of analysis of more than 20 cm². The optical image of the sample and elemental images of Mg and Pb are shown. These results illustrate that the resolution capabilities of the LIBS imaging method and annual growth layers can be observed. It is worth mentioning that other elements, such as Cu, Ag, Sr, Ba, Al, and Si, were also detected. The interpretation of their results is currently being explored by specialists.

	LIBS-imaging can also be applied on soft biological tissues. Our first work in this area was conducted on a mouse kidney to image the biodistribution of NPs before administration to animals (6,24,25). These studies were followed by several articles reporting biological studies conducted on murine animals, where LIBS was used as a complement to more conventional techniques such as transmission electron microscopy (TEM) or fluorescence microscopy (17,18,26). More recently, studies were also conducted with human specimens of medical interest (16,27). In health infrastructures, the possibility of coupling LIBS imaging with conventional histopathology could be highly valuable for making some medical diagnoses. Because LIBS imaging is compatible with the light microscopy used in pathology laboratories, it may provide important support for uncovering the elemental composition of tissues, particularly in the case of numerous pathologies relating to exposure to metals (such as respiratory diseases, dermatological diseases, and cancers).

	As an example, Figure 5 shows an image obtained for a rat kidney section. The organ was sampled 1 h after intravenous administration of gold nanoparticles (NPs). This experiment was performed with a resolution of 10 µm, and represents 2 million individual pixels. Three elements are represented with a false color scale: Au in yellow, P in blue (the P distribution reflects the histological structure of the tissue), and Fe in red (in this case, the threshold allows for visualizing Fe from hemoglobin, and by extent, blood vessels). As observed in several previous work (6,26), small NPs are mainly eliminated by the renal route with partial retention in the cortex region of the organ (peripheral region), indicating a normal elimination process. More generally, the possibility of imaging NP distribution at the entire organ scale and without any labeling is highly relevant for studies involving the development of new drugs based on metals.

	LIBS-based imaging is a mature technique for various classes of materials. These materials consist of a homogeneous matrix in which it is desired to observe the possible heterogeneity of minor or trace elements (both examples shown above in Figures 4 and 5). The advantage of such a simple material is that matrix effects can be avoided, greatly facilitating the calibration of raw intensities. In addition, since we are interested in minor or trace elements, line saturation phenomena (such as self-absorption and detector saturation) are nonexistent, and the intensity of the measured lines is linearly dependent on the concentration.

	However, difficulties arise when considering complex materials with multiple matrices. These difficulties are inherent to all elemental imaging methods. Regardless of the analytical technique, the interaction between the primary beam and the material will always depend on the sample properties. The recorded signal will then be dependent on the material (nature, hardness, composition, and so on), so it is difficult to generate a quantitative analysis on such complex material, regardless of the technique used. Currently, the classical way to calibrate the intensities is to use reference samples for each phase of the material, knowing that in practice these standards do not always exist. An interesting possibility for the future is to implement a calibration-free approach. Such a method does not require external references because it only relies on the physical analysis of the spectra. By considering a given plasma model and measuring its parameters (such as temperature or electronic density, for example), it is possible to extract the concentrations of the various detected elements, regardless the type of the matrix. Such a method could be highly promising for the development of the technique, and several groups are currently working on its implementation.

	An additional difficulty related to emission spectroscopy comes from the possible complexity of the signal. Some elements, such as transition metals, have a large number of emission lines from the ultraviolet (UV) to near infrared (NIR) range. This is the case in particular for iron or titanium. If one or more of these elements are a matrix constituent, the observed spectra will be spectrally dense, and the desired lines (generally from minor or trace elements) will likely be masked. In this case, the performance in terms of sensitivity will deteriorate.

	Managing the processing of data obtained on complex samples also represents a challenge and limits the fields of possible exploration of the technique. It must be kept in mind that the data to be processed are very large (up to several tens of gigabytes), and the software used cannot be supervised. It is indeed impractical to interpret spectra one by one for a megapixel analysis. A promising possibility that is currently being explored is the use of statistical analysis methods (such as chemometrics). In general, these methods allow the dimensionality of a series of data to be reduced by looking for correlations among the different variables, thus making the information less redundant and easier to handle. We recently published an article detailing the implementation of principal component analysis (28), and several groups in the LIBS community are currently exploring other methods, such as the possibility of using machine learning approaches.

	LIBS imaging is currently experiencing significant development, with an increase in the number of published papers during the last three years. Due to the feasibility studies conducted in laboratories, this approach is becoming increasingly mature, and is currently a good candidate to become a reference approach for spatially resolved elemental characterization in the near future. LIBS has unique specificities, such as the possibility of reaching micrometer-scale resolution with limits of detection down to ppm, operation at ambient conditions with tabletop, and all optical instrumentation. Regarding the acquisition rate, recent work has demonstrated the possibility of access to kHz acquisition rates (29), and higher operating speeds should be achieved very soon.

	The future will tell if LIBS-based imaging can be generalized and used outside academic laboratories. What is certain is that, at some point, its development will have to be supported by large industrial groups. In the meantime, researchers still have some work to explore. In particular, the analysis of complex materials (those composed of several matrices) appears challenging, both for quantification and data processing. Further studies dedicated to these issues need to be conducted, and the use of calibration-free LIBS and advanced chemometric tools represent interesting opportunities.	 

	This work was partially supported by the French region Rhônes Alpes Auvergne (Optolyse, CPER2016) and the French ANR (projects MediLIBS ANR-17-CE18-0028-04 and Imazinc ANR-16-CE29-0024). In addition, we gratefully acknowledge Florian Trichard from Ablatom S.A. S, Nicolo Omenetto from the University of Florida, Ludovic Duponchel from the University of Lille, and Christophe Dujardin from ILM for fruitful discussions.

L. Jolivet, M. Leprince, S. Moncayo, L. Sorbier, C.P. Lienemann, and V. Motto-Ros, 

B. Busser, S. Moncayo, J.-L. Coll, L. Sancey, and V. Motto-Ros, 

J.O. Cáceres, F. Pelascini, V. Motto-Ros, S. Moncayo, F. Trichard, G. Panczer, A. Marín-Roldán, J.A. Cruz, I. Coronado, and J. Martín-Chivelet, 

L. Sancey, V. Motto-Ros, B. Busser, S. Kotb, J.M. Benoit, A. Piednoir, F. Lux, O. Tillement, G. Panczer, and J. Yu, 

G. Alombert-Goget, F. Trichard, H. Li, C. Pezzani, M. Silvestre, N.S. Barthalay, V. Motto-Ros, and K. Lebbou, 

F. Trichard, L. Sorbier, S. Moncayo, Y. Blouët, C.-P. Lienemann, and V. Motto-Ros, 

K. Janssens, W. de Nolf, G. van der Snickt, L. Vincze, B. Vekemans, R. Terzano, and F.E. Brenker, 

K.L. Moore, E. Lombi, F.-J. Zhao, and C.R.M. Grovenor, 

D. Pozebon, G.L. Scheffler, V.L. Dressler, and M.A.G. Nunes, 

F. Trichard, F. Gaulier, J. Barbier, D. Espinat, B. Guichard, C.-P. Lienemann, L. Sorbier, P. Levitz, and V. Motto-Ros, 

V. Motto-Ros, E. Negre, P. Pelascini, G. Panczer, and J. Yu, 

V. Motto-Ros, S. Moncayo, F. Trichard, and F. Pelascini, Spectrochim. Act. B 

G. Munteanu, E. Kammenthaler, J. Mantenant, and C. Rico, 

V. Motto-Ros, L. Sancey, Q. L. Ma, F. Lux, X. S. Bai, X. C. Wang, J. Yu, G. Panczer, and O. Tillement, 

S. Moncayo, F. Trichard, B. Busser, M. Sabatier-Vincent, F. Pelascini, N. Pinel, I. Templier, J. Charles, L. Sancey, and V. Motto-Ros, 

S. Moncayo, L. Duponchel, N. Mousavipak, G. Panczer, F. Trichard, B. Bousquet, F. Pelascini, and V. Motto-Ros, 

K. Rifai, F. Doucet, L. Özcan, and F. Vidal, 

are with the Institut Lumière Matière at the Université de Lyon, in Villeurbanne, France. 

are with the Institute for Advanced Biosciences at the Univ. Grenoble Alpes, in Grenoble, France. 

is also with Grenoble University Hospital, in Grenoble, France. 

is with the Université de Bordeaux, in Pessac, France. 

is with the Université de Franche-Comté, in Besançon, France. 

is with Cetim Grand Est, in  Illkirch-Graffenstaden, France. Direct correspondence to: 

Articles

LIBS-based imaging has a broad range of applications. Here, we demonstrate those capabilities with examples from paleoclimate research and toxicology studies.

B. Busser

Articles

LIBS-Based Imaging: Recent Advances and Future Directions