Ultrafast Opto-magnetic Effects in the Extreme Ultraviolet Spectral Range

Coherent light-matter interactions mediated by opto-magnetic phenomena like the inverse Faraday effect (IFE) are expected to provide a non-thermal pathway for ultrafast manipulation of magnetism on timescales as short as the excitation pulse itself. As the IFE scales with the spin-orbit coupling strength of the involved electronic states, photo-exciting the strongly spin-orbit coupled core-level electrons in magnetic materials appears as an appealing method to transiently generate large opto-magnetic moments. Here, we investigate this scenario in a ferrimagnetic GdFeCo alloy by using intense and circularly polarized pulses of extreme ultraviolet radiation. Our results reveal ultrafast and strong helicity-dependent magnetic effects which are in line with the characteristic fingerprints of an IFE, corroborated by ab initio opto-magnetic IFE theory and atomistic spin dynamics simulations.


Introduction
Light-driven ultrafast magnetic phenomena occurring on ever faster and smaller time and length scales are at the core of the modern magnetism research 1,2 .A long sought-after and yet to be realized phenomenon is the coherent and deterministic control of a macroscopically ordered spin ensemble on the sub-cycle timescales of the photo-exciting laser field 3,4 .In this context, one of the most prominent effects that could lead to large, ultrafast changes of magnetization is the inverse Faraday effect (IFE).
Discovered originally in paramagnetic solids 5 , the IFE was shown to generate an effective magnetic field or, equivalently, an induced magnetic moment in a medium upon photoexcitation with circularly polarized light pulses.According to classical IFE theory 5 , the orientation and the magnitude of the lightinduced magnetic moment are respectively determined by the helicity of circularly polarized light (or its angular momentum) and the combined effect of the magneto-optical susceptibility of the photoexcited medium and the incident light intensity.In case of metallic ferro-or antiferromagnetic materials, a fully ab initio quantum-mechanical treatment of the IFE has revealed that its action on these materials is significantly different compared to the originally discovered effect, predicting a complex dependence of the IFE-induced magnetization on both helicity and wavelength, which is not covered anymore by the conventional theory [6][7][8] .However, the high absorption of the optical excitation in these materials typically leads to strong heat-induced demagnetization which hinders the observation of any potential non-thermally induced magnetization 9 .Thus, the actual influence of an IFE on the sub-picosecond spin dynamics observed in such systems as, for example, all-optical switching phenomena in ferrimagnetic rare-earth transition-metal alloys like GdFeCo [10][11][12][13] , remains highly debated and elusive so far 14 .
The magnitude of the IFE depends on the light intensity and its polarization, as well as on the spin-orbit coupling (SOC) of the involved electronic states via the wavelength-and helicity-dependent optomagnetic constant IFE ()    6,8 : where () E  denotes the light electric field, σ± its helicity, and ind M the induced magnetization along the k vector of the incident circularly polarized light beam.The IFE, which is an electronic Raman process, has thus far been investigated theoretically 6,8,[15][16][17][18] , but a direct comparison with measurements is lacking.A potential way of generating strong opto-magnetic interactions in metallic magnets is the excitation of core-level electrons, involving states with much higher SOC compared to the valence band (e.g., 1.1 eV vs. 65 meV in case of Fe 3p and 3d electrons, respectively 19 ) that can be accessed by short wavelength radiation in the extreme ultraviolet (XUV) or soft x-ray range.The availability of freeelectron lasers (FEL) providing femtosecond pulses of very high brilliance with control over the light polarization has enabled studies of x-ray induced magnetization dynamics, extending the capabilities from element-specific probing, employing, e.g., the x-ray magnetic circular dichroism (XMCD), towards element-specific excitation of highly non-equilibrium states [20][21][22][23][24] .However, besides first theoretical predictions, there is, to the best of our knowledge, no experimental evidence available so far on the existence of opto-magnetic effects like the IFE at these XUV wavelengths.
In this work, we study the ultrafast magnetization dynamics of ferrimagnetic GdFeCo alloy induced by photoexcitation with intense femtosecond circularly polarized pulses in the XUV spectral range generated by the free-electron laser FERMI.Comparing the helicity-dependent dynamics upon on-and off-resonant excitation with respect to the Fe M3,2 resonance, we find strong dynamic helicity-and wavelength-dependent effects that resemble the expected characteristics of an IFE.Analyzing the magnitude and wavelength-dependence of the observed effect, we can rule out a thermal origin caused by dichroic absorption due to the XMCD.This finding is corroborated by atomistic spin dynamics (ASD) simulations, revealing that even in case of fully on-resonant excitation, the expected XMCD response is up to an order of magnitude smaller compared to the observed effect.Instead, we find qualitative agreement with ab initio calculations of the IFE in the XUV spectral range.Our results thus identify the existence of an IFE in the XUV spectral range.

