Direct visualization of the complete evolution of femtosecond laser-induced surface structural dynamics of metals.

Despite extensive studies of femtosecond laser-material interactions, even the simplest morphological responses following femtosecond pulse irradiation have not been fully resolved. Past studies have revealed only partial dynamics. Here we develop a zero-background and high-contrast scattered-light-based optical imaging technique through which we capture, for the first time, the complete temporal and spatial evolution of the femtosecond laser-induced morphological surface structural dynamics of metals from start to finish, that is, from the initial transient surface fluctuations, through melting and ablation, to the end of resolidification. We find that transient surface structures first appear at a delay time on the order of 100 ps, which is attributed to ablation driven by pressure relaxation in the surface layer. The formation dynamics of the surface structures at different length scales are individually resolved, and the sequence of their appearance changes with laser fluence is found. Cooling and complete resolidification, observed here for the first time, are shown to occur more slowly than previously predicted by two orders of magnitude. We examine and identify the mechanisms driving each of these dynamic steps. The visualization and control of morphological surface structural dynamics not only are of fundamental importance for understanding femtosecond laser-induced material responses but also pave the way for the design of new material functionalities through surface structuring.

Introduction

Femtosecond (fs) laser-material interactions have been extensively studied in the past[1, 2, 3, 4, 5]. However, a clear understanding of even the simplest morphological responses is still lacking; direct visualization of the complete ultrafast structural dynamics in metals has never been achieved. Interestingly, fs laser-induced morphological changes often drastically alter a material’s physical properties. For example, surface nano- and microstructuring induced by fs laser pulses has led to the creation of various materials with novel effects, such as black silicon, black and colored metals, and superhydrophilic, superhydrophobic, and multifunctional surfaces[2, 3]. These highly functionalized surfaces have found a wide range of applications throughout the physical and biological sciences[1, 2, 3, 4, 5].

Ultrafast morphological and crystalline structural dynamics have been studied using various techniques, including X-ray diffraction[6, 7], electron diffraction[8], diffractive XUV imaging[9], and optical imaging[10, 11, 12, 13, 14, 15, 16]. Ultrafast X-ray and electron diffraction studies are more suitable for investigating long-range disorder phenomena in crystalline materials, such as melting[8]. In a more recent work on X-ray diffraction combined with photon-correlation spectroscopy, real-space images were extracted through phase retrieval at delay times of up to 15 ps, but such retrieval becomes increasingly more difficult and less reliable at longer time delays[17]. Diffractive imaging using time-resolved coherent XUV scattering in a transmission geometry can also be applied for indirect imaging, in which a real image is reconstructed from a diffraction pattern[9]. Optical techniques are more suitable for real-space imaging. However, optical imaging is usually achieved by collecting specular reflections[10, 11, 12, 13, 14, 15, 16] and is more suitable for resolving relatively large changes in reflection, such as Newton-ring formation[11, 12, 13, 14, 15]. The strong background of specular reflections makes it difficult to discern the small changes in reflection that are often associated with nano/microstructure formation. Here we develop an ultrafast optical imaging technique to study the structural dynamics of metals by utilizing scattered light instead of specular reflections. Scattered-light imaging has been used to study biological and organic materials because scattering is a prevalent form of light interaction with inhomogeneous media, but scattered-light imaging has not been explored for studying the structural dynamics of solid surfaces[18]. The scattered-light imaging provides a near-zero background and high contrast, allowing us to obtain high-quality real-space images of the structural dynamics of metals and resolve, for the first time, the complete temporal and spatial evolution of fs laser-induced morphological surface structure formation from start to finish in metals.

Materials and methods

Ultrafast pump-probe imaging

