## Symmetry-enforced band crossings in trigonal materials: Accordion states and Weyl nodal lines

Y.-H. Chan, 1, 2, ∗ Berkay Kilic, 3 Moritz M. Hirschmann, 3 Ching-Kai Chiu, Leslie M. Schoop, 5 Darshan G. Joshi, 3 and Andreas P. Schnyder 3, †

4

1 Institute of Atomic and Molecular Sciences, Academia Sinica, Taipei 10617, Taiwan 2 Materials Sciences Division, Lawrence Berkeley National Laboratory, Berkeley, CA 94720, USA 3 Max-Planck-Institut f¨r Festk¨rperforschung, Heisenbergstrasse 1, D-70569 Stuttgart, Germany u o 4 Kavli Institute for Theoretical Sciences, University of Chinese Academy of Sciences, Beijing 100190, China 5 Department of Chemistry, Princeton University, Princeton, NJ 08544, USA

(ΩDated: August 5, 2019)

Nonsymmoprhic symmetries, such as screw rotations or glide reflections, can enforce band crossings within high-symmetry lines or planes of the Brillouin zone. When these band degeneracies are close to the Fermi energy, they can give rise to a number of unusual phenomena: e.g., anomalous magnetoelectric responses, transverse Hall currents, and exotic surface states. In this paper, we present a comprehensive classification of such nonsymmorphic band crossings in trigonal materials with strong spin-orbit coupling. We find that in trigonal systems there are two different types of nonsymmorphic band degeneracies: (i) Weyl points protected by screw rotations with an accordionlike dispersion, and (ii) Weyl nodal lines protected by glide reflections. We report a number of existing materials, where these band crossings are realized near the Fermi energy. This includes Cu SrSnS 2 4 and elemental tellurium (Te), which exhibit accordion Weyl points; and the telluriumsilicon clathrate Te 16 Si 38 , which shows Weyl nodal lines. The ab-initio band structures and surface states of these materials are studied in detail, and implications for experiments are briefly discussed.

## I. INTRODUCTION

The famous non-crossing theorem by Wigner and von Neumann [1] states that Bloch states with the same symmetry cannot be degenerate at a generic point in the Brillouin zone (BZ), which prevents the formation of band crossings. Instead, when two Bloch bands of the same symmetry approach each other at a generic momentum, they start to hybridize and undergo an avoided level crossing. However, the non-crossing theorem does not apply to bands with non-trivial wavefunction topology, which can form topologically protected band degeneracies [2-6]. When these band crossings are in the vicinity of the Fermi energy, they lead to a range of interesting phenomena, such as arc and drumhead surface states [713], transverse topological currents [4, 5, 14], and anomalous magnetoelectric responses [15], which could potentially be utilized for device applications [16, 17].

other hand, arise in the presence of nonsymmorphic symmetries and are globally stable, i.e., they cannot be removed even by large deformations of the Hamiltonian. In other words, these band crossings are required to exist due to nonsymmorphic symmetries alone, independent of the chemical composition and other details of the material. This fact allows us to construct the following strategy to discover new topological semimetals, which consists of three steps: (i) identify the space groups (SGs) whose nonsymmorphic symmetries enforce the desired band crossings, (ii) perform a database search for materials in these SGs, and (iii) compute the electronic band structure of these materials to check whether the band crossings are near the Fermi energy. Previously, we applied this strategy to discover new topological semimetals with hexagonal symmetries [38]. Here, we extend this analysis to trigonal systems.

There are two different types of topological band crossings, namely, accidental band crossings and symmetryenforced band crossings. The former are protected by symmorphic crystal symmetries and are only perturbatively stable [18-20]. That is, they can be adiabatically removed by large symmetry-preserving deformations of the Hamiltonian, for example, through pair annihilation. Examples of accidental band crossings include Dirac points and Dirac lines that are protected by spacetime inversion, reflection, or rotation symmetry [21-26]. Symmetry-enforced band crossings [15, 27-38], on the

∗ yanghao@umich.edu

† a.schnyder@fkf.mpg.de

There has been a number of previous works, which considered specific types of nonsymmorphic band crossings [15, 33-35] and proposed candidate materials, where these are realized [39-42]. However, a comprehensive classification of all possible nonsymmorphic bandcrossings in trigonal materials has not been performed before. We find that among the 25 trigonal SGs, six support Weyl points with accordion-like dispersions and three support Weyl nodal lines (see Table I). We identify several materials, where these topological band crossings are realized. This includes elemental tellurium (Te) with accordion Weyl points, and the trigonal high-temperature form of Te 16 Si 38 with Weyl nodal lines. As a byproduct of our analysis, we also obtain the tight filling constraints for the existence of band insulators, which are in agreement with the recent literature [43, 44]. Furthermore, we find that in trigonal materials, as opposed to hexagonal

TABLE I. Classification of nonsymmorphic band crossings in trigonal materials with time-reversal symmetry and strong spinorbit coupling. The first column lists the trigonal space groups that exhibit topological band crossings. The second and third columns indicate the high-symmetry lines and planes in which Weyl nodal points and Weyl nodal lines appear, respectively. The definition of the coordinate system and high-symmetry labels in the Brillouin zones are given in Fig. 1(a) and Fig. 3. The electron fillings ν for which a band insulator is possible, are given in the fourth column, where m N denotes the set { m, 2 m, 3 m, · · · } . Example materials which realize the predicted band crossings are listed in the last column.

| Space group     | Weyl nodal points   | Weyl nodal lines   | Filling constraint   | Materials    |
|-----------------|---------------------|--------------------|----------------------|--------------|
| 144 ( P 3 1 )   | Γ∆A (6)             | -                  | 6 N                  | Cu 2 SrSnS 4 |
| 145 ( P 3 2 )   | Γ∆A (6)             | -                  | 6 N                  | Cu 2 SrGeS 4 |
| 151 ( P 3 1 12) | Γ∆A (6)             | -                  | 6 N                  | Ag 2 HPO 4   |
| 152 ( P 3 1 21) | Γ∆A (6)             | -                  | 6 N                  | Te, Se       |
| 153 ( P 3 2 12) | Γ∆A (6)             | -                  | 6 N                  | -            |
| 154 ( P 3 2 21) | Γ∆A (6)             | -                  | 6 N                  | -            |
| 158 ( P 3 c 1)  | -                   | k x k z -plane     | 4 N                  | -            |
| 159 ( P 31 c )  | -                   | k x = π            | 4 N                  | -            |
| 161 ( R 3 c )   | -                   | k x = - k y        | 4 N                  | Te 16 Si 38  |

materials [38], there exist no band crossings protected by off-centered symmetries [35-37]. We note that our analysis takes into account the full connectivity of the bands in the entire BZ, which goes beyond the results obtained through symmetry indicators [45-47].

## A. Symmetry eigenvalues

The remainder of this paper is organized as follows. In Sec. II we classify Weyl points protected by screw rotation symmetries and discuss how symmetries constrain the multiplicities of the Weyl points. Section III is devoted to the study of Weyl nodal lines protected by glide reflection. In Sec. IV we discuss the filling constraints imposed by the nonsymmorphic symmetries. In Sec. V we present the example materials that realize the predicted band crossings. Our summary and conclusions are given in Sec. VI. In Appendix A we present additional figures of the band structure calculations. The calculated surface states of tellurium are shown in Appendix B.

## II. NON-SYMMORPHIC NODAL POINTS

In this section we show how three-fold screw-rotation symmetries in trigonal space groups lead to protected Weyl points. To derive this we employ two different approaches: (i) Using symmetry analysis we show how the evolution of eigenvalues of the symmetry operator along a line in the BZ enforces band crossings (Sec. II A). (ii) Using crystallographic group theory we show how the compatibility relations between irreducible representations (irreps) at different high-symmetry points of the BZ [48-50] lead to enforced band crossings (Sec. II B). We also give a brief discussion on how time-reversal and SG symmetries lead to higher multiplicities of the Weyl points (Sec. II C).

In general, a symmetry operator G = { g glyph[vector] | t } of a given space group is a combination of a point group symmetry g and a lattice translation glyph[vector] t . In case of a nonsymmorphic symmetry the lattice translation glyph[vector] t is some fraction of the Bravais-lattice vector. Here, we are interested in the three-fold screw rotation, which has the following form:

$$\text{train} \\ \text{s de} \\ & \text{glide} \\ \text{s im} \\ V \text{ we} \quad \text{emh that } n \, \leq \, 1 \, 9 \, \mu \, \text{ } \text{Within the } a \, \text{invariant anacos of}$$

such that p ∈ { 1 2 , } . Within the g -invariant space of the BZ, i.e. where g glyph[vector] k = glyph[vector] k , the Bloch states | ψ m ( glyph[vector] k ) 〉 can be constructed to be simultaneous eigenfunctions of G and the Hamiltonian, satisfying the following eigenvalue equation:

$$C _ { 3, p } | \psi _ { m } ( \vec { k } ) \rangle = e ^ { i \pi ( 2 m + 1 ) / 3 } e ^ { - i p k _ { z } / 3 } | \psi _ { m } ( \vec { k } ) \rangle \,, \quad ( 2. 2 )$$

with m ∈ { 0 1 2 , , } . The above equation follows from the fact that C 3 3 ,p = -p glyph[vector] R , where glyph[vector] R is a Bravais lattice vector (here glyph[vector] R = ˆ) and the minus sign arises due to z 2 π rotation of the electron spin, relevant for spin-orbit coupling. From Eq. (2.2), we observe that the eigenvalues of C 3 ,p evolve along the line segment glyph[vector] k = (0 0 , , k z ) with k z ∈ [0 , π ], i.e., along the path Γ -∆ -A in Fig. 1, which is invariant under C 3 rotation.

