We propose a hybrid magnonic-oscillator system based on the combination of a spin transfer auto-oscillator and a magnonic waveguide to open new perspectives for spin-wave based circuits. The system is composed of a spin transfer oscillator based on a vortex state which is dipolarly coupled to a nanoscale spin-wave waveguide with longitudinal magnetization. In its auto-oscillating regime, the oscillator emits coherent spin waves with tunable and controllable frequencies, directions, and amplitudes into the waveguide. We demonstrate the principle of this method using micromagnetic simulations and show that reconfiguration of the system is possible by changing the chirality and polarity of the magnetic vortex. Spin waves are emitted into the waveguide with high non-reciprocity and the preferred direction depends on the core polarity of the vortex. In contrast, different vortex chiralities lead to different amplitudes of the emitted waves. Our findings open up a novel way to design an agile spintronic device for the coherent and tunable generation of propagating spin waves.
I. INTRODUCTION
The investigation of spin waves goes a long way back to the twentieth century.1,2 Since then, research in these fundamental magnetic excitations has been gaining increasing attention. Spin waves are studied due to fundamental interest,3,4 as well as for their large potential for information processing applications, such as Boolean spin-wave logic devices or alternative computing schemes.5–8 Perhaps, some of the most important advantages of spin waves are their ability to transport information without generating ohmic losses,9–11 as well as the possibility of being used in nanoscale spintronic signal-processing devices.6,8,12,13 When it comes to the realization of spin-wave devices, some of the most important ingredients are ways to efficiently excite spin-waves. The most common method for this task is direct spin-wave excitation by microwave currents in antennas.14–17 Further methods involve the utilization of the spin-transfer torque (STT)18–20 or ultrashort light pulses.21–23 Another developed method of spin-wave excitation is proposed to use the driven rotational dynamics of the formed magnetic vortex cores in an antiferro-magnetically coupled trilayer.24 The vortex configuration is a region of in-plane curling magnetization at the edges and an out-of-plane magnetization at the center. It is characterized by two degrees of freedom, namely, polarity (P) and chirality (C). The polarity describes the direction of the out-of-plane magnetization at the vortex core, with positive/negative polarities referring to upward/downward directions of magnetization. The chirality describes the direction of the in-plane magnetization, with positive/negative chiralities referring to clockwise/counterclockwise directions of magnetization.25 When STT is used to drive the auto-oscillation of the magnetic vortex, we refer to these devices as spin-transfer torque vortex oscillators (STVOs).26,27 STVOs are magnetic oscillators in the GHz range which are driven by a DC current and thus constitute very promising and energy efficient sources for GHz excitations.28–30 Due to their intrinsic non-linearity and their ability to synchronize their oscillation frequencies, STVOs are already used in novel logic concepts such as brain inspired and neuromorphic computing31–34 or vowel recognition.35 On the other hand, magnons, the eigen-excitations of the spin system, are proposed as data carriers for future beyond-CMOS logic. Several logic building blocks such as magnon transistors,36,37 majority gates,38 or nonlinear directional magnon couplers39 have been realized. Based on this, we propose to merge these two important sub-fields of spintronics, namely, the field of spin-transfer torque vortex oscillators (STVOs)40 and the field of magnonics.41–43 Here, we demonstrate a novel hybrid device concept which unifies the particular strength of both fields, which we refer to as the hybrid magnonic-oscillator system (HMOS). The proposed device consists of a vortex-based spin-transfer torque oscillator, which is dipolarly coupled to an adjacent magnonic waveguide. A coupled system where spin waves propagating in a waveguide are coupled to spin-wave modes in a vortex state has been previously studied in Ref. 44. In this case however, no auto-oscillation and vortex core dynamics have been involved. Here, we numerically study this system using the open source micromagnetic solver MuMax 3.0 package,45 as well as the software platform Aithericon.46 We demonstrate that, when carefully designed, this hybrid system can be used as a tunable source for the emission for coherent spin waves in the magnonic waveguide. In the following, we investigate the spin-wave dispersion relation in the bare magnonic waveguide. We then numerically investigate the coupled HMOS system and the emission of spin-waves by the vortex oscillator. In particular, we study the effects of changing the polarity (P) and chirality (C) of the vortex on the emitted signal intensity and directionality. Finally, we study the evolution of the spin-wave frequency, intensity, and non-reciprocity as a function of the current density.
II. SIMULATION SETUP
The proposed hybrid magnonic-oscillator system is schematically depicted in Fig. 1. It consists of a magnonic waveguide made of yttrium iron garnet (YIG) with a uniform magnetization directed along its long axis (x-direction) and in the absence of an external magnetic bias field. Adjacent to the waveguide is a cobalt–iron–boron (CoFeB) disk, operating as a spin-transfer vortex oscillator (STVO), which can be driven to auto-oscillation if an electric current is passed through it. The eigenfrequencies of the two sub-systems can be designed to match the desired functionality of the hybrid structure.
Schematic representation of the hybrid magnonic-oscillator system (HMOS). The uniform magnetization along the waveguide is pointing to the right. As an example, the vortex is shown here with clockwise chirality and negative polarity. To illustrate the dipolar coupling, the dipolar stray fields are shown as emanating from a magnetic dipole at the position of the oscillator.
Schematic representation of the hybrid magnonic-oscillator system (HMOS). The uniform magnetization along the waveguide is pointing to the right. As an example, the vortex is shown here with clockwise chirality and negative polarity. To illustrate the dipolar coupling, the dipolar stray fields are shown as emanating from a magnetic dipole at the position of the oscillator.
On the one hand, the gyrotropic frequency of the vortex oscillation () for thin disks in the absence of an external magnetic field depends on the thickness (), the radius () of the disk, as well as on the saturation magnetization .47 For , is approximated as
On the other hand, the dispersion relation for spin waves in this waveguide geometry is given by48
where = , = , is the exchange length, is the exchange constant, is the saturation magnetization, is the total in-plane wave number with its component along the propagation direction, and is the effective wave number across the waveguide width. The tensor describes the dynamical magneto-dipolar interaction, is the static internal magnetic field, is the vacuum permeability, and is the gyromagnetic ratio.
Thus, besides selecting particular materials for the oscillator and the waveguide, one can tune the individual resonance frequencies by changing the thickness and width of the waveguide, as well as the diameter and thickness of the disk to reach a resonant coupling between the STVO and the magnonic waveguide. To ensure realistic conditions, we have used standard parameters of nanometer-thick YIG films49 and circular nanopillars with diameters of 100 nm for the oscillator. At the oscillator, the magnetic configuration of the free layer is a vortex, with the possibility to tune the core polarities and chiralities.50,51 We focus here on the gyrotropic oscillation mode because the gyration of a vortex core through spin-torque in a single pillar is experimentally well proven and understood.52 Analytical descriptions of these dynamics have been shown to be quantitatively consistent with experiments.53 In addition, vortex oscillators exhibit low phase noise and have been experimentally shown to synchronize through dipolar coupling.54 In the numerical simulations, we assume that the magnetization of the polarizing magnetic layer is fixed and points out-of-plane and that the magneto-static field emanating from it is negligible. The magnetic parameters used for the simulations are displayed in Table I.
The parameters for micromagnetic modeling of the HMOS in MuMax3.
Parameter . | Disk . | Waveguide . |
---|---|---|
Msat (A/m) | 1.4 × 106 | 1.4 × 105 |
Aex (J/m) | 1.6 × 10−11 | 3.5 × 10−12 |
α | 2 × 10−3 | 2 × 10−4 |
Parameter . | Disk . | Waveguide . |
---|---|---|
Msat (A/m) | 1.4 × 106 | 1.4 × 105 |
Aex (J/m) | 1.6 × 10−11 | 3.5 × 10−12 |
α | 2 × 10−3 | 2 × 10−4 |
In the waveguide, the Gilbert damping parameter is set to a value of . In addition, to avoid spin-wave reflections, absorbing boundary conditions were modeled by exponentially increasing at the ends of the waveguide to a value of 0.5. The mesh size was set to (two layers in the z-direction). The waveguide has a total length of , a width of and a thickness of . The oscillator was placed in a distance of nm from the waveguide. The oscillator has a thickness of nm, and a disk diameter of . To excite the vortex oscillation via the spin-transfer torque effect, a DC current injected through the disk in out-of-plane direction is modeled (Simulation parameters: Spin polarization , degree of non-adiabaticity , Slonczewski secondary STT term = 0). We note that in the simulation, only this spin-polarized current injected into the free layer is explicitly modeled. In the experiment, such a current is generated, for example, by a fixed magnetic layer, which is not explicitly modeled in our simulations. Also, the Oersted field generated by the charge current is computed analytically and plugged into Mumax. The y-component of the magnetization of each cell was collected over a time period of 200 ns and recorded in intervals of 0.08 ns. The spin-wave dispersion was calculated by performing a fast Fourier transform (FFT) in space and time of in the waveguide, and the frequency of the oscillator is extracted from the fast Fourier transform of in the disk. Here, it is worth noting that we chose the parameters of CoFeB for the oscillator since this material has a high saturation magnetization which increases the frequency of the oscillator pushing it closer to the spin-wave dispersion in the YIG waveguide.
III. RESULTS AND DISCUSSION
We illustrate the waveguide dispersion relation using broadband excitation by means of an antenna in the center of the waveguide. The resulting mode spectrum of the waveguide is shown in Fig. 2(a). It contains two spin-wave modes. The low frequency mode is the fundamental waveguide mode with uniform oscillation across the width, while the higher frequency mode is the first width mode of the waveguide. This result reflects a typical waveguide dispersion curve for a width confined nanostripe geometry.48 For comparison, the fundamental mode was calculated analytically (green curve). Figure 2(b) shows the spin-wave spectrum in the waveguide, when the vortex oscillator is driven to auto-oscillation by a constantly injected spin current. The vortex oscillation couples dipolarly to the waveguide and excites spin waves at a frequency of GHz. Figure 2(c) shows the FFT of the component of the disk, showing that the auto-oscillation frequency of GHz is equal to the frequency of the emitted spin waves. Adjusting the current allows this frequency to be changed, as will be shown later.
Dispersion relation of spin waves propagating in the waveguide: (a) spin-wave spectrum excited by an antenna by applying a sinc-function field at the middle of the waveguide. The two lowest waveguide modes are visible. The green line shows the analytical dispersion relation for the uniform mode which has been calculated using the same parameters as in the simulation55 showing excellent agreement. (b) Spin waves excited by the spin transfer oscillator based on the vortex state (the electrical current density is A/m). One recognizes a dominant excitation of spin waves with propagating to the right at a frequency of . (c) FFT spectra of of the oscillator showing that . (d) Snapshots of the spatial profile of (upper subfigure) and integrated over the waveguide width (lower subfigure), as a function of the distance x from the oscillator along the waveguide, respectively.
Dispersion relation of spin waves propagating in the waveguide: (a) spin-wave spectrum excited by an antenna by applying a sinc-function field at the middle of the waveguide. The two lowest waveguide modes are visible. The green line shows the analytical dispersion relation for the uniform mode which has been calculated using the same parameters as in the simulation55 showing excellent agreement. (b) Spin waves excited by the spin transfer oscillator based on the vortex state (the electrical current density is A/m). One recognizes a dominant excitation of spin waves with propagating to the right at a frequency of . (c) FFT spectra of of the oscillator showing that . (d) Snapshots of the spatial profile of (upper subfigure) and integrated over the waveguide width (lower subfigure), as a function of the distance x from the oscillator along the waveguide, respectively.
Figure 2(d) shows the propagation of spin waves along the waveguide in the respective positive direction (to the right), as well as the distribution of the magnetization along the longitudinal axis. The dominant propagation of spin waves in the positive direction as shown in Fig. 2 for the studied vortex chirality and waveguide magnetization introduced in Fig. 1 indicates that the emission is unidirectional and that it is possible to control the generated spin waves by simply changing the two degrees of freedom of the vortex, i.e., its polarity and chirality.
In order to demonstrate the reconfigurability of the emission direction and amplitude, we present results of micromagnetic simulations for different vortex configurations in Fig. 3 for a fixed injected charge current density of . All four different combinations of vortex chirality/polarity are displayed. Figure 3(a) refers to the case of a vortex with clockwise chirality and a positive polarity (CP). In this case, the vortex core gyrates counterclockwise . The generated spin waves propagate to the left with a relatively small amplitude. When changing only the polarity (CP) and the direction of the current (required to keep the vortex auto-oscillation with STT), the vortex core gyrates clockwise and the spin waves propagate to the right with an amplitude that is nearly equal to that of case (a) [see Fig. 3(b)]. In both cases, the spin waves have a frequency of GHz and a wave vector rad/m. When changing the chirality, as displayed in (c) (CP) and (d) (CP), the emitted spin waves have a larger amplitude compared to cases (a) and (b), yet the directional dependence on the polarity is the same. However, the frequency and wave vector of the emitted spin waves in the cases (c) and (d) are different GHz and rad/m.50 It can thus be noted that changing the vortex polarity leads to changing the direction of the emitted spin waves along the waveguide. On the other side, the direction of the chirality relative to the waveguide magnetization affects the amplitude of the generated spin waves.
Two-dimensional maps of spin-wave emission showing different vortex polarity and chirality combinations and their respective effect on the direction and amplitude of the emitted spin-waves. Four different cases of different polarity and chirality combinations were presented by combining the polarities P [P: upward (a) and (c) and P: downward (b) and (d)] and the chiralities (C: clockwise and C: counterclockwise). The figure shows the out-of-plane spatial distributions of and line profiles represent integrated over the waveguide width along the longitudinal axis (lower sub-figure) , as obtained from simulation for fixed times. The injected current density is A/m for P and A/m for P.
Two-dimensional maps of spin-wave emission showing different vortex polarity and chirality combinations and their respective effect on the direction and amplitude of the emitted spin-waves. Four different cases of different polarity and chirality combinations were presented by combining the polarities P [P: upward (a) and (c) and P: downward (b) and (d)] and the chiralities (C: clockwise and C: counterclockwise). The figure shows the out-of-plane spatial distributions of and line profiles represent integrated over the waveguide width along the longitudinal axis (lower sub-figure) , as obtained from simulation for fixed times. The injected current density is A/m for P and A/m for P.
A schematic representation demonstrating the directions and amplitudes of the emitted spin waves throughout the HMOS is presented in Fig. 4. The effect of the vortex chirality can be understood by considering the characteristics of dipolarly coupled magnonic nanostructures. Different chiralities lead to an anti-parallel or parallel orientation of the local magnetization of the oscillator and the waveguide. Since the dipolar coupling is stronger for the anti-parallel alignment,55 the emitted spin-wave amplitude is larger in the presented case for a counterclockwise vortex chirality. The situation is inversed if the static magnetization in the waveguide is inversed.
Schematic representation of the HMOS emission behavior, based on the results of the micromagnetic simulations. While the polarity P controls the emission direction, the chirality C defines amplitude and frequency of the emitted spin waves.
Schematic representation of the HMOS emission behavior, based on the results of the micromagnetic simulations. While the polarity P controls the emission direction, the chirality C defines amplitude and frequency of the emitted spin waves.
A noticeable characteristic of the generated spin waves is their unidirectionality which appears despite the total reciprocity of the spin-wave dispersion for the presented waveguide geometry.56 As seen in Fig. 3, the generated spin waves propagate in the waveguide in a single direction with the maximum amplitude, and an almost vanishing amplitude in the opposite direction. To illustrate the origin of this complex behavior, we performed a vortex-only simulation to obtain the dynamic fields generated by the rotating vortex at the position of the waveguide. In Fig. 5, we present the amplitude of , which is the main component of the dynamic dipolar field exerting a torque on the wave guide’s magnetization, as a function of its position along the waveguide. Displayed is the evolution of for half of a vortex gyration cycle ( to ). Figure 5(a) refers to the case of the vortex with positive polarity which results in a counterclockwise vortex core gyration (the sense of the vortex core gyration must not be confused with the chirality C of the static vortex) . In this case, the spin waves are emitted into the waveguide in the negative direction of the x axis (refer back to Fig. 3). As can be seen from Fig. 5(a), the dynamic dipolar field is not similar to a standing wave but it has characteristics of a wave moving from right to left. This can be understood considering the direction of the vortex core movement when the latter is closest to the guideline in its upper half of the oscillation circle. Consequently, the direction of the vortex core movement when close to guideline also agrees with the direction of the spin-wave emission. Figure 5(b) presents the alternative case of a vortex with negative polarity where the core’s gyration trajectory is clockwise . Here, the spin waves are emitted into the waveguide in the positive direction of the x axis (refer back to Fig. 3). Due to the opposite sense of the vortex core gyration, the dipolar field has characteristics of a wave moving from left to right () which again agrees with the vortex core movement when close to guideline.
Dynamic magnetic field generated by a vortex at the location of the waveguide. The dynamic field’s strength is presented for the two cases (a) vortex with positive polarity, (b) vortex with negative polarity as a function of its waveguide position for different oscillation periods (). The respective generated spin waves traverse the waveguide in opposite directions for /. For more clarity, see the attached animation Fig5annimation.mp4. in the supplementary material.
Dynamic magnetic field generated by a vortex at the location of the waveguide. The dynamic field’s strength is presented for the two cases (a) vortex with positive polarity, (b) vortex with negative polarity as a function of its waveguide position for different oscillation periods (). The respective generated spin waves traverse the waveguide in opposite directions for /. For more clarity, see the attached animation Fig5annimation.mp4. in the supplementary material.
Finally, we investigate the spin waves’ characteristics as a function of the current density (in A/) (see Fig. 6) for (magnetization pointing along the positive axis of the waveguide. The vortex polarity and chirality are set as (CP) to ensure strong spin-wave emission into the positive x-direction. In Fig. 6(a), the evolution of the frequency and wave vector are shown to increase from to GHz and from to rad/m, respectively, as a function of ranging from to A/m. This well-known increase in the system’s frequency is accompanied by an initial increase in the oscillator’s amplitude as well as an increase in the intensity of the generated spin waves [see Fig. 6(b)]. With the steady increase in frequency, the vortex gradually reaches a state of nearly constant amplitude at a frequency around 2.35 GHz (current density of A/m). Beyond this point, the reached amplitude is then maintained by the oscillator despite the increasing frequency and current density since the gyration of the vortex has reached the largest possible radius and a further increase is suppressed by nonlinear effects. In this case, the increase in frequency translates to an increase in the vortex core’s gyration speed.57 Interestingly, the intensity of the generated spin waves drops at the time the oscillator reaches a constant amplitude for frequencies exceeding 2.45 GHz. The energy flow carried by the waves is proportional to the product of intensity and the wave group velocity. Concerning the emitted energy, the drop in intensity is at least partially compensated by an increase in the group velocity of the spin waves [see Fig. 6(c)]. Figure 6(d) is a plot of the spin waves’ non-reciprocity C as a function of , where C represents the ratio between the spin wave intensities emitted to the right and to the left. It gradually increases for increasing current density at values beyond that at which vortex dynamics are excited. This trend could be related to the decrease of the emitted wavelength with increasing current which changes the dipolar coupling of the wave to the oscillator, but further studies are needed to quantify this effect.
Characteristics of spin-wave emission as a function of current density (in A/) injected in the oscillator for the case of a vortex with negative polarity and chirality (CP). (a) Evolution of the frequencies in GHz (red) and wave vector in rad/m (blue) as a function of . (b) Evolution of the spin-wave intensity (red) and oscillator amplitude (blue) as a function of . (c) The spin-wave group velocity as a function of . (e) Non-reciprocity as a function of .
Characteristics of spin-wave emission as a function of current density (in A/) injected in the oscillator for the case of a vortex with negative polarity and chirality (CP). (a) Evolution of the frequencies in GHz (red) and wave vector in rad/m (blue) as a function of . (b) Evolution of the spin-wave intensity (red) and oscillator amplitude (blue) as a function of . (c) The spin-wave group velocity as a function of . (e) Non-reciprocity as a function of .
IV. SUMMARY AND CONCLUSIONS
The continuous development of spintronic devices for information processing technologies aims to design and realize spintronic-based hardware building blocks with large nonlinearity for unconventional computing concepts like neuromorphic computing. Here, we propose a new hybrid device realizing tunable spin-wave emitters which are characterized by merging two subfields, namely, the field of spintronic oscillators and the field of nano-magnonics. The main advantages of the hybrid magnonic-oscillator system (HMOS) is its nanoscale size, its operation with simple DC currents and its potential to use wave-based effects like frequency-multiplexing to establish efficient connections. We propose a system that utilizes the spatial dynamics of dipolar fields generated by a vortex-based spin transfer oscillator. Coupling these fields to a magnonic waveguide leads to the emission of spin waves. We have verified the functionality of the HMOS via micromagnetic simulations by choosing a set of parameters, including different materials, the element sizes and the gap. It was shown that the HMOS can be used as tunable source for the emission of coherent, propagating spin-waves. This proof-of-concept idea was modeled using simple geometries, no requirement for an external field, and parameters of commonly available materials. Many further modifications can be obtained, e.g., by changing the external magnetic field, the geometric parameters or the static direction of magnetization of the waveguide. Moreover, the wave vector, as well as the frequencies and the amplitudes of the spin waves are tuned by the current density injected into the spin transfer oscillator. In addition, the hybrid system offers the possibility of tuning the direction of emission by the sense of the vortex gyration given by the vortex polarity. The strength of the emission can be changed by the vortex chirality. Combining all these features together, we believe that this device will open a way for the development of new and efficient signal-processing devices for unconventional computing.
SUPPLEMENTARY MATERIAL
See the supplementary material for the dynamic magnetic field generated by a vortex at the location of the waveguide for the two cases of polarities.
ACKNOWLEDGMENTS
This research was funded by the European Research Council within the Starting Grant No. 101042439 “CoSpiN” and by the Deutsche Forschungsgemeinschaft (DFG, German Research Foundation) within the Transregional Collaborative Research Center TRR 173268565370 [Spin X(project B01).]
AUTHOR DECLARATIONS
Conflict of Interest
The authors have no conflicts to disclose.
Author Contributions
A. Hamadeh: Conceptualization (equal); Formal analysis (equal); Investigation (equal); Supervision (equal); Validation (equal); Visualization (equal); Writing – original draft (lead); Writing – review & editing (equal). D. Breitbach: Conceptualization (supporting); Formal analysis (supporting); Writing – review & editing (equal). M. Ender: Software (equal); Visualization (equal); Writing – review & editing (supporting). A. Koujok: Software (equal); Writing – review & editing (supporting). M. Mohseni: Writing – review & editing (supporting). F. Kohl: Writing – review & editing (supporting). J. Maskill: Writing – review & editing (supporting). M. Bechberger: Writing – review & editing (supporting). P. Pirro: Conceptualization (supporting); Formal analysis (supporting); Funding acquisition (lead); Investigation (supporting); Supervision (lead); Validation (supporting); Writing – review & editing (equal).
DATA AVAILABILITY
The data that support the findings of this study are available from the corresponding author upon reasonable request.