Twisted magnon beams carrying orbital angular momentum

Low-energy eigenmode excitations of ferromagnets are spin waves or magnons that can be triggered and guided in magnonic circuits without Ohmic losses and hence are attractive for communicating and processing information. Here we present new types of spin waves that carry a definite and electrically controllable orbital angular momentum (OAM) constituting twisted magnon beams. We show how twisted beams emerge in magnonic waveguides and how to topologically quantify and steer them. A key finding is that the topological charge associated with OAM of a particular beam is tunable externally and protected against magnetic damping. Coupling to an applied electric field via the Aharanov-Casher effect allows for varying the topological charge. This renders possible OAM-based robust, low-energy consuming multiplex magnonic computing, analogously to using photonic OAM in optical communications, and high OAM-based entanglement studies, but here at shorter wavelengths, lower energy consumption, and ready integration in magnonic circuits.

R ecently, the spatial structuring of photonic, electronic, or neutron beams has been demonstrated 1-5 enabling to encode additional information in the beam spatial distribution. For instance, a photonic beam can be spatially modulated to have a helical or twisted wavefront embodying a welldefined, internal, meaning origin-independent, orbital angular momentum. The external orbital angular momentum (OAM) is origin dependent and is obtained from the photonic angular momentum density L γ = r γ × P γ , where r γ is measured with respect to the beam center and P γ is the linear momentum density 6 . Similar arguments apply to matter fields [3][4][5] . The intense research on this topic is fueled not only by fundamental interest but also by the new possibilities and functionalities accomplished by structuring the spatial distributions of the waves. For instance, OAM can take very large values that can be exploited for multiplex communications and quantum information based on entangled photons with large OAM [1][2][3][4][5][7][8][9] . Here we uncover the existence of twisted magnon beams, which are low-energy spin waves that carry a well-defined OAM. We demonstrate how they are triggered in a magnonic waveguide (such as in Fig. 1), and how to quantify their topological nature by an associated topological charge, which is found to be related to the OAM of a particular twisted beam. It is shown that this topological charge is electrically tunable and protected against damping. Hence, twisted magnonic beams are ideal candidates for robust, low-energy cost multiplex magnonic computing using the versatile toolbox for material and spin waves engineering [10][11][12][13][14][15][16][17] . In addition, we find an OAM-dependent magnonic Hall effect that depends on the value of OAM and, hence its value can be steered electrically. We also inspect the nature of thermal magnetic fluctuations around the equilibrium-ordered phase in a magnetic sample supporting twisted modes and find an anisotropic dependence of the correlation length.

