Next Article in Journal
Investigating the Electronic Properties and Stability of Rh3 Clusters on Rutile TiO2 for Potential Photocatalytic Applications
Next Article in Special Issue
Versatile Porphyrin Arrangements for Photodynamic Therapy—A Review
Previous Article in Journal
Enhancing Optical and Thermal Stability of Blue-Emitting Perovskite Nanocrystals through Surface Passivation with Sulfonate or Sulfonic Acid Ligands
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

(Sub-)Picosecond Surface Correlations of Femtosecond Laser Excited Al-Coated Multilayers Observed by Grazing-Incidence X-ray Scattering

by
Lisa Randolph
1,2,*,
Mohammadreza Banjafar
2,3,
Toshinori Yabuuchi
4,5,
Carsten Baehtz
6,
Michael Bussmann
6,7,
Nicholas P. Dover
8,
Lingen Huang
6,
Yuichi Inubushi
4,5,
Gerhard Jakob
9,
Mathias Kläui
9,
Dmitriy Ksenzov
1,
Mikako Makita
2,
Kohei Miyanishi
5,
Mamiko Nishiuchi
10,
Özgül Öztürk
1,
Michael Paulus
11,
Alexander Pelka
6,
Thomas R. Preston
2,
Jan-Patrick Schwinkendorf
2,6,
Keiichi Sueda
5,
Tadashi Togashi
4,5,
Thomas E. Cowan
3,6,
Thomas Kluge
6,
Christian Gutt
1 and
Motoaki Nakatsutsumi
2,*
add Show full author list remove Hide full author list
1
Department Physik, University of Siegen, 57072 Siegen, Germany
2
European XFEL, 22869 Schenefeld, Germany
3
Fakultät Physik, TU Dresden, 01069 Dresden, Germany
4
Japan Synchrotron Radiation Research Institute (JASRI), Sayo 679-5198, Hyogo, Japan
5
RIKEN SPring-8 Center, Sayo 679-5148, Hyogo, Japan
6
Helmholtz-Zentrum Dresden-Rossendorf, 01328 Dresden, Germany
7
Center for Advanced Systems Understanding (CASUS), 02826 Görlitz, Germany
8
The John Adams Institute for Accelerator Science, Imperial College London, London SW7 2BW, UK
9
Institute of Physics, Johannes Gutenberg-University, 55099 Mainz, Germany
10
Kansai Photon Science Institute, National Institutes for Quantum Science and Technology, Kyoto 619-0215, Japan
11
Fakultät Physik/DELTA, TU Dortmund, 44221 Dortmund, Germany
*
Authors to whom correspondence should be addressed.
Nanomaterials 2024, 14(12), 1050; https://doi.org/10.3390/nano14121050
Submission received: 30 April 2024 / Revised: 24 May 2024 / Accepted: 6 June 2024 / Published: 19 June 2024

Abstract

:
Femtosecond high-intensity laser pulses at intensities surpassing 1014 W/cm2 can generate a diverse range of functional surface nanostructures. Achieving precise control over the production of these functional structures necessitates a thorough understanding of the surface morphology dynamics with nanometer-scale spatial resolution and picosecond-scale temporal resolution. In this study, we show that single XFEL pulses can elucidate structural changes on surfaces induced by laser-generated plasmas using grazing-incidence small-angle X-ray scattering (GISAXS). Using aluminium-coated multilayer samples we distinguish between sub-picosecond (ps) surface morphology dynamics and subsequent multi-ps subsurface density dynamics with nanometer-depth sensitivity. The observed subsurface density dynamics serve to validate advanced simulation models representing matter under extreme conditions. Our findings promise to open new avenues for laser material-nanoprocessing and high-energy-density science.

1. Introduction

