Controlling acoustic waves using magneto-elastic Fano resonances Controlling acoustic waves using magneto-elastic Fano resonances O. S. Latcham,1 Y. I. Gusieva,2 A. V. Shytov,1 O. Y. Gorobets,2 and V. V. Kruglyak1, a) 1)University of Exeter, Stocker Road, Exeter, EX4 4QL, United Kingdom 2)Igor Sikorsky Kyiv Polytechnic Institute, 37 Prosp. Peremohy, Kyiv, 03056, Ukraine (Dated: 5 May 2020) We propose and analyze theoretically a class of energy-efficient magneto-elastic devices for analogue signal processing. The signals are carried by transverse acoustic waves while the bias magnetic field controls their scattering from a magneto-elastic slab. By tuning the bias field, one can alter the resonant frequency at which the propagating acoustic waves hybridize with the magnetic modes, and thereby control transmission and reflection coefficients of the acoustic waves. The scattering coefficients exhibit Breit-Wigner/Fano resonant behaviour akin to inelastic scattering in atomic and nuclear physics. Employing oblique incidence geometry, one can effectively enhance the strength of magneto- elastic coupling, and thus countermand the magnetic losses due to the Gilbert damping. We apply our theory to discuss potential benefits and issues in realistic systems and suggest routes to enhance performance of the proposed devices.

