Tracking picosecond strain pulses in heterostructures that exhibit giant magnetostriction

We combine ultrafast X-ray diffraction (UXRD) and time-resolved Magneto-Optical Kerr Effect (MOKE) measurements to monitor the strain pulses in laser-excited TbFe2/Nb heterostructures. Spatial separation of the Nb detection layer from the laser excitation region allows for a background-free characterization of the laser-generated strain pulses. We clearly observe symmetric bipolar strain pulses if the excited TbFe2 surface terminates the sample and a decomposition of the strain wavepacket into an asymmetric bipolar and a unipolar pulse, if a SiO2 glass capping layer covers the excited TbFe2 layer. The inverse magnetostriction of the temporally separated unipolar strain pulses in this sample leads to a MOKE signal that linearly depends on the strain pulse amplitude measured through UXRD. Linear chain model simulations accurately predict the timing and shape of UXRD and MOKE signals that are caused by the strain reflections from multiple interfaces in the heterostructure.


I. INTRODUCTION
The generation, propagation, and detection of laser-induced strain waves contain rich physics that has been studied extensively since the seminal work of Thomsen et al., which exploited photoelasticity for detection. 1,2 Ingenious all-optical probing schemes for GHz to THz phonons have since then been used to investigate strain waves in multiple materials, [3][4][5][6] vibrational modes of nanoparticles, 7 shear waves, 8 nonlinear propagation effects, 9,10 and acoustic solitons. 11,12 Strain waves that originate from coherent phonon excitation can attain transient stresses on the order of GPa, which have been shown to interact with other phenomena such as phase transitions, 13,14 quantum well bandgaps, 15 piezo-16 /ferroelectricity, 17 and magnetism. [18][19][20][21][22] Such interactions are not only of fundamental interest but may also become relevant for applications as soon as the understanding allows for controllability. 23,24 In order to study the response to pure strain pulses, it is beneficial to spatially separate the laser excited transducer from the probed layer. The inherent limitation given by the finite optical penetration of the visible light in the transducer is often circumvented by backside probing schemes. 10,21 The development of (sub)-picosecond hard X-ray diffraction has opened the possibility to directly obtain the timeresolved strain amplitude [25][26][27][28] with penetration depths in the few mm regime. The separation of the layer peaks in reciprocal space allows for material specific probing of the energy flow 29,30 and strain evolution 31,32 in nanoscopic, crystalline heterostructures.
The envisioned manipulation of the polarization and magnetization states in ferroic materials via strain relies on a strong coupling between spin or electronic degrees of freedom and the atomic lattice. 23,33 In this regard, rare-earth-based alloys such as Terfenol (TbFe 2 ) have attracted attention due to the discovery of "giant magnetostriction," 34,35 i.e., lattice strains in excess of 10 À3 caused by magnetization change. Among the binary rare-earth alloys, TbFe 2 exhibits the largest magnetostriction. 36 It combines the large exchange interaction of the 3d orbital in Iron (Fe) with the large magnetic moment of 9l B per Tb atom and the large spin-orbit coupling associated with 4f orbitals. The exchange coupling results in a ferrimagnetic alignment of the Fe and Tb moments with the Curie point (T C % 700 K) considerably above room temperature. 37 TbFe 2 crystallizes in a cubic C 15 Laves phase structure, where the h111i-direction is the magnetic easy axis, which can be modeled by the cubic crystalline anisotropy constants K 1 ¼ À1.2 Â 10 8 erg/cm 3 and K 2 ¼ 2.08 Â 10 7 erg/cm 3 . 38 The resulting high coercivity is often reduced for application purposes by introducing Dy (h001i easy axis) to obtain the ternary alloy Terfenol-D (Tb x Dy 1Àx Fe 2 ). The desired low coercivity with large magnetostriction that is favorable for magneto-acoustic transduction applications can be tailored by different ratios x, where x ¼ 0.27 is found to be optimal at room temperature. 39 Despite the potentially rich, coupled magnetization, and lattice dynamics in this magnetostrictive ferrimagnet, there have been only a few reports 8,23 that aim at quantifying and correlating the strain evolution and its coupling to the magnetization by time-resolved measurements in binary rare-earth alloys.
Here, we display the different capabilities of table-top ultrafast Xray diffraction (UXRD) and all-optical methods to probe the strain propagation and evolution in a layered magnetostrictive heterostructure. Femtosecond laser pulses are used to excite the rare-earth alloy TbFe 2 , which serves as a transducer for strain waves into adjacent layers. UXRD measurements observe the arrival and shape of the strain waves in a thin, buried detection layer. From this, we extract the stress profile that generates the strain wave in the inhomogeneously excited TbFe 2 layer. The timings of the observed experimental features are rationalized by modeling the strain propagation in this multilayer sample using a 1-dimensional linear chain model of masses and springs. The modeling is shown to be particularly useful when the transducer is capped by a transparent layer so that the conventional symmetric bipolar strain pulse is split into an asymmetric bipolar pulse travelling into the transducer and a unipolar strain pulse that is reflected at the sample-surface after a time determined by the transparent layer thickness. Complementary to the UXRD data, we employ an all-optical polarization sensitive measurement that probes the strain propagation in a transparent silica (SiO 2 ) acoustic delay line. Timeresolved magneto-optical-Kerr-effect (MOKE) measurements are shown to be a very sensitive probe for the arrival of the multiple strain echoes at the top of the laser-excited TbFe 2 layer, which can be used to complement the bulk sensitive UXRD.