Intense, ultrashort laser–solid interactions at intensities ranging from 1013 to 1016 W/cm2 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 ( μ m ) [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 ( Q z ), 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 Q y 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
C ( R ) = < h ( 0 ) h ( R ) > = σ 2 exp ( R / ξ ) 2 H
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
< h j ( 0 ) h k ( R ) > = 1 2 σ k σ j C j ( R ) + σ j σ k C k ( R ) × exp z j z k ξ .
Here, C ( R ) 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, 4.5   n m ) and copper nitride (Cu3N, 8.5   n m ) capped with a 200   n m thick Aluminum top layer was irradiated by an optical laser with a central wavelength of 800   n m , intensities of 8 × 1014 and 8 × 1015 W/cm2, pulse duration of 40   f s 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 8.81   k e V . The X-ray beam with a pulse energy of ∼0.1 mJ/pulse was incident under a grazing incidence angle of 0.75 deg 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, Q specular = 1.16   n m 1 ) 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 4   μ m 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 ( Q y = 0 ) as a function of Q z 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 ( Q z ) and surface-parallel ( Q x ) directions. As a result, the scattering signal reflects vertical density correlations along the z-direction ( Q z ) of the ML structure, roughness correlations along the surface plane ( Q x ) direction, and cross-correlations between different interfaces. For clarity, only the Q z values are plotted here, as Q x Q z . The in-plane signal at Q z > Q specular is closely associated with the specular reflectivity curve [10], which we characterized ex situ. The spacing of the intense Bragg-like peaks at Q z = 1.6 , 2.1 and 2.6   n m 1 corresponds to the typical length scale, i.e., 13   n m thickness of each Ta/Cu3N 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 12   p s . 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 12   p s , 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 Q z = 1.6 , 2.1 and 2.6   n m 1 . Figure 2a illustrates the temporal evolution of the peak at Q z = 2.1   n m 1 following laser irradiation with an intensity of 8 × 1014 W/cm2. While only the peak at Q z = 2.1   n m 1 is displayed here, qualitatively similar behaviour is observed for the peaks at Q z = 1.6 and 2.6   n m 1 . 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
Intensity = A · erf ( x t 0 ) Δ t B
resulting in values for t 0 = 11.4 ± 0.3   p s and a width of Δ t = 2.1 ± 0.4   p s . Here, A represents a stretching factor, t 0 denotes the time at which the error function decreased to half of its initial value, Δ t 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 × 1015 W/cm2, 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 t 0 = 2.8 ± 0.3   p s and a width of Δ t = 1.0 ± 0.5   p s . 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
Q z time delay delay c
we result in an exponent c = 2.0 ± 0.4 for the lower intensity and c = 2.5 ± 0.7 for the higher intensity, respectively.
The dynamic diffraction effects at Q z < Q specular 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— Q z = 0.75   n m 1 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 Q z = 0.75   n m 1 exhibits a progressive decrease after laser irradiation, which persists for at least at t = 12   p s before almost disappearing at t = 14   p s , indicative of the ongoing presence of the solid density Al cover layer. The surface ablation speed can be approximated by the speed of sound, C s = ( γ e Z m e a n k B T e + γ i k B T i ) / M i , where T e and T i are electron and ion temperature, respectively, Z m e a n is the mean ionization, γ e 1 and γ i 3 are adiabatic index of electrons and ions, respectively and M i = 26 a.u. is the ion mass for Al [17]. In our temperature range ( T e = 5–10 eV, T i = 1–5 eV) the ablation of the 200   n m 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 ( Q z = 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 Q z shown in Figure 2b. In Figure 3a, two additional Yoneda peaks are also visible at higher Q z , corresponding to the Yoneda peaks for Cu 3 N ( Q z = 0.86   n m 1 ) and Ta ( 0.95   n m 1 ), respectively. Given that these materials are embedded within a 200   n m 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 14   p s 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 Cu3N 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 Q z -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 Q z 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 Q z -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 0   p s while the grey dots depict the signal at a delay of 1   p s . 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 2.1   n m 1 and the surface-sensitive Yoneda region, respectively. According to the model, the observed variation between the blue ( 0   p s ) and red ( 1   p s ) lines in Figure 4a–c can be attributed to the decrease in the Hurst parameter of the Al surface layer from 0.6 to 0.15 , 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 Q y at different Q z positions. Figure 4g illustrates a lineout along Q y at the Al Yoneda peak. The decay starting at Q y = 0.01   n m 1 is modeled via
Intensity Q y m
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 Q y at larger Q z = 2.1   n m 1 (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 1   p s , the vertical correlation length needs to be increased from 100 to 150   n m , and the surface roughness RMS σ should be slightly increased from 2.3 to 2.6   n m . 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 ( T e ) around the Fermi temperature ( T F ).
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 ( I L = 8 × 10 14   W / c m 2 ), the intensity of the Bragg-like peak remains constant (Figure 2b). As confirmed by MULTI-fs simulations, it appears that between 1 and 8   p s 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 Q z -position of the Bragg peaks.
Figure 5c,d are the same simulation but performed at the higher laser intensity I L = 8 × 10 15   W / c m 2 . 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 2   p s . 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 ( k B T i / M i 1 nm/ps at ∼1 eV thermal temperature, where k B is the Boltzmann constant, T i is the ion temperature and M i 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 I L = 8 × 10 15   W / c m 2 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 I L = 8 × 10 14   W / c m 2 , 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 I L = 8 × 10 15   W / c m 2 , significant heating occurs deep within the ML much before the arrival of the pressure wave.

4. Conclusions

In conclusion, our investigation of surface and subsurface dynamics in high-intensity laser-excited multi-layer (ML) systems using ultrafast grazing-incidence small-angle X-ray scattering (GISAXS) has provided important insights into the complex physical processes occurring in laser-matter interactions on a picosecond (ps) time scale.
By analyzing the GISAXS pattern, we distinguish between ultrafast sub-ps surface dynamics and subsequent multi-ps subsurface dynamics. The ultrafast modifications in surface properties, induced by strong localized surface heating, are evident in Yoneda region at low Q, where the scattering exit angle is close to the external total reflection. In contrast, the slower subsurface dynamics manifest as changes in intensity and position of the intense Bragg-like peak at Q z = 2.1   n m 1 . Comparison of the temporal evolution of the density dynamics with 1D hydrodynamic and 1D particle-in-cell simulations reveals qualitative agreement across different time regimes observed in the experiment. Notably, the observed time scales align well with experimental observations, indicating the arrival of the pressure wave at the ML surface embedded in the 200 nm thick aluminium within a range of 3 to 10 ps and the intermixing of individual layers.
The presented GISAXS experiment was limited in time resolution to ∼1.2 ps due to the large footprint in grazing-incidence. However, this can be improved by using a smaller X-ray focal spot, e.g., 100   n m , which is readily available at XFEL facilities. This would imply only ∼30 fs smearing at 0.75° grazing-incidence, allowing experimental studies on surface and subsurface dynamics with nm spatial and femtosecond temporal resolution to observe the ultrafast changes on the surface or using laser intensities above 1016 W/cm2 to study, e.g., the effect of density oscillations [22].

Author Contributions

Conceptualization, M.N. (Motoaki Nakatsutsumi) and C.G.; software, L.R., M.B. (Mohammadreza Banjafar), L.H. and M.B. (Michael Bussmann); validation, D.K. and Ö.Ö.; formal analysis, L.R. and M.B. (Mohammadreza Banjafar); investigation, L.R., M.B. (Mohammadreza Banjafar), C.B., N.P.D., M.M., M.N. (Mamiko Nishiuchi), A.P., T.R.P., J.-P.S., T.T., T.K., C.G. and M.N. (Motoaki Nakatsutsumi); resources, T.Y., Y.I., K.M., K.S., G.J., M.K. and M.P.; writing—original draft preparation, L.R., M.B. (Mohammadreza Banjafar), C.G. and M.N. (Motoaki Nakatsutsumi); writing—review and editing, all authors; supervision, M.N. (Motoaki Nakatsutsumi), C.G., T.K. and T.E.C.; funding acquisition, G.J., M.K., C.G. and M.N. (Motoaki Nakatsutsumi). All authors have read and agreed to the published version of the manuscript.

Funding

C.G. and M.N. (Motoaki Nakatsutsumi) acknowledge funding by Deutsche Forschungsgemeinschaft (DFG, German Research Foundation) Project No. 436265460 (GU 535/6-1). G.J. and M.K. acknowledge the support by the DFG Project No. 268565370 (SFB TRR173 Projects A01 and B02) by TopDyn and the ForLab MagSens.

Data Availability Statement

The data are available upon reasonable request.

Acknowledgments

The XFEL experiments were performed at the BL2 of SACLA with the approval of the Japan Synchrotron Radiation Research Institute (JASRI) (Proposal No. 2019B8076). We thank E. Brambrink, S. Göde, S.V. Rahul, C. Rödel, A. Kon, J. Koga and Y. Sentoku for various discussions and advices.

Conflicts of Interest

The authors declare no conflicts of interest.

Appendix A

Appendix A.1. Experimental Setup Details, X-ray and Laser Parameters, Timing Synchronization

The experiment was performed at the SACLA XFEL facility in Japan at the Experimental Hutch 6, which provides a high-intensity optical laser ( λ L = 800   n m central wavelength, maximum 10   J / pulse with 40   f s duration in full-width half-maximum (FWHM)) combined with ultrashort intense X-ray pulses [13]. The X-ray pulses had a photon energy of 8.81   k e V (with 43   e V FWHM bandwidth), ∼100 μ J/pulse, and a pulse duration of 7   f s in FWHM. The X-rays were focused to a 4   μ m FWHM spot on sample by a set of compound refractive lenses placed 3   m upstream from the sample. The scattered X-ray signal was recorded on an MPCCD area detector with 50   μ m pixel size [14] placed at a distance of 1.24   m from the sample and shielded by a 50   μ m thick Al foil, to remove plasma-induced bremsstrahlung background. The incident angle was fixed at 0.75° to be slightly larger than the critical angle of total external reflection α c n e for all materials, Ta ( α c _ Ta = 0.46 ), Cu ( α c _ Cu = 0.35 ) and Al ( α c _ Al = 0.22 ).
The samples were placed in a vacuum chamber and then irradiated by a high-intensity optical laser, impinging on the sample at 17 incident angle from the surface normal with p-polarization. In order to cover the X-ray footprint on sample ( 4   μ m FWHM for 0.75 grazing incidence yields 300   μ m ), the optical laser beam was defocused to a diameter of ∼500  μ m yielding an average laser intensity of about 8 × 10 14   W / c m 2 and 8 × 10 15   W / c m 2 depending on the attenuation. The synchronization between the X-ray and the laser was measured before starting of the experiment [13].

Appendix A.2. Multi-fs

The 1D MULTI-fs hydrodynamic simulation code is tailored for modeling short (≤ps) pulse high-intensity laser-solid interactions at intensities below <1017 W / cm 2 [23,24]. MULTI-fs directly solves the Maxwell equations, capturing ultrashort laser plasma dynamics with sharp plasma density gradients. 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). In the MULTI-fs simulation setup, the simulation box comprises 1478 cells distributed as follows: 400 cells for the 200 nm thick Al top layer, 678 cells for the ML, and 400 cells for a 200 nm thick Al substrate. To mitigate numerical artifacts, finer cell resolutions were employed near layer interfaces. We employed the EOS and ionization tables for aluminium (Al) as provided by the MULTI-fs package. For tantalum (Ta) and copper (Cu), the EOS and ionization data were generated using the FEOS code [25]. The radiation transfer module was disabled in the simulation, and the free streaming limiting factor was set to f = 0.6 . To maintain consistency with the experimental setup, the laser angle of incidence was fixed at 17 from the sample normal. In our simulations, Al was chosen as the substrate material instead of silicon (Si) used in the experiment. This choice was motivated by well-established EOS of aluminium. Given that the energy transfer occurs from the ML to the substrate, we anticipate that the choice between these two materials will have minimal impact on the dynamics of the ML.

Appendix A.3. 1D PICLS

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 ( T e ) around the Fermi temperature ( T F ). To accurately simulate the microscopic particle collisions at the atomic scale, the spatial resolution was configured with a cell size of Δ x = 2 Å corresponding to a time–step of Δ t = 6.67 × 10 19 s. Each computational cell accommodated 30 virtual ion particles with initial charge states of 2, 1, and 3 for Ta, Cu, and Al, respectively. Additionally, a fourth-order particle shape was employed, along with distinct particle weightings for different ion species. The ion number densities are set to realistic values for Ta ( n Ta = 31.8 n c ), Cu ( n Cu = 48.7 n c ), and Al ( n Al = 34.6 n c ), where n c = m e ω L 2 / ( 4 π e 2 ) = 1.742 × 10 21 cm 3 represents the critical plasma density at the laser wavelength of λ L = 800 nm. Here, m e and e denote the electron mass and charge, respectively, and ω L is the laser angular frequency. Ionization dynamics are modeled employing field and direct-impact ionization models. The laser incident angle is set as normal to the surface, as oblique incidence is not supported in 1D PICLS. Our test simulations using the MULTI-fs have revealed a negligible difference, approximately 1%, in laser absorption between normal incidence and a 17 angle of incidence.

References

  1. Bonse, J.; Gräf, S. Maxwell Meets Marangoni—A Review of Theories on Laser-Induced Periodic Surface Structures. Laser Photonics Rev. 2020, 14, 2000215. [Google Scholar] [CrossRef]
  2. Lutey, A.H.; Gemini, L.; Romoli, L.; Lazzini, G.; Fuso, F.; Faucon, M.; Kling, R. Towards laser-textured antibacterial surfaces. Sci. Rep. 2018, 8, 10112. [Google Scholar] [CrossRef]
  3. San-Blas, A.; Martinez-Calderon, M.; Buencuerpo, J.; Sanchez-Brea, L.M.; del Hoyo, J.; Gómez-Aranzadi, M.; Rodríguez, A.; Olaizola, S.M. Femtosecond laser fabrication of LIPSS-based waveplates on metallic surfaces. Appl. Surf. Sci. 2020, 520, 146328. [Google Scholar] [CrossRef]
  4. Li, C.Y.; Duan, S.; Yi, J.; Wang, C.; Radjenovic, P.M.; Tian, Z.Q.; Li, J.F. Real-time detection of single-molecule reaction by plasmon-enhanced spectroscopy. Sci. Adv. 2020, 6, eaba6012. [Google Scholar] [CrossRef]
  5. Bonse, J.; Kirner, S.V.; Griepentrog, M.; Spaltmann, D.; Krüger, J. Femtosecond laser texturing of surfaces for tribological applications. Materials 2018, 11, 801. [Google Scholar] [CrossRef]
  6. Rudenko, A.; Abou-Saleh, A.; Pigeon, F.; Mauclair, C.; Garrelie, F.; Stoian, R.; Colombier, J. High-frequency periodic patterns driven by non-radiative fields coupled with Marangoni convection instabilities on laser-excited metal surfaces. Acta Mater. 2020, 194, 93–105. [Google Scholar] [CrossRef]
  7. Roth, S.V. A deep look into the spray coating process in real-time—The crucial role of X-rays. J. Phys. Condens. Matter 2016, 28, 403003. [Google Scholar] [CrossRef]
  8. Su, B.; Körstgens, V.; Yao, Y.; Magerl, D.; Song, L.; Metwalli, E.; Bernstorff, S.; Müller-Buschbaum, P. Pore size control of block copolymer-templated sol–gel-synthesized titania films deposited via spray coating. J. Sol-Gel Sci. Technol. 2017, 81, 346–354. [Google Scholar] [CrossRef]
  9. Randolph, L.; Banjafar, M.; Preston, T.R.; Yabuuchi, T.; Makita, M.; Dover, N.P.; Rödel, C.; Göde, S.; Inubushi, Y.; Jakob, G.; et al. Nanoscale subsurface dynamics of solids upon high-intensity femtosecond laser irradiation observed by grazing-incidence X-ray scattering. Phys. Rev. Res. 2022, 4, 033038. [Google Scholar] [CrossRef]
  10. Holý, V.; Kuběna, J.; Ohlídal, I.; Lischka, K.; Plotz, W. X-ray reflection from rough layered systems. Phys. Rev. B 1993, 47, 15896–15903. [Google Scholar] [CrossRef]
  11. Sinha, S.K.; Sirota, E.B.; Garoff, S.; Stanley, H.B. X-ray and neutron scattering from rough surfaces. Phys. Rev. B 1988, 38, 2297–2311. [Google Scholar] [CrossRef]
  12. Schlomka, J.P.; Tolan, M.; Schwalowsky, L.; Seeck, O.H.; Stettner, J.; Press, W. X-ray diffraction from Si/Ge layers: Diffuse scattering in the region of total external reflection. Phys. Rev. B 1995, 51, 2311–2321. [Google Scholar] [CrossRef]
  13. Yabuuchi, T.; Kon, A.; Inubushi, Y.; Togahi, T.; Sueda, K.; Itoga, T.; Nakajima, K.; Habara, H.; Kodama, R.; Tomizawa, H.; et al. An experimental platform using high-power, high-intensity optical lasers with the hard X-ray free-electron laser at SACLA. J. Sync. Rad. 2019, 26, 585–594. [Google Scholar] [CrossRef]
  14. Kameshima, T.; Ono, S.; Kudo, T.; Ozaki, K.; Kirihara, Y.; Kobayashi, K.; Inubushi, Y.; Yabashi, M.; Horigome, T.; Holland, A.; et al. Development of an X-ray pixel detector with multi-port charge-coupled device for X-ray free-electron laser experiments. Rev. Sci. Instrum. 2014, 85, 033110. [Google Scholar] [CrossRef]
  15. Kiessig, H. Untersuchungen zur Totalreflexion von Röntgenstrahlen. Ann. Phys. 1931, 10, 769. [Google Scholar] [CrossRef]
  16. Yoneda, Y. Anomalous Surface Reflection of X-rays. Phys. Rev. 1963, 131, 2010–2013. [Google Scholar] [CrossRef]
  17. Drake, R.P. High-Energy-Density Physics; Graduate Texts in Physics; Springer: Cham, Switzerland, 2018. [Google Scholar]
  18. Pospelov, G.; Van Herck, W.; Burle, J.; Carmona Loaiza, J.M.; Durniak, C.; Fisher, J.M.; Ganeva, M.; Yurov, D.; Wuttke, J. BornAgain: Software for simulating and fitting grazing-incidence small-angle scattering. J. Appl. Crystallogr. 2020, 53, 262–276. [Google Scholar] [CrossRef]
  19. Salditt, T.; Metzger, T.H.; Peisl, J. Kinetic Roughness of Amorphous Multilayers Studied by Diffuse X-ray Scattering. Phys. Rev. Lett. 1994, 73, 2228–2231. [Google Scholar] [CrossRef]
  20. Nguyen, L.; Hashimoto, T.; Zakharov, D.N.; Stach, E.A.; Rooney, A.P.; Berkels, B.; Thompson, G.E.; Haigh, S.J.; Burnett, T.L. Atomic-Scale Insights into the Oxidation of Aluminum. ACS Appl. Mater. Interfaces 2018, 10, 2230–2235. [Google Scholar] [CrossRef] [PubMed]
  21. Sentoku, Y.; Kemp, A. Numerical methods for particle simulations at extreme densities and temperatures: Weighted particles, relativistic collisions and reduced currents. J. Comput. Phys. 2008, 227, 6846–6861. [Google Scholar] [CrossRef]
  22. Paschke-Bruehl, F.L.; Banjafar, M.; Garten, M.; Huang, L.; Marré, B.E.; Nakatsutsumi, M.; Randolph, L.; Cowan, T.E.; Schramm, U.; Kluge, T. Heating in Multi-Layer Targets at ultra-high Intensity Laser Irradiation and the Impact of Density Oscillation. New J. Phys. 2023, 25, 043024. [Google Scholar] [CrossRef]
  23. Eidmann, K.; Meyer-ter Vehn, J.; Schlegel, T.; Hüller, S. Hydrodynamic simulation of subpicosecond laser interaction with solid-density matter. Phys. Rev. E 2000, 62, 1202–1214. [Google Scholar] [CrossRef]
  24. Ramis, R.; Eidmann, K.; Meyer-ter Vehn, J.; Hüller, S. MULTI-fs—A computer code for laser–plasma interaction in the femtosecond regime. Comput. Phys. Commun. 2012, 183, 637–655. [Google Scholar] [CrossRef]
  25. Faik, S.; Tauschwitz, A.; Iosilevskiy, I. The equation of state package FEOS for high energy density matter. Comput. Phys. Comm. 2018, 227, 117–125. [Google Scholar] [CrossRef]
Figure 1. Schematics of the experimental setup to investigate the surface and subsurface solid-density plasma dynamics with GISAXS using single femtosecond X-ray FEL pulses. (a) An ML sample consists of 5 repetitions of tantalum (Ta) and copper nitride (Cu3N) of 4.5 and 8.5   n m thickness covered by a 200   n m aluminium capping layer. The samples were irradiated by an optical laser with a central wavelength of 800   n m , intensity of ( 8 ± 0.2 ) × 10 14 or ( 8 ± 0.2 ) × 10 15   W / c m 2 with 40   f s pulse duration. The laser is irradiated by an incident angle of 17 from the surface normal in p polarization. The X-ray pulses with 8.81   k e V photon energy, 7   f s FWHM duration are irradiated on the sample at the grazing-incidence angle of α i = 0.75 , i.e., slightly above the critical angle of external total reflection of the layer materials. The laser beam is defocused to obtain a ∼500 μ m spot diameter to cover the X-ray footprint of ∼300 μ m at the surface. Scattered photons are recorded by 2 modules of the MPCCD area detector placed around the specular direction. The strong specular peak at Q z = 1.16   n m 1 is blocked. (b) In-plane signal along Q z (at Q y = 0 ) for different time delays between 0 and 14 ps after the laser intensity peak.
Figure 1. Schematics of the experimental setup to investigate the surface and subsurface solid-density plasma dynamics with GISAXS using single femtosecond X-ray FEL pulses. (a) An ML sample consists of 5 repetitions of tantalum (Ta) and copper nitride (Cu3N) of 4.5 and 8.5   n m thickness covered by a 200   n m aluminium capping layer. The samples were irradiated by an optical laser with a central wavelength of 800   n m , intensity of ( 8 ± 0.2 ) × 10 14 or ( 8 ± 0.2 ) × 10 15   W / c m 2 with 40   f s pulse duration. The laser is irradiated by an incident angle of 17 from the surface normal in p polarization. The X-ray pulses with 8.81   k e V photon energy, 7   f s FWHM duration are irradiated on the sample at the grazing-incidence angle of α i = 0.75 , i.e., slightly above the critical angle of external total reflection of the layer materials. The laser beam is defocused to obtain a ∼500 μ m spot diameter to cover the X-ray footprint of ∼300 μ m at the surface. Scattered photons are recorded by 2 modules of the MPCCD area detector placed around the specular direction. The strong specular peak at Q z = 1.16   n m 1 is blocked. (b) In-plane signal along Q z (at Q y = 0 ) for different time delays between 0 and 14 ps after the laser intensity peak.
Nanomaterials 14 01050 g001
Figure 2. In-plane scattering signal along Q z around the Bragg-like peak at 2.1   n m 1 . (a) Circular dots represent the intensity of experimental data for various time delays, while solid lines correspond to a Gaussian model fit. (b) Peak intensity derived from the Gaussian model plotted against time delay. Circular dots denote the maximum intensity values from (a), while the dashed line represents a refinement using an error function. Blue color denotes the lower laser intensity of 8 × 10 14   W / c m 2 while red indicates the higher laser intensity of 8 × 10 15   W / c m 2 . (c) Q z -position of the Bragg-like peak. The dashed lines represent refinements using an exponential model.
Figure 2. In-plane scattering signal along Q z around the Bragg-like peak at 2.1   n m 1 . (a) Circular dots represent the intensity of experimental data for various time delays, while solid lines correspond to a Gaussian model fit. (b) Peak intensity derived from the Gaussian model plotted against time delay. Circular dots denote the maximum intensity values from (a), while the dashed line represents a refinement using an error function. Blue color denotes the lower laser intensity of 8 × 10 14   W / c m 2 while red indicates the higher laser intensity of 8 × 10 15   W / c m 2 . (c) Q z -position of the Bragg-like peak. The dashed lines represent refinements using an exponential model.
Nanomaterials 14 01050 g002
Figure 3. In-plane signal along Q z around the dynamical diffraction area between 0.6 and 1.0   n m 1 . (a) In-plane signal for different time delays at lower laser intensity case ( 8 × 10 14   W / c m 2 ). (b) The circular dots display the integrated signal between Q z = 0.65–0.79 nm−1 as a function of time delay. The dashed line represents a guide to the eyes. Blue color stands for the lower laser intensity of 8 × 10 14   W / c m 2 and red denote the higher laser intensity 8 × 10 15   W / c m 2 .
Figure 3. In-plane signal along Q z around the dynamical diffraction area between 0.6 and 1.0   n m 1 . (a) In-plane signal for different time delays at lower laser intensity case ( 8 × 10 14   W / c m 2 ). (b) The circular dots display the integrated signal between Q z = 0.65–0.79 nm−1 as a function of time delay. The dashed line represents a guide to the eyes. Blue color stands for the lower laser intensity of 8 × 10 14   W / c m 2 and red denote the higher laser intensity 8 × 10 15   W / c m 2 .
Nanomaterials 14 01050 g003
Figure 4. GISAXS signals and corresponding retrieved real-space electron density profile. (a) Lineout at 0   p s (white circular dots) and 1   p s (grey circular dots) after the laser intensity peak. Solid lines represent models using the program BornAgain (https://bornagainproject.org/). (b) Zoom into the area of the Bragg-like peak at Q z = 2.1   n m 1 and (c) zoom into the area of dynamical diffraction around Q z = 0.75   n m 1 . (d) Retrieved real-space effective electron density profile as a function of the depth (Z). The laser irradiates the sample from the left. (e) Generated surface model for H = 0.6 , σ = 2.3   n m and (f) H = 0.15 , σ = 2.6   n m . (g) Lineout along Q y at the Al Yoneda peak ( Q z = 0.75   n m 1 ). Dashed lines indicate a fit of the decay. (h) Lineout along Q y at Q z = 2.1   n m 1 .
Figure 4. GISAXS signals and corresponding retrieved real-space electron density profile. (a) Lineout at 0   p s (white circular dots) and 1   p s (grey circular dots) after the laser intensity peak. Solid lines represent models using the program BornAgain (https://bornagainproject.org/). (b) Zoom into the area of the Bragg-like peak at Q z = 2.1   n m 1 and (c) zoom into the area of dynamical diffraction around Q z = 0.75   n m 1 . (d) Retrieved real-space effective electron density profile as a function of the depth (Z). The laser irradiates the sample from the left. (e) Generated surface model for H = 0.6 , σ = 2.3   n m and (f) H = 0.15 , σ = 2.6   n m . (g) Lineout along Q y at the Al Yoneda peak ( Q z = 0.75   n m 1 ). Dashed lines indicate a fit of the decay. (h) Lineout along Q y at Q z = 2.1   n m 1 .
Nanomaterials 14 01050 g004
Figure 5. MULTI-fs hydrodynamic simulations irradiated by a p-polarized laser pulse with a duration of 50 fs, a wavelength of 800 nm, and a laser intensity of (a,b) 8 × 10 14   W / c m 2 and (c,d) 8 × 10 15   W / c m 2 , respectively. Time delays are (a,c) 0.5 ps, (b) 7.75 ps and (d) 2.7 ps, respectively. The ML structure is the same as the experiment. (e,f) 1D Particle-in-cell PICLS simulation results at time delays of (e) 1.5 and (f) 2.7 ps, respectively. In all figures, the yellow, grey, and red regions represent the laser-induced pressure wave (in Mbar or Mbar / 10 ), mass density (in g / cm 3 ), and electron temperature (in eV), respectively. The dashed line denotes the initial sample density. The red arrow in (a) indicates the direction of the incoming laser.
Figure 5. MULTI-fs hydrodynamic simulations irradiated by a p-polarized laser pulse with a duration of 50 fs, a wavelength of 800 nm, and a laser intensity of (a,b) 8 × 10 14   W / c m 2 and (c,d) 8 × 10 15   W / c m 2 , respectively. Time delays are (a,c) 0.5 ps, (b) 7.75 ps and (d) 2.7 ps, respectively. The ML structure is the same as the experiment. (e,f) 1D Particle-in-cell PICLS simulation results at time delays of (e) 1.5 and (f) 2.7 ps, respectively. In all figures, the yellow, grey, and red regions represent the laser-induced pressure wave (in Mbar or Mbar / 10 ), mass density (in g / cm 3 ), and electron temperature (in eV), respectively. The dashed line denotes the initial sample density. The red arrow in (a) indicates the direction of the incoming laser.
Nanomaterials 14 01050 g005
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Randolph, L.; Banjafar, M.; Yabuuchi, T.; Baehtz, C.; Bussmann, M.; Dover, N.P.; Huang, L.; Inubushi, Y.; Jakob, G.; Kläui, M.; et al. (Sub-)Picosecond Surface Correlations of Femtosecond Laser Excited Al-Coated Multilayers Observed by Grazing-Incidence X-ray Scattering. Nanomaterials 2024, 14, 1050. https://doi.org/10.3390/nano14121050

AMA Style

Randolph L, Banjafar M, Yabuuchi T, Baehtz C, Bussmann M, Dover NP, Huang L, Inubushi Y, Jakob G, Kläui M, et al. (Sub-)Picosecond Surface Correlations of Femtosecond Laser Excited Al-Coated Multilayers Observed by Grazing-Incidence X-ray Scattering. Nanomaterials. 2024; 14(12):1050. https://doi.org/10.3390/nano14121050

Chicago/Turabian Style

Randolph, Lisa, Mohammadreza Banjafar, Toshinori Yabuuchi, Carsten Baehtz, Michael Bussmann, Nicholas P. Dover, Lingen Huang, Yuichi Inubushi, Gerhard Jakob, Mathias Kläui, and et al. 2024. "(Sub-)Picosecond Surface Correlations of Femtosecond Laser Excited Al-Coated Multilayers Observed by Grazing-Incidence X-ray Scattering" Nanomaterials 14, no. 12: 1050. https://doi.org/10.3390/nano14121050

APA Style

Randolph, L., Banjafar, M., Yabuuchi, T., Baehtz, C., Bussmann, M., Dover, N. P., Huang, L., Inubushi, Y., Jakob, G., Kläui, M., Ksenzov, D., Makita, M., Miyanishi, K., Nishiuchi, M., Öztürk, Ö., Paulus, M., Pelka, A., Preston, T. R., Schwinkendorf, J. -P., ... Nakatsutsumi, M. (2024). (Sub-)Picosecond Surface Correlations of Femtosecond Laser Excited Al-Coated Multilayers Observed by Grazing-Incidence X-ray Scattering. Nanomaterials, 14(12), 1050. https://doi.org/10.3390/nano14121050

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop