Research ArticlePHYSICS

Giant nonlinear damping in nanoscale ferromagnets

See allHide authors and affiliations

Science Advances  25 Oct 2019:
Vol. 5, no. 10, eaav6943
DOI: 10.1126/sciadv.aav6943


Magnetic damping is a key metric for emerging technologies based on magnetic nanoparticles, such as spin torque memory and high-resolution biomagnetic imaging. Despite its importance, understanding of magnetic dissipation in nanoscale ferromagnets remains elusive, and the damping is often treated as a phenomenological constant. Here, we report the discovery of a giant frequency-dependent nonlinear damping that strongly alters the response of a nanoscale ferromagnet to spin torque and microwave magnetic field. This damping mechanism originates from three-magnon scattering that is strongly enhanced by geometric confinement of magnons in the nanomagnet. We show that the giant nonlinear damping can invert the effect of spin torque on a nanomagnet, leading to an unexpected current-induced enhancement of damping by an antidamping torque. Our work advances the understanding of magnetic dynamics in nanoscale ferromagnets and spin torque devices.


Nanoscale magnetic particles are the core components of several emerging technologies such as nonvolatile spin torque memory (1), spin torque oscillators (27), targeted drug delivery, and high-resolution biomagnetic imaging (8). Control of magnetic damping holds the key to improving the performance of many nanomagnet-based practical applications. In biomagnetic characterization techniques, such as magnetic resonance imaging, relaxometry, and magnetic particle imaging, magnetic damping affects nanoparticle relaxation times and image resolution. In spin torque memory and oscillators, magnetic damping determines the electrical current necessary for magnetic switching (1) and generation of auto-oscillations (9) and thereby determines energy efficiency of these technologies. The performance of nanomagnet-based microwave detectors is also directly affected by the damping (1012). Despite its importance across multiple disciplines, magnetic damping in nanoparticles is poorly understood and is usually modeled as a phenomenological constant (6, 9).

Here, we experimentally demonstrate that a ferromagnetic nanoparticle can exhibit dynamics qualitatively different from those predicted by the harmonic oscillator model. We show that nonlinear contributions to the damping can be unusually strong, and the effective damping parameter itself can exhibit resonant dependence on frequency. Our work demonstrates that nonlinear damping in nanomagnets is qualitatively different from that in bulk ferromagnets. We show both experimentally and theoretically that this nonlinear damping originates from multimagnon scattering in a magnetic system with a discrete spectrum of magnons induced by geometric confinement.

We also discover that the nonlinear damping markedly alters the response of a nanomagnet to spin torque. Spin torque arising from injection of spin currents polarized opposite to the direction of magnetization acts as negative damping (2). We find, however, that the effect of this antidamping spin torque is reversed, leading to an enhanced dissipation due to the nonlinear magnon scattering. This counterintuitive behavior should have substantial impact on the operation of spin torque–based memory (1), oscillators (27), and microwave detectors (1012).


Spin wave spectroscopy

We study nonlinear spin wave dynamics in nanoscale elliptical magnetic tunnel junctions (MTJs) that consist of a CoFeB free layer (FL), an MgO tunnel barrier, and a synthetic antiferromagnet (SAF) pinned layer (13). Spectral properties of the FL spin wave modes are studied in a variety of MTJs with both in-plane and perpendicular-to-plane equilibrium orientations of the FL and SAF magnetization. We observe strong nonlinear damping in both the in-plane and the perpendicular MTJs, which points to the universality of the effect.

We use spin torque ferromagnetic resonance (ST-FMR) to measure magnetic damping of the FL spin wave modes. In this technique, a microwave drive current Iac sin(2π ft) applied to the MTJ excites oscillations of magnetization at the drive frequency f. The resulting magnetoresistance oscillations Rac sin (2π ft + ϕ) generate a direct voltage Vmix. Peaks in ST-FMR spectra Vmix(f) arise from resonant excitation of spin wave eigenmodes of the MTJ (1418). To improve signal-to-noise ratio, the magnitude of external magnetic field H applied parallel to the FL magnetization is modulated, and a field-derivative signal Vmix(f)=dVmix(f)/dH is measured via lock-in detection technique (13). Vmix(f) can then be obtained via numerical integration (note S1).

