## Abstract

In recent years, transition metal dichalcogenides (TMDs) have garnered great interest as topological materials. In particular, monolayers of centrosymmetric β-phase TMDs have been identified as 2D topological insulators (TIs), and bulk crystals of noncentrosymmetric γ-phase MoTe2 and WTe2 have been identified as type-II Weyl semimetals. However, angle-resolved photoemission spectroscopy and STM probes of these semimetals have revealed huge, arclike surface states that overwhelm, and are sometimes mistaken for, the much smaller topological surface Fermi arcs of bulk type-II Weyl points. In this Letter, we calculate the bulk and surface electronic structure of both β- and γ-MoTe2. We find that β-MoTe2 is, in fact, a Z4-nontrivial higher-order TI (HOTI) driven by double band inversion and exhibits the same surface features as γ-MoTe2 and γ-WTe2. We discover that these surface states are not topologically trivial, as previously characterized by the research that differentiated them from the Weyl Fermi arcs but, rather, are the characteristic split and gapped fourfold Dirac surface states of a HOTI. In β-MoTe2, this indicates that it would exhibit helical pairs of hinge states if it were bulk insulating, and in γ-MoTe2 and γ-WTe2, these surface states represent vestiges of HOTI phases without inversion symmetry that are nearby in parameter space. Using nested Wilson loops and first-principles calculations, we explicitly demonstrate that, when the Weyl points in γ-MoTe2 are annihilated, which may be accomplished by symmetry-preserving strain or lattice distortion, γ-MoTe2 becomes a nonsymmetry-indicated, noncentrosymmetric HOTI. We also show that, when the effects of spin-orbit coupling are neglected, β-MoTe2 is a nodal-line semimetal with Z2-nontrivial monopole nodal lines (MNLSM). This finding confirms that MNLSMs driven by double band inversion are the weak-spin-orbit coupling limit of HOTIs, implying that MNLSMs are higher-order topological semimetals with flat-band-like hinge states, which we find to originate from the corner modes of 2D "fragile" TIs.

Original language | English (US) |
---|---|

Article number | 186401 |

Journal | Physical review letters |

Volume | 123 |

Issue number | 18 |

DOIs | |

State | Published - Oct 28 2019 |

## All Science Journal Classification (ASJC) codes

- Physics and Astronomy(all)

## Access to Document

## Other files and links

## Fingerprint Dive into the research topics of 'Higher-Order Topology, Monopole Nodal Lines, and the Origin of Large Fermi Arcs in Transition Metal Dichalcogenides XTe2 (X=Mo, W)'. Together they form a unique fingerprint.

## Cite this

- APA
- Author
- BIBTEX
- Harvard
- Standard
- RIS
- Vancouver

*Physical review letters*,

*123*(18), [186401]. https://doi.org/10.1103/PhysRevLett.123.186401

**Higher-Order Topology, Monopole Nodal Lines, and the Origin of Large Fermi Arcs in Transition Metal Dichalcogenides XTe2 (X=Mo, W)**. In: Physical review letters. 2019 ; Vol. 123, No. 18.

}

*Physical review letters*, vol. 123, no. 18, 186401. https://doi.org/10.1103/PhysRevLett.123.186401

**Higher-Order Topology, Monopole Nodal Lines, and the Origin of Large Fermi Arcs in Transition Metal Dichalcogenides XTe2 (X=Mo, W).** / Wang, Zhijun; Wieder, Benjamin J.; Li, Jian; Yan, Binghai; Bernevig, B. Andrei.

Research output: Contribution to journal › Article › peer-review

TY - JOUR

T1 - Higher-Order Topology, Monopole Nodal Lines, and the Origin of Large Fermi Arcs in Transition Metal Dichalcogenides XTe2 (X=Mo, W)

AU - Wang, Zhijun

AU - Wieder, Benjamin J.

AU - Li, Jian

AU - Yan, Binghai

AU - Bernevig, B. Andrei