Note that the points Γ and A are time-reversal invariant momenta, where in the presence of time-reversal symmetry bands pair-up with their respective Kramers partner. In terms of the eigenvalues in Eq. (2.2) this means that a band with real eigenvalue pairs with itself, while a band with complex eigenvalue pairs with the band whose eigenvalue is the complex conjugate. Using the label m to identify the bands and corresponding eigenvalues, we find the following pairing scenario: At the Γ point ( k z = 0) there are two Kramers pairs of bands, namely, (1 1) and (0 2), irrespective of the value of , , p .

TABLE II. This table lists how the Bloch bands | ψ m ( glyph[vector] k ) 〉 pair up into Kramers partners at the A point of the BZ in materials with a three-fold screw rotation C 3 ,p . The pairing depends on p , i.e., the fractional lattice translation part of C 3 ,p , see Eq. (2.2). The right column indicates the eigenvalue labels m of the Kramers pairs at the A point.

|   p | Pairings of m   |
|-----|-----------------|
|   1 | (0,0), (1,2)    |
|   2 | (0,1), (2,2)    |

At the A point ( k z = π ), on the other hand, Kramers pairs depend on the value of p and pair up as indicated in Table II. Since the bands pair up differently at Γ and A, they show, in the absence of additional symmetries, a nontrivial connectivity along the path Γ -∆ -A , as displayed in Fig. 1(b). We find that for both p = 1 and p = 2 there are six bands forming a connected group, which must cross at least two times, leading to two Weyl point degeneracies with an accordion-like dispersion.

We note that the presence of additional symmetries such as inversion or mirror, which either commute or anticommute with C 3 ,p , may lead to additional degeneracies at the high-symmetry points in the BZ. This can remove the above discussed Weyl points. In particular, out of the eighteen trigonal space groups with hexagonal lattice system (P-trigonal) there are only six space groups without inversion symmetry (SGs with Nos. 144, 145, 151, 152, 153, and 154). These are the only SGs which host symmetry-enforced Weyl points with accordion-like dispersion, see Table I.

In passing we note that the seven trigonal space groups with rhombohedral lattice system (R-trigonal) do not have a nonsymmorphic three-fold screw-rotation symmetry, but just three-fold symmorphic rotation symmetries (i.e. p = 0 in C 3 ,p ). Consequently, the eigenvalues of these symmorphic rotation symmetries are momentum independent and hence there are no symmetry-enforced Weyl points for R-trigonal space groups.

## B. Compatibility relations between irreps

Next, we explain how the existence of symmetryenforced Weyl points can be derived form the compatibility relations of irreducible representations (irreps) at different high-symmetry points/lines of the BZ [4850]. Here, we give a brief outline of this derivation for SG No. 144 (P3 ). 1 Additional details can be found in Ref. [38].

Bloch bands | ψ m ( glyph[vector] k ) 〉 restricted to a high-symmetry line (or point) can be labeled by the irreducible representations (irreps) of the symmetry operators that leave this line (or point) invariant. An irrep of a given symmetry group can be viewed as a set of matrices which form a group that is homomorphic to the symmetry group, with the property that there does not exist a unitary

TABLE III. Double-valued irreps for SG No. 144 (P3 1 ) without time-reversal symmetry at the points Γ and A and at the line ∆. For the labelling of the irrpes we use the same convention as in Ref. [50].

| Irreps/Sym.       | Element   | C 3                                                |
|-------------------|-----------|----------------------------------------------------|
| ¯ Γ 4 ¯ Γ 5 ¯ Γ 6 |           | - 1 e - iπ/ 3 e iπ/ 3                              |
| ¯ ∆ 4 ¯ ∆ 5 ¯ ∆ 6 |           | e ik z / 3 e i ( k z - π ) / 3 e i ( k z + π ) / 3 |
| ¯ A 4 ¯ A 5 ¯ A 6 |           | e - i 2 π/ 3 +1 e i 2 π/ 3                         |

transformation that brings all these matrices into block diagonal form. While these sets of matrices are not unique, their trace (i.e., their character) is unchanged under any unitary transformation. Therefore, for a given high-symmetry line (or point) in the BZ one can define a character table, which lists the irrpes and their characters for all the symmetry operations that leave this line (or point) invariant.

Now, as we move from a high-symmetry point to a high-symmetry line in the BZ, the symmetry is lowered, and hence, an irrep A σ of the high-symmetry point is, in general, no longer irreducible for the high-symmetry line. That is the irrep A σ of the high-symmetry point decomposes, in general, into smaller irreps B i,σ of the high-symmetry line. However, under such a decomposition the character χ of each relevant symmetry must be preserved, i.e.,

$$\chi ( A _ { \sigma } ) = \sum _ { i = 1 } ^ { n } \chi ( B _ { i, \sigma } ) \,. \quad \text{(2.3)} \\ \text{equation defines the commatifilitu relationships be}$$

The above equation defines the compatibility relations between the irreps A σ at a high-symmetry point and the irreps B i,σ at a high-symmetry line [51, 52].

We shall now use these compatibility relations to study how the irreps decompose in SG No. 144 (P3 ) for the 1 path Γ -∆ -A , which is left invariant under the screw rotation C 3 1 , . Since we are interested in electronic band structures with spin-orbit coupling, we need to consider the double-valued irreps of SG No. 144 (P3 ) (see Ta1 ble III), as the identity operator is not equivalent to a rotation by 2 π but rather by 4 π on the Bloch sphere. At the time-reversal invariant momenta (TRIMs), Bloch bands form Kramers pairs due to time-reversal symmetry. This means that at the TRIMs pairs of complex conjugate irreps combine together to form a time-reversal symmetric irrep [51, 52]. For the 1 d irreps of Table III we can construct these time-reversal symmetric irreps by inspection. At the Γ point we find that the time-reversal symmetric irreps are given by the direct product of Γ ¯ 5 with ¯ Γ 6 and ¯ Γ 4 with itself. Similarly, at the A point the

FIG. 1. (a) Bulk BZ (black lines) for trigonal space groups with hexagonal lattice systems (P-trigonal). The high-symmetry points and lines are indicated in red. The green lines show the surface BZ for surfaces perpendicular to the (-110) direction. (b), (c) Kramers pairings and band connectivity diagrams for the Γ -∆ -A line of the hexagonal BZ, which is left invariant by the three-fold screw rotation C 3 ,p . Here, we focus on the case p = 1, e.g., on SG No. 144 (P3 1 ). (Similar band crossings are also realized for p = 2.) The left panels in (b) and (c) indicate how the Bloch bands | ψ m ( glyph[vector] k ) 〉 pair up into Kramers partners at the time-reversal invariant momenta Γ and A , see Tables II and III. The right panels in (b) and (c) show the band connectivity diagrams along the path Γ -∆ -A . Six bands form a connected group with at least two symmetry-enforced Weyl points, as explained in the text. In (b) the bands are marked by their C 3 ,p eigenvalue label m , while in (c) they are marked by their double-valued irreps ∆ . ¯ i

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00901v1-figure-26e0e567-086a-4173-9ed4-73bf8bd269b6.jpg)

time-reversal symmetric irreps are given by the direct product of ¯ A 4 with ¯ A 6 and ¯ A 5 with itself, see Fig. 1(c). As we move from the TRIMs Γ and A to a point on the high-symmetry line ∆, the two-dimensional time-reversal symmetric irreps decompose into one-dimensional irreps, according to Eq. (2.3), see left part of Fig. 1(c). Connecting the respective irreps that split away from the points Γ and A we find that the bands connect in a non-trivial way, with a minimum number of two Weyl point crossings, as shown in the right part of Fig. 1(c). This is in exact agreement with the symmetry analysis presented above. Similar analysis can be done for the remaining five space groups of Table I and we find perfect agreement therein.

we assume in the following that the Weyl points at the Γ -∆ -A line have positive chirality +1.

## C. Weyl point multiplicities

- a. SG Nos. 144, 145, 151, and 153.Band structures in SG Nos. 144 and 145 possess only two symmetries: the three-fold screw rotation along the k z axis C 3 ,p and time-reversal symmetry. In order to cancel the +1 chirality of the Weyl point at the Γ -∆ -A line, we need to place a negative chirality Weyl point somewhere else in the BZ. If we put it at an arbitrary position, the screw rotation C 3 ,p would generate two more copies with the same chirality, leading to an overcompensation. Instead, in order to construct a Weyl point configuration with minimal multiplicity, we need to put the negative chirality Weyl point at one of the other two C 3 ,p -invariant lines, i.e., either on K -P -H or on K ′ -P ′ -H ′ . With this, the chiralities cancel out and we obtain a minimal Weyl point multiplicity of four, see Fig. 2(a).

So far we have only discussed the appearance of Weyl points on the Γ -∆ -A line. However, it follows from symmetry arguments together with the fermion doubling theorem [53] that these cannot be the only Weyl points in the BZ. This is because (i) there has to be another Weyl point away from the Γ -∆ -A line with opposite chirality, since the chiralities of all Weyl points formed by one pair of bands must cancel out to zero [53]; and (ii) time-reversal and SG symmetries lead to additional copies of the Weyl points that are located away from the Γ -∆ -A line. Using these arguments, we now determine for each SG the minimum number of Weyl points that is compatible with the symmetries and the fermion doubling theorem. To do so, we first note that we can divide the BZ into two parts, an upper part with k z ∈ [0 , π ] and a lower part with k z ∈ -[ π, 0], which are related by time-reversal symmetry. Hence, it is sufficient to analyze the Weyl point configuration in the upper part; the configuration in the lower part can be obtained by applying time-reversal symmetry (see Fig. 2). For concreteness,

SG Nos. 151 and 153 possess an additional two-fold screw rotation along the k y axis, which however does not affect the above argument. Hence, also in these SGs the minimal Weyl point multiplicity is four.