Figure 1A shows ST-FMR spectra Vmix(f) measured as a function of out-of-plane magnetic field H for an elliptical 52 nm by 62 nm perpendicular MTJ device (sample 1). Three spin wave eigenmodes with nearly linear frequency–field relation fn,m(H) are visible in the spectra. Micromagnetic simulations (note S2) reveal that these modes are the three lowest frequency spin wave eigenmodes of the FL (fig. S1). The lowest frequency (quasi-uniform) mode ∣0,0〉 is nodeless and has spatially uniform phase. Each of the two higher-order modes ∣n, m〉 has a single node at the FL center that is either perpendicular ∣0,1〉 or parallel ∣1,0〉 to the ellipse long axis (Fig. 1B).

Fig. 1 Spin wave spectra in a nanoscale MTJ.

(A) Normalized ST-FMR spectra V˜mix (f) of spin wave eigenmodes in a perpendicular MTJ device (sample 1) measured as a function of out-of-plane magnetic field. Resonance peaks arising from three low frequency modes of the MTJ FL ∣0,0〉, ∣0,1〉, and ∣1,0〉 are observed. (B) Schematic of the three-magnon confluence process, denoted as ψ0,1, and the inverse process of three-magnon splitting, denoted as ψ0,1*. Mode profiles are shown. (C) Spectral linewidth of the quasi-uniform ∣0,0〉 spin wave mode as a function of out-of-plane magnetic field. Strong linewidth enhancement is observed in the three-magnon scattering regime at H1 and H2.

Although the spectra shown in Fig. 1A (especially for mode ∣0,1〉) strongly resemble avoided crossing, the physics giving rise to these spectral shapes is entirely different. An avoided crossing presumes two distinct modes, which eigenfrequencies would cross as a function of an external parameter (in our case, magnetic field) in the absence of mode interaction. In the variety of investigated MTJ samples, no such modes exist. Neither spin wave modes localized in the FL nor SAF modes approach each other in energy as magnetic field is varied (Fig. 1A). Without these crossing modes, the observed lineshape irregularities of the three FL spin wave modes cannot be attributed to avoided crossing, and instead intrinsic nonlinearity of the system must be considered.

The spectral linewidth (Fig. 1C) of the resonances in Fig. 1A can be used for evaluation of the mode damping. The quasi-uniform mode ∣0,0〉 resonance visibly broadens at two magnetic field values: H1 = 0.74 kOe (4 GHz) and H2 = 1.34 kOe (6 GHz). Near H1, the mode ∣0,1〉 resonance also broadens and exhibits splitting, the same behavior is observed for the mode ∣1,0〉 at H2. At these fields, the higher-order mode frequency is twice that of the quasi-uniform mode fn,m = 2f0,0. This shows that three-magnon confluence (1924) is the mechanism of the quasi-uniform mode damping increase (Fig. 1B): Two magnons of the quasi-uniform mode ∣0,0〉 merge into a single magnon of the higher-order mode ∣n, m〉.

The most notable feature of the quasi-uniform mode resonance near H1 is its split-peak shape with a local minimum at the resonance frequency. Such a lineshape cannot be fit by the standard Lorentzian curve with symmetric and antisymmetric components (13). We therefore use a double-peak fitting function (note S1) to quantify the effective linewidth Δf0,0 of the resonance profile. For applied fields sufficiently far from H1, the ST-FMR curve recovers its single-peak shape, and Δf0,0 is determined as half width of the standard Lorentzian fitting function (13). Figure 1C shows Δf0,0 as a function of H and demonstrates a large increase of the linewidth near the fields of the three-magnon scattering regime, H1 and H2. The stepwise increase of Δf0,0 near H1 is a result of the ST-FMR curve transition between the split-peak and single-peak shapes. For fields near H2, the resonance profile broadens but does not develop a visible split-peak lineshape. As a result, Δf0,0(H) is a smooth function in the vicinity of H2.

