Spin-1 Dirac half-metal, spin-gapless semiconductor, and spin-polarized massive Dirac dispersion in transition metal dihalide monolayers
11 Spin-1 Dirac half-metal, spin-gapless semiconductor, and spin-polarized massive Dirac dispersion in transition metal dihalide monolayers
Muhammad Nadeem † , Xiaolin Wang Institute for Superconducting and Electronic Materials, Australian institute for innovative materials, and ARC Centre of Excellence in Future Low-Energy Electronics Technologies (FLEET), University of Wollongong, Australia, NSW 2500, Australia.
Spin-1 condensed matter systems characterized by the combination of a Dirac-like dispersion and flat bands are ideal for realizing high-temperature electronics and spintronic technologies in the absence of external magnetic field. In this study, we propose a three-band tight binding model, with spin-polarized Haldane-like next-nearest-neighbour tunnelling, on dice lattice and show that spin-1 Dirac half-metal, spin-1 Dirac spin-gapless semiconductor, and spin-polarized spin-1 massive Dirac dispersion with nontrivial topology can exist in two-dimensional ferromagnetic condensed matter systems with electron spin polarization P = 1. The proposed spin-polarized spin-1 phases can be realized in ferromagnetic transition metal dihalides MX monolayers effectively. By using first principle calculations, we show that a small compressive strain leads MX monolayers to be spin-one Dirac half-metal for M = Fe and X = Br, Cl while spin-one Dirac spin-gapless semiconductor for M = Co and X = Br, Cl. Spin-one Dirac spin-gapless semiconductors CoBr and CoCl embeds flat band ferromagnetism where spin-orbit coupling opens a topologically non-trivial Dirac gap between dispersing valance and conduction band while leaving flat band unaffected. The intrinsic flat-band ferromagnetism in spin-polarized spin-1 massive Dirac dispersion plays key role in materializing quantum anomalous Hall state with Chern number C = -2. I. INTRODUCTION
Discovery of exotic topological phases hosting novel low energy quasiparticles is one of the key research areas in condensed matter physics. In recent years, Weyl semimetals, Dirac semimetals, and superconducting heterostructures have been discovered where Dirac, Weyl, and Majorana fermions exist as low-energy excitations in these condensed matter systems . Moreover, unlike high energy physics where the constraint of Lorentz invariance classifies the fermions in spin S=1/2 Dirac, Weyl, and Majorana particles only, spin S=1 and S=3/2 low energy fermion-like excitations are also allowed in condensed matter systems where unconventional three, four, six and eight-fold band crossings are protected by space group symmetries . Such condensed matter systems with multi-band crossings have unique surface states and transport properties . Among these exotic three dimensional topological phases, condensed matter systems with S=1 chiral fermions characterized by the combination of a Dirac-like dispersion and flat bands are ideal for realizing high-temperature electronics and spintronic technologies in the absence of external magnetic field. This unique feature of S=1 condensed matter systems is indebted to the presence of following two key features: (i) topologically trivial flat bands and (ii) the nontrivial overall bulk band topology. These features have also been studied in two dimensional pseduspin-1 optical lattices . While Lieb lattice incorporates both of these features, a simple dice lattice lacks the later one. However, if the bulk band dispersion can be made topologically nontrivial, dice lattice would have advantage of having large topological charge than Lieb lattice. This problem can be addressed for dice lattice in half-metallic or spin-gapless ferromagnetic condensed systems with electron spin polarization P = 1, unlike regular ferromagnetic materials where 0 < P < 1. Keeping this microscopic feature of intrinsic magnetization in view, spin-polarized S=1 Dirac dispersion in dice lattice can be modelled through a three-band tight binding model (1). When only nearest-neighbour hoping is considered, single-particle band structure represents either a S=1 Dirac spin-gapless semiconductor or S=1 Dirac half-metal depending upon the on-site energies. When spin-polarized Haldane-like next-nearest-neighbour tunnelling is also considered, flat band at the Fermi level remains unaffected while a topological gap opens between dispersing valance and conduction bands. As a result, a spin-polarized S=1 massive Dirac dispersion gives topologically non-trivial bulk band dispersion with Chern number -2. Because of the presence of flat band, we refer these S=1 phases as flat band Dirac half-metals (FB-DHM) and flat band spin-gapless semiconductor (FB-SGS). While Dirac half-metal (DHM) and spin-gapless semiconductors (SGS) host spin-polarized S=1/2 Dirac dispersion, FB-DHM and FB-SGS hosts spin-polarized S=1 Dirac dispersion where a spin-polarized S=1/2 Dirac dispersion is accompanied by a single spin flat band. FB-SGS, where spin polarized flat band is lying exactly at the Fermi level, provides an ideal platform for the generalization of following two very important concepts previously studied in S=1/2 condensed matter systems – intrinsic ferromagnetism in two-dimensional materials and quantum anomalous Hall (QAH) effect in SGS – to S=1 condensed matter systems. In realistic condensed matter systems, transition metal dihalides MX monolayers, where transition metal cation M = V, Cr, Mn, Fe, Co, Ni is surrounded by six halogen anions X=Br, Cl, I in octahedral environment, are much simpler two-dimensional systems that can host FB-SGS, FB-DHM and spin-polarized S=1 massive Dirac dispersion. It has been repotred that VX , CrX , and MnX monolayers are antiferromagnetic while FeX , NiX , CoBr and CoCl monolayers are ferromagnetic. Moreover, spin-resolved density of states shows that M +2 (Fe +2 (3d ), Co +2 (3d ), Ni +2 (3d )) ions in MX monolayers exhibit high-spin configuration: one of the spin sectors (say spin up) is distributed over whole d -manifold while the other spin sector (spin-down) fills the low-lying t orbitals with one, two, and three electrons in FeX , CoX , and NiX monolayers respectively. Such a spin configuration of d -manifold renders FeX monolayers as conventional ferromagnetic half-metals , CoBr and CoCl monolayers as ferromagnetic semiconductors with a very small band gap , while NiX monolayers as ferromagnetic insulators . FeX and CoX monolayers are of particular interest among these MX family members for S=1 physics, mainly because of partially filled low-lying t orbitals. Here, by using first principle calculations, we show that a small compressive strain leads FeBr and FeCl monolayers to FB-DHM while CoBr and CoCl monolayers to FB-SGS. We further investigate the intrinsic ferromagnetic and topological bulk properties in CoBr and CoCl monolayers and found that: (i) the ferromagnetic ground state of CoBr and CoCl monolayers embeds flat band ferromagnetism . (ii) spin-orbit coupling opens a topologically non-trivial Dirac gap between dispersing valance and conduction band while leaving flat band unaffected. This is in complete agreement with the effects of Haldane-like NNN tunnelling on the single-particle band structure obtained from three-band tight binding model on dice lattice. The combined effect of intrinsic flat-band ferromagnetism and nontrivial bulk electronic topological properties leads to quantum anomalous Hall state with Chern number C = -2. II.
THREE-BAND TIGHT BINDING MODEL
The dice lattice structure, with three-sites per unit cell, is composed of two inter penetrating honeycomb lattices, figure (1). While honeycomb lattice hosts S=1/2 Dirac dispersion, dice lattice is featured by S=1 Dirac dispersion. With ferromagnetic ground state, coupling between itinerant electron’s spin and magnetization leads the dice lattice to host spin-polarized S=1 Dirac dispersion. FB-DHM, FB-SGS and spin-polarized S=1 massive Dirac band dispersion can be modelled through three-band tight binding Hamiltonian on the dice lattice † † † †0 ( ) i i i j ij i j i i ii ij iij H M c c t c c i c c m c c (1)
Here † ( ) i i c c is the creation (annihilation) electron operator on site i , ( , , ) (1, 0, 1) i a c b represents the sign convention for staggered sublattice potential, and ij ik kj v d d connects sites i and j on sublattice A (and B) via the unique intermediate site k on sublattice C. Here ik d and kj d are the nearest-neighbour bond vectors connecting A (and B) sublattices with sublattice C. In the model Hamiltonian (1), first term incorporates both on-site energy ε and Zeeman exchange interaction M > 0, second term is the nearest neighbour hoping while third term is modified Kane and Mele type SO coupling for S=1 system. By setting electron spin polarization s z = -1, Kane and Mele type SO coupling term effectively reduces to Haldane-like term for next-nearest-neighbour tunnelling . The last term represents the staggered sublattice potential. In the long wavelength limit, effective Hamiltonian and low energy single-particle band dispersion in the vicinity of spin-polarized S=1 Dirac points are M m atqH K q atq M atqatq M m (2) ( ) E K q M (3a) ( ) 2 3 3 F E K q M v q m (3b) where F v at is Fermi velocity, x y q q q and
1( 1) is the valley index representing K ( K '). In the absence of SO coupling and staggered sublattice potential, three-fold degenerate crossing at K and K' represents S=1 spin-polarized Dirac points. That is, the Hamiltonian (1) represent FB-SGS and FB-DHM for ( ) 0 M and ( ) 0 M respectively, as shown in figure (1). When SO coupling and sublattice potential are activated, a Dirac gap of ( ) 2 3 3 g E m is opened between dispersing valance and conduction bands in FB-SGS. As shown in figure (2), spin-polarized S=1 massive Dirac dispersion captures the salient features of Haldane model for S=1/2 systems: the gap remains always positive at valley K ', however, it decreases at valley K with increasing the strength of staggered potential m, closes when m , and reopens for m . Band dispersion of semi-infinite nanoribbon confirms the non-zero Chern number, C = -2. In next section, we show that electronic band structure obtained from three-band tight-binding model for infinite dice sheet qualitatively reproduces the electronic properties of strained CoBr and CoCl monolayer predicted by first principle calculations and hence reveal the underlying topological physics effectively. Figure 1:
Lattice and single-particle band dispersion for three-band tight binding model.
A two-dimensional network of two intermingled honeycomb lattices, known as dice structure with three-sites per unit cell, hosting S=1 Dirac dispersion (a) . When m , three-band tight binding model (1) gives FB-SGS for ε – M = 0 (b) and FB-DHM for ε – M > 0 (c) . Figure 2:
Effects of SO coupling and staggered sublattice potential on FB-SGS. ( Upper panel ) Bulk band structure for infinite sheet of dice lattice, showing topological phase transition from QAH phase to trivial insulating phase. (
Lower panel ) Band structure for semi-infinite sheet with zigzag edges. For m , edge state crossing Fermi level shows nontrivial behaviour, C = -2. Here we set and ε – M = 0. III.
MATERIAL REALIZATION
It is worthy to shed some light on microscopic ferromagnetic phenomena for bridging between three-band tight binding model on dice lattice and first principle caculations for MX monolayers presented in this section. In a conventional ferromagnetic system, where electron spin polarization is larger than 0 but smaller than 1; 0 < P < 1, a realistic Zeeman-type exchange splitting could not completely distangle the two oppositely spin polarized bands near Fermi level. That is, if we include both spin sectors in Hamiltonian (1), S=1 spin-up(spin-down) Dirac points will be lying ε + M above ( ε - M below) the Fermi level and there would be band crossing between oppositely spin polarized bands at Fermi level which is located at the energy zero. As a result, when Zeeman coupling M is less than the bandwidth W of topologically non-trivial dispersing bands, even a spin-orbit coupled ferromagnetic insulator would turn into a ferromagnetic metal – a partial concern for the realization of QAH phase in conventional ferromagnetic systems. This problem can be naturally addressed in ferromagnetic spin-gapless semiconductors, ferromagnetic half-metals, or spin-orbit coupled ferromagnetic insulators where electron spin polarization is P = 1; oppositely spin-polarized sectors are already well separated in energy and only one of the spin-polarized bands reside at/near Fermi level. With ferromagnetic ground state, by setting all localized Ising spins directed out-of-plane (along z-axis) and considering itinerant electrons carry spin-down, coupling between itinerant electron’s spin and magnetization provided by localized spins leads spin-polarized band dispersion. Secondly, Hamiltonian (1) effectively describes the FB-DHM, FB-SGS and spin-polarized S=1 massive Dirac band dispersion with characteristics similar as we find in first principle calculations for transition metal dihalide monolayers. However, in realistic transition metal dihalides, the edges of conduction and valence bands near Fermi level are predominantly contributed by all three low-lying t orbitals of transition metal ions. Hence, a more accurate tight binding model should be constructed based on all three orbitals along with the consideration of symmetries of the monolayers. Moreover, for 3 d orbitals involving multiple bands, possible electron-electron interactions are required to be considered for complete and quantitative understanding of flat band magnetism. A. First principle calculations
The presence of FBSGS, FB-DHM, and spin-polarized S=1 massive Dirac dispersion with nontrivial bulk topology requires large intrinsic magnetization and strong SO coupling. We show that these features are embedded in transition metal dihalides MX monolayers. We study the magnetic and electronic properties of MX (M = Fe, Co; X=Br, Cl) monolayers with space group P-3m1 using the projected augmented plane-wave method (PAW) as implemented in the Vienna ab initio simulation package (VASP) . The Perdew-Burke-Ernzerhof functional of the generalized gradient approximation (GGA) is employed for the account of exchange-correlation interactions among electrons . The cutoff energy for the plane-wave-basis set is chosen to be 400 eV for all calculations. A 2×2×1 supercell is used to calculate the exchange energy of ferromagnetic and antiferromagnetic (with periodic boundary conditions) configurations of Co spins. A 35×35×1 k-point grid is employed to sample the Brillouin zone. The structure is obtained by fully relaxing the atomic positions until the Hellmann-Feynman forces on each ion are less than 0.01 eV/A -1 . Convergence of total energies and full optimization of the atomic coordinates yields the hexagonal lattice constants for FeX and CoX as shown in table 1. These values are in close agreement as previously reported . To study the massive Dirac dispersion, SO coupling is included in the self-consistent loop. MX Pristine Strained FeBr Table 1:
Calculated lattice parameter (a = b, in Å) for pristine and strained MX monolayers and the corresponding ferromagnetic electronic phases. Figure 3:
First principle calculation for FeBr , FeCl , CoBr , and CoCl monolayers. Red (blue) curves represent spin-down(spin-up) bands.
Top and middle panels show band dispersion for pristine and strained FeBr , FeCl , CoBr , and CoCl monolayers respectively. A small compressive strain leads FeBr and FeCl monolayers to S=1 Dirac half-metallic phase while CoBr and CoCl monolayers to S=1 spin-polarized Dirac gapless state. Lower panel shows partial density of states (PDOS) for strained FeBr , FeCl , CoBr , and CoCl monolayers. B. FB-SGS and FB-DHM
In terms of electronic structure, consistent with the previous results , first principle calculations show that FeX monolayers are conventional ferromagnetic half-metals while CoBr and CoCl monolayers are ferromagnetic semiconductors with a very small band gap . However, a compressive strain leads FeBr and FeCl monolayers to FB-DHM while CoBr and CoCl monolayers to FB-SGS. As shown in figure (3), the band at the Fermi level remains flat near the corners of Brillion zone (high symmetry points K and K') for CoBr monolayer while the band at the Fermi level remains flat along the boundary of Brillion zone (high symmetry line K-M-K') for CoCl monolayer. Like FeBr and FeCl , pristine FeI monolayer is also a conventional ferromagnetic half metal. However, a very large compressive strain is required to attain FB-DHM in FeI monolayer. Spin-polarized S=1 Dirac band dispersions for FB-DHM and FB-SGS can easily be attained from electronic band dispersion of pristine monolayers under very small and experimentally accessible compressive strains of ε = - 4.4%, ε = - 5%, ε = - 1.8%, and ε = - 1.1% for FeBr , FeCl , CoBr and CoCl monolayers respectively . C. Flat band ferromagnetism
We calculate the ground state energy for both ferromagnetic and periodic antiferromagnetic (stripe order) spin configuration for M X cell as shown in figure (4). The large positive values for exchange energy confirm that the ground state of FeBr , FeCl , CoBr and CoCl monolayers are ferromagnetic with a high Curie temperature. Ferromagnetic ordering is consistent with Goodenough–Kanamori–Anderson (GKA) rules . While the localized magnetic moments arise through the competition between ionic and covalent bond strength between M and X ions, the magnetic ordering originates from the interplay between superexchange and direct exchange interaction. Although both direct and superexchange leads to antiferromagnetic ordering in principle, GKA rules indicate that the dominating superexchange interaction need not to favour antiferromagnetic ordering in all cases: If d -orbitals on neighbouring M ions hybridize with orthogonal p orbitals on X ions such that M-X-M bond angles are 90 , the superexchange gives ferromagnetic ordering. Spin-resolved density of states shown in figure (3) reflect two key features of FeBr , FeCl , CoBr and CoCl monolayers: (i) With seven (six) 3 d electrons, Co +2 (Fe +2 ) ions exhibit high-spin configuration, five with spin-up and two (one) electrons with spin-down state. The calculated total magnetic moment of B ( B ) per unit cell, with most of it carrying by Co ≈ 2.48 (Fe ≈ 3.45) ions, confirms that the ground state favours high spin-configuration. In the octahedral crystal field environment, where five 3 d orbitals split into two e g and three t orbitals, five spin up electrons are distributed over whole 3 d manifold while two (one) spin-down electrons reside in low energy t orbitals. This spin configuration of Co +2 (Fe +2 ) electrons leaves CoBr and CoCl (FeCl and FeBr ) monolayers in a spin-polarized gapless (half metal) state where only spin-down electrons in t orbitals reside near the Fermi energy. (ii) A sharp peak of density of states at Fermi level indicates the presence of spin-polarized flat band and hence flat-band ferromagnetism in CoBr and CoCl monolayers. This prominent feature is the result of specific inter-site electron hoping in MX monolayers which is composed of two intermingled honeycomb lattices with unequal number of M and X sites . The effect of on-site correlations on magnetic ordering and spin-configuration of Co +2 (Fe +2 ) ions is studied through LDA+ U . The flat band ferromagnetism , for any positive value of on-site Coulomb correlations U ≥ 0, in CoBr and CoCl monolayers can be justified in the non-interacting or weakly-interacting limit: when flat band at the Fermi level is saturated with fully spin-polarized electrons, the Pauli exclusion of double occupancy minimizes the Coulomb interaction and renders ground state of CoBr and CoCl monolayer ferromagnetic. Figure 4:
Magnetic configuration and effects of SO coupling in CoBr and CoCl monolayers. (a, b) Ferromagnetic (a) and periodic antiferromagnetic (b) spin configuration for MX . Blue and green arrows represent the mechanism for ferromagnetic superexchange interaction with 90 of M-X-M bond angle and antiferromagnetic direct M-M exchange interaction between localized moments respectively. (c, d) Band structure of CoBr (c) and CoCl (d) monolayers with SO coupling. SO coupling open a gap between dispersing valance and conduction band in whole Brillion zone for both CoBr and CoCl monolayers. D. Spin-polarized S=1 massive Dirac dispersion - Effect of SO coupling
Having established spin-polarized S=1 Dirac character of CoBr and CoCl monolayers with intrinsic flat band ferromagnetism, we explore the possible effects of SO coupling on both Dirac dispersing and flat band as shown in figure (4). In the non-interacting limit, a gap is opened between dispersing valance and conduction bands by the activation of SO coupling. However, spin-polarized flat band lying at Fermi level remains unaffected. That is, flat band gets detached from both dispersing valance and conduction bands but remains at the Fermi level. This is in complete agreement with the single-particle band structure obtained from three-band tight binding model on dice lattice. Simultaneous presence of flat band ferromagnetism and topological spin-orbit coupling gap in FB-SGS is promising for the exploration of QAH phase when Fermi energy E F is positioned within the Dirac gap between dispersing bands. Possibility of QAH phase, with a trademark of dissipationless chiral edge states, in two-dimensional S=1 spin-polarized Dirac systems has an obvious advantage over QAH phase in S=1/2 magnetic doped topological insulators (TIs) . The idea of achieving ferromagnetic ordering in TIs by inducing magnetic impurities , has already led to the experimental realization of QAH phase in Cr and V doped thin films of tetradymite semiconductors (Bi,Sb) Te . However, the highest temperature at which QAH phase is observed in magnetic doped TIs is only few Kelvin. Several experimental challenges limit the high temperature realization of QAH phase in magnetic doped TIs. For example, presence of spin scattering centres due to non-uniform distribution of magnetic impurities and parallel dissipative channels at the one-dimensional edges are likely to induce several KΩ longitudinal resistance. Moreover, as doped magnetic impurities are required to generate a sufficiently strong magnetic moment without effecting the band topology of the system, difficulty in simultaneous achievement of high T c and a large topological gap opening in magnetically doped TIs limit the applications of high temperature QAH phase . On the other hand, FB-SGS with intrinsic flat band ferromagnetism, sustainable against thermal fluctuations at high temperature in the two-dimensional limit, in combination with topologically non-trivial bulk band structure can offer QAH phase without requiring random magnetic impurity doping. With quantized transverse resistance, extremely low longitudinal resistance, and large topological charge could allow possible future dissipationless spintronics and low-energy electronics technologies in two-dimensional FB-SGS. IV.
DISCUSSION
We proposed a three-band tight binding model, with spin-polarized Haldane-like next-nearest-neighbour tunnelling, for ferromagnetic condensed matter systems with electron spin polarization P = 1. By using first principle calculations, we further show that FeBr and FeCl monolayers are FB-DHM while CoBr and CoCl monolayers are FB-SGS with spin-polarized S=1 Dirac dispersion. When SO coupling is activated, like the effects of Haldane-like NNN tunnelling on the single-particle band structure obtained from three-band tight binding model on dice lattice, a topologically non-trivial Dirac gap is opened between dispersing valance and conduction band while leaving flat band unaffected. FB-SGS embeds two interesting phenomena and has an obvious advantage over S=1/2 SGS; the presence of spin-polarized flat band at the Fermi level leads to two-dimensional intrinsic flat band ferromagnetism and topologically nontrivial electronic properties with large topological charge. Since the physical realization of quantum anomalous Hall effect is indebted to the existence of intrinsic ferromagnetism, this study greatly extends the search for material classes of quantum anomalous Hall phase from S=1/2 to S=1 condensed matter system. Moreover, in the presence of strong SO coupling, intrinsic flat band ferromagnetism can serve as guiding principle for the experimental realization of other time reversal symmetry breaking topological phases such as topological superconductors. Finally, intrinsic flat band ferromagnetism and topologically nontrivial electronic properties in two-dimensional S=1 condensed matter system is promising for playing vital role in the realization of dissipationless spintronics and electronics technologies.
1 Lv, B. Q. et al.
Experimental Discovery of Weyl Semimetal TaAs.
Physical Review X , 031013, doi:10.1103/PhysRevX.5.031013 (2015). 2 Lv, B. Q. et al. Observation of Weyl nodes in TaAs.
Nature Physics , 724, doi:10.1038/nphys3426 (2015). 3 Xu, S.-Y. et al. Discovery of a Weyl fermion semimetal and topological Fermi arcs.
Science , 613 (2015). 4 Yang, L. X. et al.
Weyl semimetal phase in the non-centrosymmetric compound TaAs.
Nature Physics et al. Discovery of a Weyl fermion state with Fermi arcs in niobium arsenide.
Nature Physics , 748, doi:10.1038/nphys3437 (2015). 6 Liu, Z. K. et al. A stable three-dimensional topological Dirac semimetal Cd3As2.
Nature Materials et al. Discovery of a Three-Dimensional Topological Dirac Semimetal, Na3Bi.
Science , 864, doi:10.1126/science.1245085 (2014). 8 Lutchyn, R. M., Sau, J. D. & Das Sarma, S. Majorana Fermions and a Topological Phase Transition in Semiconductor-Superconductor Heterostructures.
Physical Review Letters , 077001, doi:10.1103/PhysRevLett.105.077001 (2010). 9 Mourik, V. et al.
Signatures of Majorana Fermions in Hybrid Superconductor-Semiconductor Nanowire Devices.
Science , 1003, doi:10.1126/science.1222360 (2012). 10 Nadj-Perge, S. et al.
Observation of Majorana fermions in ferromagnetic atomic chains on a superconductor.
Science , 602, doi:10.1126/science.1259327 (2014). 11 He, Q. L. et al.
Chiral Majorana fermion modes in a quantum anomalous Hall insulator–superconductor structure.
Science , 294, doi:10.1126/science.aag2792 (2017). 12 Fu, L. & Kane, C. L. Superconducting proximity effect and majorana fermions at the surface of a topological insulator.
Phys Rev Lett , 096407, doi:10.1103/PhysRevLett.100.096407 (2008). 13 Wan, X., Turner, A. M., Vishwanath, A. & Savrasov, S. Y. Topological semimetal and Fermi-arc surface states in the electronic structure of pyrochlore iridates.
Physical Review B , 205101, doi:10.1103/PhysRevB.83.205101 (2011). 14 Young, S. M. et al. Dirac semimetal in three dimensions.
Phys Rev Lett , 140405, doi:10.1103/PhysRevLett.108.140405 (2012). 15 Bradlyn, B. et al.
Beyond Dirac and Weyl fermions: Unconventional quasiparticles in conventional crystals.
Science (2016). 16 Wieder, B. J., Kim, Y., Rappe, A. M. & Kane, C. L. Double Dirac Semimetals in Three Dimensions.
Physical Review Letters , 186402, doi:10.1103/PhysRevLett.116.186402 (2016). 17 Lv, B. Q. et al.
Observation of three-component fermions in the topological semimetal molybdenum phosphide.
Nature , 627, doi:10.1038/nature22390 (2017). 18 Chang, G. et al.
Unconventional Chiral Fermions and Large Topological Fermi Arcs in RhSi.
Physical Review Letters , 206401, doi:10.1103/PhysRevLett.119.206401 (2017). 19 Tang, P., Zhou, Q. & Zhang, S.-C. Multiple Types of Topological Fermions in Transition Metal Silicides.
Physical Review Letters , 206402, doi:10.1103/PhysRevLett.119.206402 (2017). 20 Takane, D. et al.
Observation of Chiral Fermions with a Large Topological Charge and Associated Fermi-Arc Surface States in CoSi.
Physical Review Letters , 076402, doi:10.1103/PhysRevLett.122.076402 (2019). 21 Yang, Y. et al.
Topological triply degenerate point with double Fermi arcs.
Nature Physics , doi:10.1038/s41567-019-0502-z (2019). 22 Sanchez, D. S. et al.
Topological chiral crystals with helicoid-arc quantum states.
Nature , 500-505, doi:10.1038/s41586-019-1037-2 (2019). 23 Rao, Z. et al.
Observation of unconventional chiral fermions with long Fermi arcs in CoSi.
Nature , 496-499, doi:10.1038/s41586-019-1031-8 (2019). 24 Bercioux, D., Goldman, N. & Urban, D. F. Topology-induced phase transitions in quantum spin Hall lattices.
Physical Review A , 023609, doi:10.1103/PhysRevA.83.023609 (2011). 25 Weeks, C. & Franz, M. Topological insulators on the Lieb and perovskite lattices. Physical Review B , 085310, doi:10.1103/PhysRevB.82.085310 (2010). 26 Goldman, N., Urban, D. F. & Bercioux, D. Topological phases for fermionic cold atoms on the Lieb lattice. Physical Review A , 063601, doi:10.1103/PhysRevA.83.063601 (2011). 27 Haldane, F. D. M. Model for a Quantum Hall Effect without Landau Levels: Condensed-Matter Realization of the "Parity Anomaly". Physical Review Letters , 2015-2018, doi:10.1103/PhysRevLett.61.2015 (1988). 28 Ishizuka, H. & Motome, Y. Dirac Half-Metal in a Triangular Ferrimagnet. Physical Review Letters , 237207, doi:10.1103/PhysRevLett.109.237207 (2012). 29 Wang, X. L. Proposal for a New Class of Materials: Spin Gapless Semiconductors.
Physical Review Letters , 156404, doi:10.1103/PhysRevLett.100.156404 (2008). 30 Gong, C. et al.
Discovery of intrinsic ferromagnetism in two-dimensional van der Waals crystals.
Nature , 265, doi:10.1038/nature22060 (2017). 31 Huang, B. et al.
Layer-dependent ferromagnetism in a van der Waals crystal down to the monolayer limit.
Nature , 270, doi:10.1038/nature22391 (2017). 32 Bonilla, M. et al.
Strong room-temperature ferromagnetism in VSe2 monolayers on van der 0 Waals substrates.
Nature Nanotechnology , 289-293, doi:10.1038/s41565-018-0063-9 (2018). 33 Shabbir, B. et al. Long range intrinsic ferromagnetism in two dimensional materials and dissipationless future technologies.
Applied Physics Reviews , 041105, doi:10.1063/1.5040694 (2018). 34 Ye, L. et al. Massive Dirac fermions in a ferromagnetic kagome metal.
Nature , 638, doi:10.1038/nature25987 (2018). 35 Wang, X.-L. Dirac spin-gapless semiconductors: promising platforms for massless and dissipationless spintronics and new (quantum) anomalous spin Hall effects.
National Science Review , 252-257, doi:10.1093/nsr/nww069 (2016). 36 He, J., Li, X., Lyu, P. & Nachtigall, P. Near-room-temperature Chern insulator and Dirac spin-gapless semiconductor: nickel chloride monolayer. Nanoscale , 2246-2252, doi:10.1039/C6NR08522A (2017). 37 Kulish, V. V. & Huang, W. Single-layer metal halides MX2 (X = Cl, Br, I): stability and tunable magnetism from first principles and Monte Carlo simulations. Journal of Materials Chemistry C , 8734-8741, doi:10.1039/C7TC02664A (2017). 38 Botana, A. S. & Norman, M. R. Electronic structure and magnetism of transition metal dihalides: Bulk to monolayer. Physical Review Materials , 044001, doi:10.1103/PhysRevMaterials.3.044001 (2019). 39 Ashton, M. et al. Two-Dimensional Intrinsic Half-Metals With Large Spin Gaps.
Nano Letters , 5251-5257, doi:10.1021/acs.nanolett.7b01367 (2017). 40 Mounet, N. et al. Two-dimensional materials from high-throughput computational exfoliation of experimentally known compounds.
Nature Nanotechnology , 246-252, doi:10.1038/s41565-017-0035-5 (2018). 41 Lieb, E. H. Two theorems on the Hubbard model. Physical Review Letters , 1201-1204, doi:10.1103/PhysRevLett.62.1201 (1989). 42 Tasaki, H. Ferromagnetism in the Hubbard models with degenerate single-electron ground states. Physical Review Letters , 1608-1611, doi:10.1103/PhysRevLett.69.1608 (1992). 43 Mielke, A. & Tasaki, H. Ferromagnetism in the Hubbard model. Examples from models with degenerate single-electron ground states. Comm. Math. Phys. , 341-371 (1993). 44 Tasaki, H. Ferromagnetism in the Hubbard Model: A Constructive Approach.
Communications in Mathematical Physics , 445-472, doi:10.1007/s00220-003-0952-z (2003). 45 Kane, C. L. & Mele, E. J. Quantum spin Hall effect in graphene.
Phys Rev Lett , 226801, doi:10.1103/PhysRevLett.95.226801 (2005). 46 Kane, C. L. & Mele, E. J. Z2 topological order and the quantum spin Hall effect. Phys Rev Lett , 146802, doi:10.1103/PhysRevLett.95.146802 (2005). 47 Blöchl, P. E. Projector augmented-wave method. Physical Review B , 17953-17979, doi:10.1103/PhysRevB.50.17953 (1994). 48 Kresse, G. & Joubert, D. From ultrasoft pseudopotentials to the projector augmented-wave method. Physical Review B , 1758-1775, doi:10.1103/PhysRevB.59.1758 (1999). 49 Kresse, G. & Hafner, J. Ab initio molecular-dynamics simulation of the liquid-metal--amorphous-semiconductor transition in germanium. Physical Review B , 14251-14269, doi:10.1103/PhysRevB.49.14251 (1994). 50 Perdew, J. P., Burke, K. & Ernzerhof, M. Generalized Gradient Approximation Made Simple. Physical Review Letters , 3865-3868, doi:10.1103/PhysRevLett.77.3865 (1996). 51 Goodenough, J. B. Theory of the Role of Covalence in the Perovskite-Type Manganites [La, M(II}]MnO3. Physical Review , 564-573, doi:10.1103/PhysRev.100.564 (1955). 52 Anderson, P. W. New Approach to the Theory of Superexchange Interactions.
Physical Review , 2-13, doi:10.1103/PhysRev.115.2 (1959). 53 Kanamori, J. Crystal Distortion in Magnetic Compounds.
Journal of Applied Physics , S14-S23, doi:10.1063/1.1984590 (1960). 1 54 Yu, R. et al. Quantized Anomalous Hall Effect in Magnetic Topological Insulators.
Science , 61, doi:10.1126/science.1187485 (2010). 55 Wang, Z. F., Liu, Z. & Liu, F. Quantum Anomalous Hall Effect in 2D Organic Topological Insulators.
Physical Review Letters , 196801 (2013). 56 Fang, C., Gilbert, M. J. & Bernevig, B. A. Large-Chern-Number Quantum Anomalous Hall Effect in Thin-Film Topological Crystalline Insulators.
Physical Review Letters , 046801 (2014). 57 Chang, C.-Z. et al.
Experimental Observation of the Quantum Anomalous Hall Effect in a Magnetic Topological Insulator.
Science , 167, doi:10.1126/science.1234414 (2013). 58 Bestwick, A. J. et al.
Precise Quantization of the Anomalous Hall Effect near Zero Magnetic Field.
Physical Review Letters , 187201, doi:10.1103/PhysRevLett.114.187201 (2015). 59 Chang, C.-Z. et al.
High-precision realization of robust quantum anomalous Hall state in a hard ferromagnetic topological insulator.
Nature Materials Acknowledgments