Topological Type-II Nodal Line Semimetal and ... - ACS Publications

Sep 19, 2017 - the Mg and Bi atoms form alternating layers of triangular sheets. Three Mg atoms situate at 1a (0, .... Interestingly, the topological ...
4 downloads 12 Views 6MB Size
Subscriber access provided by University of Sussex Library

Letter

Topological Type-II Nodal Line Semimetal and Dirac Semimetal State in Stable Kagome Compound MgBi 3

2

Xiaoming Zhang, Lei Jin, Xuefang Dai, and Guodong Liu J. Phys. Chem. Lett., Just Accepted Manuscript • DOI: 10.1021/acs.jpclett.7b02129 • Publication Date (Web): 19 Sep 2017 Downloaded from http://pubs.acs.org on September 19, 2017

Just Accepted “Just Accepted” manuscripts have been peer-reviewed and accepted for publication. They are posted online prior to technical editing, formatting for publication and author proofing. The American Chemical Society provides “Just Accepted” as a free service to the research community to expedite the dissemination of scientific material as soon as possible after acceptance. “Just Accepted” manuscripts appear in full in PDF format accompanied by an HTML abstract. “Just Accepted” manuscripts have been fully peer reviewed, but should not be considered the official version of record. They are accessible to all readers and citable by the Digital Object Identifier (DOI®). “Just Accepted” is an optional service offered to authors. Therefore, the “Just Accepted” Web site may not include all articles that will be published in the journal. After a manuscript is technically edited and formatted, it will be removed from the “Just Accepted” Web site and published as an ASAP article. Note that technical editing may introduce minor changes to the manuscript text and/or graphics which could affect content, and all legal disclaimers and ethical guidelines that apply to the journal pertain. ACS cannot be held responsible for errors or consequences arising from the use of information contained in these “Just Accepted” manuscripts.

The Journal of Physical Chemistry Letters is published by the American Chemical Society. 1155 Sixteenth Street N.W., Washington, DC 20036 Published by American Chemical Society. Copyright © American Chemical Society. However, no copyright claim is made to original U.S. Government works, or works produced by employees of any Commonwealth realm Crown government in the course of their duties.

Page 1 of 18

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry Letters

Topological Type-II Nodal Line Semimetal and Dirac Semimetal State in Stable Kagome Compound Mg3Bi2 Xiaoming Zhang, Lei Jin, Xuefang Dai, and Guodong Liu∗ School of Materials Science and Engineering, Hebei University of Technology, Tianjin 300130, China. E-mail: [email protected]

Abstract Topological type-II nodal line semimetal (NLS) was proposed quite recently and exhibits distinct properties comparing with conventional type-I NLS. Up to date, no ambient-condition stable candidate material has been reported. Here, we propose that a stable Kagome compound Mg3 Bi2 can host a type-II nodal line state with the protection of time reversal and spatial inversion symmetries. Similar to type-I NLSs, the type-II nodal line in Mg3 Bi2 is characterized by the drumhead surface states, which has not been observed in the previous type-II NLSs. The nodal line in Mg3 Bi2 can open a minor gap and a pair of three-dimensional (3D) Dirac points occur when SOC was included. The SOC-induced gap around the nodal line is quite small and the formation of 3D Dirac points are independent with the nodal line. So, the Mg3 Bi2 compound is expected to play a good candidate to investigate the exotic properties of both type-II NLS and 3D Dirac semimetal states.

∗ To

whom correspondence should be addressed

1 ACS Paragon Plus Environment

The Journal of Physical Chemistry Letters

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

TOC figure

2 ACS Paragon Plus Environment

Page 2 of 18

Page 3 of 18

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry Letters

Topological semimetals have received broad research interests in current condensed matter physics and materials science. 1,2 These materials are characterized with nontrivial band-crossings between conduction and valence bands in momentum space. Around the band-crossings, the quasiparticles exhibit distinct behavior from the usual Schrödinger-type fermions. For instance, a Dirac semimetal features with fourfold degenerate Dirac points (DPs), around which the low-energy quasiparticles are described by massless Dirac fermions. 3–7 With broken time reversal or inversion symmetry, a DP could split into two Weyl points (WPs) that are charaterized by the pair of Weyl fermions with opposite chiralities. 8–13 In momentum space, both Dirac and Weyl cones can be tilted, by controlling the tilting term in their Hamiltonian. According to the degree of tilting, two types of nodal points (usually termed as type-I and type-II) can be classified. For a type-I nodal point, the cone is only slightly tilted so that the hole-like states and electron-like states are well separated by the energy level of the nodal point. 3,8,12,14–21 However, for a type-II nodal point, the cone is completely tipped over so that the hole-like and electron-like states can coexist in certain energy levels. 11,21–24 The type-II fermions are proposed to possess drastically different properties from the type-I fermions. For example, in a type-II Weyl semimetal, the coexistence of hole and electron states could give rise to novel properties such as direction-dependent chiral anomaly, 25 modified anomalous Hall conductivity, 26 momentum space Klein tunneling and novel quantum oscillations, 27–29 which are detectable in transport and optical experiments. Therefore, type-II fermions attract increasing interests nowadays. In practice, both type-I and type-II fermions have been predicted in several materials, 3–5,7–11,22–24 and some are evidenced by experiments. 13–17,19–21 Under certain crystalline symmetry, the crossing points between bands could form one-dimensional (1D) nodal lines in the lattice momentum space, which is drastically different from the isolate nodal points in Dirac and Weyl semimetals. Such nodal line semimetals have received considerable interests since its first proposal in carbon allotropes. 30 Currently, several nodal line materials are proposed, with exotic physical properties revealed. 31–39 By analogy with the nodal point materials, the nodal line materials could also be classified into type-I and type-II categories, depending on the tilting degree of the band dispersion. Up to now, type-I nodal line materials have been intensively

3 ACS Paragon Plus Environment

The Journal of Physical Chemistry Letters

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

studied, 31–39 while there is only one report on type-II nodal line, proposed in compound K4 P3 quite recently. 40 To be noted, previous K4 P3 seems to be not a good candidate material for further experimental investigations and practical applications. Firstly, as an intermediate product of other K-P phases, K4 P3 is highly unstable under air and moisture conditions, and even explosive when exposed to the light. 41 Secondly, K4 P3 is quite different from a rigorous nodal line semimetal, because its band-crossing factually forms a pair of nodal loops (not closed nodal lines) traversing the whole Brillouin zone. Thirdly, K4 P3 is a topological metal rather than a semimetal, so that the characteristics of the nodal lines are most hidden by nearby trivial band intertwining. For example, the drumhead-like surface states, a typical signature of nodal line materials, are not evident in K4 P3 . Therefore, there is urgent need to explore stable topological materials with clear signature of type-II nodal line state. In current work, based on first-principles calculations, we report that Mg3 Bi2 , an ambientstable Kagome compound, 42 presents a clear type-II nodal line state close to Fermi level when spin-orbit coupling (SOC) is absent. Along the nodal line, the type-II band dispersion is observed. Unlike the type-II nodal loops which traverse the whole Brillouin zone in former K4 P3 , the bandcrossing in Mg3 Bi2 factually forms a closed type-II nodal line nesting inside the Brillouin zone. Benefiting from the clean band structure around the nodal line, one can observe the drumhead surface states either inside or outside the nodal line on the (001) surface depending on surface termination. When SOC is included, the nodal line is weakly gapped, and a pair of three-dimensional (3D) DPs that are independent with the nodal line appears. Therefore, we expect the topological semimetal Mg3 Bi2 could provide a good platform to study both properties of type-II nodal line and 3D Dirac semimetal states. To investigate the topological properties of Mg3 Bi2 compound, first-principles calculations are carried out in the framework of density functional theory (DFT), 43 as formulated in the Vienna ab initio simulation package (VASP). 44 The projector augmented wave (PAW) method is adopted for the ionic potentials. Exchange-correlation potential is treated by the generalized gradient approximation (GGA) of Perdew-Burke-Ernzerhof (PBE) functional. 45 Plane waves with a kinetic energy

4 ACS Paragon Plus Environment

Page 4 of 18

Page 5 of 18

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry Letters

Figure 1: Crystal structure of Kagome compound Mg3 Bi2 in (a) top and (b) side views. (c) Bulk Brillouin zone and the projected surface Brillouin zones of the (001) and (010) surfaces with highsymmetry points labelled. cutoff of 600 eV were adopted as basis set. The Brillouin zone is sampled with 11×11×5 Γcentered k-mesh for the structural optimization and with 19×19×9 k-mesh for the self-consistent calculations. The energy and force convergence criteria are set to be 10−6 eV and 0.001 eV/Å, respectively. The electronic structure calculations are performed both without and with SOC. The topological features of surface states are studied by constructing the maximally localized Wannier functions, 46,47 making use of the WANNIER-TOOLS package by Wu, 48 with combination of an iterative Green’s function approach. 49 According to the previous experimental studies, 42 the Mg3 Bi2 compound is a semimetal phase and stable under ambient conditions. It possesses a layered Kagome lattice structure with the space group of P3m1 (No. 164). As shown in Figure 1a and 1b, the Mg and Bi atoms form alternating layers of triangular sheets. Three Mg atoms situate at 1a (0, 0, 0) and 2d (1/3, 2/3, z) Wyckoff sites, and two Bi atoms occupy 2d’ (1/3, 2/3, z’) sites, respectively. The optimized lattice constants are a=b= 4.702Å, c= 7.436Å, and well match the experimental ones (a=b= 4.666Å, c= 7.401Å) 42 with the deviation less than 1%. The values of z and z’ from lattice relaxation are 0.6286 and 0.2212, respectively. The optimized lattice structure is adopted for the following calculations. We start with the electronic band structure of Mg3 Bi2 compound in the absence of SOC. As shown in Figure 2a, we find that the material manifests a semimetal band structure: conduction and valence bands overlap along the Γ-M and K-Γ directions, and are well separated around the 5 ACS Paragon Plus Environment

The Journal of Physical Chemistry Letters

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Figure 2: (a) Calculated electronic band structure of Mg3 Bi2 compound without SOC. (b) Enlarged band structure along K-Γ and Γ-M paths. (c) Illustration of the type-II nodal line in the bulk Brillouin zone. Noted, the nodal line dose not exactly situate at the kz =0 plane. (d) The Fermi surface of Mg3 Bi2 compound, where electron and hole pockets are denoted as red and blue colours, respectively. (e) The top view of the Fermi surface.

6 ACS Paragon Plus Environment

Page 6 of 18

Page 7 of 18

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry Letters

Fermi level when they are away from these directions. On considering the fact that GGA usually underestimates the band gap width of semimetal and semiconductor, here we check the electronic band structure by using a combination of modified Becke-Johnson exchange potential with local density approximation (MBJLDA), 50 as it predicts band gap and band order in topological materials with favorable accuracy similar to more expensive hybrid functional and GW calculations. 51–54 The results show a similar semimetal signature with the band crossing along the K-Γ direction (see Figure S1 in Supporting Information). Such semimetal character of Mg3 Bi2 has already been evidenced by former electrical conductivity and soft x-ray absorption experiments. 55,56 Because of the same topological nature by the two computations and sufficient experimental evidences, we use electronic structure from GGA for the discussions in the following. By a careful investigation to the band overlaps in Figure 2b, we find that there exist a linear band-crossing along the Γ-K direction and a tiny band gap along the Γ-M direction. Considering the slope of band dispersion, the linear crossing point along Γ-K direction is factually a type-II DP. Since the system possesses both time reversal (T) and spatial inversion (P) symmetries, under which the spinless Hamiltonian is always real valued, 30 this type-II DP cannot be isolated but belongs to a nodal line. After a careful scan of electronic band structure, we find the crossing between the conduction and valence bands indeed produces a closed type-II nodal line around the Γ point, as schematically shown in Figure 2c. To be noted, except for the crossing points along Γ-K direction, the obtained nodal line in Mg3 Bi2 does not exactly lie on the mirror plane of kz =0, but slightly wrinkles in the 3D Brillouin zone. Therefore, the nodal line does not arise from the mirror symmetry, but is merely protected by the coexistence of T and P symmetries. The shape of the wrinkled nodal line can be easily traced from the isoenergy contour by setting an energy level slightly above the nodal line. 38,57 A contour at the Fermi level and the obtained Fermi surface are shown in Figure 2d and 2e for Mg3 Bi2 . The Fermi surface clearly consists of two parts: six ellipsoid-shaped hole pockets and one bagel-shaped electron pocket centered at the Γ point. It can be observed that the bagel-shaped Fermi surface is slightly twisted along kz direction, which reflects the wrinkled nodal line centering the Γ point. It is also worth noting that, Mg3 Bi2 factually possesses a closed type-II

7 ACS Paragon Plus Environment

The Journal of Physical Chemistry Letters

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

nodal line, drastically different from that in previous K4 P3 , where band-crossing forms a pair of unclosed type-II nodal loops traversing through the whole Brillouin zone. 40 Topological nodal line materials usually manifest drumhead surface states either outside or inside the projected nodal line. 30–39 However, such important signature is not evident in previous type-II nodal line material K4 P3 , 40 mainly because it inherently possesses an metallic electronic structure so that the nature of nodal line is mostly hidden by nearby trivial band intertwining. Now we turn to the surface states of Mg3 Bi2 compound. For the Mg-terminated (001) surface, it clearly manifests drumhead surface states nestling inside the projected nodal line, as shown in Figure 3a. Its surface band has small dispersion with a bandwidth of 50 meV, being comparable with that of typical type-I nodal line in carbon allotropes. 30 Interestingly, the topological surface character of Bi-terminated (001) surface is drastically different, where the drumhead surface bands situate outside the projected nodal line, as displayed in Figure 3b. Such dependence between topological surface states and surface terminations is previously described to originate from the unique topological properties of 1D atomic chains under PT symmetry. 57 The above results clearly indicate that the type-II nodal-line materials also feature with drumhead surface character, similar to the previous type-I category. In general, SOC can shift the topological properties in materials. For instance, SOC usually drives the nodal line materials into the other topological phases, such as topological insulators, 58 Dirac and Weyl semimetals, 34,59 and also other kinds of nodal line materials. 60,61 For Mg3 Bi2 compound, we display the calculated electronic band structure under SOC in Figure 4a and 4b. It can be seen that a small indirect gap occurs around the type-II nodal line under SOC. Simultaneously a pair of type-I 3D DPs along the Γ-A direction is generated as shown in Figure 4b and 4c, which indicates Mg3 Bi2 compound is a 3D Dirac semimetal. The existence of the 3D DPs under SOC is also confirmed by the MBJLDA computation (see Figure S2 in Supporting Information). The properties of 3D Dirac semimetals have been mostly studied in the well known Na3 Bi 3 and Cd3 As2 materials. 5 However, on considering the facts that Na3 Bi is not stable in air condition and Cd3 As2 naturally contains toxic arsenic element, we argue that current Mg3 Bi2 could serve as a

8 ACS Paragon Plus Environment

Page 8 of 18

Page 9 of 18

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry Letters

Figure 3: (a) Calculated (001) surface band structure without SOC for (a) Mg-terminated and (b) Bi-terminated surface. Drumhead surface states are denoted by arrows.

9 ACS Paragon Plus Environment

The Journal of Physical Chemistry Letters

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Figure 4: (a) Calculated electronic band structure of compound Mg3 Bi2 with SOC. (b) Enlarged band structure along and K-Γ and Γ-A paths. (c) The position of the pair of 3D DPs in the bulk Brillouin zone. good 3D Dirac material for further experimental investigations and practical applications. Next, we turn to the surface states of Mg3 Bi2 compound under SOC. When projected on (001) surface, the whole Γ-A path superposes into Γ point on the surface (see Figure 1c). In this case, the projected bulk Dirac cone that situates at the Γ point is visible (see Figure 5a), while we cannot draw the configuration of its nontrivial surface states because the Fermi surface at the Dirac cone (at E= 0.251eV) is just a point (see Figure 5b). Then we consider (110) surface projection, where e Ze path in the side surface Brillouin zone (see Figure 1c). The the bulk Γ-A path projects to Γprojected surface band structure and Fermi surface at E= 0.251eV are shown in Figure 5c and 5d, where the exact positions of DPs are indicated. Although the projection of Dirac cone is hidden by the other bulk states (see Figure 5c), the profile of two half-circle-shaped Fermi arcs could be faintly visible in the given Fermi surface (see Figure 5d). Before closing, we would like to make further discussions on the topological properties of Mg3 Bi2 compound. Because of the inherent type-II band dispersion, SOC opens an indirect gap

10 ACS Paragon Plus Environment

Page 10 of 18

Page 11 of 18

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry Letters

Figure 5: (a) Calculated (001) surface band structure with SOC and (b) corresponding Fermi surface at the chemical potential of 0.251 eV. (c) and (d) are similar with (a) and (b), but for the case of side (110) surface. The position of bulk DP is denoted by the cyan points. In (d), the arrows point the profile of two half-circle-shaped Fermi arcs connecting the DPs. around the nodal line, which is different from the direct gap in previous type-I nodal line materials. Importantly, we note that the SOC-induced gap is only 36 meV along Γ-K direction (see Figure 4b), quite smaller than that of most typical nodal line materials such as antiperovskite Cu3 NPd (>60meV), 33,34 hexagonal pnictide CaAgBi ( 80meV), 51 and BaSn2 category (50-160meV). 50 Benefiting from the small SOC-induced gap, the type-II nodal line semimetal state and its unique properties are promising to be observed in Mg3 Bi2 compound. The SOC strength may be further diminished by doping lighter elements such as Sb or As, potentially making the type-II nodal-line state more prominent. Tuning the SOC strength by alloying with lighter elements has already been proved effectively in former topological Weyl and Dirac semimetals . 62,63 The pair of 3D DPs along the Γ-A direction is found to be independent with the existence of type-II nodal line, therefore both type-II nodal line and 3D Dirac states are expected to be achievable in current Mg3 Bi2 system. In summary, on the basis of first-principles calculations, we predict the existence of a type-II nodal line state in a stable Kagome compound Mg3 Bi2 when SOC is absent. The nodal line is

11 ACS Paragon Plus Environment

The Journal of Physical Chemistry Letters

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 12 of 18

formed by continuous type-II nodal points and is protected by the time reversal and spatial inversion symmetries. Different from previous unclosed type-II nodal loops in the unstable compound K4 P3 , we firstly report a closed type-II nodal line in ambient-condition stable Mg3 Bi2 material. Depending on surface termination, the drumhead surface states either inside or outside the nodal line were observed on the (001) surface projection of Mg3 Bi2 . The nodal line is weakly gapped and a pair of 3D DPs emerges when SOC is included. So, Mg3 Bi2 compound is also a 3D Dirac semimetal. On considering the facts that the SOC-induced gap around the nodal line is much smaller than typical nodal line materials and the formation of 3D Dirac points is independent with the nodal line, the Mg3 Bi2 is promising to serve as a good candidate material to study both properties of type-II nodal line semimetal and 3D Dirac semimetal states.

Acknowledgement This work is supported by the Natural Science Foundation of Hebei Province (No. E2016202383), the Program for Leading Talents in Science and Technology Innovation of Chongqing City (No. cstckjcxljrc19), and the Basic and Frontier Research Project of Chongqing City (No. cstc2014jcyja50005)

Supporting Information Available Electronic band structures with modified Becke-Johnson potential for Mg3 Bi2 ; electronic band structures for Mg3 Sb2 and Mg3 As2 ; some discussions and references.

References (1) Burkov A. A. Topological Semimetals. Nat. Mater. 2016, 15, 1145-1148. (2) Yan, B.; Felser, C. Topological Materials: Weyl Semimetals. Annu. Rev. Condens. Matter Phys. 2017, 8, 337-354.

12 ACS Paragon Plus Environment

Page 13 of 18

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry Letters

(3) Wang, Z.; Sun, Y.; Chen, X.-Q.; Franchini, C.; Xu, G.; Weng, H.; Dai, X.; Fang, Z. Dirac Semimetal and Topological Phase Transitions in A3 Bi (A=Na, K, Rb). Phys. Rev. B 2012, 85, 195320. (4) Young, S. M.; Zaheer, S.; Teo, J. C. Y.; Kane, C. L.; Mele, E. J.; Rappe, A. M. Dirac Semimetal in Three Dimensions. Phys. Rev. Lett. 2012, 108, 140405. (5) Wang, Z.; Weng, H.; Wu, Q.; Dai, X.; Fang, Z. Three-Dimensional Dirac Semimetal and Quantum Transport in Cd3 As2 . Phys. Rev. B 2013, 88, 125427. (6) Yang, B.-J.; Nagaosa, N. Classification of Stable Three-dimensional Dirac Semimetals with Nontrivial Topology. Nat. Commun. 2014, 5, 4898. (7) Sheng, X.-L.; Wang, Z.; Yu, R.; Weng, H.; Fang, Z.; Dai X. Topological Insulator to Dirac Semimetal Transition Driven by Sign Change of Spin-orbit Coupling in Thallium Nitride. Phys. Rev. B 2014, 90, 245308. (8) Wan, X.; Turner, A. M.; Vishwanath, A.; Savrasov, S. Y. Topological Semimetal and Fermi-arc Surface States in the Electronic Structure of Pyrochlore Iridates. Phys. Rev. B 2011, 83, 205101. (9) Murakami, S. Phase Transition Between the Quantum Spin Hall and Insulator Phases in 3D: Emergence of a Topological Gapless Phase. New J. Phys. 2007, 9, 356. (10) Sheng, X.; Yu, Z.; Yu, R.; Weng, H.; Yang, S. A. d Orbital Topological Insulator and Semimetal in the Antifluorite Cu2 S Family: Contrasting Spin Helicities, Nodal Box, and Hybrid Surface States. J. Phys. Chem. Lett. 2017, 8, 3506-3511. (11) Chang, G.; Xu, S.-Y.; Sanchez, D. S.; Huang, S.-M.; Lee, C.-C.; Chang, T.-R.; Bian, G.; Zheng, H.; Belopolski, I.; Alidoust, N.; et al. A Strongly Robust Type II Weyl Fermion Semimetal State in Ta3 S2 . Sci. Adv. 2016, 2, e1600295. (12) Burkov, A. A; Balents, L. Weyl Semimetal in a Topological Insulator Multilayer. Phys. Rev. Lett. 2011, 107, 127205. 13 ACS Paragon Plus Environment

The Journal of Physical Chemistry Letters

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

(13) Nielson, H. B.; Ninomiya, M. The Adler-Bell-Jackiw Anomaly and Weyl Fermions in a Crystal. Phys. Lett. B 1983, 130, 389-396. (14) Liu, Z. K.; Zhou, B.; Zhang, Y.; Wang, Z. J.; Weng, H. M.; Prabhakaran, D.; Mo, S.-K.; Shen, Z. X.; Fang, Z.; Dai, X.; et al. Discovery of a Three-Dimensional Topological Dirac Semimetal, Na3 Bi. Science 2014, 343, 864-867. (15) Xu, S.-Y.; Liu, C.; Kushwaha, S. K.; Sankar, R.; Krizan, J. W.; Belopolski, I.; Neupane, M.; Bian, G.; Alidoust, N.; Chang, T.-R.; et al. Observation of Fermi Arc Surface States in a Topological Metal. Science 2015, 347, 294-298. (16) Neupane, M.; Xu, S.-Y.; Sankar, R.; Alidoust, N.; Bian, G.; Liu, C.; Belopolski, I.; Chang, T.R.; Jeng, H.-T.; Lin, H.; et al. Observation of a Three-Dimensional Topological Dirac Semimetal Phase in High-Mobility Cd3 As2 . Nat. Commun. 2014, 5, 3786. (17) Yang B.-J.; Nagaosa, N. Classification of Stable Three-Dimensional Dirac Semimetals with Nontrivial Topology. Nat. Commun. 2014, 5, 4898. (18) Weng, H.; Fang, C.; Fang, Z.; Bernevig, B. A.; Dai, X. Weyl Semimetal Phase in Noncentrosymmetric Transition-Metal Monophosphides. Phys. Rev. X 2015, 5, 011029. (19) Huang, S.-M.; Xu, S.-Y.; Belopolski, I.; Lee, C.-C.; Chang, G.; Wang, B.; Alidoust, N.; Bian, G.; Neupane, M.; Zhang, C.; et al. A Weyl Fermion Semimetal with Surface Fermi Arcs in the Transition Metal Monopnictide TaAs Class. Nat. Commun. 2015, 6, 7373. (20) Xu, S.-Y.; Belopolski, I.; Alidoust, N.; Neupane, M.; Bian, G.; Zhang, C.; Sankar, R.; Chang, G.; Yuan, Z.; Lee, C.-C.; et al. Discovery of a Weyl Fermion Semimetal and Topological Fermi arcs. Science 2015, 349, 613-617. (21) Lu, L.; Wang, Z.; Ye, D.; Ran, L.; Fu, L.; Joannopoulos, J. D.; Soljai, M. Experimental Observation of Weyl Points. Science 2015, 349, 622-624.

14 ACS Paragon Plus Environment

Page 14 of 18

Page 15 of 18

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry Letters

(22) Chang, T.-R.; Xu, S.-Y.; Chang, G.; Lee, C.-C.; Huang, S.-M.; Wang, B.; Bian, G.; Zheng, H.; Sanchez, D. S.; Belopolski, I.; et al. Arc-Tunable Weyl Fermion Metallic State in MoxW1xTe2. Nat. Commun. 2016, 7, 10639. (23) Soluyanov, A. A.; Gresch, D.; Wang, Z.; Wu, Q.; Troyer, M.; Dai, X.; Bernevig, B. A. Type-II Weyl Semimetals. Nature 2015, 527, 495-498. (24) Chang, T.-R.; Xu, S.-Y.; Sanchez, D. S.; Tsai, W.-F.; Huang, S.-M.; Chang, G.; Hsu, C.H.; Bian, G.; Belopolski, I.; Yu, Z.-M. et al. Type-II Symmetry-Protected Topological Dirac Semimetals. Phys. Rev. Lett. 2017, 119, 026404. (25) Udagawa M.; Bergholtz, E. J. Field-Selective Anomaly and Chiral Mode Reversal in Type-II Weyl Materials. Phys. Rev. Lett. 2016, 117, 086401. (26) Zyuzin A. A.; Tiwari, R. P. Intrinsic Anomalous Hall Effect in Type-II Weyl Semimetals. JETP Lett. 2016, 103, 717-722. (27) O’Brien, T. E.; Diez, M.; Beenakker, C. W. J. Magnetic Breakdown and Klein Tunneling in a Type-II Weyl Semimetal. Phys. Rev. Lett. 2016, 116, 236401. (28) Koshino, M. Cyclotron Resonance of Figure-of-Eight Orbits in a Type-II Weyl Semimetal. Phys. Rev. B 2016, 94, 035202. (29) Yu, Z.-M.; Yao, Y.; Yang A. S. Predicted Unusual Magnetoresponse in Type-II Weyl Semimetals. Phys. Rev. Lett. 2016, 117, 077202. (30) Weng, H.; Liang, Y.; Xu, Q.; Yu, R.; Fang, Z.; Dai, X.; Kawazoe, Y. Topological Node-Line Semimetal in Three Dimensional Graphene Networks. Phys. Rev. B 2015, 92, 045108. (31) Fang, C.; Weng, H. M.; Dai, X.; Fang, Z. Topological Nodal Line Semimetals. Chin. Phys. B 2016, 25, 117106. (32) Chen, Y.; Xie, Y.; Yang, S. A.; Pan, H.; Zhang, F.; Cohen, M. L.; Zhang, S. Nanostructured Carbon Allotropes with Weyl-like Loops and Points. Nano Lett. 2015, 15, 6974-6978. 15 ACS Paragon Plus Environment

The Journal of Physical Chemistry Letters

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

(33) Kim, Y.; Wieder, B. J.; Kane, C. L.; Rappe, A. M. Dirac Line Nodes in Inversion-Symmetric Crystals. Phys. Rev. Lett. 2015, 115, 036806. (34) Yu, R.; Weng, H.; Fang, Z.; Dai, X.; Hu, X. Topological Node-Line Semimetal and Dirac Semimetal State in Antiperovskite Cu3 PdN. Phys. Rev. Lett. 2015, 115, 036807. (35) Fang, C.; Chen, Y.; Kee, H. Y.; Fu, L. Topological Nodal Line Semimetals with and without Spin-Orbital Coupling. Phys. Rev. B 2015, 92, 081201. (36) Xie, L. S.; Schoop, L. M.; Seibel, E. M.; Gibson, Q. D.; Xie, W.; Cava, R. J. A New Form of Ca3 P2 with a Ring of Dirac Nodes. APL Mater. 2015, 3, 083602. (37) Chan, Y.-H.; Chiu, C.-K.; Chou, M. Y.; Schnyder, A. P. Ca3 P2 and other Topological Semimetals with Line Nodes and Drumhead Surface States. Phys. Rev. B 2016, 93, 205132. (38) Li, R. H.; Ma, H.; Cheng, X.; Wang, S.; Li, D.; Zhang, Z.; Li, Y.; Chen, X.-Q. Dirac Node Lines in Pure Alkali Earth Metals. Phys. Rev. Lett. 2016, 117, 096401. (39) Zhang, X.; Yu, Z.-M.; Sheng, X.-L.; Yang, H. Y.; Yang, A. S. Coexistence of Four-Band Nodal Rings and Triply Degenerate Nodal Points in Centrosymmetric Metal Diborides. Phys. Rev. B 2017, 95, 235116. (40) Li, S.; Yu, Z.-M.; Liu, Y.; Guan, S.; Wang, S.-S.; Zhang, X.; Yao, Y.; Yang A. S. Type-II Nodal Loops: Theory and Material Realization. Phys. Rev. B 2017, 96, 081106(R). (41) Schnering, H. G. von; Hartweg, M.; Hartweg, U.; Hnle, W. K4 P3 , a Compound Containing the Radical Anion P. Angew. Chem. Int. Ed. Engl. 1989, 28, 56-58. (42) Watson, L. M.; Marshall, C. A. W.; Cardosot, C. P. On the Electronic Structure of the Semiconducting Compounds Mg3 Bi2 , and Mg3 Sb2 . J. Phys. F: Met. Phys. 1984, 14, 113-121. (43) Kresse, G.; Joubert, D. From Ultrasoft Pseudopotentials to the Projector Augmented-Wave Method. Phys. Rev. B 1999, 59, 1758-1775. 16 ACS Paragon Plus Environment

Page 16 of 18

Page 17 of 18

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry Letters

(44) Kresse, G.; Hafner, J. Ab Initio Molecular Dynamics for Liquid Metals. Phys. Rev. B 1993, 47, 558-561. (45) Perdew, J. P.; Burke, K.; Ernzerhof, M. Generalized Gradient Approximation Made Simple. Phys. Rev. Lett. 1996, 77, 3865-3868. (46) Marzari, N.; Vanderbilt, D. Maximally Localized Generalized Wannier Functions for Composite Energy Bands. Phys. Rev. B 1997, 56,12847. (47) Mostofi, A. A.; Yates, J. R.; Lee, Y.-S.; Souza, I.; Vanderbilt, D.; Marzari, N. Wannier90: a Tool for Obtaining Maximally-Localised Wannier Functions. Comput. Phys. Commun. 2008, 178, 685-699. (48) Wu, Q. S.; Zhang, S. N.; Song, H.-F.; Troyer, M.; Soluyanov, A. A. WannierTools: An Open-Source Software Package for Novel Topological Materials. ArXiv e-prints: 1703.07789 2017. (49) Sancho, M. P. L.; Sancho, J. M. L.; Rubio, Quick Iterative Scheme for the Calculation of Transfer Matrices: Application to Mo (100). J. J. Phys. F 1984, 14, 1205-1215. (50) Becke, A. D.; Johnson, E. R. A Simple Effective Potential for Exchange. J. Chem. Phys. 2006 124, 221101. (51) Tran, F.; Blaha, P. Accurate Band Gaps of Semiconductors and Insulators with a Semilocal Exchange-Correlation Potential. Phys. Rev. Lett. 2009, 102, 226401. (52) Feng, W.; Xiao, D.; Zhang, Y.; Yao Y. Half-Heusler Topological Insulators: a First-Principles Study with the Tran-Blaha Modified Becke-Johnson Density Functional. Phys. Rev. B 2010, 82, 235121. (53) An, W.; Wu, F.; Jiang, H.; Tian, G.-S.; Li, X.-Z. Systematic Investigation on Topological Properties of Layered GaS and GaSe Under Strain. J. Chem. Phys. 2014, 141, 084701.

17 ACS Paragon Plus Environment

The Journal of Physical Chemistry Letters

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

(54) Quan, Y.; Yin, Z. P.; Pickett W. E. Single Nodal Loop of Accidental Degeneracies in Minimal Symmetry: Triclinic CaAs3 . Phys. Rev. Lett. 2017, 118, 176402. (55) Ferrier, R.P.; Herrell, D.J. Conduction in Amorphous Mg-Bi and Mg-Sb Alloys. J. Non-Cryst. Solids 1970, 2, 278-283. (56) Slowik J. H. Core Exeitons in Amorphous Magnesium Alloys. Phys. Rev. B 1974, 10, 416431. (57) Huang,H.; Liu, J.; Vanderbilt, D.; Duan, W. Topological Nodal-Line Semimetals in AlkalineEarth Stannides, Germanides, and Silicides. Phys. Rev. B 2016, 93, 201114(R). (58) Yamakage, A.; Yamakawa, Y.; Tanaka, Y.; Okamoto, Y. Line-Node Dirac Semimetal and Topological Insulating Phase in Noncentrosymmetric Pnictides CaAgX (X = P, As). J. Phys. Soc. Jpn. 2016, 85, 013708 (59) Weng, H.; Fang, C.; Fang, Z.; Dai, X. Coexistence of Weyl Fermion and Massless Triply Degenerate Nodal Points. Phys. Rev. B 2016, 94, 165201. (60) Bian, G.; Chang, T.-R.; Sankar, R.; Xu, S.-Y.; Zheng, H.; Neupert, T.; Chiu, C.-K.; Huang, S.-M.; Chang, G.; Belopolski I.; et al. Topological Nodal-Line Fermions in Spin-Orbit Metal PbTaSe2 . Nat. Commun. 2016, 7, 10556. (61) Bian, G.; Chang, T.-R.; Zheng, H.; Velury, S.; Xu, S.-Y.; Neupert, T.; Chiu, C.-K.; Huang, S.M.; Sanchez,D. S.; Belopolski, I.; et al. Drumhead Surface States and Topological Nodal-Line Fermions in TlTaSe2 . Phys. Rev. B 2016, 93, 121113(R). (62) Wang, Z.; Vergniory, M. G.; Kushwaha, S.; Hirschberger, M.; Chulkov, E. V.; Ernst, A.; Ong, N. P.; Cava, R. J.; Bernevig, B. A. Time-Reversal-Breaking Weyl Fermions in Magnetic Heusler Alloys. Phys. Rev. Lett. 2016, 117, 236401. (63) Narayan, A.; Sante, D. D.; Picozzi, S.; Sanvito S. Topological Tuning in Three-Dimensional Dirac Semimetals. Phys. Rev. Lett. 2014, 113, 256403. 18 ACS Paragon Plus Environment

Page 18 of 18