II. EXPERIMENTAL DETAILS
We investigate laser-excited samples that consist of (110) oriented Terfenol (TbFe 2 ) layers grown by MBE on (1121) oriented Sapphire (Al 2 O 3 ) with a buried Niobium (Nb) (110) buffer layer as previously described. 38,40 The basic sample structure is only capped by a 2 nm thin protective Titanium (Ti) layer, which does not significantly contribute to the experimental transients. We therefore refer to sample 1 as uncapped. The second sample was instead capped with an 885 nm thick amorphous silica (SiO 2 ) layer. The UXRD measurements are carried out at a laser-driven, plasma-based diffraction setup (PXS) that supplies 200 fs X-ray pulses at Cu K a -energy. 41 The tabletop laser-pump X-ray-probe setup uses p-polarized excitation pulses at a central wavelength of 800 nm, with a 1 kHz repetition rate and a full width at half maximum spot size of a 2-dimensional Gaussian function of 1.4 mm Â 1.5 mm for the laser pulses and 0.3 mm Â 0.3 mm spot size of the X-ray pulses. Using the top-hat approximation with the 1/e width for the laser excitation profile and the pulse energy, we calculate the fluence for the TbFe 2 and Nb experiments, respectively. The pump-fluence at the Nb angle is approximately 4% larger compared to the TbFe 2 experiments since the 1.6 larger diffraction angle leads to a smaller laser footprint, whereas the Fresnel reflection coefficient for the p-polarized laser light decreases by approximately 1.4%.
A representative reciprocal-space map (RSM) of the uncapped sample structure obtained at the PXS alongside the static X-ray diffraction curve and the temporal evolution of the material specific Bragg peaks are displayed in Fig. 1. In the probed RSM volume, we find three separated peaks with their maximum intensity at the out-of-plane reciprocal space coordinate q z ¼ 2.42 Å À1 , 2.64 Å À1 , and 2.69 Å À1 , which are attributed to TbFe 2 (220), Al 2 O 3 (1121), and Nb (110), respectively, according to their bulk lattice plane spacings. The layer thicknesses set by the sample growth are 500 nm TbFe 2 on top of 50 nm Nb and 330 mm Al 2 O 3 as schematically depicted in Fig. 1(e).
The TbFe 2 diffraction peak is significantly broadened along the in-plane reciprocal space coordinate q x compared to the instrument function limited Al 2 O 3 substrate peak. This is a hallmark for microcrystalline domains that in this case exhibit a large mosaic spread of 1.5 around the bulk diffraction angle, 40 very similar to previously reported UXRD experiments on ferroelectric samples. 42,43 The presence of such structural imperfections in the TbFe 2 becomes evident by comparison to the diffraction peak of the Nb layer, which exhibits a much smaller width in q x . The blue and orange lines in Fig. 1(b) represent the reciprocal space slices that are probed in our setup for two fixed angles of incidence (AOI) x that are chosen to be selectively sensitive to the TbFe 2 and Nb lattice strains, respectively.
The combination of an optic that focuses X-rays onto the sample with a convergence of Dx % 0.3 (Montel optic from Layertec) and an X-ray area pixel detector (Dectris PILATUS-100k) allows for swift data acquisition that avoids time-consuming mesh scans of the AOI (x) and the diffraction angle (h). Each pixel of the X-ray area detector is mapped to reciprocal space coordinates q x and q z using the mapping routine described in a previous publication, 42 which is applicable in the thin film regime. Using this fixed angle detection scheme, the flux of 10 6 photons/s incident on the sample is sufficient to probe the evolution of the material specific diffraction peaks and the laser excited heterostructure with subpicosecond time resolution within few hours. The X-ray diffraction curves of the unexcited sample at the Nb and TbFe 2 AOI are indicated by orange and blue solid lines in Fig. 1 to minor modification of the UXRD signals consistent with a slightly increased sound velocity.