Time-resolved pump -probe measurements
The experimental concept employing the XUV pump − magneto-optics probe technique to trigger and measure the IFE is depicted in Fig. 1a.The sample was excited at normal incidence using ≈90 fs (full width at half maximum, FWHM) XUV pulses of variable polarization [circular right/left (σ±) and linear horizontal (lin.hor.)], tuned to photon energies around the Fe M3,2 resonance.The XUV-induced dynamics were probed under 45° by linearly polarized ≈90 fs (FWHM) optical pulses at 400 nm wavelength, measuring both the magneto-optical Faraday rotation in transmission and the Kerr rotation in reflection, obtaining magnetic contrast by flipping an out-of-plane saturating magnetic field of ±8 mT applied to the sample.Fig. 2 shows the temporal evolution of magnetization after resonant excitation by linearly and circularly polarized XUV pulses at a photon energy of 54.1 eV, as measured by the Faraday signal; for the Kerr measurements leading to similar results, see Supplementary Note 3. The data was fitted using a doubleexponential fit function convoluted with a Gaussian (solid lines in Fig. 2), taking into account the ultrafast demagnetization, the subsequent relaxation process and the experimental time resolution of ≈280 fs.For clarity, the pump-probe scans are displayed only for two incident XUV fluences, showing the dynamics upon low and high fluence excitation.The data reveals a clear influence of the XUV polarization on the demagnetization amplitudes, emerging upon excitation with increasing incident fluence (compare 1.3 and 6.7 mJ/cm 2 ).Above 1.3 mJ/cm 2 , the XUV pump helicity drives either an enhancement or a decrease of the demagnetization effect that takes place within the first hundreds of femtoseconds; this helicity-dependent difference in demagnetization magnitudes lasts for tens of picoseconds, i.e., also during the magnetization relaxation process, emphasizing the strength and robustness of the transient XUV-induced magnetization.Such a helicity-dependent behavior is at variance with the case of linearly polarized XUV excitation that generates demagnetization based on purely thermal effects, i.e., demagnetization driven by ultrafast electronic heating.

