## Universal prethermal dynamics in Heisenberg ferromagnets

Saraswat Bhattacharyya, Joaquin F. Rodriguez-Nieva, and Eugene Demler Department of Physics, Harvard University, Cambridge, MA 02138, USA

We study the universal far from equilibrium dynamics of magnons in Heisenberg ferromagnets. We show that such systems exhibit universal scaling in momentum and time of the quasiparticle distribution function, with the universal exponents distinct from those recently observed in Bose-Einstein condensates. This new universality class originates from the SU(2) symmetry of the Hamiltonian, which leads to a strong momentum-dependent magnon-magnon scattering amplitude. We compute the universal exponents using the Boltzmann kinetic equation and incoherent initial conditions that can be realized with microwave pumping of magnons. We compare our numerical results with analytic estimates of the scaling exponents and demonstrate the robustness of the scaling to variations in the initial conditions. Our predictions can be tested in quench experiments of spin systems in optical lattices and pump-probe experiments in ferromagnetic insulators such as yttrium iron garnet.

Introduction. Understanding the emergence of universal dynamics in interacting quantum systems far from thermodynamic equilibrium is a central challenge in theoretical physics. A large body of theoretical works have proposed that isolated quantum systems, such as quarkgluon plasma after heavy ion collision, the early universe after inflation, or cold atoms after a quench, can exhibit universal relaxation of quasiparticles as they evolve far from equilibrium[1-13]. Such behavior is manifested in universal scaling in time and momenta of correlation functions in the prethermal regime, with the universal exponents independent of microscopic details or initial conditions. The universal relaxation can be attributed to the existence of non-thermal fixed points in the system's phase space (Fig.1), with the non-equilibrium state inheriting its universal exponents[14]. Manifestations of universal relaxation were recently observed experimentally for the first time in cold atomic gases[15-17]. This experimental feat creates new challenges, both in classifying all the possible universality classes as well as devising new tabletop experiments to explore them.

presses collisions between quasiparticles and the condensate that arises due to symmetry breaking, preventing the renormalization of the quadratic dispersion of quasiparticles [unlike a BEC where Goldstone modes have linear dispersion]. These two features lead to distinct universal exponents in a broad range of wavevectors.

Here we uncover a new universality class arising in the relaxation dynamics of magnons in the Heisenberg model. Similar to the case observed experimentally with an interacting Bose-Einstein condensate (BEC) [15-17], the prethermal quasiparticle distribution becomes selfsimilar in momentum k and time t :

$$n _ { k } ( t ) = t ^ { \alpha } f ( t ^ { \beta } | k | ),$$

where α and β are universal exponents independent of microscopic details or initial conditions, and f ( x ) is a universal function. Key properties that determine α and β are dimensionality, quasiparticle dispersion and the nature of interactions. While a ferromagnet at low energy can be described as an interacting Bose gas after a Holstein-Primakoff transformation[18], the SU(2) symmetry of the Hamiltonian sets the Heisenberg model apart from a conventional BEC in two important ways. First, interaction between quasiparticles are strongly constrained by SU(2) symmetry giving rise to 'soft' collisions [Eq.(4) below]. Second, the SU(2) symmetry sup-

Crucially, we find that the key requirement for the observation of the universal scaling in Eq.(1) is that a sufficiently large population of magnons is pumped into the system-the details of the initial conditions are unimportant. In terms of the experimentally controllable density of magnons ρ and pumping frequency ω , the occupation number of modes at frequency ω is given by n ω ≈ ρa d ( JS/ω ) d/ 2 , where we assumed that magnons are pumped in an energy window ∆ ω ∼ ω and have dispersion ω k = JS ( | k | a ) 2 ( a : lattice constant, J : exchange coupling, d : system dimension, S : spin number). The condition n ω ∼ 1 defines a boundary below which we do not observe self-similarity and the system evolves directly into the thermal fixed point, see Fig.1. Furthermore, when pumping quasiparticles close to the bottom of the band, we expect interactions that break SU(2) symmetry, e.g. , dipolar or anisotropic exchange interactions, to be important (Fig.1). In this case, we expect different