Results
Analytical model for twisted magnon beams. What are the requirements for the emergence of OAM-carrying waves in magnetic systems and what is their relevance? In magnetically ordered systems key low-energy carriers of information are spin waves or magnons, which are the quanta of low-energy magnetic excitations. Identifying twisted magnonic beams that carry OAM would open a subfield in magnonics since OAM in principle can have an unbounded value. Thus, let us consider a generic ferromagnetic (FM) cylindrical waveguide modeled by localized magnetic moments at site i described by the magnetic moment operator M i or the corresponding (dimensionless) spin operator S i = −M i /(γħ) (γ is the gyromagnetic ratio). The (Heisenberg) Hamiltonian reads is the exchange coupling between nearest-neighbor sites (ij), and K s > 0 is a magnetic anisotropy contribution including the shape anisotropy (the demagnetizing factor tensor D ij at any given direction of the cylindrical nanowire are approximately zero, except for D xx = D yy ≃ 1/2 18,19 ). B = (0, 0, B z ) is an external magnetic field along the waveguide (this direction is taken as the zaxis). For magnetically ordered systems longitudinal excitations, meaning changes in the expectation value of M i , caused for example by changes in J ij , are much higher in energy than transversal ones. The latter correspond to a precession of the unit vector field m i = M i /|M i |, while |M i | = M = constant and their energy is typically set by K s , which is much smaller than J ij 20 . Hence, we inspect the transversal spin dynamics governed by the Heisenberg equation of motion dm i =dt ¼ Ài= h½m i ; H. The dynamics of the magnetization unit vector field proceeds as (we introduce the We are interested in weak transversal fluctuations around the ordered magnetic phase (weak means slow variation in m i within a lattice constant a). Furthermore, we will consider sufficiently large spins where quantum/thermal fluctuations are subsidiary and magnetic reversals are rare. Thus, it is useful to go over to a continuous magnetic texture vector field m i (t) → m(r i , t) and interpret m(r, t) as the local, time-dependent expectation value of the field operator m(r, t). The exchange contribution in Eq.
(1) appears as the divergence of the spin current 21 , where ψ(r, t) = m x (r, t) − im y (r, t) and where p = −iħ∇ is the momentum operator and m ? J ¼ h=ð2γJMa 2 Þ is an effective mass. The wave function ψ(r, t) vanishes outside the waveguide and, otherwise, is subject to the potential V(r) = −ħγ(K s M + B z ). In an extended cylindrical waveguide and in cylindrical coordinates r → (r, ϕ, z), Eq.  Fig. 1 Proposal for generating a twisted magnon beam exploiting the Aharanov-Casher effect. Incident spin waves (impinging from left) that propagate in a cylindrical magnetic insulator waveguide traverse a region with a linear charge density λ e (marked yellow). As demonstrated here, a magnonic twisted beam emerges to the right and can be quantified by the topological charge ' / λ e , which is associated with a monopole-like vector potential, A e ðrÞ $ λ e r b e ϕ that is intimately related to the Dirac phase for magnons ref. 22 ) ψðr; tÞ / J ' ðk ? rÞ expði'ϕ þ ik z zÞ expðÀiωtÞ; ð5Þ with J ' ðxÞ is the Bessel function of the first kind with order '. The total energy is Note that similar to optics, conventional (meaning ' ¼ 0 in Eq. (5)) Bessel modes 22,23 also exist, but those do not carry any OAM. A striking consequence is that Bessel beams do not interact with electric fields, in contrast to twisted magnon beams and also do not allow for '-based operations. This fact is decisive when it comes to constructing OAM-sensitive magnonic circuits (similar to ref. 11 ) that are controlled by electric fields offering so new functionality, as discussed in ref. 24 . A further important aspect is the difference between our propagating OAM modes and static/ driven magnetization vortices [25][26][27] . The latter are localized topological excitations that are much higher in energy than a magnon wave (cf. Eq. (5)), which extends over the whole waveguide. The much lower energy set by E(k) is tunable by J; m ? J and V. Both magnetization vortices and propagating, OAM-carrying magnons have well-defined topological features that can be quantified. To do that in our case, we consider the z-component of the magnon current J m z , which is nothing but the probability current density coiling around the z-axis meaning that the spin waves spiral around the axis of the waveguide, similarly to the case of photon, electron, or neutron twisted beams [1][2][3][4][5] . In fact, the analogy runs deeper: let us introduce pseudo-electric and -magnetic fields as e E ¼ m x þ im y and e B ¼ iðm x À im y Þ: Physically, e B and e E describe the (paraxial) spin wave modes rotating clockwise and anticlockwise, respectively. The magnon density is W m ¼ 1 2 j e Ej 2 þ j e Bj 2 À Á , and the m z polarized magnon current density J m z can be related to the pseudo-Poynting-like vector 28 It is straightforward to infer that P has the proper time and space symmetry of a canonical momentum. Noting that our twisted beam is localized within the waveguide, we can introduce the canonical OAM of the magnon twisted beam as L ¼ r P, which is extrinsic and dependent on the choice of the coordinate origin. However, the integral orbital angular momentum over the cross-section is indeed intrinsic, amounting to with k ' k z b e z being the mean wave vector of the spin wave. We have then hLi Á hPi=P ¼ ', in full analogy with paraxial optical twisted beams 28 . This ' can be identified with the topological charge associated with a particular twisted beam. Generally, our magnons are characterized by the integral momentum hP i and the orbital angular moment 〈L〉 as two independent properties, meaning that in addition to the intrinsic spin angular moment ħ along the b e z , the magnon in the cylindrical nanowire may carry a longitudinal and intrinsic orbital angular moment L z $ h' as well, which formally can take arbitrarily large values and can be utilized in magnonic and spintronics applications.
Undamped topological charge of twisted beams. An important feature of twisted magnons is the robustness of the associated topological charge to magnetic damping. This fact is crucial when it comes to utilization of the OAM of twisted magnons as information carriers, analogously to OAM-based optical communications 7,8 . Generally, the (Gilbert) damping of magnetization procession 11,17 is governed by a term of the form αm × ∂m/∂ t . In Eq. (3), this amounts to substituting ∂ t → (1 + iα)∂ t , which leads to the time-decaying magnon density ρ ' ðtÞ ¼ ρ ' ðt ¼ t 0 Þ expðÀ2αωðt À t 0 ÞÞ when propagating from time t 0 to t. The magnon OAM ' derives, however, by an averaging of 〈L〉 over the cross-section of the waveguide S, Thus, damping leads to a decaying magnon density (its photon analog is the light intensity), which is however independent of ' and hence of OAM. In other words, the OAM (') is not affected by damping and is conserved during the magnetization evolution. Another point of view is that the damping term would dominate the interfacial spinpumping current by m ∂m=∂t / = ψðr; tÞ∂ t ψ ? ðr; tÞ ½ ¼ ωρ ' , which clearly does not depend on the internal phase structure of the twist excitation (and hence on '). These conclusions are supported by the full-fledge numerical simulations presented below.
Simulating an experimental realization. Having laid down the principle existence of twisted OAM and the robustness of the associated topological charge to perturbations that subsume into damping of the precessional dynamics, the experimental feasibility as well as the role of further magnetic interactions that emerge in a realistic setting, have to be clarified. Therefore, we resort to full-fledge numerical micromagnetic simulations (see Methods) for an experimentally realistic waveguide made of insulating yttrium iron garnet (YIG).
The cylindrical YIG waveguide has a diameter of 0.4 μm and a length of 2.0 μm. The saturation magnetization is M s = 1.4 × 10 5 A/m, the exchange stiffness constant is A ex = 3 × 10 −12 J/m, the uniaxial magnetocrystalline anisotropy along z is K u = 5 × 10 3 J/m -3 , and the Gilbert damping α = 0.01 24 . The ground state equilibrium magnetization is oriented along the +z-axis, parallel to the main axis of the waveguide. To launch the spin waves, we excite locally at one end z = 0 with a twisted magnetic B field (see Methods) that is linearly polarized in the x-direction and having the amplitude B max = 10 mT and the frequency ω B = 5 GHz. Figure 2a shows a snapshot of the x-component of magnetization excitations propagating after 2 ns. More details of excitation mode are gained from slices along the tube in Fig. 2b, which demonstrates the generation of propagating twisted magnon beams. An animated version of m x at the transversal slice z = 1 μm during the magnon propagation is provided as Supplementary Movie 1. The finite damping leads to a decaying magnon density along the YIG tube with the amplitude being normalized with respect to the maximum excitation at z = 0. Remarkably, when evaluating the OAM ' from the numerically calculated pseudo-Poynting vector, we indeed find that ' is independent of the magnon density and is conserved during the spin wave propagation along the tube (cf. Fig. 2c). These full numerical calculations endorse the above formal analysis as well as the experimental feasibility and relevance of twisted magnon beams.
Aharonov-Casher effect and Landau levels. The simulations shown in Fig. 2 as well as our formal analysis evidence that the twisted magnon beams embody a non-vanishing z-polarized spin current density allowing for a coupling to an external electric field E through the Aharonov-Casher (AC) effect 29 . For magnonbased information transfer, this fact is crucial as it allows to control the flow of OAM-carrying magnons with electric means. In the presence of an electric field, the exchange interaction in the Heisenberg Hamiltonian becomes anisotropic according to where L j'j n are the associated Laguerre polynomials, and n = 0, 1, 2, … is the quantized radial number. The Landau energy levels read Here w e ¼ 2 ffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffi h=q m jBj p is a characteristic width that depends on the amplitude of the applied electric field, Ω e ¼ q m σjBj=2m ? J is the Larmor frequency (σ = ±1 correspond to the clockwise or anticlockwise rotation, respectively). Thus, the electric field E confines the magnons to form paraxial LG beams with The azimuthal component is directly related to the kinetic orbital angular momentum L ¼ r P e as hL z i ¼ À hhi∂ ϕ þ q m rA e i.
Unlike the Bessel beams with uniform rotational direction described by the topological charge ', the direction of hL z i of the LG beams depends on ' and the twist vector potential A e as well. Note, the hL z i radial structure depends strongly on σ, that is, the direction of the applied electric field in terms of the first radial maximum of the LG mode r j'j ¼ w e ffiffiffiffiffiffiffiffiffi ffi j'j=2 p : For 'σ>0, we have a uniform, rotationally azimuthal current as usual. However, for 'σ<0 the currents in the region of r<r ' and r>r ' are countercirculating and are dominated by the topological charge ' and the twist vector potential A e , respectively. For an estimate, we take: JM 2 = 0.1 μeV, a = 10 Å, and E ¼ 1 V/nm, which give w e ≈ 1.4 μm and ħ|Ω e | ≈ 0.33 μeV. Increasing the strength of the applied electric field, higher electric Larmor frequencies Ω e are achievable but the electric Landau radius w e is smaller.
OAM-tunable magnon Hall effect. Having clarified the role of an external electric field, one may wonder how the twisted magnon beam is affected by magnetic perturbations. A hallmark of a charge wave-packet propagation in a magnetic field is the transverse deflection, that is, the Hall effect, which in its simplest variant is assigned to the Lorentz force. A similar effect for uncharged particles, such as photons, phonons, or magnons, may thus appear surprising. In our FM case, on the other hand, the Hall effect is also sensitive to the residual magnetization, a phenomenon termed the anomalous Hall effect. Indeed, the Hall effect for magnons has already been observed 34 . Is this effect OAM dependent and can we control it by tuning the OAM? A positive answer would be an important advance, since OAM can be controlled externally, as demonstrated in Fig. 2. To approach this issue, we note that LG modes have well-defined azimuthal and radial wavefront distributions (quantified by ' and n, respectively) and they form an orthogonal and complete basis in terms of which an arbitrary function can be represented. So let us consider a magnonic wave packet centered at (p c , r c ) (similar to electronic wave packets 35,36 ). From the above, we conclude that the trajectory spirals around the z-axis. Projecting the local coordinate frame along the curved trajectory, we end up with a noncommutative geometry and a covariant derivative in momentum space [37][38][39][40] : ;n j∂=∂ p jψ L ';n i is the Berry connection. The corresponding Berry gauge field (curvature), B g ¼ ∂=∂ p A g ðpÞ gives rise to a deflection of the wave packet center. Given that the z-axis is now locally directed along p, the vector ðb e x þ ib e y Þ is orthogonal to p and moves on the unit sphere with p/p under variation of p, which results in a magnetic-monopole-type field 41,42 , B g ¼ À'p=p 3 for the orbital motion of magnons (since expði'ϕÞ ¼ ðb e x þ ib e y Þ ' ). As a result, we infer the Hamiltonian equations of motion 35,39 dr The anomalous velocity term h' _ p c p c =p 3 is perpendicular to p c causing a transversal motion of the wave packet. The corresponding transverse shift δr c is cast as δr c ¼ h' R L ðp c dp c Þ=p 3 c and is determined by the geometry of the contour L in momentum space and '. This is in so far remarkable, as arbitrarily large values of ' are possible giving rise to a large deflection at small driving force. For concreteness, let us consider a weak spatially varying, perpendicular magnetic field B z , say along the x-direction 43,44 . The driving force is F x ¼ À∂ x B z ðrÞb e x and thus _ p c ¼ F x and δ_ r c ¼ h'F x p c =p 3 c . Assuming paraxial magnonic propagation, meaning p c ' p z b e z , the transverse deflection reads (an exact but cumbersome analytical expression is available) δr c ' δy / h'F x z p 3 c b e y : Different LG modes with different topological charges ' split at the waveguide end (cf. Fig. 3a). Similar to photon 41 and electron 42 twisted beams, we predict a magnonic OAM-dependent Hall effect, which is qualitatively different from the usual magnon Hall effect 10,34 . For a cylindrically inhomogeneous waveguide, for instance due to inhomogeneity in the demagnetization field 18,19 , we may formally write the potential in Eq. (4) as VðrÞ The displacement direction is along b e ϕ , meaning that LG modes are rotated uniformly clockwise or anticlockwise depending on ' (cf. Fig. 3b), a result in line with the optical Magnus effect 6,45 .
Another aspect arises when the waveguide is an FM metallic wire, in which case the electric field vanishes inside the wire. However, we may still have nontrivial magneto-electric effects. Starting from a general symmetrical but inhomogeneous electric potential V e (r), the energy shift due to the AC effect is given to the first order in the perturbation V e (r) by δE ' q m m ?
∂r i: For a step-function profile of the wire electric potential V e (r) = UΘ(r − R) (here Θ(r) is the Heaviside function and R is the radius of the metallic wire), we find ∂V e /∂r = Uδ(r − R). No electric field is present in the nanowire except for the surface. Nonetheless, the twisted magnon beam acquires an additional phase upon traveling a distance z along the nanowire φð'; zÞ ¼ Uρ ' ðRÞ q m hc ' k z z R : This phase accumulation may result in an interference between the different LG modes, for example, ψ À' þ ψ ' À Á / cos½j'jϕ þ φðj'j; zÞ.
Anisotropic correlation length of magnetic fluctuations. While the control of twisted magnons via external electromagnetic fields is crucial for applications, one may wonder about the stability of the FM state to twisted excitations. To explore this aspect we have to study the influence of vorticity on the magnetic fluctuations and in particular the correlation length. Both are indicators on how magnetic ordering reacts to excitations 46 . Let us start with the Ginzburg-Landau free energy density of the FM waveguide at temperature T (in unit of 1/M 2 ), F ¼ Ja 2 2 ð=mÞ 2 þ ðk B T À K s 2 Þm 2 z þ K 4 m 4 À B z m z (with k B is Boltzmann constant and |m z | ≤ 1). The quartic term with K 4 > 0 stabilizes the long-range FM ordering below T c ≈ K s /2k B . Because of the uniaxial anisotropy and the cylindrical demagnetizing field, the SU(2) spin symmetry is reduced to SO(2) Z 2 . Low-energy Goldstone modes are due to the continuous SO(2) symmetry.
For T < T c the saddle point approximation delivers the meanfield value of m z in the absence of an external magnetic field: For the paramagnetic region T > T c , the quartic K 4 term is not relevant and the mean-field equilibrium solution is m z ¼ B z =k B ðT À T c Þ. For transverse fluctuations around the saddle point solution m z , we write m ' ffiffiffiffiffiffiffiffiffiffiffiffi ffi 1 À ρ 2 p b e z þ ρb e r ; with 0 ≤ ρ ≤ 1. Substituting in δF and neglecting constants and higher order terms in ρ leads to δF ρ ¼ Ja 2 2 ð=ρÞ 2 À αðTÞρ 2 þ βρ 4 with α (T) = [k B T − (K s + B z )/2] and β = B z /8. Interestingly, transverse fluctuations are most probable not at the boundary ρ = 0 but at , the global minimum of the free energy δF ρ . For cylindrical waveguide the transverse fluctuations are subject to a non-zero restoring force and the correlation length is 1 Ja 2 q : Thus, the transverse scattering density is Lorentzian. Generation of twisted magnon beams. In Fig. 2 we demonstrated how to trigger twisted magnons via magnetic fields. In fact, there are various other ways to accomplish this goal: one may use, for instance, an engineered magnetic spiral phase plate or fabricate a hologram such as a pitch fork for scattering of incoming magnons before entering the waveguide. Another possibility is the following (cf. also Fig. 1): let us consider a uniformly charged thin wire with a linear charge density λ e located in the region z ∈ [d 1 , d 2 ] on the axis of a FM cylindrical waveguide. The radial electric field EðrÞ ¼ λ e 2πϵr b e r in the FM nanowire within z ∈ [d 1 , d 2 ] (ϵ is the electric permittivity of the FM). For z-polarized magnons, the interaction of magnons with the radial electric field E involves the vector potential A e ðrÞ ¼ ÀE b e z =c ¼ λ e 2πϵcr b e ϕ ; which is intimately related to the so-called Dirac phase for magnons moving along a contour C, C A e ðrÞ Á dr; producing a phase structure expði'ϕÞ with the quantized number ' ¼ À q m λ e 2π hϵc . Similar to twisted electron beams 47 , the twisted magnons can be produced by considering a transition of a conventional magnon with (' ¼ 0) in the region z ≤ d 1 without any flux (λ e = 0) to the region z ∈ [d 1 , d 2 ] with the non-zero flux (λ e ≠ 0). For an estimate, we have q m λ e 2π hϵc $ 5:6λ e 10 À15 with the linear charge density λ e in unit of e/m. For 10 22 cm −3 carrier density in copper and magnon wavelength (governed by the exchange interaction)~μm, the linear density λ ẽ 10 14 -10 16 e/m is experimentally feasible. Demonstration of OAM magnon beams evolution. As evidenced, twisted magnon beams are controlled in numerous ways including external electric and magnetic fields, geometry design and holographic means. For instance, with suitably designed holograms, different superpositions of twisted beams can be created such as the OAM-balanced superposition. This is determined by two LG eigenmodes with the same radial number n but opposite topological charges ± ', meaning ψ 1 ¼ ψ L À';n þ ψ L ';n : This superposition has zero net canonical OAM 〈L z 〉 = 0 exhibiting a flower-like symmetry pattern: jψ 1 j 2 / cos 2 ð'ϕÞ at z = 0, as shown in Fig. 4a. The OAM-unbalanced superposition consists of two LG eigenmodes with the same radial number n but different topological charge 0 and ' means that ψ 2 ¼ ψ L 0;n þ ψ L ';n : This superposition carries a finite net canonical OAM, 〈L z 〉 ∝ l, and is characterized by a pattern with ' off-axis vortices (cf. Fig. 4b).
In the paraxial approximation with a relatively small transverse kinetic energy and considering that Ω e / σjBj, we infer the allowed wave number as k z ≃ k + Δk z , with Δk z ¼ À½σ' þ ð2n þ j'j þ 1Þ=z e . Here z e ¼ w e ffiffiffiffiffiffiffiffiffiffiffiffiffiffiffi ffi E= hjΩ e j p is a characteristic longitudinal scale. On the propagation distance z, the electric field modifies the longitudinal wave vector resulting in an additional phase φ = Δk z z, which strongly depends on the direction (σ) of the applied electric field. As for 'σ>0, the induced phase φ shows a linear dependence on the topological charge '. In contrast, for 'σ<0, φ = −(2n + 1)z/z e is independent of '. These phenomena are caused by partial cancellation of counter-circulating azimuthal currents produced by the beam with expði'ϕÞ and by the twist vector potential A e .
The propagation of twisted magnon beams in the waveguide can give rise to a Faraday effect similar to the case of electron vortex states 32,33 . As for the OAM-balanced superpositions, the induced phases for the modes ψ L ± ';n reads φ ± ¼ Çσ'z=z e , which results in rotation of the interference pattern by the angle Δφ = σz/z e , as displayed in Fig. 4a. For the OAM-unbalanced superpositions, the mode ψ L ';n acquires an additional phase Δφ ¼ Àðσ' þ j'jÞz=z e , as compared to the ψ L 0;n mode. The superposition with σ'>0 shows a rotation of the interference pattern by the angle Δφ = 2σz/z e , as demonstrated by the lower panels in Fig. 4b with σ = −1, whereas the superposition with σ'<0, has no rotation at all (Δφ = 0, as evidenced by the σ = −1 panels in Fig. 4b.
All these electric field tunable phases, including the mixed cases with ψ 3 ¼ ψ L À';n þ ψ L 0;n þ ψ L ';n , as shown in Fig. 4c, indicate rich magneto-electric patterns of the interference intensity, which also may serve as a marker for identifying the twisted magnon beam.

Discussion
Imprinting an orbital angular momentum on a magnon beam puts a new twist on magnonics, as OAM can be functionalized as an independent robust, meaning damping resistant parameter. Together with the fact that OAM can be tuned to very large values, the twisted magnon beams offer new opportunities for reliable multiplex communications. The susceptibilities of the OAM-carrying magnons to electric and magnetic fields are thereby a key advantage, as they offer versatile tools to extract and steer the flow of information. Further immediate consequences of the topological nature of the discovered twisted beams are that the internal global phases, related to OAM, result in a controllable interference pattern when such beams interfere. This we may employ to imprint spatially modulated magnon density on a length scale well below the beam's extensions. The parameters of the twisted beams are directly reflected in the magnon density spatio-temporal pattern. Thus, in addition to the well-established, energy-saving applications of magnonics and magnon-based spintronics, this new offspring is useful for OAM-based, electrically controlled functionalities in magnonic-based information transfer and processing. Future studies are focused on the scattering of twisted magnon beams from a magnetization vortex core 48 and its possible application in a long distance (at short wavelength 49 ) electrically controlled multiplex communication channel between different vortex cores.

Methods
Micromagnetic simulations. We used the open-source, graphical processing unit (GPU)-accelerated software package MuMax 350 for the micromagnetic simulations. The (0.4 × 0.4 × 2.0) μm 3 large system is discretized into (100 × 100 × 500) cubic cells with a size of (4 nm) 3 . The spatial profile of the excitation field applied at the z = 0 layer reads B ' ðρ; ϕ; tÞ ¼ <êB max ð ρ w 0 Þ j'j exp j'j 2 1 À ρ 2 The waist parameter w 0 denotes the distance of maximum field strength B max ,ê corresponds to the polarization vector, which is equal toê x in our case, whereas w 0 = 75 nm. The time evolution is computed using a Dormand-Prince algorithm with adaptive time-steps. To avoid artifacts due to finite size simulations, we apply open boundary conditions (disregard demagnetizing fields).

Code availability
For obtaining the results plotted in Fig. 2 as well as the movie file, we used the opensource, GPU-accelerated software package MuMax 350 for the micromagnetic simulations. For reproducing the data, one has to use the paramter for YIG, as given in the text and discretize the (0.4 × 0.4 × 2.0) μm 3 large cylinder into (100 × 100 × 500) cubic cells with a size of (4 nm) 3 . Full details on the code are in ref 50 .