France



The mutual synchronization of spin-torque nano-oscillators (STNOs) arranged in a large array improves the individual properties of each STNO, increasing the output-power and reducing the linewidth of the emitted signal. These improvements are required to implement STNOs in wireless communication schemes [1] and furthermore open applications for different oscillator based neuromorphic computing schemes [2]. Many theoretical and experimental studies considered single STNO devices [3], while experimentally the synchronization of a large number of STNOs remains a challenging task [4]. This is because the synchronized state and the phase-locking depend on the geometry, device resistance, input signals, coupling mechanism, among others [5].&lt;br/&gt;&lt;br/&gt;As a first step to guide experiments on the synchronization behavior of STNO arrays, we have undertaken a theoretical study, combining numerical and analytical approaches based on the spin-wave formalism [3]. The STNO devices considered here are of circular cross-section (100nm in diameter) with a perpendicular polarizer and an in-plane magnetized free layer under a strong out-of-plane field. This field is larger than the saturation field. The free layer magnetization is orientated out-of-plane leading to an out-of-plane precession mode that is advantageous to induce strong dipolar coupling between STNOs, which is the coupling mechanism considered here. We derived analytically the coupled equations of amplitudes and phases where the latter are similar to the Kuramoto-phase equation, but limited to a finite number of devices. Solving for the stationary states, this model provides analytical expressions for the locking-range, the stationary amplitudes, and the frequency of oscillations of the synchronized state. These solutions depend on the array geometry, number N of STNOs in the array, separation between STNOs, distribution of the DC current, as well as the different magnetic material parameters, and the external magnetic field. In this work, we investigate the solution of these equations for two specific array geometries: a straight line and a ring structure (N=3,6 and 10), shown in figure 1b.&lt;br/&gt;&lt;br/&gt;Our results show that the locking-range depends on the separation among the STNOs (because this scales the coupling strength), but it is independent of the angle among them. Furthermore, we show that the locking-range is larger when the STNOs are placed in a ring structure (or triangle for N=3) as compared to a straight-line structure. Finally, the most interesting result is that for the ring structure, we identify two dynamic modes that represent a standing wave for the STNO phases along the ring structure (observed for N&amp;gt;5). These two dynamics modes are shown in Fig. 2 and are called here (a) the in-phase mode, with identical phases for all STNOs, and (b) the out-of-phase mode, where the phase difference between two adjacent STNOs is 2π/N. We study how these states and their existence range are modified upon variation of the current density through the different STNOs (identical and non-identical currents), the STNO separation, and the number of STNOs in the ring. The analytical results are in good agreement to the numerical solutions of the Landau-Lifshitz-Gilbert (LLG) equation, as long as the separation among the STNOs is above a certain critical value. For smaller separations, the dipolar interactions become too strong and non-linear effects would require to go beyond the first-order approximations used to solve the amplitude and phase equations. The results reported here represent a first step to validate the phase equations and to define their limitations, and can be used to address now the Kuramoto equations to find solutions for large scale structures (N very large).&lt;br/&gt;&lt;br/&gt;We acknowledge financial support in Chile from FONDECYT 1200867, and Financiamiento Basal para Centros Científicos y Tecnológicos de Excelencia AFB 180001. D. M.-A. acknowledges Postdoctorado FONDECYT 2018, 3180416. M. A. Castro acknowledges Conicyt-PCHA/Doctorado Nacional/2017-21171016. The project was supported in part by ERC Advanced Grant MAGICAL No. 669204

**References:**

[1] T. Chen, et al., Proceedings of the IEEE., 104, 1919-1945(2016).&lt;br/&gt;[2] J. Grollier, et al., Nat. Electron. 3, 360 (2020).&lt;br/&gt;[3] A. Slavin and V., Tiberkevich, IEEE Trans. Magn., 45,1875(2009).&lt;br/&gt;[4] S. Tsunegi, et al., Scientific Reports., 8, 13475(2018).&lt;br/&gt;[5] R. Lebrun, et al., Nat. Commun. 8, 15825 (2017).

**Images:**