Effect of spin torque

In MTJs, direct bias current Idc applied across the junction exerts spin torque on the FL magnetization, acting as antidamping for Idc > 0 and as positive damping for Idc < 0 (15, 25). The antidamping spin torque increases the amplitude of the FL spin wave modes (15, 26) and decreases their spectral linewidth (27). We can use spin torque from Idc to control the amplitude of spin wave eigenmodes excited in ST-FMR measurements and thereby study the crossover between linear and nonlinear regimes of spin wave resonance. Increasing ac drive through the MTJ at zero direct current can also, in principle, be used to increase the spin wave amplitude and measure nonlinear damping. However, we found that this approach is not suitable for studying the linear-to-nonlinear transition because varying ac drive affects signal-to-noise ratio and often causes dielectric breakdown of the MgO barrier.

Figure 2 shows the dependence of ST-FMR resonance curve of the ∣0,0〉 mode Vmix(f) on Idc for a 50 nm by 110 nm elliptical in-plane MTJ (sample 2). For in-plane magnetic field values far from three-magnon scattering, the amplitude of ST-FMR resonance curve Vmix(f) shown in Fig. 2A monotonically increases with increasing antidamping spin torque, as expected. At the three-magnon scattering field H = H1, the antidamping spin torque has a radically different and rather unexpected effect on the resonance curve. As illustrated in Fig. 2C, increasing antidamping spin torque first broadens the resonance at H = H1 and then transforms a single-peak resonance lineshape into a split-peak lineshape with a local minimum at the resonance frequency f = f0,0. The data in Fig. 2 demonstrate that the unusual split-peak lineshape of the resonance is only observed when (i) the three-magnon scattering of the quasi-uniform mode is allowed by the conservation of energy and (ii) the amplitude of the mode is sufficiently high, confirming that the observed effect is resonant in the frequencies of two modes undergoing three-magnon scattering and nonlinear in nature.

Fig. 2 Effect of spin torque on spin wave resonance lineshape.

(A and B) Spin wave resonance lineshapes for different values of direct bias current Idc, far from three-magnon scattering regime H > H1. (C and D) Spin wave resonance lineshapes in the three-magnon regime at H = H1. (A and C) Measured ST-FMR spectra (sample 2). (B and D) Solutions of Eqs. 3 and 4. Identical current values Idc, displayed in (A), are used in all four figure panels.

Figure 2C reveals that antidamping spin torque can increase the spectral linewidth and the effective damping of the quasi-uniform spin mode if the mode undergoes three-magnon scattering. Figure 3 further illustrates this counterintuitive effect. It shows the linewidth of the quasi-uniform mode of a 50 nm by 110 nm elliptical in-plane MTJ (sample 3) measured as a function of bias current. In Fig. 3, blue open squares show the linewidth measured at an in-plane magnetic field sufficiently far from the three-magnon scattering. At this field, the expected quasi-linear dependence of the linewidth on Idc is observed for currents well below the critical current for the excitation of auto-oscillatory magnetic dynamics. Near the critical current, the linewidth increases due to a combination of the fold-over effect (28, 29) and thermally activated switching between the large- and small-amplitude oscillatory states of the fold-over regime (15). The red open circles in Fig. 3 show the linewidth measured in the three-magnon scattering regime at H = H1. Here, in contrast, the linewidth increases with increasing ∣Idc∣ for both current polarities. Furthermore, the maximum linewidth is measured for the antidamping current polarity.

Fig. 3 Effect of spin torque on linewidth.

Linewidth of the quasi-uniform spin wave mode as a function of applied direct bias current (sample 3): red open circles—in the three-magnon scattering regime H = H1; blue open squares—far from the three-magnon scattering regime, HH1. Lines are numerical fits using Eqs. 3 and 4.

Theoretical model