Fluence-and wavelength-dependence
In Fig. 3 3b) show that the effect is wavelength-dependent, clearly changing its polarity along the resonance.The data further reveals that the effect is strongly fluence-dependent; within the experimental uncertainty, two fluence regimes become apparent (illustrated by the white/grey shaded background in Fig. 3b) where the magnitude of ∆M initially increases with XUV fluence and then either reaches a saturated state at a constant ∆M value or decreases again as the fluence is further elevated.In case of below-and above-resonant excitation (51.0 and 64.0 eV), ∆M scales linearly with the fluence until the saturated state is achieved.It thereby reaches up to (32.2±0.7)% difference in pump-induced change with respect to the equilibrium magnetization (largest ∆M observed for 5.3 mJ/cm 2 at 64.0 eV).In case of resonant excitation (54.1 and 56.1 eV), the data reveals a more complex fluence-dependence, indicating either a threshold-like behavior (54.1 eV) or even a fluence-dependent sign change of ∆M (56.1 eV).For excitation between 54.1 and 64.0 eV, the demagnetization reaches amplitudes of up to ≈100% at high fluences (depending on the helicity), corresponding to full demagnetization of the magnetic layer.This leads to the observed saturation or even quenching of ∆M in the high fluence regime, as the increasing thermal load and demagnetization starts to counteract the helicity-dependent effect (grey shaded regions).
The most striking observation is the presence of strong helicity-dependent differences even for the offresonant excitation at photon energies where the XUV spectroscopy shows very small or even zero XMCD (compare Fig. 1b for 51.0 and 64.0 eV).In general, the XMCD provides a potential source of helicity-dependent dynamics due to the interaction of circularly polarized XUV pulses and the angular momentum of the excited electrons leading to a polarization-dependence of the absorbed fluence.This dichroic absorption effect would manifest in a different amount of thermally induced demagnetization, depending on the helicity of the XUV excitation pulse.However, due to the lack of sizeable dichroic absorption in case of 51.0 and 64.0 eV (compare Fig. 1b), this mechanism can already be ruled out as the driving force behind the observed helicity-dependence at these two photon energies, which strongly indicates the existence of an opto-magnetic effect such as the IFE.This interpretation is further corroborated by the linear scaling with the fluence (XUV intensity) observed at 51.0 and 64.0 eV (white shaded regions in Fig. 3b), resembling a typical IFE feature, since the magnetization induced by the IFE depends linearly on the incident light intensity via the opto-magnetic, wavelength-dependent constant (see Eq. 1 and Battiato et al. 8 ).A possible explanation for the more complex fluence-dependence and the deviation from a purely linear scaling observed at 54.1 and 56.1 eV can be given by the steep slope of the XAS and XMCD spectra (see Fig. 1b) at these photon energies, significantly changing the response of the system when the spectrum is slightly reshaped or shifted while the system is driven out of equilibrium, depending on the strength of the excitation 25,26 .
The generation of helicity-dependent effects that are longer-lived compared to the duration of the driving pump pulse (see Fig. 2) have been previously reported in various experimental studies of the IFE on different materials 9,13,15,[27][28][29] .Although the IFE is active only during the pump pulse duration, the subsequent magnetization evolution is affected on much longer time scales, especially, when the impulsively IFE-induced magnetization is large (here we reach a maximum of ≈30% of the ground-state M0 values) and, consequently, the system needs longer times to accommodate such large magnetic changes.Particularly for high fluence excitation, the helicity-dependent effects therefore persist until the end of the pump-probe delay range (-2.2 to 30 ps) scanned in the experiment (for a full set of time traces as a function of fluence, see Supplementary Figure 4 and Supplementary Note 2).

Ab initio calculations of the inverse Faraday effect
To compare our experimental results to theoretical predictions of an opto-magnetic effect in the XUV spectral range, we carried out first ab initio calculations of the IFE 6,8 in vicinity of the Fe M3,2 resonance.
To reduce the computational complexity, the system was modeled as a stoichiometric GdFe2Co unit cell, focusing on the static opto-magnetic response for left (σ−) and right (σ+) circularly polarized excitation.
The difference between the calculated opto-magnetic constants is shown in Fig. 4 (right axis) in comparison to the largest helicity-dependent effects that were experimentally observed for the studied excitation photon energies (left axis, see also Fig. 3b).Note that this comparison is purely qualitative, as it is based on two non-equivalent quantities, namely the difference between the static opto-magnetic constants vs. the transient demagnetization amplitudes obtained at a certain time after excitation.A fully quantitative comparison cannot be made yet, as the constants IFE ()    are obtained from stationary IFE calculations, so far not taking into account the fluence-dependent, nonequilibrium state present upon femtosecond excitation.To obtain a qualitative agreement with the experimental data, the calculated spectrum of the IFE response had to be shifted by +2.5 eV, which can be attributed to the effect of the core holes that is not included in the calculations [30][31][32] (see also the Methods section).The calculations support our experimental findings regarding the existence of an optomagnetic effect at these wavelengths by predicting a highly wavelength-dependent IFE in vicinity of the M3,2 resonance.The calculated ∆IFE reverses its sign going through the resonance, which is again in qualitative agreement to the different polarity of ∆M observed in the experiment for below-and aboveresonant excitation.This sign reversal occurs because the IFE spectra As mentioned before, even a small but ultrafast spectral shift of the Fe M3,2 resonance, which has been experimentally observed 25,26 for laser-driven non-equilibrium states in similar magnetic materials, can lead to a sign reversal of the IFE depending on the excitation fluence.Furthermore, it has to be emphasized that the calculations show a finite response at the photon energies 51.0 and 64.0 eV, corroborating the existence of a helicity-dependent interaction at photon energies where a thermal mechanism due to the XMCD effect could already be ruled out.However, the calculations are so far not able to quantitatively explain the experimental results regarding the magnitude of the observed helicity-

