Hourglass Dirac chain metal in rhenium dioxide

Nonsymmorphic symmetries, which involve fractional lattice translations, can generate exotic types of fermionic excitations in crystalline materials. Here we propose a topological phase arising from nonsymmorphic symmetries—the hourglass Dirac chain metal, and predict its realization in the rhenium dioxide. We show that ReO2 features hourglass-type dispersion in the bulk electronic structure dictated by its nonsymmorphic space group. Due to time reversal and inversion symmetries, each band has an additional two-fold degeneracy, making the neck crossing-point of the hourglass four-fold degenerate. Remarkably, close to the Fermi level, the neck crossing-point traces out a Dirac chain—a chain of connected four-fold-degenerate Dirac loops—in the momentum space. The symmetry protection, the transformation under symmetry-breaking, and the associated topological surface states of the Dirac chain are revealed. Our results open the door to an unknown class of topological matters, and provide a platform to explore their intriguing physics.

Although the essential band-crossings are solely determined by the space group for which theoretical analysis has offered valuable guidelines, the search for realistic materials that exhibit them at low energy is still challenging. This is because the bands in real materials typically have complicated 3D dispersions, such that the , which consists of one (red) loop in kz = π plane and one (green) loop in kx = π plane. There is another isolated Dirac point (orange dot) on T-Y. These crossings are fourfold degenerate and correspond to the neck crossing-point of the hourglass-type dispersion. For example, (b) shows the schematic band dispersion along a path on the kx = π plane connecting R and P (an arbitrary point on U-X). Each band is two-fold degenerate, and the neck point (green dot) is fourfold degenerate.
crossing point that we are chasing may be far away from the Fermi energy. The situation could be even worse for nodal loops, since the points on the loop are not guaranteed to have the same energy, there might be large energy variation around the loop. So far, the proposed nonsymmorphic topological metals are still limited, therefore, it is urgent to discover more suitable candidate materials to expedite studies of their intriguing properties.
In this paper, based on first-principles calculations and symmetry analysis, we report on hourglass-type essential band-crossings in an existing material-ReO 2 . We show that due to time reversal (T ) and inversion (P) symmetries, the hourglass here is actually doubled, and the neck crossing-point here becomes a Dirac point with fourfold degeneracy [ Fig. 1(b)]. Remarkably, close to Fermi level, the neck point traces out a Dirac chain-a chain of connected (four-fold-degenerate) Dirac loops-in the momentum space, as schematically shown in Fig. 1(a), hence the state may be dubbed as an hourglass Dirac chain metal. This chain is essential, robust against SOC, and dictated by two orthogonal glide mirror planes com-arXiv:1705.01424v1 [cond-mat.mes-hall] 3 May 2017 bined with time-reversal and inversion symmetries. In addition, there is another pair of bulk hourglass Dirac points on a symmetry line [see Fig. 1(a)]. We clarify the protection of these exotic band-crossings, and discuss their transformations under symmetry-breaking as well as the associated topological surface states. Our findings provide an exciting platform for explore the novel topological fermions from nonsymmorphic symmetries.
Single crystal ReO 2 is observed with two structures denoted as α and β [48,49]. β-ReO 2 is energetically more stable, and is found to be a stable paramagnetic metal at ambient conditions [50]. Hence we focus on β-ReO 2 here. It adopts the PbO 2 -type orthorhombic crystal structure with space group No. 60 (P bcn) [48]. As shown in Fig. 2(a), the structure is characterized by zigzag chains of Re atoms running along the c-axis, and each Re atom is contained in a slightly distorted octahedron of six surrounding O atoms. The space group of the structure may be generated by the following symmetry operations that will be important in our discussion: the inversion P, and two glide mirror planes involving half Here the tilde above a symbol indicates that it is a nonsymmorphic symmetry. One also notes that combining all three operations leads to a third glide mirror M y : (x, y, z) → (x, −y, z + 1 2 ). We performed first-principles calculations based on the density functional theory (DFT). SOC was included, and possible correlation effect of Re(5d) orbitals was tested. The calculation details are in the Supplemental Material [51]. The experimental values of the lattice parameters (a = 4.809Å, b = 5.643Å, c = 4.601Å) [48] were used in the calculation.
In octahedral crystal field, Re(5d) orbitals are split into t 2g and e g groups, with the latter at higher energy. For Re 4+ with 3 valence electrons, the Re-t 2g orbitals will be half-filled, resulting in a metallic state. Figure 2(c) shows the calculated band structure of ReO 2 along with the projected density of states (PDOS). Indeed, one observes a metallic phase with fairly dispersive bands around Fermi level, and the low-energy states are dominated by the Re-t 2g orbitals. Understanding that each band is at least two-fold degenerate due to the presence of T and P, two interesting type of band features can be observed from Let's first investigate feature (i) regarding the four-fold degeneracy along the three high-symmetry lines. Consider the U-X line at k x = π and k y = 0 (in unit of the inverse of respective lattice parameter). It is an invariant subspace of M x , so each Bloch state |u there can be chosen as an eigenstate of M x . Since the M x eigenvalue g x must be ±i on U-X. Here T 010 denotes the translation along y by one unit cell, and E is the 2π spin rotation. The commutation relation between M x and P given by means that { M x , P} = 0 on U-X. Consequently, each state |u and its Kramers-degenerate partner PT |u must share the same M x eigenvalue. For example, assume |u has g x = +i (denoted as | + i ), then where in the second step we used the fact that T is an anti-unitary operator. Same result holds for a state with g x = −i. On the other hand, U-X is invariant under another anti-unitary symmetry M z T , which also generates a Kramers-like degeneracy since ( M z T ) 2 = −1 on U-X. Note that  (2) and M z P = T 111 P M z , one finds that so the Kramers pair |u and PT |u at any k-point on U-R share the same (g x , g z ) eigenvalues. In addition, points R and U are time-reversal invariant momenta. At R = (π, π, π), (g x , g z ) = (±1, ±i), hence if |u has eigenvalues (g x , g z ), its Kramers partner T |u must have (g x , −g z ).
Focusing on the eigenvalue g x , the analysis shows that the four states in the degenerate quartet (may be chosen as {|u , T |u , P|u , PT |u }) at R all have the same g x (+1 or −1); whereas at point U, they consist of two states with g x = +i and two other states with g x = −i. Hence there has to be a switch of partners between two quartets along U-R, during which the eight bands must be entangled to form the hourglass-type dispersion. The situation is schematically shown in Fig. 3(a). It is important to note that the four-fold-degenerate neck crossingpoint (denoted as D on U-R) is protected because the two crossing doubly-degenerate bands have opposite g x [with each degenerate pair sharing the same g x , as shown in Eq. (5) and illustrated in Fig. 3(a)].
Furthermore, since the whole k x = π plane is invariant under M x , g x is well defined for any state on this plane. Hence the above argument applies to any path lying on the k x = π plane and connecting points U and R, which should feature an hourglass spectrum with fourfold-degenerate crossing-point in between. The crossingpoint must trace out a closed Dirac loop L 1 on this plane, as indicated in Fig. 3(b). One also notes that not only U, actually any point P on U-X has four-fold degeneracy with two g x = +i and two g x = −i, as we analyzed before. Thus hourglass pattern is guaranteed to appear on any path connecting R to an arbitrary point on U-X [ Fig. 3(b)].
Similar analysis as in the last two paragraphs applies to the k z = π plane, with the role played by M x replaced by M z . It shows that hourglass pattern appears on any path connecting U to an arbitrary point on Z-T or T-R, and the neck point of the hourglass traces out a second Dirac loop L 2 , as illustrated in Fig. 3(c). Interestingly, L 1 and L 2 are orthogonal to each other, and they touch at the point D on the U-R line. Thus they constitute a Dirac chain in the momentum space, as shown in Fig. 1(a). Figure 3(d,e) shows the locations of the Dirac loops obtained from DFT calculations, which are consistent with our symmetry analysis. The chain is close to the Fermi level and has small energy variation (< 0.2 eV). We stress that the presence of such band-crossing pattern is solely determined by the space group (plus T ). However, whether such crossings could manifest around Fermi level and have relatively small energy variation will depend on the specific material.
Up to now, one may wonder whether there exists a third loop on the k y = π plane, given that M y is also a symmetry. It turns out not to be the case. Consider any state |g y on k y = π plane with M y eigenvalue g y , one can show that M y (PT |g y ) = −g y (PT |g y ).
Thus each Kramers pair |u and PT |u have opposite g y , which is in contrast with Eq. The hourglass dispersion and the Dirac chain are dictated by symmetry. They must be kept as long as the space group symmetry is maintained. In Fig. 4(a,b), we demonstrate that when we distort the crystal lattice while maintaining the symmetry, the shape and the size of the chain can change, but it cannot be destroyed. In contrast, if we break the symmetry, e.g., by varying the angle between a and b axis away from 90 • (corresponding to some shear strain) to change the lattice from orthorhombic to monoclinic, the chain will lose (part of) its protection. In this case, the distortion breaks M x but still preserves M z and P, thus the Dirac loop on the k z = π plane is still protected [Fig. 4(c)], whereas the loop on the k x = π plane and the Dirac point on T-Y are removed. These are confirmed by the DFT calculation.
Nodal loops could feature topological drumhead-like surface states [16]. We find similar phenomena for the Dirac chain here. For example, on the (001) surface, the projected loop L 2 is centered around X point, around which one indeed observes drumhead-like surface states emanating from the projected bulk band-crossing point [ Fig. 5(a,b)]. The pair of isolated Dirac points on T-Y also generates surface Fermi arcs. As shown in Fig. 5(c) for the (010) surface, the arcs connect the surface-projections of the bulk Dirac points, similar to the Dirac semimetals Na 3 Bi and Cd 3 As 2 [6,7].
A few remarks are in order before closing. First, to our knowledge, this is the first time that an hourglass Dirac chain metal is reported for a real material. The Dirac chain revealed in ReO 2 represents an essential bandcrossing: it is robust against SOC and dictated by the crystalline symmetry. Compared with the Weyl chain studied in Ref. 47, the Dirac chain here acquires additional degeneracy due to the presence of an inversion center. This important difference poses more stringent condition regarding the symmetry protection. Indeed, if there were no such degeneracy in the current case, the (missing) loop on the k y = π plane would be well protected.
The Dirac chain, the hourglass dispersion, and the surface states are close to the Fermi level. They could be detected in spectroscopic experiment such as ARPES [45]. Since a nodal loop has intrinsic anisotropy, the magnetic response could be very different for different magnetic field orientations. It has been shown that for essential nodal loops, there could also be pronounced anomaly in the longitudinal magnetotransport when the fields are aligned in the loop plane [47].
Finally, it has been argued that drumhead surface states could lead to huge surface density of states, which may provide a route towards high-temperature superconductivity [52]. In this aspect, nodal chain metals could be ideal since the orthogonal loops dictates the presence of drumhead surface states on multiple surfaces. The size of the chain and hence the surface density of states can be tuned, e.g., by strain, which could offer an additional control of the possible phase transition.