![](https://assets.underline.io/uploads/markdown_image/1/image/67e3f596fc435dd303beec991462b132.jpg)
Fig.1: (a) Schematics of an STNO device and (b) two configurations for a coupled STNO array.

![](https://assets.underline.io/uploads/markdown_image/1/image/fdc26dc3ca23f437e028887120de5d0d.jpg)
Fig. 2: Schematics of the two dynamic modes obtained for the ring structure. The snapshots show the instantaneous distribution of the x- and y- components of the magnetization in three consecutive time steps t&lt;span style=&quot;font-size:10.8333px&quot;&gt;&lt;sub&gt;1&lt;/sub&gt; &lt; t&lt;sub&gt;2 &lt;/sub&gt;&lt; t&lt;sub&gt;3&lt;/sub&gt;&lt;/span&gt;, for (a) the in-phase and (b) the out-of-phase modes.

INTERMAG 2021

 Mutual synchronization of an array of spin-torque oscillators with perpendicular polarizer

#### Welcome to the 67th Annual Conference on Magnetism and Magnetic Materials, MMM 2022!

The Conference was held from October 31 to November 4 and offers both in-person and on-demand content.

MMM 2022 is sponsored jointly by AIP Publishing and the IEEE Magnetics Society. Members of the international scientific and engineering communities interested in recent developments in fundamental and applied magnetism are invited to attend and contribute to the technical sessions. The technical program will include invited and contributed papers in oral and poster sessions, invited symposia, a tutorial, and several special sessions. This Conference provides an outstanding opportunity for worldwide participants to meet their colleagues and collaborators and discuss developments in all areas of magnetism research.

<iframe width="560" height="315" src="https://s3.amazonaws.com/pf-user-files-01/u-59356/uploads/2022-11-02/un13zp7/MMM%20Welcome%20Video%20v5.mp4" title="YouTube video player" frameborder="0" allow="accelerometer; autoplay; clipboard-write; encrypted-media; gyroscope; picture-in-picture" allowfullscreen></iframe>

For tickets to the event, please register at the following link: https://magnetism.org/registration

Please register for the MMM 2022 conference to join the event.

MMM 2022

This Conference provides an outstanding opportunity for worldwide participants to meet their colleagues and collaborators and discuss developments in all areas of magnetism research.

The mutual synchronization of spin-torque nano-oscillators (STNOs) arranged in a large array improves the individual properties of each STNO, increasing the output-power and reducing the linewidth of the emitted signal. These improvements are required to implement STNOs in wireless communication schemes [1] and furthermore open applications for different oscillator based neuromorphic computing schemes [2]. Many theoretical and experimental studies considered single STNO devices [3], while experimentally the synchronization of a large number of STNOs remains a challenging task [4]. This is because the synchronized state and the phase-locking depend on the geometry, device resistance, input signals, coupling mechanism, among others [5]. As a first step to guide experiments on the synchronization behavior of STNO arrays, we have undertaken a theoretical study, combining numerical and analytical approaches based on the spin-wave formalism [3]. The STNO devices considered here are of circular cross-section (100nm in diameter) with a perpendicular polarizer and an in-plane magnetized free layer under a strong out-of-plane field. This field is larger than the saturation field. The free layer magnetization is orientated out-of-plane leading to an out-of-plane precession mode that is advantageous to induce strong dipolar coupling between STNOs, which is the coupling mechanism considered here. We derived analytically the coupled equations of amplitudes and phases where the latter are similar to the Kuramoto-phase equation, but limited to a finite number of devices. Solving for the stationary states, this model provides analytical expressions for the locking-range, the stationary amplitudes, and the frequency of oscillations of the synchronized state. These solutions depend on the array geometry, number N of STNOs in the array, separation between STNOs, distribution of the DC current, as well as the different magnetic material parameters, and the external magnetic field. In this work, we investigate the solution of these equations for two specific array geometries: a straight line and a ring structure (N=3,6 and 10), shown in figure 1b. Our results show that the locking-range depends on the separation among the STNOs (because this scales the coupling strength), but it is independent of the angle among them. Furthermore, we show that the locking-range is larger when the STNOs are placed in a ring structure (or triangle for N=3) as compared to a straight-line structure. Finally, the most interesting result is that for the ring structure, we identify two dynamic modes that represent a standing wave for the STNO phases along the ring structure (observed for N&gt;5). These two dynamics modes are shown in Fig. 2 and are called here (a) the in-phase mode, with identical phases for all STNOs, and (b) the out-of-phase mode, where the phase difference between two adjacent STNOs is 2π/N. We study how these states and their existence range are modified upon variation of the current density through the different STNOs (identical and non-identical currents), the STNO separation, and the number of STNOs in the ring. The analytical results are in good agreement to the numerical solutions of the Landau-Lifshitz-Gilbert (LLG) equation, as long as the separation among the STNOs is above a certain critical value. For smaller separations, the dipolar interactions become too strong and non-linear effects would require to go beyond the first-order approximations used to solve the amplitude and phase equations. The results reported here represent a first step to validate the phase equations and to define their limitations, and can be used to address now the Kuramoto equations to find solutions for large scale structures (N very large). We acknowledge financial support in Chile from FONDECYT 1200867, and Financiamiento Basal para Centros Científicos y Tecnológicos de Excelencia AFB 180001. D. M.-A. acknowledges Postdoctorado FONDECYT 2018, 3180416. M. A. Castro acknowledges Conicyt-PCHA/Doctorado Nacional/2017-21171016. The project was supported in part by ERC Advanced Grant MAGICAL No. 669204

**References:**

[1] T. Chen, et al., Proceedings of the IEEE., 104, 1919-1945(2016). [2] J. Grollier, et al., Nat. Electron. 3, 360 (2020). [3] A. Slavin and V., Tiberkevich, IEEE Trans. Magn., 45,1875(2009). [4] S. Tsunegi, et al., Scientific Reports., 8, 13475(2018). [5] R. Lebrun, et al., Nat. Commun. 8, 15825 (2017).

**Images:**

![](https://assets.underline.io/uploads/markdown_image/1/image/67e3f596fc435dd303beec991462b132.jpg)
Fig.1: (a) Schematics of an STNO device and (b) two configurations for a coupled STNO array.

![](https://assets.underline.io/uploads/markdown_image/1/image/fdc26dc3ca23f437e028887120de5d0d.jpg)
Fig. 2: Schematics of the two dynamic modes obtained for the ring structure. The snapshots show the instantaneous distribution of the x- and y- components of the magnetization in three consecutive time steps t1 < t2 < t3, for (a) the in-phase and (b) the out-of-phase modes.

technical paper

INTERMAG 2021 Welcome Video



	
<figure class="video_container">
 <iframe src="https://screencast-o-matic.com/player/crflqFVnUpZ?v=6&ff=1&title=0&controls=1" width=640 height=360 frameborder="0" allowfullscreen="true"> </iframe>
</figure>

Want to buy a ticket? 

Register for INTERMAG 2021 here

[![](https://assets.underline.io/uploads/markdown_image/1/image/17ce76e7960db633258024eef8d6da7f.png)](https://intermag.org/register)

Got your ticket, but still can’t access the event?

Contact us and our team will look into it.

Magnetoelectric multiferroics, which simultaneously stabilize magnetic and polar static orders, have been predicted to present unusual dynamical magnetoelectric effects (ME). The relatively large magnitude of this dynamical ME response as well as its strong interac-tion with light is extremely attractive for potential applications among others are multifunctional spin-based memories or optical switches. The dynamical ME effects under consideration in this presentation, are characterized by the existence of an electro-active magnetic excitation called electromagnon. The electroactivity of electromagnons was first observed in type II multiferroics, where ferroelectricity is magnetically induced. In orthorhombic RMnO3 and RMn2O5(R=Tb, Y) families, both THz [1] and Raman [2] spectroscopies showed that this mode can be excited by means of the electric field of the light. The demonstration of the complementary magnetic nature of this excitation has been the prerogative of inelastic neutron scattering (INS) [3]. Sincethen, electromagnons have also been reported in Type I multiferroics hexagonal RMnO3 [4], in BiFeO3 and even in a paraelectric compound Ba2Mg2Fe12O22. To describe this unusual property, different mechanisms were proposed, in particular for the manganites families. These models basically rely on the very same magnetoelectric coupling which accounts for the static properties: lattice vibrations alter the exchange couplings or anisotropies which in turn affect the magnon spectrum. In helical magnetic systems for instance, where the multiferroicity originates from Dzyaloshinskii-Moriya Interac-tion (DMI), the same DMI has been proposed to explainthe origin of the electromagnon [5]. Although later refuted in DyMnO3, such DMI based scenario seems to be at play in CuO. Other propositions are based on Exchange-Striction (ES) [6], and couple a polarlattice mode to the magnon spectrum [1]. So far, these models have been consistent with experimental results in TbMn2O5 [3]. These interpretations tend to reduce the definition of the electromagnon to a hybrid magnon-phonon mode, or, as proposed in ErMnO3, to a mixing between a magnon and an electro-active rareearth crystal field transition [4]. In this presentation we question this issue and support the idea that electroactivity can be a purely electronic phenomenon. The RMn2O5 multiferroic compounds crystallize in the polar Pm space group [7] but present quasi-symmetries of the paraelectric Pbam space group. Along the c direction, the structure is composed of chains of Mn4+O6 octahedra, separated by layers of R3+ or Mn3+ ions. In the (a,b) plane, zig-zag chains of Mn4+O6 octahedra and Mn3+O5 square-pyramids run along the a axis and are stacked along b. GdMn2O5, like the other members of the family, shows an Incommensurate Magnetic ordering (ICM) belowTN=38 K, with a magnetic propagation wave vector k1= (½,0,0.28). Then, below Tc=32 K, the magnetic order becomes commensurate with a propagation wave vector k2= (½,0,0), accompanied by the largest electric polarization known in type II multiferroics [8]. This strong electric polarization has been ascribed to two cooperative ES mechanisms [6]. As all members of the series, the electric polarization is mainly oriented along the b axis [8] even if few studies have recently suggested the existence of small additional components along the a and even the c directions. The dynamical properties of this compound have been studied by means of IR and INS on high quality single crystals. In this presentation, we will mainly focus on the demonstration and characterization of an electromagnon present in the compound GdMn2O5. More precisely, the IR measurements allowed us to demonstrate the existence of an excitation below the TN, at an energy of 35cm-1 at 5K, which is electroactive because it is sensitive to the direction of the electric field of light (see fig. 1). The temperature tracking shows that this excitation disappears above the TN, which shows the magnetic origin of this excitation. Moreover, INS measurements have allowed us to measure the excitation in the center of the magnetic zone (½,0,0) (see fig. 2), but also to measure its dispersion, which we have been able to interpret through simulations based on the «Holstein-Primakov» representation, which shows that this excitation is an electroactive spin wave, and thus an electromagnon. Polarized INS measurements (see fig. 2) allowed us to discriminate between the nuclear and magnetic contributions, and to show that this excitation is purely magnetic, thus excluding the hypothesis of a coupling with a phonon. Finally, the presence of a crystalline field level was excluded by both INS and IR measurements, and by ab-initio calculations. Consequently, the electroactivity of this electromagnon is of a different nature than those observed in the past, and seems to have a purely electronic origin.

**References:**

a[1] A. B. Sushkov, R. V. Aguilar and S. W. Cheong, Phys. Rev. L 98, 3 (2007) [2] P. Rovillain, Cazayous and A. Sacuto, Phys. Rev. B 86, 014437 (2012) [3] S. Petit, V. Balédent and P. Foury-Leylekian, Phys. Rev. B 87, 1 (2013) [4] L. Chaix P. Roy and V. Simonet, Phys. Rev. L 112, 1 (2014) [5] H. Katsura, A. V. Balatsky and N. Nagaosa, Phys. Rev. L 98, 1 (2007) [6] V. Balédent, M. B. Lepetit, and P. Foury-Leylekian, Phys. Rev. B 97, 5 (2018) [7] V. Balédent, S. Chattopadhyay and P. Foury-Leylekian, Phys. Rev. L 114, 117601 (2015) [8] B. Khannanov, E. I. Golovenchits and V. A. Sanina, J.Phys.: Conf. Series 572, 12046 (2014)

**Images:**

![](https://assets.underline.io/uploads/markdown_image/1/image/bd2e503f7f47b740d4f6a1478408d52f.jpg)
Inelastic polarized neutron scattering at Q=(3/2,0,0) performed at THALES@ILL for 15 K. NSFx, SFx SFy and SFz stands for Non-Spin-Flip, Spin-flip with polarisation parallel to a, c and b crystallographic axis respectively.

![](https://assets.underline.io/uploads/markdown_image/1/image/bd2ea5b97909165cb1bda30c587ddab0.jpg)
Absorbance at 5 K with a reference at 40 K in Gd for e//a h//b & e//b & h//a (electric and magnetic field of the light) with respect to crystallographic axis.

 Study of a new type of excitation in the multiferroic magnetoelectric compound GdMn2O5, the electromagnon

Integrated magnetic devices show great potential for applications of radiofrequency and microwave integrated circuits, such as chip-scale magnetic inductors, bandpass filters, phase shifters, etc. [1, 2]. Recently, many researchers have been vastly concentrating on ferromagnetic/insulator sandwich structure for high-frequency applications. This thin film stack reduces the eddy-current as well as the out-of-plane anisotropy [3, 4], reducing the energy dissipation in form of heat. The FeGaB has soft magnetic nature, a large magnetostriction coefficient, and saturated at lower magnetic fields. Due to the low loss property of Al2O3 in microwave frequency, we have attempted to obtain low damping (α) and narrow linewidth through the FeGaB/Al2O3 thin film stack and characterize the ferromagnetic resonance (FMR) by utilizing Cryo-FMR. In this article, Ta (5)/ [FeGaB (15)/ Al2O3 (3)] n/Ta (3) / Al2O3 (2) thin-film stack are deposited by DC and RF sputtering on Si substrate, where n indicates the number of multilayer (n=1, 5, and 10) and parenthesis indicates thickness in nanometers. Fig.1 (a) shows normalized magnetic loop of the thin film stack, which reveals that the magnetic saturation moment has increased for increasing layer thickness. The coercive filed (Hc) has a low value below 3 Oe as indicated in Fig.1 (b). The dynamical magnetic properties have been studied using NanOsc Instruments Cryo-FMR in VersaLab system with temperature variation 300-100 K and exciting frequency from 2 to 20 GHz. The ferromagnetic absorption spectrum is shown in Fig.2 (a). The absorption spectrum contains the perpendicular standing spin wave (PSSW) modes for n=5 and n=10 in addition to FMR mode. The linewidth of absorption spectrum is fitted with dH= dH0 + αf/γ, (d=delta) to obtain damping factor (α), and inhomogeneous line width (dH0)[5]. The thin film stack with n=1 has a narrow linewidth and a low damping factor compared to other two thin film stacks, as shown in Fig.2 (b). When the temperature decreases, the damping factor has decreased for n=1, but increased for n=5 and n=10 due to defects arise for lowering the temperature. The Fig.2 (c) shows temperature dependent inhomogeneous linewidth for thin film stacks. As the temperature lowering to 100 K, the dH0 has increased for n=1 and n=10 thin film stacks, whereas decreased for n=5. The dH0 is lower for n=1 and larger value for n=5. The frequency dependence of FMR fields (Hres) plot provides effective magnetization (Meff) and magnetic anisotropy, when the plot has fitted to the Kittle equation. The Meff has a lower value for n=10 than other two stacks. The Meff has increased in all thin film stack as the temperature decreases. Since the absorption of the spectrum shows the PSSW modes in additional to FMR mode, the exchange magnetic fields (Hex) have been calculated by the subtraction method [6]. The plot of Hex and number of PSSW mode has provided the spin-wave exchange stiffness (D). The D values are 13.6615 and 10.7371 x10-18 Tm-2 for n=5, and n=10, respectively at 300 K. As the temperature decreases to100 K, the D has decreased for n= 5 thin film stack whereas increased for n=10. The temperature-dependent stiffness is following the T2 model in a linear form, and it suggests the existence of the electron–magnon interaction in this thin film stack. In order to calculate the constant of spin-wave exchange stiffness (A), saturation magnetization of thin film stack has multiplied with the exchange stiffness. At 300K, the A values are 7.1108, and 5.3213 pJ/m, respectively. In summary, we have deposited multilayer film stack of Ta (5)/ [FeGaB (15)/ Al2O3 (3)] n/Ta (3) / Al2O3 (2) onto Si substrates, and studied magnetodynamic (ferromagnetic resonance) properties with multilayer thickness and temperature dependence. The magnetic properties of thin film stack reveal that saturation magnetic moment has enhanced for increasing thickness of thin film stack, and Hc has obtained below 3 Oe. Therefore, the multilayer thin film stacks have the most significant for microwave applications. The dynamic magnetic properties suggest that the lowest damping and narrow linewidth are obtained for n=1 thin film stack compared to other thin film stacks. The FMR absorption spectrum of thin film stack with n=5, and n=10 show the PSSW modes in addition to the FMR mode. Consequently, the spin-wave exchange stiffness has been extracted. The temperature-dependent stiffness has followed the T2 model, which provides the strong presence of electron–magnon interaction in materials. The constant of spin-wave exchange stiffness has calculated with help of saturation magnetic moment and the exchange stiffness. Hence, the multilayer thin film stacks can be employed for microwave device applications.

**References:**

[1] Y. He, et.al,” IEEE Trans. on Magn., vol. 54, no. 9, pp. 1-4, 2018. [2] T. Wu, et.al, IEEE Magn. Lett., vol. 2, pp. 6000104-6000104, 2011. [3] X. Xing, et.al, IEEE Trans. on Magn., vol. 47, no. 10, pp. 3104-3107, 2011. [4] T. Wu, et.al, Appl. Phys. Lett., vol. 98, no. 1, pp. 012504, 2011. [5] K. Yadagiri, & T. Wu,J. Magn. Magn Mater., vol. 515, pp. 167277, 2020/12/01/, 2020. [6] K. Yadagiri, et.al, IEEE Trans. on Magn., vol. 57, no. 1, pp. 1-7, 2021.

**Images:**

![](https://assets.underline.io/uploads/markdown_image/1/image/958b8bb25594b09daba2dc26ecee751d.jpg)
Fig.1 (a) Normalized M-H loop for Ta (5)/ [FeGaB (15)/ Al2O3 (3)] n/Ta (3) / Al2O3 (2); n=1, n=5 and n=10, inset of Fig.1 (a) shows the Hc of thin film stacks, (b) coercive field HC (Oe) and Mr/Ms as a function of number of multilayer (n).

![](https://assets.underline.io/uploads/markdown_image/1/image/f01c6f6422c349af3f382d094dd43e8f.jpg)
Fig.2 (a) The absorption spectra of Ta (5)/ [FeGaB (15)/ Al2O3 (3)] n/Ta (3) / Al2O3 (2); n=10 at 300 K with exciting frequencies 2-20 GHz; inset of figure shows absorption spectra at 15 GHz with FMR mode and PSSW modes, (b) the damping factor (α) versus termeapture for thin film stack, and (c) the inhomogeneous linewidth (deltaH0) as a function of temperature for Ta (5)/ [FeGaB (15)/ Al2O3 (3)] n/Ta (3) / Al2O3 (2).

 Magnetodynamic properties of FeGaB/Al2O3 multilayer thin film stack for microwave applications

Magnetic materials are at the hearts of billions of high frequency communication devices and on-chip integration of magnetic microwave devices is of paramount importance for efficient utilization of the huge electromagnetic (EM) spectrum beyond 5G. Magnonics – a sub-field of spintronics – employs magnetic waves and oscillations in magnetic thin films and patterned nanostructures for processing GHz to sub-THz EM signal [1-2]. Magnetic ferrites are conventionally used in such magnetic microwave devices. However, due to the need for external bias magnetic field and incompatibility with complementary metal oxide semiconductor (CMOS) process, ferrites are not suitable for miniaturized microwave components. On the other hand, metallic magnetic thin films and nanostructures have great potential for efficient engineering of microwaves. Although, the need for bias field remains a bottleneck. We have been working on solutions where the requirement of the bias-field is absent. We have successfully demonstrated several strategies for achieving such device architecture by utilizing self-biased nanomagnets with reprogrammable remanent states which are associated with distinct magnetization dynamics [2-11]. To illustrate reconfigurable microwave properties, we have shown in Fig 1 as an example, a pair of dipolarly coupled rhomboid nanomagnets (RNMs). Two different remanent states (↑↑ and ↑↓) are obtained by utilizing a simple field initialization (HI) scheme: applying a saturating field along short and long axes of the RNMs with subsequent removal of the field (Fig. 1(a)). Due to different distributions of the stray fields for the two remanent states, the effective fields (Fig. 1(b)) are different which manifest in the shift of the microwave spectra for the two remanent states. Measurement of microwave properties from a single isolated nanomagnetic pair was obtained by using micro-focused Brillouin light scattering (micro-BLS) experiment (Fig. 1(c)) and the micro-BLS spectra are shown in Fig. 1(d). This strategy for reconfigurable magnonic responses is expanded to 2D arrays of the RNMs for different magnetic configurations: antiferromagnetic-type and ferrimagnetic-type artificial magnonic crystals. The maximum amount of the shift in the microwave responses is achieved by using multilayer isolated as well as dipolar-coupled networks of RNMs. Beyond, RNM shape, we have also demonstrated reconfigurable microwave operation based on other strategic shapes with multiple remanent states, such as arrow-, C-, L- and S-type nanomagnets and zigzag nanowires. Furthermore, we emphasize that such self-biased nanomagnets are also suitable for magnonic device applications where the propagation and gating of magnons – quanta for spin waves – are of primary interest. Conventionally, a large bias magnetic field is used in order to satisfy an orthogonal orientation between spin wave propagation and magnetization direction for efficient transmission characteristics (known as surface geometry or Damon-Eshbach mode). The need of such bias field hinders device integration on-chip for magnonic waveguides. Therefore, we have proposed and demonstrated a self-biased magnonic waveguide based on dipolar-coupled but physically separated chain of RNMs (Fig. 2). A prototype of such magnonic devices is shown in Fig. 2(a) where a microwave antenna (input) excites the spin waves and the spin waves are detected using micro-BLS technique with a focused laser spot (output). In order to gate the spin wave propagation, we have flipped a RNM in the waveguide and it is marked as ‘gate’ in the scanning electron microscopy (SEM) image (Fig. 2(b)). Such waveguide design enables one to have two different remanent states: ferromagnetically ordered (FO) and FO with a defect (FO*) when field initialized (same scheme as mentioned above) along the width and the length of the waveguide, respectively. FO and FO* remanent magnetic states are shown in Fig. 2(b) using magnetic force microscopy (MFM) images. Gating of the spin wave propagation has been shown by plotting 2D spatial profiles of the spin wave intensities using micro-BLS technique which clearly indicate gated magnon flow for the FO* state in comparison to the FO state. Thus, one can also achieve binary gating of output signal amplitude using such self-biased nanomagnetic waveguide which is remarkable as it also eliminates the need for bias field. In this presentation, we would like to emphasize the need for bias-free reconfigurable magnonics along with recent developments and future roadmap.

**References:**

[1] A. Haldar, C. Tian and A. O. Adeyeye, Sci. Adv. 3, e1700638 (2017). [2] A. Haldar, D. Kumar and A. O. Adeyeye, Nat. Nanotech. 11, 437 (2016). [3] A. Haldar and A. O. Adeyeye, ACS Nano. 10, 1690 (2016). [4] A. Haldar and A. O. Adeyeye, Appl. Phys. Lett. 108, 192404 (2016). [5] A. Haldar and A. O. Adeyeye, Appl. Phys. Lett. 108, 022405 (2016). [6] A. Haldar and A. O. Adeyeye, J. Appl. Phys. 123, 243901 (2018). [7] K. Begari and A. Haldar, J. Phys. D: Appl. Phys. 51, 275004 (2018). [8] K. Begari and A. Haldar, J. Phys. D: Appl. Phys. 52, 335003 (2019). [9] K. Begari and A. Haldar, J. Phys. D: Appl. Phys. 53, 455005 (2020). [10] A. Haldar and A.O. Adeyeye, Appl. Phys. Lett. 116, 162403 (2020). [11] A. Haldar and A. O. Adeyeye, J Appl. Phys. 128, 240902 (2020).

**Images:**

![](https://assets.underline.io/uploads/markdown_image/1/image/d84bb5df54b8e122952d84e6e7237690.jpg)
Figure 1: (a) Reprogrammable remanent states (X- and Y-state) of a RNM network. (b) Stray field distributions showing difference in the effective field for the two remanent states. (c) Measurement scheme from a single RNM network using micro-BLS technique. (d) Reconfigurable microwave responses for the two different remanent states.

![](https://assets.underline.io/uploads/markdown_image/1/image/408cd838fb91485f6c26689a83c74bc1.jpg)
Figure 2: (a) Propagating spin wave measurement scheme from a RNM chain using micro-BLS technique. (b) Logic device based on dipolar coupled RNM chain is shown in the SEM image. MFM images show the two different magnetic states (FO and FO*). Logic operation is demonstrated by using 2D spatial maps of the spin wave propagation for the FO and FO* states.

 Reconfigurable magnonics using self-biased reprogrammable nanomagnetic structures INVITED

Magnonics where magnons, the associated quasiparticles of spin waves are used to perform advanced logical or analogue functions represents a promising paradigm for new means of information processing for beyond complementary metal oxide semiconductor (CMOS) based architectures. Furthermore, in view of the emerging nano-magnonics which can eventually lead to a large-scale integration of magnonic devices new approaches for a compact, reconfigurable operation of magnonic device are required. Although such a reconfigurability is key to any efficient application such as magnonic filters, it has hitherto only theoretically [1] but not experimentally been shown on a sub-micron scale. On the other hand, multiferroics inherently offer compact means to implement voltage-control and reconfigurability in such microscopic magnonic systems [2]. In this work, we combine multiferroics with ferromagnets by employing a fully epitaxial heterostructure to achieve such voltage-controlled and reconfigurable magnonic system. Frequency filtering is most often based on magnonic crystals (MC) where a periodic modulation of the spin wave’s energy landscape results in frequency bandgaps in the magnon spectrum [3,4]. Here, we employ a thin-film multiferroic bismuth ferrite BiFeO3 (BFO)- lanthanum strontium manganite La2/3Sr1/3MnO3 (LSMO) heterostructure with a low Gilbert damping ~6x10-3 for the LSMO layer. Instead of modulating the sample shape, we create a MC by imprinting a periodic remnant electrical polarization in the multiferroic BFO layer which modulates the effective magnetic field in the ferromagnetic layer. We characterize the spin wave propagation spectra in this artificial, voltage-induced magnonic crystal and demonstrate the occurrence of a robust magnonic bandgap with &gt;20 dB rejection, i.e. filtering. The imprinted ferroelectric domains can be reconfigured repeatedly and are long-lived. Our results open a new path for the transduction and configurable filtering of spin waves within a single magnonic device. In general, bridging between between the scientific fields of functional oxides and magnonics, it opens new perspectives for the development of beyond CMOS based technologies [5].

**References:**

[1] Q. Wang , A. V. Chumak, L. Jin et al., Phys. Rev. B. 95, 134433 (2017) [2] N. Spaldin and R. Ramesh, Nat. Mat. 18, 203-212 (2019) [3] P. Frey , A.A. Nikitin, D. Bozhko et al., Commun Phys 3, 17 (2020) [4] A. V. Chumak, A. A. Serga and B. Hillebrands, J. Phys. D: Appl. Phys. 50, 244001 (2017) [5] H. Merbouche, I. Boventer, V. Haspot et al., submitted

**Images:**

![](https://assets.underline.io/uploads/markdown_image/1/image/243806faa5791ce2503db35d16d621e6.jpg)
Schematics of the experimental setup shown for the periodic ferroelectric domain state for the BFO layer.

![](https://assets.underline.io/uploads/markdown_image/1/image/4e0581eb0df6af24b4bcc64b618bde47.jpg)
Evidence for the bandgap for the BFO/LSMO heterostructure with periodic domains in the BFO in the amplitude representation in logarithmic scale. The data is normalized by the mean amplitude. At the bandgap frequency, the amplitude is dropping by >20 dB. A reconfiguration of the periodic domain state to the uniform down state erases the gap and reconstitutes the initial uniform state's response. Pn,m label the different domain states where n is the cycle number (1 cycle : down (virgin)-> up-> periodic-> ) and m the direction of the ferroelectric polarization.

 Voltage-Controlled Magnonic Crystal at the Sub-Micron Scale

In the last decades, the investigation of microwave spin waves (SW) in thin magnetic films has been an attractive field of research due to their very short wavelength reaching down to nanometers at GHz frequencies. SW also have a rich dispersion relation, that depends on their propagation direction with respect to a magnetic field. Hence, SW -based devices are promising candidates for microwave information processing, and eventually for overcoming the limitations encountered with CMOS-devices such as power consumption. In magnonic devices, a key challenge is to achieve long propagation distances of the spin waves. This requires to work with magnetic materials with small magnetic damping such as Yttrium Iron Garnet (YIG). (YIG), has by far the lowest magnetic damping which allows the spin waves propagation to be spread over millimeter-scale distances making it a reference material for spin-wave dynamics studies. Kajiwara and al. [1] confirmed that an electric current injection in a YIG film can be converted into spin-waves. In fact, when placing a non-magnetic metal with a high spin orbit coupling such as platinum on top of a magnetic insulator such as YIG, flow of charge current in the metal leads to its conversion into a pure spin current (via the Spin Hall Effect), which is then injected in the magnetic insulator where it adds up as an additional torque on the YIG’s magnetization. The above-mentioned mechanism of STT generation originating from pure spin-current obtained through the spin-orbit interaction of a heavy metal is called spin orbit torque (SOT). Recent studies [2], [3] with Platinum stripes placed on top of a 20 nm thick YIG waveguide (Pt/YIG) showed that SOT could generate a full compensation of the damping, leading to auto-oscillations of the magnetization above a critical injected current density. When excited with a microwave field, spin waves propagation length was increased by a factor of 10 in the bilayer. However, the possibility to achieve an amplification of propagating spin-waves was not observed yet, due to the onset of nonlinear dissipative processes above the critical current. More recently, materials having perpendicular magnetic anisotropy have demonstrated a full damping compensation along with SW emission, such behavior is accomplished by substituting Bismuth in YIG [4]. Here, we report on a study of SW propagation driven by SOT in a Pt/Bi-YIG waveguide where SW are detected using micro-focused Brillouin Light Scattering spectroscopy (BLS). The 20 nm thick and 21 µm long Bi-YIG waveguide is covered with 7 nm of Platinum. The experiments show the lossless propagation of a pulse of spin-waves excited by an rf-field generated by an antenna placed on top of the waveguide, under zero or positive damping during the course of 200 ns. Here, we show the full scale micromagnetic simulations done to emulate this experiment. The result shows a verification of the lossless propagation of SWs as seen in the experiment. By including a realistic Gilbert damping parameter and the Slonczewski torque to the Bi-YIG structure to compensate the magnetic losses, the results reproduce not only the steady state dynamics but also transitional dynamics when the rf-excitation is turned on and off at the ns time-scale. In summary, the inclusion of STT in the micromagnetic simulations showed an accordance with the experimental results, which verifies the ability to obtain a lossless propagation of SWs in the Bi-YIG waveguide and their non-suppression even after a hundred of nanoseconds after cutting off the excitation field.

**References:**

[1] Y. Kajiwara et al., “Transmission of electrical signals by spin-wave interconversion in a magnetic insulator,” Nature, vol. 464, no. 7286, Art. no. 7286, Mar. 2010 [2] M. Collet et al., “Generation of coherent spin-wave modes in yttrium iron garnet microdiscs by spin–orbit torque,” Nat. Commun., vol. 7, no. 1, Art. no. 1, Jan. 2016 [3] M. Evelt et al., “High-efficiency control of spin-wave propagation in ultra-thin yttrium iron garnet by the spin-orbit torque,” Appl. Phys. Lett., vol. 108, no. 17, p. 172406, Apr. 2016 [4] L. Soumah et al., “Ultra-low damping insulating magnetic thin films get perpendicular,” Nat. Commun., vol. 9, no. 1, Art. no. 1, Aug. 2018

**Images:**

![](https://assets.underline.io/uploads/markdown_image/1/image/2a26b6d8f9040dc1bef07c68aad030d7.jpg)
Fig.1: Two dimensional colormap showing the spin waves intensity in the simulated Bi-YIG waveguide after 80 ns from turning the rf-excitation on.

![](https://assets.underline.io/uploads/markdown_image/1/image/ea7d4882c8bc05ce1b677bad5a334f51.jpg)
Fig.2: Two dimensional colormap showing a lossless propagation of the SW intensity under SOT amplification, 100 ns after emission. The extracted group velocity is 55 m/s. The intensity of z-component of the magnetization is identical to that of Fig.1 recorded at a timestamp of 80 ns.

 Micromagnetic Simulations of Spin Waves propagation by SOT in a Bi-YIG waveguide

Conversion of heat into spin current has recently drawn a lot of attention due to its potential usage on heat recovery via the spin Seebeck effect (SSE) [1] as well as fundamental interests including generation of spin super currents [2] and Bose-Einstein condensation of magnons [3]. The underlying mechansim for SSE is understood as a combination of interfacial and bulk (intrinsic) contriubtion. The interfacial SSE is driven by the tempearture drop at the interface between metal and magnetic layer while the bulk SSE is driven by diffusive thermal magnons. In the nonlocal spin transport measurements, the long range exponential decay of thermally generated spin signal was explained solely by the intrinsic SSE [4] while there are several studies showing that there is another decay mechanism in subhundred nanometer scale [5-7]. Here, we study the spatial decay of thermally excited magnon currents inside a thin magnetic insulator by focusing on the short-range behavior in the nonlocal geometry. We compare SSE signal (Δ1,2) on the same device before and after adding a nonmagnetic Al capping layer. The Al capping alters the thermal profile significantly near the heat source, which results in that the nonequilibrium thermal magnon profile deviates from an exponential decay and shows two sign reversals within 1 μm from the heat source. From simulated temperature profiles, we find the vertical temperature gradient also reverses twice within the same length scale in the Al deposited case. The correlation suggests that thermal magnons behave locally at this length scale. Using a phenomenological spin transport model, we calculate the spatial profile of nonequilibrium spin density. The observed short-range behavior can be accounted if the relavant decay length is on the order of subhundred nanometers. Our result shows the existance of the shorter decay length in the nonlocal magnon transport experiments. This suggests that the long-range SSE signal in previous measurements should be attributed to magnons with lower energy.

**References:**

[1] S. R. Boona, R. C. Myers, and J. P. Heremans, Energy & Environmental Science, 7, 885 (2014). [2] D. A. Bozhko, A. A. Serga, P. Clausen, et. al., Nature Physics, 12, 1057 (2016). [3] Y. Tserkovnyak,S. A. Bender, R. A. Duine, et. al., Physical Review B, 93, 100402. (2016). [4] L. J. Cornelissen, J. Liu, R. A. Duine, et. al., Nature Physics, 11, 1022 (2015). [5] A. Kehlberger, U. Ritzmann, D. Hinzke, et. al., Physical Review Letters, 115, 096602. (2015). [6] A. Prakash, B. Flebus, J. Brangham, et. al., Physical Review B, 97, 020408. (2018). [7] J. S. Jamison, Z. Yang, B. L. Giles, et. al., Physical Review B, 100, 134402. (2019).

**Images:**

![](https://assets.underline.io/uploads/markdown_image/1/image/71fcbac32cc3cda92f6d18b1228e05ae.jpg)
Fig.1: Comparison of the experimentally measured local and non-local SSE signals (Δ1,2) at 1 and 2 mA injection and the numerically calculated chemical potential profiles with different input values of decay length λ for YIG/Pt/Si3N4 (case A) on a symmetric log scale. The insets show the schematic of the sample structure (bottom left) and temperature profile (top right).

![](https://assets.underline.io/uploads/markdown_image/1/image/898460b298a266d784aeaeab1dbbca16.jpg)
Fig.2: Same plots as Fig.1 for YIG/Pt/Si3N4/Al (case B). The result with λ = 0.1 µm reproduces the measured double crossing as shown, revealing the short range thermal magnon diffusion. The second crossing does not appear anymore when increasing λ to 0.5 or 1 µm. The insets show the schematic of the sample structure (bottom left) and temperature profile (top right).

 Short-range Thermal Magnon Diffusion in Magnetic Garnet

Spin waves are collective excitations of a spin system, where the magnetic moments precess with a spatially periodic phase shift that determines the wavelength [1]. From a technological point of view, the use of spin waves as signal carriers in future spintronic logic and memory devices has the advantages of a potentially lower power consumption and improved miniaturization prospects compared to the present charge-based CMOS technology [2]. For achieving nanoscale integration, topological spin textures such as magnetic vortex cores can be driven by alternating magnetic fields to excite coherent spin waves with short, sub-micrometer wavelengths, as shown recently for both synthetic ferrimagnetic systems and single ferromagnetic layers [3, 4]. In this contribution, we will demonstrate the versatility of this excitation mechanism by showing that the application of a static magnetic bias field makes the vortex core a manipulable spin-wave source that can be continuously translated in space without altering the resulting spin-wave dispersion relation [5]. We employ time-resolved scanning transmission x-ray microscopy (TR-STXM) to directly image vortex-core spin textures as well as the emission and propagation of spin waves from them. By applying a static in-plane magnetic bias field, we can continuously displace the vortex core from its equilibrium position in the center of a structured element. For a synthetic ferrimagnet (Co/Ru/Ni81Fe19), it is possible to shift the core relatively close to the edge while maintaining a radial, isotropic spin-wave emission pattern (see Figure 1). The wavelengths of the excited spin waves are not noticeably changed by the applied magnetic bias field, but remain tunable by the driving frequency in the same way as without bias field. For a single Ni81Fe19 layer with a thickness of the order of 100 nm, we observe a similar behavior for relatively small magnetic bias fields but when the vortex core approaches the edge of the disk, it also expands (see Figure 2). This expansion can be seen as a transformation from a zero-dimensional point source to a one-dimensional curved object where its shape and extent is controlled by the magnitude of the applied magnetic bias field. Such an expansion leads to a directional emission of the spin waves that propagate away from the core in mainly two opposite directions rather than radially as observed for non-expanded cores. For certain combinations of magnetic bias fields and driving frequencies, focusing effects appear. Such effects provide promising means for a controlled and directional propagation of spin waves without the need for additional patterning or waveguides.

**References:**

[1] A.G. Gurevich and G.A. Melkov, Magnetization Oscillations and Waves. Vol. 1, CRC Press (1996) [2] A.V. Chumak et al., Nat. Phys. Vol.11, p.453-461 (2015) [3] S. Wintz et al. Nat. Nanotechnol. Vol.11, p.948-953 (2016) [4] G. Dieterle et al., Phys. Rev. Lett. Vol.122, p.117202 (2019) [5] S. Mayr et al., Nano Lett. accepted (2021)

**Images:**

![](https://assets.underline.io/uploads/markdown_image/1/image/1ab356e578b624eb1b59b9aeeb5b4f0a.jpg)
Fig. 1: TR-STXM snapshots of a synthetic ferrimagnet, taken with out-of-plane magnetic sensitivity at the Fe L3 edge at a magnetic bias field of 5 mT. (a) Absolute absorption image showing the black vortex core shifted from the center of the structure (green circle). (b) Normalized contrast displaying the radial spin-wave emission from the core for a 1 GHz alternating field excitation.

![](https://assets.underline.io/uploads/markdown_image/1/image/2d254de94235b1183f441f680c4fed5a.jpg)
Fig. 2: TR-STXM images with out-of-plane sensitivity displaying the response of the vortex in a Ni81Fe19 layer to a magnetic bias field of 22.5 mT. (a) Absolute absorption image showing the expansion of the vortex core (encircled in green). (b) Normalized image displaying only the temporal changes with respect to the average magnetization, showing the spin waves emitted from the expanded vortex core in two directions for an excitation frequency of 8.6 GHz.

 Spin-wave emission from vortex cores under static magnetic bias fields

Synthetic antiferromagnets (SyAFs) are an excellent platform to explore novel spintronic and magnonic concepts with coupled magnetic moments[1,2]. Unlike a homogeneously magnetised single-layer ferromagnet, coupled magnetic layers are able to offer rich magnetic states where the competition between interlayer exchange coupling, external-field-induced Zeeman interaction as well as other magnetic anisotropy terms plays a role. This can expand into their dynamic regimes as the coupled moment nature inherently provides two eigenmodes, acoustic and optical modes[3-7], where time-dependent components of two coupled moments are oscillating in-phase (acoustic) and out-of-phase (optical)in a spin-flop (canted) regime. When we excite spin-waves by oscillating magnetic fields, the spatial symmetry/profile of microwave excitation determines which spin-wave modes are excited. In the simplest case, a uniform distribution of microwave excitation fields across a magnet can excite the uniform spin-wave mode (wavevector k = 0) as well as higher-order spin-wave resonance modes (k &gt; 0) with odd index numbers since the spatial profile of the microwave excitation and spin-wave mode amplitude (with phase) matches to each other in terms of symmetry. This in turn suggests that it should be possible to control the selection rules of spin-wave excitations by designing the spatial profile of excitation fields. Accessing hidden spin-wave states as well as tailoring spin-wave excitation efficiency by this approach has not been much explored in the past. We report in this work the current-induced-torque excitation of acoustic and optical modes in Ta/NiFe/Ru/NiFe/Ta synthetic antiferromagnet stacks grown on SiO2/Si substrates(Fig 1(a)). The two Ta layers serve as spin torque sources with the opposite polarisations both in spin currents and Oersted fields acting on their adjacent NiFe layers. This can create the odd symmetry of spatial spin torque distribution across the growth direction, allowing us to observe different spinwave excitation efficiency from synthetic antiferromagnets excited by homogeneous torques. We analyse the torque symmetry by in-plane angular dependence of symmetric and anti-symmetric lineshape amplitudes for their resonance (as shown in Fig 2(a-b)) and confirm that the parallel (perpendicular) pumping nature for the acoustic (optical) modes in our devices, which is in stark difference from the modes excited by spatially homogeneous torques. The terms parallel (perpendicular) pumping here signify that the torque (hence spin-wave excitation efficiency) is maximised when the oscillating fields and dc magnetic field are collinear(perpendicular) to each other respectively. This is evident from Fig 2(c-d) that strongly suggests the torque symmetry for the acoustic(optical) mode excitation being of form sin Φ(cos Φ) . We also generalize macrospin model for this particular spin-torque excitation geometry, which excellently supports our experimental observation. Our results offer capability of controlling spin-wave excitations by local spin-torque sources and we can explore further spin-wave control schemes based on this concept. We envisage that the control of spin-wave excitations in STT nanodevices will be useful for future spintronic and magnonic nano-devices.

**References:**

[1] R. A. Duine, K.-J. Lee, S. S. P. Parkin, and M. D. Stiles, Nat. Phys. 14, 217 (2018). [2] A. V. Chumak, V. I. Vasyuchka, A. A. Serga, and B. Hillebrands, Nat. Phys. 11, 453461 (2015). [3] F. Keffer and C. Kittel, Phys. Rev. 85, 329 (1952). [4] J. J. Krebs, P. Lubitz, A. Chaiken, and G. A. Prinz, J. Appl. Phys. 67, 5920 (1990). [5] S. M. Rezende, A. Azevedo, and R. L. Rodrguez-Surez, J. Appl. Phys. 126, 151101 (2019). [6] A. Sud, Y. Koike, S. Iihama, C. Zollitsch, S. Mizukami, H. Kurebayashi, arXiv:2012.12999 . [7] A. Sud, C. W. Zollitsch, A. Kamimaki, T. Dion, S. Khan, S. Iihama, S. Mizukami, and H. Kurebayashi, Phys. Rev. B 102, 100403(R) (2020).

**Images:**

![](https://assets.underline.io/uploads/markdown_image/1/image/b406e148e73969fbd07f012bcc683fa7.jpg)
FIG. 1: (a) A schematic of the sample geometry used for STT-FMR measurements in our study. (b) V obtained at different frequencies for acoustic and optical modes in our device for Φ = 55ο (the values for some frequencies have been scaled to show them properly). (c) A 2D colorplot of V as a function of applied field and frequency, measured for Φ = 55ο . (d) Numerical results calculated for the experimental conditions as in (c).

![](https://assets.underline.io/uploads/markdown_image/1/image/aaeaf82e3f75f0759e29b390a7923bff.jpg)
FIG. 2: (a-b) Angular dependence of symmetric (Vsym) and antisymmetric (Vasy) components of voltage for (a) acoustic mode measured at 8 GHz and (b) optical mode measured at 8.5 GHz. The solid curves are obtained using theoretical model. (c-d) The symmetry of torques obtained by dividing the Voltage by (c) sin 2Φ for acoustic mode and (d) cos 2Φ for optical mode. The dominant sin Φ (cos Φ) dependence of torques confirms the parallel(perpendicular) pumping configuration.

 Parity-controlled spin-wave excitations in synthetic antiferromagnets

Spin wave (SW) based computing i.e. magnonics relies on propagating SW as information carriers. Being quasiparticles, magnons (the SW quanta) have a fine lifetime characterizing the exponential decay of the magnons population. Consequently, up to now, all magnonic devices operate within a short time-window and on short propagation distances. Finding a SW amplification paradigm is hence a prerequisite for the development of magnonics as a credible CMOS alternative. Two schemes are currently envisioned for SW amplification: Parametric pumping1 and spin-orbit-torque amplification (SOT)2. Here we will discuss SOT amplification in BiYIG/Pt bilayer3 (Fig.1). Using micro Brillouin light scattering spectroscopy (µ-BLS) on a 500 nm wide Pt/BiYIG waveguides, we take advantage of the strong BLS signal of BiYIG (nearly two orders of magnitude the one of YIG3) to extensively study the influence of SOT on coherent propagating magnon modes. By passing current in the Pt layer, the spin accumulation resulting from the spin-Hall effect induces a positive torque that compensates the Gilbert damping of BiYIG when the magnetic field orientation is perpendicular to the waveguide axis. At subcritical values, the current has a linear effect that lowers the effective damping, thus significantly increasing the decay length of the propagating magnons, similar to what was observed previously in a YIG/Pt waveguide4. However, by increasing the excitation frequency of these propagating magnons for a constant external field, we observe a gradual but strong decrease of the modulation efficiency of their lifetime in our BiYIG/Pt waveguide. Once the current in the Pt layer reaches a threshold value, the vanishing effective damping allows for the onset of the auto-oscillation regime5 where the thermal magnons get incoherently amplified and the magnetization auto-oscillates in the main uniform mode. However, coherent propagating magnons strongly interact with the auto-oscillations6, preventing the observation of their propagation at continuous supercritical current. To circumvent these non-linear interactions, we study the auto-oscillation transient regime for the first time by performing time-resolved µ-BLS measurements. When the SOT is switched on, the uniform auto-oscillation grows exponentially, in agreement with the Slavin & Tiberkevitch non-linear auto-oscillator model7 in the linear regime limit. The onset of non-linearities only occurs when the auto-oscillation reaches its stable orbital, after few hundreds of nanoseconds typically. Using a rf and DC pulse scheme, we propagate rf-excited magnons before the onset of non-linearities at supercritical currents and observe the amplification of propagating magnons (Fig. 2-a). Similarly to what was observed in the continuous regime at subcritical currents, the modulation efficiency of the propagating magnons lifetime strongly decreases with frequency (Fig. 2-b). This experiment sheds new light on the influence of SOT on propagating magnons and pave the way toward a magnonic spin-wave amplifier using DC current.

**References:**

1 T. Brächer, P. Pirro, and B. Hillebrands, Phys. Rep. 699, 1 (2017). 2 V.E. Demidov, S. Urazhdin, A. Anane, V. Cros, and S.O. Demokritov, J. Appl. Phys. 127, 170901 (2020). 3 L. Soumah, N. Beaulieu, J. Ben Youssef, and A. Anane, Nat. Commun. 9, 3355 (2018). 4 M. Evelt, V.E. Demidov, and A. Anane, Appl. Phys. Lett. 108, 172406 (2016). 5 M. Collet, V. Cros, A. Anane, G. de Loubens, and O. Klein, Nat. Commun. 7, 10377 (2016). 6 B. Divinskiy, S. Urazhdin, S.O. Demokritov, and V.E. Demidov, Nat. Commun. 10, 5211 (2019). 7 A. Slavin and V. Tiberkevich, IEEE Trans. Magn. 45, 1875 (2009).

**Images:**

![](https://assets.underline.io/uploads/markdown_image/1/image/7c591b0f58bbf93bc6d6866a53fdf516.jpg)
Fig. 1 : Schematic representation of the 0.5µm-wide BiYIG/Pt waveguide. The DC current is applied in the Pt along the waveguide axis and a 300nm-wide strip rf-antenna is used to excite propagating magnons. The µ-BLS laser spot is focused from the back of the sample.

![](https://assets.underline.io/uploads/markdown_image/1/image/b8aea6e40370d823a1c720f70ddf0448.jpg)
Fig. 2 : Spin-wave intensity as a function of the distance to the antenna. The experimental curves (solid lines) are linearly fitted (dotted lines) in log scale for 4 frequencies (only 2 shown here) 4.925 GHz (a), 4.95 GHz, 4.975 GHz and 5.0 GHz (b) and 4 DC voltages corresponding to one subcritical current (1.7V) and three supercritical currents (2.2, 2.3 and 2.4V). The FMR is constant during the pulse with value 4.885 ± 0.01 GHz. 
The positive slope at supercritical currents at 4.925 GHz indicates an amplification of the propagating spin-waves (negative effective damping) due to SOT. The damping compensation efficiency gradually decreases as the frequency is increased, and the spin-waves are no longer amplified at 5.0 GHz for the applied voltages.

 Insight on the effect of SOT on propagating magnons: Spin-wave amplification and frequency dependent efficiency

Artificial neural networks (NN) running on traditional computing hardware tend to be relatively slow and power-hungry. One possible solution to this issue is to create hardware with neurons that spike in a manner similar to biological neurons. For example, Intel [1] and IBM [2] have recently fabricated neuromorphic chips based on spiking NN. Still, neurons in these chips are fabricated with multiple silicon transistors, which limits power efficiency of such an approach. The development of neuromorphic chips would benefit from a nanoscale, single element synthetic neuron that behaves in a manner that closely resembles a biological neuron. Such a synthetic neuron can be fabricated with nanometer dimensions from an antiferromagnetic (AFM) material mounted on platinum [3]. These “AFM neurons” are spin Hall oscillators driven by a sub-threshold spin current that can produce ultra-short (~ 5 ps) spikes in response to a weak external stimulus. In traditional spiking NN, synapses that connect artificial neurons are formed by complex active circuits. A much simpler approach would be to use passive memristors as artificial synapses [4]. However, memristors create symmetric coupling between neurons, which allows a spike produced by one neuron to travel to all neighboring neurons, both forwards and backwards. Unfortunately, this bi-directional signal propagation makes use of standard machine learning algorithms impossible. Therefore, it would be interesting to design NN in which unidirectional spike propagation is possible even with symmetrically coupled neurons. Here, we demonstrate that the unidirectional signal propagation is realizable in symmetrically-coupled NNs based on AFM neurons. This regime is possible due to an effective inertia that originates from exchange coupling between the two magnetic sublattices [3] and creates a relatively long refractory period - an interval of time after a neuron fires when it is unable to fire again. We show, that it is possible to adjust synapse weights such that each neuron is still in its refractory period when it receives a backward propagating signal from its neighbors. Then, a simple chain of neurons receiving an input at one end carry the signal to the opposite end without any reflections. However, for such a simple chain of AFM neurons, an input can occur at either end, and thus spikes would be able to travel in both directions. It would be beneficial to create an ‘isolator’ neuromorphic circuit, which would allow the propagation of spikes in only one direction. Figure 1 depicts the schematic of a possible ‘isolator’ circuit composed of AFM neurons. This figure uses circles to represent AFM neurons. The isolator is composed of neurons 4, 5, and 6. In the figure, double arrows represent symmetrical synapses with arrow width proportional to the coupling strength. This isolator allows spikes to be carried “forwards” from neuron 1 to 9, but stops signals carried “backwards” from neuron 9 to 1, as illustrated by the results of numerical simulations shown in Fig. 2. The isolator functions as follows. When the spike is initiated in neuron 1 (Fig. 2(a)), it travels along the chain to neuron 4. Then, the spike generated by neuron 4 travels to both neuron 5 and neuron 6. Neuron 5 also generates a spike, which is carried to neuron 6. The synaptic weights connecting to neuron 6 are small enough that a single input is not sufficient to induce a spike in the AFM neuron. However, when both spikes combine at neuron 6, they provide enough stimulus for neuron 6 to fire. The spike can then continue to propagate from neuron 6 to 9. When traveling backwards (Fig. 2(b)), a spike begins at neuron 9 and travels unidirectionally to neuron 6. Neuron 6 will generate a spike, but the synaptic weights from neuron 6 to 4/5 are insufficient to induce an additional spike. This, in effect, prevents the backward travel of spikes through the isolator. Using the proposed isolator NN, it is possible to design a number of other neuromorphic circuits supporting unidirectional information propagation, for example, neuromorphic combiners and unidirectional memory rings. The creation of such circuits with symmetric synapses shows hope for AFM neurons to be used to create more realistic spiking NN. The delay caused by the exchange inertia in the AFM neurons plays a key role in any neuromorphic circuit constructed with these neurons. These physically realistic neurons will increase the possibility for more functional spiking neural network hardware.

**References:**

[1] M. Davies et al., Loihi: A Neuromorphic Manycore Processor with On-Chip Learning, in IEEE Micro, vol. 38, no. 1, pp. 82-99, January/February 2018, doi: 10.1109/MM.2018.112130359. [2] M. V. DeBole et al., TrueNorth: Accelerating From Zero to 64 Million Neurons in 10 Years, in Computer, vol. 52, no. 5, pp. 20-29, May 2019, doi: 10.1109/MC.2019.2903009. [3] R. Khymyn, I. Lisenkov, J. Voorheis, O. Sulymenko, O. Prokopenko, V. Tiberkevich, J. Akerman, and A. Slavin, Ultra-fast artificial neuron: generation of picosecond-duration spikes in a current-driven antiferromagnetic auto-oscillator, Scientific Reports 8, 15727 (2018). [4] Sung Hyun Jo, Ting Chang, Idongesit Ebong, Bhavitavya B. Bhadviya, Pinaki Mazumder, and Wei Lu Nanoscale Memristor Device as Synapse in Neuromorphic System, Nano Lett. 2010, 10, 4, 1297–1301

**Images:**

![](https://assets.underline.io/uploads/markdown_image/1/image/eb70901938b9e7ce0e61457aefc96cde.jpg)
Fig 1: Chain of spiking neurons, configured as an isolator circuit. Gray circles represent neurons, and black arrows represent coupling between neurons. Thick black arrows represent symmetric synaptic weight tuned such that a spike cannot flow backwards. Thin black double arrows represent symmetric synaptic weight weak coupling.

![](https://assets.underline.io/uploads/markdown_image/1/image/7460602315b4c8ef8afb70fa1a01a721.jpg)
 
Fig 2: Results of simulation of AFM neurons configured as an isolator circuit. (a) Initial spike is in neuron 1. The spike travels through the chain, to neuron 9. (b) The initial spike is at neuron 9. The spike travels to neuron 6, then is unable to continue.

Downloads

Next from INTERMAG 2021

Study of a new type of excitation in the multiferroic magnetoelectric compound GdMn2O5, the electromagnon

Similar lecture

Synchronization of Spin Torque Nano

.css-70qvj9{display:-webkit-box;display:-webkit-flex;display:-ms-flexbox;display:flex;-webkit-align-items:center;-webkit-box-align:center;-ms-flex-align:center;align-items:center;}Downloads

Next from INTERMAG 2021

Study of a new type of excitation in the multiferroic magnetoelectric compound GdMn2O5, the electromagnon

Similar lecture

Synchronization of Spin Torque Nano

Downloads