1. Introduction
Intense, ultrashort laser–solid interactions at intensities ranging from 10
13 to 10
16 W/cm
2 are of paramount importance for laser nanoprocessing to achieve functional surfaces. Among these interactions, one of the most notable examples is the generation of Laser-Induced Periodic Surface Structures (LIPSS) [
1], which find applications, e.g., in antibacterial coatings, optical devices, chemical sensing, and tribology [
2,
3,
4,
5]. Understanding the dynamics of both surface and sub-surface phenomena at nanometer (nm) or even atomic scales, within the requisite temporal frame, is pivotal for comprehending the underlying physics responsible for creating desired surface structures in a controlled manner. Despite notable advancements in theoretical frameworks and models in recent years, aimed at elucidating intricate mechanisms involved in the self-organization of nanostructures under ultrashort laser irradiation within relevant temporal scales [
6], direct experimental visualization of this process remains scarce. Thus far, the majority of experimental findings regarding surface manufacturing have relied on post-mortem analyses, which lack temporal dynamics information. Time-resolved experiments employing optical lasers, such as optical reflectometry, interferometry, and spectroscopy, suffer from limited spatial resolution and a lack of bulk sensitivity. Laser processing involves a complex chain of various physical processes occurring across different temporal and spatial scales. Initially, ultrashort laser interactions with metals excite electrons within the surface skin layer, typically spanning tens of nm. These excited electrons then propagate into the bulk at Fermi velocity before rapid thermalization occurs within sub-picoseconds (ps) timescales via collisions. Depending on the excitation strength and material properties, this process triggers electron-lattice/ion thermalization, coherent phonon oscillation, and lattice heating in ps timescales. This can lead to subsequent thermal or non-thermal melting, followed by ablation through spallation in sub-nanosecond(ns) timescales. These phenomena collectively contribute to later crystallization or amorphization, leading to specific surface nanostructures. Therefore, understanding above-mentioned early time phenomena in (sub-)ps dynamics is crucial. To address this need, here, we demonstrate experimental capabilities of visualizing physical processes at both surface and subsurface levels with ps and nm resolutions using an X-ray Free Electron Laser (XFEL) operating in a grazing-incidence geometry. Grazing-incidence X-ray small-angle scattering (GISAXS) is a well-established technique for probing lateral structures and correlations along the surface normal. It has been extensively utilized over the past decades at synchrotron X-ray facilities to characterize material structures across length scales ranging from sub-nm to micrometers (
) [
7,
8], albeit limited to the millisecond timescale due to photon accumulation requirements.
Recently, we demonstrated the applicability of GISAXS to XFEL to track nanometric multilayer (ML) dynamics at unprecedented ps scales, which is six orders of magnitude faster than previously achievable [
9]. This breakthrough is enabled by the XFEL’s ability to deliver an immense number of photons, comparable to those typically accumulated over a second using third-generation synchrotron X-rays, all within a single pulse lasting only ∼10 femtoseconds (fs). Despite inherent time smearing due to grazing-incidence, providing different parts of the X-ray arriving at the surface at different times, achieving approximately a ps time resolution is feasible.
Utilizing the in-plane scattering signal, characterized by the wavevector transfer along the sample depth (
), which correlates closely with specular reflectivity [
10], our experiments revealed the compressed, heated, intermixed, and ablated ML dynamics after the laser irradiation. On the other hand, out-of-plane diffuse scattering along the
direction provides insights into ultrafast changes in the lateral distribution of surface roughnesses. This can be described with the help of the height–height correlation function
where
R is the spatial separation of two points,
is the root mean square (RMS) roughness,
is the lateral correlation length and
H the Hurst parameter [
11]. The correlation function along the sample depth can be written as
Here,
is the height–height correlation function,
z denotes the depth inside the sample and
the cross-correlation length. The indices
j and
k denote different interfaces of the ML sample [
12].
Additionally, by incorporating a low-Z capping layer on top of a high-Z ML, GISAXS facilitates measurements of energy transport several hundred nm within the material. Upon ultrashort laser irradiation, strongly confined temperature gradients induce compression (shock) waves, where the progression of temperature gradients depends on the mean-free-path of excited electrons which is a function of both temperature and density. Our results show that increasing laser excitation strength transitions the dominant energy transfer mode from compression wave dominance to the electron–thermal wave dominant regime.
2. Experimental Methods and Results
The experiment was performed at the BL2 EH6 station at the SACLA XFEL facility in Japan [
13]. A metallic ML sample, consisting of 5 repetitions of tantalum (Ta,
) and copper nitride (Cu
3N,
) capped with a
thick Aluminum top layer was irradiated by an optical laser with a central wavelength of
, intensities of 8 × 10
14 and 8 × 10
15 W/cm
2, pulse duration of
and a focus size of ∼500
m (
Figure 1a). After a variable delay time, the surface state of the sample was probed via surface-sensitive X-ray scattering employing ultrafast XFEL pulses of X-ray photons of energy
. The X-ray beam with a pulse energy of ∼0.1 mJ/pulse was incident under a grazing incidence angle of
and the diffusely scattered intensity was recorded by a multi-port CCD (MPCCD) area detector [
14] with the intense specular peak (incident angle equals exit angle,
) being blocked by a beam stop. More details about the experimental setup can be found in the
Appendix A.
Figure 1a displays a typical single-pulse GISAXS pattern. The shallow angle-of-incidence, in combination with the beam size of
leads to a large footprint (∼300
m) of the beam on the sample surface, limiting the temporal resolution of the detected X-ray signal to ∼1 ps.
Figure 1b illustrates the in-plane scattering signal (
) as a function of
for various delay times after laser irradiation. Because of the peculiar geometry of GISAXS, the in-plane scattering signal contains momentum transfers with components in both the normal (
) and surface-parallel (
) directions. As a result, the scattering signal reflects vertical density correlations along the z-direction (
) of the ML structure, roughness correlations along the surface plane (
) direction, and cross-correlations between different interfaces. For clarity, only the
values are plotted here, as
. The in-plane signal at
is closely associated with the specular reflectivity curve [
10], which we characterized ex situ. The spacing of the intense Bragg-like peaks at
,
and
corresponds to the typical length scale, i.e.,
thickness of each Ta/Cu
3N double layer. The Kiessig fringes [
15], represented as smaller peaks between the Bragg peaks, are a fingerprint of the number of double-layer repeats in the sample. Upon laser irradiation, substantial changes in the in-plane scattering signal are observed during the first
. The number and intensity of Kiessig fringes remain relatively stable during the initial ∼9 ps after laser irradiation, thereafter exhibiting a gradual decrease, concomitant with a simultaneous broadening and reduction of the intense Bragg-like peaks. After
, the Kiessig fringes are nearly gone, leaving behind only a broad residual of the Bragg-like peaks.
The most notable structural changes inside the ML sample are manifested by the change of intensity and position of the intense peaks at
,
and
.
Figure 2a illustrates the temporal evolution of the peak at
following laser irradiation with an intensity of 8 × 10
14 W/cm
2. While only the peak at
is displayed here, qualitatively similar behaviour is observed for the peaks at
and
. The solid lines represent refinements using a Gaussian function. The peak amplitude as a function of time delay, extracted from the Gaussian refinement, is depicted in
Figure 2b), with blue dots, which shows an ultrafast increase within the first ps. Subsequently, the intensity remains constant, followed by a gradual decrease at around ∼10 ps associated with a loss of structural correlation between the layers. This decay is modelled using an error function
resulting in values for
and a width of
. Here,
A represents a stretching factor,
denotes the time at which the error function decreased to half of its initial value,
is the duration of the decrease, and
B stands for a constant offset. In comparison, the same analysis for the higher-laser-intensity case of 8 × 10
15 W/cm
2, is depicted with red circular dots. Similar to the lower-intensity case, the ultrafast increase in scattering intensity within the first ∼1 ps is evident. However, for the time delays ≥1 ps the amplitude quickly begins to decline again. The refined error function yields a time constant of
and a width of
. On top of the change in the Bragg intensity, we observe a steady increase in the Q-position with time delay for both laser intensities, as summarized in
Figure 2c. This shift towards larger Q-values implies compression of the double-layer structure as the typical length scale decreases. Assuming that
we result in an exponent
for the lower intensity and
for the higher intensity, respectively.
The dynamic diffraction effects at
yield insight into the dynamics of the top Al layer. In particular, the scattering signals around exit angles equal to the critical angle for the total external reflection—
=
for Al—are particularly surface-sensitive, originating from an evanescent X-ray wave travelling parallel to the surface. This so-called Yoneda peak [
16], primarily originated from interference within the topmost surface layer, serves as a sensitive marker of its surface structure. A close-up of this Yoneda peak is summarized in
Figure 3a. The peak at
exhibits a progressive decrease after laser irradiation, which persists for at least at
before almost disappearing at
, indicative of the ongoing presence of the solid density Al cover layer. The surface ablation speed can be approximated by the speed of sound,
, where
and
are electron and ion temperature, respectively,
is the mean ionization,
and
are adiabatic index of electrons and ions, respectively and
a.u. is the ion mass for Al [
17]. In our temperature range (
= 5–10 eV,
= 1–5 eV) the ablation of the
layer takes 10–20 ps, which aligns with the timescale for the disappearance of scattering from Al observed in
Figure 3a. The integrated intensity of the Yoneda peak (
= 0.65–0.79 nm
−1) is depicted in
Figure 3b for both low- (blue) and high- (red) intensity cases. In both cases, a rapid initial decay in intensity is observed within the first ps, followed by a period of quasi-constant signal intensity extended up to ∼12 ps and ∼5 ps before subsequent decrease, respectively. This temporal behaviour aligns qualitatively with that observed at the Bragg-like peak at higher
shown in
Figure 2b. In
Figure 3a, two additional Yoneda peaks are also visible at higher
, corresponding to the Yoneda peaks for
(
) and Ta (
), respectively. Given that these materials are embedded within a
thick Al layer, the time-dependency of their intensity is not as pronounced as that of the Al Yoneda peak, except during the period between 12 and
when all peaks suddenly disappear. As we will discuss below, aided by plasma simulations, this phenomenon can be attributed to the arrival of a compression wave to the ML, which initiates significant modulation of the ML structure. An additional noteworthy feature in the figure is the presence of high-frequency fringes superimposed on the Ta and Cu
3N Yoneda regions. These fringes arise from the interference between waves scattered from the Al surface and those from the Al-ML interface. Following laser excitation and the subsequent expansion of the top Al surface, these fringes rapidly vanish as a result of blurring of the distinct interface between vacuum and Al.
Both Bragg-like peaks and the Yoneda peak indicate the presence of three distinct time regimes. Immediately after laser irradiation ∼1 ps, we observe an ultrafast increase in the intensity of the Bragg-like peaks without any change in their -positions, accompanied by the decay of the Yoneda peak. As we will discuss further below, this phenomenon is attributed to a modification in the surface structural properties. In the second time frame, spanning from 1 to ∼10 ps or to ∼2 ps for low and high laser intensities, respectively, the intensities of the Bragg-like peaks and Yoneda peak remain constant, while the position of the Bragg-like peaks steadily increases. This suggests that the embedded ML undergoes compression, although ML structure remains constant. In the third time frame, the intensities of the Bragg and Yoneda peaks are reduced, accompanied by a continued shift of the -position to even higher values. This indicates further compression of the entire ML, as well as significant modulation of the ML structure, resulting in the loss of correlation and subsequent decrease in the X-ray scattering intensity.
3. Discussion
To facilitate a more quantitative discussion of the ultrafast surface dynamics, we employed the BornAgain [
18] GISAXS analysis program (v 1.19) to model our experimental observations. The white circular dots in
Figure 4a represent the experimental data at
while the grey dots depict the signal at a delay of
. The solid lines represent the BornAgain model corresponding to the effective electron density profiles shown in
Figure 4d.
Figure 4b,c provide enlarged views of the intense peak at
and the surface-sensitive Yoneda region, respectively. According to the model, the observed variation between the blue (
) and red (
) lines in
Figure 4a–c can be attributed to the decrease in the Hurst parameter of the Al surface layer from
to
, indicating an increase in the spatial frequency of roughness, as illustrated schematically in
Figure 4e,f. This observation is supported by analyzing the out-of-plane signal along
at different
positions.
Figure 4g illustrates a lineout along
at the Al Yoneda peak. The decay starting at
is modeled via
where
m defines the slope of the decay which is proportional to the Hurst parameter
H [
19]. Refining the lineouts for 0 and 1 ps indeed reveals a reduction of the Hurst parameter as
m has decreased. Conversely, the lineout along
at larger
(
Figure 4h) shows no change in the Hurst parameter, indicating that the ultrafast change is localized to the surface rather than at the ML interfaces. Furthermore, it appears that in order to align the model with the experimental observation from 0 to
, the vertical correlation length needs to be increased from 100 to
, and the surface roughness RMS
should be slightly increased from
to
. We hypothesize that this ultrafast change in surface roughness properties to be attributed to the presence of a thin aluminium oxide layer on the surface, which typically forms within minutes when exposed to air [
20]. This oxide layer likely possesses a distinct surface morphology compared to the intrinsic surface morphology parameters of bare aluminium, leading to a reduction in spatial frequency of roughness as well as a decrease in vertical correlation length. Given that this layer is typically extremely thin, in the order of a nm, it evaporates instantaneously (<1 ps) after laser excitation revealing the intrinsic surface properties of the initial Al layer.
To substantiate our experimental findings further, we also conducted two advanced simulations: a 1D MULTI-fs hydrodynamic simulation and 1D PICLS Particle-in-Cell (PIC) simulation. The 1D MULTI-fs hydrodynamic simulation code is tailored for modeling short pulse high-intensity laser–solid interactions which directly solves the Maxwell equations. The code incorporates a temperature-dependent collision frequency, thermal conductivity ranging from metallic solids to high-temperature ideal plasmas, and a distinct equation-of-state (EOS) for electrons and ions (two-temperature model). The details of the simulation setups are summarized in
Appendix A section. The 1D PICLS simulation incorporates the collisional 1d3v (one-dimensional in space and three-dimensional in velocity) [
21]. An interpolated collision frequency akin to the MULTI-fs method was integrated into the code to address Angstrom-scale collisions occurring at electron temperatures (
) around the Fermi temperature (
).
Figure 5a,c display the density profile simulated using MULTI-fs at a delay of 0.5 ps. The localized heating induced by the laser within the skin depth instantaneously elevates the surface electron temperature to a few tens of eV and generates pressures exceeding >10 Mbar, prompting immediate surface expansion. This observation is aligned with the experiment demonstrating ultrafast changes in surface properties.
During the second time frame, spanning from 1 to ∼9 ps in the case of lower laser intensity (
), the intensity of the Bragg-like peak remains constant (
Figure 2b). As confirmed by MULTI-fs simulations, it appears that between 1 and
the compression wave is only in the Al layer (
Figure 5b), highlighted in yellow). Primarily, only the Al layer undergoes modulation during this interval, while the ML remains static with slight compression. Subsequently, after ∼8 ps, the simulation reveals a strong pressure wave compressing the ML. This aligns with the experimental findings (
Figure 2c) that demonstrated a further increase in the
-position of the Bragg peaks.
Figure 5c,d are the same simulation but performed at the higher laser intensity
. Here, the pressure peak arrives at the ML surface already after ∼3 ps, which encounters strong modulation subsequently. On the contrary, the experiment (
Figure 2b) shows the significant decrease in Bragg intensity already started at
. It implies that the significant modulation of the ML occurs before the arrival of the strong compression wave to the ML surface, which is not seen in the MULTI-fs simulation.
It is known that Lagrangian hydrodynamic simulations, such as MULTI-fs, are incapable of accurately simulating the atomic mixing of adjacent layers, which inhibits particle interdiffusion between Lagrangian cells used in the code [
9]. Given the ion thermal velocity (
1 nm/ps at ∼1 eV thermal temperature, where
is the Boltzmann constant,
is the ion temperature and
is the ion mass), particles are expected to penetrate into adjacent layers over timescales of several ps. To address this, we conducted 1D-PICLS kinetic simulation for the
case, summarized in
Figure 5e,f. While the simulation requires significantly higher computational cost, it properly includes the aforementioned kinetic effects. The simulation shows that the embedded ML undergoes significant layer intermixing from a thermal wave, even before the arrival of the strong compression wave. This phenomenon begins already at 1.5 ps delay as shown in
Figure 5e. Even though the pressure peak is still in the middle of the Al layer, the ML has already started to undergo modulation. This leads to reductions in the Bragg peak in the X-ray scattering. By 2.7 ps delay, when both MULTI-fs and PICLS indicate the arrival of the pressure peak (
Figure 5f), PICLS demonstrates that the entire ML is intermixed, and the original ML structure is entirely disrupted. Such a structure cannot produce any pronounced Bragg-like peaks due to the absence of a statistically relevant periodic structure. Note that the cooling is overestimated in PICLS (compare
Figure 5d,f, highlighted in the yellow-shaded area) due to the absence of an opacity package in the code, which reabsorbs Bremsstrahlung radiation, as discussed in [
9]. This results in reduced particle kinetic energy and is likely responsible for the density peaks observed at approximately 170 nm inside the Al layer in
Figure 5f.
In summary, for the low-intensity case , modifications to the ML structure are primarily driven by the compression wave, as the electron temperature is low enough that the thermal wave cannot significantly alter the ML before the arrival of the pressure wave. Conversely, for the high-intensity case , significant heating occurs deep within the ML much before the arrival of the pressure wave.