Paper The following article is Open access

Topological Weyl semimetals in the chiral antiferromagnetic materials Mn3Ge and Mn3Sn

, , , , and

Published 20 January 2017 © 2017 IOP Publishing Ltd and Deutsche Physikalische Gesellschaft
, , Focus on Topological Semimetals Citation Hao Yang et al 2017 New J. Phys. 19 015008 DOI 10.1088/1367-2630/aa5487

1367-2630/19/1/015008

Abstract

Recent experiments revealed that Mn3Sn and Mn3Ge exhibit a strong anomalous Hall effect at room temperature, provoking us to explore their electronic structures for topological properties. By ab initio band structure calculations, we have observed the existence of multiple Weyl points in the bulk and corresponding Fermi arcs on the surface, predicting antiferromagnetic Weyl semimetals in Mn3Ge and Mn3Sn. Here the chiral antiferromagnetism in the Kagome-type lattice structure is essential to determine the positions and numbers of Weyl points. Our work further reveals a new guiding principle to search for magnetic Weyl semimetals among materials that exhibit a strong anomalous Hall effect.

Export citation and abstract BibTeX RIS

Original content from this work may be used under the terms of the Creative Commons Attribution 3.0 licence. Any further distribution of this work must maintain attribution to the author(s) and the title of the work, journal citation and DOI.

1. Introduction

Recent discovery of Weyl semimetals (WSMs) [13] in realistic materials has stimulated tremendous research interest in topological semimetals, such as WSMs, Dirac semimetals, and nodal line semimetals [49], as a new frontier of condensed matter physics after the discovery of topological insulators [10, 11]. The WSMs are of particular interest not only because of their exotic Fermi-arc-type surface states but also because of their appealing bulk chiral magneto-transport properties, such as the chiral anomaly effect [1214], nonlocal transport [15, 16], large magnetoresistance, and high mobility [17]. Currently discovered WSM materials can be classified into two groups. One group breaks crystal inversion symmetry but preserves time-reversal symmetry (e.g., TaAs-family transition-metal pnictides [18, 19] and WTe2- and MoTe2-family transition-metal dichalcogenides [2026]). The other group breaks time-reversal symmetry in ferromagnets with possible tilted moments (e.g., magnetic Heusler GdPtBi [27, 28] and YbMnBi2 [29]). An antiferromagnetic (AFM) WSM compound has yet to be found, although Y2Ir2O7 with a noncoplanar AFM structure was theoretically predicted to be a WSM candidate [5].

In a WSM, the conduction and valence bands cross each other linearly through nodes called Weyl points. Between a pair of Weyl points with opposite chiralities (sink or source of the Berry curvature) [4], the emerging Berry flux can lead to the anomalous Hall effect (AHE) [30], as observed in GdPtBi [27, 28], and an intrinsic spin Hall effect (SHE), as predicted in TaAs-type materials [31], for systems without and with time-reversal symmetry, respectively. Herein, we raise a simple recipe to search for WSM candidates among materials that host strong AHE or SHE.

Recently, Mn3X (where ${X}=\mathrm{Sn}$, Ge, and Ir), which exhibit noncollinear antiferromagetic (AFM) phases at room temperature, have been found to show large AHE [3235] and SHE [36], provoking our interest to investigate their band structures for possible WSMs. In this work, we report the existence of Weyl fermions for Mn3Ge and Mn3Sn compounds and the resultant Fermi arcs on the surface by ab initio calculations, awaiting experimental verifications. Dozens of Weyl points exist near the Fermi energy in their band structure, and these can be well understood with the assistance of lattice symmetry.

2. Methods

The electronic ground states of Mn3Ge and Mn3Sn were calculated by using density-functional theory (DFT) within the Perdew–Burke–Ernzerhof-type generalized-gradient approximation (GGA) [37] using the Vienna ab initio simulation package (vasp) [38]. The $3{{\rm{d}}}^{6}4{{\rm{s}}}^{1}$, $4{{\rm{s}}}^{2}4{{\rm{p}}}^{2}$, and $5{{\rm{s}}}^{2}5{{\rm{p}}}^{2}$ electrons were considered as valance electrons for Mn, Ge, and Sn atoms, respectively. The primitive cell with experimental crystal parameters $a=b=5.352$ and c = 4.312 Å for Mn3Ge and $a=b=5.67$ and c = 4.53 Å for Mn3Sn were adopted. Spin-orbit coupling (SOC) was included in all calculations.