Our ultrafast pump-probe imaging setup is depicted in Figure 1a. We use an amplified Ti:sapphire laser system that generates 65-fs pulses with a pulse energy of 1 mJ, a repetition rate of 1 kHz, and a central wavelength of 800 nm. A fast electromechanical shutter is used to select a single pulse. A pump beam induces the formation of surface structures on Zn. Surface nano/microstructures are produced by a single laser pulse at a fluence in the range of 0.1–1.0 J cm−2, under which extensive nano/microstructure formation is observed[19]. The incidence angle of the pump beam is 36°. The Zn sample is first mechanically polished using 0.05-μm Al2O3 powder, washed in methanol, and finally cleaned using Ar plasma. The probe beam is used to provide stroboscopic illumination of the surface structures induced by the pump beam. An optical delay line is used to provide various time delays, ranging from 0 to 408 ns, between the pump and probe pulses. A delay of 408 ns is achieved by allowing the beam to propagate over 120 m through multiple reflections. A zero delay is achieved by using the autocorrelation technique using a BBO crystal. After passing the delay line, the probe pulse is directed onto the BBO crystal to generate second harmonics at 400 nm. A blue band-pass filter is used to block the fundamental light after the BBO crystal. The 400-nm probe pulse is incident on the sample at an angle of 18°. A 50-mm (numerical aperture=0.25) long-working-distance objective is used to image the sample surface onto the sensor of a CCD camera (1600 × 1200 pixels). The CCD camera is positioned normal to the sample surface. A 400-nm narrow band-pass filter is placed in front of the imaging objective to block the pump light. The imaging optics we use are capable of resolving both nanoclusters (microscale size) and microstructures. At a fixed laser fluence, we acquire a sequence of three surface images. The first snapshot is taken 10 s before the pump laser pulse irradiation for the identification of initial structures, if any, that may be present on the surface. The second snapshot of the surface is taken at a preset time delay. The third snapshot is taken 10 s after the pump laser pulse irradiation to record the permanent surface structures produced by the laser pulse. By recording such sequences of surface images at different time delays and laser fluences, we can track the formation dynamics of the transient surface structures at different delays and identify the effect of the laser fluence on their evolution, including transient surface structural fluctuations and solidification. The scattered-light images may contain speckle, but the speckle is not pronounced for femtosecond pulse imaging because femtosecond pulses have a very short coherence length.

Figure 1

(a) Ultrafast pump-probe imaging setup. (b–e) SEM images of surface nano/microstructures of Zn in the center of the irradiated spot at various laser fluences, as indicated in the figure.

Reflectance measurements

The optical reflection from a surface will change when the surface is irradiated with a high-intensity laser beam, and the surface reflectance at a given fluence is crucial for our study. Therefore, we also measure the surface reflectance at different fluences by using a hemiellipsoidal reflector to collect both the specular and diffuse light reflection[20]. The studied sample is placed at the internal focal point of the reflector. The reflected light is collected at the external focal point of the reflector, and its energy ER is measured using a joulemeter. The sample is tilted at the same incidence angle as in our imaging setup with respect to the pump beam (36°). A second joulemeter is used to measure the incident pulse energy, EI. The reflectance R=ER/EI is obtained from the measured ER and EI. To calibrate our reflectance measurements, we compare the reflectance under low-intensity laser light (below the damage threshold) with that measured at 800 nm using a Perkin–Elmer Lambda 900 spectrophotometer with an integrating sphere.

Characterization of surface structures

A scanning electron microscope (SEM) and a three-dimensional laser scanning microscope are used to examine the laser-induced surface nano/microstructures after the resolidification of the surface.

Results and discussion

The single-pulse damage threshold of our zinc sample was measured to be 22 mJ cm−2. This threshold was determined based on the first appearance of surface modification observed under an optical microscope. Figure 1b–1e are SEM images of some typical surface nano/microstructures produced in the center of the ablated spot at various laser fluences, where one can see that nanostructures dominate at lower fluences, whereas microstructures dominate at higher fluences. Figure 2 shows the transient structures of Zn at different delays for a laser fluence of 0.1 J cm−2, which is approximately four times larger than the damage threshold. For each delay, the final image of the same spot is also shown for comparison. It is seen that surface structures begin to appear at a time delay of 400 ps (Figure 2c), but a comparison with the structures observed a long time after the laser pulse (t=∞) (Figure 2d) reveals that the two images show nearly no resemblance, indicating that the material in the irradiated spot is not yet solidified, and therefore, the surface structures observed at this delay are transient. A comparison of Figure 2e and 2f reveals some similarity at the periphery of the irradiated spot, which indicates that the transient structures begin to solidify at 527 ps and that interestingly, solidification first occurs at the periphery. With increasing delay time, the solidified area increases. Figure 2g and 2h shows that ~90% of the surface structures are frozen at t=50 ns. Finally, Figure 2i and 2j shows that the structures become identical to their final configuration at 124 ns, indicating that the surface structures are completely solidified at this time. Figure 1b shows that for a fluence of 0.1 J cm−2, the surface structures are dominated by extensive nanostructures with dimensions between 30 and 500 nm, and some microstructures. Note that the optical imaging system used here, operating with 400 nm light, resolves nanoclusters at the micron scale but not individual nanostructures. To fully characterize the surface structural evolution dynamics, three characteristic times are important: (1) the time at which transient surface structures first appear, t1; (2) the time at which solidification begins, t2; and (3) the time at which the surface structures become completely solidified, t3. These three characteristic times as found in our study for different pump fluences are summarized in Table 1. Previously, fs laser-excited surface structure formation has been studied for semiconductors[12, 13]; in these studies, the time required for complete solidification in silicon was found to be ~75 ns, much shorter than what we observe here for Zn.

