# The family of topological Hall effects for electrons in skyrmion crystals

###### Abstract

Hall effects of electrons can be produced by an external magnetic field, spin orbit-coupling or a topologically non-trivial spin texture. The topological Hall effect (THE) — caused by the latter — is commonly observed in magnetic skyrmion crystals. Here, we show analogies of the THE to the conventional Hall effect (HE), the anomalous Hall effect (AHE), and the spin Hall effect (SHE). In the limit of strong coupling between conduction electron spins and the local magnetic texture the THE can be described by means of a fictitious, ‘emergent’ magnetic field. In this sense the THE can be mapped onto the HE caused by an external magnetic field. Due to complete alignment of electron spin and magnetic texture, the transverse charge conductivity is linked to a transverse spin conductivity. They are disconnected for weak coupling of electron spin and magnetic texture; the THE is then related to the AHE. The topological equivalent to the SHE can be found in antiferromagnetic skyrmion crystals. We substantiate our claims by calculations of the edge states for a finite sample. These states reveal in which situation the topological analogue to a quantized HE, quantized AHE, and quantized SHE can be found.

## 1 Introduction

If an electric current is applied to an electrical conductor in the presence of a static magnetic field in perpendicular direction, a voltage across the sample is observed. This ‘Hall effect’ of electrons was discovered almost 140 years ago Hall (1879). In the following and especially in the last 20 years it has been found that the external magnetic field is not mandatory to produce a Hall effect and can be exchanged for example by spin-orbit coupling Nagaosa et al. (2010) or a chiral magnetic texture Bruno et al. (2004) that produces an effective magnetic field, that is a ‘Berry curvature’ Berry (1984). The family of Hall effects for electrons in solids was extended by the quantum Hall effect Landau (1930); Hofstadter (1976); Klitzing et al. (1980), anomalous Hall effect Nagaosa et al. (2010), quantum anomalous Hall effect Liu et al. (2008); Chang et al. (2013), spin Hall effect D’yakonov and Perel (1971); Kato et al. (2004), quantum spin Hall effect Kane and Mele (2005); Bernevig et al. (2006), and topological Hall effect Neubauer et al. (2009); Schulz et al. (2012); Kanazawa et al. (2011); Lee et al. (2009); Li et al. (2013); Bruno et al. (2004); Hamamoto et al. (2015); Yin et al. (2015); Lado and Fernández-Rossier (2015); Göbel et al. (2017, 2017); Ndiaye et al. (2017).

The Hall resistivity

(1) |

is the off-diagonal element of the resistivity tensor , that connects electric field and current density via Ohm’s law . It is given by contributions from the ordinary Hall effect (HE), the anomalous Hall effect (AHE) and eventually the topological Hall effect (THE) in magnetic systems. The HE can appear quantized as quantum Hall effect (QHE) and is due to an external magnetic field that couples to the charge of the conduction electrons; the electrons are deflected by a Lorentz force in a semiclassical picture. The AHE arises due to spin-orbit coupling (SOC) with intrinsic (Berry curvature) and extrinsic (skew-scattering and side jump) contributions Nagaosa et al. (2010), all of which deflect electrons with opposite spin into opposite directions. In a ferromagnet, the imbalance of majority and minority conduction electrons leads to both spin and charge transport. In non-magnetic materials the charge accumulation at both sides of the sample is compensated, leading to a pure spin Hall effect (SHE).

This Paper addresses the THE of electrons in topologically non-trivial magnetic textures in two dimensions. The probably most prominent example is the skyrmion Skyrme (1962); Bogdanov and Yablonskii (1989); Bogdanov and Hubert (1994); Rößler et al. (2006); Mühlbauer et al. (2009). Its winding is quantified by the topological charge density

(2) |

that integrates to the integer topological charge (or skyrmion number) Nagaosa and Tokura (2013)

(3) |