Optical and, more generally, wave-based computing paradigms gain momentum on a promise to replace and com- HB plement the traditional semiconductor-based technology.1 The I energy savings inherent to non-volatile memory devices has spurred the rapid growth of research in magnonics,2,3 in which spin waves4 are exploited as a signal or data carrier. Yet, the θ Tω progress is hampered by the magnetic loss (damping).5,6 In- deed, the propagation distance of spin waves is rather short in Rω M ferromagnetic metals while low-damping magnetic insulators are more difficult to structure into nanoscale devices. In con- trast, the propagation distance of acoustic waves is typically much longer than that of spin waves at the same frequencies.7 Hence, their use as the signal or data carrier could reduce the propagation loss to a tolerable level. Notably, one could con- NM M NM trol the acoustic waves using a magnetic field by coupling them to spin waves within magnetostrictive materials.8–10 To FIG. 1: The prototypical magneto-elastic resonator is a thin minimize the magnetic loss, the size of such magneto-acoustic magnetic slab (M) of width δ, biased by an external field HB, functional elements should be kept minimal. This implies and embedded into a non-magnetic (NM) matrix. The coupling propagating acoustic waves to confined spin wave acoustic wave with amplitude I incident at angle θ induces modes of finite-sized magnetic elements. As we show below precession of the magnetisation vector M via the this design idea opens a route towards hybrid devices combin- magneto-elastic coupling. As a result, the wave is partly 2,3 ing functional benefits of magnonics with the energy effi- transmitted and reflected, with respective amplitudes Tω 7,11,12 ciency of phononics. and Rω . The phenomena resulting from interaction between coher- ent spin and acoustic waves have already been addressed in the research literature: the spin wave excitation of prop- acoustic devices in which the signal is carried by acoustic agating acoustic waves7,13–15 and vice versa,8,16–18 acous- waves while the magnetic field controls its propagation via tic parametric pumping of spin waves,19–21 magnon-phonon the magnetoelastic interaction in thin isolated magnetic in- coupling in cavities22–24 and mode locking,25 magnonic- clusions as shown in Fig. 1. By changing the applied mag- phononic crystals,26,27 Bragg scattering of spin waves from a netic field, one can alter the frequency at which the incident surface acoustic wave induced grating,28–30 topological prop- acoustic waves hybridize with the magnetic modes of the in- erties of magneto-elastic excitations,15,31 acoustically driven clusions. Thereby, one can control the acoustic waves by the arXiv:1906.07297v2 [physics.app-ph] 4 May 2020 spin pumping and spin Seebeck effect,32,33 and optical excita- resonant behaviour of Breit-Wigner and Fano resonances in tion and detection of magneto-acoustic waves.34–40 However, the magnetic inclusion.41 We find that the strength of the res- studies of the interaction between propagating acoustic waves onances is suppressed by the ubiquitous magnetic damping and spin wave modes of finite-sized magnetic elements, which in realistic materials, but this can be mitigated by employing are the most promising for applications, have been relatively oblique incidence geometry. To compare magneto-acoustic scarce to date.10,34,36,39 materials for such devices, we introduce a figure of merit. The Here, we explore theoretically the class of magneto- magneto-elastic Fano resonance is identified as most promis- ing in terms of frequency and field tuneability. To enhance res- onant behaviour, we explore the oblique incidence as a means by which to enhance the figure of merit. a)Electronic mail: [email protected] We consider the simplest geometry in which magneto- Controlling acoustic waves using magneto-elastic Fano resonances 2 elastic coupling can affect sound propagation. A ferromag- For the magnetization precession in the magnetic layer driven netic slab ("magnetic inclusion") of thickness δ, of the or- by the acoustic wave, we thus obtain der of 10 nm, is embedded within a non-magnetic medium γB(ωk + iω k ) (Fig.1). The slab is infinite in the Y − Z plane, has satu- ω,y ey ω,x mx = 2 U, (6) ω − ωxωy ration magnetization Ms, and is biased by the applied field e e HB = HBzˆ. Due to the magneto-elastic coupling, this equilib- iγB(ωexkω,y + iωkω,x) my = 2 U, (7) rium configuration is perturbed by shear stresses in the xz- and ω − ωexωey yz planes associated with the incident acoustic wave. where we have denoted = (H + N M ) and To derive the equations of motion, we represent the mag- ωx(y) γµ0 B x(y) s netic energy density F of the magnetic material as a sum of the ωex(y) = ωx(y) − iωα. The complex-valued wave number kω,x 42 is given by the dispersion relation magneto-elastic FME and purely magnetic FM contributions. Taking into account the Zeeman and demagnetizing energies, ρ 2 2
2  2 γB2  µ 2 2 ω ω − ωxωy − k ω − ωxωy + ωx we write F = −µ H M + 0 (N M + N M ), where N C e e ω,y e e MsC e M 0 B 2 x x y y x(y) k2 = , ω,x h 2 i are the demagnetising coefficients, Nx + Ny = 1, M is the ω2 − ω ω + γB ω ex ey MsC ey magnetization and µ0 is the magnetic permeability. In a crys- tal of cubic symmetry, the magnetoelastic contribution takes (8) the form43 where kω,y is equal to that of the incident wave, and the branch with Imkω,x > 0 describes a forward wave decaying 0 B B 2 into the slab. Eq. (8) describes the hybridization between FME = 2 ∑ MiMjui j + 2 ∑Mi uii, i, j = x,y,z, (1) acoustic waves and magnetic precession at frequencies close Ms i6= j Ms i to ferromagnetic resonance (FMR) at frequency ωFMR, with linewidth Γ . The frequency at which the precession am- B0 B FMR where and are the linear isotropic and anisotropic plitudes (Eqs. (6) and (7)) diverge is given by the condition magneto-elastic coupling constants, respectively.44 The strain (ω + iΓ /2)2 = ω ω . In the limit of small α, this tensor is u = 1 (∂ U + ∂ U ), where U are the displace- FMR FMR ex ey jk 2 j k k j j yields ω = ω ω and Γ = α(ω + ω ). Away from ment vector components. To maximize the effect of the cou- FMR x y FMR x y the resonance, Eq. (8) gives the linear dispersion of acous- pling B, we consider a transverse acoustic plane wave incident tic waves. In the non-magnetic medium (B = 0), one finds on the slab from the left and polarized along the bias field, so 2 2 k0 = ω ρ0/C0. Here and below, the subscript ’0’ is used to that Ux = Uy = 0, Uz = U(x,y,t). The non-vanishing com- 1 1 mark quantities pertaining to the non-magnetic matrix. ponents of the strain tensor are uxz = ∂xU and uyz = ∂yU, 2 2 To calculate the reflection and transmission coefficients, Rω and FME is linear in both M and U: and Tω , for a magnetic inclusion, we introduce the mechanical B impedance as Z = iσxz/ωUω . Solution of the wave matching FME = (Mxuxz + Myuyz). (2) problem can then be expressed via the ratio of load (ZME) and Ms source (Z0) impedances. For impedances in the forward (F) The magnetization dynamics in the slab is due to the effec- and backward (B) directions in the magnetic slab, we find tive magnetic field, µ H = −δF/δM. We define m as  k ,y  0 eff 2 ω Ck γB ωey ∓ iω k the small perturbation of the magnetic order, i.e. |m|  Ms. Z(F/B) = ω,x 1 + ω,x . (9) 4 ω,ME  2  Linearizing the Landau-Lifshitz-Gilbert equation, we write ω CMs ω − ωexωey

