Magnon SpinMomentum Locking: Various Spin Vortices and Dirac Magnons in Noncollinear Antiferromagnets
Abstract
We generalize the concept of the spinmomentum locking to magnonic systems and derive the formula to calculate the spin expectation value for onemagnon states of general twobody spin Hamiltonians. We give nogo conditions for magnon spin to be independent of momentum. As examples of the magnon spinmomentum locking, we analyze a onedimensional antiferromagnet with the Néel order and twodimensional kagome lattice antiferromagnets with the 120 structure. We find that the magnon spin depends on its momentum even when the Hamiltonian has the axis spin rotational symmetry, which can be explained in the context of a singular band point or a symmetry breaking. A spin vortex in momentum space generated in a kagome lattice antiferromagnet has the winding number , while the typical one observed in topological insulator surface states is characterized by . A magnonic analogue of the surface states, the Dirac magnon with , is found in another kagome lattice antiferromagnet. We also derive the sum rule for by using the PoincaréHopf index theorem.
 PACS numbers

75.30.Ds,75.50.Ee,72.20.i, 75.76.+j,85.70.w
pacs:
Introduction.—The physics of magnons chumak (); magnonics (), the quanta of spin wave excitations, is enriched by their multiband nature. Although magnons have no internal degrees of freedom other than spin, such as atomic orbitals, the presence of chemical and magnetic sublattices allows magnonic systems to exhibit nontrivial band structures. In particular, recent studies have generalized many concepts in topological band theory hasan (); xlq (); niu (), established in multiband electron systems, to magnonic systems, e.g. the magnon Hall effect fujimoto (); katsura (); onose (); owerre3 (), magnon topological insulators shindou (); lzhang (); chisnell (); owerre (), and Weyl fyli () (Dirac fransson ()) magnons.
Another interesting feature is spin angular momentum carried by one magnon. Magnon spin in multiband systems, which has been studied for a long time kittel (), has attracted renewed attention rcheng2 (); zyuzin (); ohnuma (); rcheng () motivated by recent developments in spintronic techniques chumak (). Except for in simple collinear ferromagnets, magnon spin generally depends on its band properties. For instance, the magnon spin Nernst effect in antiferromagnets, which is a magnonic analogue of the spin Hall effect, is interpreted as the two copies of the magnon Hall effect for magnons with opposite spins rcheng2 (); zyuzin (). This example shows the importance of considering the spin direction of each magnon mode in magnon spintronics chumak ().
In this Letter, we generalize the concept of spinmomentum locking to magnonic systems. The conventional spinmomentum locking, in which electron spin depends on its momentum, is described by a noninteracting Hamiltonian without rotational symmetry in spin space, such as the Dirac Hamiltonian of topological insulator surface states hasan (); xlq (). We define the magnon spin for each band and prove a nogo theorem which states that magnon spin is momentum independent for several conditions. By performing a numerical calculation for kagome lattice antiferromagnets with a 120 structure, we find a spinmomentumlocked magnon band characterized by the winding number winding () [Fig. 1]. Our results for spin Hamiltonians with axis spin rotational symmetry demonstrate for the first time that spinmomentum locking can be generated through spontaneous symmetry breaking.
Definition of dependent magnon spin.—A general twobody spin interaction Hamiltonian is given by
(1) 
where is the spin operator at each site, denote the magnetic lattice vectors, denote the magneticsublattice indices, and is the number of sites in a magnetic unit cell. To rewrite the spin Hamiltonian (1) in terms of spin excitations (magnons) around a classical ground state, we introduce the HolsteinPrimakoff boson operators (),
(2) 
where is the size of the spin and is the set of the basis vectors of the rotating frame in which is in the direction of classical spin at . The corresponding quadratic form of boson Hamiltonian is given by
(3) 
where is the crystal momentum, , and is a bosonic Bogoliubovde Gennes Hamiltonian. We ignore magnonmagnon interactions in Eqs. (2) and (3) by assuming magmageffect (). The eigenenergy problem of Eq. (3) can be solved by the Bogoliubov transformation colpa (),
(4) 
where are paraunitary matrices and . () are the magnon field operators, which satisfy
(5) 
Using Eqs. (2) and (5), the total spin operator is given by
(6) 
Thus, the dependent magnon spin carried by a onemagnon state , where is the Fock vacuum of (), is given by
(7) 
Nogo conditions for ordered magnets.—In noninteracting electron systems, spinmomentum locking is forbidden for Hamiltonians with rotational symmetries rotsym (). In magnonic systems, however, it does not hold, because of spontaneous symmetry breaking. In the following, we write nogo conditions for magnon spin to be independent of momentum. We first consider twobody spin Hamiltonians with rotational symmetry in spin space such as isotropic Heisenberg models. Suppose that is a ground state with a spontaneous symmetry breaking. Using , we can rewrite in terms of up to the second order supplemental (),
(8) 
where () are field operators of massless NambuGoldstone (NG) modes associated with the spontaneous symmetry breaking, denotes the index of the independent massless NG modes, is the total number of the massless NG modes, and are complex numbers. For collinear magnets with the symmetry breaking: , is an eigenstate of the unbroken generator , which does not include (). A onemagnon state is also an eigenstate of :
(9) 
Using the notation that is an eigenstate of with an eigenvalue M instead of , we obtain
(10) 
where satisfy . Equation (10) shows that take quantized values. Because quantized spin components are not changed under a small momentum change , we cannot expect the dependent magnon spin in isotropic Heisenberg models with collinear ground states. Note that there is a exception. The above statement assumes the smoothness of on the magnon band . However, when we cannot avoid a singularity such as a band crossing point in a onedimensional system in the adiabatic deformation , can be changed across the singular region. For noncollinear and noncoplanar systems with symmetry breaking: , where is the identity element, we can also expect the dependent magnon spin since there is no unbroken generator such that triangular ().
Next, we apply a similar argument to the Hamiltonians with the rotational symmetry around the axis such as models. Since and , Eq. (8) holds only for the component spin. When there is no massless mode, and in Eq. (8) are zero. Thus, and are eigenstates of , and we cannot expect the dependent magnon spin. It is important to note that the above argument does not hold in the presence of the singularity discussed above and symmetry breaking: in which there is one massless NG mode watanabe (), or equivalently, the states are no longer the eigenstates of . We construct examples for both cases in the following parts.
Trivial example in 1D antiferromagnet.—To gain some insight into the dependent magnon spin, we first consider a simple model of a onedimensional antiferromagnet, which has been studied in the context of the spin wave field effect transistor rcheng (),
(11) 
where is the nearestneighbor exchange coupling, is the strength of the DzyaloshinskiiMoriya (DM) interaction, and is the easyaxis anisotropy. Although the DM interaction and the anisotropy breaks the symmetry, they preserve the symmetry around the axis.
For sufficiently small , the classical ground state is the Néel state with two sublattices, with up spin and with down spin [Fig. 2(a)]. Using the HolsteinPrimakoff transformation around the Néel state,
(12) 
we can rewrite Eq. (11) in terms of magnons as
(13) 
where is the onedimensional momentum, , , and . We set the lattice constant . By using a standard Bogoliubov transformation technique colpa (), we can find paraunitary matrices and that diagonalize Eq. (13). By performing numerical calculations, we plot magnon energies and the component magnon spin defined by Eq. (7) for in Figs. 2 (b) and (c). The band structure has two splitted bands with two crossing points and a finite energy gap. As shown in the previous section, the absence of the symmetry breaking ensures that is an eigenstate of . Thus, the component magnon spin is quantized, while the contribution from each sublattice does not have to be. The fact that each contribution can be over comes from the quantum nature of the antiferromagnetic magnon quantumcorrection (). In the upper and lower bands, the sign of is changed across the band crossing points, which does not conflict with our discussion above. This can be interpreted as a example of the spinmomentum locking with a collinear spin structure in momentum space.
To explore noncollinear spin structures in momentum space such as in the topological insulator surface state, we should consider classical ground states with noncollinear spin structures in real space [see Eq. (7)] noncoplanar ().
Magnon spin texture in momentum space.—As an example of a noncollinear structure, we consider the 120 structure in kagome lattice antiferromagnets, which have not only magnetic but also chemical sublattices. We analyze the following Hamiltonian:
(14) 
where are the nearestneighbor exchange couplings, and is the DM vector defined in Fig. 1. We consider the two interesting limits: (a) , and (b) , , both of which have classical ground states with the 120 structure groundstate () and preserve the symmetry around the axis. We here choose the ground state shown in Fig. 1(b). By mapping Eq. (14) to the magnon Hamiltonian and performing the numerical Bogoliubov transformation maestro (), we plot the magnon band dispersions and dependent magnon spin for each case in Fig. 3 skkim ().
The band structure for the case (a) has the finiteenergy flat band, which is reminiscence of the zero energy flat band in the classical spin liquid phase of the isotropic Heisenberg model. There is one massless NG mode associated with the symmetry breaking , and we can observe noncollinear spin structure in momentum space, as shown in Fig. 3(a). The norm of magnon spin is no longer quantized due to the absence of any spin rotational symmetries. The most striking feature is that spin structures in the highest and flat bands have the winding number defined in a closed curve around a point, while the original spinmomentum locking in electron systems is characterized by . Although the vector plot can depend on the choice of the ground state, all plots for degenerated ground states can be identified up to overall rotation in spin space, which preserves the winding number of the vortexlike spin structures. This model is thought to be realized in KFe(OH)(SO) grohol (); matan () except for some terms that slightly modify canted (). A similar model has also been investigated in terms of a topological thermal Hall effect owerre2 ().
The band structure for the case (b) also has the finiteenergy flat band and one massless NG mode for the same reasons. In addition, there are two Dirac points with a finite energy in the and points, as shown in Fig. 3(b). For each Dirac cone, a noncollinear spin structure characterized by is realized, which is a magnonic analogue of the topological insulator surface state. Note that the magnon spinmomentum locking does not require the relativistic effect, the DM interaction, while the conventional one has been found only in systems with strong relativistic effect, the spinorbit interaction. In magnonic systems, interesting physics can occur even in the absence of the DM interaction. For instance, Owerre showed that the topological thermal Hall effect occurs in such a situation owerre2 ().
Before ending this section, we remark on the relation between a noncollinear spin texture in the twodimensional Brillouin zone and a mathematical theorem. Although spin is not quantized, there is a mathematical way to characterize such spin structures, i.e., the PoincaréHopf index theorem poincare (), which states that
(15) 
where is a compact differentiable manifold, is the Euler characteristic of , which is a typical topological invariant of the manifold poincare (); euler (), are the winding numbers around isolated zero points of a vector field, and the sum of indices is over all isolated zero points. Using this theorem for the vector field on the twodimensional Brillouin zone (the twodimensional torus ), we obtain the sum rule for spinmomentum locking,
(16) 
where we use . In the spinmomentumlocked band with discussed above, the isolated zero points are , , with , respectively.
.—We here briefly discuss the detection of magnon’s spin angular momentum. A current of magnon with finite spin, known as magnon spin current, can be detected by the spin pumping sandweg (); kovalev1 () and the spin Seebeck caloritronics (); caloritronics2 (); kovalev2 () measurements. In conventional antiferromagnets without an external magnetic field, net spin current in the bulk vanishes due to the degeneracy between the up and down bands ohnuma (). The antiferromagnetic examples in Fig. 3, on the other hand, have no band degeneracy except for the crossing regions, and we can expect finite spin Seebeck signals. However, there is no established experimental method to detect directly the momentumdependent magnon spin, while the spin and angleresolved photoemission spectroscopy hasan (); xlq (); hoesch () enables us to detect the momentumdependent electron spin. To this end, we theoretically propose a setup to detect the spinmomentum locking with [Fig. 1(b)]. In spin structure, the spin direction of the magnon with is opposite to that with . This property would be observed as the spin flip under the magnon propagation in the specular reflection setup. The polarized inelastic neutron scattering chatterji () is another possibility. By investigating a change in neutron spin before and after the scattering, it is possible, in principle, to detect magnon spinmomentum locking.
In summary, we presented a theory of the magnon spinmomentum locking. We gave conditions for magnon spin to be independent of momentum and constructed examples of spinmomentum locking by avoiding such conditions. We find the first example of spinmomentum locking induced by spontaneous symmetry breaking.
.—We acknowledge many fruitful discussions with Masao Ogata, Masatoshi Imada, Hiroyasu Matsuura, Yohei Yamaji, Yusuke Kousaka, Tomonari Mizoguchi, Yohei Ema, and Yuta Kikuchi. N. O. is supported by the Japan Society for the Promotion of Science (JSPS) through Program for Leading Graduate Schools (MERIT). N. O. is also supported by JSPS KAKENHI (Grant No. 16J07110).
Footnotes
 preprint: APS/123QED