The spins of the conduction electrons tend to align with the magnetic texture, thereby accumulating a Berry phase Everschor-Sitte and Sitte (2014). If skyrmions are arranged periodically in a skyrmion crystal (SkX) this Berry phase can be translated to a Berry curvature in reciprocal space that deflects traversing electrons and leads to nonzero transverse Hall coefficients for charge and spin. Skyrmion crystals have first been observed in non-centrosymmetric B20 bulk-compounds like MnSi Mühlbauer et al. (2009) and FeGe Yu et al. (2011). They also appear at interfaces, like in Fe/Ir Heinze et al. (2011).

This Paper is dedicated to characterize precisely the THE and to update the classification of electronic Hall effects. The THE is commonly mentioned in line with the other Hall effects. We think it is more sensible to speak of topological versions of the conventional, the anomalous and the spin Hall effect because the THE manifests itself differently in dependence on the strength of the Zeeman interaction of electron spin and magnetic texture. We suggest a nomenclature (Fig. 1) for the variety of THEs and tell which effect can be expected in the different scenarios.

In the strong-coupling limit, the Hamiltonians of THE and a spin-polarized QHE can be mapped onto each other (Sect. 4). This locally fully spin-polarized quantized topological Hall effect (QTHE) is the extreme case of a quantized topological anomalous Hall effect (QTAHE). In the weak-coupling limit spin and charge transport coefficients are decoupled, what corresponds to the typical anomalous Hall scenario in its topological version (TAHE; Sect. 5). The topological analogue to the spin Hall effect Göbel et al. (2017); Buhl et al. (2017); Akosa et al. (2017) (TSHE), i. e., pure spin currents, is found in antiferromagnetic skyrmion crystals and completes the topological Hall trio (Sect. 6; see Fig. 1). By calculating edge states of a finite sample we illustrate when these topological Hall effects show up in their quantized version.

## 2 Model and methods

In this Paper we focus on the topological contribution to the Hall effect that arises due to the nonzero topological charge density of the magnetic texture formed by localized magnetic moments . Therefore, we consider the tight-binding Hamiltonian

(4) |

that features nearest-neighbor hopping of conduction electrons ( creation and annihilation operators at site , amplitude ), and Hund’s coupling (quantified by ), i. e., the Zeeman interaction of the electrons’ spins with the fixed magnetic texture [visualized in Fig. 1(e)]. The magnetic texture is thought to be produced by localized electrons that are not explicitly featured in the calculation.

To segregate the effect of the magnetic texture, we neglect (i) SOC which would contribute to the AHE and (ii) the influence of external magnetic fields that are commonly needed to stabilize a SkX phase; they would contribute to the conventional or quantized Hall effect.

Having diagonalized the Hamiltonian, the eigenenergies and the eigenvectors , both depending on the wave vector , allow to calculate the Berry curvature Berry (1984)

(5) |

whose integral over the Brillouin zone (BZ) yields the Chern number

(6) |

of band and the topological Hall conductivity Nagaosa et al. (2010)

(7) |

Concerning spin transport in a ferromagnet, one velocity has to be replaced by the spin current operator , with that accounts for spin-up and spin-down orientation with respect to the global quantization axis Gradhand et al. (2012). The noncollinear texture of a skyrmion requires to consider the local quantization axis (at each site) by means of the operator . We obtain the spin Berry curvature

(8) |

and the spin Hall conductivity for spin polarization parallel to the local quantization axes

(9) |

## 3 Influence of the Hund’s coupling on the band structure

Before presenting results for the transverse charge and spin conductivities, the effect of the Hund’s coupling on the electron spin and on the band structure is prerequisite. The strength of the interaction is quantified by .

Aiming at the simplest, yet most distinct manifestation of topological properties of the SkX in band structure and transport coefficients, we consider only s-electrons on a square structural lattice (lattice constant ). The respective zero-field band structure (band structure of the structural lattice for ) has one band,

(10) |

with a maximum at , a minimum at , and two saddle points at (van Hove singularity, VHS). On top of this, it is spin-degenerate.