∂mx ∂U ∂my Here, the ‘-’ and ‘+’ signs correspond to (F) and (B), re- − = γµ (H + NyM )my + γB + α , (3) ∂t 0 B s ∂y ∂t spectively. For the non-magnetic material, Eq. (9) recov- 45 √ ∂my ∂U ∂mx ers the usual acoustic impedance Z0 = cosθ ρ0C0. Due = γµ0 (HB + NxMs)mx + γB + α , (4) (F/B) ∂t ∂x ∂t to magnon-phonon hybridization, ReZω,ME diverges at ωFMR and vanishes at a nearby frequency ωME. For α = 0, the latter where γ is the gyromagnetic ratio and α is the Gilbert is given by damping constant. To describe the acoustic wave, we in- s (ME) 2 clude the magneto-elastic contribution to the stress, σ jk = γB ωME = ωxωy − ωy. (10) δFME /δu jk, into the momentum balance equation: MsC

2     Reflection Rω and transmission Tω coefficients are then ∂ U ∂ ∂U B ∂ ∂U B 45 ρ 2 = C + mx + C + my , found via the well-known relations as ∂t ∂x ∂x Ms ∂y ∂y Ms (5) (η + 1)(1 − η )sin(k δ) R = eω ω ω,x , C = c ω where 44 is the shear modulus and ρ is the mass density. (ηeω ηω + 1)sin(kω,xδ) + i(ηω + ηeω )cos(kω,xδ) The non-magnetic medium is described by Eq.(5) with B = 0. (11) Since the values of C, B, and Nx,y are constant within each i(ηω + ηeω ) individual material, we shall seek solutions of the equations in Tω = , (ηω ηω + 1)sin(kω,xδ) + i(ηω + ηω )cos(kω,xδ) the form of plane waves U,mx(y) ∝ exp[i(kω,xx +kω,yy−ωt)]. e e From herein, we consider all variables in the Fourier domain. (12) Controlling acoustic waves using magneto-elastic Fano resonances 3 where δ is the thickness of the magnetic inclusion, ηω = acoustic waves are damped in our model only due to the cou- (F) (B) 46 pling with spin waves. ZME/Z0 and ηeω = ZME/Z0. In close proximity to the res- onance, the impedances changes rapidly. Expanding Eq. (11) To consider how the magneto-elastic resonance is affected −3 near ωME in the limit kω δ  1, we obtain by the damping, we also plot the response for α of 10 and 10−2, red and blue curves in Fig.2, respectively. An increase of α from 10−4 to 10−3 significantly suppresses and broad- iΓR/2 iφ ens the resonant peak. For a more common, realistic value of Rω = e + R0, (13) −2 (ω − ωME) + iΓR/2 10 the resonance is quenched entirely. A stronger magne-   toelastic coupling (i.e. high values of B) could, in principle, C rωx φ = −2 arctan tanθ , countermand this suppression. This, however, is also likely C ω 0 y to enhance the phonon contribution to the magnetic damping, leading to a correlation between B and α observed in realistic where R0 represents a smooth non-resonant contribution due 49 to elastic mismatch at the interfaces, while φ represents a res- magnetic materials. onant phase, which is non-zero for finite θ and approaches π To characterise the strength of the Fano resonance, we note rapidly. In a system with no magnetic damping, the hybridiza- that the fate of the magnon excited by the incident acoustic wave is decided by the relation between the emission rate ΓR, tion yields a resonance of finite linewidth ΓR, see Eq. (14), and absorption rate ΓFMR. Hence, we introduce 2  2  the respective figure of merit as ϒ = ΓR/ΓFMR. This quan- γB p 2 C 2 Γ = ρ C ωycos θ + ωxsin θ δ. tity depends upon the material parameters, device geometry, R 2M C2 cosθ 0 0 C2 s 0 and bias field. As seen from the first terms on the l.h.s. of (14) Eqs. (6) and (7), the relation between the dynamic magnetisa- The origin of this linewidth can be explained as follows. Due tion components mx,y are determined by the quantities ωx and to the magneto-elastic coupling incident propagating acoustic p modes can be converted into localised magnon modes. These ωy. Equating these terms, one finds mx ∝ my ωy/ωx, i.e. the 50 modes in turn either decay due to the Gilbert damping or are precession of m is highly elliptical, due to the demagnetis- re-emitted as phonons. The rates of these transitions are pro- ing field along x. This negatively affects the phonon-magnon coupling for normal incidence (k = 0): the acoustic wave portional to ΓFMR and ΓR, respectively, and the total decay y couples only to m , as given by the second term in Eqs. (6) rate is Γ = ΓR + ΓFMR. This is similar to resonant scattering x 47 and (7). One way to mitigate this is to increase H , mov- in quantum theory , such that ΓR and ΓFMR are analogous to B ing the ratio ω /ω closer to 1 and thus improving the figure the the elastic (Γe) and inelastic (Γi) linewidths respectively. y x of merit. To compare different magneto-elastic materials, the When α = 0, ΓFMR vanishes, and Γ = ΓR. Acoustic waves in the geometry of Fig. 1 can be scattered dependence on the layer thickness δ and elastic properties of via several channels. E.g. in a non-magnetic system (B = 0), the non-magnetic matrix (i.e. ρ0 and C0) can be eliminated by elastic mismatch can yield Fabry-Pérot resonance due to the calculating a ratio of the figures of merit for the compared ma- quarter wavelength matching of δ and the acoustic wave- terials. The comparison can be performed either at the same length. However, this occurs at very high frequencies, which value of the bias field, or at the same operating frequency. The we do not consider here. To understand the resonant magneto- latter situation is more appropriate for a device application, elastic response, it is instructive to consider first the case of but to avoid unphysical parameters, we present our results for normal incidence (θ = 0), when the demagnetising energy the same µ0HB. An example of such comparisons for yttrium takes a simplified form due to the lack of immediate inter- iron garnet (YIG), cobalt (Co) and permalloy (Py) is offered in Table I. faces to form surface poles in y the direction, so that Nx = 1 Another way to improve ϒ is to employ the oblique inci- and Ny = 0. Including magneto-elastic coupling (B 6= 0), we dence (θ 6= 0), in which the acoustic mode is also coupled to plot the frequency dependence of Rω and Tω using Eq. (11) and (12) in Fig.2. To gain a quantitative insight, we analysed the magnetisation component my. The latter is not suppressed a magnetic inclusion made of cobalt (ρ = 8900kgm−3, B = by the demagnetisation effects if Ny  1. The resulting en- 10MPa, C = 80GPa, γ = 176GHzT−1, M = 1MAm−1), em- hancement in ϒ is reflected in the full equation by the inclu- −3 sion of ω and ω from Γ , bedded into a non-magnetic matrix (ρ0 = 3192kgm ,C0 = x y R 298GPa). To highlight the resonant behaviour, we first sup-  2  −4 2 C 2 press α to 10 . The reflection coefficient exhibits an asym- 2 HBcos θ + 2 Mssin θ ΓR γδB p C0 metric non-monotonic dependence, shown as a black curve in ϒ = = ρ0C0 2 2 , (15) ΓFMR 2 αC M cosθ Fig.2(a), characteristic of Fano resonance.27,41 This line shape s can be attributed to coupling between the discrete FMR mode where ωx  ωy and HB  Ms is assumed. For small θ, the of the magnetic inclusion and the continuum of propagating approximation Nx ' 1 and Ny ' 0 still holds. As a result, non- acoustic modes in the surrounding non-magnetic material.41 zero θ increases peak reflectivity, as seen in Fig.3. The evolu- If the two materials had matching elastic properties, Rω would tion of the curves in Fig.3 with θ is explained by the variation exhibit a symmetric Breit-Wigner lineshape.47 The transmis- of the phase φ of the resonant scattering relative to that of the sion shown in Fig.2(b) exhibits an approximately symmetric non-resonant contribution R0. The latter changes its sign at in- 48 2 2 ◦ dip near the resonance. The absorbance |Aω | = 1−|Rω | − cidence angle of about 30 , which yields a nearly symmetric 2 |Tω | , shown in Fig.2(c) exhibits a symmetric peak, since the curve (blue), and an inverted Fano resonance at larger angles Controlling acoustic waves using magneto-elastic Fano resonances 4