References
 V. V. Kruglyak, S. O. Demokritov, and D. Grundler, J. Phys. D , 264001 (2010).
 A. V. Chumak, V. I. Vasyuchka, A. A. Serga, and B. Hillebrands, Nature Phys. , 453 (2015).
 M. Z. Hasan and C. L. Kane, Rev. Mod. Phys. , 3045 (2010).
 X.L. Qi and S.C. Zhang, Rev. Mod. Phys. , 1057 (2011).
 D. Xiao, M.C. Chang, and Q. Niu, Rev. Mod. Phys. , 1959 (2010).
 S. Fujimoto, Phys.Rev.Lett. , 047203 (2009).
 H. Katsura, N. Nagaosa, and P. A. Lee, Phys.Rev.Lett. , 066403 (2010).
 Y. Onose, T. Ideue, H. Katsura, Y. Shiomi, N. Nagaosa, and Y. Tokura, Science , 297 (2010).
 S. A. Owerre, J. Phys.: Condens. Matter , 03LT01 (2017).
 R. Shindou, R. Matsumoto, S. Murakami, and J. I. Ohe, Phys. Rev. B , 174427 (2013).
 L. Zhang, J. Ren, J. S. Wang, and B. Li, Phys. Rev. B , 144101 (2013).
 R. Chisnell, J. S. Helton, D. E. Freedman, D. K. Singh, R. I. Bewley, D. G. Nocera, and Y. S. Lee, Phys. Rev. Lett. , 147201 (2015).
 S. A. Owerre, Journal of Applied Physics , 043903 (2016).
 F. Y. Li, Y. D. Li, Y. B. Kim, L. Balents, Y. Yu, and G. Chen, Nature Commun. , 12691 (2016).
 J. Fransson, A. M. BlackSchaffer, and A. V. Balatsky, Phys. Rev. B , 075401 (2016).
 C. Kittel, (John Wiley & Sons, New York, 1986).
 R. Cheng, S. Okamoto, and D. Xiao, Phys. Rev. Lett. , 217202 (2016).
 V. A. Zyuzin and A. A. Kovalev, Phys. Rev. Lett. , 217203 (2016).
 Y. Ohnuma, H. Adachi, E. Saitoh, and S. Maekawa, Phys. Rev. B , 014423 (2013).
 R. Cheng, M. W. Daniels, J.G. Zhu, and D. Xiao, Sci. Rep. , 24223 (2016).
 Although it is not a mathematically rigorous definition, we would like to give the intuitive meaning of the winding number below. Let be a vector field with an isolated zero . The winding number is the total number of counterclockwise turns experienced by after completing the counterclockwise motion along the closed curve around . Positive winding number means the counterclockwise turns of and conversely, negative winding number means the clockwise turns. The winding number can be calculated by the formula: , where , and is the inplane projection of .
 J. H. P. Colpa, Physica A , 327 (1978).
 In this paper, we treat the magnets in the semiclassical picture (). It would be an interesting future work to include interaction effects in systems such as nonperturbative damping discussed in A. L. Chernyshev and P. A. Maksimov, Phys. Rev. Lett. ,187203 (2016).
 In SU(2) symmetric systems, we can not define the spinmomentum locking due to the spin degeneracy. In U(1) symmetric systems, all Bloch states in the Brillouin zone have the common quantized axis, and spin structures can not be interesting.
 See the Supplemental Material for the details of calculation.
 Note that spins on different magnetic sublattices in spin Hamiltonians without chemical sublattices such as a triangular lattice antiferromagnet contribute to with the equal weight, and cancel out each other.
 In this simple symmetry breaking, the number of massless NG bosons is given by dim U(1)/=1. The formula to calculate is given in H. Watanabe and H. Murayama, Phys. Rev. Lett. , 251602 (2012).
 In usual context, this is known as the quantum correction to the classical antiferromagnetic ground state. See P. Fazekas, (World Scientific, Singapore,1999).
 Although we focus on the noncollinear spin structure, there is no reason to forbid the noncoplanar spin structure in momentum space, which could be realized in noncoplanar magnets.
 Strictly speaking, there is another ground state, structure, for the case (b). In the presence of small but finite DM interaction, which exists in realistic materials, we do not have to consider this degeneracy.
 A. D. Maestro and M. Gingras, J. Phys. Cond. Matt. , 3399 (2004).
 Note that the concept of magnon spin is different from the hvector, which is closely related to the topological property, plotted such as in S. K. Kim, H. Ochoa, R. Zarzuela, and Y. Tserkovnyak, Phys. Rev. Lett. , 227201 (2016).
 D. Grohol, K. Matan, J.H Cho, S.H. Lee, J. W. Lynn, D. G. Nocera and Y. S. Lee, Nat. Mater. , 323 (2005).
 K. Matan, D. Grohol, D. G. Nocera, T. Yildirim, A. B. Harris, S. H. Lee, S. E. Nagler, and Y. S. Lee, Phys. Rev. Lett. , 247201 (2006).
 As mentioned in Refs. grohol () and matan (), the ground state spins are canted slightly out of the kagome plane. Thus, the magnons have small zcomponent spin angular momentum accordingly.
 S. A. Owerre, Phys. Rev. B , 014422 (2017).
 V. Guillemin and A. Pollack, (American Mathematical Society, 2010).
 M. Nakahara, (A. Hilger, 1990).
 C. W. Sandweg, Y. Kajiwara, A. V. Chumak, A. A. Serga, V. I. Vasyuchka, M. B. Jungfleisch, E. Saitoh, and B. Hillebrands, Phys. Rev. Lett. , 216601 (2011).
 A. A. Kovalev, V. A. Zyuzin, and B. Li, Phys. Rev. B ,165106 (2017).
 A. A. Kovalev and V. A. Zyuzin, Phys. Rev. B , 161106 (R) (2016).
 G. E. W. Bauer, A. H. MacDonald, and S. Maekawa (Eds), , (Elsevier, 2010).
 G. E. W. Bauer, E. Saitoh, and B. J. van Wees, Nature materials , 391 (2012).
 M. Hoesch, M .Muntwiler, V. N. Petrov, M. Hengsberger, L. Patthey, M. Shi, M. Falub, T. Greber, and J. Osterwalder, Phys. Rev. B , 241401(R) (2004).
 T. Chatterji, (Elsevier, Amsterdam, 2006).