FIG. 1. The presence of a non-thermal fixed point in phase space can induce a long-lived prethermal state. Appearance of the prethermal state starting from an incoherent state ˆ ρ inc depends on the average density ρ of magnons pumped into the system and the frequency ω at which they are pumped, see discussion in the main text. Here J is the exchange coupling, g is the strength of the interaction that breaks SU(2) symmetry and a is the lattice constant.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00554v2-figure-f6615ff3-7ae2-485f-ac0b-eb192bedff97.jpg)

universal exponents which will not be discussed here.

The prethermal state lives for a wide and experimentally accessible temporal window given the slow relaxation of long wavelength Goldstone modes in a ferromagnet. In particular, the thermalization time is determined by terms that break SU(2) symmetry, e.g. , dipolar interactions. For example, in the context of BoseEinstein condensation of magnons in yttrium iron garnet (YIG) after microwave pumping[19-21], it was found a thermalization time τ ∼ 10 -100 ns that is consistent with dipolar interactions, which are much slower than microscopic timescales associated to the exchange coupling J ∼ 100 meV.

Beside of its fundamental appeal, our predictions are relevant to various experiments. For example, we argue that the universal scaling exponents can be accessed in YIG [19-23] given its negligible magnetic anisotropy and despite its ferrimagnetic order[24]. In particular, there is a large energy window in which quasiparticles can be pumped such that: (i) the collisions rate due to exchange coupling is much faster than the collision rate due to terms that break SU(2) symmetry and (ii) a single parabolic spin wave mode is populated ( ω glyph[lessmuch] 100 meV). Cold atom platforms are also promising because the system can be effectively isolated from the environment and the exchange interaction can be engineered using various mechanisms, e.g. Feshbach resonances, dipolar interactions or lattice shacking[25-31].

Microscopic model. Focusing only on a single magnon band, we consider a two-dimensional Heisenberg ferromagnet on a square lattice with nearest neighbour exchange and Zeeman field

$$\hat { \mathcal { H } } = - J \sum _ { \langle j j ^ { \prime } \rangle } \hat { S } _ { j } \cdot \hat { S } _ { j ^ { \prime } } + h _ { z } \sum _ { j } \hat { S } _ { j } ^ { z }, \quad \ ( 2 ) \quad \ \frac { \text{Equ} } { \text{vani} }$$

with J &gt; 0 and 〈 jj ′ 〉 denoting nearest neighbors. We assume that the system has N lattice sites, each containing a spin S degree of freedom, and periodic boundary conditions in each spatial direction. The spin operators satisfy the commutation relations [ ˆ S , S z j ˆ ± j ′ ] = ± δ jj ′ ˆ S ± j and [ ˆ S + j , S ˆ -j ′ ] = 2 δ jj ′ ˆ S z j , with ˆ S ± j = ˆ S x j ± iS ˆ y j . We also assume that the lattice is at a small temperature such that magnon-phonon interactions can be neglected, which is the case for the prethermal timescales of interest t &lt; ∼ 1 µ s [32, 33]. We include a Zeeman field h z , which is present in many relevant experiments and seems to break the SU(2) symmetry, to illustrate that h z has no effect on the relaxation dynamics.

We proceed to build an effective theory valid when the density of quasiparticles is small, ρa 2 glyph[lessmuch] S . We recall that one magnon states | k 〉 = ˆ S + k | F are exact eigenstates 〉 of ˆ with energies H