b. SG Nos. 152 and 154.SG Nos. 152 and 154 contain, besides the three-fold screw rotation C 3 ,p , a two-fold rotation C 2 along the (110) axis, which maps ( k , k x y , k z ) → ( k , x -k , y -k z ). This C 2 rotation puts extra constraints on the multiplicity of the Weyl points. To find these constraints we combine the C 2 rotation with time-reversal symmetry, yielding the mirror symmetry ( k , k x y , k z ) → -( k , k x y , k z ) [dashed lines in Fig. 2(b)]. If we now put a negative chirality Weyl point on the K -P -H line, the mirror symmetry would generate another negative chirality Weyl point at the K ′ -P ′ -H ′ line, leading to an overcancellation. Instead we place positive chirality Weyl points both at K -P -H and K ′ -P ′ -H ′ , which we then compensate by negative chirality Weyl points on the three mirror planes, as shown in Fig. 2(b). This configuration of twelve Weyl points has minimal multiplicity, satisfies all the symmetries, and has

FIG. 2. Weyl point configuration with minimal multiplicity for the six trigonal SGs with symmetry-enforced Weyl points (see Table I). The red and blue hexagons represent the top view of the upper and lower parts of the hexagonal BZ, which are related by time-reversal symmetry. The brown and green dots indicate the positions of the Weyl points with positive and negative chirality, respectively. (a) In SG Nos. 144, 145, 151, and 153 the minimal Weyl point multiplicity is four. Weyl points with the same chirality are related by time-reversal symmetry. (b) In SG Nos. 152 and 154 the minimal Weyl point multiplicity is twelve. Weyl points with the same chirality are related by time-reversal symmetry and the mirror symmetries indicated by the dashed lines.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00901v1-figure-1586c5ca-5bae-4213-905e-3b4ccf0e9d3b.jpg)

zero total chirality.

## III. NON-SYMMORPHIC NODAL LINES

Using the techniques discussed in the previous section we will now prove that glide-reflection symmetries lead to symmetry-enforced nodal lines in materials with trigonal space groups. As before, this derivation is base on (i) the momentum-dependence of the symmetry eigenvalues and (ii) the compatibility relations between irreps.

## A. Symmetry eigenvalues

We consider glide mirror symmetries M with a fractional shift of one half of the lattice constant along the z axis. Applying these symmetries twice we obtain M 2 = -ˆ, where the minus sign is due to the 2 z π rotation of the electron spin. From this it follows that the eigenvalues of the Bloch wavefunctions in the mirror planes are given by

$$M | \psi _ { \pm } ( \vec { k } ) \rangle = \pm i e ^ { - i k _ { z } / 2 } | \psi _ { \pm } ( \vec { k } ) \rangle \,. \quad \ ( 3. 1 ) \quad \ \overline { \text{TR} } \text{add}$$

The mirror planes contain either two or four TRIMs at which the Bloch bands are Kramers degenerate (blue and red dots in Fig. 3). It follows from Eq. (3.1) that the M eigenvalues of the Bloch bands at the TRIMs A L , ,

FIG. 3. (a) Band connectivity diagram for a path within the k k x z -plane of the hexagonal BZ of SG No. 158 connecting Γ to A [see panel (b)]. The color indicates the glide-mirror eigenvalues (3.1) of the Bloch bands. (b)-(d) Weyl nodal lines (purple solid and dashed curves) protected by glide-mirror symmetries for (b) SG No. 158, (c) SG No. 159 and (d) SG No. 161 are shown in the respective bulk BZs. The TRIMs are indicated by the blue and red dots and the mirror planes are shown in green and orange. For SG Nos. 158 and 161 there are two nodal lines separating the blue from the red TRIMs. Depending on the dimensionality of the irreps at the A point ( Z point) these nodal lines either touch (dashed curves) or do not touch (solid curves) the A point ( Z point). In SG No. 159 there is a nodal line in the k x = π plane enclosing a single TRIM. Note that there are different possibilites for the shape and connectivity of the nodal lines; here we only show the simplest one.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00901v1-figure-d9d78792-ac5a-4088-8a0f-a08848881656.jpg)

and Z are ± 1, while at the TRIMs Γ, M , and F they are ± i . Due to the anti-unitary time-reversal symmetry, at Γ, M , and F the Kramers pairs are formed between bands with opposite M eigenvalues (blue dots), while at A , L , and Z they are formed between bands with the same eigenvalues (red dots). Hence, as we move from a blue TRIM to a red TRIM the Kramers pairs must switch partners. In the absence of additional symmetries this leads to a non-trivial band connectivity with a symmetry-enforced band crossing [see Fig. 3(a)]. Since this is true for any path connecting a blue TRIM to a red TRIM, the Weyl band crossings occur along a onedimensional line separating the blue TRIMs from the red TRIMs [see Figs. 3(b)-(d)]. We note that the presence of additional symmetries, in particular inversion, leads to additional degeneracies, thereby invalidating the above argument. Thus, symmetry-enforced nodal lines only occur in SGs with a mirror-glide symmetry but no inversion symmetry. Among the trigonal SGs there are only three

TABLE IV. Double-valued irreps of SG No. 158 ( P 3 1) withc out time-reversal symmetry at the TRIMs Γ, A M , , and L , and the plane D (i.e., the k k x z -plane), see Fig. 3(b). We use the same convention as in Ref. [50] for labelling the irreps.

| Irreps/Sym.       | Element                      | M 010      |
|-------------------|------------------------------|------------|
| ¯ Γ 4 ¯ Γ 5 ¯ Γ 6 | - i + i ( 0 e - e - i 2 π/ 3 | iπ/ 3 0 )  |
| ¯ A 4 ¯ A 5 ¯ A 6 | +1 - 1 ( 0 e                 | iπ/ 6 )    |
| ¯ M 3 ¯ M 4       | e - iπ/ 6 - i + i            | 0          |
| ¯ D 3 ¯ D         | - ie ik z + ik z             | / 2 ie / 2 |
| 4 ¯ L 3           | +1                           |            |
| ¯                 |                              |            |
| L 4               |                              | - 1        |

SGs satisfying these criteria, namely Nos. 158, 159, and 161.

SG No. 158 ( P c 3 1) has hexagonal lattice system (Ptrigonal) and contains the glide-mirror symmetry

$$\lambda \quad \lambda _ { 0 } \colon ( x, y, z ) \to ( x, x - y, z + \frac { 1 } { 2 } ) \otimes ( - \frac { i \sqrt { 3 } } { 2 } \sigma _ { x } - \frac { i } { 2 } \sigma _ { y } ) \,, \quad \text{$pmea$} \\ M _ { 0 1 0 } \colon ( x, y, z ) \to ( x, x - y, z + \frac { 1 } { 2 } ) \otimes ( - \frac { i \sqrt { 3 } } { 2 } \sigma _ { x } - \frac { i } { 2 } \sigma _ { y } ) \,, \quad \text{$pmea$} \\ \text{which have to use to $k$ $k$ $_$n_{lano}$ invariant $Followin\n_{uart}$ to $bv$ $th$}$$

which leaves the k k x z -plane invariant. Following the above arguments we find that this glide-mirror symmetry leads to two nodal loops in the k k x z -plane, which separate the red TRIMs from the blue TRIMs, as shown in Fig. 3(b). If the irrep at the A point has dimension four (instead of two), these nodal lines touch the A point [dashed curve in Fig. 3(b)]. We note that these nodal loops must be symmetric under time-reversal symmetry and the point-group operations of SG No. 158, but are otherwise free to move within the mirror plane.

SG No. 159 ( P 31 ) has also hexagonal lattice system c (P-trigonal) and contains the glide-mirror symmetry

$$M _ { 2 1 0 } \colon ( x, y, z ) \to ( - x, - x + y, z + \frac { 1 } { 2 } ) \otimes ( \frac { i } { 2 } \sigma _ { x } + \frac { i \sqrt { 3 } } { 2 } \sigma _ { y } ) \,, \quad \, & L, \, \text{the} \\ \text{which no race both the } l \, - \, 0 \text{ and } \text{the } l \, - \, \pi \text{ balance} \quad \text{here } t$$

which leaves both the k x = 0 and the k x = π planes invariant, see Fig. 3(c). Both mirror planes have two TRIMs, between which the Kramers pairs change partners. Hence, according to the above arguments, there should be nodal lines in both mirror planes, which encircle either the Γ/ M or A L / points. However, at the A point there exist four-dimensional irreps, leading to a four-fold band degeneracy. For this reason, symmetryenforced nodal lines generically occur only in the k x = π plane, but not in the k x = 0 plane [54].

SG No. 161 ( R c 3 ) has rhombohedral lattice system (R- trigonal) and contains the glide mirror symmetry