To identify the Weyl points with the monopole feature, we calculated the Berry curvature distribution in momentum space. The Berry curvature was calculated based on a tight-binding Hamiltonian based on localized Wannier functions [39] projected from the DFT Bloch wave functions. Chosen were atomic-orbital-like Wannier functions, which include Mn-spd and Ge-sp/Sn-p orbitals, so that the tight-binding Hamiltonian is consistent with the symmetry of ab initio calculations. From such a Hamiltonian, the Berry curvature can be calculated using the Kubo-formula approach [40]

Equation (1)

where ${{\rm{\Omega }}}_{n}^{\gamma }(\vec{k})$ is the Berry curvature in momentum space for a given band n, ${\hat{v}}_{\alpha (\beta ,\gamma )}=\tfrac{1}{{\hslash }}\tfrac{\partial \hat{H}}{\partial {k}_{\alpha (\beta ,\gamma )}}$ is the velocity operator with $\alpha ,\beta ,\gamma =x,y,z$, and $| {u}_{n}(\vec{k})\rangle $ and ${E}_{n}(\vec{k})$ are the eigenvector and eigenvalue of the Hamiltonian $\hat{H}(\vec{k})$, respectively. The summation of ${{\rm{\Omega }}}_{n}^{\gamma }(\vec{k})$ over all valence bands gives the Berry curvature vector ${\boldsymbol{\Omega }}\,({{\rm{\Omega }}}^{x},{{\rm{\Omega }}}^{y},{{\rm{\Omega }}}^{z})$.

In addition, the surface states that demonstrate the Fermi arcs were calculated on a semi-infinite surface, where the momentum-resolved local density of states (LDOS) on the surface layer was evaluated based on the Green's function method. We note that the current surface band structure corresponds to the bottom surface of a half-infinite system.

3. Results and discussion

3.1. Symmetry analysis of the AFM structure

Mn3Ge and Mn3Sn share the same layered hexagonal lattice (space group $P{6}_{3}/{mmc}$, No. 194). Inside a layer, Mn atoms form a Kagome-type lattice with mixed triangles and hexagons and Ge/Sn atoms are located at the centers of these hexagons. Each Mn atom carries a magnetic moment of 3.2 μB in Mn3Sn and 2.7 μB in Mn3Ge. As revealed in a previous study [41], the ground magnetic state is a noncollinear AFM state, where Mn moments align inside the ab plane and form 120° angles with neighboring moment vectors, as shown in figure 1(b). Along the c axis, stacking two layers leads to the primitive unit cell. Given the magnetic lattice, these two layers can be transformed into each other by inversion symmetry or with a mirror reflection (My) adding a half-lattice ($c/2$) translation, i.e., a nonsymmorphic symmetry $\{{M}_{y}| \tau =c/2\}$. In addition, two other mirror reflections (Mx and Mz) adding time reversal (T), ${M}_{x}T$ and ${M}_{z}T$, exist.

Figure 1.

Figure 1. Crystal and magnetic structures of Mn3X (where ${X}=\mathrm{Sn}$ or Ge) and related symmetry. (a) Crystal structure of Mn3X. Three mirror planes are shown in purple, corresponding to {${M}_{y}| \tau =c/2$}, MxT, and MzT symmetries. (b) Top view along the c axis of the Mn sublattice. Chiral AFM with an angle of 120° between neighboring magnetic moments is formed in each Mn layer. The mirror planes that correspond to MxT and {${M}_{y}| \tau =c/2$} are marked by dashed lines. (c) Symmetry in momentum space, My, MxT, and MzT. If a Weyl point appears at $({k}_{x},{k}_{y},{k}_{z})$, eight Weyl points in total can be generated at $(\pm {k}_{x},\pm {k}_{y},\pm {k}_{z})$ by the above three symmetry operations. For convenience, we choose the ${k}_{y}=\pi $ plane for My here.

Standard image High-resolution image

In momentum space, we can utilize three important symmetries, ${M}_{x}T$, ${M}_{z}T$, and My, to understand the electronic structure and locate the Weyl points. Suppose a Weyl point with chirality χ (+ or −) exists at a generic position ${\bf{k}}\,({k}_{x},{k}_{y},{k}_{z})$. Mirror reflection reverses χ while time reversal does not and both of them act on ${\bf{k}}$. Further, mirror reflection ${M}_{\gamma }$ preserves the Berry curvature ${{\rm{\Omega }}}^{\gamma }$ while time reversal reserves it. The transformation is as follows:

Equation (2)