$$\varepsilon _ { k } = h _ { z } + J S ( \gamma _ { 0 } - \gamma _ { k } ), \ \gamma _ { k } = \sum _ { \tau } e ^ { i k \cdot \tau }. \quad ( 3 ) \quad ( \omega = \quad \text{mag}$$

Here | F 〉 = | ↓↓ . . . ↓〉 denotes the ferromagnetic ground state and ˆ S + k denotes ˆ S + k = 1 √ N ∑ j e -i k r · j ˆ S + j . Two magnon states | k p , 〉 = 1 2 S ˆ S + k ˆ S + p | F , 〉 however, are not eigenstates of ˆ [18, H 34]. The interaction between magnons can be obtained from the matrix elements

$$\stackrel { \sim } { \mathfrak { e } } _ { \mathfrak { e } } _ { \mathfrak { e } t } \quad \hat { \mathcal { H } } | k, p \rangle = ( \varepsilon _ { k } + \varepsilon _ { p } ) | k, p \rangle + \sum _ { q } G _ { k, p } ^ { q } | k + q, p - q \rangle, \\ \intertext { a n a t o } \text{$G_{k,p}^{q}=\frac{J}{N} \left(\gamma_{q-p}+\gamma_{q+k}-\gamma_{q}-\gamma_{k+q-p})$}, \\ \text{an} _ { u } \quad \text{such that one magnon states are counted via momentum}.$$

such that one magnon states are coupled via momentumconserving collision G q k p , . When the incoming magnons have long wavelength, the collision term takes the simple form G q k p , ≈ -Ja 2 N ( k · p ) for all values of S . This characteristic ( k · p ) interaction, which is independent of h z , arises from the global SU(2) symmetry of the exchange coupling and justifies why magnons propagate ballistically when | k | → 0. Hard core collisions may arise if the SU(2) symmetry is broken.

A long-wavelength description that captures the features of the SU(2) symmetric parent Hamiltonian is

$$\intertext { s y - } \intertext { t and } \intertext { r i o u s } \intertext { \frac { \hat { \mathcal { H } } } { } = \sum _ { k } \varepsilon _ { k } \hat { a } _ { k } ^ { \dagger } \hat { a } _ { k } - \frac { J a ^ { 2 } } { N } \sum _ { k p q } ( k \cdot p ) \hat { a } _ { p + q } ^ { \dagger } \hat { a } _ { k - q } ^ { \dagger } \hat { a } _ { p } \hat { a } _ { k } + h. c., } \\ \intertext { e r a c - } \intertext { \text{ show a. i.o. b o c o n o u n o d f o u n d a t o w e t i o l t i o n } }$$

where ˆ a k is a bosonic operator defined after a HolsteinPrimakoff transformation of the spin operators. In contrast to the Bogoliubov theory of weakly interacting Bose particles, Eq.(5) explicitly shows the absence of anomalous terms ( a † -k a † k ) which describe coherent scattering between finite energy quasiparticles and the condensate. Equation (5) also shows that scattering of quasiparticles vanishes as their momentum approaches zero, which precludes the formation of a magnon condensate at k = 0.

Initial conditions. We consider an initial incoherent population of magnons in a narrow band of energies centered around ω . Such initial condition can be achieved via transverse microwave pumping h ⊥ ( ) at fret quency Ω = 2 ω [35]. The amount of magnons pumped into the systems can be controlled by the strength of h ⊥ , giving two independent knobs to control which k modes are excited and their respective population n k . Although parametric pumping of magnons also create anomalous correlations 〈 ˆ a † -k ˆ a † k 〉 , these decohere rapidly since pairs with different wavevectors k oscillate with different frequencies. We also note that this protocol leads to no net spin texture, 〈 ˆ a k 〉 = 0. We parametrize the initial condition as

$$n _ { k } ( t = 0 ) = n _ { * } \exp \left [ - \frac { ( | { \mathbf k } | - { k } _ { * } ) ^ { 2 } } { \Gamma ^ { 2 } } \right ], \quad \ \ ( 6 )$$

where k ∗ is the wavevector at which magnons are pumped ( ω = k / m 2 ∗ 2 ), n ∗ parametrizes the occupation number of magnons at k ∗ , and Γ determines the initial width of

the distribution (its value depends on the details of the pump pulse, e.g. , its duration). We note that, although n ∗ can be much larger than 1, n ∗ and Γ need to satisfy ρa 2 ≈ k ∗ Γ a n 2 ∗ glyph[lessmuch] S for Eq.(5) to be valid.

Kinetic equation. The measurable quantity of interest is the magnon population n k ( ) t = 〈 ˆ a † k ( )ˆ t a k ( ) t 〉 as a function of time. In ferromagnetic materials, such quantity can be measured via Brillouin scattering[20, 21]. An alternative technique is spin qubit magnetometry[36], which has been used to measure (steady-state) magnon population[37, 38] as well as imaging single spins[39], but also has been proposed to access a variety of elementary excitations in ferromagnets[40, 41], spin ice[42], spin chains[43], and spin liquids[44]. In cold atom experiments, it is possible to use snapshots of local spin measurements 〈 ˆ S S x i ˆ x j 〉 for the different spin pairs in order to compute n k .

The time evolution of n k ( ) at intermediate timescales t (after decoherence of anomalous terms has occured) can be described using the kinetic equation ∂ n t k = I k , with I k the collision integral:

$$I _ { k } = J ^ { 2 } a ^ { 8 } & \int _ { p } \int _ { q } ( k \cdot p ) ^ { 2 } \left [ n _ { k } n _ { p } ( 1 + n _ { k + q } ) ( 1 + n _ { p - q } ) \quad \text{FIG.} \quad \text{ing fr} \\ & - ( 1 + n _ { k } ) ( 1 + n _ { p } ) n _ { k + q } n _ { p - q } \right ] \delta ( \varepsilon _ { i } - \varepsilon _ { f } ). \quad \text{$k_{*}$, w}$$

Here ε i = ε k + ε p and ε f = ε k + q + ε p -q the energies of the initial and final states, respectively, and ∫ p ≡ ∫ d 2 p (2 π ) 2 . The kinetic equation is valid in the weak coupling regime, ∫ k | I k | ∫ k ω n k k glyph[lessmuch] 1, which is true if ρa 2 glyph[lessmuch] S . We note that the lack of condensate formation at k = 0 justifies the kinetic equation at long times because, otherwise, condensate formation would give rise to non-perturbative corrections to the kinetic equation.

Universal exponents from kinetic simulations. We first focus on prethermalization in the large pumping regime, n ∗ glyph[greatermuch] 1, where the collision term (7) defines a timescale τ ∗ = J 2 ω n 3 3 ∗ . As shown in Fig.2(a), on a timescale O ( τ ∗ ), the details of the initial conditions are lost and the distribution function at intermediate momenta, | k | ∼ k ∗ , acquires a self-similar form governed by Eq.(1). We find that the distribution function n k can be fitted by Eq.(1) in a broad range of momenta (between one and two decades) using the parameters

$$\alpha = 0. 6 5 \pm 0. 0 5, \ \beta = 0. 3 0 \pm 0. 0 5, \quad \ ( 8 ) \quad _ { \ e a a l ; } ^ { a }$$

and the universal function f ( x ) ∼ 1 /x 2 3 . . The uncertainty in Eq.(8) is obtained by initializing the simulation from qualitatively distinct initial conditions (but same energy and particle number) and computing the variations in ( α, β ). Figure 2(a) also shows the lack of scattering of magnon states at | k | ≈ 0. For small momenta, the distribution function evolves as ∂ n t k ∼ | t k | 2 . Such behavior is cutoff by interactions that break SU(2) symmetry. We emphasize that the exponents ( α, β ) in Eq.(8) are

FIG. 2. Evolution of the occupation number n k starting from an initial incoherent pump at wavevector | k | = k ∗ , with occupation (a) n ∗ = 100 and (b) n ∗ = 1 [see Eq.(6)]. Indicated with dashed-dotted lines is the distribution function once the details of the initial conditions are lost. The inset (a) illustrates the collapse of the data points using a self-similar distribution function in Eq.(1), with α and β defined in Eq.(8). The inset (b) exhibits no collapse of the data points. Solid lines indicate the initial (black) and final Bose-Einstein (light blue) distribution of the magnon gas. The distributions are plotted at times t/τ ∗ = 0 0 01 0 02 0 09 0 11 0 13 0 16 0 18 0 22 0 28 , . , . , . , . , . , . , . , . , . , ∞ for (a) and t/τ ′ ∗ = 0 0 01 0 02 0 5 0 1 0 15 0 2 0 25 , . , . , . , . , . , . , . , ∞ for (b), with decreasing tones of blue [ τ ∗ = J /ω n 2 3 3 ∗ and τ ′ ∗ = J /ω n 2 3 2 ∗ ]. In both panels, we normalized n k with the BoseEinstein distribution at k = 0 consistent with the energy and particle number of the initial state. Parameters used: Γ = 0 2 . k ∗ , momentum cutoff Λ = 4 k ∗ (for n ∗ = 100 and ω/J = 0 01, . ρa 2 ≈ 0 1). .

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00554v2-figure-fe528c6d-b3de-42ae-951f-154355d91326.jpg)