III. ANALYSIS AND DISCUSSION OF THE RESULTS A. Signatures from the sample without SiO 2 capping
Before discussing the transient strain of the SiO 2 capped Terfenol (TbFe 2 ) structures, it is instructive to rationalize the signals seen in the UXRD experiment on the uncapped sample 1 that is schematically depicted in Fig. 1(e) for a fluence of 12.7 mJ/cm 2 and 13.3 mJ/cm 2 for the TbFe 2 and Nb, respectively. At first we discuss the strain evolution in the directly excited, approximately 500 nm thick TbFe 2 layer shown in Fig. 2(a). The blue data points show the experimentally obtained strain from Gaussian fits to the diffraction curves. The representative fits and raw time-resolved data from Figs. 1(c) and 1(d) are shown in Fig. S1 of the supplementary material. Beyond 20 ps, we observe an expansion that manifests in a shift of the diffraction peaks to smaller q z . Within the first 20 ps, one observes a transient shift of the majority of the Bragg peak to larger q z , which coincides with the appearance of a shoulder at smaller q z . Between 40 and 80 ps, we detect a pronounced, triangular shaped strain increase and subsequent decrease in addition to an overall rising background.
No background is observed in the strain response of the 50 nm thin, buried Nb layer displayed in Fig. 2(b). The strain in the Nb layer is close to zero up to 4 ns (not shown), except for the very pronounced, nearly symmetric bipolar strain pulse that starts with a contraction at approximately 90 ps, reverses sign at 117 ps, and ceases at approximately 150 ps. The diffraction peak evolution of the thin Nb layer displayed in Fig. 1(c) exhibits a large peak shift that corresponds to a   The interpretation of the UXRD data from the thick TbFe 2 layer is based on the insights into ultrafast lattice response of photoexcited thin films studied by UXRD, which were previously discussed by Schick et al. 44 The laser illumination leads to the excitation of coherent and incoherent phonons, which superimpose in the strain response of the absorbing layer. 2 The strain pulse composed of the coherent excitation of phonons subsequently propagates at the longitudinal acoustic phonon velocity whereas the thermal energy leaves the excited layer by a slower diffusion process.
Our experiment represents the limiting case of an inhomogeneously excited transducer since the TbFe 2 layer thickness is approximately 25 times larger than the 19.9 nm light intensity penetration depth at 800 nm, which we obtain from ellipsometry measurements using a commercial setup and analyzing software (SENTECH), as discussed in Sec. III of the supplementary material. The appearance of a marked shoulder in the TbFe 2 diffraction signal on the lower q z side for the main diffraction in Fig. 1(d) signals the existence of a highly strained surface layer on top of the nearly unperturbed TbFe 2 . Schick et al. 44 have analyzed in detail that the exponential stress profile originating from inhomogeneous laser heating leads to an initial compression of the majority of the layer. Since the Gaussian fit is most sensitive to the central region of the diffraction peak, we observe the leading compressive strain front as a shift of the Bragg peak maximum to larger angles in the first 15 ps. When the free surface expansion propagates into the material at the speed of sound, the strong expansive component finally shifts the Bragg peak maximum to smaller angles. We attribute the remaining slope to the heat transport that equilibrates the inhomogeneous temperature profile within the TbFe 2 layer on a timescale of hundreds of picoseconds to several nanoseconds.
The strain response of the Nb layer seen in Fig. 2(b) confirms that the thermal transport occurs mainly within the TbFe 2 layer since we observe no thermal expansion that would appear as a background within our 4 ns measurement window. The bipolar strain pulse marks the delayed passage of the coherently excited phonon wave packet, which is launched at the sample-air interface, through the buried Nb layer. The detected diffraction peak shift of the 50 nm thick detection layer thus shows a background-free signal of the strain pulse, consisting of a compressive leading edge, which is followed by an expansive trailing edge as it is known from previous picosecond acoustic investigations. 2,45 The smaller layer thickness leads to higher average strain signals and sharper features as compared to the strain detected in the thick transducer layer.
The signature of the exit of the bipolar strain pulse from the probed TbFe 2 layer is an increase in the average layer strain followed by a decrease back to the thermal expansion background since the leading compressive edge exits while the trailing expansive part is still in the layer. In our experiment, we observe a pronounced delay between the exit of the strain wave from the TbFe 2 layer, at approximately 40 ps and its arrival in the adjacent Nb at 90 ps. This 50 ps delay of the signatures can only be rationalized if the strain pulse traverses a TbFe 2 layer that does not contribute significantly to the X-ray diffraction signal. Using v sound ¼ 3.94 nm/ps, known for polycrystalline TbFe 2 , 46 this corresponds to a layer with a thickness of approximately 187 nm TbFe 2 that has a considerably different texture. Structural inhomogeneities are in-line with the mosaic peak broadening and the comparably small X-ray diffraction intensity of the TbFe 2 peak. The existence of a structurally imperfect interface layer at the TbFe 2 -Nb interface is further supported by Atomic Force microscopy measurements that observed that the rare-earth alloy layer growth proceeds first as separated 3-dimensional islands that only coalesce to form a continuous film for thicknesses on the order of 100 nm and above. 38 The large in-plane lattice-constant mismatch of 11.6% between the underlying Nb template and the TbFe 2 is reduced by a thin FeNb layer but is probably the origin for the large mosaicity and for the limited coherence length along the growth direction, 47 which amounts to 50 nm in the present TbFe 2 film.
Modeling the excitation and propagation of picosecond acoustic strain pulses is achieved by solving the partial differential equation for the time-dependent local strain in which the given spatio-temporal stress profile acts as source term. 2,44,48 The thin film geometry reduces this to a 1-dimensional problem as the laser excitation spot is much larger than the film thickness of the nanostructure and the probed X-ray spot. Numerical solutions for the strain evolution in nanoscopic layers are frequently applied in nanoscopic heterostructure geometries where multiple interface reflections complicate analytical solutions. 32,48,49 In Figs. 2(a) and 2(b), we compare the UXRD data to simulation results obtained with the udkm1Dsim toolbox package that we used to calculate the time-resolved strain response based on a linear chain model of masses and springs. Although details of the software are given in the reference publication, 50 we briefly outline the workflow of the modeling. Upon input of the thermophysical material properties and the known sample geometry that are listed in Table I of the supplementary material, we first calculate the absorbed optical energy density and temperature profiles according to the heat diffusion equation 51 with unit cell resolution. The resulting spatio-temporal temperature profile represents the thermoelastic stress that drives a linear chain of masses and springs, where the masses represent individual unit cells. In the last simulation step, the obtained time-resolved strains are used as an input for the computation of the timedependent dynamical X-ray diffraction signal. 52 The resulting diffraction peaks are fitted with a Gaussian line profile to yield the strain signal displayed as lines in Figs. 2(a) and 2(b).
The dashed line in Fig. 2(a) shows the modeled strain of a structurally perfect 436 nm thick TbFe 2 layer on top of a 50 nm Nb layer attached to an Al 2 O 3 substrate. The simulation data represented by the solid line assume only 249 nm structurally perfect TbFe 2 on top of a 187 nm TbFe 2 layer with substantial disorder. The improved fit of the model regarding the triangular feature beginning at 40 ps substantiates the evidence for a structurally different TbFe 2 layer at the Nb interface. The total TbFe 2 layer thickness is determined by the arrival time of the bipolar strain pulse in Nb using the directionally averaged speed of sound of v sound ¼ 3.94 nm/ps for polycrystalline TbFe 2 , 46 due to the lack of exact elastic constants for single-crystalline TbFe 2 . Despite the agreement between the simulated and experimental strains in the TbFe 2 layer, the simulation substantially overestimates the bipolar strain pulse amplitude in the Nb layer. This may be accounted for by taking into account a slowly rising stress profile in TbFe 2 and acoustic damping as well as scattering of the coherent phonons at the interface. [53][54][55] Structural Dynamics ARTICLE scitation.org/journal/sdy