Figure 2

(a–j) A comparison of transient surface structures observed at various delay times and final solidified structures following a pump pulse at a fluence of 0.10 J cm−2.

Table 1

Measurements and estimates regarding surface structural dynamics

F (J cm−2)	t1 (ps)	t2 (ns)	t3 (ns)	R	Tm (K)	tcool (ns)	Ablation mechanism
0.022	—	—	—	0.57	890	—	Damage threshold
0.10	400	0.527	123.8	0.57	2970	32	Spallation
0.14	300	0.86	322	0.57	4080	41
0.19	200	1.3	—	0.57	5430	50
0.26	100	2.8	—	0.56	7320	62
0.34	50	5	—	0.57	9480	74	Phase explosion
0.54	—	9.3	—	0.57	14880	101
0.78	—	50	—	0.55	22330	130
1.0	—	124	—	0.56	28580	156
1.1	—	—	—	0.5	—	—	Solid-density plasma regime
1.8	—	—	—	0.33	—	—

First, we study the mechanism of the transient surface structure formation and its associated start time, t1. Previous studies have shown that the characteristic time for surface melting following the ultrafast laser irradiation of Zn is ~3 ps at laser fluences of 0.1–0.15 J cm−2 (Ref. 21). However, our data show that the onset of transient surface structure formation occurs at ~400 and ~300 ps at laser fluences of 0.1 and 0.14 J cm−2, respectively. This time difference suggests that the irradiated surface does not undergo a morphological change immediately after melting; instead, there is a time delay between the melting process and the hydrodynamic motion of the melted material that results in the transient morphological fluctuations observed on the surface. In general, transient structural fluctuations can be induced by the flow of the melted material due to a temperature gradient (the Marangoni effect) and the flow due to laser-induced stresses in the melted surface layer. To identify whether the Marangoni effect plays the key role here, we estimate the characteristic timescale tM of the Marangoni effect as follows[22]:

where η is the dynamic viscosity of the melt, L is the typical radial dimension, h is the average melt depth, Tm is the melting point of the metal, and |γT|=|dγ/dT| is the absolute value of the temperature coefficient of the surface tension γ. Using values of η=2.5 × 10−3 Pa·s (Ref. 23), |γT|=0.25 × 10−3 N m−1 K−1 (Ref. 24), Tm=693 K, and L=200 nm (determined from Figure 1b) and assuming h=30 nm, we obtain tM≈20 ns, which is two orders of magnitude larger than the observed values of t1 (which are on the ps timescale). Therefore, the Marangoni effect does not appear to be the primary cause of the onset of the transient surface structures. However, this does not exclude the Marangoni effect as a factor affecting the formation of surface structures at longer timescales.

Following high-intensity ultrashort laser pulse heating, a highly pressurized melted surface layer can form as a result of isochoric heating, as has been theoretically predicted[25, 26]. Pressure relaxation in this melted layer will induce hydrodynamic motion, resulting in morphological fluctuations on the surface. The timescale of these morphological fluctuations is governed by the speed of sound in the melted material[26] and can be estimated as tP≈L/vsl, where vsl is the speed of sound in liquid Zn. The speed of sound in solid Zn is 3.8 × 103 m s−1 (Ref. 27). Typically, the speed of sound in liquid metals is 1.5 times smaller than that in solids. Therefore, assuming vsl=2.3 × 103 m s−1 and again using the typical radial dimension L=200 nm found from Figure 1b, we obtain tP≈90 ps, which is on the same timescale as our measured t1 values. In general, melting occurs at laser fluences near the damage threshold, and ablation occurs at higher fluences. The laser fluence we use in our study is more than four times greater than the damage threshold. Therefore, we believe that the start time of transient surface structure formation, t1, as observed in our experiment is associated with ablation driven by pressure relaxation in the surface layer. Figure 3 presents t1 and t2 as functions of the laser fluence. We can see that t1 decreases with laser fluence, which indicates an earlier onset of ablation due to higher pressure in the surface layer.

Figure 3

Time at which transient surface structures appear, t1, and start time of solidification, t2, as functions of laser fluence.