different from those found in cold atom experiments[1517], a consequence of different scattering dynamics and different dimensionality (for instance, α ≈ β ≈ 0 1 in a . 1D Bose-Einstein condensate in Ref.[16]).

As shown in Fig.2(b), we do not observe self-similar scaling in our numerical results at small occupation number ( n ∗ = 1). In particular, we observe that the distribution at intermediate/large momenta relax to the thermal form n k ∝ e -ε k /T without exhibiting self-similarity (small momenta states, | k | glyph[lessmuch] k ∗ , still scales as n k ∝ t | k | 2 ), in agreement with previous results on different models[5, 13].

Importantly, the universal exponents are independent of the details of the initial condition. Figure 3(a) shows the evolution of the quasiparticle distribution after an

FIG. 3. (a) Robustness of the self-similarity under different initial conditions. Shown is the evolution of the distribution function using a two-peak initial condition. After a short time, the distribution function becomes self-similiar like in Fig.2. (b) Thermalization in the presence of SU(2) symmetry breaking terms which gives rise to scattering at small momenta. We do not observe a noticeable change in the ( α, β ) values from those found in Eq.(8). Parameters used are described in the main text.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00554v2-figure-f406685f-e8c5-4b35-86ce-efd173d53eb1.jpg)

initial pump at two frequencies. Similarly to Fig.2, the details of the initial conditions are lost in a time scale O ( τ ∗ ) and the system evolves in a self-similar fashion with the same universal exponents in Eq.(8).