The square magnetic unit cell of the SkX, with a period of , leads to back-folding of the bands into the magnetic Brillouin zone and thus to bands. A nonzero coupling breaks time reversal symmetry, and the spin degeneracy is lifted. This is accompanied by a shift of the bands to lower (spin antiparallel to the magnetic texture) and higher (spin parallel to the magnetic texture) energies. For (half the band width of ), the band structure is separated into two blocks, the electron spin is almost completely aligned with the magnetic texture, and small band gaps show up. These findings are best understood in the strong-coupling limit .

## 4 Quantized topological Hall effect in the strong-coupling limit

### 4.1 Transformation to the emergent field

In the limit the electron spin completely aligns with the texture to minimize the system’s energy. In this case the Zeeman term of the Hamiltonian can be diagonalized by applying a set of unitary transformations,

(11) |

so that the spin at site points into the direction of the local coordinate system. The new electron operators read , while the new hopping strength is now a non-diagonal (22) matrix. The full transformation is discussed in Ref. Hamamoto et al., 2015.

In the strong-coupling limit, spin-parallel electrons cannot align antiparallel to the magnetic texture and vice versa. For the parallel aligned electrons, only the element

(12) |

of the hopping is relevant, which depends on the azimuthal and the polar angle of the local magnetic moment at site . It can be transformed into an effective hopping strength

(13) |

with the angle between the moments at sites and ; the phase is given by

(14) |

With these substitutions the Hamiltonian becomes

(15) |

which describes fully spin-aligned electrons (equivalent to spinless electrons; creation and annihilation ). The constant Zeeman term has been dropped. The cosine in (13) merely scales the hopping strength and converges to for large skyrmions (); the THE has then been mapped onto the QHE described by the Hamiltonian

(16) |

In the quantum Hall effect spinless electrons interact with a uniform magnetic field , that is oriented out of plane (spin-polarized electrons can be described likewise). This field couples to the charge, rather than to the spin of the conduction electrons and affects their kinetic energy. In the tight-binding description this leads to complex hopping amplitudes with phase factors

(17) |

whose phase is determined by the vector potential , with . The above integral is along the hopping path.

The phase for the topological Hall effect can be identified with the phase from the quantum Hall effect. The local skyrmion density can therefore be identified with a fictitious magnetic field that is collinear but inhomogeneous.

Application of the unitary transformations is exact; however, taking only a diagonal element is strictly valid only in the strong-coupling limit . In this limit the Hund’s coupling can be expressed as an interaction of the skyrmion’s emergent field with the charge of the conduction electrons. Considering the element instead of the element, one would have described the lower block of bands, in which the electron spin is aligned antiparallel to the magnetic texture. In that case, the emergent field as well as the topological Hall conductivity change sign.

Describing the QHE on a lattice restricts the magnetic field to discrete values, so that the magnetic flux per magnetic unit cell is an integer multiple of the flux quantum . Therefore, a large magnetic unit cell has to be considered to account for the appropriate phase . The magnetic flux per plaquette then reads ; and are coprime integers Hofstadter (1976). The emergent field of the skyrmion induces an integer number of flux quanta per unit cell as well. Both systems differ by the fact that the emergent field is not homogeneous; still one can find a quantum Hall system with a magnetic field equal to the average emergent field of the SkX. The average flux per plaquette then reads .

This equivalence of THE and QHE in the limit allows to carry over all the known results from one effect to the other. Bands in the THE can, for example, be interpreted as dispersive Landau levels which carry a Chern number of [] for the upper [lower] block. One result is a quantized transverse conductivity, provided the Fermi level lies within a band gap. Furthermore, behaves peculiarly for Fermi energies near VHSs, as is discussed in the next section.

### 4.2 Conductivity

The analogy of THE and QHE holds strictly speaking in the limit . Nevertheless, the emergent-field description is reasonable as long as block separation in the band structure is given (). We emphasize that the computations were performed with the exact Hamiltonian (4); the numerical results are only interpreted by means of the emergent field.

The transverse conductivity [red in Fig. 2(a)] depends strongly on the position of the Fermi energy. The two blocks in the band structure, easily recognized in the spectrum, show opposite sign, which is explained by the opposite spin alignment with respect to the magnetic texture, giving emergent fields of opposite sign. This leads to deflection of electrons into opposite transverse directions in a semi-classical picture. It is sufficient to discuss only the upper block.