Nonlinear interactions among spin wave eigenmodes of a ferromagnet give rise to a number of spectacular magnetodynamic phenomena such as Suhl instability of the uniform precession of magnetization (30, 31), spin wave self-focusing (32), and magnetic soliton formation (33, 34). In bulk ferromagnets, nonlinear interactions generally couple each spin wave eigenmode to a continuum of other modes via energy- and momentum-conserving multimagnon scattering (30). This kinematically allowed scattering limits the achievable amplitude of spin wave modes and leads to broadening of the spin wave resonance (30). It cannot explain the observed split-peak lineshape of the resonance. In nanoscale ferromagnets, geometric confinement discretizes the spin wave spectrum and thereby generally eliminates the kinematically allowed multimagnon scattering. This suppression of nonlinear scattering enables persistent excitation of spin waves with very large amplitudes (35) as observed in nanomagnet-based spin torque oscillators (2). Tunability of the spin wave spectrum by external magnetic field, however, can lead to a restoration of the energy-conserving scattering (21). To describe the nonlinear spin wave resonance in the nanoscale ferromagnet and to derive the theory of nonlinear damping, we start with a model Hamiltonian that explicitly takes into account the nonlinear scattering between the quasi-uniform mode and a higher-order spin wave mode (in reduced units with ℏ ≡ 1)H=ω0,0aa+ωn,mbb+Ψ0,02aaaa+Ψn,m2bbbb+(ψn,maab+ψn,m*aab)+ζ{exp (i ωt)a+exp (i ωt)a}(1)where a, a and b, b are the magnon creation and annihilation operators for the quasi-uniform mode ∣0,0〉 with frequency ω0,0 and for the higher-order spin wave mode ∣n,m〉 mode with frequency ωn,m, respectively. The nonlinear mode coupling term proportional to the coupling strength parameter ψn,m describes the annihilation of two ∣0,0〉 magnons and creation of one ∣n,m〉 magnon, as well as the inverse process. The Hamiltonian is written in the resonant approximation, where small nonresonant terms such as aab and aaa are neglected. The terms proportional to Ψ0,0 and Ψn,m describe the intrinsic nonlinear frequency shifts (36) of the modes ∣0,0〉 and ∣n,m〉. The last term describes the excitation of the quasi-uniform mode by an external ac drive with the amplitude ζ and frequency ω.

We further define classically a dissipation function Q, where α0,0 and αn,m are the intrinsic linear damping parameters of the modes ∣0,0〉 and ∣n,m〉 (37)Q=dadtdadt(α0,0+η0,0aa)+dbdtdbdt(αn,m+ηn,mbb)(2)

For generality, Eq. 2 includes intrinsic nonlinear damping (9) of the modes ∣0,0〉 and ∣n,m〉 described by the nonlinearity parameters η0,0 and ηn,m. However, our analysis below shows that the split-peak resonance lineshape is predicted by our theory even if η0,0 and ηn,m are set equal to zero.

Equations describing the nonlinear dynamics of the two coupled spin wave modes of the system follow from Eqs. 1 and 2 viaidadt=Ha+Q(da/dt)(3)idbdt=Hb+Q(db/dt)(4)

It can be shown (note S3) that these equations have a periodic solution a=āexp (i ωt) and b=b̄exp (i 2ωt), where ā and b̄ are the complex spin wave mode amplitudes. For this periodic solution, Eqs. 3 and 4 are reduced to a set of two nonlinear algebraic equations for absolute values of the spin wave mode amplitudes ā and b̄, which can be solved numerically. Since the ST-FMR signal is proportional to ā2 (note S1), the calculated ā2(ω) function can be directly compared to the measured ST-FMR resonance lineshape.

We use the solution of Eqs. 3 and 4 to fit the field dependence of the quasi-uniform mode linewidth in Fig. 1C. In this fitting procedure, the resonance lineshape ā2(ω) is calculated, and its spectral linewidth Δω0,0 is found numerically. The resonance frequencies ω0,0 and ωn,m are directly determined from the ST-FMR data in Fig. 1A. The intrinsic damping parameters α0,0 and αn,m near H1 and H2 are found from linear interpolations of the ST-FMR linewidths Δf0,0 and Δfn,m measured at fields far from H1 and H2. We find that Δω0,0 weakly depends on the nonlinearity parameters Ψ and η, and thus, these parameters are set to zero (fig. S2). We also find that the calculated linewidth Δω0,0 depends on the product of the drive amplitude ζ and mode coupling strength ψn,m but is nearly insensitive to the individual values of ζ and ψn,m as long as ζ · ψn,m = const (fig. S3 and note S4). Therefore, we use ζ · ψn,m as a single fitting parameter in this fitting procedure. The solid line in Fig. 1C shows the calculated field dependence of the quasi-uniform mode linewidth on magnetic field. The agreement of this single-parameter fit with the experiment is excellent.