N1 - Funding Information: The authors thank Barry Bradlyn, Jennifer Cano, Yichen Hu, Fan Zhang, and Youngkuk Kim for helpful discussions. Z.W., B.J.W., and B.A.B. were supported by the Department of Energy Grant No. DE-SC0016239, the National Science Foundation EAGER Grants No. DMR 1643312 and NSF-MRSEC No. DMR-1420541, Army Research Office Grant No. ARO MURI W911NF-12-1-0461, Simons Investigator Grant No. 404513, ONR Grant No. N00014-14-1-0330, the Packard Foundation, the Schmidt Fund for Innovative Research, and a Guggenheim Fellowship from the John Simon Guggenheim Memorial Foundation. Z.W. additionally acknowledges support from the National Thousand-Young-Talents Program and the CAS Pioneer Hundred Talents Program. This work was also supported by the National Natural Science Foundation of China (Grants No. 11504117 and No. 11774317). B.Y. acknowledges support from the Willner Family Leadership Institute for the Weizmann Institute of Science, the Benoziyo Endowment Fund for the Advancement of Science, the Ruth and Herman Albert Scholars Program for New Scientists, and the European Research Council (ERC) under the European Union Horizon 2020 Research and Innovation Programme (Grant No. 815869). Funding Information: X Te 2 ( X = Mo , W ) Wang Zhijun 1,2,3 https://orcid.org/0000-0003-2540-6202 Wieder Benjamin J. 3 Li Jian 4 Yan Binghai 5 Bernevig B. Andrei 3,6,7 Beijing National Laboratory for Condensed Matter Physics, and Institute of Physics, 1 Chinese Academy of Sciences , Beijing 100190, China 2 University of Chinese Academy of Sciences , Beijing 100049, China Department of Physics, 3 Princeton University , Princeton, New Jersey 08544, USA School of Science, 4 Westlake University , 18 Shilongshan Road, Hangzhou 310024, China Department of Condensed Matter Physics, 5 Weizmann Institute of Science , Rehovot 7610001, Israel Dahlem Center for Complex Quantum Systems and Fachbereich Physik, 6 Freie Universität Berlin , Arnimallee 14, 14195 Berlin, Germany 7 Max Planck Institute of Microstructure Physics , 06120 Halle, Germany 28 October 2019 1 November 2019 123 18 186401 16 September 2019 2 July 2018 © 2019 American Physical Society 2019 American Physical Society In recent years, transition metal dichalcogenides (TMDs) have garnered great interest as topological materials. In particular, monolayers of centrosymmetric β -phase TMDs have been identified as 2D topological insulators (TIs), and bulk crystals of noncentrosymmetric γ -phase MoTe 2 and WTe 2 have been identified as type-II Weyl semimetals. However, angle-resolved photoemission spectroscopy and STM probes of these semimetals have revealed huge, arclike surface states that overwhelm, and are sometimes mistaken for, the much smaller topological surface Fermi arcs of bulk type-II Weyl points. In this Letter, we calculate the bulk and surface electronic structure of both β - and γ - MoTe 2 . We find that β - Mo Te 2 is, in fact, a Z 4 -nontrivial higher-order TI (HOTI) driven by double band inversion and exhibits the same surface features as γ - MoTe 2 and γ - WTe 2 . We discover that these surface states are not topologically trivial, as previously characterized by the research that differentiated them from the Weyl Fermi arcs but, rather, are the characteristic split and gapped fourfold Dirac surface states of a HOTI. In β - MoTe 2 , this indicates that it would exhibit helical pairs of hinge states if it were bulk insulating, and in γ - MoTe 2 and γ - WTe 2 , these surface states represent vestiges of HOTI phases without inversion symmetry that are nearby in parameter space. Using nested Wilson loops and first-principles calculations, we explicitly demonstrate that, when the Weyl points in γ - MoTe 2 are annihilated, which may be accomplished by symmetry-preserving strain or lattice distortion, γ - MoTe 2 becomes a nonsymmetry-indicated, noncentrosymmetric HOTI. We also show that, when the effects of spin-orbit coupling are neglected, β - MoTe 2 is a nodal-line semimetal with Z 2 -nontrivial monopole nodal lines (MNLSM). This finding confirms that MNLSMs driven by double band inversion are the weak-spin-orbit coupling limit of HOTIs, implying that MNLSMs are higher-order topological semimetals with flat-band-like hinge states, which we find to originate from the corner modes of 2D “fragile” TIs. U.S. Department of Energy 10.13039/100000015 DE-SC0016239 National Science Foundation 10.13039/100000001 DMR 1643312 DMR-1420541 Army Research Office 10.13039/100000183 W911NF-12-1-0461 Simons Investigator 404513 Office of Naval Research 10.13039/100000006 N00014-14-1-0330 David and Lucile Packard Foundation 10.13039/100000008 Schmidt Fund for Innovative Research John Simon Guggenheim Memorial Foundation 10.13039/100005851 Chinese Academy of Sciences 10.13039/501100002367 National Natural Science Foundation of China 10.13039/501100001809 11504117 11774317 Weizmann Institute of Science 10.13039/501100001735 Benoziyo Endowment Fund Ruth and Herman Albert Scholars Program H2020 European Research Council 10.13039/100010663 815869 Over the past decade, the number of topological insulating and semimetallic phases identified in real materials has grown immensely. With the discovery of increasingly intricate topological insulators (TIs) and semimetals (SMs), such as higher-order TIs (HOTIs) [1–12] and unconventional fermion SMs [13–18] , previously overlooked, but readily accessible compounds, including bismuth [19] and chiral B20 [20–24] crystals, have been experimentally verified as topologically nontrivial. In this Letter, we extend the theory and experimental applicability of higher-order topology by recognizing that the X Te 2 ( X = Mo , W ) family [25–30] of transition metal dichalcogenides (TMDs), a large, well-studied, and readily synthesizable class of materials, are HOTIs, and not topologically trivial. X Te 2 TMDs, which exhibit a trivial magnetoelectric polarizability [11,31–38] , therefore, provide an intriguing experimental platform for examining topological response effects beyond the magnetoelectric effect [39–43] . We also draw connections between an exotic class of nodal-line SMs (NLSMs) [44–46] and the recently introduced notion of “fragile” topology [11,47–51] , leading to the discovery of fractionally charged corner modes in fragile TIs. Our findings further establish the seemingly esoteric concept of higher-order topology as crucial for characterizing topological transport and response effects in everyday materials. All of the spinful TIs discovered to date represent the gapped, spin-orbit coupled (strong-SOC) limits of gapless topological (SM) phases without SOC. This intrinsic link between gapped and gapless phases has defined topological condensed matter physics since the recognition that graphene [52–54] and HgTe gap into Z 2 TIs [55,56] under the introduction of SOC [57] . As the number of known topological SMs has increased [13–18,58–70] , the number of known topological (crystalline) insulators realized by gapping them has kept pace [7,71–76] . In one particularly simple example, an SM phase with a ring of linearly dispersing degeneracies, known as a “nodal line” (NL), can occur in a weak-SOC crystal with only inversion ( I ) and time-reversal ( T ) symmetries [44,77,78] . NLs may be created and annihilated at single, time-reversal-invariant (TRIM) points in the Brillouin zone (BZ) by inverting bands with opposite parity ( I ) eigenvalues [44] , such that the number of NLs is constrained by the same Fu-Kane parity criterion [79,80] that indicates 3D TI phases in strong-SOC crystals. This recognition has driven the rapid identification of candidate NLSMs, including Ca 3 P 2 [81] , Cu 3 ( Pd , Zn ) N [44,77] , and 3D graphene networks [82] , all of which exhibit characteristic nearly flat-band “drumhead” surface states. Crucially, it also implies that weak-SOC NLSMs gap directly into 3D TIs upon the introduction of I -symmetric SOC [44] . Very recently, fundamentally distinct I - and T -symmetric SMs and insulators have been proposed that escape this paradigm. In [45] , Fang et al. introduced a second kind of weak-SOC NL, which, unlike the previous example, can only be removed by pairwise annihilation. Though the mechanisms underpinning the protection and identification of these “monopole-charged” NLs (MNLs) have been explored in detail [46,83,84] ; MNLs have, thus far, only been proposed in magnonic systems [85] and 3D graphdiyne [46,86] . Recent works have also identified 3D “higher-order” TIs [1–12] stabilized by only I and T symmetries [3,19,87–90] . Notably, 3D HOTIs exhibit gapped 2D surfaces and gapless 1D hinges with characteristic helical modes [2–4] , which represent the domain wall states between 2D faces with oppositely gapped fourfold Dirac fermions [7,89] . By enumerating the parity eigenvalues of trivial (atomic) insulators, whose occupied bands define “elementary” band representations (EBRs) [48,91–96] , it can be shown that the Z 2 Fu-Kane criterion should be promoted to a Z 4 index that captures both TIs and HOTIs [3,19,87–89] . Using EBRs [93] , HOTI phases have been identified in systems with double band inversion (DBI) (Fig. 1 ), most notably in rhombohedral bismuth crystals [19] . 1 10.1103/PhysRevLett.123.186401.f1 FIG. 1. When two pairs of degenerate bands with positive parity eigenvalues and two pairs with negative parity eigenvalues are inverted at a TRIM point [3,19] , the occupied bands cannot be expressed as a linear combination of EBRs [48,91–96] and the Z 4 topological index [87–90] is changed by 2. In a T -symmetric crystal with vanishing SOC, this process may nucleate a pair of Dirac nodal lines with nontrivial monopole charge (MNLs) [46] (dashed lines in left panel). On the 1D hinges, the projections of the MNLs are spanned by nearly flat bands (an explicit model is provided in the Supplemental Material, Sec. A [97] ). These hinge states represent an example of higher-order topology in a bulk-gapless system: they are the d - 2 -dimensional generalization of drumhead surface states, and are the spinless analogs of the hinge states recently predicted in spinful Dirac SMs [50,119] . When I -symmetric SOC is introduced, the MNLSM will necessarily gap into a HOTI if all other bands are uninverted, and the flat-band hinge states will open into helical pairs spanning the bulk and surface gaps. HOTIs driven by this “double band inversion” (DBI) include bismuth [19] and, as shown in this Letter, β - MoTe 2 [Fig. 2(d) ]. Employing this Z 4 index and Wilson loops [120–124] , we identify the TMD [25] β - (1T’-) MoTe 2 [space group (SG) 11 P 2 1 / m ] as a HOTI with large, gapped, arclike surface states, and explicitly show that, in the absence of SOC, it forms an NLSM with MNLs. This cements the suggestion, introduced in [125] , that monopole nodal-line semimetals (MNLSMs) formed from DBI are the weak-SOC limit of HOTIs, in analogy to the earlier recognition [44] that monopole-trivial NLSMs are the weak-SOC limit of 3D TIs. It has been shown that the noncentrosymmetric γ - (Td-) phases of X Te 2 ( X = Mo , W ) (SG 31 P m n 2 1 ), previously identified as type-II (tilted) Weyl (semi)metals [68–70] , exhibit the same large surface states as β - MoTe 2 . These states were previously identified as topologically trivial [70,126–137] [Fig. 4(d) ], as the actual topological Fermi arcs from the Weyl points are considerably shorter [68,70,126,127,129,130] . However, using nested Wilson loops and first-principles calculations (Supplemental Material, Sec. B4 [97] ), we explicitly demonstrate that γ - MoTe 2 , which exhibits the same band ordering and surface states as β - MoTe 2 , transitions into a nonsymmetry-indicated HOTI when its narrowly separated Weyl points are annihilated. Therefore, the large surface states in γ - X Te 2 are not trivial but, rather, are vestiges of a nearby HOTI phase, and originate from DBI, like those in β - MoTe 2 . TMDs are a class of readily synthesizable [25–30] layered materials. Originally highlighted for the semiconducting band gap of exfoliated monolayers [138] , TMDs have recently been recognized as topological materials—quasi-2D samples of β -phase TMDs have been identified as 2D TIs [139–141] , and 3D samples of γ - X Te 2 have been identified in theory [68–70] and experiment [142–144] as type-II Weyl SMs. First, we focus on MoTe 2 , and then generalize our findings to the isostructural phases of WTe 2 . MoTe 2 can crystallize in two distinct structures at room temperature: the hexagonal α ( 2 H ) phase (SG 194 P 6 3 / m m c ) and the distorted monoclinic β phase [25–27,135] [Fig. 2(a) ]. When β - MoTe 2 is further cooled below 250 K, it transitions into the noncentrosymmetric γ phase [25,135,145] . Using first-principles calculations detailed in the Supplemental Material, Sec. B1 [97] , we calculate the electronic structure of β - MoTe 2 with and without the effects of SOC incorporated [Figs. 2(d) and 2(b) , respectively]. β - MoTe 2 exhibits DBI (Fig. 1 ) at Γ as a consequence of the β -phase lattice distortion [Fig. 2(a) ]. When SOC is neglected, a T -reversed pair of topological NLs [44,45] forms, intersecting Y Γ in an irregular, 3D shape with a significant pucker in the k y direction [schematically depicted in Fig. 2(c) ]. The NLs [red dots in Fig. 2(b) ] represent the only crossing points between the bands at E F (taking the direct gap to lie above N = 28 spin-degenerate pairs of bands). As prescribed in [46,83,84] , we surround each NL with a closed surface and calculate the Wilson loop (holonomy) matrix [120–124] over the lower N bands as a function of the polar momentum k θ [Fig. 2(c) and Supplemental Material, Sec. B2 [97] ]. This Wilson spectrum exhibits the characteristic winding of a MNL [46,83,84] . 2 10.1103/PhysRevLett.123.186401.f2 FIG. 2. (a) The monoclinic lattice of β - (1T’-) MoTe 2 [25] in SG 11 P 2 1 / m . (b) and (d) Bulk bands of β - MoTe 2 calculated without and with the effects of SOC incorporated, respectively (details provided in the Supplemental Material, Sec. B [97] ). DBI occurs about the Γ point, as indicated by the parity eigenvalues in (d). (c) When SOC is neglected, a T -reversed pair of irregularly shaped NLs forms at E F , intersecting Y Γ [red dots in (b)] between k y = ± 0.06 ( 2 π / b ) and k y = ± 0.19 ( 2 π / b ) , where b = 4.369 Å is the lattice spacing along the b → lattice vector [146] . We surround one of the NLs with a closed, tetragonal prism and calculate the Wilson loop around k y -normal squares as a function of the polar momentum k θ (exact coordinates provided in the Supplemental Material, Sec. B2 [97] ); we observe Z 2 -nontrivial winding [inset panel in (b)], indicating a nontrivial monopole charge [46,83,84] . (d) When SOC is introduced, a gap near E F develops at all crystal momenta (gray shaded region) with the Z 4 parity index (Table I ) of a HOTI [87–90] . We also explore the topology of the gapped regions between the MNLs by calculating the z -directed Wilson loop W 1 ( k x , k y ) [Fig. 3(a) ] over the lower N bands in the absence of SOC. In BZ planes indexed by k y away from the MNLs, W 1 ( k x , k y ) exhibits gaps at θ 1 / 2 π ≈ ± 0.25 [Figs. 3(c) and 3(d) and Supplemental Material, Sec. B3 [97] ], allowing us to calculate a nested Wilson loop matrix W 2 ( k y ) for which det [ W 2 ( k y ) ] = exp ( i γ 2 ) , where γ 2 is the nested Berry phase [1–3] . In bulk-gapped k y -indexed planes, det [ W 2 ( k y ) ] is quantized at ± 1 [Fig. 3(b) ], indicating that γ 2 = π (0) below (above) the MNL, implying that the Hamiltonians H k y ( k x , k z ) of planes in the two regions are topologically distinct. This quantization can be understood from two perspectives: the bulk and the Wilson loop. From a bulk perspective, H k y ( k x , k z ) is invariant under a local spinless time-reversal symmetry I × T ˜ that preserves the signs of k x , z and squares to + 1 . H k y ( k x , k z ) , therefore, lies in Class AI of the Altland-Zirnbauer classification [147,148] with codimension [44,149] D mod 8 = 6 , implying a Z 2 topology. This topology can be diagnosed by considering the Wilson-loop perspective. In the Supplemental Material, Sec. B3 [97] , we show that I × T ¯ , which acts on W 1 ( k x , k y ) as an antiunitary particle-hole symmetry Ξ ˜ that preserves the signs of k x , y [7,122,124] , enforces det [ W 2 ( k y ) ] = ± 1 when it is evaluated over any Ξ ˜ -symmetric grouping of Wilson bands [Figs. 3(c) and 3(d) ], including (but not limited to) the Wilson bands between θ 1 / 2 π ≈ ± 0.25 in Figs. 3(c) and 3(d) . Crucially, building on [1,2,50] , which employ a different definition of γ 2 reliant on fourfold rotation symmetry, the π shift in γ 2 indicates that H k y ( k x , k z ) above and below an MNL are equivalent to topologically distinct 2D magnetic atomic limits [50,93] (or trivialized fragile phases [11,47–51] ) that differ by the presence or absence of topological corner (hinge) modes (Supplemental Material, Sec. A [97] ). This implies that MNLSMs are higher-order topological SMs [50,119] with flat-band-like hinge states (Fig. 1 and Supplemental Material, Sec. A [97] ). Thus, the jump in γ 2 as k y passes through an MNL [Fig. 3(b) ] represents a new example of a topological “descent relation,” analogous to the jump in Berry phase as the line on which it is calculated passes through a Dirac point in 2D and an NL in 3D [44] . As in a Weyl SM [150,151] , the winding of the Wilson loop evaluated on a closed surface around the MNL [Fig. 2(b) ] captures the difference in topology between the gapped planes above and below it [46] , which, here, is the gapless point in W 1 (which is well defined when W 1 is evaluated on a slightly distorted path that avoids the MNLs). 3 10.1103/PhysRevLett.123.186401.f3 FIG. 3. (a) Bulk Wilson loops calculated for β - MoTe 2 from first principles in the absence of SOC (Supplemental Material, Sec. B3 [97] ). At values of k y away from the MNLs in Fig. 2 , there is a large gap in the z -directed Wilson loop spectrum W 1 ( k x , k y ) between θ 1 / 2 π ≈ ± 0.25 ; representative examples are shown in (c) and (d) for k y = 0 , π , respectively. (b) The determinant of the nested Wilson matrix W 2 ( k y ) calculated over the Wilson bands between θ 1 = ± π / 2 [1–3,50] is quantized at ± 1 by the antiunitary symmetry ( I × T ˜ ) 2 = + 1 (Supplemental Material, Sec. B3a [97] ), and jumps as it passes over an MNL, indicating that k y -indexed planes above and below the MNL are topologically distinct [1–3] . When SOC is introduced, β - MoTe 2 , though remaining metallic, develops a gap at all crystal momenta [Fig. 2(d) ]. Calculating the parity eigenvalues (Table I ), we find that, though the Fu-Kane Z 2 index [79,80] is trivial, the occupied bands, nevertheless, cannot be expressed as a sum of EBRs, indicating an overall nontrivial topology [93] . Specifically, while every EBR in an I -symmetric space group exhibits Kramers pairs of parity eigenvalues [93,94] for which ∑ n - ( k → ) mod 4 = 0 , the DBI in β - MoTe 2 induces an insulator with ∑ n - ( k → ) mod 4 = 2 (Table I ). Alternatively, this defines a Z 4 index [87–90] which is nontrivial. From both perspectives, β - MoTe 2 carries the parity eigenvalues of a HOTI [3,19] . Therefore, like bismuth [19,152,153] , β - MoTe 2 is a 2D TI when viewed as a quasi-2D system [139–141] , but is actually a HOTI when taken to be fully 3D. I 10.1103/PhysRevLett.123.186401.t1 TABLE I. The number of Kramers pairs with - 1 parity eigenvalues n - ( k → ) at each TRIM point in β - MoTe 2 (Supplemental Material, Sec. B1 [97] ). The Z 4 index [87–90] ∑ n - ( k → ) mod 4 = 2 . Along with the trivial weak indices [79,80,154] and the absence of fourfold and sixfold rotation symmetries [2–4,76,87,89] , this indicates that β - MoTe 2 is a HOTI. TRIM Γ X Y Z S T U R n - ( k → ) 12 14 14 14 14 14 14 14 Unlike β - MoTe 2 , β - WTe 2 , while stabilizable as a monolayer [27,141] , is unstable as a bulk crystal [155] . Nevertheless, the calculated electronic structure of artificial β - WTe 2 also exhibits DBI at the Γ point [155] , indicating that it would also be a HOTI if it could be stabilized. However, shortly, we will see that remnants of this HOTI phase are still observable in γ - WTe 2 . In Figs. 4(d) and 4(e) , we plot the (001) surface states of β - MoTe 2 calculated from first-principles (Supplemental Material, Sec. B1 [97] ). We observe large, arclike surface states around the projection of the Γ point [white arrows in Fig. 4(d) ], as well the projections of bulk states at the k y = 0 surface TRIM points. Crucially, we determine that the surface states are, in fact, gapped at all values of k x , y . This can be understood by considering the symmetry and topology consequences of the bulk DBI at Γ . In the absence of SOC, each band inversion nucleates one drumhead surface state (per spin) around the surface projection of the Γ point (Supplemental Material, Sec. A [97] ). In the absence of additional surface wallpaper group symmetries, such as mirror or glide [7,156] , these drumhead states hybridize and gap. When SOC is reintroduced, the four hybridized drumhead states (two per spin) open into two hybridized twofold surface TI cones [Figs. 4(a)–4(c) ]. Therefore, rather than being trivial Fermi arcs, the surface states of β - MoTe 2 are, in fact, the characteristic split and gapped fourfold Dirac-cone states of a HOTI [3,7,19] . Unlike the gapless surface states of 3D TIs [79,80] , the gapped Fermi arcs only appear at low energies because of the interplay of SOC and band dispersion in β - MoTe 2 and, theoretically, could be moved away from E F without changing the bulk or surface topology (Supplemental Material, Sec. A [97] ). 4 10.1103/PhysRevLett.123.186401.f4 FIG. 4. (a)–(c) Schematic surface state evolution of a HOTI driven by DBI. (a) Two bulk bands inverted at the same energy (blue dashed lines) realize a fourfold surface Dirac fermion (purple lines) [19] . (b) In the absence of specific glide reflection symmetries, this fermion is unstable [7] , and will split into two, twofold surface fermions, which may be stabilized by either a surface mirror [topological crystalline insulator (TCI)] [72,157] , glide (hourglass TCI) [7,75] , or C 2 z × T symmetry (rotation anomaly TCI) [11,76,158] . (c) In the absence of surface reflection or rotation symmetries, the twofold cones [yellow circles in (b), dashed lines in (c)] hybridize and gap, realizing the surface of a HOTI [3] . (d) Spectral weight at E F of states on the (001) surface of β - MoTe 2 plotted as a function of k x , y , and (e) along k x = 0 as a function of energy (Supplemental Material, Sec. B1 [97] ). Each of the two band inversions at the bulk Γ point [Fig. 2(d) ] nucleates a topological twofold surface cone centered at k x = k y = 0 (purple); the cones then repel each other in energy and merge with the projections of the bulk states (b). As depicted in (c), the surface bands from these cones [white arrows in (d)] hybridize and gap [yellow dashed lines in (e)] to form a narrowly avoided crossing. In γ - X Te 2 , these hybridized cones also appear as surface states [68–70,127–137] , but their gap is spanned along k x by small, topological Fermi arcs from bulk type-II Weyl points below E F [black arrows in (d)]. In the Supplemental Material, Sec. B4 [97] , we calculate the (001) surface states of γ - MoTe 2 gapped with symmetry-preserving distortion, and observe gapped HOTI surface states nearly identical to those in β - MoTe 2 [(d) and (e)]. This observation solves a longstanding mystery in γ - X Te 2 . In theoretical predictions [68–70] and bulk experimental probes [142–144] , both γ - MoTe 2 and γ - WTe 2 exhibit narrowly separated type-II Weyl points in the vicinity of doubly inverted bands. Nevertheless, as measured both directly by angle-resolved photoemission spectroscopy [127,128,130,131,136] and through quasiparticle interference in STM probes [133–135] , γ - X Te 2 crystals also exhibit huge, arclike surface states that largely overwhelm possible signatures of topological Weyl Fermi arcs. Previous works determined these large surface arcs to be topologically trivial [68–70,127–135] . However, in light of our previous analysis of similar nontrivial surface states in β - MoTe 2 , we recognize this determination to be incomplete. By explicitly calculating the surface states and bulk topology of γ - MoTe 2 when its Weyl points are gapped by slight distortion (Supplemental Material, Sec. B4 [97] ), we, instead, discover that the large Fermi arcs in γ - X Te 2 represent the split surface Dirac cones of a nonsymmetry-indicated HOTI phase that is nearby in parameter space and driven by DBI (Fig. 4 ). Given the small separation of the bulk Weyl points, which is quite sensitive to experimental conditions [135,136,159,160] , this HOTI phase may be accessible via symmetry-preserving distortion or strain, and may already be realized in existing samples. In this Letter, we have demonstrated higher-order topology in both β - and γ - X Te 2 . When SOC is neglected in β - MoTe 2 , we observe a pair of MNLs at E F . While SOC cannot be neglected in β - MoTe 2 , other centrosymmetric materials with lighter atoms and DBI are likely to also exhibit MNLs [46] . When the effects of SOC are incorporated, β - MoTe 2 develops a direct gap with the parity eigenvalues of a HOTI. Though β - MoTe 2 is, in fact, metallic, and thus, does not host the projected hinge gap required to observe its characteristic helical hinge modes [19] , it is possible that a TMD with more favorable band dispersion could be engineered by intercalation or chemical substitution. Finally, we observe that both γ - MoTe 2 and γ - WTe 2 exhibit the same large topological surface arcs as β - MoTe 2 , resolving an outstanding puzzle in TMDs, and presenting a new venue for investigating nonsymmetry-indicated higher-order topology. The authors thank Barry Bradlyn, Jennifer Cano, Yichen Hu, Fan Zhang, and Youngkuk Kim for helpful discussions. Z. W., B. J. W., and B. A. B. were supported by the Department of Energy Grant No. DE-SC0016239, the National Science Foundation EAGER Grants No. DMR 1643312 and NSF-MRSEC No. DMR-1420541, Army Research Office Grant No. ARO MURI W911NF-12-1-0461, Simons Investigator Grant No. 404513, ONR Grant No. N00014-14-1-0330, the Packard Foundation, the Schmidt Fund for Innovative Research, and a Guggenheim Fellowship from the John Simon Guggenheim Memorial Foundation. Z. W. additionally acknowledges support from the National Thousand-Young-Talents Program and the CAS Pioneer Hundred Talents Program. This work was also supported by the National Natural Science Foundation of China (Grants No. 11504117 and No. 11774317). B. Y. acknowledges support from the Willner Family Leadership Institute for the Weizmann Institute of Science, the Benoziyo Endowment Fund for the Advancement of Science, the Ruth and Herman Albert Scholars Program for New Scientists, and the European Research Council (ERC) under the European Union Horizon 2020 Research and Innovation Programme (Grant No. 815869). Z. W. and B. J. W. contributed equally to this work. Publisher Copyright: © 2019 American Physical Society.

