Efficient spin excitation via ultrafast damping-like torques in antiferromagnets

TL;DR

This paper reveals that ultrafast damping significantly influences spin dynamics in antiferromagnets, enabling highly efficient and rapid control of magnetic states for spintronic applications.

Contribution

It demonstrates that ultrafast damping effects are crucial in antiferromagnets, leading to large-amplitude magnetic modulations during impulsive excitations, unlike in ferromagnets.

Findings

Ultrafast damping causes immediate spin canting in antiferromagnets.

Exchange and anisotropy interactions amplify spin canting effects.

This mechanism enables ultrafast manipulation of antiferromagnetic order.

Abstract

Damping effects form the core of many emerging concepts for high-speed spintronic applications. Important characteristics such as device switching times and magnetic domain-wall velocities depend critically on the damping rate. While the implications of spin damping for relaxation processes are intensively studied, damping effects during impulsive spin excitations are assumed to be negligible because of the shortness of the excitation process. Herein, we show that, unlike in ferromagnets, ultrafast damping plays a crucial role in antiferromagnets because of their strongly elliptical spin precession. In time-resolved measurements, we find that ultrafast damping results in an immediate spin canting along the short precession axis. The interplay between antiferromagnetic exchange and magnetic anisotropy amplifies this canting by several orders of magnitude towards large-amplitude modulations of the antiferromagnetic order parameter. This leverage effect discloses a highly efficient route towards the ultrafast manipulation of magnetism in antiferromagnetic spintronics.

Figures

Tables

Table 1. Table 1: Fit parameters. Values for λ = J N ​ ( g ​ μ B ) 2 𝜆 𝐽 𝑁 superscript 𝑔 subscript 𝜇 B 2 \lambda=\frac{J}{N\left(g\mu_{\mathrm{B}}\right)^{2}} and D z = K z N ​ ( g ​ μ B ) 2 subscript 𝐷 𝑧 subscript 𝐾 𝑧 𝑁 superscript 𝑔 subscript 𝜇 B 2 D_{z}=\frac{K_{z}}{N\left(g\mu_{\mathrm{B}}\right)^{2}} were calculated from the exchange interaction J 𝐽 J and the out-of-plane anisotropy constant K z subscript 𝐾 𝑧 K_{z} given in the respective references.

Parameter	Unit	HoMnO3	YMnO3
$T_N$	[$𝖪$]	69.2(7)	74.7(1)
$T_{\mathrm{SR}}$	[$𝖪$]	33.5(2)	—
$\lambda$	[${𝖳}^{\mathrm{𝟤}}{\mathrm{𝖼𝗆}}^{\mathrm{𝟥}}{𝖩}^{-\mathrm{𝟣}}$]	70.845	71.125
$D_z$	[${𝖳}^{\mathrm{𝟤}}{\mathrm{𝖼𝗆}}^{\mathrm{𝟥}}{𝖩}^{-\mathrm{𝟣}}$]	11.045	13.925
$D_y\left(T\right)$	[${𝖳}^{\mathrm{𝟤}}{\mathrm{𝖼𝗆}}^{\mathrm{𝟥}}{𝖩}^{-\mathrm{𝟣}}$]	$\kappa \left(T_{\mathrm{SR}}-T\right)$	0.0251(1)
$\kappa$	[${𝖳}^{\mathrm{𝟤}}{\mathrm{𝖼𝗆}}^{\mathrm{𝟥}}{𝖩}^{-\mathrm{𝟣}}{𝖪}^{-\mathrm{𝟣}}$]	0.0093(6)	0
$\Delta \cdot {\left(g{\mu }_B\right)}^2$	[${𝖳}^{\mathrm{𝟤}}{\mathrm{𝖼𝗆}}^{\mathrm{𝟥}}{𝖩}^{-\mathrm{𝟣}}$]	0.0354(9)	0
$𝒥$	[${𝖳}^{\mathrm{𝟤}}{\mathrm{𝖼𝗆}}^{\mathrm{𝟥}}{𝖩}^{-\mathrm{𝟣}}$]	117.2$+\left|D_y\right|$	120.6
$𝒟$	[${𝖳}^{\mathrm{𝟤}}{\mathrm{𝖼𝗆}}^{\mathrm{𝟥}}{𝖩}^{-\mathrm{𝟣}}$]	$\Delta M_0^2+\left|D_y\right|$	0.0251(1)
$A=\sqrt{𝒥/𝒟}$	—	Fig. 3e	69.3