Atomistic spin dynamics simulations
In case of fully resonant excitation at 54.1 and 56.1 eV, the observed helicity-dependence could further result from the combined action of an IFE together with the dichroic absorption caused by the finite XMCD present at these photon energies (compare Fig. 1b).To quantify the actually expected effect of the XMCD on the helicity-dependent dynamics, atomistic spin dynamics (ASD) simulations were carried out 11,[33][34][35] , simulating the influence of the dichroic absorption on the helicity-dependent demagnetization amplitudes (see Fig. 4).For best comparability, the excitation fluence was calibrated between experiment and simulation by matching the demagnetization amplitudes obtained upon linearly polarized excitation.While the absence of any helicity-dependence in the XMCD-driven ASD simulations is obvious for photon energies without XMCD (51.0 and 64.0 eV), the simulation shows that even in the fully resonant case (54.1 and 56.1 eV), the expected dichroic demagnetization effects would be smaller by at least a factor of six (comparing ∆Msim and ∆Mexp at 54.1 eV) and by almost an order of magnitude compared to the largest effect observed across the spectrum.Thus, a significant contribution of a purely thermal mechanism (via dichroic absorption) on the observed helicitydependent dynamics can be ruled out, strengthening our conclusions regarding the existence of an optomagnetic effect like the IFE for both on-and off-resonant excitation in vicinity of the Fe M3,2 resonance.

Conclusions
In conclusion, we studied the impact of femtosecond XUV photoexcitation of a ferrimagnetic GdFeCo x-rays 20,36,37 .The observation of large opto-magnetic effects in the XUV spectral range also implies a significant, yet to be explored impact on ultrafast magnetization switching phenomena in addition to the MCD-induced helicity-dependence known from the visible-light regime 38 (regarding an IFE in the visible-light regime, see also Supplementary Note 4).Moreover, the different nature of the Gd 4f and Fe 3d magnetic moments with respect to their SOC and predominantly localized vs. itinerant character suggests that comparing element-specific excitation of the rare-earth and transition metal sublattices will contribute to fully revealing the microscopic mechanisms behind the XUV and soft x-ray IFE.
From a more general perspective, our findings reveal an efficient method to transiently generate large macroscopic magnetization on ultrafast time scales.Microscopically rooted in the stimulated Raman scattering process, the XUV-induced IFE demonstrated here can be seen, depending on the helicity of the exciting XUV pulse, to either generate or annihilate a transient magnetization within a range of up to ≈30% of the ground-state magnetization value of ferrimagnetic GdFeCo.Given its strength and robustness, such an opto-magnetic core-level effect is expected to be of relevance for the fields of ultrafast magnetism and spintronics, coherent magnetization control and nonlinear x-ray science.
Moreover, touching upon the generality of the effect, the IFE via core-level excitations could essentially be extended to any nonmagnetic materials with strong spin-orbit coupling in order to transiently generate a macroscopic magnetic moment.