PY - 2019/10/28

Y1 - 2019/10/28

N2 - In recent years, transition metal dichalcogenides (TMDs) have garnered great interest as topological materials. In particular, monolayers of centrosymmetric β-phase TMDs have been identified as 2D topological insulators (TIs), and bulk crystals of noncentrosymmetric γ-phase MoTe2 and WTe2 have been identified as type-II Weyl semimetals. However, angle-resolved photoemission spectroscopy and STM probes of these semimetals have revealed huge, arclike surface states that overwhelm, and are sometimes mistaken for, the much smaller topological surface Fermi arcs of bulk type-II Weyl points. In this Letter, we calculate the bulk and surface electronic structure of both β- and γ-MoTe2. We find that β-MoTe2 is, in fact, a Z4-nontrivial higher-order TI (HOTI) driven by double band inversion and exhibits the same surface features as γ-MoTe2 and γ-WTe2. We discover that these surface states are not topologically trivial, as previously characterized by the research that differentiated them from the Weyl Fermi arcs but, rather, are the characteristic split and gapped fourfold Dirac surface states of a HOTI. In β-MoTe2, this indicates that it would exhibit helical pairs of hinge states if it were bulk insulating, and in γ-MoTe2 and γ-WTe2, these surface states represent vestiges of HOTI phases without inversion symmetry that are nearby in parameter space. Using nested Wilson loops and first-principles calculations, we explicitly demonstrate that, when the Weyl points in γ-MoTe2 are annihilated, which may be accomplished by symmetry-preserving strain or lattice distortion, γ-MoTe2 becomes a nonsymmetry-indicated, noncentrosymmetric HOTI. We also show that, when the effects of spin-orbit coupling are neglected, β-MoTe2 is a nodal-line semimetal with Z2-nontrivial monopole nodal lines (MNLSM). This finding confirms that MNLSMs driven by double band inversion are the weak-spin-orbit coupling limit of HOTIs, implying that MNLSMs are higher-order topological semimetals with flat-band-like hinge states, which we find to originate from the corner modes of 2D "fragile" TIs.