Interpreting the bands as dispersive Landau levels with a Chern number of tells that the conductivity decreases stepwise upon increasing . Due to the inhomogeneity of the emergent field, the bands are ‘deformed’ and have a nonzero band width. The spectrum consists therefore of energy regions in which is quantized [insulating system; QTHE; Fig. 1(d)] or not quantized [metallic system; THE; Fig. 1(a)].

The energies of the Landau levels are determined by the constraint that, like for the QHE, for each Landau level the number of states is identical (Onsager’s quantization scheme). At these energies, the corresponding Fermi line of the zero-field band structure encloses an area of

(18) |

in the Brillouin zone; is the area of the (smaller) magnetic Brillouin zone. The curvature of these Fermi lines dictates both the fermion character and the sign of the transverse conductivity: electron-like (hole-like) orbits give a negative (positive) sign for energies below (above) . The sign of the transverse conductivity changes at the VHS Göbel et al. (2017, 2017, 2018); Hatsugai et al. (2006); Arai and Hatsugai (2009), because it separates electron- from hole-like Fermi lines. The particular Landau level closest to the VHS thus carries a large Chern number of which compensates the Chern number of all other bands (for details see Ref. Göbel et al., 2017). The zero-field band structure dictates therefore the shape of : different lattices lead to different spectra [cf. a triangular lattice in Ref. Göbel et al., 2017 and a honeycomb lattice in Ref. Göbel et al., 2017].

### 4.3 Edge states and higher order skyrmions

The bulk-boundary correspondence connects the sum of Chern numbers of the occupied bands with the number of topologically protected edge states Hatsugai (1993a, b), the latter being computed by Green function renormalization for the semi-infinite system Henk and Schattke (1993); Bödicker et al. (1994).

Figure 3(a) illustrates the edge states of a skyrmion crystal in the strong-coupling limit . Due to the large degree of spin polarization we show only the lower block. The number of edge states changes by for high and low energies. The edge states change from left to right propagation at , which complies with the sign change of the conductivity.

In Fig. 3(b) the system has been manipulated by renormalizing the hopping [ according to Eq. (13)]. Additionally, the skyrmion texture has been tilted so that the emergent field becomes almost homogeneous. This way, the topological charge of a skyrmion remains invariant but the SkX mimics a quantum Hall system with dispersionless Landau levels [cf. Figs. 3(b) and (c)] even better. The dispersion relations of edge states in the quantum Hall system and the SkX are very similar, especially their slope. Keep in mind that the SkX magnetic unit cell comprises sites, while for the QH system it comprises sites (with equal to the number of atoms in the magnetic unit cell). The X direction of the QH system is times as long as that of the SkX.

The equivalence of the spin polarized QTHE and the QHE holds also for larger skyrmion numbers. The ratio dictates the ratio for the related quantum Hall system. Figs. 3(d) and (e) show the results for the ratio , i. e., for a skyrmion with . The expansion

(19) |

tells that the band structure shows bundles that are connected by surface states like in the case (Refs. Chang and Niu, 1995, 1996). Each of these bundles consists of three subbands with a total Chern number of , except for the bundle at the VHS, which consists of four bands with a total Chern number of , compensating the other Chern numbers. This behavior is well understood for the QHE and can be carried over to the topological version. The Diophantine equation Dana et al. (1985); Thouless et al. (1982); Chang and Niu (1996) tells the Chern numbers of the subbands. Furthermore, the essential physics are represented by Hofstadter butterflies Hofstadter (1976); such figures depict the band energies versus magnetic flux, which is characterized by the ratio for the QH case and can now be understood as the magnetic flux of the emergent field given by .

### 4.4 Spin conductivity

The transverse spin conductivity of the SkX in the strong-coupling limit is shown in Fig. 2 (panel a; black curve). Charge and spin transport appear inseparably linked, which is attributed to the high degree of spin polarization. In the case of non-overlapping blocks () the results are well resembled by the emergent field picture. The spin conductivity is given by for all Fermi energies in the upper block. For the lower block the charge conductivity is sign-reversed but the spin conductivity is not. Due to the opposite sign of the emergent field, electrons with opposite spins are deflected into opposite directions.