$$\frac { 1 } { 2 } \frac { 1 } { 2 } \sigma _ { y } } { 2 } \frac { 1 } { 2 } \sigma _ { y } } { 2 } \frac { 1 } { 2 } \sigma _ { y } } { 2 } \frac { 1 } { 2 } \sigma _ { y } } { 2 } \frac { 1 } { 2 } \sigma _ { y } } { 2 } \frac { 1 } { 2 } \sigma _ { y } } { 2 } \frac { 1 } { 2 } \sigma _ { y } } { 2 } \frac { 1 } { 2 } \sigma _ { y } } { 2 } \frac { 1 } { 2 } \sigma _ { y } } { 2 } \frac { 1 } { 2 } \sigma _ { y } } { 2 } \frac { 1 } { 2 } \sigma _ { y } } { 2 } \frac { 1 } \sigma _ { y } } { 2 } \frac { 1 } { 2 } \sigma _ { y } } { 2 } \frac { 1 } { 2 } \sigma _ { y } } { 2 } \frac { 1 } { 2 } \sigma _ { y } } { 2 } \frac { 1 } { 2 } \sigma _ { y } } { 2 } \frac { 1 } { 2 } \sigma _ { y } } { 2 } \frac { 1 } { 2 } \sigma _ { y } } { 2 } \frac { 1 } { 2 } \sigma _ { y } } { 2 } \frac { 1 } { 2 } \sigma _ { y } } { 2 } \frac { 1 } { 2 } \sigma _ { y } } { 2 } \frac { 1 } { 2 } \sigma _ { y } } \frac { 1 } { 2 } \sigma _ { y } } { 2 } \frac { 1 } { 2 } \sigma _ { y } } { 2 } \frac { 1 } { 2 } \sigma _ { y } } { 2 } \frac { 1 } \sigma _ { y } } { 2 } \frac { 1 } \sigma _ { y } } { 2 } \frac { 1 } \sigma _ { y } } { 2 } \frac { 1 } \sigma _ { y } } { 2 } \frac { 1 } \sigma _ { y } } { 2 } \frac { 1 } \sigma _ { y } } { 2 } \frac { 1 } \sigma _ { y } } { 2 } \frac { 1 } \sigma _ { y } } { 2 } \frac { 1 } \sigma _ { y } } { 2 } \frac { 1 } \sigma _ { y } } { 2 } \frac { 1 } \sigma _ { y } } { 2 } \frac { 1 } \sigma _ { y } } { 2 } \frac { 1 } \sigma _ { y + 1 } } { 2 } \frac { 1 } \sigma _ { y + 1 } } { 2 } \frac { 1 } \sigma _ { y + 1 } } { 2 } \frac { 1 } \sigma _ { y + 1 } } { 2 } \frac { 1 } \sigma _ { y + 1 } } { 2 } \frac { 1 } \sigma _ { y + 1 } } { 2 } \frac { 1 } \sigma _ { y + 1 } } { 2 } \frac { 1 } \sigma _ { y + 1 } } { 2 } \frac { 1 } \sigma _ { y + 1 } } { 2 } \frac { 1 } \sigma _ { y + 1 } } { 2 } \frac { 1 } \sigma _ { y + 1 } } { 2 } \frac { 1 } \sigma _ { y + 1 } \frac { 1 } \sigma _ { y + 1 } } { 2 } \frac { 1 } \sigma _ { y + 1 } } { 2 } \frac { 1 } \sigma _ { y + 1 } } { 2 } \frac { 1 } \sigma _ { y + 1 } } { 2 } \frac { 1 } \sigma _ { y + 1 } } { 2 } \frac { 1 } \sigma _ { y + 1 } } { 2 } \frac { 1 } \sigma _ { y + 1 } } { 2 } \frac { 1 } \sigma _ { y + 1 } } { 2 } \frac { 1 } \sigma _ { y + 1 } } { 2 } \frac { 1 } \sigma _ { y + 1 } } { 2 } \frac { 1 } \sigma _ { y + 1 } } { 2 } } \frac { 1 } \sigma _ { y + 1 } } { 2 } \frac { 1 } \sigma _ { y + 1 } } { 2 } \frac { 1 } \sigma _ { y + 1 } } { 2 } \frac { 1 } \sigma _ { y + 1 } } { 2 } \frac { 1 } \sigma _ { y + 1 } } { 2 } \frac { 1 } \sigma _ { y + 1 } } { 2 } \frac { 1 } \sigma _ { y + 1 } } { 2 } \frac { 1 } \sigma _ { y + 1 } } { 2 } \frac { 1 } \sigma _ { y + 1 } } { 2 } \frac { 1 } \sigma _ { y + 1 } } { 2 } \frac { 1 } \sigma _ { y + 1 } } { 2 } { 2 } \frac { 1 } \sigma _ { y + 1 } } { 2 } \frac { 1 } \sigma _ { y + 1 } } { 2 } \frac { 1 } \sigma _ { y + 1 } } { 2 } \frac { 1 } \sigma _ { y + 1 } } { 2 } \frac { 1 } \sigma _ { y + 1 } } { 2 } \frac { 1 } \sigma _ { y + 1 } } { 2 } \frac { 1 } \sigma _ { y + 1 } } { 2 } \frac { 1 } \sigma _ { y + 1 } } { 2 } \frac { 1 } \sigma _ { y + 1 } } { 2 } \frac { 1 } \sigma _ { y + 1 } } { 2 } \frac { 1 } \sigma _ { y + 1 } } \frac { 1 } \sigma _ { y + 1 } } { 2 } \frac { 1 } \sigma _ { y + 1 } } { 2 } \frac { 1 } \sigma _ { y + 1 } } { 2 } \frac { 1 } \sigma _ { y + 1 } } { 2 } \frac { 1 } \sigma _ { y + 1 } } { 2 } \frac { 1 } \sigma _ { y + 1 } } { 2 } \frac { 1 } \sigma _ { y + 1 } } { 2 } \frac { 1 } \sigma _ { y + 1 } } { 2 } \frac { 1 } \sigma _ { y + 1 } } { 2 } \frac { 1 } \sigma _ { y + 1 } } { 2 } \frac { 1 } \sigma _ { y + 1 } } { 2 ] \sinh { 1 } } { 2 } \frac { 1 } \sinh { 1 } } { 2 } \frac { 1 } \sinh { 1 } } { 2 } \frac { 1 } \sinh { 1 } } { 2 } \frac { 1 } \sinh { 1 } } { 2 } \frac { 1 } \sinh { 1 } } { 2 } \frac { 1 } \sinh { 1 } } { 2 } \frac { 1 } \sinh { 1 } } { 2 } \frac { 1 } \sinh { 1 } } { 2 } \frac { 1 } \sinh { 1 } } { 2 } \frac { 1 } \sinh { 1 } } { 2 } \frac { 1 } \sinh { 1 } } { 2 } \frac { 1 } \sinh { 1 } } { 2 } \frac { 1 } \sinh { 1 } } { 2 } \frac { 1 } \sinh { 2 } } { 2 } \frac { 1 } \sinh { 1 } } { 2 } \frac { 1 } \sinh { 1 } } { 2 } \frac { 1 } \sinh { 1 } } { 2 } \frac { 1 } \sinh { 1 } } { 2 } \frac { 1 } \sinh { 1 } } { 2 } \frac { 1 } \sinh { 1 } } { 2 } \frac { 1 } \sinh { 1 } } { 2 } \frac { 1 } \sinh { 1 } } { 2 } \frac { 1 } \sinh { 1 } } { 2 } \frac { 1 } \sinh { 1 } } { 2 } \frac { 1 } \sinh { 1 } } { 2 } \frac { 1 } \sinh { 1 } } { 2 } \frac { 1 } \sinh { 1 } } { 2 ] \sinh { 1 } } { 2 } \frac { 1 } \sinh { 1 } } { 2 } \frac { 1 } \sinh { 1 } } { 2 } \frac { 1 } \sinh { 1 } } { 2 } \frac { 1 } \sinh { 1 } } { 2 } \frac { 1 } \sinh { 1 } } { 2 } \frac { 1 } \sinh { 1 } } { 2 } \frac { 1 } \sinh { 1 } } { 2 } \frac { 1 } \sinh { 1 } } { 2 } \frac { 1 } \sinh { 1 } } { 2 } \frac { 1 } \sinh { 1 } } { 2 } \frac { 1 } \sinh { 1 } } { 2 } \frac { 1 } \sinh { 1 } } { 2 } \frac { 0 } \sinh { 1 } } { 2 } \frac { 1 } \sinh { 1 } } { 2 } \frac { 1 } \sinh { 1 } } { 2 } \frac { 1 } \sinh { 1 } } { 2 } \frac { 1 } \sinh { 1 } } { 2 } \frac { 1 } \sinh { 1 } } { 2 } \frac { 1 } \sinh { 1 } } { 2 } \frac { 1 } \sinh { 1 } } { 2 } \frac { 1 } \sinh { 1 } } { 2 } \frac { 1 } \sinh { 1 } } { 2 } \frac { 1 } \sinh { 1 } } { 2 } \frac { 1 } \sinh { 1 } } { 2 } \frac { 1 } \sinh { 1 } } { 2 } \frac { 1 } \sinh \sinh { 1 } } { 2 } \frac { 1 } \sinh { 1 } } { 2 } \frac { 1 } \sinh { 1 } } { 2 } \frac { 1 } \sinh { 1 } } { 2 } \frac { 1 } \sinh { 1 } } { 2 } \frac { 1 } \sinh { 1 } } { 2 } \frac { 1 } \sinh { 1 } } { 2 } \frac { 1 } \sinh { 1 } } { 2 } \frac { 1 } \sinh { 1 } } { 2 } \frac { 1 } \sinh { 1 } } { 2 } \frac { 1 } \sinh { 1 } } { 2 } \frac { 1 } \sinh { 1 } } { 2 } \frac { 1 } \sinh { 1 } } { 2 } \frac { 1 } \sinh { 0 } } { 2 } \frac { 1 } \sinh { 1 } } { 2 } \frac { 1 } \sinh { 1 } } { 2 } \frac { 1 } \sinh { 1 } } { 2 } \frac { 1 } \sinh { 1 } } { 2 } \frac { 1 } \sinh { 1 } } { 2 } \frac { 1 } \sinh { 1 } } { 2 } \frac { 1 } \sinh { 1 } } { 2 } \frac { 1 } \sinh { 1 } } { 2 } \frac { 1 } \sinh { 1 } } { 2 } \frac { 1 } \sinh { 1 } } { 2 } \frac { 1 } \sinh { 1 } } { 2 } \frac { 1 } \sinh { 1 } } { 2 } \frac { 1 } \sinh { 1 } } { 2} \frac { 1 } \sinh { 1 } } { 2 } \frac { 1 } \sinh { 1 } } { 2 } \frac { 1 } \sinh { 1 } } { 2 } \frac { 1 } \sinh { 1 } } { 2 } \frac { 1 } \sinh { 1 } } { 2 } \frac { 1 } \sinh { 1 } } { 2 } \frac { 1 } \sinh { 1 } } { 2 } \frac { 1 } \sinh { 1 } } { 2 } \frac { 1 } \sinh { 1 } } { 2 } \frac { 1 } \sinh { 1 } } { 2 } \frac { 1 } \sinh { 1 } } { 2 } \frac { 1 } \sinh { 1 } } { 2 } \frac { 1 } \sinh { 1 } } { 2 } \frac { 1. \sinh { 1 } } { 2 } \frac { 1 } \sinh { 1 } } { 2 } \frac { 1 } \sinh { 1 } } { 2 } \frac { 1 } \sinh { 1 } } { 2 } \frac { 1 } \sinh { 1 } } { 2 } \frac { 1 } \sinh { 1 } } { 2 } \frac { 1 } \sinh { 1 } } { 2 } \frac { 1 } \sinh { 1 } } { 2 } \frac { 1 } \sinh { 1 } } { 2 } \frac { 1 } \sinh { 1 } } { 2 } \frac { 1 } \sinh { 1 } } { 2 } \frac { 1 } \sinh { 1 } } { 2 } \frac { 1 } \sinh { 1 } } { 2 } \frac { 1 } \sinh { { 1 } } { 2 } \frac { 1 } \sinh { 1 } } { 2 } \frac { 1 } \sinh { 1 } } { 2 } \frac { 1 } \sinh { 1 } } { 2 } \frac { 1 } \sinh { 1 } } { 2 } \frac { 1 } \sinh { 1 } } { 2 } \frac { 1 } \sinh { 1 } } { 2 } \frac { 1 } \sinh { 1 } } { 2 } \frac { 1 } \sinh { 1 } } { 2 } \frac { 1 } \sinh { 1 } } { 2 } \frac { 1 } \sinh { 1 } } { 2 } \frac { 1 } \sinh { 1 } } { 2 } \frac { 1 } \sinh { 1 } } { 2 } \frac { 1 } \sinh { 1 } } \sinh { 1 } } { 2 } \frac { 1 } \sinh { 1 } } { 2 } \frac { 1 } \sinh { 1 } } { 2 } \frac { 1 } \sinh { 1 } } { 2 } \frac { 1 } \sinh { 1 } } { 2 } \frac { 1 } \sinh { 1 } } { 2 } \frac { 1 } \sinh { 1 } } { 2 } \frac { 1 } \sinh { 1 } } { 2 } \frac { 1 } \sinh { 1 } } { 2 } \frac { 1 } \sinh { 1 } } { 2 } \frac { 1 } \sinh { 1 } } { 2 } \frac { 1 } \sinh { 1 } } { 2 } \frac { 1 } \sinh { 1 } } { 2 } \frac { 1 } \sinh{ 1 } } { 2 } \frac { 1 } \sinh { 1 } } { 2 } \frac { 1 } \sinh { 1 } } { 2 } \frac { 1 } \sinh { 1 } } { 2 } \frac { 1 } \sinh { 1 } } { 2 } \frac { 1 } \sinh { 1 } } { 2 } \frac { 1 } \sinh { 1 } } { 2 } \frac { 1 } \sinh { 1 } } { 2 } \frac {$$

which leaves the plane k x = -k y invariant, see Fig. 3(d). As in SG No. 158, this glide-mirror symmetry gives rise to two nodal lines separating the two red TRIMs from the two blue TRIMs. If the bands transform at the Z point under an irrep with dimension four (instead of two), then the nodal lines touch the Z point [dashed curve in Fig. 3(d)].

## B. Compatibility relations between irreps

The existence of symmetry-enforced nodal lines can also be derived from the compatibility relations between irreps. Here, we present this derivation for SG No. 158 ( P c 3 1). The derivation for SG Nos. 159 and 161 proceeds in a very similar way.

For SG No. 158 the relevant high-symmetry plane is the plane D (i.e., the k k x z -plane), which is left invariant under the glide-mirror symmetry M 010 , Eq. (3.2). This mirror plane contains the four TRIMs Γ, A L , , and M . To determine the connectivity of the bands we first need to derive the irreps at these four TRIMs and at the D plane. The double-valued irreps without time-reversal symmetry are listed in Table IV. We find that these irreps are all complex with the exception of Γ , which is ¯ 6 pseudoreal. Pseudoreal irreps are time-reversal invariant by themselves, while complex irreps need to be paired up to construct time-reversal invariant irreps [51, 52] (cf. Fig. 4). Note that at the A point there are both twodimensional and four-dimensional time-reversal invariant irreps. The decomposition of these time-reversal invariant irreps, as we move from a TRIM to a nearby point in the mirror plane D , can be inferred from Eq. (2.3), or alternatively from the program DCOMPREL on the Bilbao Crystallographic Server [50]. The resulting compatibility relations are depicted in Figs. 4(a), 4(c), and 4(e). Using these compatibility relations, we can now construct the band connectivity diagrams for a path within the D plane that connects two TRIMs, see Figs. 4(b), 4(d), and 4(f). We observe that for any path connecting Γ/ M to L , there must be at least one band crossing. The same is true for paths connecting Γ to A or M to A , albeit here the band crossings can be pinned at the A point, if the corresponding irrep is four dimensional. Hence, the D plane must contain at least two Weyl nodal lines enclosing the Γ and M points, as shown in Fig. 3(b).

## IV. FILLING CONSTRAINTS

According to band theory, a noninteracting band insulator can only exist if the electron filling ν is an even integer, i.e., ν ∈ 2 ◆ . Vice versa, materials with ν / ∈ 2 ◆

FIG. 4. (a), (c), (e) Kramers pairings and compatibility relations for SG No. 158 ( P 3 1) between the irreps at the TRIMs c and at a point in the mirror plane D (i.e., a point in the k k x z -plane). (b), (d), (f) Band connectivity diagrams for SG No. 158 for paths within the D plane connecting two TRIMs. Due to these nontrivial band connectivities the D plane must contain at least two Weyl nodal lines, see Fig. 3(b).

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00901v1-figure-c2b7ca2a-4911-4a0b-9afb-9d0cbb3924db.jpg)

