United States

In recent years, the research of long-distance transport of magnon in magnetic insulator showed a negative voltage in thermal excited magnons signal [1,2]. Subsequent research investigated the negative voltage and find its dependence on the thickness of YIG, which raising a challenge for traditional magnon diffusion model. Similar phenomenon of viscous electrons’ backflow in graphene inspired the magnon study. In graphene at low temperature where the scale of electron-electron scattering is much smaller than the scale of electron mean free path, the momentum-conserving process dominated will give make electrons behave as a viscous liquid and exhibit hydrodynamic phenomena as well as classical liquids [3], which could present a negative voltage at two edge of system. In magnons system, same momentum-conserving process would happen, even though magnon as a quasiparticle, people don’t emphasize the conservation laws of particular numbers. The conservation of magnon momentum will lead whole system relaxed to an collective equilibrium N(λ)=1/exp[(ε&lt;sub&gt;k&lt;/sub&gt;-k●λ)/k&lt;sub&gt;B&lt;/sub&gt;T] , which is called displaced Planck distribution [4]. The collective velocity λ of the ordered motion as a function of coordinates must satisfy a Navier-Stokes type hydrodynamic equation. The backflow of viscous magnon naturally produce the negative signal in magnetic insulator system, as similar behaviors in graphene.&lt;br/&gt; &lt;br/&gt; In this work, based on Boltzmann method and considering momentum-conserving dominated process, we developed a collective hydrodynamic diffusion equation of magnons’ viscous fluid, without vanishing k-space scattering. We also found nonlinear properties of hydrodynamics magnon system. Two- dimensional calculations of magnon flux showed the viscosity of magnon fluid created two eddies near the heat source. And our result could not only fit experimental data precisely, but also reconstruct physical scattering process.

2022 Joint MMM-INTERMAG

Collective Hydrodynamic Magnon Transport Study in Magnetically Ordered Insulator Based on Boltzmann Approach

The 15th Joint MMM-INTERMAG Conference (2022 Joint) 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 plenary session, and an evening session, with about 1500 presentations overall. 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 2022 Joint Conference will offer both in-person and prerecorded on-demand content.

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

Please register for the Joint MMM-Intermag conference to join the event.

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

In recent years, the research of long-distance transport of magnon in magnetic insulator showed a negative voltage in thermal excited magnons signal [1,2]. Subsequent research investigated the negative voltage and find its dependence on the thickness of YIG, which raising a challenge for traditional magnon diffusion model. Similar phenomenon of viscous electrons’ backflow in graphene inspired the magnon study. In graphene at low temperature where the scale of electron-electron scattering is much smaller than the scale of electron mean free path, the momentum-conserving process dominated will give make electrons behave as a viscous liquid and exhibit hydrodynamic phenomena as well as classical liquids [3], which could present a negative voltage at two edge of system. In magnons system, same momentum-conserving process would happen, even though magnon as a quasiparticle, people don’t emphasize the conservation laws of particular numbers. The conservation of magnon momentum will lead whole system relaxed to an collective equilibrium N(λ)=1/exp[(εk-k●λ)/kBT] , which is called displaced Planck distribution [4]. The collective velocity λ of the ordered motion as a function of coordinates must satisfy a Navier-Stokes type hydrodynamic equation. The backflow of viscous magnon naturally produce the negative signal in magnetic insulator system, as similar behaviors in graphene. In this work, based on Boltzmann method and considering momentum-conserving dominated process, we developed a collective hydrodynamic diffusion equation of magnons’ viscous fluid, without vanishing k-space scattering. We also found nonlinear properties of hydrodynamics magnon system. Two- dimensional calculations of magnon flux showed the viscosity of magnon fluid created two eddies near the heat source. And our result could not only fit experimental data precisely, but also reconstruct physical scattering process.

poster