lengthscales in the relaxation of spin models is a future direction to explore.

Universal exponents from dimensional analysis. The scaling exponents can be analytically estimated using arguments based in wave turbulence[45, 46]. In the presence of two positive conserved quantities (particle number and energy), a dual cascade typically occurs in two different regions of k -space, one with uniform flux of energy and the other with uniform flux of particles. If we assume that a particle cascade dominates at wavevectors | k | ∼ k ∗ , then magnon number conservation, ρ = ∫ d d k (2 π ) d n k ( ) = t t α -βd ∫ d d κ f ( | κ | ), implies α = dβ ( κ = t β k ). Alternatively, if an energy cascade dominates at | k | ∼ k ∗ , then energy conservation, E = ∫ d d k (2 π ) d ε k n k ( ) t = t α -β d ( +2) ∫ d d κ f ( | κ | ), implies α = β d ( + 2). A second relation between ( α, β ) is obtained from the kinetic equations assuming that n k glyph[greatermuch] 1: the collision term (7) can be rescaled as I k = t 3 α -4 β -2 dβ +2 β I κ , and ∂ n t k can be rescaled as ∂ n t k = t α -1 [ αf ( | κ | ) + βκf ′ ( | κ | )]. Matching the coefficients of t of these two terms results in 2( d +1) β -2 α = 1.

