## The hybrid-order topology of weak topological insulators

Sander H. Kooi, 1 Guido van Miert, 1 and Carmine Ortix 1, 2

1 Institute for Theoretical Physics, Center for Extreme Matter and Emergent Phenomena,

Utrecht University, Princetonplein 5, 3584 CC Utrecht, the Netherlands

2 Dipartimento di Fisica 'E. R. Caianiello', Universit` di Salerno I-84084 Fisciano (Salerno), Italy a

(Dated: August 4, 2020)

We consider weak topological insulators with a twofold rotation symmetry around their 'dark' direction, and show that these systems can be endowed with the topological crystalline structure of a higher-order topological insulator protected by rotation symmetry. These hybrid-order weak topological insulators display surface Dirac cones on all surfaces. Translational symmetry breaking perturbations gap the Dirac cones on the side surfaces leaving anomalous helical hinge modes behind. We also prove that the existence of this topological phase comes about due to a novel crystalline topological invariant of quantum spin-Hall insulators that can neither be revealed by symmetry indicators nor using Wilson loop invariants. Considering the minimal symmetry requirements, we anticipate that our findings could apply to a large number of weak topological insulators.

Introduction The essence of a free-fermion topological insulator is that it cannot be adiabatically deformed to a trivial atomic insulator, whose nature can be understood considering electrons as localized point particles. Put differently, topological insulators do not admit a representation in terms of exponentially localized Wannier functions (WFs). This obstruction to Wannier representability is, in turn, reflected in the presence of anomalous gapless boundary modes. Examples include the chiral (helical) edge modes in quantum (spin) Hall insulators [1-6], as well as the surface Dirac cones of three-dimensional topological insulators (TI) [7]. In crystalline systems with an additional set of spatial symmetries, additional topological phases can arise [8]. These topological crystalline insulators (TCI) cannot be represented in terms of WFs respecting the spatial symmetries of the system, and feature, by the bulk-boundary correspondence, anomalous surface states violating a stronger version of the fermion doubling theorem [9] on surfaces that are left invariant under the protecting symmetry. Mirror Chern insulators [10, 11], for instance, are characterized by the presence of gapless surface Dirac cones pinned to mirror planes. Similarly, higher-order topological insulators (HOTI) [12-15] feature anomalous gapless one-dimensional modes at the hinges connecting two surfaces related by the protecting crystalline symmetry [16].

however, certain topological crystalline phases that are neither characterizable by symmetry indicators, nor by the tenfold way [32]. In two-dimensional systems these phases have recently started to be classified [33].

The topologies related to the internal and spatial symmetries do not necessarily exclude each other and can also coexist. This occurs, for instance in different 'dual' topological materials [17-20], which have the topological structure of both a weak TI and a mirror Chern insulator. Likewise, it has been recently proposed that certain topological superconductors can concomitantly feature both surface cones and Majorana hinge modes [21-23]. In all these systems, the topological crystalline structure can be diagnosed using the spatial symmetry content of the electronic bands [24-28] while the topology due to the internal symmetry is uniquely determined by the 'tenfold-way' invariants [29-31]. There exist,

The question that immediately arises is whether crystalline topologies without symmetry indicators can be embedded in a topological non-trivial insulating phase protected by an internal symmetry. In this work, we provide an affirmative answer by showing that twodimensional topological insulators in the wallpaper group p 2 - where time-reversal symmetry (TRS) guarantees the complete absence of symmetry indicators - can be characterized by a set of three crystalline topological Z 2 indices. They correspond to two quantized partial Berry phases [34, 35] and one additional novel topological index that cannot be diagnosed even from the Wilson loop. We subsequently use this new invariant to show that weak TIs possessing a twofold rotation symmetry around the [ ν , ν 1 2 , ν 3 ] direction, ν 1 2 3 , , indicating the so-called weak topological indices [36], can be in a non-trivial topological crystalline phase. It is characterized by the presence of anomalous unpinned Dirac cones at the surfaces whose Miller indices (modulo 2) are identical to the weak topological indices, i.e. the so-called 'dark' surfaces of weak TIs where surface Dirac cones protected by TRS are absent. This topological crystalline phase corresponds to a form of hybrid-order topology since the system can be switched to a HOTI with helical hinge modes using translational symmetry breaking perturbations.

Crystalline topological invariants in quantum spin-Hall insulators We start out by developing a scheme that is able to capture the full crystalline topology of quantum spin-Hall insulators (QSHI) in systems with a twofold rotation symmetry C 2 . To do so, we first recall that for atomic insulating phases, the crystalline topology is fully determined by the gauge-invariant charge centers [9, 33, 37] of time-reversal symmetric Wannier functions that respect the symmetries of the crystal. The construction of such symmetric Kramers pairs of Wannier

Figure 1. (a) Schematic drawing of the Brillouin zone of C 2 symmetric crystal, spanned by reciprocal lattice vectors G 1 and G 2 . The partial polarizations are calculated along the green and blue line, and a typical Wilson loop contour is shown in red. (b) Wilson loop spectrum of a QSHI. The winding reflects the topological non-trivial nature of the insulating phase. The quantized value of the Wilson loop for k 1 = 0 corresponds to the quantized partial polarizations γ I 2 .

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00879v2-figure-97104fd3-a990-45fb-9706-3a13d2ee3722.jpg)

functions requires the construction of two time-reversed channels [5] of Bloch waves | Ψ I,II n ( q ) 〉 that are separately C 2 symmetric, where n is a band index running from one to N / F 2 and N F the total number of occupied bands. The Bloch waves | Ψ I,II n ( q ) 〉 need not be individual eigenstates of the Hamiltonian, but are still basis states spanning the eigenspace corresponding to the N F occupied bands. Importantly, the construction of symmetric Wannier functions requires a smooth, periodic, and symmetric gauge for the | Ψ I,II n ( q ) 〉 Bloch waves. Since we want to study crystalline topology in non Wannier representable QSHI, we relax these constraints on the gauge by demanding its smoothness, periodicity and symmetry modulo a U ( N / F 2) ⊗ U ( N / F 2) gauge degree of freedom, with these two residual gauges acting in the two time-reversed and C 2 symmetric channels. In other words, we require a smooth, periodic and symmetric set of projectors ρ I ( II ) ( q ) = ∑ n | Ψ I ( II ) n ( q ) 〉 〈 Ψ I ( II ) n ( q ) . | In the Supplemental Material we show how to construct such a gauge assuming for simplicity there are no degeneracies in the band structure other than those required by time-reversal. Since within each sector we have not demanded a continuous gauge, it follows that the channels described by the Bloch waves | Ψ I,II n ( q ) 〉 can be characterized by non-vanishing but opposite Chern numbers C I,II . Furthermore, the twofold rotation symmetry endows the two channels with Z indices that correspond to the multiplicities of the rotation eigenvalues m I ± i ≡ m II ∓ i at the high-symmetry points in the Brillouin zone (BZ), i.e. m = Γ , X, Y, M [see Fig. 1(a)].

We will now show that these integer crystalline indices and the Chern numbers of the channels can be used to construct four Z 2 invariants that fully characterize the topology of C 2 and time-reversal symmetric insulators in two dimensions. Two Z 2 invariants can be immediately identified in the quantized partial polarizations [34] on the C 2 symmetric lines of the BZ k 1 2 , ≡ 0. They correspond to the centers of charge of one-dimensional hybrid Wannier functions and are diagnosed [33] by the Wilson loop spectra ν k ( 1 2 , ) [see Fig. 1(b)]. These quantized partial polarizations can be expressed in terms of the crystalline indices m I ± i as γ I 1(2) ≡ [ Γ I i + X I i ( Y I i ) ] mod 2 (see the Supplemental Material). The third Z 2 invariant corresponds to the Fu-Kane-Mele invariant that characterizes QSHI and can be expressed (see the Supplemental Material) in terms of the crystalline indices as ν FKM = ( Γ I i + X I i + M I i + Y I i ) mod 2. To define a fourth Z 2 invariant, notice that the additional combination of eigenvalues ν I 1 d = Γ ( I -i -X I -i -Y I -i + M I -i ) / 2 mod 2 is linearly independent from the previously defined Z 2 indices. For an atomic insulating phase, ν I 1 d corresponds to the parity of the time-reversed pairs of symmetric Wannier functions centered at the corner of the unit cell with coordinates 1 d = 1 2 1 2 . The fact that { / , / } in a QSHI the two time-reversed channels I, II are characterized by an odd Chern number immediately yields a semi-integer value ν I 1 d = ± 1 / 2. However, and this is key, we can still define a Z 2 number reading

$$\underset { \text{li-} } { \mathbb { d } } \quad \gamma _ { 3 } ^ { I } = & \frac { 1 } { 2 } \left [ C ^ { I } + \left ( \Gamma _ { - i } ^ { I } - X _ { - i } ^ { I } - Y _ { - i } ^ { I } + M _ { - i } ^ { I } \right ) \, \right ] \mod 2.$$

Being independent of the partial polarizations, this new integer cannot be diagnosed by the Wilson loop spectrum but still represents a well-defined and gauge-invariant crystalline topological number. In fact, γ I 3 is manifestly gauge-invariant under intrachannel U ( N / F 2) transformations since it is made out of a Chern number and the twofold rotation symmetry eigenvalues. Furthermore, γ I 3 is also invariant under interchannel gauge transformations, which correspond to the swapping of the channels ( I ↔ II ) for isolated pairs of bands. These transformations concomitantly change the sign of the Chern numbers of the channels and the multiplicities of the C 2 symmetry eigenvalues, and therefore do not change γ I 3 . We have thus identified three gauge-invariant Z 2 crystalline topological indices, which together with the Fu-KaneMele invariant yield a Z 4 2 classification in agreement with a recent K theory study [38].

We finally emphasize that the gauge-invariant γ I 3 is different in nature from the 'spin Chern numbers' existing in systems with a mirror symmetry M z . In this situation, the two time-reversed and C 2 symmetric channels I, II can be taken to be the spin eigenstates |↑〉 , |↓〉 , such that C I ≡ C ↑ . However, this does not determine the value of γ I 3 , as the spin Chern number does not determine ν I 1 d . Thus one can find both γ I 3 = 0 1 for the , same spin Chern number.

Hybrid-order weak TIs Next, we exploit the existence of the novel crystalline topology of γ I 3 in threedimensional bulk crystals with a C 2 z rotational symmetry. To do so, let us consider the three-dimensional Bril-

Figure 2. Schematic of a hybrid-order weak topological insulator. At the top and bottom surfaces left invariant under the protecting C 2 z symmetry a single pair of surface Dirac cones exist. On the side surfaces an even number of Dirac cones pinned to time-reversal invariant surface momenta are mandated by the weak topological invariants. When breaking C 2 symmetry, the topological crystalline surface Dirac cones at the top and bottom surfaces can be gapped out leaving these surfaces completely dark. By breaking the translational symmetry, i.e. doubling the unit cell, the time-reversal symmetry protected Dirac cones gap out, and the topological crystalline Dirac cones are then connected by helical hinge states.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00879v2-figure-2716a416-e09e-4a7c-85f3-98f66ce5e020.jpg)

louin zone of our time-reversal invariant system as a collection of two-dimensional momentum cuts parametrized by the momentum k z parallel to the twofold rotation axis. At the time-reversal invariant two-dimensional planes k z = 0 , π we consider the system to be a topological non-trivial QSHI. As a result, the bulk three-dimensional crystal will be a three-dimensional topological insulator of the weak class. In principle, we could choose the two Z 2 topological crystalline indices corresponding to the quantized partial polarization of the k z = 0 , π QSHI to be different. This, however, would imply that in the triad of 'weak' topological invariants [36] ( ν , ν 1 2 , ν 3 ), ν 1 and/or ν 2 are different from zero. Hence, the three-dimensional system would feature an even number of surface Dirac cones protected by time-reversal at the (001) and (00 ¯ ) 1 surfaces that are left invariant under the C 2 z rotation symmetry. As a result, any physical consequence of the crystalline topology cannot manifest itself: it would be completely obscured by the internal, time-reversal, symmetry topology.

However, we can choose the two Z 2 topological crystalline invariants at the time-reversal invariant planes to be equal, thus constraining the weak invariants to be (0 0 1). , , The time-reversal symmetric topology now guarantees the existence of an even number of massless Dirac cones appearing at time-reversal invariant (100) and (010) surface momenta [c.f. Fig. 2], while the C 2 z invariant (001) surfaces are completely gapped. A nontrivial crystalline topology, which can thus only arise from a difference in γ I 3 at the k z = 0 , π planes, will then be in full force and lead to the appearance of a single pair of surface Dirac cones [c.f. Fig. 2] at unpinned sur- face momenta related by the twofold rotation symmetry. This pair of surface Dirac cones realizes the rotational anomaly discussed in Ref. 9, and can be only removed by breaking the protecting C 2 z and/or Θ symmetry [c.f. Fig. 2]. We point out that the existence of this rotation anomaly cannot be diagnosed by considering the flow of gauge-invariant Wannier centers between the k z = 0 , π planes as in Ref. 9. This is because at k z = 0 , π our system is a topological insulator, and therefore cannot be represented in terms of localized Wannier functions. The appearance of the unpinned Dirac surface cones is instead detected by considering the k z -directed Wilson loop [see the Supplemental Material] in agreement with Ref. [39], although the stability of the surface Dirac cones cannot be inferred from the Wilson loop that consequently cannot be used to derive a 'topological index'. We dub this new three-dimensional insulating phase a hybrid-order weak topological insulator: it is by itself a first-order topological insulator in d = 3 dimensions with d -1 gapless boundary modes, but it can be switched using unit cell doublings in the ˆ-direction, and thus without breakz ing any protecting symmetry, to a second-order topological crystalline insulator with anomalous gapless hinge modes [c.f. Fig. 2] and C 2 rotation anomaly [9], and reminiscent of the surface cones one predicted to appear in α -Bi Br 4 4 and a family of Zintl compounds [40, 41].

Stacked Kane-Mele model Having established the existence of the hybrid-order weak topological insulator, we now present an explicit model based on stacked KaneMele systems realizing this phase. Let us consider a tightbinding model for spin-1 / 2 electrons on AA stacked honeycomb lattices. In momentum space the Bloch Hamiltonian can be written as:

$$\underset { \ell } { \text{e} } \text{of} \\ / \text{or} \\ \overset { \text{e} } { \text{of} } \text{$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$ \quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\ \text{of}$$

where the τ i 's and s i 's are the Pauli matrices acting in sublattice and spin space respectively. The first two terms in the Hamiltonian above correspond to intralayer spin-independent nearest-neighbor hopping processes, and the corresponding coefficients are d 1 ( k ) = -t [1 + cos x 1 +cos x 2 ] and d 2 ( k ) = -t [sin x 1 +sin x 2 ]. Here we have introduced the hopping amplitude t while x 1 2 , = k · a 1 2 , , a 1 2 , being the Bravais lattice vectors. The third term in the Hamiltonian Eq. 1 corresponds to spin-orbit interaction which involves intralayer spindependent second-neighbor hopping. We take the corresponding coefficient d 5 ( k ) = 2 t 2 sin ( x 1 ), with t 2 the hopping strength, thus explicitly breaking the threefold rotation symmetry. Finally, the last term in the Hamiltonian involves interlayer spin-dependent hopping amplitudes and the corresponding coefficient reads d 4 ( k ) = -2 t 3 sin ( k z ). We introduce this term to explicitly break the effective 'in-plane' time-reversal symmetry [42] to allow for the possibility of a change of (crystalline) topology in the two time-reversal symmetric planes k z = 0 , π .

Figure 3. (a) Bulk bandstructure of the stacked Kane-Mele model Eq. 1 by further accounting for a Rashba spin-orbit coupling term of strength λ/t = 0 1. . The strength of the modified intralayer spin-orbit coupling term has been fixed to t 2 /t = 0 7, whereas the interlayer spin-orbit coupling strength . has been fixed to t 3 /t = 0 4. (b) Energy spectrum in a slab ge-. ometry with open boundary conditions along the ˆ direction. y The (010) and (0 ¯ 0) surfaces exhibit an even number of Dirac 1 cones pinned at time-reversal invariant surface momenta as required by the weak topological invariants. (c) Surface energy spectrum along the stacking ˆ direction. z There are two pairs of surface Dirac cones localized at the (001) and the (00 ¯ ) 1 surface. The Dirac points are found at unpinned surface momenta related by the C 2 symmetry. The zoom-in (d) shows that the Dirac cones at opposite surface are located at different momenta due to the lack of inversion symmetry. All energies have been measured in unit of the hopping strength t .

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00879v2-figure-e373a297-3bcf-4350-a965-cd1d474a811a.jpg)

Figure 4. (a) Bulk band structure of the stacked Kane-Mele model with a translational breaking perturbation. The parameter set is the same as in Fig. 3. Moreover the translational symmetry breaking parameters have been fixed to glyph[epsilon1]/t = 0 1 and . δ/t = 0 2. (b) Surface energy spectrum showing . the gapping of the time-reversal symmetry protected Dirac cones. (c) Surface energy spectrum along the stacking direction that still feature the C 2 -protected Dirac cones at unpinned surface momenta. (d) Energy spectrum in a ribbon geometry with periodic boundary conditions only along the stacking direction. Within the surface energy gap we find gapless anomalous helical hinge modes, colored in red.

Since the Hamiltonian Eq. 1 preserves bulk inversion symmetry, we can immediately obtain the strong and weak topological indices and thus obtain ( ν 0 ; ν , ν 1 2 , ν 3 ) = (0; 0 , 0 1). , In this form, however, Eq. 1 does not model a hybrid-order weak topological insulator: it can be adiabatically connected to a stack of uncoupled QSHI and consequently its (001) surface does not feature gapless modes. To endow the system with a non-trivial crystalline topology we instead modify the intralayer spinorbit coupling as d 5 ( k ) → cos( k z ) d 5 ( k ). This modification keeps the strong and weak topological indices intact but changes the crystalline topology of the system. Note that also the inversion eigenvalues remain unchanged, thus implying that the hybrid-order phase cannot be diagnosed by inversion symmetry indicators.

To show this, we have computed the bulk bandstructure [see Fig. 3(a)] and the surface energy spectra [see Fig. 3(b),(c),(d)] of this modified model by further accounting for an intralayer Rashba spin-orbit coupling term [43] of strength λ that explicitly breaks inversion symmetry. At the side surfaces we observe the conventional surface Dirac cones of a weak topological insulators [c.f. Fig. 3(b)]. More importantly, diagonalization of the Hamiltonian with open boundary conditions along the stacking direction [c.f. Fig. 3(c)] reveals the presence of two C 2 symmetry protected surface Dirac cones thus verifying that our model realizes a hybrid-order weak topological insulator. Note that the pairs of Dirac cones at the (001) and (00 ¯ ) surface are found at different surface 1 momenta in agreement with the lack of inversion symmetry.

We have also verified that our model can be switched to a higher-order topological insulator by suitable translational symmetry breaking perturbations. Specifically we have introduced an interlayer staggered chemical potential of strength glyph[epsilon1] that provides the required doubling of the unit cell and further introduced an interlayer coupling in the enlarged unit cell of the form -δτ z s x . Fig. 4(a) shows the corresponding bulk bandstructure that is still characterized by a substantial gap. At the (010) [(0 ¯ 0)] 1 surface the time-reversal symmetry protected Dirac cones are gapped out [see Fig. 4(b)] while the twofold rotation symmetry-protected Dirac cones at the (001) [(00 ¯ )] sur1 face are preserved [see Fig. 4(c)]. Notice that Dirac cone

pairs localized at opposite surface are connected by helical hinge states [see Fig. 4(d)] as expected for a helical higher-order topological insulator protected by a twofold rotation symmetry.

Conclusions To sum up, we have shown that weak topological insulators with an additional twofold rotation symmetry around the dark direction can feature a pair of Dirac cones on their dark surfaces, which are protected by the rotation symmetry. This hybrid-order weak topological insulator can be turned into a higher-order topological insulator with protected helical hinge modes by translational symmetry breaking perturbations. We have shown that the existence of such a topological phase comes about due to a third Z 2 topological invariant characterizing quantum spin-Hall insulators in C 2 -symmetric crystals, that can be read off neither from symmetry indicators nor from the properties of the Wilson loop spectrum. Considering the minimal symmetry requirements and the fact that the C 2 protected surface Dirac cones appear at unpinned points in the surface Brillouin zone, we anticipate that our findings could apply to a large number of weak topological insulators.

C.O. acknowledges support from a VIDI grant (Project 680-47-543) financed by the Netherlands Organization for Scientific Research (NWO). This work is part of the research programme of the Foundation for Fundamental Research on Matter (FOM), which is part of the Netherlands Organization for Scientific Research (NWO). S.K. acknowledges support from a NWO-Graduate Program grant.

- [1] K. v. Klitzing, G. Dorda, and M. Pepper, Phys. Rev. Lett. 45 , 494 (1980).
- [2] D. J. Thouless, M. Kohmoto, M. P. Nightingale, and M. den Nijs, Phys. Rev. Lett. 49 , 405 (1982).
- [3] C. L. Kane and E. J. Mele, Phys. Rev. Lett. 95 , 226801 (2005).
- [4] C. L. Kane and E. J. Mele, Phys. Rev. Lett. 95 , 146802 (2005).
- [5] L. Fu and C. L. Kane, Phys. Rev. B 74 , 195312 (2006).
- [6] B. A. Bernevig, T. L. Hughes, and S.-C. Zhang, Science 314 , 1757 (2006).
- [7] L. Fu, C. L. Kane, and E. J. Mele, Phys. Rev. Lett. 98 , 106803 (2007).
- [8] L. Fu, Phys. Rev. Lett. 106 , 106802 (2011).
- [9] C. Fang and L. Fu, Science Advances 5 , eaat2374 (2019).
- [10] T. H. Hsieh, H. Lin, J. Liu, W. Duan, A. Bansil, and L. Fu, Nat. Comm. 3 , 982 (2012).
- [11] T. H. Hsieh, J. Liu, and L. Fu, Phys. Rev. B 90 , 081112(R) (2014).
- [12] W. A. Benalcazar, B. A. Bernevig, and T. L. Hughes, Science 357 , 61 (2017).
- [13] F. Schindler, A. M. Cook, M. G. Vergniory, Z. Wang, S. S. Parkin, B. A. Bernevig, and T. Neupert, Sci. Adv. 4 , eaat0346 (2018).
- [14] F. Schindler, Z. Wang, M. G. Vergniory, A. M. Cook, A. Murani, S. Sengupta, A. Y. Kasumov, R. Deblock, S. Jeon, I. Drozdov, et al. , Nat. Phys. 14 , 918 (2018).
- [15] G. van Miert and C. Ortix, Phys. Rev. B 98 , 081110(R) (2018).
- [16] S. H. Kooi, G. van Miert, and C. Ortix, Phys. Rev. B 98 , 245102 (2018).
- [17] M. Eschbach, M. Lanius, C. Niu, E. Mglyph[suppress]y´czak, l n P. Gospodariˇ c, J. Kellner, P. Sch¨ffelgen, M. Gehlmann, u S. D¨ring, E. Neumann, o et al. , Nature Comm. 8 , 14976 (2017).
- [18] N. Avraham, A. Norris, Y. Sun, Y. Qi, L. Pan, A. Isaeva, A. Zeugner, C. Felser, B. Yan, and H. Beidenkopf, arXiv:1708.09062 (2017).
- [19] I. Rusinov, T. Menshchikova, A. Isaeva, S. Eremeev, Y. M. Koroteev, M. Vergniory, P. Echenique, and E. V. Chulkov, Sci. Rep. 6 , 20734 (2016).
- [20] J. I. Facio, S. K. Das, Y. Zhang, K. Koepernik, J. van den Brink, and I. C. Fulga, Phys. Rev. Materials 3 , 074202 (2019).
- [21] N. Bultinck, B. A. Bernevig, and M. P. Zaletel, Phys. Rev. B 99 , 125149 (2019).
- [22] S. A. A. Ghorashi, X. Hu, T. L. Hughes, and E. Rossi, Phys. Rev. B 100 , 020509(R) (2019).
- [23] S. A. A. Ghorashi, T. L. Hughes, and E. Rossi, arXiv preprint arXiv:1909.10536 (2019).
- [24] J. Zak, Phys. Rev. B 26 , 3010 (1982).
- [25] B. Bradlyn, L. Elcoro, J. Cano, M. G. Vergniory, Z. Wang, C. Felser, M. I. Aroyo, and B. A. Bernevig, Nature 547 , 298 (2017).
- [26] J. Cano, B. Bradlyn, Z. Wang, L. Elcoro, M. G. Vergniory, C. Felser, M. I. Aroyo, and B. A. Bernevig, Phys. Rev. Lett. 120 , 266401 (2018).
- [27] E. Khalaf, H. C. Po, A. Vishwanath, and H. Watanabe, Phys. Rev. X 8 , 031070 (2018).
- [28] H. C. Po, A. Vishwanath, and H. Watanabe, Nature Communications 8 , 50 (2017).
- [29] A. Altland and M. R. Zirnbauer, Phys. Rev. B 55 , 1142 (1997).
- [30] A. Kitaev, in AIP Conference Proceedings , Vol. 1134 (AIP, 2009) pp. 22-30.
- [31] J. Kruthoff, J. de Boer, J. van Wezel, C. L. Kane, and R.-J. Slager, Phys. Rev. X 7 , 041069 (2017).
- [32] Z. Song, L. Elcoro, N. Regnault, and B. A. Bernevig, arXiv e-prints , arXiv:1905.03262 (2019).
- [33] S. H. Kooi, G. van Miert, and C. Ortix, Phys. Rev. B 100 , 115160 (2019).
- [34] A. Lau, J. van den Brink, and C. Ortix, Phys. Rev. B 94 , 165164 (2016).
- [35] G. van Miert and C. Ortix, Phys. Rev. B 96 , 235130 (2017).
- [36] L. Fu and C. L. Kane, Phys. Rev. B 76 , 045302 (2007).
- [37] Z. Song, Z. Fang, and C. Fang, Phys. Rev. Lett. 119 , 246402 (2017).
- [38] J. Kruthoff, J. de Boer, and J. van Wezel, Phys. Rev. B 100 , 075116 (2019).
- [39] L. Fidkowski, T. Jackson, and I. Klich, Phys. Rev. Lett. 107 , 036601 (2011).
- [40] C.-H. Hsu, X. Zhou, Q. Ma, N. Gedik, A. Bansil, V. M. Pereira, H. Lin, L. Fu, S.-Y. Xu, and T.-R. Chang, 2D Mat. 6 , 031004 (2019).
- [41] T. Zhang, C. Yue, T. Zhang, S. Nie, Z. Wang, C. Fang, H. Weng, and Z. Fang, Phys. Rev. Research 1 , 012001 (2019).

[42] A. Lau, C. Ortix, and J. van den Brink, Phys. Rev. B 91 , 085106 (2015).

$$\text{ev. } B \quad [sin(k _ { 1 } ) \sqrt { 3 } / 2 ] \tau _ { x } \otimes s _ { y }. \\ \quad [44 ] \text{ } T. \, L. \, Hughes, E. \, Prodan, \text{ and } B. \, A. \, Bernevig, Physical }.$$

[43] The term is [ -1 / 2 + cos( k 2 ) -cos( k 1 )] τ y ⊗ s x + [sin( k 2 ) -sin( k 1 ) / 2] τ x ⊗ s x +[ √ 3 2(cos( / k 1 ) -1)] τ y ⊗ σ y +

T. L. Hughes, E. Prodan, Review B 83 , 245132 (2011)

.

## SUPPLEMENTAL MATERIAL FOR 'THE HYBRID-ORDER TOPOLOGY OF WEAK TOPOLOGICAL INSULATORS'

## Constructing C 2 -symmetric channels for pairs of isolated Kramers pairs

Our goal is to numerically construct, for systems with C 2 and time-reversal symmetry, a gauge that divides our system in two channels that are by themselves C 2 symmetric, and are related to each other by time-reversal Θ, and consequently C 2 Θ symmetry. Since we are assuming that there are no other degeneracies than those required by time-reversal symmetry, it is sufficient to consider a single Kramers pair of bands. We then want to find a gauge such that we have two states ψ I k and ψ II k that satisfy

$$\mathcal { C } _ { 2 } \Theta \psi _ { k } ^ { I } = \psi _ { k } ^ { I I },$$

$$\mathcal { C } _ { 2 } \psi _ { k } ^ { I, I I } = \psi _ { - k } ^ { I, I I }.$$

To construct such a gauge, let us first see how to construct locally, at each k -point, a gauge that satisfies Eq. (2). To do so we start by numerically diagonalizing the Hamiltonian to find the two occupied eigenstates Ψ I,II k . We then calculate the unitary matrix M ij = Ψ i † k C 2 θ Ψ j k and diagonalize it by a transformation ψ i k = U ij k Ψ where j k , U ij k consists of the eigenvectors of M ij . The matrix M then has the form

$$M = \begin{pmatrix} e ^ { i \rho _ { 1 } } & 0 \\ 0 & e ^ { i \rho _ { 2 } } \end{pmatrix}.$$

After a further unitary transformation ψ 1 2 , k → e iρ 1 2 , / 2 ψ 1 2 , k , M will be the identity matrix. A final unitary transformation

$$\begin{pmatrix} \psi ^ { I } _ { \tilde { I } _ { I } } \\ \psi ^ { I I } _ { k } \end{pmatrix} \rightarrow \frac { 1 } { \sqrt { 2 } } \begin{pmatrix} - 1 & 1 \\ 1 & 1 \end{pmatrix} \begin{pmatrix} \psi ^ { I } _ { \tilde { I } _ { I } } \\ \psi ^ { I I } _ { k } \end{pmatrix},$$

then ensures M is completely off-diagonal, and thus we have found states ψ I k and ψ II k that satisfy Eq. (2).

In order to construct a gauge satisfying Eqs. (2) and (3) across the Brillouin zone, let us start by noting that for an isolated Kramers pair of bands, the local gauge freedom is U (1) ⊗ SU (2). That is, we are free to do a rotation between the basis states, and can then change each basis state by a U (1) phase. Now the conditions of Eqs. (2) and (3) completely fix the SU (2) part of the gauge freedom, up to the unitary operation of relabeling the states ψ I k ←→ ψ II k . There will thus be a remaining U (1) degree of freedom within each C 2 symmetric channel.

We can now start at the Γ point in the Brillouin zone at k = (0 0), and follow the steps described to make the , matrix M off-diagonal. Since Γ is a time-reversal invariant point, Kramers theorem guarentees that the time-reversal operator θ is off-diagonal. Rotating such that C θ 2 is off-diagonal thus also ensures that this basis is C 2 diagonal.

Once we have found the states at the Γ point, we take a point close to it at k = ( glyph[epsilon1], 0) for a small glyph[epsilon1] and repeat the procedure to make M off-diagonal. We then consider the overlap matrix ψ i (0 0) , ψ j ( glyph[epsilon1], 0) , which will either be almost completely diagonal, or almost completely off-diagonal. If it is off-diagonal, this means that we need to switch the labels I, II by a unitary transformation

$$\binom { \psi _ { k } ^ { I } } { \psi _ { k } ^ { I I } } \rightarrow \binom { 0 } { 1 } \binom { \psi _ { k } ^ { I } } { \psi _ { k } ^ { I I } } \,.$$

We then define ψ i -k = C 2 ψ i k . Continuing in this way along half the Brillouin zone also ensures that the second condition Eq. (2) is satisfied.

The gauge constructed in this way fixes the SU (2) part of the gauge degrees of freedom, since we ensured that the overlap matrix is diagonal. Within each channel, there is still a U (1) phase degree of freedom which we have not fixed, and for which there is an obstruction to smoothness if the band has a non-zero Chern number. For our purposes we do not need to fix this phase, as neither the eigenvalues nor the Chern number can change under this U (1) gauge transformation

## Topological invariants of systems with two-fold rotational symmetry

Using the procedure of the previous section, we can divide the occupied bands of a system into two channels that are tow-fold rotation symmetric, and have broken time-reversal symmetry. In this section we show how one can compute the topology of the system once such a gauge has been found.

First consider the partial polarizations, which are quantized along the high-symmetry lines of the BZ by C 2 and time-reversal symmetry. Using the constructed gauge, the two C 2 -symmetric channels give us an easy way to compute the partial polarizations. Since the two channels are related by time-reversal symmetry, the partial polarization is equal to the polarization in one channel. This means we have to calculate the Berry phase of one channel along highsymmetry lines in the BZ. On these contours, each channel is effectively a one-dimensional system with C 2 symmetry. For such systems, the Berry phase γ is quantized and can directly be expressed in terms of its C 2 -eigenvalues as [44]

$$\gamma = & ( - i ) \prod _ { i } ^ { N _ { F } } \log [ - \xi _ { i } ( 0 ) \xi _ { i } ( \pi ) ] \mod 2 \pi,$$

where ξ i ( k ) is the C 2 eigenvalue of band i at momentum k . Since the eigenvalues can be either ± i , the term inside the logarithm can either be ± and hence each term inside the product is either 0 or π . This allows us to rewrite

$$\gamma = & \sum _ { i } ^ { N _ { F } } - i \log [ - \xi _ { i } ( 0 ) \xi _ { i } ( \pi ) ] \mod 2 \pi \\ = & \sum _ { i } ^ { N _ { F } } - i ( \log [ i \xi _ { i } ( 0 ) ] + \log [ i \xi _ { i } ( \pi ) ] ) \mod 2 \pi \\ = & - i \sum _ { i } ^ { N _ { F } } \log [ i \xi _ { i } ( 0 ) ] - i \sum _ { i } ^ { N _ { F } } \log [ i \xi _ { i } ( \pi ) ] ) \mod 2 \pi \\ = & \pi ( m _ { i } ^ { k = 0 } + m _ { i } ^ { k = \pi } ) \mod 2 \pi, \\ \text{the multiplicity of eigenvalue} \ i \text{ at } k = 0 \text{ and the last equality} \text{ to}$$

Where m k =0 i denotes the multiplicity of eigenvalue i at k = 0 and the last equality follows since whenever an eigenvalue is -i , the log evaluates to zero.

Applying this to the contours Γ -X and Γ -Y we find γ I 1(2) ≡ [ Γ I i + X I i ( Y I i ) ] mod 2 as presented in the main text.

The Fu-Kane-Mele invariant is equal to the Chern number in one of the channels modulo 2. This Chern number can in turn be expressed as the difference of the berry phases of the contours Γ -Y and X -M [38]. Hence the FKM-invariant can be expressed as

$$\nu _ { \text{FKM} } = & ( \Gamma _ { i } ^ { I } + X _ { i } ^ { I } + Y _ { i } ^ { I } + M _ { i } ^ { I } ) \mod 2$$

So far we have found three invariants constructed out of four eigenvalue multiplicities. For atomic insulators, a fourth linearly independent crystalline invariant is given by [15]

$$\nu _ { 1 d } = & \frac { 1 } { 2 } ( \Gamma _ { - i } - X _ { - i } - Y _ { - i } + M _ { - i } ),$$

where m -i is the multiplicity of eigenvalue -i at high-symmetry point m = Γ , X, Y, M .

This invariant is related to the number of Wannier functions localized at Wyckoff position 1 d in real space by [15]

$$\nu _ { 1 d } = - \, N _ { 1 d ; i } + N _ { 1 d ; - i }$$

where N 1 d i ; is the number of symmetric Wannier functions localized at Wyckoff position 1 d with two-fold rotation eigenvalue i . When considering non-atomic insulators, one cannot find localized symmetric Wannier functions, however

we can still calculate ν 1 d . For systems with an odd Chern number ν 1 d will be half-integer, since there will be an odd number of i eigenvalues, and an odd number of -i eigenvalues. Hence we can construct an invariant that is always integer by the combination

$$\gamma _ { 3 } = ( \mathcal { C } / 2 + \nu _ { 1 d } ) \mod 2.$$

## Sums of quantum spin Hall insulators

Here we explicitly show how adding together and coupling two QSHI can result in either a fragile topological phase or a trivial phase. To this end we consider a Kane-Mele model,

$$\mathcal { H } _ { K M } \left ( t, t _ { 2 } \right ) & = - t \left [ 1 + \cos \left ( k _ { 1 } \right ) + \cos \left ( k _ { 2 } \right ) \right ] \tau _ { x } \otimes s _ { 0 } \\ & - t \left [ \sin \left ( k _ { 1 } \right ) + \sin \left ( k _ { 2 } \right ) \right ] \tau _ { y } \otimes s _ { 0 } + 2 t _ { 2 } \sin \left ( k _ { 1 } \right ) \tau _ { z } \otimes s _ { z },$$

where t and t 2 are the nearest-neighbor and intrinsic spin-orbit coupling amplitudes, and τ i and s i are Pauli matrices acting in sub-lattice and spin-space. We also add a Rashba coupling

$$\mathcal { H } _ { R } = & \lambda \Big \{ [ - 1 / 2 + \cos ( k _ { 2 } ) - \cos ( k _ { 1 } ) ] \tau _ { y } \otimes s _ { x } + [ \sin ( k _ { 2 } ) - \sin ( k _ { 1 } ) / 2 ] \tau _ { x } \otimes s _ { x } + [ \sqrt { 3 } / 2 ( \cos ( k _ { 1 } ) - 1 ) ] \tau _ { y } \otimes \sigma _ { y } \\ + & [ \sin ( k _ { 1 } ) \sqrt { 3 } / 2 ] \tau _ { x } \otimes s _ { y } \Big \},$$

where λ is the Rashba amplitude. Note that the intrinsic spin-orbit coupling term only acts along k 1 such that threefold rotation symmetry is broken.

We now take two copies of this model, and couple them by adding a term -δiσ x ⊗ τ x ⊗ s y , where δ is the amplitude of the coupling and σ x is a Pauli matrix acting in 'copy' space. Let us denote the parameters of the two copies by t i and t i 2 , where i = 1 2 denotes the copy. To examine the phase that results from coupling the two QSHI, we consider , its Wilson loop spectrum along k 1 (see main text). Let us first consider the parameters of the two Kane-Mele models to be equal and such that the system is insulating. This leads to a Wilson loop spectrum as in Fig. 5(a), showing this is an atomic insulator with two Kramers pairs at unpinned points in the BZ. Changing the sign of t 2 between the copies, such that t 1 2 = -t 2 2 and keeping the sign of t the same, results in a fragile topological phase [see Fig. 5(b)]. Note that when turning off Rashba coupling, such that we restore inversion symmetry, the two copies have the same inversion eigenvalues, meaning this phase or γ 3 cannot be diagnosed by inversion. To obtain this phase, we can also instead change the sign of t between the copies such that t 1 = -t 2 , while keeping the sign of the other parameters equal. Changing both the sign of t and t 2 between the two copies again results in the trivial phase.

Finally, we note that when taking the hybrid-order topology model and adding translational-symmetry breaking perturbation (as presented in the main text), the resulting system at k z = 0 will be fragile topological insulator with Wilson loop Fig. 5(b), while at k z = π it will be a trivial insulator with Wilson loop Fig. 5(a).

## Spectral flow in the hybrid-order topological insulator

In general protected edge states can be detected by studying the Wilson loop spectrum along the momentum direction perpendicular to the surface [39]. For the hybrid-order topological insulator, the C 2 -related cones on the top and bottom surfaces can be detected by considering the Wilson loop along k z , as a function of k 1 and k 2 . This is plotted in Fig. 6 for the hybrid-order topological insulator presented in the main text. Here we see two protected gap-closing points that are related by C 2 symmetry.

## Channel decomposition of the hybrid-order topological insulator

In Table I provide the multiplicities of the C 2 eigenvalues in the two time-reversed and rotation-symmetric channels for our pristine hybrid order topological phase. Then we consider a dimerization in the z -direction that trivializes the time-reversal invariant weak topological indices but does not change the value of γ I 3 (Table II). Finally, in Table III we consider a dimerization in the x -direction, which does not change the weak topological indices, nor the Z 2 index governing the appearance of unpinned Dirac cones at the bottom and top surfaces given by [ γ I 3 ( k z = π ) -γ I 3 ( k z = 0)] mod 2.

Figure 5. (a) Wilson loop spectrum of two coupled Kane-Mele models for t 2 = t 1 , t 1 2 /t 1 = t 2 2 /t 1 = 0 7, . λ/t 1 = 0 2 and . δ/t 1 = 0 4. . (b) Wilson loop spectrum for the same model as in a), changing only t 1 2 /t 1 = -t 2 2 /t 1 = 0 7. .

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00879v2-figure-83e2ac0c-d0c0-40ed-9556-4b9b8f647443.jpg)

Figure 6. Wilson loop spectrum along k z plotted as a function of k 1 and k 2 of the hybrid-order topological insulator presented in the main text. Two C 2 related band crossings indicate gapless modes on the top and bottom surfaces.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00879v2-figure-024c9b9c-a93b-42ab-ba71-767f9f8e164c.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00879v2-figure-2fb99b80-18e3-4721-a529-ec59c0d2f1c8.jpg)

Table I. Tables showing the channel eigenvalues and channel Chern number of the two channels I and II at k z = 0 , π of the hybrid-order topological insulator.

Table II. Tables showing the channel eigenvalues and channel Chern number of the two channels I and II for the hybrid-order topological insulator with a translational-symmetry perturbation in the z -direction at k z = 0 , π .

|         | C I   | Γ I          | X I          | Y I          | M I          |
|---------|-------|--------------|--------------|--------------|--------------|
| k z = 0 | 2     | {- i, + i }  | {- i, + i }  | {- i, + i }  | {- i, + i }  |
| k z = π | 0     | { + i, + i } | { + i, + i } | { + i, + i } | {- i, - i }  |
|         | C II  | Γ II         | X II         | Y II         | M II         |
| k z = 0 | - 2   | {- i, + i }  | {- i, + i }  | {- i, + i }  | {- i, + i }  |
| k z = π | 0     | {- i, - i }  | {- i, - i }  | {- i, - i }  | { + i, + i } |

Table III. Tables showing the channel eigenvalues and channel Chern number of the two channels I and II for the hybrid-order topological insulator with a translational-symmetry perturbation in the x -direction at k z = 0 , π .

|         | C I   | Γ I          | X I          | Y I          | M I          |
|---------|-------|--------------|--------------|--------------|--------------|
| k z = 0 | 1     | { + i, + i } | {- i, + i }  | {- i, + i }  | {- i, + i }  |
| k z = π | 1     | {- i, - i }  | { + i, + i } | { + i, + i } | {- i, - i }  |
|         | C II  | Γ II         | X II         | Y II         | M II         |
| k z = 0 | - 1   | {- i, - i }  | { + i, - i } | { + i, - i } | { + i, - i } |
| k z = π | - 1   | { + i, + i } | {- i, - i }  | {- i, - i }  | { + i, + i } |
Sander H. Kooi Guido van Miert Carmine Ortix
The hybrid-order topology of weak topological insulators
We consider weak topological insulators with a twofold rotation symmetry
around the dark direction, and show that these systems can be endowed with the
topological crystalline structure of a higher-order topological insulator
protected by rotation symmetry. These hybrid-order weak topological insulators
display surface Dirac cones on all surfaces. Translational symmetry breaking
perturbations gap the Dirac cones on the side surfaces leaving anomalous
helical hinge modes behind. We also prove that the existence of this
topological phase comes about due to a hidden crystalline topological invariant
of quantum spin-Hall insulators that can neither be revealed by symmetry
indicators nor using Wilson loop invariants. Considering the minimal symmetry
requirements, we anticipate that our findings could apply to a large number of
weak topological insulators.