1. Introduction
Recently, the use of plasma for producing nanostructures of different materials are gaining growing interest. One of the most powerful methods for producing nanoparticles (NPs) with plasma is the use of laser pulse, namely, pulsed laser ablation in liquid (PLAL) or laser-assisted synthesis in solution (LASiS) [
1,
2]. In this technique, the material ejection is obtained by focusing a laser pulse on a specific target immersed in liquid. The laser-induced plasma (LIP) and cavitation bubble generation are the main processes occurring after the laser-matter interaction. During the plasma and cavitation bubble evolution, the production of NPs, generally with the same composition of the target, occurs. One of the most attractive features of such a tool for NP production is that with a very simple set-up, it is possible to obtain high quality NPs, in terms of purity and size distribution. Moreover, PLAL does not require stabilizers or chemical compounds for their growth and stabilization in solution.
In the last decades, several kinds of nanostructures have been produced with PLAL including metallic NPs [
1], various carbon NPs and nanotubes [
3], alloy NPs, core-shell NPs [
4], etc.
Despite the ease of the instrumental set up, the process is very complex to investigate for several reasons: the plasma produced with laser−matter interaction in water has a maximum dimension comparable to the laser spot size and a lifetime at microsecond scale; the emission of the plasma contains poor information on the ablated species because it is dominated by continuum radiation. A few hundreds of ns after the laser pulse, the cavitation is induced around the plasma and it evolves very fast in the order of hundreds of microseconds changing, in turn, the refraction index. Finally, during the PLAL processes, the observation of the NPs’ evolution is very difficult with visible sources, since the NPs themselves have a smaller size than the visible light wavelength. The previously mentioned issues imply the necessity of using high temporally and spatially resolved techniques to retrieve information on the generation of NPs during the PLAL processes. Optical tools like optical emission spectroscopy (OES), shadowgraph and laser scattering techniques are generally employed for the diagnostics of the processes occurring during PLAL, while most of information is obtained with theoretical calculation [
5,
6]. In this paper, a general overview of these techniques including the advantages and disadvantages is discussed.
2. Materials and Methods
The employed laser was ns pulsed laser (Quanta System, PILS-GIANT, 6 ns pulse duration) operating with the fundamental (1064 nm) and the second harmonic (532 nm), a repetition rate varying between 1 and 10 Hz and with energy varying between 50 mJ and 217 mJ, depending on the experiment performed. The processes of laser ablation in liquid were studied during the production of NPs with high temporal and spatial resolution techniques, as described in further paragraphs. The laser ablation in liquid was performed by focusing the laser beam on a target immersed in deionized water using a stainless steel chamber with a set of optics for laser transmission and the collection of the plasma emission and/or shadowgraph detection. For the optical emission spectroscopy, a spectrometer (Jobin Yvon TRIAX 550) equipped with an intensified charge-coupled device (ICCD) as the detector (Jobin Yvon i3000) was employed. For the shadowgraph technique a collimated light source obtained by a set of lenses and a fast ICCD detection system (New iStar, Andor) coupled with a pulse generator (Stanford inc. DG 535) for the laser pulse-ICCD synchronization were employed. The laser scattering measurements during the PLAL processes were performed with a continuum laser DPSS at 532 nm (Cobolt Samba, 25 mW, CW) by employing the same spectroscopic system used for the shadowgraph experiment. For each presented experiment, the laser focal position onto the target surface and the water were changed after each measurement to always have a flat target and fresh water, respectively. The laser crater on the target had a diameter of 1 ± 0.2 mm.
The targets were titanium (99.995% pure, 1” diameter × 0.125” thick, Kurt J. Lesker) and silver (99.99% pure, 1” diameter × 0.125” thick, Kurt J. Lesker).
More details on the experimental set-ups can be found in Refs. [
7,
8,
9]. Each figure caption reports detailed experimental conditions relative to the shown measurements.
3. Temporal sequence of PLAL processes
The PLAL process can be grouped into four main stages [
10,
11,
12], as represented in
Figure 1:
- (a)
laser ablation and plasma production,
- (b)
plasma expansion and shockwave propagation,
- (c)
cavitation bubble expansion and compression,
- (d)
NPs releasing in solution.
PLAL starts with the interaction of the laser pulse with the target. The laser pulse energy is converted into electron excitation in the material and the generation of charged particles (electrons and ions) occurs by multiphoton ionization. Because of these processes, the target material is ejected in a few tens of fs. In the case of ns irradiation, the laser pulse continues to interact with the ablated material and in particular with the free electrons, releasing further energy and heating the electron and, consequently, inducing the plasma.
The plasma reaches a temperature of several thousands of Kelvin and consequently, it transfers part of its energy to the surrounding liquid. In these conditions, a thin layer of vapor at high temperature is suddenly formed. As a consequence of the high temperature, the vapor pressure is very high and this layer starts to expand at supersonic speed inducing the cavitation bubble. When the bubble reaches its equilibrium with the surrounding liquid, the expansion stops and the liquid mass accumulated at the border of the bubble compresses the vapor inducing the collapse stage. If enough energy has been stored in the bubble, the latter rebounds, inducing further expansion/compression stages. Recently, it has been shown that the compression stage plays an important role in the release of the formed NPs into the solution because the electrostatic pressure of the NP cloud inside the bubble becomes greater than the pressure of the liquid at the border [
8].
4. OES for the Laser Induced Plasma Investigation
As mentioned in the previous paragraph, in the case of ns-laser ablation, the source of material for NP generation is the plasma. For this reason, retrieving information on the plasma system is of great importance to understand the PLAL processes. As a consequence of the high ionization degree and the high temperature, the plasma is characterized by an intense emission of light that can be investigated with the optical emission spectroscopy, OES. The instrumental set up for OES is sketched in
Figure 2, where the integrated plasma emission intensity is collected on the slit of the spectrometer through a lens and the spectrometer is coupled with an ICCD to obtain the temporal evolution acquisition of the plasma emission.
In the case of laser-induced plasma in liquid, the plasma expansion is notably confined, not only for the low compressibility of the liquid phase, but especially because the fast energy transfer from the plasma to the liquid induces the vaporization of the surrounding liquid and the production of a water vapor thin layer [
9]. The latter, that is the starting point of the cavitation bubble [
13], is characterized by high temperature and high pressure and consequently expands in all directions. Therefore, the high pressure of the cavitation at the early stage of bubble formation contains the plasma and confines it in a small volume. The main effect of this strong confinement is that the matter in the plasma phase holds a high number density for the entire plasma lifetime [
9,
12]. This phenomenon can be directly appreciated with the spectral features of the plasma emission. The emission spectrum of the laser-induced plasma is characterized by the absence of discrete emission lines and by a continuum spectrum resembling the Planck-like distribution of photons, as shown in
Figure 3a, where a typical example of temporal evolution of the continuum radiation spectra of LIP is shown. The absence of the discrete emission lines that generally characterize laser-induced plasma in a gas and vacuum environment, is related to the high-density effects. According to Fermi or Debye theories, the cut-off of the levels occurs when excited levels cannot be achieved as a consequence of a sharp lowering of the ionization energy [
13]. This happens because of the high density (10
20–10
22 cm
−3) that constrains the excited atom’s expansion (Fermi)—the distance between particles is smaller than the size of the excited particle—or because of the high electron number density (10
19–10
21 cm
−3) that affects the Coulomb potential of the particle (Debye−Huckel effect) [
14]. For this condition, the radiative recombination is the main contribution to the spectrum as discussed in Ref. [
15]. As a result of this high density effect, the spectral distribution is proportional to the Planck’s distribution [
15].
In the visible range of the spectrum, the Wien’s approximation can be assumed so that [
7]:
In Equation (1), h is the Planck constant, k the Boltzmann constant, c the speed of the light and G(T, Ne) represents a function depending on detection efficiency, on the electron number density (Ne) and plasma temperature (T), but not on the wavelength. The Planck distribution can be plotted in the logarithmic form to determine the temperature of the plasma from the inverse of the slope as shown in
Figure 3b:
Note that the intensity, I
λ, used in the Planck plot is integrated over a small wavelength range (e.g., corresponding to the spectral resolution of the detection system). As an example of the result obtained with temporally resolved OES,
Figure 3c shows the temporal temperature evolution determined from the continuum spectra of
Figure 3a.
Another interesting measurement of the plasma evolution can be performed by the fast imaging technique. The imaging of the total emission of plasma can also be simply detected with an ICCD coupled to an objective for studying the dynamics of the plasma [
11]. Then, plasma emission imaging and OES can be matched for retrieving temporally resolved spatial maps of temperature. This can be obtained as sketched in
Figure 4: the plasma emission image is split and focused on different areas of the ICCD. Along the optical paths, two narrowband filters at two different wavelengths are located in the optical paths of the split image, therefore the simultaneous acquisition of two plasma images at different wavelengths is collected. In this way, it is possible to employ Equation (2) in the form of an intensity ratio for the estimation of the temperature trend across the plasma. Although this approach is not so accurate as the one using a Planck-like distribution obtained with the wavelength range of the whole visible spectrum, if the two filter wavelengths are optically far enough (>100 nm) from each other, the investigation of the temperature of the plasma with a micrometric spatial resolution can be done as shown in
Figure 5. Ref. [
7] describes in detail both the experimental procedures and the used equations for retrieving the plasma temperature maps temporally and spectrally resolved.
As an example,
Figure 5 shows the temperature maps obtained with laser ablation of a titanium target immersed in water. In particular,
Figure 5 reports the results at three delay times from the laser pulse and the corresponding axial temperature profile along the expansion axis of the plasma. It is evident that the temperatures were not homogeneous in the whole plasma volume since the map presented the inner zone with lower values with respect to those at the edge of the plasma (30% hotter). At longer delay times, the temperatures tended to reach homogeneous values all over the plasma volume. On the other hand, from the plasma images, the plasma volume as a function of time could be estimated. This information can be very important if the ablated mass can be measured independently, for example, by measuring the ablated volume of the crater because it allows the calculation of the ablated mass from the density of the target [
17]. When the total ablated mass and the plasma volume are known, the initial number density can be determined. Moreover if the initial number density and the plasma temperature is known, the Saha equation [
18] can be applied and the temporal and spatial evolution of the number density of the species in the plasma can be determined. The details of these calculations are reported in Ref. [
7]. The traditional OES on plasma emission, although it allows good temporal resolution, does not permit spatial resolution and, consequently, the estimated temperature is spatially integrated along the all emitting volume. Since it does not allow the observation of several aspects of the plasma structure, the coupling of plasma OES and emission imaging can introduce a deeper inspection of the plasma characteristics during its expansion dynamics in water. As a matter of fact, the temperature and number density spatial distributions can reveal important information about the formation mechanisms of NPs during the ablation in water, as reported in Ref. [
7].
5. Shadowgraph Technique for the Cavitation Bubble Investigation
The laser-induced cavitation bubble dynamics can be studied by employing a high temporally resolved shadowgraph technique [
9,
19,
20,
21].
In
Figure 6, a typical sketch for the shadowgraph technique is shown.
The bubble is generated by the fast energy exchanges between the plasma and the liquid and its evolution strictly depends on laser pulse characteristics as wavelength, output energy and focusing condition, pulse duration, external pressure [
11,
22]. As an example, at atmospheric external pressure, the ns laser-induced bubble duration is usually in the order of hundreds of microseconds, that is orders of magnitude longer than the plasma duration itself [
23]. The bubble dynamics consists of an early expansion stage, a maximum expansion, reached when the vapor pressure inside the bubble equals the surrounding water pressure, and a further shrinking stage where the bubble collapses. After the first collapse, a bubble rebound can occur with a higher order bubble formation, depending on the initial laser energy [
24]. When higher order bubbles are induced, the expansion/compression stages resemble a damping oscillator, until all energy of the cavitation bubble is dissipated into the surrounding liquid.
The evolution of the radius of a cavitation bubble can be experimentally determined from the shadowgraph images, as reported in
Figure 6b where Ry is employed to determine the bubble radius, Rx is employed to determine the bubble front expansion, z is the shadowgraph detection direction and x is the laser direction. The determination of the radii can be directly employed for the determination of the vapor pressure inside the bubble during the bubble dynamics.
From the experimental point of view,
Figure 7 shows the temporal evolution of a laser-induced bubble on Ag target immersed in water.
Figure 8 shows the experimental bubble radius determined from the shadowgraph images reported in
Figure 7.
As a matter of fact, by determining the bubble radius at each observation time, the estimation of pressure and of the temperature inside the bubble with the hardcore Van der Waals model can be performed [
8,
9]. In Ref. [
25], the Rayleigh–Plesset equation was employed starting from the experimental parameters. In particular, to take into account the liquid compressibility, the Keller–Miksis formulation [
24] was used, starting from the initial experimental radius (in this case, equal to the laser spot size, but the first experimental radius determined with the shadowgraph image can also be employed) and assuming an adiabatic Van der Waals equation of state for the gas in the spherical bubble. With this approach, the radius, and in turn the volume, are related to the vapor pressure and temperature. The theoretical estimation of the bubble radius, pressure and temperature evolution of the vapor inside the bubble can also be done as in Ref. [
26] where the extreme conditions of the pressure reached during the collapse stage were revealed.
Once the estimation of vapor pressure inside the bubble is determined, a comparison with the electrostatic pressure due to the NPs inside the bubble volume can be investigated in order to obtain an evaluation of the number of the NPs ejected in the water during the bubble dynamics. When the bubble begins to expand, the NPs produced in the plasma stage occupy the whole bubble volume following its dynamics. Immediately after the plasma stage, the NPs have absorbed enough charge on their surface to produce an electrostatic pressure whose behavior follows the bubble radius evolution. When the NPs’ electrostatic pressure is higher than the pressure at the water/vapor boundary, the NPs escape from the cavitation bubble and are released in solution. In Ref. [
8], by the comparison between the bubble water pressure and the NPs’ electrostatic pressure and by an estimation of the number of the AgNPs, the releasing of the NPs from the cavitation bubble to the liquid has been modeled.
In
Figure 8, the same procedures described in [
8] have been followed for estimating the percentages of ejected AgNPs during the bubble evolution of
Figure 7. In this case, most of AgNPs ejection occurs immediately after the beginning of the expansion phase and, to a greater extent, between the collapsing stage and the bubble rebound.
6. NP Visualization during the PLAL Processes
In the last years, part of the literature on the PLAL was dedicated to the NP formation mechanism during the different PLAL steps occurring before the NPs releasing in solution. The study of the PLAL processes coupled with the investigation of NP location and the mechanisms of release in solution is a very important goal for the control and the systematic production of NPs. The NPs before the release in solution are trapped inside the cavitation bubble following, of course, the bubble dynamics as discussed in the previous section. The visualization of the NPs inside the cavitation bubble during its evolution by the use of optical diagnostics is extremely difficult just because the bubble surrounds the NPs and light scattering is not so efficient when scatters are smaller than the incident wavelength.
One of the most powerful techniques for obtaining information on the size and mass of particles inside the bubble is the small-angle X-ray scattering (SAXS), usually obtained by synchrotron-radiation [
27]. A simpler way [
10] to check qualitatively the location of the NPs in the bubble is to shoot a second laser pulse focused inside the bubble produced by the first laser pulse in order to produce a secondary plasma from the material ablated by the first pulse. This technique, that is double pulse laser ablation in liquid (DP-LAL), can be efficiently used as a diagnostic tool for understanding the distribution of NPs in the bubble volume. Finally, also the light-scattering-based technique [
9,
28] has been used for the inspection of the materials inside the cavitation bubble during its evolution.
In the following, a brief description of the techniques mentioned above is given.
In Ref. [
29], small-angle X-ray scattering was used during the cavitation bubble evolution to identify the species inside the bubble generated by laser ablation of a silver target in water with temporal and spatial resolution. Since the X-rays can penetrate the bubble interior, X-ray scattering can analyze the nanometer-sized particle distributions.
Figure 9a shows a sketch of the employed experimental set-up, and
Figure 9b, the main results obtained with this technique. The authors of Ref. [
29] show in detail how, in the bubble volume, during the bubble evolution, two different groups of particle size distributions can be identified. One contains particles around 8–10 nm size and the other one contains bigger agglomerates around 40–60 nm size. Further research has been done by the same group, improving the temporal resolution of the acquisition [
30] and the quantification as well as a survey of the crystallinity of NPs during the release in solution [
31].
In
Figure 10a, a sketch of the experimental set-up for the DP-LAL technique is shown [
32,
33]. Using the double pulse technique as a diagnostic of the materials inside the cavitation bubble, the first laser pulse is employed for the PLAL, while the second laser pulse is fired collinearly to the first one and focused on the target at different interpulse delays between the two laser shots. By simultaneously acquiring the shadowgraph of the bubble and the images of the plasma generated by the second pulse, it is possible to visualize the second plasma at a different lifetime of the bubble. Depending on the position of the NPs inside the bubble, the second plasma is generated. In part on the target, in part on the materials that occupy the bubble.
Figure 10b clearly shows this phenomenon, described in detail in Ref. [
10], where it is possible to discriminate between both the NPs inside the bubble, located immediately below the expansion bubble front, and the materials ejected in solution after the bubble collapse.
Figure 11a shows the sketch of the experimental set-up for the laser scattering technique. The continuum laser employed for scattering is perpendicular to the bubble axis propagation, and the images of the scattering are acquired during the cavitation bubble evolution. Of course, the scattering can easily come from the bubble surface instead of penetrating the bubble itself. To overcome this drawback, usually more than one technique can be used, to retrieve some useful information. Spectroscopy, shadowgraph and laser scattering measurements with high temporal resolution can be applied simultaneously to analyze respectively the laser-induced plasma, the bubble and the delivered particles.
As a matter of fact, as described in [
9], employing the three high temporally resolved techniques, it was seen that the scattering signal suggested that the produced particles, after the bubble collapse, were lying close to the target surface.
Figure 11b shows an example of laser scattering coupled with the shadowgraph technique, where the images obtained during the laser scattering were acquired with the same experimental conditions of the shadowgraph images reported in the upper part of
Figure 11b. Finally, in Ref. [
34], the laser scattering was coupled with simultaneous spectroscopy of the materials inside the bubble.
To probe the materials inside the bubble employing laser scattering without the previously described drawbacks, Ref. [
35] proposed a variation of the classical experimental set-up, where the probe laser was introduced into the ablation bubble through a flat and transparent window created by the pump laser, after the ablation of a thin Ag film on a glass substrate.
Figure 12a shows a sketch of the experimental set-up where it is clearly visible that the laser probe enters the bubble from the bottom side of the bubble expansion direction, therefore it can persist inside the bubble during its evolution stages (expansion and shrinking). This NP visualization method is based on the different polarization dependences of the NP clouds from the reflection inside the bubble, that allows the distinction between the scattered signal of the NPs and the background reflection.
Figure 12b shows the results of this experiment where it is visible that the dim spots, generated on the materials inside the bubble, are all inside the bubble volume and along the probe laser path. More details can be found in Ref. [
35].
7. Conclusions
Several optical techniques for the diagnostics of the PLAL process have been presented and discussed. A summary of the advantages and shortcomings is reported in
Table 1.
The use of OES allows the determination of the plasma parameters that may be extremely useful for investigating the source of matter of the all PLAL processes. Moreover, OES can be coupled with imaging in order to obtain a temporally and spatially resolved map of the plasma parameters.
On the other hand, shadowgraph imaging permits the investigation of the cavitation bubble dynamics. The cavitation stage has focused the attention of many researchers because it is related with the production of the final NP form and with the releasing of the nanoparticle in solution. Shadowgraph imaging enables the determination of the bubble shape and radius. These are crucial pieces of information for investigating the cavitation thermodynamics and the rate of the expansion and compression, as the radius and in turn the volume can be related with the thermodynamic parameters with different real gas-state equations.
Finally, the visualization of the produced nanostructure during the PLAL process can be investigated with scattering techniques or double pulse techniques. Among the scattering techniques, SAXS provides the best information, as light scattering has poor efficacy with scatters smaller than the employed wavelength. The double pulse technique can be a good option for a qualitative investigation of the position of the material in the bubble as well as for visualizing the released NPs.