B. Signals in the SiO 2 capped sample
The central experimental result of this work is summarized in Fig. 3. Here, we combine the experimental results from a sample capped by an amorphous SiO 2 layer as sketched in the inset of Fig.  3(e). Figure 3(a) shows the spatio-temporal strain profile that is obtained within the linear chain model for the second sample structure with the stacking sequence 882 nm SiO 2 / 342 nm TbFe 2 / 50 nm Nb / Al 2 O 3 substrate, which is solved by the udkm1Dsim toolbox. 50 Red and blue colors correspond to regions of expansive and compressive strain, respectively. One observes that the expansion of the laserheated region slowly spreads as the heat diffuses within the TbFe 2 layer. Now, the bipolar strain pulse that is launched towards the substrate is clearly asymmetric, where a large-amplitude leading compressive part is followed by a smaller expansive tail. In addition, a unipolar compression pulse propagates in the SiO 2 capping layer towards the surface where it is converted into an expansion. The simulated strain in Fig. 3(a) clarifies that for a very thin capping layer, the expansive unipolar wave reflected at the surface would superimpose with the small expansive tail of the asymmetric bipolar wave to form the symmetric bipolar wave observed in sample 1, see Fig. 2. Multiple reflections of the strain pulses occur at the material interfaces indicated by horizontal dashed lines where the reflection and transmission arise due to the acoustic impedance mismatch. 2,56,57 A direct comparison of the simulation results is presented in Sec. IV of the supplementary material.
The occurrence of multiple unipolar strain pulse echoes that traverse the TbFe 2 and Nb layers at different timings is readily seen in the UXRD data presented in Fig. 3(b). The modeled average strain shown as solid lines accurately predicts the timing and shape of the observed features but the amplitude of the Nb strain is substantially overestimated. This may indicate a finite electron-phonon-coupling time and scattering of the coherent phonons in TbFe 2 from structural imperfections, which are both not captured in the current modeling. The presented UXRD data were obtained under identical excitation conditions as the experiments on the uncapped sample except for the larger pump-fluence 24.6 mJ/cm 2 , which leads to an increase in the detected strain amplitudes.
In the following, we discuss the results of time-resolved MOKE measurements, which probe the change of the polarization state of the probe light upon reflection due to the permanent magnetization of the sample. The measurement displayed in Fig. 3(c) was carried out close to the polar MOKE geometry with an external out-of-plane magnetic field of l 0 H ¼ 800 mT using 200 fs laser pulses at a central wavelength of approximately 800 nm, a repetition rate of 250 kHz, and a pump fluence of approximately 2.7 mJ/cm 2 . In these measurements, the difference of the polarization changes for opposite external field orientations [S(H up ) À S(H down )] is probed using the reflection of 800 nm probe-light-pulses analyzed by a half-wave plate in combination with a Wollaston-prism and a balanced photo-diode. Lock-in detection using an acousto-optical modulation of the pump beam intensity at 50 kHz was employed. The resulting polar MOKE signal displayed as a solid grey line essentially probes the out-of-plane magnetization component of the TbFe 2 layer within the 19.9 nm optical penetration depth and is probably sensitive to the Fe sub-lattice. 58 Subtraction of the slowly varying thermal background approximated by a doubleexponential decay (red line) from the MOKE signal (grey line) reveals multiple sharp peaks in the residual black curve. By comparison with the linear chain model results in Fig. 3(a), it becomes obvious that the observed features occur at the time when the longitudinal strain pulse echoes traverse the top few nanometers of the TbFe 2 layer given by the penetration depth of the probe pulse. Note that the sign of the peaks correlates with the sign of the (unipolar) strain-pulse echoes and that even the small reflections from the TbFe 2 /Nb interface produce observable MOKE signatures at around 180, 485, and 790 ps. Figure 3(d) displays the time-resolved polarization analysis signal of the reflected 800 nm probe beam independent of the magnetization state, which is obtained from the sum signal [S(H up ) þ S(H down )] of the balanced detection. Similar time-resolved Brillouin scattering experiments have shown that the observed oscillations originate from the interference of the reflected light from the traveling strain pulse in the transparent SiO 2 medium and the static interfaces. 59,60 Pronounced phase jumps in this Brillouin signal occur when the strain pulses invert their sign due to the reflection at the SiO 2 /air interface. 56,61 Consequently, the strain propagation as modeled by the 1dimensional-linear chain model accurately predicts the timings of all