Fig. 2 (B and D) illustrates that Eqs. 3 and 4 not only describe the field dependence of ST-FMR linewidth but also qualitatively reproduce the spectral lineshapes of the measured ST-FMR resonances and the effect of the antidamping spin torque on the lineshapes. Figure 2B shows the dependence of the calculated lineshape ā2(ω) on antidamping spin torque for a magnetic field H far from the three-magnon scattering fields. Here, the increasing antidamping spin torque induces the fold-over of the resonance curve (28) without resonance peak splitting. The dependence of ā2(ω) on antidamping spin torque for H = H1 is shown in Fig. 2D. At this field, the resonance peak in ā2(ω) first broadens with increasing antidamping spin torque and then splits, in qualitative agreement with the experimental ST-FMR data in Fig. 2C. Our calculations (fig. S2) reveal that while the nonlinearity parameters Ψ0,0, η0,0, Ψn,m, and ηn,m have little effect on the linewidth Δω0,0, they modify the lineshape of the resonance. Given that the nonlinearity parameter values are not well known for the systems studied here, we do not attempt to quantitatively fit the measured ST-FMR lineshapes.

Equations 3 and 4 also quantitatively explain the observed dependence of the quasi-uniform mode linewidth Δω0,0 on direct bias current Idc. Assuming antidamping spin torque linear in bias current (27, 38, 39), α0,0α0,0(1Idc/Ic0,0), αn,mαn,m(1Idc/Icn,m), where Icn,m>Ic0,0 are the critical currents, we fit the measured bias dependence of ST-FMR linewidth in Fig. 3 by solving Eqs. 3 and 4. The solid lines in Fig. 3 are the best numerical fits, where ζ · ψn,m and Ic are used as independent fitting parameters. The rest of the parameters in Eqs. 3 and 4 are directly determined from the experiment, following the procedure used for fitting the data in Fig. 1C. Theoretical curves in Fig. 3 capture the main feature of the data at the three-magnon scattering field H1, which is an increase of the linewidth with increasing antidamping spin torque.


Quantitative understanding of the mechanisms of the nonlinear spin wave resonance peak splitting and broadening by antidamping spin torque can be gained by neglecting the intrinsic nonlinearities Ψn,m and ηn,m of the higher-order mode ∣n,m〉. Setting Ψn,m = 0 and ηn,m = 0 in Eqs. 3 and 4 allows us to reduce the equation of motion for the quasi-uniform mode amplitude ā to the standard equation for a single-mode damped driven oscillator (notes S5 and S6), where a constant damping parameter α0,0 is replaced by an effective frequency-dependent nonlinear damping parameter α0,0effα0,0eff=α0,0+[η0,0+4αn,mψn,m2(2ωωn,m)2+4αn,m2ω2]a¯2(5)and the resonance frequency is replaced by an effective resonance frequencyω0,0eff=ω0,0+[Ψ0,0+2ψn,m2(2ωωn,m)(2ωωn,m)2+4αn,m2ω2]a¯2(6)