Each of the above three operations doubles the number of Weyl points. Thus, eight nonequivalent Weyl points can be generated at $(\pm {k}_{x},+{k}_{y},\pm {k}_{z})$ with chirality χ and $(\pm {k}_{x},-{k}_{y},\pm {k}_{z})$ with chirality $-\chi $ (see figure 1(c)). We note that the ${k}_{x}=0/\pi $ or ${k}_{z}=0/\pi $ plane can host Weyl points. However, the ${k}_{y}=0/\pi $ plane cannot host Weyl points, because My simply reverses the chirality and annihilates the Weyl point with its mirror image if it exists.

In addition, the symmetry of the 120° AFM state is slightly broken in the materials, owing to the existence of a tiny net moment (∼0.003 μB per unit cell) [4143]. Such weak symmetry breaking seems to induce negligible effects in the transport measurement. However, it gives rise to a perturbation of the band structure, for example, shifting slightly the mirror image of a Weyl point from its position expected, as we will see in the surface states of Mn3Ge.

3.2. The anomalous Hall effect

The intrinsic anomalous Hall conductivity ${\sigma }_{\gamma }(\gamma =x,y,z)$ can be calculated by integrating the Berry curvature ${{\rm{\Omega }}}^{\gamma }$ over the whole Brillouin zone [40, 44]

Equation (3)

According to equation (2), ${{\rm{\Omega }}}^{x}$ and ${{\rm{\Omega }}}^{z}$ are odd with respect to ${M}_{x}T$ and ${M}_{y}T$, respectively. Thus, corresponding ${\sigma }^{x}$ and ${\sigma }^{z}$ are zero. Since ${{\rm{\Omega }}}^{y}$ is even with respect to the My mirror plane, corresponding ${\sigma }^{y}$ is nonzero. This is also consistent with the distribution of Weyl points in the k-space. As shown in figure 1(c), only '+' Weyl points appear on one side of the My plane and only '–' Weyl points locate on the other side of My plane. Then there are net Berry flux (starting from '+' to '–' Weyl points) ${{\rm{\Omega }}}^{y}$ crossing the My plane, resulting in the nonzero anomalous Hall conductivity ${\sigma }^{y}$. In contrast, an equal number of '+' and '–' Weyl points appear on each side of Mx (Mz) planes. Consequently, the net Berry flux of ${{\rm{\Omega }}}^{x}$ (${{\rm{\Omega }}}^{z}$) should be zero, giving rise to vanishing ${\sigma }^{x}$ (${\sigma }^{z}$). According to recent numerical calculations [36], ${\sigma }^{y}=330(133)\,{{\rm{S}}}^{-1}\,{\mathrm{cm}}^{-1}$ for Mn3Ge (Mn3Sn).

In the measurement of AHE, an external magnetic field is usually applied to uniform different magnetic domains. Further, for Mn3Ge and Mn3Sn, the triangular spins can be rotated inside the xy plane even by a very weak magnetic field due to the residual magnetic moment [43]. The rotation of an arbitrary angle can break the My and ${M}_{x}T$ symmetry, showing nonzero ${\sigma }^{y}$ and ${\sigma }^{x}$. However, ${\sigma }^{z}$ is still zero due to the ${M}_{z}T$ symmetry. As observed for both compounds in experiment [34, 35], ${\sigma }^{x,y}$ are indeed very large and ${\sigma }^{z}$ is negligible. The in-plane anomalous Hall conductivity is about 500 (100) ${{\rm{S}}}^{-1}\,{\mathrm{cm}}^{-1}$ for Mn3Ge (Mn3Sn) at low temperature, which are in the same order of magnitude as the calculations [36].

3.3. Weyl points in the bulk band structure

The bulk band structures are shown along high-symmetry lines in figure 2 for Mn3Ge and Mn3Sn. It is not surprising that the two materials exhibit similar band dispersions. At first glance, one can find two seemingly band degenerate points at Z and K points, which are below the Fermi energy. Because of ${M}_{z}T$ and the nonsymmorphic symmetry {${M}_{y}| \tau =c/2$}, the bands are supposed to be quadruply degenerate at the Brillouin zone boundary Z, forming a Dirac point protected by the nonsymmorphic space group [4547]. Given the slight mirror symmetry breaking by the residual net magnetic moment, this Dirac point is gapped at Z (as shown in the enlarged panel) and splits into four Weyl points, which are very close to each other in k space. A tiny gap also appears at the K point. Nearby, two additional Weyl points appear, too. Since the Weyl point separations are too small near both Z and K points, these Weyl points may generate little observable consequence in experiments such as those for studying Fermi arcs. Therefore, we will not focus on them in the following investigation.

Figure 2.

Figure 2. Bulk band structures for (a) Mn3Sn and (b) Mn3Ge along high-symmetry lines with SOC. The bands near the Z and K (indicated by red circles) are expanded to show details in (a). The Fermi energy is set to zero.