ARTICLE
scitation.org/journal/sdy the experimental signals we observed in this heterostructure. We have employed a single temperature model for the driving stress on the lattice. This certainly oversimplifies the equilibration process of the electron-, lattice-, and spin-subsystems to occur instantaneously. A detailed analysis of the dynamics prior to the equilibration is beyond the scope of the current investigation as it requires the knowledge of the thermophysical properties for each sub-system as well as the coupling constants of this largely unexplored material.

C. Experimental results from the buried detection layer
In the henceforth presented data analysis, we put the focus on the qualitative and quantitative information that can be directly extracted from the UXRD signal in the buried Nb detection layer. In Fig. 4(a), we see that the normalized, background-subtracted MOKE signal originating from TbFe 2 matches the normalized Nb strain when shifted by 84 ps, which is the longitudinal acoustic propagation time through the TbFe 2 layer. This agreement proves a linear relation between the lattice strain and the observed MOKE signal. The slight discrepancies at 485 and 790 ps probably originate from the fact that the MOKE signal results from a superposition of the strain pulses reflected at the surface and at the TbFe 2 /Nb interface, which traverse the top TbFe 2 layer simultaneously [see Fig. 3(a)]. Since only part of the reflection at the TbFe 2 /Nb interface is again reflected at the TbFe 2 /SiO 2 interface, the Nb layer senses a different strain wave composition. It will be important for future experiments investigating the interaction of shear waves with the magnetization to accurately identify also the small longitudinal acoustic pulse echoes in such multilayered structures. Previous picosecond acoustic investigations in magnetic samples have observed that strain pulses can exert a torque on the sample magnetizationM via a transient modification of the crystalline anisotropy, often resulting in a damped precessional motion ofM. 18,62,63 Although a torque on the magnetization by the strain pulse is expected, the absence of precessional oscillations challenges the theoretical interpretation of the observed MOKE signal based on the Landau-Lifshitz-Gilbert model. Crystalline defects in the TbFe 2 , resulting in small magnetic domains, 63 in combination with a magneto-crystalline anisotropy and damping could drastically suppress the coherent precessional signal. The signal might have contributions from a modulation of the reflectivity driven by the photoelastic effect; 21 however, the reflectivity signal does not exhibit significant spikes at the echo positions. In any case, the striking resemblance of the detected strain pulses in the Nb layer to extracted features in the MOKE measurements demonstrates a high sensitivity of MOKE for probing strain pulses arriving at the TbFe 2 surface.
In Fig. 4(b), we compare the initial bipolar strain pulses from the capped and uncapped samples, normalized to their compressive part. We observe that the leading, compressive parts coincide, whereas the expansive parts in the SiO 2 capped sample 2 are strongly reduced. The black dashed line indicates a single exponential fit to the falling edge of the compressive strain with a time constant of 4.6 ps, which translates to a spatial extension of approximately 18 nm using v sound ¼ 3.94 nm/ps. This value provides an estimation of the spatial extension of the driving stress profile. 3,64 This matches the optical penetration depth obtained from ellipsometry, which shows that potential hot-electron diffusion does not substantially increase the excitation profile length.
The evolution of the unipolar strain pulse within the SiO 2 capping is shown in Fig. 4(c), where we compare the first and second unipolar strain echoes that traverse the Nb layer at 396 ps and 695 ps to the initial pulse launched into the SiO 2 capping. The latter is extracted from the difference between the bipolar strain pulses observed on the capped and uncapped samples. In this analysis, we assume that the laser generated stress profile is identical in the capped and uncapped TbFe 2 samples and that the difference of the initially detected bipolar strain pulse seen in Fig. 4(b) originates solely from partial reflection of the expansion at the top TbFe 2 interface. The FWHM of the detected strain signals increases from 9 to 16 and 21 ps. The modeling does not include any broadening mechanisms and reports echoes with a constant width of 8 ps FWHM. Anharmonic interactions in the lattice potential have been shown to change the shape and broaden high amplitude coherent phonon wavepackets. 65 Contributions from the SiO 2 surface roughness should also be taken into account.