Figure 4 shows a set of time-resolved surface images acquired at various delays following pump pulse irradiation at a higher fluence of 1.0 J cm−2. Somewhat counter- intuitively, we can see that transient surface structures first appear at the edge of the irradiated spot and form a ring, in contrast to what occurs under ablation at lower laser fluences, for which the transient structures first appear in the center (Figure 2c). Over time, the transient surface structures move toward the center. At 5–9.3 ns, the center of the irradiated spot is covered with structures. The solidification of the transient surface structures begins at ~124 ns, as indicated by the resemblance to the final image. At the longest delay investigated in our experiments (408 ns), solidification is not fully complete, as a difference still exists between Figure 4k and 4l. A comparison of Figures 2 and 4 shows that the spatial sequencing of the transient structures before solidification varies with the laser fluence. To better understand this dependence, we consider surface images acquired at a fixed delay time of 527 ps but at different laser fluences, as shown in Figure 5. We can see that a transition from first emergence of the transient structures in the center to first emergence at the edge occurs at 0.54 J cm−2. The SEM images presented in Figure 1 show that at this laser fluence, the surface structures in the center also change from being mostly nanostructures to being mostly microstructures. At this and higher laser fluences, nanostructures can still be produced at the edge of the ablated spot because of the Gaussian profile of the laser beam. Therefore, the spatial evolution of the ring surface pattern observed between 300 ps and 2.8 ns in Figure 4 is attributed to the fact that nanostructures first form at the periphery and the center region is later populated with microstructures; as shown earlier with regard to the timescale of morphological fluctuations (tP≈L/vsl), the time required to form nanostructures (smaller L) is shorter than that required to form microstructures (larger L). Another possible explanation for the central dark spot is the formation of highly light-absorptive matter in the central part of the interaction area, possibly due to the formation of an ablation plasma that masks the surface structures. Over time, such an ablation plasma would become more transparent, and the surface structures in the central part of the ablated area would become more visible. We note that the absorption of laser light by an ablation plasma has been previously observed in double-pulse ablation studies[28, 29].

Figure 4

(a–l) CCD images of a Zn surface at various delay times following a pump pulse at a fluence of 1.0 J cm−2.

Figure 5

(a–h) CCD images of Zn surfaces at a fixed delay time of 527 ps following pump pulses at various laser fluence values.

Next, we study the mechanisms of the evolution of the transient surface structures from their first appearance (t1) up through the onset of solidification (t2). Several possible ablation mechanisms exist[25, 26, 30, 31, 32], that is, evaporation[31], thermomechanical fragmentation/spallation[25, 26] and phase explosion[25, 32]. To identify the ablation mechanisms that are responsible to our observations, we estimate the maximum lattice temperature Tm (Refs. 21, 33):

where T0 is room temperature, c is the specific heat of the lattice, α is the light attenuation coefficient, g is the electron–phonon coupling factor, k is the thermal conductivity, and FA=(1−R)F is the absorbed laser fluence, where R is the reflectance of Zn at the wavelength of the laser. The reflectance of a metal surface depends on its condition after polishing, the angle of the incident light and the laser fluence[20]. The key to obtaining a reliable estimate of the surface temperature is to accurately determine the value of the sample reflectance under our experimental conditions. To do so, we perform a rigorous surface reflectance measurement using a hemiellipsoidal collector to collect both the specular and diffuse reflected light, as first demonstrated in Ref. 20. The measured reflectance in a fluence range of 0.01–1.0 J cm−2 is R=0.55–0.57, as shown in Table 1. Table 1 also shows the surface temperatures estimated using T0=300 K, ci=2.78 × 106 J m−3 K−1 (Ref. 34), α=7.5 × 107 m−1 (Ref. 34), g=6 × 1016 W m−3 K−1 (Ref. 34), and k=116 W m−1 K−1. As we can see, at the damage threshold laser fluence of 0.022 J cm−2, Tm is estimated to be 890 K. This value is only slightly higher than the melting point of Zn (693 K), indicating that Equation (2) provides a very good estimate of the surface temperature. At F=0.1 J cm−2, the surface temperature is estimated to be 2970 K, exceeding the boiling point of zinc (1180 K), and this laser fluence is clearly in the ablation regime. At laser fluences of 0.14 J cm−2 or above, the estimated surface temperatures are above the thermodynamic critical point Tc (3190 K), suggesting that the laser-excited material becomes a supercritical fluid (with no distinction between the liquid and vapor phases). The surface temperature for ablation at F=0.1 J cm−2 can lead to either vaporization[31] or thermomechanical fragmentation/spallation[25, 26]. However, it is believed that the predominant contribution to ablation is from thermomechanical fragmentation/spallation, whereas the role of vaporization is not essential[35]. Therefore, we believe that our results at a fluence of 0.1 J cm−2 are mostly related to ablation through the spallation mechanism. In the spallation regime, an ultrashort laser pulse induces fast isochoric heating of the surface layer, which causes a buildup of high stress, cavitation, nanofoam formation and thermomechanical spallation (ablation) of a liquid layer on the surface due to the relaxation of tensile stresses developed in the surface layer[25, 26, 36, 37, 38]. For a laser fluence of higher than 0.14 J cm−2, the surface temperature is estimated to be significantly higher than 0.9 Tc, which is the characteristic threshold temperature for phase explosion[32]. Phase explosion occurs at laser fluences below those for substantial ionization and plasma formation[38]. It is known that plasma formation causes a significant decrease in reflectance[20, 39]. Our reflectance measurements show that a significant decrease in the reflectance of the Zn surface occurs at laser fluences above 1.1 J cm−2; see Table 1. This indicates that in the laser fluence range of 0.14–1.1 J cm−2, the ablation is primarily caused by phase explosion. Therefore, the transient structural dynamics between t1 and t2 originate from either the spallation mechanism or the phase explosion mechanism, depending on the laser fluence. As shown in Figure 3, both t1 and t2 depend on the laser fluence. The start time of solidification, t2, increases with laser fluence, indicating that a longer time is required for solidification to occur at a higher laser fluence.