Standard image High-resolution image

Mn3Sn and Mn3Ge are actually metallic, as seen from the band structures. However, we retain the terminology of Weyl semimetal for simplicity and consistency. The valence and conduction bands cross each many times near the Fermi energy, generating multiple pairs of Weyl points. We first investigate the Sn compound. Supposing that the total valence electron number is Nv, we search for the crossing points between the ${N}_{{\rm{v}}}^{\mathrm{th}}$ and ${({N}_{{\rm{v}}}+1)}^{\mathrm{th}}$ bands.

As shown in figure 3(a), there are six pairs of Weyl points in the first Brillouin zone; these can be classified into three groups according to their positions, noted as W1, W2, and W3. These Weyl points lie in the Mz plane (with W2 points being only slightly off this plane owing to the residual-moment-induced symmetry breaking) and slightly above the Fermi energy. Therefore, there are four copies for each of them according to the symmetry analysis in equation (2). Their representative coordinates and energies are listed in table 1 and also indicated in figure 3(a). A Weyl point (e.g., W1 in figures 3(b) and (c)) acts as a source or sink of the Berry curvature ${\boldsymbol{\Omega }}$, clearly showing the monopole feature with a definite chirality.

Figure 3.

Figure 3. Surface states of Mn3Sn. (a) Distribution of Weyl points in momentum space. Black and white points represent Weyl points with − and + chirality, respectively. (b) and (c) monopole-like distribution of the Berry curvature near a W1 Weyl point. (d) Fermi surface at EF = 86 meV crossing the W1 Weyl points. The color represents the surface LDOS. Two pairs of W1 points are shown enlarged in the upper panels, where clear Fermi arcs exist. (e) Surface band structure along a line connecting a pair of W1 points with opposite chirality. (f) Surface band structure along the white horizontal line indicated in (d). Here p1 and p2 are the chiral states corresponding to the Fermi arcs.

Standard image High-resolution image

Table 1.  Positions and energies of Weyl points in first Brillouin zone for Mn3Sn. The positions (kx, ky, kz) are in units of π. Energies are relative to the Fermi energy EF. Each type of Weyl point has four copies whose coordinates can be generated from the symmetry as $(\pm {k}_{x},\pm {k}_{y},{k}_{z}=0)$.

Weyl point kx ky kz Chirality Energy (meV)
W1 −0.325 0.405 0.000 86
W2 −0.230 0.356 0.003 + 158
W3 −0.107 0.133 0.000 493

In contrast to Mn3Sn, Mn3Ge displays many more Weyl points. As shown in figure 4(a) and listed in table 2, there are nine groups of Weyl points. Here ${{\rm{W}}}_{\mathrm{1,2},\mathrm{7,9}}$ lie in the Mz plane with W9 on the ky axis, W4 appears in the Mx plane, and the others are in generic positions. Therefore, there are four copies of ${{\rm{W}}}_{\mathrm{1,2},\mathrm{7,4}}$, two copies of W9, and eight copies of other Weyl points. Although there are many other Weyl points in higher energies owing to different band crossings, we mainly focus on the current Weyl points that are close to the Fermi energy. The monopole-like distribution of the Berry curvature near these Weyl points is verified; see W1 in figure 4 as an example. Without including SOC, we observed a nodal-ring-like band crossing in the band structures of both Mn3Sn and Mn3Ge. SOC gaps the nodal rings but leaves isolating band-touching points, i.e., Weyl points. Since Mn3Sn exhibits stronger SOC than Mn3Ge, many Weyl points with opposite chirality may annihilate each other by being pushed by the strong SOC in Mn3Sn. This might be why Mn3Sn exhibits fewer Weyl points than Mn3Ge.

Figure 4.

Figure 4. Surface states of Mn3Ge. (a) Distribution of Weyl points in momentum space. Black and white points represent Weyl points with − and + chirality, respectively. Larger points indicate two Weyl points ($\pm {k}_{z}$) projected into this plane. (b) and (c) monopole-like distribution of the Berry curvature near a W1 Weyl point. (d) Fermi surface at EF = 55 meV crossing the W1 Weyl points. The color represents the surface LDOS. Two pairs of W1 points are shown enlarged in the upper panels, where clear Fermi arcs exist. (e) Surface band structure along a line connecting a pair of W1 points with opposite chirality. (f) Surface band structure along the white horizontal line indicated in (d). Here p1 and p2 are the chiral states corresponding to the Fermi arcs.

Standard image High-resolution image