IV. CONCLUSION
In this work, we have combined multiple techniques to follow the trajectory of strain pulses that are generated by femtosecond laser pulses exploiting the giant magnetostriction material Terfenol (TbFe 2 ) as a transducer. MOKE measurements in TbFe 2 are shown to provide a surface sensitive method to probe strain pulses at the top of the metallic TbFe 2 that is complementary to bulk sensitive X-ray diffraction.
Probing the strain pulse in a thin, buried detection layer adjacent to an optically opaque transducer via UXRD allows for a characterization of the coherent strain pulse separately from heat expansion without frontside-pump backside-probe schemes. By detecting the strain pulse in the buried and perfect Nb layer, we demonstrate how UXRD can clearly locate the structurally imperfect fraction of the TbFe 2 layer in this opaque heterostructure. Contrary to all-optical methods, UXRD provides a quantitative measure of the average lattice strain that does not require detailed knowledge of photo elastic coefficients and optical properties. In combination with 1dimensional-linear chain models, UXRD can provide quantitative information on realistic strains with unit cell resolution, which can be

ARTICLE
scitation.org/journal/sdy used as input for modeling strain-assisted magnetization switching approaches. 23 The obtained maximum strain amplitude in TbFe 2 of 2& is well below the deterministic switching limit of a few percent strain, but the local heating substantially lowers the anisotropy 66 as it is known from heat-assisted magnetic recording schemes. The combination of MOKE and UXRD outlines a potential path towards an experimental calibration of the magneto-elastic coefficient, i.e., the magnetization change per strain amplitude for picosecond acoustic pulses. Such a quantity is not only relevant for testing fundamental research that models magneto-elastic couplings but also represents a valuable input for application-oriented research.
We believe that probing the strain-pulse in a buried detection layer is a versatile method for studying the stress generation profile as it separates coherent from incoherent phonon excitations. It will be especially useful in situations where multiple mechanisms with different spatial or temporal characteristics superimpose in the strain generation process as it is the case in (anti-)ferromagnetic 67,68 materials. The use of an acoustic delay line further introduces the possibility to study the evolution of the strain pulse shape and to calibrate the magnetization response to unipolar compression and expansion pulses. This will support important future steps towards a full understanding of the demagnetization process especially in high-anisotropy, giant magnetostriction materials. A combination of time-resolved probes that monitor different degrees of freedom within the same experiment will foster the understanding of the intricate couplings between electron-, spin-, and lattice systems in solids, which forms the basis for many useful devices.

SUPPLEMENTARY MATERIAL
See supplementary material for the details of the time evolution of the TbFe 2 and Nb Bragg peaks (S1) and the material parameters used in the modeling (S2). Furthermore, we provide the complex index of refraction and the resulting optical properties that were extracted from spectroscopic ellipsometry on the uncapped TbFe 2 sample (S3) as well as a section that compares the modeling results for the strain propagation in the capped and uncapped samples (S4).