The universal exponents found numerically [Eq.(8)] are modestly close to those associated with a particle cascade ( α p = 1 , β p = 0 5); however, we do not see in our sim-. ulations a sharp separation between energy and particle cascade. The lack of a dual cascade appears to be linked to a combination of two effects: (i) absence of an emergent lengthscale in Eq.(5) ( i.e. , the only lengthscale in the problem is k ∗ defined in the initial condition) and (ii) the arrested dynamics at intermediate | k | due to absence of scattering at k → 0. This is distinct from a weakly interacting Bose gas where there is an interaction lengthscale 1 /λ = √ mgρa 2 (which defines a boundary between energy and particle cascades[13]; m : mass] and the unconstrained particle cascade towards k → 0. Incorporating real space dynamics to understand the emergent

SU(2) symmetry breaking terms. Interactions that break SU(2) symmetry allow to populate the k ≈ 0 modes and give rise to different universal exponents if they dominate over Heisenberg exchange. One common interaction is exchange anisotropy, ˆ H z = δJ z ∑ 〈 i,j 〉 ˆ S S z i ˆ z j , which reduces the symmetry of the Hamiltonian from SU(2) to U(1) and effectively gives rise to hard-core collisions ˆ H z = δJ z N ∑ k p q , , ˆ a † k + q ˆ a † p -q ˆ a k ˆ a p . Dipolar interaction, in addition, breaks the remaining U(1) symmetry and gives rise to anomalous terms in the Hamiltonian. So long as quasiparticles are pumped at energy scales in which Heisenberg exchange dominate over interactions that break SU(2) symmetry, we expect the exponents ( α, β ) in Eq.(8) to hold.

To confirm that this is indeed the case, Fig.3(b) shows the magnon relaxation in the presence of a weak exchange anisotropy δJ z = 0 05 . J k a ( ∗ ) 2 and n ∗ glyph[greatermuch] 1. Contrary to the results above, the k ≈ 0 modes are populated. However, the universal exponents at intermediate momenta, where Heisenberg exchange dominates, remain within the values found in Eq.(8).

Summary &amp; outlook. We showed that the Heisenberg model hosts universal quasiparticle relaxation after an incoherent pump. One direction to explore is whether other universal regimes are possible under different classes of initial conditions, e.g. , spin textures [47, 48] or in the presence of orbital degrees of freedom[49]. Other open questions are whether the same self-similar scaling survives in the large magnon density regime, i.e. as we approach criticality, and its connection with coarsening dynamics and ageing[50, 51]. Such studies need to go beyond the kinetic equation, for instance, using Truncated Wigner Approximation for spin systems[52-55]. On the experimental front, probing the predicted non-thermal fixed point in ferromagnets is within grasp of ongoing experiments, namely driven ferromagnetic insulators and quenches of spins in optical lattices.

Acknowledgements. We thank C. Du, G. Falkovich, B. Halperin, J. Marino, A. Pi˜eiro Orioli, D. Podolsky, n A. M. Rey, A. A. Rosch, D. Sels, A. Yacoby, T. Zhou for enlightening discussions. We acknowledge support from Harvard-MIT CUA, NSF Grant No. DMR-1308435 and AFOSR-MURI: Photonic Quantum Matter (award FA95501610323).