Equation 5 shows that the damping parameter of the quasi-uniform mode itself becomes a resonant function of the drive frequency with a maximum at half the frequency of the higher-order mode (ω=12ωn,m). The amplitude and the width of this resonant behavior of α0,0eff(ω) are determined by the intrinsic damping parameter αn,m of the higher-order mode ∣n,m〉. If αn,m is sufficiently small, the quasi-uniform mode damping is strongly enhanced at ω=12ωn,m, which leads to a decrease of the quasi-uniform mode amplitude at this drive frequency. If the drive frequency is shifted away from 12ωn,m to either higher or lower values, the damping decreases, which can result in an increase of the quasi-uniform mode amplitude ā. Therefore, the amplitude of the quasi-uniform mode ā(ω) can exhibit a local minimum at ω=12ωn,m. Because of its nonlinear origin, the tendency to form a local minimum in ā(ω) at 12ωn,m is enhanced with increasing ā. Since ā is large near the resonance frequency ω0,0, tuning ω0,0 to be equal to 12ωn,m greatly amplifies the effect of local minimum formation in ā(ω). This qualitative argument based on Eq. 5 explains the data in Fig. 2—the split-peak nonlinear resonance of the quasi-uniform mode is only observed when external magnetic field tunes the spin wave eigenmode frequencies to the three-magnon scattering condition ω0,0=12ωn,m.

Equation 6 reveals that the nonlinear frequency shift of the quasi-uniform mode also presents a resonant behavior with respect to the drive frequency. In contrast to the nonlinear damping resonance described by Eq. 5, the frequency shift is an antisymmetric function of (ω12ωn,m). The nonlinear shift is negative for ω<12ωn,m and thus causes a fold-over toward lower frequencies, while it is positive for ω>12ωn,m, causing fold-over toward higher frequencies. At the center of the resonance profile, the three-magnon process induces no frequency shift. This double-sided fold-over also contributes to the formation of the split-peak lineshape of the resonance shown in Fig. 2 (C and D) and to the linewidth broadening. As with the nonlinear damping enhancement, the antisymmetric nonlinear frequency shift and the double-sided fold-over become greatly amplified when the spin wave mode frequencies are tuned near the three-magnon scattering ω0,0=12ωn,m. Note that the mode ∣0,1〉 at the three-magnon scattering field H = H1 in Fig. 1A presents a qualitatively similar split-peak profile. When the MTJ is driven by microwave at frequencies near f = f0,1, the mode ∣0,1〉 decays into the ∣0,0〉 mode via three-magnon splitting process. Because of the discrete spectrum of magnons, the effective damping of the spin wave modes ∣n,m〉 presents a resonant character as well.

Equations 5 and 6 also shed light on the origin of the quasi-uniform mode line broadening by the antidamping spin torque. The antidamping spin torque increases the quasi-uniform mode amplitude ā via transfer of angular momentum from spin current to the mode. Since the nonlinear damping and the nonlinear frequency shift are both proportional to ā2 and both contribute to the line broadening, the antidamping spin torque can indeed give rise to the line broadening. Equation 5 reveals two competing effects of the antidamping spin torque on the quasi-uniform mode damping parameter α0,0eff: Spin torque from Idc decreases the linear component of the damping parameter α0,0α0,0(1Idc/Ic0,0) and increases the nonlinear component via increased ā2. Whether the antidamping spin torque decreases or increases the spectral linewidth of the mode depends on the system parameters. Our numerical solution of Eqs. 3 and 4 shown in Fig. 3 demonstrates that the antidamping spin torque can strongly increase the linewidth of the quasi-uniform mode when the three-magnon scattering condition ω0,0=12ωn,m is satisfied.

Here, we find that the three-magnon processes in strongly confined magnets are cardinally different from multimagnon processes known in bulk (19, 30). A prominent example of three-magnon processes in bulk systems is Suhl instability, where a magnon of uniform precession decays in two magnons of nonuniform precession. In our case, we observe an inverse process, where two uniform magnons merge into one nonuniform magnon. This process (24) does not exhibit a threshold behavior characteristic for Suhl instabilities; we find no threshold behavior upon increasing amplitude (fig. S2) or decreasing intrinsic damping. Yet, there is another fundamental difference of the three-magnon process studied here from multimagnon processes in bulk ferromagnets. In a nanomagnet, translational invariance is broken. The three-magnon processes thus require only conservation of energy, as opposed to extended systems with continuous magnon spectra where momentum must be conserved as well. The three-magnon process in an extended system is described by ω(0)=ω(k)+ω(k). The three-magnon process studied here is not allowed in extended systems, where it would require merging of two k = 0 (uniform mode) magnons into a magnon with k ≠ 0, which violates conservation of momentum. The pronounced resonant character (40) of the nonlinearity (resonant behavior with respect to the frequencies of the modes undergoing three-magnon scattering), the split-peak spectra, and the inversion of the antidamping effect are thus unique to confined systems such as nanomagnets.