Finally, we study the mechanism of complete solidification and its associated time, t3. After ablation, rapid cooling due to high electronic heat conduction causes solidification of the melted surface layer, giving rise to permanent surface structures. The cooling and resolidification of the melted surface layer in the spallation regime has previously been modeled using the atomistic molecular dynamics method combined with a continuum-level two-temperature model[36, 37]. In these works, the time of complete resolidification was predicted to be 2.5 ns for Al (Ref. 36) and 1.4–1.7 ns for Ag (Ref. 37). These timescales are approximately two orders of magnitude smaller than our observed complete solidification time, t3, for Zn (Table 1). It is generally believed that the cooling of an ultrafast laser-heated metal surface is dominated by thermal conduction[21, 33]. Under these conditions, the cooling time tcool is estimated to be[21]

where χ is the thermal diffusivity and IA is the absorbed laser intensity. The estimated values of tcool for Zn are shown in Table 1, where we see that tcool is significantly shorter than the observed t3. Thus, both the model predictions[36, 37] and the estimation based on Equation (3) predict significantly faster cooling than that observed here in our experiments, indicating that the current understanding of the rapid cooling process is incomplete. To explain this discrepancy, we speculate that the slower cooling process is due to an enhanced thermal coupling that has previously been observed and modeled[40, 41], where it was found that a significant amount of thermal energy remains in a metal following fs laser ablation because of the transfer of thermal energy from the ablation plume back to the sample on a timescale that is much longer than the laser pulse. Our speculation is supported by previous studies of residual thermal coupling in the ablation of Zn (Ref. 40), which indicates that the onset of the enhanced thermal coupling occurs at F=0.45 J cm−2 and is correlated with an abrupt increase in the start time of solidification, t2, above 0.54 J cm−2, as seen in Figure 3. Furthermore, it has been shown that the heating of a sample through the transfer of energy from the ablation plasma causes ultradeep drilling of the material[42]. In addition to the enhanced thermal coupling, surface cooling can also be slowed by the latent heat of resolidification and by exothermic chemical reactions induced during laser ablation[43]. In one example of the large contribution of chemical sources of energy, this contribution has been shown to reach up to 30% of the incident laser energy during the nanosecond laser ablation of zinc in air[43].

Conclusions

In conclusion, for the first time, we have captured and visualized the complete temporal and spatial evolution of the fs laser-induced surface structural dynamics of metals by utilizing a time-resolved optical imaging technique. We find that transient surface structures first appear at a delay time on the order of 100 ps, and this is attributed to ablation driven by pressure relaxation in the surface layer. At lower laser fluences that favor nanostructure formation, the transient surface structures first appear in the central portion of the ablated spot. At higher laser fluences that favor microstructure formation, nanostructures first emerge at the periphery, and the center region is later populated with microstructures. For zinc, the predicted cooling time is two orders of magnitude shorter than the time found in our observations, and we speculate that this slower cooling process is due to an enhanced thermal coupling phenomenon that has been previously observed. The visualization and control of surface structural dynamics not only are of fundamental importance for understanding the fs laser-induced responses of materials but also pave the way for the design of new material functionalities through surface structuring.