Equations

$$\frac{d{\mathbf{M}}_i}{dt}=-\gamma {\mu }_0{\mathbf{M}}_i\times {\mathbf{H}}_i-\gamma {\mu }_0\frac{\alpha }{M_0}{\mathbf{M}}_i\times \left({\mathbf{M}}_i\times {\mathbf{H}}_i\right),$$

$$\mathcal{F}=\lambda \sum _{⟨i,j⟩}{\mathbf{M}}_i\cdot {\mathbf{M}}_j+D_z\sum _i{M}_{i,z}^2+D_y\sum _i{M}_{i,y}^2+\Delta \sum _i{M}_{i,x}^2{M}_{i,y}^2,$$

$$\mathbf{M}_{i}(+0)=M_{0}\left(\begin{array}[]{c}1\\ \gamma\mu_{0}H\theta\\ \alpha\gamma\mu_{0}H\theta\end{array}\right).$$

$$\omega$$

$$\tau$$

$$\tan {\varphi }_0$$

$$\Delta \mathcal{F}=\mathcal{F}(+0)-\mathcal{F}(-\infty )\approx 3(A^2{\alpha }^2+1)\mathcal{D}{\mathbf{M}}_y^2,$$

$$\tan {\varphi }_0\approx A\left(\alpha +\frac{\delta }{\gamma {\mu }_{0}H\theta M_0}\right).$$

$$\Phi (t)=A_0+A_1{e}^{-t/t_1}+A_2{e}^{-t/\tau }\sin \left(\omega t+{\varphi }_0\right).$$

Full text

Efficient spin excitation via ultrafast damping-like torques in antiferromagnets

Christian Tzschaschel

[ctzschaschel@fas.harvard.edu

Current address: Department of Chemistry and Chemical Biology, Harvard University, Cambridge 02138, USA](mailto:ctzschaschel@fas.harvard.edu)

Department of Materials, ETH Zürich, 8093 Zurich, Switzerland

  

Takuya Satoh

Department of Physics, Tokyo Institute of Technology, Tokyo 152-8551, Japan

  

Manfred Fiebig

Department of Materials, ETH Zürich, 8093 Zurich, Switzerland

I Abstract

Damping effects form the core of many emerging concepts for high-speed spintronic applications. Important characteristics such as device switching times and magnetic domain-wall velocities depend critically on the damping rate. While the implications of spin damping for relaxation processes are intensively studied, damping effects during impulsive spin excitations are assumed to be negligible because of the shortness of the excitation process. Herein, we show that, unlike in ferromagnets, ultrafast damping plays a crucial role in antiferromagnets because of their strongly elliptical spin precession. In time-resolved measurements, we find that ultrafast damping results in an immediate spin canting along the short precession axis. The interplay between antiferromagnetic exchange and magnetic anisotropy amplifies this canting by several orders of magnitude towards large-amplitude modulations of the antiferromagnetic order parameter. This leverage effect discloses a highly efficient route towards the ultrafast manipulation of magnetism in antiferromagnetic spintronics.

II Introduction

Antiferromagnets (AFMs) are a promising material class for spintronic applications 1, 2. Their robustness against external magnetic fields, along with the possibility of picosecond antiferromagnetic switching, could help to substitute semiconductor-based information technologies with antiferromagnetic spintronic devices 3, 4. In fact, an electrically switched antiferromagnetic bit has been presented recently 5 and synthetic AFMs have substantially improved the speed and data density of racetrack memories 6. Efficient manipulation of AFMs relies on spin currents acting as effective magnetic fields generating staggered (anti-) damping-like torques 7. Such electrical control, however, has two disadvantages. On the one hand, the required current densities are in the range of $\mathsf{\sim 10^{6}-10^{8}\,A\,cm^{-2}}$, which involves a considerable amount of Joule heating. On the other hand, timescales accessible with electrical currents are typically too slow to explore the aforementioned picosecond dynamics. In contrast, optical laser pulses provide the time-resolution required to track THz-spin dynamics. The inverse Faraday effect (IFE), i.e. the optical generation of effective sub-picosecond magnetic field pulses with circularly polarised light, allows to impulsively create a spin canting in AFMs 8, 9, 10, 11, 12, 13, 14, 15.