Experimental techniques and data acquisition
The time-resolved studies were performed at the DiProI end-station of the FERMI FEL 39 , using a pumpprobe geometry as schematically shown in Fig. 1a.For simultaneously measuring both the magnetooptical Faraday rotation in transmission and the Kerr rotation in reflection, two polarization-sensitive balanced photodetection schemes were employed.The magnetic contrast was obtained upon flipping an out-of-plane magnetic field of ±8 mT that was magnetically saturating the sample, also restoring the initial magnetization state after each pump-probe cycle.Taking the difference between the two traces recorded for opposite magnetization states leads to a quantity that is directly proportional to the magnetization, as the non-magnetic background cancels out.The experimental time resolution was evaluated to be ≈280 fs, using the method described by Ziolek et al. 40 , accounting for the enclosed angle of 45° between pump and probe beams as well as their spatial footprint on the sample.The FEL and optical laser system were operating at 50 Hz repetition rate.The effective repetition rate of the XUV pump was reduced to 25 Hz by seeding only every second electron bunch in the FEL, allowing for an interleaved measurement to record both pumped and unpumped states of the sample.
Supplementary Figure 1 exemplarily shows the data as recorded, i.e., the individually acquired magnetooptical Faraday signals as a function of pump-probe delay.Sorting the data after FEL pulse energy, which was shot-by-shot recorded by an I0 gas monitor detector, allows retrieving the fluence-dependence (compare Supplementary Figure 2).For a more extensive technical background about the experimental methods and data acquisition, see Supplementary Notes 1 and 2.

Sample design and spectroscopy
The studied sample was a 20 nm thin amorphous film of a ferrimagnetic Gd24Fe67Co9 alloy with an outof-plane magnetic anisotropy.For static absorption spectroscopy, it was deposited on a 30 nm Si3N4 membrane, allowing transmission of the XUV radiation.For the time-resolved measurements, a 1 mm thick glass substrate was used in order to allow simultaneous detection of the Faraday and Kerr signals.
The static XUV absorption and XMCD of the sample shown in Fig. 1b was characterized employing the ALICE reflectometer 41 at the PM3 beamline of BESSY II 42 .The XMCD spectrum thereby corresponds to the asymmetry (A) of the transmitted circularly polarized XUV radiation (T↑,↓) measured for opposite magnetic field directions perpendicular to the sample plane, calculated by A = (T↑ − T↓) / (T↑ + T↓).For obtaining the XAS (B↑,↓), the transmitted XUV intensity was normalized to the beamline spectrum (T0) by B↑,↓ = −ln(T↑,↓ / T0).For comparison, this characterization was also carried out for a similar Gd24Fe76 sample, i.e., without the Co ingredient but otherwise identical composition.The comparison is shown in Supplementary Figure 8, showing the influence of the small fraction of Co atoms on the XAS and XMCD.To correct for the limited degree of circular polarization approaching the lower photon energy limits of the beamline, the magnetic part ∆β of the absorptive refractive index was calculated from the measured XMCD for the effective thickness of the Fe content within the alloy and normalized to reference XMCD measurements on pure Fe systems 43 .The corrected XMCD magnitude at the Fe M3,2 resonance agrees with the values reported in the literature for comparable GdFe systems 44,45 .