The antiferromagnetic (AFM) medium is an ideal candidate for emitting terahertz (THz) electromagnetic (EM) waves, which have extraordinary potential application in information technology and biological medicine [1–4]. In this work, we report the theoretical and numerical investigation about the emission of terahertz (THz) spin waves from a moving AFM DW under the magnetic anisotropy energy gradient (dEa/dx). The DW is accelerated by dEa/dx that can be generated using the technique displayed in Fig. 1: A DC voltage applied across the wedge-shaped insulating layer generates dEa/dx, which triggers the AFM DW to move towards the end with a lower magnetic anisotropy energy [5]. Based on an atomistic spin model, we have found that the DW velocity exceeds the relativistic limit and keeps increasing, which is accompanied by the broadening of the DW width. This is different from a traditional relativistic AFM DW motion with an upper limit of the DW velocity and Lorentz contraction [6]. The DW acceleration effectuates the DW precession. When the frequency of the DW precession is greater than the frequency gap of the spin wave propagating in the AFM medium, a continuous THz spin wave is resonantly emitted. The frequency of the spin wave keeps increasing with the decrease of the anisotropy energy constant, which seems to be contrary to the AFM spin wave obeying the relativistic-like dynamics. Our results demonstrate that obeying the relativistic dynamics is not necessary for the spin-wave emission from a moving AFM DW.