In the strong-coupling limit, the THE becomes a spin-polarized quantized topological Hall effect, which is related to a spin-polarized QHE brought about by a magnetic field [Fig. 1(a)]. The spin-polarized QHE is also the extreme case of a QAHE. Therefore the spin-polarized THE can also be labeled TAHE [Fig. 1(a)] and the quantized version, for Fermi energies located in a band gap, can also be labeled QTAHE [Fig. 1(d)]. The Hall conductivity and the transverse spin conductivity of both compared systems agree well.

Note, that we calculate a spin conductivity in a local coordinate system. The measured spin accumulation (in a global system) is determined by this quantity but the texture at the sample’s edge has to be taken into account [cf. accumulated spins in Fig. 1(a)].

Commonly, in an AHE scenario charge and spin transport are (strictly) uncoupled. This happens in the weak-coupling limit, in which complete spin-polarization is not given for all Fermi energies.

## 5 Topological anomalous Hall effect in the weak-coupling limit

In materials with weaker Hund’s coupling , the separation into non-overlapping blocks of bands is not given, hence the emergent-field picture becomes invalid.

In collinear ferromagnets the Hall effect can be treated well in a two-channel model, where each channel deals with one kind of either spin parallel or spin antiparallel electrons. This picture breaks down if the two spin species hybridize, for example, due to spin-orbit coupling or, as in our case, due to a noncollinear magnetic texture. When the electron spin does not follow the texture adiabatically, complete local spin polarization, i. e., parallel or antiparallel alignment, is not given. Still, we utilize the two channel model and compare its prediction to the actual results.

Figures 2(b) and (c) show charge and spin conductivities for and , respectively. Both conductivities change sign at the energy of the zero-field VHS . For energies, at which the two zero-field bands do not overlap () both curves are well explained by the two-channel model, because the electron spin is almost completely aligned with the magnetic texture.

For energies at which the blocks overlap, the two-channel picture describes qualitatively the spectrum but fails to predict the amplitude, especially for small . This is due to the fact that the spin is locally not aligned parallel or antiparallel to the texture. The bands become dispersive, and the eigenstates of the Hamiltonian show contributions from both spin species.

For the charge conductivity can exceed the spin Hall conductivity in magnitude (in units of their respective quanta) if spin-antiparallel hole-like states hybridize with spin-parallel electron-like states for energies [cf. for and in Figs. 2(b) and (c)]. In the two-channel approximation, then, electrons with opposite spins are deflected into the same direction and the spin conductivity can drop to zero.

For the spin conductivity can exceed the charge conductivity. This can happen for , where spin-parallel electron-like states mix with spin-antiparallel electron-like states, and for , where spin-parallel hole-like states mix with spin-antiparallel hole-like states. In the two-channel approximation, electrons with opposite spins are deflected into opposite directions [this corresponds to the scenario depicted in Fig. 1(b)]. This increases the spin conductivity and decreases the charge conductivity compared to the case of non-overlapping zero-field bands. Due to the missing alignment the actual transverse conductivities are always reduced in magnitude compared to the results expected for the two-channel model.

The scenario of mixed zero-field blocks corresponds to the AHE in the stricter sense [cf. Fig. 1(b)], i. e., where both types of spin carriers dominate the transport effects. In this case the emergent field interpretation breaks down and bands do not behave like Landau levels. They become dispersive and band gaps disappear. Therefore, edge states are visible but are always superposed by bulk bands and the quantized version of the TAHE for mixed carriers cannot be observed.

For Fermi energies , where only one spin species is present, the local spin alignment is given and bands become merely flat. This QTAHE is again equivalent to a spin-polarized version of the QHE as in the strong-coupling limit (cf. Sect. 4).

## 6 Topological spin Hall effect in antiferromagnetic skyrmion crystals