Atomistic spin dynamics simulations
Atomistic spin dynamics (ASD) simulations were carried out 11 to simulate the influence of the dichroic absorption on the helicity-dependent demagnetization amplitudes.The magnitude of the dichroic absorption was set according to the XMCD asymmetry obtained from the static XUV and XMCD spectroscopy (see Fig. 1b), i.e., (−6.2±0.5)% in case of 54.1 eV, (−3.0±0.5)% in case of 56.1 eV and (0.5±0.5) % in case of 51.0 eV excitation, respectively (note that for 64.0 eV, the XMCD is fully zero, thus no dichroic absorption is expected).The resulting simulated fluence-dependencies of the maximum demagnetization amplitudes upon σ±-polarized excitation are shown in Supplementary Figure 9 together with their difference ∆M.
Simulating also the dynamics excited by linearly polarized radiation allows better comparability The ASD simulations are done using the VAMPIRE software package 33,34 based on the Landau-Lifshitz-Gilbert (LLG) equation.The system Hamiltonian is given by: where Si, Sj are the normalized spin vectors on i, j sites, Jij is the exchange constant and ku is the uniaxial magnetocrystalline anisotropy energy constant per site.We used similar parametrization as reported in literature 12,35 : µFe = 1.92µB and µGd = 7.63µB, the anisotropy is taken as ku = 8.07246 • 10 -24 J and the exchange interactions are JFe-Fe = 2.835 • 10 -21 J, JGd-Gd = 1.26 • 10 -21 J and JFe-Gd = −1.09• 10 -21 J.We incorporate the rapid change in thermal energy of a system under the influence of a femtosecond laser pulse using the two-temperature model 46 : where Ce = 225 Jm -3 K -2 , Cph = 3.1 • 10 6 Jm -3 K -1 and Ge-ph = 2.5 • 10 17 Wm -3 K -1 .The spin system is coupled to the electron temperature Te and the quantity Pabs(t) corresponds to the laser power density absorbed by the electronic system and depends on laser fluence and light polarization.The XMCD effect is included by modifying Pabs(t) proportional to the percentage obtained from the experiments.

Ab initio calculations of the inverse Faraday effect
The ab initio calculations of the IFE were carried out using the second-order response theory formalism as derived by Berritta et al. 6 and Battiato et al. 8 .The calculations were performed with the full-potential, all-electron code WIEN2k 47 , with spin-orbit interaction included.
The system was modeled as a stoichiometric GdFe2Co unit cell in the AuCu3 structure, i.e., as GdFe3, but with one Fe atom replaced by Co.The lattice parameters were first optimized using the VASP code 48 , after which the electronic structure and IFE was computed with the WIEN2k program.To capture the strong electron correlations in the Gd 4f shell, we employed the GGA+U method, in the atomic limit version, with parameters U = 7 eV and J = 1 eV, in combination with the generalized gradient approximation (GGA) 49 to the exchange-correlation potential.The calculated band structure is shown in Supplementary Figure 11.The occupied Gd 4f states are located at −6 eV and the unoccupied 4f states start at +4 eV.The spin-polarized Fe and Co 3d states that are accessed by XUV photon excitation of the 3p states are in the energy window of 0 to 4 eV.
To compute the IFE, the exchange splitting of the 3p semi-core states of Fe and Co atoms was included as well 50 .The computed opto-magnetic IFE response is shown in Supplementary Figure 12.The photon energy-and helicity-dependent IFE constants are given for the total optically induced magnetization Mind, i.e., they contain both spin and orbital contributions (see Berritta et al. 6 ).The IFE constants for σ− and σ+ circular polarizations are clearly distinct in the photon energy range of the Fe and Co M absorption edges while for energies outside of this range the distinction between σ− and σ+ polarizations vanishes.The bottom panel of Supplementary Figure 12 shows ∆IFE, i.e., the difference between the opto-magnetic constants for the two helicities.To understand this spectrum, we can first note that, in the absence of spin magnetization in the alloy, the IFE spectra for σ− and σ+ polarizations would be identical but with opposite sign.This feature can be approximately recognized in Supplementary Figure 12.When there is spin magnetization present, the exchange splitting of the 3d valence bands leads to an energydependent shift in the individual IFE spectra, such that where the 3p→3d resonance occurs for Fe in the calculations.Around 58−60 eV, a similar difference is expected to appear near the Co 3p absorption edge.However, there will be an overlap of the M edges of Fe and Co in this photon energy range which can partially cancel out each other.In addition, we mention that previously, it was found that the ab initio computed XUV spectra of 3d transition metals compared well with measured spectra, but the computed energy onset of the M3,2 edge differs from the measured onset position by about 2.5 eV [30][31][32] .This difference occurs because the ab initio calculated energy positions of the semi-core levels can deviate from the real energy positions and, in addition, there is an effect of the core hole in the XUV excitation which can lead to a shift in the binding energy of the semicore levels 30,50 .The ∆IFE spectrum shown in Fig. 4 has therefore been shifted by +2.5 eV to align with the experimental onset of the Fe M edge, which precisely corresponds to the shift that has been used in previous work 32 .We remark further that the computed IFE spectrum is valid in the quasi-static approximation, i.e., assuming that the computed IFE response is modulated only by the envelope of the XUV pulse.The employed second-order response theory is furthermore expected to be valid for reasonable, but not too high laser fluences.