- [1] R. Micha and I. I. Tkachev, Phys. Rev. Lett. 90 , 121301 (2003).
- [2] R. Micha and I. I. Tkachev, Phys. Rev. D 70 , 043538 (2004).
- [3] J. Berges, S. Bors´nyi a and C. Wetterich, Phys. Rev. Lett. 93 , 142002 (2004).
- [4] J. Berges, D. Gelfand, and J. Pruschke, Phys. Rev. Lett. 107 , 061301 (2011).
- [5] J. Berges, A. Rothkopf, and J. Schmidt, Phys. Rev. Lett. 101 , 041603 (2008).
- [6] J. Berges and D. Sexty, Phys. Rev. D 83 , 085004 (2011).
- [7] J. Berges, K. Boguslavski, S. Schlichting, R. Venugopalan, Phys. Rev. D 89 , 074011 (2014).
- [8] A. Kurkela and G. D. Moore, Phys. Rev. D 86 , 056008 (2012).
- [9] M. C. Abraao York, A. Kurkela, E. Lu, and G. D. Moore, Phys. Rev. D 89 , 074036 (2014).
- [10] S. Schlichting, Phys. Rev. D 86 , 065008 (2012).
- [11] J. Berges, K. Boguslavski, S. Schlichting, and R. Venugopalan, Phys. Rev. D 89 , 114007 (2014).
- [12] A. Kurkela and E. Lu, Phys. Rev. Lett. 113 , 182301 (2014).
- [13] A. Pi˜eiro Orioli, K. Boguslavski, and J. Berges, Phys. n Rev. D 92 , 025041 (2015).
- [14] J. Berges, arXiv:1503.02907.
- [15] M. Pr¨fer, u et al. , Nature 563 , 217 (2018).
- [16] C. Eigen, et al. , Nature 563 , 221 (2018).
- [17] S. Erne, R. B¨ ucker, T. Gasenzer, J. Berges, and J. Schmiedmayer, Nature 563 , 225 (2018).
- [18] D. C. Mattis, The Theory of Magnetism Made Simple (World Scientific, 2006).
- [19] S. O. Demokritov, et al. , Nature 443 , 430 (2006).
- [20] V. E. Demidov, et al. , Phys. Rev. Lett. 99 , 037205 (2007).
- [21] A. V. Chumak, et al. , Phys. Rev. Lett. 102 , 187205 (2009).
- [22] A. A. Serga, et al. , Nature Communications 5 , 3452 (2014).
- [23] A. J. E. Kreil, et al. , Phys. Rev. Lett. 121 , 077203 (2018).
- [24] A. J. Princep, R. A. Ewings, S. Ward, S. T´th, C. Dubs, o D. Prabhakaran, and A. T. Boothroyd, npj Quantum Materials 2 , 63 (2017).
- [25] L.-M. Duan, E. Demler, and M. D. Lukin, Phys. Rev. Lett. 91 , 090402 (2003).
- [26] S. Trotzky, et al. , Science 319 , 295 (2008).
- [27] B. Yan, et al. , Nature 501 , 521 (2013).
- [28] C. V. Parker, L.-C. Ha, and C. Chin, Nature Physics 9 ,

769 (2013).