For ferromagnetic excitations, the magnetic anisotropy typically leads to a circular or only slightly elliptical spin precession 16. In AFMs, however, the competition between magnetic anisotropy and the $\gtrsim$100 times stronger exchange interaction results in strongly elliptical trajectories. An oscillating net magnetisation occurs along the minor axis of the ellipse, whereas the antiferromagnetic order parameter is modulated due to a collective in-phase spin canting along the major axis 17, 18, 15. Therefore, the efficiency of impulsive spin excitations in antiferromagnets depends strongly on the direction of the initial spin canting. In particular, an impulsively generated net magnetisation in an AFM would be highly efficient in triggering large-amplitude changes of the antiferromagnetic order.

Antiferromagnetic spin dynamics are typically described by the Landau-Lifshitz-Gilbert equation 19

[TABLE]

where $\mathbf{M}_{i}$ denotes the magnetisation of the $i$-th antiferromagnetic sublattice with magnitude $M_{0}$; $\mathbf{H}_{i}$ is the field acting on $\mathbf{M}_{i}$, and $\alpha$ is the Gilbert damping parameter. Further, $\gamma$ and $\mu_{0}$ are the gyromagnetic ratio and magnetic permeability, respectively. The first term of Equation (1) is commonly referred to as the field-like torque $\mathcal{T}_{\mathrm{FL}}$, whereas the second term corresponds to the damping-like torque $\mathcal{T}_{\mathrm{DL}}$ 7. For fully compensated AFMs excited by the IFE, $\mathcal{T}_{\mathrm{FL}}$ modulates the antiferromagnetic order parameter, but does not induce a net magnetisation at $t=0$ (see Supplementary Movie 1). Therefore, the impulsive generation of a net magnetisation by the effective magnetic field of the IFE can only occur via an ultrafast damping-like torque $\mathcal{T}_{\mathrm{DL}}$ (see Supplementary Movie 2). So far, however, damping has been considered to describe only the spin-relaxation after the excitation 20. The ultrafast damping-like torque of the IFE was neglected for the excitation process itself based on the shortness of the interaction and the suspected smallness of the torque.

Here, we highlight the crucial role of the ultrafast damping-like torque during the excitation process for efficient impulsive spin excitations in AFMs. We detect such damping in the initial phase of a coherent spin precession of hexagonal YMnO3 and HoMnO3. We tune the spin damping strength via temperature variations and study the influence of magnetic order–order and order–disorder transitions on the damping. Our phenomenological theory considers both damping and magnetic anisotropy effects during the excitation, and thus develops and quantifies a consistent picture of impulsively induced coherent spin dynamics. Strikingly, we find that the hitherto neglected ultrafast spin damping during the impulsive spin excitation can even provide the dominant spin excitation mechanism.

III Results

III.1 Impulsive spin excitations in damped h-RMnO3