must necessarily by (semi-)metals. In materials with nonsymmorphic symmetries, however, these filling constraints are more rigid. This is because nonsymmorphic symmetries enforce band crossings, that lead to groups of more than two connected bands. Thus, using the analysis form Secs. II and III we can derive the more rigid filling constraints. For example, for SG No. 144 (P3 ) we 1 observe that along the Γ -∆ -A line six bands form a connected group. Thus, in an insulator with SG No. 144 this group of six bands must be fully filled, i.e., the electron filling ν must be an element of 6 ◆ . Using these arguments we have derived the filling constraints for all the nonsymmorphic trigonal space groups, see fourth column of Table I. These results are in full agreement with the analysis of Ref. [44].

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00901v1-figure-87624179-a744-4a6f-adaa-13818f33b75c.jpg)

## V. EXAMPLE MATERIALS

Having identified the SGs whose nonsymmorphic symmetries enforce Weyl points and Weyl nodal lines, we can now perform a database search for materials in these SGs. We look for suitable compounds in the ICSD database [55], the AFLOW database [56], and the Materials Project database [57, 58], which yields four materials with accordion Weyl points and one material with Weyl nodal lines.

For each of these example materials we perform electronic band structure calculations with the Vienna ab initio simulation package (VASP) [59, 60], using the projector augmented wave (PAW) method [61, 62] and the PBE [63] and mBJ [64, 65] exchange-correlation functionals. As input for the ab initio calculations we use the experimental crystal structures of Refs. [66-70]. In the following we discuss the band structures along those high-symmetry lines/planes where the predicted band crossings occur. The full band structures along all highsymmetry directions of the BZ are presented in Appendix A. For two examples we also compute surface states using tight-binding models constructed from Wannier functions [71].

FIG. 5. Electronic band structure of Cu SrSnS 2 4 in SG No. 144 ( P 3 ) computed with the mBJ functional [64, 65]. 1 The band crossings along the Γ-∆-A line are symmetryenforced by the screw rotation C 3 1 , , see Sec. II.

## A. Materials with Weyl nodal points

Here, we present four materials which exhibit Weyl nodal points with accordion-like dispersions.

a. Cu SrSnS 2 4 and Cu SrGeS 2 4 The chalcogenide Cu SrSnS , which crystallizes in SG No. 144 ( 2 4 P 3 ) [66, 1 72], is an example of a material with Weyl points along the Γ-∆A line, as discussed in Sec. II. The electronic band structure computed with the mBJ functional, see Figs. 5 and 10(a), shows a band gap of about 1.1 eV, which is somewhat smaller than the experimentally reported value of 1.78 eV [72]. We observe in Fig. 5 that there are groups of 6 n connected bands between Γ and A , which form a large number of Weyl points. This is in full agreement with the band connectivity diagrams of Figs. 1(b) and 1(c). The Weyl points of Cu SrSnS 2 4 are symmetry protected by the screw rotation C 3 1 , , Eq. (2.1). Their topological stability is ensured by quantized Chern numbers. Due to the bulk boundary correspondence, these Chern numbers lead to arc surface states. Unfortunately, due to the large band gap of ∼ 1 8 eV, it will .

FIG. 6. (a) First-principles band structure of silver hydrogen phosphate Ag HPO 2 4 in SG No. 151 ( P 3 12) computed 1 with the mBJ functional. The band crossings along the Γ-∆-A line are Weyl points that are symmetry enforced by the screw rotation C 3 1 , (see Sec. II).

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00901v1-figure-e1406601-8274-4b9f-9c7b-71adf460bb69.jpg)

be very challenging to probe the bulk Weyl points or arc surface states of Cu SrSnS 2 4 using photoemission, tunneling spectroscopy, or transport experiments. Perhaps, one could try to grow Cu SrSnS 2 4 as a thin film on a conducting substrate, which might allow to probe the Weyl points using photoemission.

Cu SrGeS 2 4 is very similar to Cu SrSnS . 2 4 It crystallizes in SG No. 145 ( P 3 ) 2 [67] and its band structure shows a large number of Weyl points along the Γ-∆A line, see Fig. 10(b) in Appendix A.