- [29] S. Hild, et al. , Phys. Rev. Lett. 113 , 147205 (2014).
- [30] C.-L. Hung, A. Gonz´lez-Tudela, J. I. Cirac, and H. J. a Kimble, Proc. Natl. Acad. Sci. 113 , E4946 (2016).
- [31] E. J. Davis, G. Bentsen, L. Homeier, T. Li, and M. H. Schleier-Smith, Phys. Rev. Lett. 122 , 010405 (2019).
- [32] J. Xiao, G. E. W. Bauer, K.-c. Uchida, E. Saitoh, and S. Maekawa, Phys. Rev. B 81 , 214418 (2010).
- [33] M. Agrawal, V. I. Vasyuchka, A. A. Serga, A. D. Karenowska, G. A. Melkov, B. Hillebrands, Phys. Rev. Lett. 111 , 107204 (2013).
- [34] F. J. Dyson, Phys. Rev. 102 , 1217 (1956).
- [35] A. G. Gurevich and G. A. Melkov, Magnetization oscillations and waves (CRC Press, 1996).
- [36] C. L. Degen, F. Reinhard, and P. Cappellaro, Rev. Mod. Phys. 89 , 035002 (2017).
- [37] T. van der Sar, F. Casola, R. Walsworth, and A. Yacoby, Nat Commun 6 , 7886 (2015).
- [38] C. Du, T. van der Sar, T. X. Zhou, P. Upadhyaya, F. Casola, H. Zhang, M. C. Onbasli, C. A. Ross, R. L. Walsworth, Y. Tserkovnyak, and A. Yacoby, Science 357 , 195 (2017).
- [39] M. S. Grinolds, S. Hong, P. Maletinsky, L. Luan, M. D. Lukin, R. L. Walsworth, and A. Yacoby, Nat Phys 9 , 215 (2013).
- [40] B. Flebus and Y. Tserkovnyak, Phys. Rev. Lett. 121 , 187204 (2018).
- [41] J. F. Rodriguez-Nieva, D. Podolsky, and E. Demler, arXiv:1810.12333.
- [42] F. K. K. Kirschner, F. Flicker, A. Yacoby, N. Y. Yao, and S. J. Blundell, Phys. Rev. B 97 , 140402 (2018).
- [43] J. F. Rodriguez-Nieva, K. Agarwal, T. Giamarchi, B. I. Halperin, M. D. Lukin, and E. Demler, Phys. Rev. B 98 , 195433 (2018).
- [44] S. Chatterjee, J. F. Rodriguez-Nieva, and E. Demler, Phys. Rev. B 99 , 104425 (2019).
- [45] S. Nazarenko, Wave Turbulence , Lecture Notes in Physics (Springer Berlin Heidelberg, 2011).
- [46] V. Zakharov, V. Lvov, and G. Falkovich, Kolmogorov Spectra of Turbulence I: Wave Turbulence , (Springer Berlin Heidelberg, 1992).
- [47] D. D. Sheka, I. A. Yastremsky, B. A. Ivanov, G. M. Wysin, and F. G. Mertens, Phys. Rev. B 69 , 054429 (2004).
- [48] Y.-Q. Li, Y.-H. Liu, and Y. Zhou, Phys. Rev. B 84 , 205123 (2011).
- [49] K. Fujimoto and M. Tsubota, Phys. Rev. A 93 , 033620 (2016).
- [50] A. J. Bray, Adv. Phys. 43 , 357 (1994).
- [51] Calabrese and Gambassi, J. Phys. A 38 , R133 (2005).
- [52] M. Babadi, E. Demler, and M. Knap, Phys. Rev. X 5 , 041005 (2015).
- [53] S. M. Davidson and A. Polkovnikov, Phys. Rev. Lett. 114 , 045701 (2015).
- [54] J. Schachenmayer, A. Pikovski, and A. M. Rey, Phys. Rev. X 5 , 011022 (2015).
- [55] B. Zhu, A. M. Rey, and J. Schachenmayer, arXiv:1905.08782.
Saraswat Bhattacharyya Joaquin F. Rodriguez-Nieva Eugene Demler
Universal prethermal dynamics in Heisenberg ferromagnets
We study the universal far from equilibrium dynamics of magnons in Heisenberg
ferromagnets. We show that such systems exhibit universal scaling in momentum
and time of the quasiparticle distribution function, with the universal
exponents distinct from those recently observed in Bose-Einstein condensates.
This new universality class originates from the SU(2) symmetry of the
Hamiltonian, which leads to a strong momentum-dependent magnon-magnon
scattering amplitude. We compute the universal exponents using the Boltzmann
kinetic equation and incoherent initial conditions that can be realized with
microwave pumping of magnons. We compare our numerical results with analytic
estimates of the scaling exponents and demonstrate the robustness of the
scaling to variations in the initial conditions. Our predictions can be tested
in quench experiments of spin systems in optical lattices and pump-probe
experiments in ferromagnetic insulators such as yttrium iron garnet.