![](https://assets.underline.io/uploads/markdown_image/1/image/b3d5ff6329dec082c7c8fd30cb4693b4.jpg)

Fig.1. Schematic of calculation model: emission of THz spin wave from a moving DW triggered by a voltage-induced gradient of magnetic anisotropy energy (the emission of the spin wave is accompanied by an increase of the DW width due to the reduction of the magnetic anisotropy energy).

Resonant spin-wave emission from the moving AFM DW without Lorentz invariance

Magnetic domain walls (DWs), as a novel information carrier, provide a reliable technique for modern information processing. In previous study, to control the movement of magnetic DWs, different external excitations, such as magnetic field, spin transfer torque, spin-orbit torque, electric field, and spin waves, etc., have been used as the driving force for the DW motion. In this work, we combine electromechanical and micromagnetic simulations to explore a promising way to manipulate the perpendicular magnetic anisotropy (PMA) DW motion in a given direction using voltage-controlled gradient strain. In particular, we show the voltage induced gradient strain driven the DW motion more complicated than simple uniform linear motion, especially for DW motion in a finite length nano-stripe. As shown schematically in Fig. 1(a), the proposed device consists of a ferromagnetic (CoFeB)/heavy metal (FM/HM) nanostrip deposited on a piezoelectric substrate. The initial spin configurations with different i-DMI coefficient are obtained by relax the DW along the y direction. Fig. 1(b) shows the snapshots of magnetization with various i-DMI coefficient when simulated 10 ns after the strain gradient applied. As shown in Fig. 1(b), the direction of the DW will parallel to the x direction due to the i-DMI and form a Neel type DW, while a Bloch type DW formed without i-DMI. From the micromagnetic simulation results shown in Fig. 1(b), along with the DW displacement, DWs are found to tilt away from the y-axis, which agrees with the experimental observations. From these magnetization images, we can also directly see an increment of DW tilting angel as the i-DMI constant increase. These results indicates that strain gradient induced PMA DW motion can be used to develop next generation artificial intelligent memory and logic devices.

![](https://assets.underline.io/uploads/markdown_image/1/image/66bfbae9388ee43210a60dde3c0ce6ad.jpg)

Fig. 1 (a) Schematic view of the investigated strain driven domain wall in a nanostripe with voltage-controlled strain gradient. (b) Snapshots of the magnetization pattern with different i-DMI coefficient.

Dynamics of Domain Wall Induced by the Voltage-controlled Strain Field Gradients

Spin-based neuromorphic computing (NC) is investigated as a power-efficient enabler for artificial intelligence [1-3]. Hence, magnetic tunnel junction (MTJ) based domain wall devices are investigated as synthetic synapses and neurons to mimic their biological counterparts [2-3]. While a biological neuron resets itself after firing, such a function has not been demonstrated experimentally in spintronic systems [3]. The self-reset feature of neurons reduces the power consumption and the complexity of NC architecture. In this work, we present a domain wall (DW) based artificial neuron that automatically resets without applying a reset pulse. At first, we simulated a synthetic antiferromagnetic (SAF) DW device where the free layer (CoFeB) is driven by spin-transfer-torque (STT) to achieve integrate function. Leakage function was achieved when the antiferromagnetic coupling (AFC) between the free and pinned layers pushes the DW towards the starting position automatically (Fig. 1). The left quarter region of the free layer is ferromagnetically coupled (FC) to the pinned layer so that as DW returns (during reset), it is stabilized at the initial position (Fig. 2(a) and Fig. 2(b)). Experimentally, the integrate and leakage functions are also achieved as illustrated in Fig. 2(d) and Fig. 2(e), suggesting that an SAF-based neuron is viable for NC.

![](https://assets.underline.io/uploads/markdown_image/1/image/4f911e30d06059db93a1e69c895660a8.jpg)

Fig.1 Proposed design of the artificial neuron. Domain wall is moved forward by STT and (after firing), reset backward by SAF coupling.


![](https://assets.underline.io/uploads/markdown_image/1/image/c297c876d46646d24593590c9f897b4d.jpg)

Fig. 2 (a and b) DW position (free layer) vs time graphs for different AFC strength in the simulations. (c, d and e) Kerr images of the fabricated SAF neuron device with magnetic field in the -z direction. (c) No current applied and device is in antiparallel configuration, (d) current is applied, pushing the DW and achieving the integrate function, and (e) current is removed and the device automatically resets, indicating the leakage function.

Synthetic antiferromagnetic domain wall device as artificial neuron for neuromorphic computing

Domain wall (DW) dynamics has been intensively studied as the base for novel spintronic devices [1-2]. These dynamics have attracted special attention in ultrathin ferromagnetic strips that show high perpendicular magnetic anisotropy (PMA) since the formation of ‘up’ or ‘down’ domains, separated by DWs, can be exploited to store bits. One of the challenges is the initial nucleation of the domains in sequence, typically performed by an Oersted magnetic field from an attached microwave ‘bit line’ orthogonal to the strip [3]. However, fabrication constraints due to the size of the bit line sets a lower bound for the bits width and, more importantly, the Oersted field that nucleates a DW can annihilate a shifted DW in the strip [3]. Solutions based on a double bit line have been proposed [3], but these still require injecting external currents and a more intricate and multilayer design. In this work, we propose a different approach based on a fluctuating DW as the source. The DW is pinned at one end at the junction of a Y-shaped strip (see red solid circles, Fig. 1) as the DW fluctuates across the strip width (see Fig. 1(b)-(d)). The pinned DW can be originally nucleated by using a bit line. Then, it is forced to oscillate by applying non-localized, pulsed currents (JHM) in the heavy metal (HM) in a controlled manner. Such oscillation injects subsequent DWs that travel along the strip. We characterize the obtained ‘domain sequence’ as a function of duration (d) and periodicity (T) of the current pulses. Particular combinations result in ‘coherent’ sequences of domains (see red dashed circles in Fig. 1(e)-(f)) and/or domain reversal, which may be useful to develop DW-based applications.

![](https://assets.underline.io/uploads/markdown_image/1/image/8e6fb080c602a3bfb788123271a304d5.jpg)

Fig.1. Nucleation of two domains (see red dashed circles in (f)) in the Y-shaped ferromagnetic nanostrip. Black and white are ‘down’ and ‘up’ domains, respectively (see green symbols), and grey represents the HM layer beneath the strip. A pulsed current (JHM) is applied in the x-direction over the entire HM. (a)-(d) Snapshots from a full fluctuation of the DW during a pulse of 0.6 ns duration, and (e) after a stabilisation period (JHM=0) of the nucleated domain (~90 nm width) at 1.4 ns. (f) Results after two pulses of period 1.4 ns.

Fast and Controlled Sequential Nucleation of Narrow Domains in a Ferromagnetic Nano-strip with Non-localised Electrical Currents

One of the tools to move a skyrmion is an injected current[1,2], which gives rise to the risk of producing excessive heat. On the other hand, spin-waves (SW) are a promising way to drive skyrmions without producing such high temperatures and hence affecting the skyrmionic device's performance. Most recent work on SW-driven skyrmion motion is based on micromagnetic simulations on single-layers with either small magnetic track width (W), high SW excitation rf field (Arf), or low Gilbert damping constant α[3,4]. However, to complement the realization of skyrmion motion experimentally, simulations should meet experimentally achievable conditions at the same time: W ≥100nm, Arf ≤100mT and α &gt;0.05. The dataset and characteristic analysis we present concern SW driven skyrmion simulations by OOMMF[5] with fixed parameters, D=1.5 mJ/m2, MS=1004 KA/m, A=3.6 pJ/m, Ku=1 MJ/m3. We varied the α=0.02-0.06, W=40-100 nm and Arf=100-800 mT. The size of the relaxed skyrmion is from 20nm to 23nm when W increases from 40nm to 100nm, so there is more than 5 nm between the skyrmion to each sample edge at initial stage. Fig.1 shows the skyrmion displacement (dsky) with time along the SW propagation direction for various W is shown. When W &gt;50 nm, the skyrmion will initially move opposite to the SW propagation direction before moving forward at longer times. Fig.2 shows that with larger α, the skyrmion will have smaller initial velocity when moving backward. If the leading edge of the skyrmion is regarded as a planar wall on the sample, the domain wall (DW) plane is predicted to rotate when linear momentum transfer between the SW and DW occurs. The angle between the DW wall plane and the out-of-sample direction is defined as DW wall plane angle (ΦDW)[6]. Also, the inset shows that ΦDW of skyrmion leading edge has similar behaviour to the dsky, but with an offset. The skyrmion backward motion gives a possible observation result in the future and guidance about finding a balance between W and magnetic properties to have skyrmions with aiming motion type.

![](https://assets.underline.io/uploads/markdown_image/1/image/20d3a9acd993b9c5372ac498eb39b2f1.jpg)

Fig.1 Time-dependent skyrmion displacement in nm for various W.

![](https://assets.underline.io/uploads/markdown_image/1/image/3a0f9f542c8a882ca49dfb2dc79008fd.jpg)

Fig.2 Skyrmion displacement with time for various α. Inset: comparison between ΦDW and skyrmion displacement when α=0.06.

Transient retrograde motion of spin wave-driven skyrmions in multilayer wave guides

Magnonic devices are based on the wave properties of spin wave (SW), so that one can transport more information in a single device if modes of SWs are freely controllable. SWs in a narrow waveguide have eigenmodes due to standing waves formed in the width direction and the thickness direction [1, 2]. However, excitation of one specific mode of propagating SW has not been realized, yet. Here, we report that high SW mode can be selectively excited by employing SW conversion. We numerically investigated SW propagation in a T-shaped structure using MuMax3 [3]. A 5-nanometer-thick Py film is assumed as the device and a magnetic field of 70 mT is applied along the y-direction (Fig. 1). By applying an rf field to the cells 1 μm away from the junction, Damon-Eshbach spin waves (DESWs) are excited and travel along the x-direction. The DESWs induce a magnetization precession in the junction, which excites SWs in the magnetostatic backward volume wave configuration (BVSWs). First, we varied the SW resonant frequency (f=12, 18 GHz) using a device having w1=500 nm and w2=100 nm. As shown in Fig. 2(a) and (b), the wavelength of the DESW (λDESW) depends on f which is well known [2]. In addition, the BVSWs clearly change their wave front. To extract the mode number, we applied a Fourier transformation (FT) along the y-direction. The FT results indicate that the 1st mode is dominantly excited at 12 GHz but the 2nd mode is dominant at 18GHz (Fig. 2(d) and (e)). We next changed the device to have w2=500 nm and repeated the simulation at 18 GHz (Fig. 2(c)). The result indicates the 12th mode was dominant at this condition (Fig.2 (f)). To further understand the high-mode-excitation, we varied the frequency and plotted the FT peak height of each mode as a function of the ratio of the width to λDESW (Fig.2 (g)). The result implies that modes can be controlled by optimizing a T-shaped waveguide.

![](https://assets.underline.io/uploads/markdown_image/1/image/29fb2eeaedfa09fa59c4d092aa944455.jpg)

Fig.1. Overview of the calculation area. SWs are excited at the orange area.

![](https://assets.underline.io/uploads/markdown_image/1/image/fc18c5114e3e73123e5e447b234cee46.jpg)

Fig.2. (a)-(c) The z-component of the magnetization and (d)-(f) their corresponding FT of BVSWs. (f, w2) is (a, d) (12 GHz, 100 nm), (b, e) (18 GHz, 100 nm), and (c, f) (18 GHz, 500 nm). (g) The peak height of FT of each mode obtained from the device with w2=500 nm.

Mode selective excitation of spin waves

Spin torque nano-oscillators (STNO) generate high frequency output voltage signals upon injection of a spin polarized DC current into a magnetic tunnel junction. For their applications in wireless communications or neuromorphic computing schemes, a central aspect is the control of their phase dynamics via injection locking to an external rf signal [1, 2]. Many interesting, non-linear effects may arise at higher excitation levels that manifest themselves in the Arnold tongue [3, 4]. For the in-plane magnetized STNO (in-plane magnetized polarizer and free layer) [2], we provide here analytical expressions for (i) the asymmetric Arnold tongue boundaries, arising from the dependence of the forcing on the STNO power; (ii) the Takens-Bogdanov bifurcation threshold [4] which leads to a reduced locking range and (iii) the hysteresis boundaries [3,4]. These effects are discussed within a common model for ‘pure’ rf signals arising from either a spin polarized current or a magnetic field. ‘Pure’ refers to a well-defined direction ui (i=x,y,z ; ux the STNO easy axis) of the spin polarization or a magnetic field [2]. In a further study we also investigate synchronization at harmonic order n=1 (rf signal frequency close to the STNO frequency) for mixed spin polarization with x- and y- components. The important and unexpected finding here is that the STNO phase does not follow an arcsine (Adler-type) behavior [5] within the locking range, but shows a discrete jump that can be as large as Pi/2. Our analysis shows that the phase jump results from a competition between locking to the ux- and uy-components of the rf signal, where locking to the ux-component at n=1 is enabled through the fact that the free running oscillation power is non-constant with frequency components at f. This demonstrates that, for mixed rf signals, injection locking can be more complex, and that the results can be of interest for applications such as the phase jump for phase shift keying concepts [6]. MIG was supported by the French Space Agency (CNES) and the CEA NUMERICS program, which has received funding from the European Union's Horizon 2020 program (No 800945). Partial support was provided by the French National Research Agency (ANR-20-CE24-0002).

Arnold tongue of injection locked spin torque nano-oscillators under high and mixed driving signals

We theoretically investigate the transmission of exchange-dominated spin waves in a one-dimensional magnonic crystal (MC) with a periodic exchange bias field. By recasting the Landau-Lifshitz-Gilbert equation into an effective Schrödinger equation and establishing spin-wave functions, it is found that MCs with upward (up) and downward (down) magnetization, respectively, correspond to the rectangular N-fold barriers and wells for magnons. The broadband transmission spectra in up and down states are systematically investigated. We show the phenomena of the magnon resonant transmission in both states and calculate the resonant transmission wave functions, which are related to the magnon density. Our results also show a transmission spectra shift effect (TSSE) between up and down states, which is found to be general in this system. Furthermore, the TSSE is useful to design a type of magnon valve, the magnonic-crystal-based magnon valve (MCMV), which has a large on/off ratio and bandwidth. By high-throughput screening, 125 000 groups of parameters of the MC are calculated, and 1948 parameter groups of high-performance MCMVs are screened out. Our work clarifies the physical details of the exchange-dominated spin-wave transmission in rectangular N-fold barriers and wells and also provides a promising route for designing novel magnonic devices.

![](https://assets.underline.io/uploads/markdown_image/1/image/06606510a3d9fad6f0a47ccb6c39f8f6.jpg)

The magnon transmission in a one-dimensional MC. The schematic diagrams of the MC in (a) the up state and (b) the down state. (c) The broadband transmission spectra in up and down states, corresponding to the black curve (T↑) and red curve (T↓) respectively. The parameters of the MC are the following: N = 14, a = 10 nm, and b = 25.2 nm. (d) The difference between T↑ and T↓ as a function of ω/2π. ω/2π is shown in logarithm scale.

Magnon valve effect and resonant transmission in a one-dimensional magnonic crystal

The study of high-frequency spin dynamics in three-dimensional (3D) magnetic nanostructures may lead to paradigm-shifting, next generation devices including high-density spintronics and neuromorphic systems. Despite remarkable progress in fabrication [1], the measurement and interpretation of spin dynamics in complex 3D magnetic structures remains exceptionally challenging. Here we present the measurement of coherent spin waves within a 3D artificial spin ice (ASI) system using Brillouin light scattering. The 3D-ASI was fabricated by using a combination of two-photon lithography and thermal evaporation. Two spin-wave modes were observed in the experiment whose frequencies showed nearly monotonic variation with the applied magnetic field strength. Numerical simulations qualitatively reproduced the observed spin-wave modes. The simulated mode profiles revealed the collective nature of the modes extending throughout the complex network of nanowires, while showing spatial quantization with varying mode quantization numbers [2]. The study shows a well-defined means to explore high-frequency spin dynamics in complex 3D spintronic and magnonic structures. &lt;!--![endif]----&gt;

![](https://assets.underline.io/uploads/markdown_image/1/image/03908b10f9b4a86f48a8879cc94464fd.jpg)

Figure 1: Schematic of BLS measurement geometry. The measurement was performed at θ = 0°. (a) A scanning electron micrograph image of the sample (scale bar: 2 µm) with applied field (H) direction, (b) BLS spectra at H = 1.8 kOe. (c) Schematic of 3D-ASI network, (d) simulated equilibrium state of magnetisation at H = 1.6 kOe and (d) mode profiles of M1 at 1.6 kOe are shown. Respective legends are shown on right side.

Study of Coherent Spin Waves in a Three-Dimensional Artificial Spin Ice Structure

Magnonic and phononic crystals are artificial nanostructures with periodic lateral variations in their magnetic and acoustic properties [1,2]. Recently, there has been growing interest in the fundamental understanding of spin waves propagation in magnonic crystals because of their huge potential in a wide range of applications such as microwave resonators, filters and spin wave logic devices [3]. On the other hand, phononic crystal are periodic arrays of nano-objects with alternating elastic properties and mass densities, they have promising applications such as acoustic microsensing or signal processing. Besides these possible applications, the wealthy physics of the acoustic waves is one of the main reasons for the attention to their propagation. Finally, a magphonic crystal is also an artificial crystal that simultaneously presents magnonic and phononic behaviors [4,5]. The phononic and magnonic dispersions of Ni80Fe20 antidot arrays deposited on a Si substrate have been investigated by using Brillouin light scattering (BLS) and ferromagnetic resonance (FMR). The periodicity of the arrays is 400 nm with a nanohole diameter of 180 nm and with a thickness of 40 nm. The FMR technique has been employed to first investigate the magnetostatic modes in the center of the first Brillouin zone. We have thus highlighted the presence of several magnetostatic modes and we conducted micromagnetic simulations (OOMMF) that allowed us to determine the spatial location of these different modes. In a second step, we have used BLS to probe the phononic and magnonic features of the magphonic crystal. We have experimentally demonstrated the simultaneous presence of magnonic and phononic bandgaps. The different dispersion curves (both magnonic and phononic) as well as the potential applications of these magphonic crystals will be presented and discussed during the conference.

Next from 2022 Joint MMM-INTERMAG

Resonant spin-wave emission from the moving AFM DW without Lorentz invariance

Similar lecture

12/10 Switched Reluctance Motor with a Hex Connection