b. Ag HPO 2 4 Silver hydrogen phosphate with the chemical formula Ag HPO 2 4 is a salt, that crystallizes in SG No. 151 [68]. According to our analysis of Sec. II, band structures with this SG symmetry exhibit Weyl points along the Γ-∆A line. Figures 6 and 10(c) display the first-principles band structure of Ag HPO , 2 4 computed with the mBJ functional. The band structure shows a band gap of about 2 8 . eV. Along the Γ-∆A direction we observe a group of 24 (= 6 × 4) connected bands, with many crossings, similar to the band connectivity diagrams of Fig. 1. These band crossings are Weyl points, whose topology is characterized by a nonzero Chern number. The nontrivial topology of these Weyl points gives rise to arc states at the surface, which connect Weyl points with opposite chirality. In Fig. 7 we demonstrate this for the Weyl point marked by the gray circle in Fig. 6. At the (-110) surface we observe an arc state that emanates out from the projected Weyl point.

While the ab-initio band structure of Ag HPO 2 4 nicely illustrates the predicted accordion states, this material is rather unsuitable for experimental investigations. Its large band gap and the fact that it cannot be doped easily makes it experimentally impossible to probe the bulk Weyl points and arc surface states.

c. TeAt ambient conditions, elemental tellurium (Te) has a trigonal crystal structure with SG No. 152

&lt;

l

a

t

e

x

i

s

h

1

\_

b

6

4

=

"

V

U

2

u

D

p

n

z

S

K

g

R

T

y

f

9

k

q

N

J

O

o

&gt;

A

B

X

c

L

E

r

Y

P

M

F

7

j

5

m

0

w

I

/

d

v

3

8

G

Q

W

Z

C

H

+

&lt;

l

a

t

e

x

i

s

h

1

\_

b

6

4

=

"

+

K

O

7

j

D

k

m

F

5

A

9

o

Q

W

&gt;

B

/

H

c

V

N

S

8

E

3

U

r

q

y

C

p

M

w

X

2

Z

d

P

f

u

G

z

0

Y

n

J

R

L

g

I

v

T

&lt;

l

a

t

e

x

i

s

h

1

\_

b

6

4

=

"

N

m

D

k

z

I

Q

L

C

u

q

n

W

o

w

0

A

3

&gt;

B

/

c

V

S

8

E

X

5

y

p

K

M

r

Y

+

g

R

F

P

v

H

2

7

G

T

d

9

Z

J

O

f

U

j

FIG. 7. Momentum-dependent surface density of states for the (-110) surface of Ag HPO 2 4 at the energy E -E F = -0 704 . eV. The crystal is terminated such that the silver atoms 'Ag3', 'Ag4', and 'Ag5' are at the top layer (see crystal structure in inset). Yellow and black correspond to high and low densities, respectively. The location of the projected Weyl point is marked by a gray circle, cf. Fig. 6. For a plot of the (-110) surface BZ see Fig. 1(a).

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00901v1-figure-dd067bda-a348-43a7-b352-5a30a8b2b3e9.jpg)

( P 3 21) 1 [69, 73]. In this structure the Te atoms form helical chains along the z direction, which are arranged in a hexagonal array [74]. Recent first principles calculation have shown that the valence and conduction bands of Te at the H point exhibit several Weyl points [75, 76]. Moreover, it was found that hydrostatic or uniaxial strain leads to a band inversion at the H point, thereby transforming Te into a strong topological insulator [77].

It follows from our analysis of Sec. II, that Te exhibits Weyl points not only at the H point, but also along the Γ-∆A line, which is left invariant by the screw rotation C 3 1 , . This is clearly visible in the ab-initio band structure calculations of Figs. 8 and 10(d), which show Weyl points with accordion-like dispersions between Γ and A . As before, we observe that the band connectivity fully agrees with the theoretical prediction of Fig. 1. By the bulk boundary correspondence, the Weyl points at the H point and on the Γ-∆A line lead to numerous arc surface states, see Appendix B. These arc surface states should be readily observable in photoemission and scanning tunneling experiments.

In closing, we note that trigonal selenium (Se), which also crystalizes in SG No. 152 ( P 3 21) [78], has a very 1 similar band structure to tellurium, albeit the spin-orbit coupling is considerably smaller [Fig. 10(e)]. Therefore, the bands forming Weyl points along the Γ-∆-A line of Se are only separated by less than ∼ 50 meV. This makes the experimental detection of the accordion Weyl points in Se more challenging than in Te.

FIG. 8. Electronic band structure of trigonal tellurium in SG No. 152 ( P 3 21), computed with the PBE functional. 1 The band crossings with accordion dispersions along the Γ∆-A line are symmetry enforced by the screw rotation C 3 1 , (see Sec. II).

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00901v1-figure-8dcb558d-f668-4339-8cad-d294aab31ff3.jpg)

## B. Weyl nodal lines in Te 16 Si 38

The tellurium-silicon clathrate Te 16 Si 38 crystallizes in two different forms. While the low-temperature form is cubic, the high-temperature form is rhombohedral with trigonal SG No. 161 ( R c 3 ) [70, 79]. According to our analysis of Sec. III, materials in this SG exhibit Weyl nodal lines that are protected by the glide mirror symmetry M 110 . In Fig. 9 we present the first principles band structure (with PBE functional) of rhombohedral Te 16 Si 38 along high-symmetry paths of the k x = -k y plane, which is left invariant by M 110 [see also Fig. 10(f)]. We observe that along any path that connects a blue TRIM to a red TRIM in the BZ of Fig. 3(d), 4 n bands form a connected group with many Weyl crossings. These Weyl crossings are part of Weyl nodal lines that separate the blue from the red TRIMs in the k x = -k y plane. The stability of these Weyl nodal lines is guaranteed by a π -Berry phase, which by the bulk-boundary correspondence gives rise to drumhead surface states [80]. In addition, the bands of rhombohedral Te 16 Si 38 carry a nonzero Berry curvature which becomes especially large close to the Weyl nodal lines.

## VI. CONCLUSIONS

In summary, we have performed a systematic analysis of nonsymmorphic band crossings in nonmagnetic trigonal materials with strong spin-orbit coupling. We have found that among the 25 trigonal space groups, six support Weyl points, while three support Weyl nodal lines (see Table I). The electronic bands at these Weyl points and Weyl lines are required to cross due to nonsymmorphic mirror or rotation symmetries alone, regardless of the chemical composition and other material details.

FIG. 9. DFT-PBE band structure of rhombohedral Te 16 Si 38 along high-symmetry paths of the k x = -k y plane of the rhombohedral BZ [see Fig. 3(d)]. The band crossings are part of Weyl nodal lines, which are protected by the glide mirror symmetry M 110 .

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00901v1-figure-ab5540a9-4cd1-4cf6-968e-eaba3456727f.jpg)

Hence, these Weyl crossings occur in all bands of any material crystallizing in the space groups listed in Table I. Using this insight, we have performed a database search and identified several existing materials that possess the predicted topological band crossings, see last column in Table I. The nontrivial topology of the band crossings in these materials leads to a number of interesting and experimentally observable phenomena, e.g., arc and drumhead surface states, anomalous magneto-transport properties, and anomalous Hall effects. Particularly interesting are the nonsymmorphic band crossings and surface states of tellurium (Figs. 8 and 11). We hope that our findings will encourage experimentalists to study the topological properties of these compounds.

## ACKNOWLEDGMENTS

The authors thank C. Ast, K. von Klitzing, A. Topp, M. G. Vergniory, and A. Yaresko for useful discussions. B.K. thanks the Max Planck Institute for Solid State Research in Stuttgart for financial support. C.-K.C. is supported by the Strategic Priority Research Program of the Chinese Academy of Sciences (Grant XDB28000000). Work at Princeton was supported by NSF through the Princeton Center for Complex Materials, a Materials Research Science and Engineering Center DMR-1420541.

## Appendix A: Additional band structure calculations

In this Appendix, we present additional band-structure calculations for the example materials studied in Sec. V of the main text. In Figs. 10(a), 10(b), and 10(c),

FIG. 10. (a), (b), (c) DFT band structures computed with the mBJ functional for Cu SrSnS , Cu SrGeS , and Ag HPO . 2 4 2 4 2 4 Along the Γ -∆ -A line, these materials exhibit Weyl points with accordion dispersions. For a zoom in of the band structure near the Weyl points see Sec. V A in the main text. (d), (e) DFT-PBE band structure of trigonal tellurium and selenium. Along the Γ -∆ -A direction there are groups of six connected bands, which cross at several Weyl points. (f) DFT-PBE band structure for Te 16 Si 38 , which exibits Weyl nodal lines within the k x = -k y plane. For a zoom in of the band structure around the Weyl nodal lines see Fig. 9 in the main text.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00901v1-figure-350eccce-c9f1-4329-a10e-973709871423.jpg)

we present the full band structures of the materials Cu SrSnS , 2 4 Cu SrGeS , 2 4 and Ag HPO , 2 4 respectively, which exhibit Weyl points along Γ -∆ -A . These Weyl points are symmetry enforced by the three-fold C 3 ,p screw rotations. In the chalcogenides Cu SrSnS 2 4 and Cu SrGeS 2 4 the band gap is about 1.5 eV, while in the salt Ag HPO 2 4 it is about 2.8 eV. This makes it challenging to observe the Weyl points in photoemission or electron tunneling experiments. However, it might be perhaps possible to grow thin films of Cu 2 SrSnS 4 and Cu SrGeS 2 4 on a conducting substrate to introduce charge carriers.

Figures 10(d) and 10(e) display the full band structures of trigonal tellurium and selenium. Along the Γ -∆ -A line there are groups of six connected bands with accordion dispersions. Since tellurium is metallic, it should be possible to observe its accordion dispersion and Weyl points using angle-resolved photoemission experiments. Tellurium has many Weyl points also away from the Γ -∆ -A line, for example, at the H point [75, 76]. These give rise to numerous arc surface states, see Appendix B.