Table 2.  Positions and energies of Weyl points in the first Brillouin zone for Mn3Ge. The positions (kx, ky, kz) are in units of π. Energies are relative to the Fermi energy EF. Each of ${{\rm{W}}}_{\mathrm{1,2,7}}$ has four copies whose coordinates can be generated from the symmetry as $(\pm {k}_{x},\pm {k}_{y},{k}_{z}=0)$. W4 has four copies at $({k}_{x}\approx 0,\pm {k}_{y},\pm {k}_{z})$ and W9 has two copies at $({k}_{x}\approx 0,\pm {k}_{y},{k}_{z}=0)$. Each of the other Weyl points has four copies whose coordinates can be generated from the symmetry as $(\pm {k}_{x},\pm {k}_{y},\pm {k}_{z})$.

Weyl point kx ky kz Chirality Energy (meV)
W1 −0.333 0.388 −0.000 57
W2 0.255 0.378 −0.000 + 111
W3 −0.101 0.405 0.097 48
W4 −0.004 0.419 0.131 + 8
W5 −0.048 0.306 0.164 + 77
W6 0.002 0.314 0.171 59
W7 −0.081 0.109 0.000 + 479
W8 0.069 −0.128 0.117 + 330
W9 0.004 −0.149 −0.000 + 470

3.4. Fermi arcs on the surface

The existence of Fermi arcs on the surface is one of the most significant consequences of Weyl points inside the three-dimensional (3D) bulk. We first investigate the surface states of Mn3Sn that have a simple bulk band structure with fewer Weyl points. When projecting ${{\rm{W}}}_{\mathrm{2,3}}$ Weyl points to the (001) surface, they overlap with other bulk bands that overwhelm the surface states. Luckily, W1 Weyl points are visible on the Fermi surface. When the Fermi energy crosses them, W1 Weyl points appear as the touching points of neighboring hole and electron pockets. Therefore, they are typical type-II Weyl points [20, 48]. Indeed, their energy dispersions demonstrate strongly tilted Weyl cones.

The Fermi surface of the surface band structure is shown in figure 3(d) for the Sn compound. In each corner of the surface Brillouin zone, a pair of W1 Weyl points exists with opposite chirality. Connecting such a pair of Weyl points, a long Fermi arc appears in both the Fermi surface (figure 3(d) and the band structure (figure 3(e)). Although the projection of bulk bands exhibit pseudo-symmetry of a hexagonal lattice, the surface Fermi arcs do not. It is clear that the Fermi arcs originating from two neighboring Weyl pairs (see figure 3(d)) do not exhibit Mx reflection, because the chirality of Weyl points apparently violates Mx symmetry. For a generic kxkz plane between each pair of W1 Weyl points, the net Berry flux points in the $-{k}_{y}$ direction. As a consequence, the Fermi velocities of both Fermi arcs point in the $+{k}_{x}$ direction on the bottom surface (see figure 3(f)). These two right movers coincide with the nonzero net Berry flux, i.e., Chern number = 2.

For Mn3Ge, we also focus on the W1-type Weyl points at the corners of the hexagonal Brillouin zone. In contrast to Mn3Sn, Mn3Ge exhibits a more complicated Fermi surface. Fermi arcs exist to connect a pair of W1-type Weyl points with opposite chirality, but they are divided into three pieces as shown in figure 4(d). In the band structures (see figures 4(e) and (f)), these three pieces are indeed connected together as a single surface state. Crossing a line between two pairs of W1 points, one can find two right movers in the band structure, which are indicated as p1 and p2 in figure 4(f). The existence of two chiral surface bands is consistent with a nontrivial Chern number between these two pairs of Weyl points.

4. Summary

In summary, we have discovered the Weyl semimetal state in the chiral AFM compounds Mn3Sn and Mn3Ge by ab initio band structure calculations. Multiple Weyl points were observed in the bulk band structures, most of which are type II. The positions and chirality of Weyl points are in accordance with the symmetry of the magnetic lattice. For both compounds, Fermi arcs were found on the surface, each of which connects a pair of Weyl points with opposite chirality, calling for further experimental investigations such as angle-resolved photoemission spectroscopy. The discovery of Weyl points verifies the large anomalous Hall conductivity observed recently in titled compounds. Our work further reveals a guiding principle to search for Weyl semimetals among materials that exhibit a strong anomalous Hall effect.

Acknowledgments

We thank Claudia Felser, Jürgen Kübler and Ajaya K Nayak for helpful discussions. We acknowledge the Max Planck Computing and Data Facility (MPCDF) and Shanghai Supercomputer Center for computational resources and the German Research Foundation (DFG) SFB-1143 for financial support.

Please wait… references are loading.