YMnO3 orders antiferromagnetically at the Néel temperature $T_{\mathrm{N}}\sim$73\text{,}\mathrm{K} with a transition from $P6_{3}cm1^{\prime}$ ($T>T_{\mathrm{N}}$) to $P6_{3}^{\prime}cm^{\prime}$ ($T<T_{\mathrm{N}}$) symmetry [21](#bib.bib21) (Fig. [1](#S3.F1)a). HoMnO3 undergoes a $P6_{3}cm1^{\prime}\rightarrow P6_{3}^{\prime}c^{\prime}m$ transition at $T_{\mathrm{N}}\sim$75\text{\,}\mathrm{K} followed by a first-order spin-reorientation transition (SRT) to $P6_{3}^{\prime}cm^{\prime}$ at $T_{\mathrm{SR}}\sim$35\text{,}\mathrm{K}$$ 22, 23. In a mean-field approach, the free energy density for both YMnO3 and HoMnO3 reads

[TABLE]

where $\mathbf{M}_{i}$ denotes the magnetisation at site $i$ ($i=1,2,3$). Furthermore $\lambda>0$ and $D_{z}>0$ parametrise the antiferromagnetic intra-plane exchange interaction and the easy-plane anisotropy, respectively. The inter-plane exchange is two orders of magnitude weaker than the intra-plane exchange; thus, we neglect it 24. $D_{y}<0$ stabilises the ground state as $P6_{3}^{\prime}c^{\prime}m$ ($\mathbf{M}_{i}$ then points along equivalent crystallographic $y$ axes), whereas $D_{y}>0$ stabilises $P6_{3}^{\prime}cm^{\prime}$ (with $\mathbf{M}_{i}$ along equivalent crystallographic $x$ axes). $\Delta\geq 0$ is a weak fourth-order in-plane anisotropy, which becomes relevant at the SRT, where $D_{y}$ crosses zero. The three-sublattice antiferromagnet exhibits three optically excitable spin precessions called $X$, $Y$, and $Z$ mode 25, 18, which relate to an oscillating net magnetisation along the $x$, $y$, or $z$ axes, respectively. We investigate the damping-like torque in relation to spin excitations via the IFE, which leaves us with the $Z$-mode excitation illustrated in Fig. 1c.

Figure 2a shows an exemplary measurement of the time-resolved Faraday rotation in YMnO3 after optical excitation with a circularly polarised pump pulse at $30\text{\,}\mathrm{K}$. Two regimes can be distinguished: (i) a pronounced peak around $0\text{\,}\mathrm{p}\mathrm{s}$ reflecting the direct interaction between pump and probe pulse and (ii) a damped sinusoidal modulation of the Faraday rotation that marks the ensuing spin precession and relaxation.

From a fit of the latter, we extract the relaxation time $\tau$, frequency $\omega$ and initial phase $\phi_{0}$ of the precession (see Methods). Figure 2b shows a magnified view of the region of the temporal overlap along with an extrapolation of the sinusoidal fit. Most strikingly, we find a finite magnetisation at $0\text{\,}\mathrm{p}\mathrm{s}$. This is unexpected because the conventional field-like torque of the IFE only perturbs the antiferromagnetic order, but does not induce a net magnetisation at the moment of the excitation 26, 27, 28, 14. We will show that the observed immediate onset of a finite net magnetisation is a consequence of spin damping during the impulsive spin excitation, which has been neglected so far.

Originally, YMnO3 is in the magnetic ground state with $\mathbf{M}_{i}(t<0)=M_{0}\mathbf{\hat{x}}_{i}$, where $M_{0}$ is the sublattice saturation magnetisation. The optomagnetic field pulse of the IFE can be modelled as $\mathbf{H}_{\mathrm{IFE}}(t)=H\theta\delta(t)\mathbf{\hat{z}}$, where $H$ is the effective magnetic field strength, $\theta$ is the laser pulse duration, and $\delta(t)$ denotes the Dirac delta function 29. Analytically integrating Equation (1) from $t=-\infty$ to $+0$ with $\mathbf{H}_{i}=\mathbf{H}_{\mathrm{IFE}}$ yields:

[TABLE]

The canting along the $z$ axis is typically significantly smaller than the $y$ canting because $\alpha\ll 1$. However, any spin canting along the magnetic hard axis ($z$) will be enhanced during the spin precession as a combined effect of antiferromagnetic exchange and magneto-crystalline anisotropy. This so-called exchange enhancement 17, 29, 14, 1, 30 is an effect that is intrinsic to AFMs, but is absent in ferromagnets. Thus, whereas damping effects during the spin excitation may be neglected for ferromagnets, the neglecting is not generally justified in AFMs. In fact, we will show that despite $\alpha\ll 1$, the damping-like contribution to the total spin excitation may even dominate over the field-like contribution.

The threefold rotational symmetry, which is preserved during the $Z$-mode precession, allows us to describe the antiferromagnetic dynamics by considering only one spin precessing in effective out-of-plane and in-plane anisotropies $\mathcal{J}=\frac{3}{2}\lambda+D_{z}+|D_{y}|$ and $\mathcal{D}=\Delta M_{0}^{2}+|D_{y}|$, respectively (Supplementary Note 1). For the $z$ component of the oscillating net magnetisation $\mathbf{M}$, this yields a damped sinusoidal oscillation with

[TABLE]

Here, $\hbar\omega_{0}=2g\mu_{\mathrm{B}}M_{0}\sqrt{\mathcal{J}\mathcal{D}}$ is the precession frequency for the undamped case with the Landé factor $g=2$ and the Bohr magneton $\mu_{\mathrm{B}}$. $A=\sqrt{\mathcal{J}\mathcal{D}^{-1}}\gg 1$ is the exchange enhancement factor 29, 1. Geometrically, $A$ corresponds to the aspect ratio of the elliptical spin precession 17, 14. Note that the initial phase $\phi_{0}$ depends on the damping during the impulsive spin excitation in direct relation to the initial phase in Fig. 2 and the ultrafast damping-like torque $\mathcal{T}_{\mathrm{DL}}$.

III.2 Ultrafast spin damping in HoMnO3

We further investigate the role of the damping-like torque by performing temperature-dependent measurements of the optically induced spin dynamics. Temperature affects both the Gilbert-damping parameter $\alpha$ (and therefore the magnitude of $\mathcal{T}_{\mathrm{DL}}$) and the magneto-crystalline anisotropy (and therefore the exchange enhancement $A$). In particular, we trace the $Z$ mode through the first-order SRT in HoMnO3 and towards the second-order antiferromagnetic-paramagnetic phase transition in YMnO3 (see Supplementary Fig. 2 for exemplary time-domain data). Conceptually, we thus investigate an order–order and order–disorder phase transition, thereby showcasing the general importance of spin damping during ultrafast antiferromagnetic spin excitations.

We first consider the case of HoMnO3. Based on measurements as in Fig. 2, we extract $\omega_{0}$, $\tau$ and $\phi_{0}$ as depicted in Figs. 3a–c, respectively. The SRT causes a dip of $\omega_{0}$ around $33\text{\,}\mathrm{K}$. The fact that the frequency does not drop to zero at the SRT highlights the importance of the fourth-order anisotropy term in Equation (2). All anisotropy parameters can be extracted from a fit of the temperature dependence of $\omega_{0}$ with a minimal model (see Methods). (The deviations below $25\text{\,}\mathrm{K}$ are caused by the incipient ordering of the Ho(4b) moments 31, 32, 24, see Supplementary Fig. 3).

Moreover, using Eqs. (4) and (5), we determine the exchange-enhanced damping $A\alpha$ and, in combination with the extracted anisotropy constants, its constituents $A$ and $\alpha$ (Figs. 3d–f). We find a significant increase in $A\alpha$ between $30\text{\,}\mathrm{K}$ and $45\text{\,}\mathrm{K}$. Together with Equation (6) we can then calculate the temperature dependence of the initial phase and compare it with the direct measurement (Fig. 3c). The agreement around the SRT and towards higher temperatures is impressive and demonstrates the importance of the previously unrecognised damping-like torque $\mathcal{T}_{\mathrm{DL}}$. More strikingly, combining Eqs. (2) and (3) yields the energy density transferred into the magnetic lattice during the excitation as

[TABLE]

where the first term originates from the ultrafast damping-like torque along the $z$ axis, whereas the second term is related to the field-like in-plane torque of the IFE. Thus, as $A\alpha$ exceeds one, spin excitation via the damping-like torque of the IFE even becomes the dominant excitation mechanism. The enhancement of $A$ around the SRT is based on changes of the magneto-crystalline anisotropy, while the anomalous increase of $\alpha$ matches the observations of a temporary small-domain state in passing the SRT 23. Such a small-domain state exhibits a high density of domain walls which will locally affect the spin precession frequency. Thus, we attribute the enhanced Gilbert damping to dephasing due to magnon scattering at the walls 33.

III.3 Ultrafast spin damping in YMnO3

The general trend of $\alpha$ is captured well by an offset Bose-Einstein distribution for a quasiparticle with an energy of $\sim 4.3\text{\,}\mathrm{m}\mathrm{e}\mathrm{V}$, suggesting that scattering with the crystal-field excitations of the Ho $4f$ electrons ($\approx$3.5\text{,}\mathrm{m}\mathrm{e}\mathrm{V}$$ 24) contributes significantly to the spin damping. This is confirmed by measurements on YMnO3 (Fig. 4), yielding relaxation times that are $\sim$20 times longer than those for HoMnO3 in line with the unoccupied $4f$ orbital. Accordingly, the extracted values for the temperature-dependent Gilbert-damping parameter $\alpha$ in YMnO3 (Fig. 5) are an order of magnitude lower than those in HoMnO3. At low temperature, $\alpha$ can be approximated by a Bose-Einstein distribution with a quasiparticle energy of $2.6\text{\,}\mathrm{m}\mathrm{e}\mathrm{V}$, which corresponds to a hybrid spin-lattice mode in YMnO334. Assuming that the frequency of this mode is, just like the $Z$ mode frequency, proportional to the magnitude of the antiferromagnetic order parameter, we can consistently describe the temperature dependence up to the Néel temperature.

Note, however, that despite good qualitative agreement, Fig. 4c reveals a difference between the measured and modelled value of the initial phase $\phi_{0}$. This indicates secondary, non-damping-based excitation mechanisms, such as an optical spin transfer torque 35 or optical orientation 36, which could become relevant because of the smallness of $\alpha$ in YMnO3 as compared with that in HoMnO3. These mechanisms involve an angular momentum transfer from the circularly polarised pump pulse to the material. As the angular momentum of the circular laser pulse is directed along the $z$ axis, we heuristically include such a contribution in our model by setting $M_{i,z}(+0)=\alpha\gamma\mu_{0}H\theta M_{0}+\delta$ in Equation (3), which changes the initial phase to

[TABLE]

With this, we restore the quantitative agreement between the measured and modelled initial phase in Fig. 4c.

IV Discussion

Summarising, we thus demonstrated the indispensable role of the previously neglected ultrafast damping-like torque during the optical excitation of AFMs, both experimentally and theoretically. The torque impulsively generates a magnetisation along the minor axis of the intrinsically highly elliptical spin precession in AFMs. The extreme ellipticity converts the small damping-related magnetisation into a large precessional canting. Because of this leverage effect, this damping-based mechanism can even facilitate the dominant spin excitation mechanism in striking contrast to its complete neglect thus far. In addition to studying the spin dynamics deep in the antiferromagnetic phase, we also traced the associated $Z$ mode in the vicinity of both first- and second-order phase transitions. This enabled us to verify the role of the ultrafast damping-like torque in distinctly different scenarios. We find that the Gilbert-damping parameter $\alpha$ remains almost constant throughout the SRT in HoMnO3, while the changes of the magneto-crystalline anisotropy driving the SRT lead to a two-fold increase in the exchange enhancement factor $A$. In contrast, as we approach the Néel temperature $T_{N}$ (here in YMnO3), the magneto-crystalline anisotropy remains unchanged, but $\alpha$ diverges. Although fundamentally different, we find an enhanced ultrafast damping-like torque in both cases, which again highlights the general importance of this hitherto neglected excitation mechanism.

More specifically, in terms of damping mechanisms we identified scattering with the $4f$ crystal-field excitations of HoMnO3 as the leading contribution to the spin damping. Therefore, the choice of the A-site ion enables us to tune the strength of the damping across several orders of magnitude. This is, on one hand, an useful degree of freedom for tuning material properties to specific spintronics applications. On the other hand, it allowed us to identify the presence of additional non-damping-based mechanisms inducing a net magnetisation. With respect to this, we note that any ultrafast mechanism that induces a net magnetisation in an AFM will benefit from the observed exchange enhancement. While field-like excitations via the IFE or using THz pulses typically induce spin-cantings of the order of $1^{\circ}$37, 15, there are various possibilities beyond increasing the pump fluence that may allow us to exceed a $90^{\circ}$ spin canting and thus induce antiferromagnetic switching. Exploiting different excitation mechanisms such as thermally and optically induced ultrafast spin-transfer torques 38, 35, 39 may drastically increase the impulsively induced net magnetisation. Increasing the exchange enhancement, e.g. by tuning the system close to a SRT, will increase the ellipticity of the spin precession. Therefore, the impulsively induced net magnetisation translates into a larger modulation of the antiferromagnetic order parameter thereby reducing the switching threshold. Ultimately, our results show that both field-like and damping-like torques — a combination that has proved to be powerful for the electrical control of magnetic order 7, 40 — are available optically, too, and can be utilized to act on the magnetic order. Thus, optically generated torques might provide the long sought-after tool enabling the efficient realisation of ultrafast coherent precessional switching of AFMs.

V Methods

V.1 Setup

Both samples are flux-grown and polished, (0001)-oriented platelets of approximately $55\text{\,}\mathrm{\SIUnitSymbolMicro m}$ thickness. Our setup is schematically shown in Fig. 1b. The fundamental light source is a regeneratively amplified laser system providing $130\text{\,}\mathrm{f}\mathrm{s}$ pulses at $1.55\text{\,}\mathrm{e}\mathrm{V}$ photon energy and $1\text{\,}\mathrm{k}\mathrm{H}\mathrm{z}$ repetition rate. We use the output of an optical parametric amplifier combined with a quarter-wave plate to generate circularly polarised pump pulses at $0.95\text{\,}\mathrm{e}\mathrm{V}$ photon energy, which create an effective magnetic field in the sample via the IFE. The pump fluence on the sample is approximately $60\text{\,}\mathrm{m}\mathrm{J}\,\mathrm{c}\mathrm{m}^{-2}$. By pumping the material far below the nearest absorption resonance at $\approx$1.6\text{,}\mathrm{e}\mathrm{V} [41](#bib.bib41), we minimise parasitic heating effects [42](#bib.bib42). We then probe the ensuing dynamics by measuring the time-resolved Faraday rotation $\Phi(t)$ of a linearly polarised probe pulse at $1.55\text{\,}\mathrm{e}\mathrm{V}$ with balanced photodetectors. The probe fluence is $<0.1\text{,}\mathrm{m}\mathrm{J},\mathrm{c}\mathrm{m}^{-2}$. Both the pump and probe beams are arranged in a quasi-collinear geometry with the beam propagation direction along the hexagonal$z$axis, which is the optic axis. Thus, the incident light experiences optical isotropy and complications arising from birefringence or other anisotropic optical effects are avoided. The sample is mounted in a cryostat for varying the temperature between$\approx5\text{\,}\mathrm{K}$$ and 300\text{,}\mathrm{K}$.

The oscillatory part of the signal (Fig. 2a) is fitted by

[TABLE]

While the first two terms model parasitic contributions to the signal arising from thermal electron dynamics, we extract the frequency $\omega(T)$, relaxation time $\tau(T)$ and initial phase $\phi_{0}(T)$ from the oscillatory third term. Error bars in Figs. 3, 4 and 5 reflect the standard error of that fit or its propagation through Eqs. (4)-(6).

V.2 Temperature dependence of $\omega_{0}$

As stated in the main text, the analytic solution of the Landau-Lifshitz-Gilbert equation yields (Supplementary Note 1): $\hbar\omega_{0}=2g\mu_{\mathrm{B}}M_{0}\sqrt{\mathcal{J}\mathcal{D}}$. According to mean-field theory, the temperature-dependent sublattice magnetisation $M_{0}(T)$ can be described as a paramagnet aligned by the exchange field of the neighbouring sublattices, i.e. $M_{0}(T)=Ng\mu_{\mathsf{B}}S\mathcal{B}_{S=2}(M_{0}(T)T_{\mathrm{N}}T^{-1})$, where $\mathcal{B}_{S=2}$ is the Brillouin function of an $S=2$ Mn spin 16, 43, 44 and $N=6/V$ with $V=$0.374\text{,}\mathrm{n}\mathrm{m}^{3}$$, being the unit cell volume, is the number density of Mn atoms. We model the SRT in HoMnO3 by introducing a linear dependence of the anisotropy parameter $D_{y}$ on temperature. The SRT temperature $T_{\mathrm{SR}}$ is defined as the point where $D_{y}$ vanishes 24, i.e. $D_{y}(T)=\kappa(T-T_{\mathrm{SR}})$. For YMnO3, assuming a temperature-independent value for $D_{y}$ yields a reasonable fit. This model is fitted to the measured temperature dependence of $\omega_{0}$. The results are summarised in Table 1.

VI Data availability

The data that support the findings of this study are available from the corresponding author on reasonable request. Source data are provided with this paper.

VII Acknowledgment

The authors thank B. A. Ivanov and A. V. Kimel for valuable discussions and R. Gmünder for experimental assistance. T.S. was supported by the Japan Society for the Promotion of Science (JSPS) KAKENHI (Nos. JP19H01828, JP19H05618, JP19K21854 and JP26103004) and the JSPS Core-to-Core Program (A. Advanced Research Networks). C.T. acknowledges support by the SNSF under fellowship P2EZP2-191801. C.T. and M.F. acknowledge support from the SNSF project 200021/147080 and by FAST, a division of the SNSF NCCR MUST.

VIII Author contributions

C.T. and T.S. conceived the project. C.T. designed and evaluated the experiments. T.S. and M.F. supervised the project. All authors discussed the results and contributed to writing the manuscript.

IX Corresponding Author

Correspondence to Christian Tzschaschel.

X Competing interests

The authors declare no competing interests.