AB - In recent years, transition metal dichalcogenides (TMDs) have garnered great interest as topological materials. In particular, monolayers of centrosymmetric β-phase TMDs have been identified as 2D topological insulators (TIs), and bulk crystals of noncentrosymmetric γ-phase MoTe2 and WTe2 have been identified as type-II Weyl semimetals. However, angle-resolved photoemission spectroscopy and STM probes of these semimetals have revealed huge, arclike surface states that overwhelm, and are sometimes mistaken for, the much smaller topological surface Fermi arcs of bulk type-II Weyl points. In this Letter, we calculate the bulk and surface electronic structure of both β- and γ-MoTe2. We find that β-MoTe2 is, in fact, a Z4-nontrivial higher-order TI (HOTI) driven by double band inversion and exhibits the same surface features as γ-MoTe2 and γ-WTe2. We discover that these surface states are not topologically trivial, as previously characterized by the research that differentiated them from the Weyl Fermi arcs but, rather, are the characteristic split and gapped fourfold Dirac surface states of a HOTI. In β-MoTe2, this indicates that it would exhibit helical pairs of hinge states if it were bulk insulating, and in γ-MoTe2 and γ-WTe2, these surface states represent vestiges of HOTI phases without inversion symmetry that are nearby in parameter space. Using nested Wilson loops and first-principles calculations, we explicitly demonstrate that, when the Weyl points in γ-MoTe2 are annihilated, which may be accomplished by symmetry-preserving strain or lattice distortion, γ-MoTe2 becomes a nonsymmetry-indicated, noncentrosymmetric HOTI. We also show that, when the effects of spin-orbit coupling are neglected, β-MoTe2 is a nodal-line semimetal with Z2-nontrivial monopole nodal lines (MNLSM). This finding confirms that MNLSMs driven by double band inversion are the weak-spin-orbit coupling limit of HOTIs, implying that MNLSMs are higher-order topological semimetals with flat-band-like hinge states, which we find to originate from the corner modes of 2D "fragile" TIs.

UR - http://www.scopus.com/inward/record.url?scp=85074882402&partnerID=8YFLogxK

UR - http://www.scopus.com/inward/citedby.url?scp=85074882402&partnerID=8YFLogxK

U2 - 10.1103/PhysRevLett.123.186401

DO - 10.1103/PhysRevLett.123.186401

M3 - Article

C2 - 31763917

AN - SCOPUS:85074882402

VL - 123

JO - Physical Review Letters

JF - Physical Review Letters

SN - 0031-9007

IS - 18

M1 - 186401

ER -