The nontopological AHE in a ferromagnet arises due to SOC and a net magnetization. In a SkX, the effect of the SOC is taken over by the magnetic texture, which also induces the block separation. Therefore, a pure topological spin Hall effect — the number of spin-parallel electrons is identical to the number of spin-antiparallel electrons at — could only be achieved when the texture is removed. Consequently, all transverse charge and spin transport phenomena would vanish as well.

The topological analogue of the SHE Buhl et al. (2017); Göbel et al. (2017); Akosa et al. (2017) has been established in antiferromagnetic (AFM) skyrmions Zhang et al. (2016a); Barker and Tretiakov (2016); Zhang et al. (2016b) and antiferromagnetic SkXs Göbel et al. (2017). These textures consist of two sublattices or layers with conventional skyrmions, but with opposite magnetic moments in the two sublattices. Hence, the topological charge and the emergent field of each sublattice have opposite signs, leading to a zero total topological charge and a vanishing THE.

In this section we briefly review the results of Ref. Göbel et al., 2017 to complete the family of topological Hall effects. First, we consider only hopping within a sublattice. The AFM-SkX then consists of two non-interacting SkXs on the two sublattices. The conduction electrons are deflected into opposite directions depending on the sublattice on which they ‘live’ [Fig. 1(c)]. Parallel spins (arbitrarily labeled as ‘up’) in one sublattice are interpreted as antiparallel spins in the other sublattice (also ‘up’). For the calculation of the spin Hall conductivity this is respected by a modified that enters Eq. (8) via .

The results can be understood as two copies of Fig. 1(a). Locally spin-parallel electrons from sublattice A (spin up) are deflected to the left, and parallel electrons from sublattice B (spin down because of the reversed magnetization) are deflected to the right. This leads to a zero charge Hall conductivity and a nonzero spin Hall conductivity, i. e., a TSHE. looks just like for the SkX, but with twice the magnitude because of the two sublattices. Within the band gaps, the TSHE is quantized [QTSHE, Fig. 1(f)].

In real materials, nearest-neighbor hopping, i. e., inter-sublattice hopping, is allowed. Nevertheless, the THE in an AFM-SkX vanishes, in contrast to a topological spin Hall effect, as shown in Ref. Göbel et al., 2017. appears modified by the altered zero-field band structure.

## 7 Conclusion

We have analyzed the electron transport in different coupling regimes in conventional and antiferromagnetic skyrmion crystals.

In the strong-coupling limit, the topological equivalent to a fully spin-polarized QHE is found and explained in the emergent-field picture. Because the electron spin follows the texture adiabatically, charge and spin transport are twined and it can therefore also be seen as a special case of a TAHE or a QTAHE.

In the weak-coupling limit, spin-mixing terms in the Hamiltonian prevent the interpretation by means of the emergent field picture. The charge conductivity is no longer quantized and the spin conductivity is not proportional to the charge conductivity. Here, one is confronted with a topological analogue to the typical anomalous Hall effect: an unequal number of electrons with their spin aligned parallel and antiparallel to the magnetic texture propagate into opposite transverse directions, resulting in a decoupling of charge and spin conductivities.

An extreme version of the anomalous Hall effect is the spin Hall effect, with identical numbers of electrons with spin parallel and antiparallel to the magnetic texture. This effect is not possible in a skyrmion crystal, but in an antiferromagnetic skyrmion crystal. We find a TSHE or QTSHE, depending on the position of the Fermi energy.

We have shown that topologically nontrivial magnetic textures lead to a topological Hall effect, which is viewed not as a separate Hall effect but as the topological version of well-known Hall effects previously discussed in ferromagnets or conventional metals.

###### Acknowledgements.

This work is supported by SPP 1666 and SFB 762 of Deutsche Forschungsgemeinschaft (DFG).## Author contribution

B. Göbel conducted all calculations. All authors discussed the results and contributed to the manuscript.

## References