0.25 (a) 1.00 (c) Damping, α : (b) 4 0.20 10− 0.3 3 0.95 10− 2 |

2 | | ) 10− ) 0.15 )

f 0.90 f f

( 0.2 ( ( T R | A | 0.10 0.85 | 0.1 0.05 0.80

0.00 0.75 0.0 7.10 7.12 7.14 7.16 7.18 7.10 7.12 7.14 7.16 7.18 7.10 7.12 7.14 7.16 7.18 Frequency, f (GHz)

FIG. 2: The frequency dependence of the absolute values of (a) reflection and (b) transmission coefficients and (c) absorbance is shown for a 20nm thick magnetic inclusion. The vertical dashed and solid black lines represent the ferromagnetic resonance frequency ωFMR and magneto-elastic resonance frequency ωME respectively. The non-magnetic and magnetic materials are assumed to be silicon nitride and cobalt, respectively, with parameters given in the text. The bias field is µ0HB = 50mT, which leads to fME ≈ 7.138 GHz.

(green). Although larger incidence angles may be hard to im- 0.06 plement in a practical device, the resonant scattering is still 0.25 1 ΓFMR − ) enhanced at smaller angles. s Υ

n , ( t i

Above, we have focused on the simplest geometry that ad- r 0.20 e 0.04 Υ R M mits full analytic treatment. To implement our idea exper- M

f F o / imentally, particular care should be taken about the acous- )

e 0.15 1 r − tic waves polarization and propagation direction relative to u 0 g

i 0.02 1 F the direction of the magnetization. Indeed, our choice max- ΓR ( 0.10 R imises magnetoelastic response. If however, the polariza- Γ tion is orthogonal to the bias field H , i.e. U = 0, the cou- B z 0.00 0.05 pling would be second-order in magnetization components 0 10 20 30 40 Angle, θ (deg.) mx,y, and would not contribute to the linearized LLG equation. Furthermore, we have neglected the exchange and magneto- dipolar fields that could arise due to the non-uniformity of the FIG. 4: Figure of merit ϒ and radiative linewidth ΓR are both > ◦ magnetization. To assess the accuracy of this approximation, enhanced in the oblique incidence geometry (θ 0 ). we note that the length scale of this non-uniformity is set by Ferromagnetic linewidth ΓFMR remains unchanged. Co is assumed with = 10−3. the acoustic wavelength λ, of about 420nm for our parame- α ters rather than by the magnetic slab thickness δ. The asso-

2 ciated exchange field is µ0Ms(klex) ' 9mT. The k-dependent contributions to the magneto-dipole field vanish at normal in- 0.06

45◦ 0.05 0◦ TABLE I: Comparison of the figure of merit ϒ for different

| 15◦ materials, assuming δ = 20nm, µ0HB = 50mT and

) 0.04 f

( C0 = 298GPa. R | 0.03 Parameters YIG Co Py 0.02 ϒ(θ = 0◦) 4.3x10−2 1.7x10−3 2.7x10−4 Γ (ns−1) 1.9x10−4 7.5x10−3 2.0x10−4 30◦ R −1 −3 0.01 ΓFMR (ns ) 4.4x10 4.3 0.74 ◦ −2 −3 −4 7.00 7.05 7.10 7.15 7.20 7.25 ϒ(θ = 30 ) 4.1x10 2.5x10 2.8x10 −1 −4 −2 −4 Frequency, f (GHz) ΓR (ns ) 1.8x10 1.1x10 2.1x10 −1 −3 ΓFMR (ns ) 4.4x10 4.3 0.74 FIG. 3: Peak R( f ) is enhanced and slightly shifted to the left fME = ωME/2π (GHz) 2.97 7.14 6.26 in the oblique incidence geometry (θ > 0◦). Coloured curves B (MJm−3) 0.55 10 -0.9 represent specific incidence angles sweeping from 0◦ to 45◦. C (GPa) 74 80 50 −3 Moderate Gilbert damping of α = 10−3 is assumed. The ρ (kgm ) 5170 8900 8720 −4 −2 −3 dashed vertical line corresponds to the magnetoelastic α 0.9x10 1.8x10 4.0x10 −1 resonance frequency. Ms (kAm ) 140 1000 760 Controlling acoustic waves using magneto-elastic Fano resonances 5 cidence but may become significant at oblique incidence, giv- 17P. G. Gowtham, T. Moriyama, D. C. Ralph, and R. A. Buhrman, “Traveling ◦ ing µ0Mskyδ ' 98mT at θ = 15 . In principle, these could surface spin-wave resonance spectroscopy using surface acoustic waves,” J. increase the resonant frequency of the slab by a few GHz but Appl. Phys. 118, 233910 (2015). 18H. Ulrichs, D. Meyer, F. Döring, C. Eberl, and H. U. Krebs, “Spectral would complicate the theory significantly. The detailed analy- control of elastic dynamics in metallic nano-cavities,” Sci. Rep. 7, 10600 sis of the associated effects is beyond the scope of this report. (2017). In summary, we have demonstrated that the coupling be- 19A. G. Gurevich, “Parametric amplification of magnetic waves in ferrites by tween the magnetisation and strain fields can be used to con- an elastic wave,” Phys. Sol. State 6, 1885 (1965). 20H. Keshtgar, M. Zareyan, and G. E. W. Bauer, “Acoustic parametric pump- trol acoustic waves by magnetic inclusions. We show that ing of spin wave,” Sol. State. Commun. 198, 30 (2014). the frequency dependence of the waves’ reflection coefficient 21P. Chowdhury, P. Dhagat, and A. Jander, “Parametric amplification of spin from the inclusions has a Fano-like lineshape, which is partic- waves using bulk acoustic waves,” IEEE Trans. Magn. 51, 1300904 (2015). ularly sensitive to the magnetic damping. Figure of merit is 22A. N. Litvinenko, A. V. Sadovnikov, V. V. Tikhonov, and S. A. Nikitov, introduced to compare magnetoelastic materials and to char- “Brillouin light scattering spectroscopy of magneto-acoustic resonances in a thin-film garnet resonator,” IEEE Magn. Lett. 6, 3200204 (2015). acterize device performance. In particular, the figure of merit 23X. F. Zhang, C. L. Zou, L. Jiang, and H. X. Tang, “Cavity magnomechan- is significantly enhanced for oblique incidence of acoustic ics,” Sci. Adv. 2, e1501286 (2016). waves, which enhances their coupling to the magnetic modes. 24C. Kong, B. Wang, Z. X. Liu, H. Xiong, and Y. Wu, “Magnetically control- We envision that further routes may be taken to transform lable slow light based on magnetostrictive forces,” Opt. Express 27, 5544 our prototype designs into working devices, such as forming (2019). 25S. Wang and T. lin Hsu, “Observation of spin-wave spectrum in an instabil- a magneto-acoustic metamaterial to take advantage of spatial ity and mode-locking experiment,” Appl. Phys. Lett 16, 534 (1970). resonance. 26S. Nikitov, Y. Filimonov, S. Vysotsky, Y. Khivintsev, and E. Pavlov, “Yt- The research leading to these results has received funding trium iron garnet based phononic-magnonic crystal,” Proc. of 2012 IEEE from the Engineering and Physical Sciences Research Coun- International Ultrasonics Symposium , 1240–1243 (2012). 27P. Graczyk, J. Kłos, and M. Krawczyk, “Broadband magnetoelastic cou- cil of the United Kingdom (Grant No. EP/L015331/1) and pling in magnonic-phononic crystals for high-frequency nanoscale spin- from the European Union’s Horizon 2020 research and inno- wave generation,” Phys. Rev. B 95, 104425 (2017). vation program under Marie Skłodowska-Curie Grant Agree- 28A. V. Chumak, P. Dhagat, A. Jander, A. A. Serga, and B. Hillebrands, ment No. 644348 (MagIC). “Reverse doppler effect of magnons with negative group velocity scattered from a moving bragg grating,” Phys. Rev. B 81, 140404 (2010). 29R. G. Kryshtal and A. V. Medved, “Influence of magnetic anisotropy on 1 D. G. Feitelson, Optical computing: A survey for computer scientists (MIT dynamic magnonic crystals created by surface acoustic waves in yttrium Press, 1988). iron garnet films,” J. Magn. Magn. Mater. 426, 666 (2017). 2 V. V. Kruglyak, S. O. Demokritov, and D. Grundler, “Magnonics,” J. Phys. 30R. G. Kryshtal and A. V. Medved, “Nonlinear spin waves in dynamic D: Appl. Phys 43, 264001 (2010). magnonic crystals created by surface acoustic waves in yttrium iron gar- 3 S. Nikitov, D. Kalyabin, I. Lisenkov, A. Slavin, Y. Barabanenkov, S. Os- net films,” J. Phys. D: Appl. Phys. 50, 495004 (2017). okin, A. Sadovnikov, E. Beginin, M. Morozova, Y. Sharaevsky, Y. Fil- 31R. Takahashi and N. Nagaosa, “Berry curvature in magnon-phonon hybrid imonov, Y. Khivintsev, S. Vysotsky, V. Sakharov, and E. Pavlov, “Magnon- systems,” Sci. Rep. 117, 217205 (2016). ics: a new research area in spintronics and spin wave electronics,” Phys. 32K. Uchida, H. Adachi, T. An, T. Ota, M. Toda, B. Hillebrands, S. Maekawa, Uspekhi 58, 1002 (2015). and E. Saitoh, “Long-range spin seebeck effect and acoustic spin pumping,” 4 A. I. Akhiezer, V. G. B. yakhtar, and S. V. Peletminskii, Spin waves (North- Nature Mater. 10, 737 (2011). Holland, Amsterdam, 1968). 33N. I. Polzikova, S. G. Alekseev, V. A. Luzanov, and A. O. Raevskiy, “Elec- 5 V. N. Krivoruchko, “Spin waves damping in nanometre-scale magnetic ma- troacoustic excitation of spin waves and their detection due to the inverse terials (review article),” Low Temp. Phys. 41, 670 (2015). spin hall effect,” Phys. Solid State 60, 2211 (2018). 6 S. Azzawi, A. T. Hindmarch, and D. Atkinson, “Magnetic damping phe- 34Y. Yahagi, B. Harteneck, S. Cabrini, and H. Schmidt, “Controlling nano- nomena in ferromagnetic thin-films and multilayers,” J. Phys. D: Appl. Phys magnet magnetization dynamics via magnetoelastic coupling,” Phys. Rev. 50, 471001 (2017). B 90, 140405(R) (2014). 7 J. H. Collins, “Short history of microwave acoustics,” IEEE Trans. Mi- 35V. Kats, T. Linnik, A. Salasyuk, A. Rushforth, M. Wang, P. Wadley, A. Aki- crowave Theory Tech. MTT-32, 1127 (1984). mov, S. Cavill, V. Holy, A. Kalashnikova, and A. Scherbakov, “Ultra- 8 C. Kittel, “Interaction of spin waves and ultrasonic waves in ferromagnetic fast changes of magnetic anisotropy driven by laser-generated coherent and crystals,” Phys. Rev. 110, 836 (1958). noncoherent phonons in metallic films,” Phys. Rev. B 93, 214422 (2016). 9 H. Bömmel and K. Dransfeld, “Excitation of hypersonic waves by ferro- 36C. Berk, Y. Yahagi, S. Dhuey, S. Cabrini, and H. Schmidt, “Controlling the magnetic resonance,” Phys. Rev. Lett 3, 83 (1959). influence of elastic eigenmodes on nanomagnet dynamics through pattern 10 L. Dreher, M. Weiler, M. Pernpeintner, H. Huebl, R. Gross, M. Brandt, and geometry,” J. Magn. Magn. Mater 426, 239 (2017). S. Goennenwein, “Surface acoustic wave driven ferromagnetic resonance in 37H. Yang, F. Garcia-Sanchez, X. Hu, S. Sievers, T. Böhnert, J. D. Costa, nickel thin films: Theory and experiment,” Phys. Rev. B 86, 134415 (2012). M. Tarequzzaman, R. Ferreira, M. Bieler, , and H. W. Schumacher, “Exci- 11 N. Li, J. Ren, L. Wang, G. Zhang, P. Hanggi, and B. W. Li, “Colloquium: tation and coherent control of magnetization dynamics in magnetic tunnel Phononics: Manipulating heat flow with electronic analogs and beyond,” junctions using acoustic pulses,” Appl. Phys. Lett 113, 072403 (2018). Rev. Mod. Phys. 84, 1045 (2012). 38M. Deb, M. H. E. Popova, N. Keller, S. Mangin, and G. Malinowski, “Pi- 12 M. Maldovan, “Sound and heat revolutions in phononics,” Nature 503, 209 cosecond acoustic-excitation-driven ultrafast magnetization dynamics in di- (2013). electric bi-substituted yttrium iron garnet,” Phys. Rev. B 98, 174407 (2018). 13 R. B. Höllander, C. Müller, J. Schmalz, M. Gerken, and J. McCord, “Mag- 39S. Mondal, M. Abeed, K. Dutta, A. De, S. Sahoo, A. Barman, and netic domain walls as broadband spin wave and elastic magnetisation wave S. Bandyopadhyay, “Hybrid magnetodynamical modes in a single magne- emitters,” Sci. Rep 8, 13871 (2018). tostrictive nanomagnet on a piezoelectric substrate arising from magnetoe- 14 S. Streib, H. Keshtgar, and G. E. W. Bauer, “Damping of magnetization lastic modulation of precessional dynamics,” ACS Appl. Mater. Interfaces dynamics by phonon pumping,” Phys. Rev. Lett 121, 027202 (2018). 10, 43970 (2018). 15 E. Thingstad, A. Kamra, A. Brataas, and A. Sudbø, “Chiral phonon trans- 40Y. Hashimoto, D. Bossini, T. H. Johansen, E. Saitoh, A. Kirilyuk, and port induced by topological magnons,” Phys. Rev. Lett 122, 107201 (2019). T. Rasing, “Frequency and wavenumber selective excitation of spin waves 16 X. Li, D. Labanowski, S. Salahuddin, and C. S. Lynch, “Spin wave gener- through coherent energy transfer from elastic waves,” Phys. Rev. B 97, ation by surface acoustic waves,” J. Appl. Phys. 122, 043904 (2017). Controlling acoustic waves using magneto-elastic Fano resonances 6

140404 (2018). 46M. Born and E. Wolf, Principles of Optics (Oxford, New York: Pergamon, 41M. F. Limonov, M. V. Rybin, A. N. Poddubny, and Y. S. Kivshar, “Fano 1964). resonances in photonics,” Nature Photon 11, 543 (2017). 47L. Landau and E. M. Lifshitz, Quantum Mechanics (Oxford, New York: 42R. L. Comstock and B. A. Auld, “Parametric coupling of the magnetization Pergamon, 1965) Chap. 145. and strain in a ferrimagnet. i. parametric excitation of magnetostatic and 48S. Klaiman, “Reflections on one dimensionl transmission,” Chem. Phys. elastic modes,” J. Appl. Phys. 34, 1461 (1963). 482, 277 (2017). 43A. Kamra, H. Keshtgar, P. Yan, and G. E. W. Bauer, “Coherent elastic 49S. Emori, B. Gray, H.-M. Jeon, J. Peoples, M. Schmitt, K. Mahalingam, excitation of spin waves,” Phys. Rev. B 91, 104409 (2015). M. Hill, M. McConney, M. Gray, U. Alaan, A. Bornstein, P. Shafer, 44E. Callen and H. B. Callen, “Magnetostriction, forced magnetostriction, A. N’Diaye, E. Arenholz, G. Haugstad, K.-Y. Meng, F. Yang, D. Li, S. Ma- and anomalous thermal expansion in ferromagnets,” Phys. Rev 139, A455 hat, D. Cahill, P. Dhagat, A. Jander, N. Sun, Y. Suzuki, and B. Howe, “Co- (1965). existence of low damping and strong magnetoelastic coupling in epitaxial 45L. M. Brekhovskikh and O. A. Godin, Acoustics of Layered Media (Berlin, spinel ferrite thin films,” Adv. Mater. 29, 1701130 (2017). Heidelberg: Springer, 1997). 50J. V. Kim, Spin Torque Oscillators, Vol. 63 (Academic Press, 2012) Chap. 4.