Figure 10(f) shows the full band structure of the tellurium-silicon clathrate Te 16 Si 38 , which exhibits Weyl nodal lines in the k x = -k y plane. We observe that along the paths Γ -L , Γ -Z , and F -Z , the bands form groups of 4 n connected bands with Weyl crossings. These Weyl crossings are part of nodal lines that separate Γ and F

from Z and L in the k x = -k y plane of the rhombohedral BZ, see Fig. 4(d).

## Appendix B: Surface states of tellurium

Trigonal tellurium has many Weyl points and associated arc surface states. In Fig. 11 we present the surface state dispersion for the (-110) surface of Te [81]. The labels of the high-symmetry points in the (-110) surface BZ are given in Fig. 1(a). We observe that there are several surface states emanating from the Weyl points near the H point. Remarkably, these surface states cross the entire surface BZ and lie within the bulk band gap, in between the valence and conduction bands. Since their energy is close to the Fermi level, they could play an important role for the (re-)interpretation of the unusual surface transport and photoemission measurements reported in Refs. [82-84].

There are also arc surface states that are associated with the accordion bands discussed in Sec. V A 0 c and Fig. 8. For example, for the topmost Weyl point in the valence bands at the energy E -E F = -1 072 (gray circle . in Fig. 8), we observe an arc state that disperses towards the Γ point. ˜

&lt;

l

a

t

e

x

i

s

h

1

\_

b

6

4

=

"

5

D

o

B

8

k

C

M

v

z

S

2

A

9

W

g

&gt;

+

H

c

V

L

N

E

O

G

j

q

U

Q

Y

3

F

Z

J

7

w

d

P

p

r

T

X

n

f

0

R

K

m

u

y

I

/

&lt;

l

a

t

e

x

i

s

h

1

\_

b

6

4

=

"

5

D

o

B

8

k

C

M

v

z

S

2

A

9

W

g

&gt;

+

H

c

V

L

N

E

O

G

j

q

U

Q

Y

3

F

Z

J

7

w

d

P

p

r

T

X

n

f

0

R

K

m

u

y

I

/

&lt;

l

a

t

e

x

i

s

h

1

\_

b

6

4

=

"

5

D

o

B

8

k

C

M

v

z

S

2

A

9

W

g

&gt;

+

H

c

V

L

N

E

O

G

j

q

U

Q

Y

3

F

Z

J

7

w

d

P

p

r

T

X

n

f

0

R

K

m

u

y

I

/

&lt;

l

a

t

e

x

i

s

h

1

\_

b

6

4

=

"

5

D

o

B

8

k

C

M

v

z

S

2

A

9

W

g

&gt;

+

H

c

V

L

N

E

O

G

j

q

U

Q

Y

3

F

Z

J

7

w

d

P

p

r

T

X

n

f

0

R

K

m

u

y

I

/

&lt;

l

a

t

e

x

i

s

h

1

\_

b

6

4

=

"

5

D

o

B

8

k

C

M

v

z

S

2

A

9

W

g

&gt;

+

H

c

V

L

N

E

O

G

j

q

U

Q

Y

3

F

Z

J

7

w

d

P

p

r

T

X

n

f

0

R

K

m

u

y

I

/

˜ ˜ ˜ ˜ ˜ FIG. 11. Energy dispersion of the surface states of trigonal tellurium for the (-110) surface along high-symmetry momentum directions, cf. Fig. 1(a). The color grading indicates the surface density of states, with yellow and black corresponding to high and low densities, respectively. The gray circle denotes the projected position of the Weyl point at the energy E -E F = -1 072 eV, see Fig. 8. .

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00901v1-figure-f3532007-5e65-45d3-900a-7fc2702a3f53.jpg)

- [1] J. von Neumann and E. Wigner, Z. Physik 30 , 467 (1929).
- [2] C.-K. Chiu, J. C. Y. Teo, A. P. Schnyder, and S. Ryu, Rev. Mod. Phys. 88 , 035005 (2016).
- [3] G. E. Volovik, Topology of quantum vacuum , Lecture Notes in Physics, Vol. 870 (Springer Berlin, 2013) pp. 343-383.
- [4] N. P. Armitage, E. J. Mele, and A. Vishwanath, ArXiv e-prints (2017), arXiv:1705.01111 [cond-mat.str-el].
- [5] A. Burkov, Annual Review of Condensed Matter Physics 9 , 359 (2018), https://doi.org/10.1146/annurevconmatphys-033117-054129.
- [6] S.-Y. Yang, H. Yang, E. Derunova, S. S. P. Parkin, B. Yan, and M. N. Ali, ArXiv e-prints (2017), arXiv:1707.04523 [cond-mat.mtrl-sci].
- [7] X. Wan, A. M. Turner, A. Vishwanath, and S. Y. Savrasov, Phys. Rev. B 83 , 205101 (2011).
- [8] A. A. Burkov, M. D. Hook, and L. Balents, Phys. Rev. B 84 , 235126 (2011).

lopolski, D. S. Sanchez, M. Neupane, N. Alidoust, C. Liu, B. Wang, C.-C. Lee, H.-T. Jeng, C. Zhang, Z. Yuan, S. Jia, A. Bansil, F. Chou, H. Lin, and M. Z. Hasan, Nature Communications 7 (2016).

- [13] G. Bian, T.-R. Chang, H. Zheng, S. Velury, S.-Y. Xu, T. Neupert, C.-K. Chiu, S.-M. Huang, D. S. Sanchez, I. Belopolski, N. Alidoust, P.-J. Chen, G. Chang, A. Bansil, H.-T. Jeng, H. Lin, and M. Z. Hasan, Phys. Rev. B 93 , 121113 (2016).
- [14] W. B. Rui, Y. X. Zhao, and A. P. Schnyder, Phys. Rev. B 97 , 161113 (2018).
- [15] T. Bzdusek, Q. Wu, A. R¨egg, M. Sigrist, u and A. A. Soluyanov, Nature 538 , 75 (2016).
- [16] Z. Shi, M. Wang, and J. Wu, Applied Physics Letters 107 , 102403 (2015), https://doi.org/10.1063/1.4930875.
- [17] J. R. Schaibley, H. Yu, G. Clark, P. Rivera, J. S. Ross, K. L. Seyler, W. Yao, and X. Xu, Nature Reviews Materials 1 (2016).
- [18] C.-K. Chiu and A. P. Schnyder, Phys. Rev. B 90 , 205136 (2014).
- [19] K. Shiozaki and M. Sato, Phys. Rev. B 90 , 165114 (2014).
- [20] Y. X. Zhao, A. P. Schnyder, and Z. D. Wang, Phys. Rev. Lett. 116 , 156402 (2016).
- [9] S.-Y. Xu, C. Liu, S. K. Kushwaha, R. Sankar, J. W. Krizan, I. Belopolski, M. Neupane, G. Bian, N. Alidoust, T.-R. Chang, H.-T. Jeng, C.-Y. Huang, W.-F. Tsai, H. Lin, P. P. Shibayev, F.-C. Chou, R. J. Cava, and M. Z. Hasan, Science 347 , 294 (2015), http://science.sciencemag.org/content/347/6219/294.full.pdf.
- [10] S.-M. Huang, S.-Y. Xu, I. Belopolski, C.-C. Lee, G. Chang, B. Wang, N. Alidoust, G. Bian, M. Neupane, C. Zhang, S. Jia, A. Bansil, H. Lin, and M. Z. Hasan, Nature Communications 6 (2015).
- [21] A. H. Castro Neto, F. Guinea, N. M. R. Peres, K. S. Novoselov, and A. K. Geim, Rev. Mod. Phys. 81 , 109 (2009).
- [22] Z. Wang, Y. Sun, X.-Q. Chen, C. Franchini, G. Xu, H. Weng, X. Dai, and Z. Fang, Phys. Rev. B 85 , 195320 (2012).
- [11] Y.-H. Chan, C.-K. Chiu, M. Y. Chou, and A. P. Schnyder, Phys. Rev. B 93 , 205132 (2016).
- [12] G. Bian, T.-R. Chang, R. Sankar, S.-Y. Xu, H. Zheng, T. Neupert, C.-K. Chiu, S.-M. Huang, G. Chang, I. Be-
- [23] C.-K. Chiu and A. P. Schnyder, Journal of Physics: Conference Series 603 , 012002 (2015).
- [24] Z. Wang, H. Weng, Q. Wu, X. Dai, and Z. Fang, Phys. Rev. B 88 , 125427 (2013).