- Hall (1879) E. H. Hall, American Journal of Mathematics 2, 287 (1879).
- Nagaosa et al. (2010) N. Nagaosa, J. Sinova, S. Onoda, A. MacDonald, and N. Ong, Rev. Mod. Phys. 82, 1539 (2010).
- Bruno et al. (2004) P. Bruno, V. Dugaev, and M. Taillefumier, Phys. Rev. Lett. 93, 096806 (2004).
- Berry (1984) M. V. Berry, in Proceedings of the Royal Society of London A: Mathematical, Physical and Engineering Sciences (The Royal Society, 1984), vol. 392, pp. 45–57.
- Landau (1930) L. Landau, Z. Phys. 64, 629 (1930).
- Hofstadter (1976) D. R. Hofstadter, Phys. Rev. B 14, 2239 (1976).
- Klitzing et al. (1980) K. v. Klitzing, G. Dorda, and M. Pepper, Phys. Rev. Lett. 45, 494 (1980).
- Liu et al. (2008) C.-X. Liu, X.-L. Qi, X. Dai, Z. Fang, and S.-C. Zhang, Phys. Rev. Lett. 101, 146802 (2008).
- Chang et al. (2013) C.-Z. Chang, J. Zhang, X. Feng, J. Shen, Z. Zhang, M. Guo, K. Li, Y. Ou, P. Wei, L.-L. Wang, et al., Science 340, 167 (2013).
- D’yakonov and Perel (1971) M. D’yakonov and V. Perel, Soviet Journal of Experimental and Theoretical Physics Letters 13, 467 (1971).
- Kato et al. (2004) Y. Kato, R. Myers, A. Gossard, and D. Awschalom, Science 306, 1910 (2004).
- Kane and Mele (2005) C. L. Kane and E. J. Mele, Phys. Rev. Lett. 95, 226801 (2005).
- Bernevig et al. (2006) B. A. Bernevig, T. L. Hughes, and S.-C. Zhang, Science 314, 1757 (2006).
- Neubauer et al. (2009) A. Neubauer, C. Pfleiderer, B. Binz, A. Rosch, R. Ritz, P. Niklowitz, and P. Böni, Phys. Rev. Lett. 102, 186602 (2009).
- Schulz et al. (2012) T. Schulz, R. Ritz, A. Bauer, M. Halder, M. Wagner, C. Franz, C. Pfleiderer, K. Everschor, M. Garst, and A. Rosch, Nature Phys. 8, 301 (2012).
- Kanazawa et al. (2011) N. Kanazawa, Y. Onose, T. Arima, D. Okuyama, K. Ohoyama, S. Wakimoto, K. Kakurai, S. Ishiwata, and Y. Tokura, Phys. Rev. Lett. 106, 156603 (2011).
- Lee et al. (2009) M. Lee, W. Kang, Y. Onose, Y. Tokura, and N. Ong, Phys. Rev. Lett. 102, 186601 (2009).
- Li et al. (2013) Y. Li, N. Kanazawa, X. Yu, A. Tsukazaki, M. Kawasaki, M. Ichikawa, X. Jin, F. Kagawa, and Y. Tokura, Phys. Rev. Lett. 110, 117202 (2013).
- Hamamoto et al. (2015) K. Hamamoto, M. Ezawa, and N. Nagaosa, Phys. Rev. B 92, 115417 (2015).
- Yin et al. (2015) G. Yin, Y. Liu, Y. Barlas, J. Zang, and R. K. Lake, Phys. Rev. B 92, 024411 (2015).
- Lado and Fernández-Rossier (2015) J. L. Lado and J. Fernández-Rossier, Phys. Rev. B 92, 115433 (2015).
- Göbel et al. (2017) B. Göbel, A. Mook, J. Henk, and I. Mertig, Phys. Rev. B 95, 094413 (2017).
- Göbel et al. (2017) B. Göbel, A. Mook, J. Henk, and I. Mertig, New J. Phys. 19, 063042 (2017).
- Ndiaye et al. (2017) P. B. Ndiaye, C. A. Akosa, and A. Manchon, Phys. Rev. B 95, 064426 (2017).
- Skyrme (1962) T. H. R. Skyrme, Nuclear Physics 31, 556 (1962).
- Bogdanov and Yablonskii (1989) A. Bogdanov and D. Yablonskii, Zh. Eksp. Teor. Fiz 95, 182 (1989).
- Bogdanov and Hubert (1994) A. Bogdanov and A. Hubert, J. Magn. Magn. Mater. 138, 255 (1994).
- Rößler et al. (2006) U. Rößler, A. Bogdanov, and C. Pfleiderer, Nature 442, 797 (2006).
- Mühlbauer et al. (2009) S. Mühlbauer, B. Binz, F. Jonietz, C. Pfleiderer, A. Rosch, A. Neubauer, R. Georgii, and P. Böni, Science 323, 915 (2009).
- Nagaosa and Tokura (2013) N. Nagaosa and Y. Tokura, Nature Nanotechnology 8, 899 (2013).
- Everschor-Sitte and Sitte (2014) K. Everschor-Sitte and M. Sitte, J. Appl. Phys. 115, 172602 (2014).
- Yu et al. (2011) X. Yu, N. Kanazawa, Y. Onose, K. Kimoto, W. Zhang, S. Ishiwata, Y. Matsui, and Y. Tokura, Nature Mater. 10, 106 (2011).
- Heinze et al. (2011) S. Heinze, K. Von Bergmann, M. Menzel, J. Brede, A. Kubetzka, R. Wiesendanger, G. Bihlmayer, and S. Blügel, Nature Physics 7, 713 (2011).
- Göbel et al. (2017) B. Göbel, A. Mook, J. Henk, and I. Mertig, Phys. Rev. B 96, 060406(R) (2017).
- Buhl et al. (2017) P. M. Buhl, F. Freimuth, S. Blügel, and Y. Mokrousov, physica status solidi (RRL)-Rapid Research Letters 11, 1700007 (2017).
- Akosa et al. (2017) C. A. Akosa, O. Tretiakov, G. Tatara, and A. Manchon, arXiv preprint arXiv:1709.02931 (2017).
- Gradhand et al. (2012) M. Gradhand, D. Fedorov, F. Pientka, P. Zahn, I. Mertig, and B. Györffy, Journal of Physics: Condensed Matter 24, 213202 (2012).
- Göbel et al. (2018) B. Göbel, A. Mook, J. Henk, and I. Mertig, arXiv preprint:1802.06411 (2018).
- Hatsugai et al. (2006) Y. Hatsugai, T. Fukui, and H. Aoki, Phys. Rev. B 74, 205414 (2006).
- Arai and Hatsugai (2009) M. Arai and Y. Hatsugai, Phys. Rev. B 79, 075429 (2009).
- Hatsugai (1993a) Y. Hatsugai, Phys. Rev. B 48, 11851 (1993a).
- Hatsugai (1993b) Y. Hatsugai, Phys. Rev. Lett. 71, 3697â3700 (1993b).
- Henk and Schattke (1993) J. Henk and W. Schattke, Computer Physics Communications 77, 69 (1993).
- Bödicker et al. (1994) A. Bödicker, W. Schattke, J. Henk, and R. Feder, Journal of Physics. Condensed Matter 6, 1927 (1994).
- Chang and Niu (1995) M.-C. Chang and Q. Niu, Phys. Rev. Lett. 75, 1348 (1995).
- Chang and Niu (1996) M.-C. Chang and Q. Niu, Phys. Rev. B 53, 7010 (1996).
- Dana et al. (1985) I. Dana, Y. Avron, and J. Zak, Journal of Physics C: Solid State Physics 18, L679 (1985).
- Thouless et al. (1982) D. Thouless, M. Kohmoto, M. Nightingale, and M. Den Nijs, Phys. Rev. Lett. 49, 405 (1982).
- Zhang et al. (2016a) X. Zhang, Y. Zhou, and M. Ezawa, Scientific Reports 6, 24795 (2016a).
- Barker and Tretiakov (2016) J. Barker and O. A. Tretiakov, Phys. Rev. Lett. 116, 147203 (2016).
- Zhang et al. (2016b) X. Zhang, Y. Zhou, and M. Ezawa, Nature Communications 7, 10293 (2016b).