We expect that the nonlinear damping discovered in this work will have strong impact on the performance of spin torque devices such as spin torque magnetic memory, spin torque nano-oscillators, and spin torque microwave detectors. Since all these devices rely on large-amplitude oscillations of magnetization driven by spin torque, the amplitude limiting resulting from the nonlinear damping is expected to have detrimental effect on the device performance.

In conclusion, our measurements demonstrate that magnetic damping of spin wave modes in a nanoscale ferromagnet has a strong nonlinear component of resonant character that appears at a discrete set of magnetic fields corresponding to three-magnon scattering. This nonlinearity can give rise to unusual spin wave resonance profiles with a local minimum at the resonance frequency in sharp contrast to the properties of the linear and nonlinear spin wave resonances in bulk ferromagnets. The nonlinearity has a profound effect on the response of the nanomagnet to spin torque. Antidamping spin torque, which reduces the quasi-uniform spin wave mode damping at magnetic fields far from the three-magnon scattering regime, can enhance the damping in the scattering regime. This inversion of the effect of spin torque on magnetization dynamics is expected to have substantial impact on the performance of nanoscale spin torque devices such as magnetic memory and spin torque oscillators.


Supplementary material for this article is available at

Note S1. Calculation of linewidth.

Note S2. Micromagnetic simulations.

Note S3. Solution of the equations of motion.

Note S4. Effect of the drive amplitude and intrinsic nonlinearitites on the resonance lineshape.

Note S5. Effective single-mode nonlinear oscillator approximation.

Note S6. Mode coupling parameter.

Fig. S1. Spatial profiles of spin wave eigenmodes.

Fig. S2. Effect of intrinsic nonlinearities on the quasi-uniform spin wave resonance lineshape.

Fig. S3. Effect of the drive amplitude on linewidth in the three-magnon regime.

References (4147)

This is an open-access article distributed under the terms of the Creative Commons Attribution-NonCommercial license, which permits use, distribution, and reproduction in any medium, so long as the resultant use is not for commercial advantage and provided the original work is properly cited.


Acknowledgments: Funding: This work was supported by the National Science Foundation through grants nos. DMR-1610146, EFMA-1641989, ECCS-1708885, and ECCS-1810541. We also acknowledge support by the Army Research Office through grant no. W911NF-16-1-0472 and Defense Threat Reduction Agency through grant no. HDTRA1-16-1-0025. A.M.G. thanks CAPES Foundation, Ministry of Education of Brazil for financial support. R.E.A. acknowledges Financiamiento Basal para Centros Cientificos y Tecnologicos de Excelencia under project FB 0807 (Chile) and Grant ICM P10-061-F by Fondo de Innovacion para la Competitividad-MINECON. B.A.I. was supported by the National Academy of Sciences of Ukraine via project no. 1/17-N and by the Program of NUST “MISIS” (grant no. K2-2017-005), implemented by a governmental decree dated 16 March 2013, no. 211. Author contributions: I.B. designed experiments, performed measurements, and carried out numerical calculations. H.K.L. and A.M.G. performed measurements. J.A.K. made the samples. Y.-J.C. and C.S. performed micromagnetic simulations. A.A.J. carried out numerical calculations. R.E.A. developed a model for the spin wave coupling parameter. B.A.I. developed the nonlinear theory. I.N.K. managed the project. All authors analyzed the data and co-wrote the paper. Competing interests: The authors declare that they have no competing interests. Data and materials availability: All data needed to evaluate the conclusions in the paper are present in the paper and/or the Supplementary Materials. Additional data related to this paper may be requested from the authors.

Stay Connected to Science Advances

Navigate This Article