Supplementary Figure 2.
Analysis of the shot-resolved FEL pulse energies recorded by the I0 gas monitor detector.The graph shows exemplary the data recorded for linearly polarized XUV radiation at a photon energy of 54.1 eV.(a) FEL pulse energies per shot as recorded for an average pulse energy of 5.0 µJ.Only every twentieth shot is plotted for better visibility.(b) Histograms of the statistical distribution of XUV pulse energies recorded during three subsequent pump-probe delay scans using different average excitation fluences.The Faraday/Kerr data was sorted by averaging only over those data points where the sample was excited by FEL shots with the same pulse energy, as defined by a 1 µJ grid (shown as dashed lines). .m /cm Supplementary Figure 4.
Difference between the magnetization transients obtained for σ− and σ+ excitation (data from Supplementary Figure 3), showing the helicity-dependent effect as a function of pump-probe delay and excitation fluence.The error bars are calculated from the error propagation of the standard error of the mean. .m /cm

Fig. 1b shows
Fig.1bshows the XUV absorption spectra (XAS) measured on the GdFeCo sample across the Fe M3,2 , the demagnetization amplitudes obtained from the double-exponential fits are shown as a function of XUV photon energy, polarization, and incident excitation fluence, evaluated for the full range of used XUV fluences.The data reveals a strong helicity-dependent effect that is present over the whole range of photon energies used for excitation, be it on-(54.1 and 56.1 eV) or off-resonant (51.0 and 64.0 eV) with respect to the Fe M3,2 resonance.The splitting between either enhanced or attenuated demagnetization amplitudes upon σ±-polarized excitation (D = 1−min[M/M0], Fig. 3a) and accordingly the magnitude and sign of the difference (∆M = D[σ−]−D[σ+], Fig.
with respect to each other at the 3p→3d resonance, generated by the exchange splitting of the 3d states.The very steep slope of the ∆IFE spectrum around the sign reversal thereby coincides with the fluence-dependent sign change of the helicity-dependent effect ∆M observed for 56.1 eV excitation.
dependence.Although the XAS show only a minor contribution due to the small amount of Co in the alloy (compare Fig. 1b at ≈62.0 eV), the IFE calculations show a finite extension of the opto-magnetic response up to the Co M3,2 resonance, which can qualitatively explain the observation of helicitydependent effects at 64.0 eV, but not yet the relatively large magnitude of the effect compared to the Fe response.
photon energy and polarization reveals a strong and long-lived helicity-dependent effect on the ultrafast between experiment and simulation by comparing the simulated XMCD-induced and experimentally obtained helicity-dependence ∆M based on the respective fluences that lead to the same amount of demagnetization in the linearly polarized case.Both the experimentally observed (yellow color) and simulated (purple color) ∆M values are shown in Supplementary Figure 10 as a function of the maximum demagnetization amplitudes induced by linearly polarized excitation of the same fluence.The comparison between experiment and ASD simulations presented in Fig. 4 is based on the largest ∆M observed in the experiment for each XUV photon energy and the interpolated value of ∆M taken from the ASD simulation resembling the respective magnitude of XMCD.

3 .
Transient magnetization dynamics of the system induced by XUV pulses of 54.1 eV photon energy, probed by the normalized magneto-optical Faraday rotation as a function of pump-probe delay, XUV polarization and incident excitation fluence.The magnetization is normalized to the equilibrium magnetization in the unexcited state (M/M0).The error bars are calculated as the standard error of the mean.