- [25] L. S. Xie, L. M. Schoop, E. M. Seibel, Q. D. Gibson, W. Xie, and R. J. Cava, APL Mater. 3 , 083602 (2015), http://dx.doi.org/10.1063/1.4926545.
- [26] A. Yamakage, Y. Yamakawa, Y. Tanaka, and Y. Okamoto, Journal of the Physical Society of Japan 85 , 013708 (2016), http://dx.doi.org/10.7566/JPSJ.85.013708.
- [27] S. M. Young, S. Zaheer, J. C. Y. Teo, C. L. Kane, E. J. Mele, and A. M. Rappe, Phys. Rev. Lett. 108 , 140405 (2012).
- [28] L. M. Schoop, M. N. Ali, C. Straßer, A. Topp, A. Varykhalov, D. Marchenko, V. Duppel, S. S. P. Parkin, B. V. Lotsch, and C. R. Ast, Nature Communications 7 (2016).
- [29] Y. X. Zhao and A. P. Schnyder, Phys. Rev. B 94 , 195109 (2016).
- [30] L. Michel and J. Zak, Phys. Rev. B 59 , 5998 (1999).
- [31] S. M. Young and C. L. Kane, Phys. Rev. Lett. 115 , 126803 (2015).
- [32] A. Alexandradinata, Z. Wang, and B. A. Bernevig, Phys. Rev. X 6 , 021008 (2016).
- [33] A. Furusaki, Science Bulletin 62 , 788 (2017).
- [34] R. Takahashi, M. Hirayama, and S. Murakami, Phys. Rev. B 96 , 155206 (2017).
- [35] B.-J. Yang, T. A. Bojesen, T. Morimoto, and A. Furusaki, Phys. Rev. B 95 , 075135 (2017).
- [36] C. Fang, Y. Chen, H.-Y. Kee, and L. Fu, Phys. Rev. B 92 , 081201 (2015).
- [37] M. Malard, P. E. de Brito, S. Ostlund, and H. Johannesson, ArXiv e-prints (2018), arXiv:1804.09736 [condmat.str-el].
- [38] J. Zhang, Y.-H. Chan, C.-K. Chiu, M. G. Vergniory, L. M. Schoop, and A. P. Schnyder, Phys. Rev. Materials 2 , 074201 (2018).
- [39] R. Chen, H. C. Po, J. B. Neaton, and A. Vishwanath, Nature Physics 14 (2017).
- [40] B. Bradlyn, L. Elcoro, J. Cano, M. G. Vergniory, Z. Wang, C. Felser, M. I. Aroyo, and B. A. Bernevig, Nature 547 , 298 (2017).
- [41] L. M. Schoop, A. Topp, J. Lippmann, F. Orlandi, L. M¨ uchler, M. G. Vergniory, Y. Sun, A. W. Rost, V. Duppel, M. Krivenkov, S. Sheoran, P. Manuel, A. Varykhalov, B. Yan, R. K. Kremer, C. R. Ast, and B. V. Lotsch, 4 (2018), 10.1126/sciadv.aar2317.
- [42] A. Topp, J. M. Lippmann, A. Varykhalov, V. Duppel, B. V. Lotsch, C. R. Ast, and L. M. Schoop, New Journal of Physics 18 , 125014 (2016).
- [43] S. A. Parameswaran, A. M. Turner, D. P. Arovas, and A. Vishwanath, Nature Physics 9 (2013).
- [44] H. Watanabe, H. C. Po, M. P. Zaletel, and A. Vishwanath, Phys. Rev. Lett. 117 , 096404 (2016).
- [45] T. Zhang, Y. Jiang, Z. Song, H. Huang, Y. He, Z. Fang, H. Weng, and C. Fang, Nature 566 , 475 (2019).
- [46] M. G. Vergniory, L. Elcoro, C. Felser, N. Regnault, B. A. Bernevig, and Z. Wang, Nature 566 , 480 (2019).
- [47] F. Tang, H. C. Po, A. Vishwanath, and X. Wan, Nature 566 , 486 (2019).
- [48] L. Michel and J. Zak, Physics Reports 341 , 377 (2001), symmetry, invariants, topology.
- [49] J. Kruthoff, J. de Boer, J. van Wezel, C. L. Kane, and R.-J. Slager, Phys. Rev. X 7 , 041069 (2017).
- [50] L. Elcoro, B. Bradlyn, Z. Wang, M. G. Vergniory, J. Cano, C. Felser, B. A. Bernevig, D. Orobengoa, G. de la Flor, and M. I. Aroyo, Journal of Applied Crys-
- tallography 50 , 1457 (2017).
- [51] C. Bradley and A. P. Cracknell, The Mathematical Theory of Symmetry in Solids: Representation theory for point groups and space groups (Clarendon Press, Oxford, 1972).
- [52] S. C. Miller and W. F. Love, Tables of Irreducible Representations of Space Groups and Co-representations of Magnetic Groups (Boulder: Pruett, 1967).
- [53] H. Nielsen and M. Ninomiya, Nuclear Physics B 185 , 20 (1981).
- [54] We note that at the A point there exist both twodimensional and four-dimensional irreps. Hence, nodal lines are also possible in the k x = 0 plane, but they are not guaranteed to exist by symmetry alone.
- [55] Https://icsd.fiz-karlsruhe.de/.
- [56] S. Curtarolo, W. Setyawan, S. Wang, J. Xue, K. Yang, R. H. Taylor, L. J. Nelson, G. L. Hart, S. Sanvito, M. Buongiorno-Nardelli, N. Mingo, and O. Levy, Computational Materials Science 58 , 227 (2012).
- [57] A. Jain, S. P. Ong, G. Hautier, W. Chen, W. D. Richards, S. Dacek, S. Cholia, D. Gunter, D. Skinner, G. Ceder, and K. a. Persson, APL Materials 1 , 011002 (2013).
- [58] 'Materials Project,' .
- [59] G. Kresse and J. Furthm¨ller, Phys. Rev. B u 54 , 11169 (1996).
- [60] G. Kresse and J. Furthm¨ller, Computational Materials u Science 6 , 15 (1996).
- [61] P. E. Bl¨chl, Phys. Rev. B o 50 , 17953 (1994).
- [62] G. Kresse and D. Joubert, Phys. Rev. B 59 , 1758 (1999).
- [63] J. P. Perdew, K. Burke, and M. Ernzerhof, Phys. Rev. Lett. 77 , 3865 (1996).
- [64] A. D. Becke and E. R. Johnson, The Journal of Chemical Physics 124 , 221101 (2006), https://doi.org/10.1063/1.2213970.
- [65] F. Tran and P. Blaha, Phys. Rev. Lett. 102 , 226401 (2009).
- [66] C. L. Teske, Zeitschrift f¨r u anorganische und allgemeine Chemie 419 , 67 (1976), https://onlinelibrary.wiley.com/doi/pdf/10.1002/zaac.19764190112.
- [67] C. L. Teske, Zeitschrift fr Naturforschung B 34 , 386 (1979).
- [68] I. Tordjman, A. Boudjada, J.-C. Guitel, and R. Masse, Acta Crystallographica Section B 34 , 3723 (1978), https://onlinelibrary.wiley.com/doi/pdf/10.1107/S056774087801196
- [69] C. Adenis, V. Langer, and O. Lindqvist, Acta Crystallographica Section C 45 , 941 (1989).
- [70] N. Jaussaud, P. Toulemonde, M. Pouchard, A. SanMiguel, P. Gravereau, S. Pechev, G. Goglio, and C. Cros, Solid State Sciences 6 , 401 (2004).
- [71] A. A. Mostofi, J. R. Yates, Y.-S. Lee, I. Souza, D. Vanderbilt, and N. Marzari, Computer Physics Communications 178 , 685 (2008).
- [72] Z. Tong, J. Yuan, J. Chen, A. Wu, W. Huang, C. Han, Q. Cai, C. Ma, Y. Liu, L. Fang, and Z. Liu, Materials Letters 237 , 130 (2019).
- [73] P. Cherin and P. Unger, Acta Crystallographica 23 , 670 (1967).
- [74] Note that depending on the chirality of these chains Te can also be in SG No. 153 ( P 3 12). 2
- [75] M. Hirayama, R. Okugawa, S. Ishibashi, S. Murakami, and T. Miyake, Phys. Rev. Lett. 114 , 206401 (2015).
- [76] K. Nakayama, M. Kuno, K. Yamauchi, S. Souma, K. Sugawara, T. Oguchi, T. Sato, and T. Takahashi, Phys. Rev. B 95 , 125204 (2017).

- [77] L. A. Agapito, N. Kioussis, W. A. Goddard, and N. P. Ong, Phys. Rev. Lett. 110 , 176401 (2013).
- [78] P. Cherin and P. Unger, Inorganic Chemistry 6 , 1589 (1967), https://doi.org/10.1021/ic50054a037.
- [79] N. Jaussaud, M. Pouchard, P. Gravereau, S. Pechev, G. Goglio, C. Cros, A. San Miguel, and P. Toulemonde, Inorganic Chemistry 44 , 2210 (2005).
- [80] C. K. Chiu, Y. H. Chan, and A. P. Schnyder, arXiv e-prints , arXiv:1810.04094 (2018), arXiv:1810.04094 [cond-mat.mes-hall].
- [81] Note that there is only one type of termination.
- [82] K. von Klitzing and G. Landwehr, Solid State Communications 9 , 2201 (1971).
- [83] L. D. Laude, R. F. Willis, and B. Fitton, Phys. Rev. Lett. 29 , 472 (1972).
- [84] T. Englert, K. von Klitzing, R. Silbermann, and G. Landwehr, physica status solidi (b) 81 , 119 (1977), https://onlinelibrary.wiley.com/doi/pdf/10.1002/pssb.2220810110.
Y. -H. Chan Berkay Kilic Moritz M. Hirschmann Ching-Kai Chiu Leslie M. Schoop Darshan G. Joshi Andreas P. Schnyder
Symmetry-enforced band crossings in trigonal materials: Accordion states and Weyl nodal lines
Nonsymmoprhic symmetries, such as screw rotations or glide reflections, can
enforce band crossings within high-symmetry lines or planes of the Brillouin
zone. When these band degeneracies are close to the Fermi energy, they can give
rise to a number of unusual phenomena: e.g., anomalous magnetoelectric
responses, transverse Hall currents, and exotic surface states. In this paper,
we present a comprehensive classification of such nonsymmorphic band crossings
in trigonal materials with strong spin-orbit coupling. We find that in trigonal
systems there are two different types of nonsymmorphic band degeneracies: (i)
Weyl points protected by screw rotations with an accordion-like dispersion, and
(ii) Weyl nodal lines protected by glide reflections. We report a number of
existing materials, where these band crossings are realized near the Fermi
energy. This includes Cu2SrSnS4 and elemental tellurium (Te), which exhibit
accordion Weyl points; and the tellurium-silicon clathrate Te16Si38, which
shows Weyl nodal lines. The ab-initio band structures and surface states of
these materials are studied in detail, and implications for experiments are
briefly discussed.