## Boundary-obstructed topological phases

Eslam Khalaf, 1 Wladimir A. Benalcazar, 2 Taylor L. Hughes, 3 and Raquel Queiroz 4

1 Department of Physics, Harvard University, Cambridge, MA 02138

2 Department of Physics, The Pennsylvania State University, University Park, PA 16802, USA

3 Department of Physics and Institute for Condensed Matter Theory,

University of Illinois at Urbana-Champaign, IL 61801, USA

4 Department of Condensed Matter Physics, Weizmann Institute of Science, Rehovot 7610001, Israel

(Dated: March 24, 2021)

Symmetry protected topological (SPT) phases are gapped phases of matter that cannot be deformed to a trivial phase without breaking the symmetry or closing the bulk gap. Here, we introduce a new notion of a topological obstruction that is not captured by bulk energy gap closings in periodic boundary conditions. More specifically, given a symmetric boundary termination we say two bulk Hamiltonians belong to distinct boundary obstructed topological phases (BOTPs) if they can be deformed to each other on a system with periodic boundaries, but cannot be deformed to each other in the open system without closing the gap at at least one high symmetry surface. BOTPs are not topological phases of matter in the standard sense since they are adiabatically deformable to each other on a torus but, similar to SPTs, they are associated with boundary signatures in open systems such as surface states or fractional corner charges. In contrast to SPTs, these boundary signatures are not anomalous and can be removed by symmetrically adding lower dimensional SPTs on the boundary, but they are stable as long as the spectral gap at high-symmetry edges/surfaces remains open. We show that the double-mirror quadrupole model of [Science, 357(6346), 2018] is a prototypical example of such phases, and present a detailed analysis of several aspects of boundary obstructions in this model. In addition, we introduce several threedimensional models having boundary obstructions, which are characterized either by surface states or fractional corner charges. Furthermore, we provide a complete characterization of boundary obstructed phases in terms of symmetry representations. Namely, two distinct BOTP phases correspond to equivalent band representations in the periodic system which become inequivalent upon restricting the symmetry group to that of the open system. This is used to shown that for a given open boundary, there is only one class of BOTPs which corresponds to a local representation of the symmetry of the open system and thus can be designated as the trivial phase. All other BOTP classes do not correspond to local representation of the open system and as a result necessarily exhibit a filling anomaly or gapless surface states.

## CONTENTS

| I.   | Introduction                                                            |   2 |
|------|-------------------------------------------------------------------------|-----|
| II.  | 2D quadrupole insulator and boundary obstructions                       |   3 |
| A.   | Review of the quadrupole insulator model                                |   3 |
| B.   | Real space picture                                                      |   5 |
| C.   | Wannier bands and edge spectrum                                         |   7 |
| D.   | Corner charge and filling anomaly                                       |   9 |
| E.   | Symmetry representations                                                |  12 |
|      | 1. Equivalence between band representations in the bulk periodic system |  12 |
|      | 2. Symmetry representations of the open system                          |  13 |
|      | 3. Wannier band representations                                         |  14 |
| III. | Boundary obstructions: Generalities                                     |  14 |
| A.   | General definition                                                      |  14 |
|      | 1. Boundary termination                                                 |  15 |
|      | 2. Boundary obstructed phases                                           |  16 |
|      | 3. Relationship to higher-order topological insulators                  |  17 |
| B.   | Band representations of BOTPs                                           |  18 |
|      | 1. Symmetry obstruction with open boundaries                            |  18 |
| C.   | Diagnosing BOTPs using the Wannier spectrum                             |  19 |
|      | 1. Wannier band representations                                         |  20 |
| IV.  | BOTP and filling anomalies in other 2D systems                          |  20 |
| A.   | Generalities on filling anomalies                                       |  20 |

|     | B. Filling anomalies in 2D systems                                |   22 |
|-----|-------------------------------------------------------------------|------|
| V.  | 3D models with boundary obstructions                              |   22 |
| A.  | Recipe for constructing 3D BOTPs                                  |   22 |
| B.  | 3D BOTPs with filling anomalies                                   |   24 |
|     | 1. A new class of C 2 nh Hamiltonians                             |   24 |
|     | 2. Real space picture and filling anomaly                         |   25 |
|     | 3. Wannier spectrum                                               |   26 |
|     | 4. A hinge-obstructed 3D BOTP                                     |   29 |
| C.  | 3D BOTPs with gapless surface states                              |   29 |
|     | 1. Dimerized weak Chern insulator                                 |   30 |
|     | 2. Dimerized weak quantum spin Hall insulator                     |   31 |
| VI. | Discussion                                                        |   33 |
|     | Acknowledgments                                                   |   34 |
| A.  | Character table of D π 2                                          |   34 |
| B.  | Real Space Proof of Anti-commuting Symmetries                     |   34 |
|     | 1. Anticommuting Mirrors in the 2D DMQI                           |   34 |
|     | 2. Anticommuting C 2 z and M z for 3D C 2 nh Models               |   35 |
| C.  | Band representations of BOTPs                                     |   35 |
|     | 1. Induction of a band representation from a local representation |   35 |
|     | 2. Movable Wyckoff positions                                      |   36 |

References

37

## I. INTRODUCTION

A symmetry protected topological (SPT) phase is a gapped phase of matter that cannot be adiabatically deformed to a trivial phase without breaking the symmetry or closing the bulk energy gap [1-8]. For free fermion systems with internal symmetries, a complete understanding of SPTs was achieved in the pioneering work of Refs. [4, 9-11]. In any given symmetry class and dimension ( d ), an SPT hosts anomalous gapless ( d -1) -dimensional surface states whose existence is tied to the non-trivial bulk topology. This connection between bulk and boundary topological properties is known as a bulkboundary correspondence.

For spatial or crystalline symmetries, the relationship between bulk and boundary topological signatures can be more subtle [12, 13]. On one hand, there are topological crystalline phases, such as mirror Chern insulators [14, 15], that exhibit a conventional bulk-boundary correspondence signaled by the appearance of 2D gapless surface states on any mirror invariant surface plane. On the other hand, less traditional types of surface states, known as higher-order surface states [13, 1625], which have a lower dimensionality, e.g., states confined to corners or hinges of the sample, are also possible in systems with crystalline symmetry. Even in the absence of any midgap surface (or higher-order boundary) states, additional signatures such as fractional charge at corners [26, 27] or defects [16, 28-30] can be used to distinguish different SPTs. Furthermore, there are some bulk SPTs protected by crystalline symmetries that are not associated with surface states or boundary fractional charges at all. For instance, two atomic insulators corresponding to filling the same Wyckoff positions, but with orbitals that transform differently under site symmetries, cannot be smoothly deformed to each other, but they cannot be distinguished by surface states or corner charges either [29, 31-33]. Thus, for SPTs protected by crystalline symmetries, bulk topological distinctions do not always imply boundary signatures in the form of gapless states or fractional charges at the boundary.

The purpose of this work is to investigate the reverse question. Rather than asking whether the distinction between two bulk SPTs can be captured by a boundary signature, we ask whether two systems which exhibit different boundary signatures, e.g., low-energy modes or fractional charges, necessarily correspond to topologically distinct systems in the bulk. There is an obvious counter-example to this statement corresponding to the case where non-trivial SPTs are placed on the boundary of a trivial bulk. For example, one can imagine gluing 2D layers on the surface of a 3D trivial insulator. If these layers are gapless in their two-dimensional bulk, then the 3D system will have surface states. If the layers are 2D topological insulators instead, then the system will exhibit gapless hinge modes. However, these boundary signatures are not associated with a bulk property. These phases have been dubbed

'extrinsic higher-order" phases and their surface states can be removed or 'peeled off' by means of symmetric boundary manipulations[13, 24]. However, this definition does not distinguish the cases where the boundary modes are completely decoupled from the bulk, e.g., a completely trivial bulk is attached to a non-trivial boundary, from those where the boundary modes are associated with the bulk in some way such that they can be created or destroyed by tuning bulk parameters.

One of the first examples of the latter phenomenon is provided by the quadrupole model of Ref. 17. This is a 2D model with gapped bulk and edges that hosts quantized fractional charge at the corners. Two symmetry variants of the model were considered: one with fourfold rotation symmetry C 4 z which satisfies ( C 4 z ) 4 = -1 , and the other with two anticommuting mirror symmetries. In the former case, the quadrupole phase is a topological phase in the standard sense: it is an obstructed atomic limit [32], which denotes an atomic insulator whose charge centers are displaced relative to the underlying positive atoms, and which is separated from the trivial atomic limit by a bulk gap-closing phase transition. However, in the latter case, the quadrupole phase is not a topological phase, in the standard sense, since the value of the quadrupole moment can be changed without closing the bulk gap (in periodic boundary conditions). The topology in this model is captured by a more subtle distinction contained in the spectrum of the Wilson loop operators (Wannier spectrum), or in the entanglement spectrum[20, 34], rather than the bulk energy spectrum. This was argued to imply that the two phases of the model are separated by an edge (rather than a bulk) phase transition when considered with open boundary conditions [17, 18]. Indeed, since both the gaps in the bulk and on the edge protect the quantized corner charge, one can imagine the corner charge delocalizing along the edges and changing values if the edge gap closes.

Despite the number of subsequent works that have studied several aspects of the quadrupole model and its generalizations, the subtle topological distinction captured by the double-mirror quadrupole insulator (DMQI) has ramifications that have been mostly overlooked. In particular, several questions regarding the nature of topological distinctions that are not captured by a bulk phase transition remain unanswered. These include: (i) how can one, in general, define a topological distinction that does not involve a gap closing phase transition in the bulk? (ii) are there other examples of models in two or three dimensions that exhibit similar phenomenology? (iii) under what conditions does an obstruction in connecting two Wannier spectra indicate a gap-closing transition at the boundary?, and (iv) how is this related to boundary signatures, e.g., fractional corner charge, or boundary states, in a general setting?

In this work, we answer these questions by introducing the concept of a boundary obstructed topological phase (BOTP) that captures distinctions between Hamiltonians that can be adiabatically connected with periodic boundary conditions, but cannot be connected for symmetric surface terminations with open boundary conditions. We provide a general definition of such distinctions, and show how they can be understood in terms of the Wannier spectra as well as real-space

symmetry representations. We show that the presence of topologically robust boundary modes in BOTPs depends on the lattice termination, and, inversely, that for any BOTP there is a choice of boundary characterized by topologically robust modes protected by a boundary gap closing.

In order to make the presentation as clear as possible, we will first start with detailed explanations of different aspects of the boundary obstruction concept that are realized in the DMQI. Then we will move on to discuss the general aspects of boundary obstructions. After reviewing the DMQI model in Sec. II A, we provide a real space framework for its boundary obstruction in Sec. II B, followed by a detailed analysis of the Wannier spectra, and how they relate to the physical edge spectra of particular surface terminations in Sec. II C. Afterwards, we establish how the boundary obstruction in the DMQI model is related to the existence of fractional corner charge via a detailed symmetry analysis of the model with open boundaries in Sec. II D. Finally, in Sec. II E we show how the notion of boundary obstructions in the DMQI can be understood in terms of restricting the bulk (periodic) symmetry representations to representations of the point group of the open system.

After the detailed study of the DMQI model, we introduce the general definition of boundary obstructed phases, and discuss their stability from a homotopic point of view in Sec. III A. In Sec. III B we present a definition of BOTPs from the point of view of the symmetry representations in the open system, and show how to diagnose them through the symmetry representations of the Wannier spectrum.

Afterwards, we discuss other possible 2D models with boundary obstructions in Sec. IV B and show that, apart from the DMQI (and variants of it), all such models require specific, complicated surface terminations to probe the boundary obstructions and the associated corner charges. In Sec. V we introduce a general recipe to generate 3D models with various types of boundary obstructions from a 2D building block. Such obstructions can be associated with fractional corner charges, fractional hinge charge density (Sec. V B), or boundary states (Sec. V C) depending on the choice of the 2D building block. For example, in the latter case we introduce two 3D models built from either a 2D Chern insulator or a 2D quantum spin Hall insulator with robust one dimensional states localized at their hinges. We also provide a complete characterization of all the 3D models using the Wannier spectra. Finally, we close with some concluding remarks and discussion in Sec. VI.

## II. 2D QUADRUPOLE INSULATOR AND BOUNDARY OBSTRUCTIONS

## A. Review of the quadrupole insulator model

We start our discussion by reviewing the DMQI model introduced in Ref. [17]. It consists of four orbitals describing spinless fermions arranged on a two-dimensional rectangular lattice with dimerized hopping amplitudes along the x and y directions, and π -fluxes threaded through each plaquette as

γ y λ y ( a ) · ( b, c ) · ν x ( π ) = 1 / 2 ν x ( π ) = 0 ν y ( π ) = 0 ν y (0) = 0 ν x ( π ) = 0 FIG. 1. Lattice model used to define the quadrupole insulator in Eq. 1. In (a), dashed lines have a relative negative sign to account for a flux of π threading each plaquette. The flux is the origin of the anticommuting reflection operations. (b) Band structure of the quadrupole insulator with γ /λ x x = 0 5 . and γ /λ y y = 0 4 . . Each energy band is twofold degenerate for a total of four bands.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00011v2-figure-72c47d89-a6ca-4f74-8a8e-a13cc742be2a.jpg)

2

2

/

•

(

d

)

/

(0) = 1

γ

x

) = 1

π

(

λ x ν x (0) = 0 ν x (0) = 1 / 2 ν y ν y ν x (0) = 0 shown in Fig. 1. The Bloch Hamiltonian is given by

$$\mathcal { H } ( k ) & = \left [ \gamma _ { x } + \lambda _ { x } \cos ( k _ { x } ) \right ] \Gamma _ { 4 } + \lambda _ { x } \sin ( k _ { x } ) \Gamma _ { 3 } \\ & + \left [ \gamma _ { y } + \lambda _ { y } \cos ( k _ { y } ) \right ] \Gamma _ { 2 } + \lambda _ { y } \sin ( k _ { y } ) \Gamma _ { 1 }, \quad ( 1 )$$

where γ x and γ y are hopping amplitudes within a unit cell along x and y respectively, and λ x and λ y are the inter-cell hopping amplitudes to nearest neighbor unit cells. The negative signs, represented by the dashed lines in Fig. 1a, correspond to our gauge choice for the π -flux threaded through each plaquette. The Γ 0 ,..., 4 matrices satisfy { Γ i , Γ j } = 2 δ i,j , and are represented explicitly by Γ 0 = σ τ 3 0 , Γ k = -σ τ 2 k , Γ 4 = σ τ , 1 0 for k = 1 2 3 , , , where σ 0 ,..., 3 ( τ 0 ,..., 3 ) denote the Pauli matrices, and the tensor product is implicit.

The model in Eq. 1 has reflection symmetries along both x and y , which for spinless electrons satisfy M 2 x = M 2 y = 1 . They are represented by

$$& \text{ar} \\ & \text{or} \\ & \text{ice} \\ & \text{$\Phi$} \quad \Phi$ \quad \Phi$ \quad \Phi$ \quad \Phi$ \quad \Phi$ \quad \Phi$ \quad \Phi$ \quad \Phi$ \quad \Phi$ \quad \Phi$ \quad \Phi$ \quad \Phi$ \quad \Phi$ \quad \Phi$ \quad \Phi$ \quad \Phi$ \quad \Phi$ \quad \Phi$ \quad \Phi$ \quad \Phi$ \quad \Phi$ \quad \Phi$ \quad \Phi$ \quad \Phi$ \quad \Phi$ \quad \Phi$ \quad \Phi$ \quad \Phi$ \quad \Phi$ \quad \Phi$ \quad \Phi$ \quad \Phi$ \quad \Phi$ \quad \Phi$ \ \Phi$ \quad \Phi$ \quad \Phi$ \quad \Phi$ \quad \Phi$ \quad \Phi$ \quad \Phi$ \quad \Phi$ \quad \Phi$ \quad \Phi$ \quad \Phi$ \quad \Phi$ \quad \Phi$ \quad \Phi$ \quad \Phi$ \quad \Phi$ \quad \Phi$ \quad \Phi$ \quad \Phi$ \quad \Phi$ \quad \Phi$ \quad \Phi$ \quad \Phi$ \quad \Phi$ \quad \Phi$ \quad \Phi$ \quad \Phi$ \quad \Phi$ \quad \Phi$ \quad \Phi$ \quad \Phi$ \quad \Phi$ \quad \Phi$ \quad \Psi$ \quad \Phi$ \quad \Phi$ \quad \Phi$ \quad \Phi$ \quad \Phi$ \quad \Phi$ \quad \Phi$ \quad \Phi$ \quad \Phi$ \quad \Phi$ \quad \Phi$ \quad \Phi$ \quad \Phi$ \quad \Phi$ \quad \Phi$ \quad \Phi$ \quad \Phi$ \quad \Phi$ \quad \Phi$ \quad \Phi$ \quad \Phi$ \quad \Phi$ \quad \Phi$ \quad \Phi$ \quad \Phi$ \quad \Phi$ \quad \Phi$ \quad \Phi$ \quad \Phi$ \quad \Phi$ \quad \Phi$ \quad \Phi$ \quad \Phi` \quad \Phi$ \quad \Phi$ \quad \Phi$ \quad \Phi$ \quad \Phi$ \quad \Phi$ \quad \Phi$ \quad \Phi$ \quad \Phi$ \quad \Phi$ \quad \Phi$ \quad \Phi$ \quad \Phi$ \quad \Phi$ \quad \Phi$ \quad \Phi$ \quad \Phi$ \quad \Phi$ \quad \Phi$ \quad \Phi$ \quad \Phi$ \quad \Phi$ \quad \Phi$ \quad \Phi$ \quad \Phi$ \quad \Phi$ \quad \Phi$ \quad \Phi$ \quad \Phi$ \quad \Phi$ \quad \Phi$ \quad \Phi$ \quad \Phi$ \$$

The model is also invariant under the combination M M x y = C 2 z , where C 2 z is a two-fold rotation symmetry. The π -flux leads to the anticommutation of the two reflection operators { M ,M x y } = 0 . As a result, C 2 z satisfies ( C 2 z ) 2 = -1 , as would be the case for spinful fermions. Other than the crystalline symmetries, the model H ( k ) , as written, also has (spinless) time-reversal symmetry given by complex conjugation T = K , as well as chiral symmetry {H , Γ 0 } = 0 . As a result, it lies in class BDI, though these internal symmetries are not crucial for the conclusions.

One important aspect to highlight here is that, due to the π fluxes, the spatial symmetries of the model as defined above already include some gauge transformations. As a result, the symmetry representation in question is projective. Although the implementations of M x and M y themselves are gauge dependent, the quantity

$$M _ { x } M _ { y } M _ { x } ^ { - 1 } M _ { y } ^ { - 1 } = - 1$$

captures the gauge invariant plaquette flux. We note that since the flux attachment does not affect the translation symmetries ( T x and T y still commute), and since the spatial symmetries are symmorphic, the full space group is given by a simple product of the point group F = D π 2 and the group of 2D translations i.e. G = F × T x × T y . The projective representation considered can be equivalently understood as a Z 2 extension of the point group D 2 defined by Eq. 3. We refer to such Z 2 extension with a π superscript, i.e., D π 2 . The point group D π 2 has four one-dimensional, irreducible representations ( A 1 , A 2 , B 1 , B 2 ), and one two-dimensional irreducible representation ( ¯ E ) (see Appendix A). We notice that the full projective representation defined by Eq. (3) is equivalent to a double (spinful) representation for the space group 25 ( Pmm 2 ) restricted to the 2D plane or equivalently the layer group 23 ( pmm 2 ). This aspect was discussed in the original works [17, 18], and investigated in more detail recently [35], but for the rest of this work we will focus on the spinless version to simplify the intuition for the degrees of freedom in the unit cell.

The energy spectrum of the Hamiltonian (1) is two-fold degenerate and gapped across the entire bulk Brillouin zone (BZ) (cf. Fig. 1) unless both | γ /λ x x | = 1 and | γ /λ y y | = 1 . Hence, we can connect the Hamiltonians at any pair of points in the ( γ /λ , γ x x y /λ y ) -plane without closing the bulk gap. This implies the absence of any bulk topological distinctions in the model. We note that the presence of C 4 z symmetry alters this conclusion since it forces λ x = λ y = λ and γ x = γ y = γ , thus changing the sign of | γ/λ | -1 is necessarily accompanied by a bulk-gap-closing phase transition, which indicates that the model has at least two distinct bulkprotected topological phases.

Despite the absence of bulk topological distinctions in the model protected by mirror symmetries, Refs. 17 and 18 have uncovered a more subtle topological distinction encoded in the topology of the Wannier bands instead. The Wannier bands along the x y ( ) direction are obtained as the eigenspectrum of a Wilson loop operator along this direction for a fixed momentum k y ( k x )[36, 37]. The Wilson loop operator is denoted W b ( k ) , where b is the direction along which the operator is taken in the Brillouin zone. W b ( k ) is a unitary operator whose eigenvalues have the form e 2 πiν b ( k ) . The values ν b ( k ) are defined modulo one, and represent the positions of the charge centers (within the unit cell) of hybrid Wannier functions that are maximally localized in the b -direction, but are delocalized Bloch waves in the perpendicular direction parameterized by the transverse momentum k [36, 37]. For the DMQI model, the two Wannier bands of the occupied states are generally gapped and symmetrically displaced away from the high-symmetry lines ν b = 0 and ν b = 1 2 / (cf. Fig. 2b). Among other things, this indicates that there is no spectral flow as would be the case, for example, for a strong topological insulator [38].

Whenever the two Wannier bands are separated by gaps from above and below we can consider the projector onto one of these two bands. This projector is effectively projecting onto the ground state of a 1D reflection-symmetric insulator with one electron per unit cell. Such an insulator can be char- acterized by a half-quantized polarization p, measured in integer units of the charge, distinguishing the cases where the charge center is at the center ( p = 0 ), or the edge ( p = 1 2 / ), of the 1D unit cell. Since there are two possible sets of Wannier spectra, i.e., one along the x -direction and one along the y -direction, there are two distinct quantized Wannier band polarizations p x,y = 0 1 2 , / and p y,x = 0 1 2 , / where p i,j denotes the polarization in the j -direction for a band taken from the Wilson loop in the i -direction. This yields a Z 2 × Z 2 invariant capturing the topology of the Wannier bands. Such invariants are protected by the gap in the Wannier spectrum rather than the energy spectrum, which means that their value can be changed without going through a bulk gap closing. As we will explain in detail later, these types of transitions can in some cases be associated with an energy gap-closing at the edge , rather than the bulk, when the system is considered with open boundary conditions in both directions. This will be the topic of Sec. II C.

For p ν = ( p x,y , p y,x ) = (1 2 1 2) / , / , the model exhibits mid-gap corner modes when considered on a rectangular geometry with open boundaries in both directions, (and when the boundary coincides an edge of the unit cell). However, these corner modes are protected by chiral symmetry which can be broken without changing the Wannier band polarization (whose quantization relies on only the mirror symmetries). Thus, these corner modes are not generically associated with the quantized Wannier band polarization. Instead we can associate the (1 / 2 1 , / 2) phase with corner charge [17, 18, 26, 29, 39] as will be discussed in detail in Sec. II D. It is worth noting that the fractional corner charge can be removed by adding edge degrees of freedom. Adding fractional corner charge to the trivial phase, will interchange which phase we identify as topologically non-trivial, but the fact that the two phases can be distinguished by the fractional corner charge remains true. This topological distinction will be discussed in Sec. II D.

Following this review of the DMQI, we will dedicate the remainder of this section to showing that the topology of the DMQI can be captured using the notion of edge topological obstructions; a new type of topological obstruction that is present only when the model is placed on certain geometries with open boundary conditions. In particular, we will provide an intuitive picture for a topological obstruction that is associated with an edge gap-closing transition rather than a bulk one. Furthermore, we will show how the Z 2 × Z 2 topological distinction encoded in the Wannier band polarizations is related to the corner charge, which is determined by only a single Z 2 invariant that distinguishes p ν = (1 2 1 2) / , / from the other three cases. We note here that picking p ν = (1 2 1 2) / , / as the non-trivial phase relies on the boundary termination convention shown in Fig. 1 which is chosen such that the vertical (horizontal) edge is invariant under M y ( M x ), and the edge lies at the border of the unit cell depicted in the figure (denoted by x, y = ± 1 / 2 ). As we will discuss later, the choice of the boundary decides which Wannier phase is distinguished from the others. In all cases, however, the resulting classification for a fixed boundary is Z 2 . This is similar to the SuSchrieffer-Heeger chain where changing the boundary termi-

high

)

R

LDOS(

low

z C 2 C 4 C 6 FIG. 2. (Color online) (a) Band structure of the two-dimensional quadrupole insulator with open boundaries in a single direction, y (left) or x (right). The blue bands correspond to the bulk band continuum, and the black bands that are slightly separated from the bulk bands are localized at the edges. (b) Wannier bands ν y ( k x ) (left) and ν x ( k y ) (right), with a quantized Wannier band polarization of p x,y = p y,x = 1 2 / . The Wannier spectra are topologically equivalent to the surface bands in (a) when we identify the gap around ν i ( k ) = 1 / 2 to represent the surface gap, see Sec. II C. (c) Local density of states at zero energy LDOS( ) = r -Imtr( H ( r ) + iη ) -1 /π for a small broadening η = 0 1 . in a system with open boundaries in x and y . The inset shows the existence of zero energy states inside the gap.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00011v2-figure-d5f0ab38-375d-493d-a5ac-f68f52df0366.jpg)

nation by half a unit cell will change the identification of the trivial and obstructed atomic limits.

## B. Real space picture

In the following, we will present an understanding of the edge topological obstruction of the DMQI in terms of a real space picture. At half-filling, the space of filled bands is Wannier representable, i.e., it is possible to find a basis of localized symmetric orbitals that span the subspace of filled bands. In a Wannier representable system, localized orbitals are labelled by their so-called Wyckoff position Q = { q i } which denotes a set (symmetry orbit) of spatial positions q i with i = 1 , ..., N Q , which is invariant as a whole under the symmetry group G . The positions q i can map to each other under some elements of G , q i = g i q for some reference q , but they are invariant under a (site symmetry) subgroup G q ⊆ G 1 . N Q denotes the size of the set of spatial positions q i , and is called the Wyckoff multiplicity of position Q . Wannier representable (atomic) insulators have electronic configurations constructed by decorating the various Wyckoff positions with orbitals that transform under the site-symmetry representations of G q . These different Wannier representable phases transform under space group representations that are induced from a local representation of G q [40, 41], and are each associated with a band representation (BR) [32, 40-43] that encodes both its transformations under the symmetry group, and its Zak-Berry phases in momentum space. A short derivation of the band repre-

1 The site symmetry groups for different elements q i in the same Wyckoff position are equivalent up to conjugation by an element in G , so we denote them by the same symbol G . q

FIG. 3. Wyckoff positions for the space group of the DMQI model, with site symmetry groups D π 2 , D 1 and C 1 , for those Wyckoff positions with multiplicity 1, 2, and 4, respectively. Figure adapted from Ref. [44].

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00011v2-figure-ed8502cf-5d2d-4787-90f9-f83bdca3117c.jpg)

sentation from the local representation of G q is presented in Sec. III B.

In the DMQI case, the full space group is a direct product of the point group F = D π 2 and the group of translations of a rectangular lattice in two dimensions, G = F × T . As previously mentioned, it can be thought of a Z 2 extension of the layer group 23 ( pmm 2 ) using Eq. 3. This group has the same set of Wyckoff positions as the original layer group: (i) four maximal Wyckoff positions with multiplicity 1 corresponding to positions lying at the intersection of the two mirrorinvariant lines, i.e., at positions (0 , 0) , (0 1 , / 2) , (1 / 2 0) , and (1 / 2 1 , / 2) (labelled 1 a , 1 b , 1 c and 1 d , respectively) with G q coinciding with the full point group F ; (ii) four Wyckoff positions with multiplicity 2 corresponding to positions lying along one of the mirror-invariant lines, ( ± x, 0) , ( ± x, 1 2) / , (0 , ± y ) , (1 / 2 , ± y ) , (labelled 2 e , 2 f , 2 g , and 2 h , respectively) with G q = D 1 ; and (iii) one general Wyckoff position 4 i with multiplicity 4 at four generic symmetry-related points ( ± x, ± y ) where G q = C 1 . The positions in the unit cell associated with these labels are shown in Fig. 3.

To illustrate the absence of a bulk obstruction in the DMQI model, consider a translationally invariant lattice with M x and M y symmetries having periodic boundary conditions 2 . For a filling of two electrons per unit cell, there are only three allowed options for placing the atomic orbitals: (i) both are placed at the same maximal Wyckoff position, (ii) they are placed at a Wyckoff position of multiplicity 2, or (iii) the two orbitals are placed at different maximal Wyckoff positions. Case (iii) has a non-vanishing polarization of 1 / 2 along x or y (or both) which rules it out for the DMQI model for which both polarizations vanish. Interestingly, cases (i) and (iii) may be adiabatically deformed into each other in some cases, i.e., they may belong to equivalent band representations. This is the case when the two electrons (per cell) transform as a two dimensional representation of the D π 2 point group, i.e., when the matrix representations of the two mirror operators anticommute when acting on the electron orbitals. To show this,

2 In general, point group operations act non-locally in real-space, and generically relate orbitals in different unit cells. However, with periodic boundary conditions one can use translation symmetry on a transformed orbital to shift it back to the initial unit cell. This allows us to represent the different point-group operations as an action on Wannier centers inside a unit cell.

we note that equivalence between two band representations can be demonstrated if there is a symmetric, adiabatic deformation between the two configurations, see Refs. [32, 45] and App. C. In Appendix B we show that in order to carry out such a continuous deformation between Wyckoff position configurations, the two mirror operators are required to anticommute, which is exactly the case for the DMQI model. The necessity for anticommutation of mirror operators was previously proven in Ref. [17, 18] using the symmetry indicators of the Wannier bands, where it was shown to be linked to the existence of gaps in the Wannier spectra.

glyph[negationslash]

glyph[negationslash]

In the following, we will present an alternative and more intuitive argument based for the equivalence of these representations based on an illustration of real space deformations. This will be backed by a more rigorous analysis in subsequent sections. In Fig. 4 we explicitly show how to deform the configuration with both electrons at 1 a to the one where both are at 1 d while preserving the symmetry and the bulk gap. There are two distinct ways to do this. The first one goes through 1 c by first moving horizontally through 2 e then vertically through 2 h, whereas the second goes through 1 b by first moving vertically through 2 g then horizontally through 2 f . Such deformations can be written in terms of deformations of the parameters of the Hamiltonian. For example, the deformation 1 a → 1 b → 1 d starts with the parameters γ x = γ y = 0 and λ x = λ y = 0 which yields Wannier orbitals centered at 1 a . Next, we gradually increase λ x and decrease γ x to get to λ x = 0 and γ x = 0 whose Wannier orbitals are centered at 1 b . Finally, we increase λ y and decrease γ y to get to γ x = γ y = 0 and λ x = λ y = 0 which yields Wannier orbitals centered at 1 a . The center of the Wannier orbitals in these cases can be obtained easily since all of these are fully dimerized limits.

glyph[negationslash]

glyph[negationslash]

It is instructive again to compare to the model with C 4 z symmetry for which the movable Wyckoff positions 2 e , 2 f , 2 g , 2 h are not allowed. In this case, it is impossible to move two electrons from 1 a to 1 d without breaking the symmetry, indicating that the two correspond to distinct bulk phases. In terms of the Hamiltonian parameters, C 4 symmetry forces λ x = λ y and γ x = γ y . Under this constraint, it is impossible to go from λ x = λ y = 0 , γ x = γ y = 0 (corresponding to 1 a ) to λ x = λ y = 0 , γ x = γ y = 0 (corresponding to 1 d ) without going through the point γ x = γ y = λ x = λ y where the bulk Hamiltonian is gapless.

glyph[negationslash]

Let us now see what happens in the presence of a boundary termination which we take, for definiteness, to be a rectangle whose horizontal (vertical) edges are parallel to the x ( y ) axis and lie at a position x, y = 1 2 / . This means that the edges of the sample coincide with the the edge rather than the center of our chosen unit cell. We then imagine how the bulk deformations explained in the previous paragraph, namely the two possible deformations connecting 1 a to 1 d either through 1 b or 1 c , are affected by the presence of the vertical and horizontal edges. Such deformations are carried out by changing the parameters λ x,y and γ x,y of the Hamiltonian while keeping the boundary fixed. For definiteness, we consider a trajectory of bulk Hamiltonians H ( ) t such that H (0) = H 1 a and H (1) = H 1 d which, upon choosing a boundary termination, induces a trajectory of Hamiltonians in the open system

FIG. 4. Deformation of the Wannier centers (open circles), in the double-mirror quadrupole model, from the 1 a to the 1 d Wyckoff position along movable Wyckoff positions. We see that from panel (b) to panel (c) we cannot preserve the symmetry at the edges.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00011v2-figure-3ce27d13-10b7-4681-998c-4018ee5362bc.jpg)

H obc ( ) t (we provide a precise description of this in Section III A ). Deep inside the bulk, H ( ) t and H obc ( ) t are the same so that the picture of sliding the Wannier centers explained above makes sense inside the bulk of the open system. In addition, if the Hamiltonian H obc ( ) t is gapped for all t , we may be tempted to assume that the position of the charge centers remain well-defined also close to the edge as we change t . However, the schematic illustration of Fig. 4 suggests that this is impossible, as we now discuss in detail.

To explain the intuitive picture behind the boundary obstruction depicted in Fig. 4, let us consider the trajectory connecting 1 a to 1 d through 1 c . We note that in the presence of the vertical boundary at x = 1 2 / , the electron filling of position 1 c on the edge is half of its filling in the bulk (since such a site is shared by two unit cells when in the bulk, but one of them is now absent on an edge). This means that there is only one electron at the 1 c position at the edge, instead of two. This makes it impossible to vertically move it to position 1 d while preserving the symmetry and energy gap. Similarly, the presence of the horizontal edge at y = 1 2 / prohibits the deformation 1 a → 1 b → 1 d as position 2 f is unavailable at the edge. Thus, with open boundaries we can pass from 1 a to 1 b or 1 c , but not to 1 d, which distinguishes 1 d from the other three configurations in the presence of open, symmetric boundaries.

To summarize, Fig. 4 suggests there are two ways to move two electrons from 1 a and 1 d : one of which goes through position 1 c and one through 1 b . Importantly, each of the two edges ( x or y ) prohibits only one of the two trajectories. Thus, if we open the boundaries along only one direction by considering the system on a cylinder, it is still possible to deform 1 a to 1 d , and thus they are not distinguished. However, once we take open boundaries along both the x and y directions (with an edge termination consistent with the unit cell), these two atomic configurations cannot be continuously connected while preserving the symmetry. This picture suggests the DMQI exhibits what we will call an edge obstruction wherein two bulk Hamiltonians which are smoothly deformable in a periodic system (i.e. in the bulk) are not continuously deformable in the open system, i.e., any trajectory connecting them will involve a gap closing on the edge or breaking the symmetry. In the following two subsections, we will expand

µ

y

∆

µ

ν

y

(

k

on the intuitive understanding presented here by presenting two rigorous formulations of the boundary obstruction in the system: one in terms of edge phase transitions, and the other in terms of a Wannier sliding argument in the open system. x ) y

## C. Wannier bands and edge spectrum µ x x

We now switch our attention to the diagnosis of the edge obstruction in terms of the Wannier spectrum. As discussed in Sec. II A, and Refs. 17 and 18, the Wannier bands of the DMQI model can be characterized by a Z 2 × Z 2 invariant p ν = 1 2 (Θ(1 -| γ /λ x x | ) , Θ(1 -| γ /λ y y | )) (where Θ( ) x is the step function). Using the correspondence between the Wannier spectrum and the edge spectrum [46], we may naively think that the DMQI has four distinct phases that cannot be smoothly connected to each other without closing the energy gap at the edge. This, however, contradicts the analysis of the previous sections as well as the result of Refs. [17, 18] which found two, rather than four, distinct phases distinguished by the corner charge (or alternatively the quadrupole moment).

The resolution to this puzzle lies in the observation that the edge spectrum and the Wannier spectrum differ in one fundamental aspect, despite being continuously deformable to each other. Namely, the Wannier spectrum is periodic, whereas the edge spectrum is not. This means that of the two gaps of the Wannier spectrum, one at ν = 0 and one at ν = 1 2 / , only one corresponds to the actual energy gap at the edge. Thus, only Wannier transitions which involve closing this particular gap correspond to actual edge obstructions even though gap closings at either ν can act to change Wannier topology; this is an important subtlety that we treat in detail below. The determination of which Wannier gap corresponds to the edge energy gap depends on details, including the surface termination as we will see below.

To understand the relationship between a surface spectrum and a Wannier spectrum, let us follow Ref. 46 in implementing the edge via the replacement

$$\mathcal { H } \to \mathcal { H } _ { \text{obc} } ( x ) = P \phi _ { M } ( x ) P + M ( 1 - P ) \quad ( 4 ) \quad \text{$\quad $\mathbb{s.h}$} \\ \text{$\quad$tered}$$

where P is the projector onto the filled bands of the periodic Hamiltonian H pbc and φ M ( x ) is a regularized linear potential given by

$$\phi _ { M } ( x ) = \begin{cases} x, & | x | < M \\ \text{sgn} ( x ) M, & | x | > M. \end{cases} \quad ( 5 ) \quad \text{res} \quad \text{a ga}$$

For a given filling specified by a chemical potential µ satisfying -M &lt; µ &lt; M , this Hamiltonian reduces to the spectrally-flattened bulk Hamiltonian with filled (empty) states at energy -M M ( ) in the bulk. Thus, H obc ( x ) has the same filled states as H pbc deep within the sample, and it implements the vacuum Hamiltonian, where all eigenstates are empty (trivial projection operator), far outside the sample. Although the existence of an edge obstruction does not depend on the details of the chosen potential φ x ( ) , the choice (5) makes the connection to Wannier bands manifest since the

1.5

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00011v2-figure-42a5a709-3712-48cb-a1a3-808d839143d6.jpg)

-

π

0

π

200 150 135 FIG. 5. Spectrum for the Hamiltonian H obc ( x ) implementing the edge termination with the linear potential φ M ( x ) (Eq. 5) truncated at M = 2 5 . . We can see that the spectrum for energies localized close to the edge matches exactly the Wannier spectrum. Bulk states are all accumulated at energies ± M .

linear potential sets a direct proportionality between distance and energy 3 .

The resulting spectrum of H obc is shown in Fig. 5. We can see that the spectrum for the states localized close to the edge consists of several copies of the Wannier spectrum shifted relative to each other by some integer. The rest of the spectrum far away from the edge accumulates close to the energies ± M and represent states in the bulk ( -M ) or outside of the sample ( + M ). It is clear from Fig. 5 that, in this potential, the edge termination is decided by choosing an energy that separates the filled states inside the sample from the empty states outside, i.e., an 'edge chemical potential." The linear form of the potential φ M ( x ) gives a direct relation between the edge chemical potential and the 'position" of the edge termination. Since the edge energy spectrum is just repeated copies of the Wannier spectrum, the edge chemical potential determines a corresponding 'Wannier chemical potential" (WCP). The key implication of this analysis is that only Wannier gap-closings that happen at a position corresponding to the WCP indicate a gap closing in the edge energy spectrum. This establishes the important link between gap closings in the Wannier spectrum and at the edge. 4

Similar to the previous subsection, we take a unit cell centered at 1 a (modulo integers) with edges lying on the lines x = 1 2 / and y = 1 2 / (positions 2 h and 2 f ) that intersect at the position 1 d where the corner is located. This termination corresponds to the WCP ν x = 1 2 / and ν y = 1 2 / . Consequently, a Wannier gap-closing at ν x = 1 2 / or ν y = 1 2 / represents a genuine edge gap closing in our convention, whereas a gap-closing at ν x = 0 or ν y = 0 , does not. Instead the latter type of gap closing indicates mixing between different states at the edge, or states at the edge mixing with the bulk, without

3 To match the units between energy and distance, we note that the Hamiltonian H obc ( x ) is a dimensionless band-flattened Hamiltonian where filled (empty) bands correspond to -M ( + M ). Similarly, the position x is measured in units of the lattice constant and can be taken to be dimensionless.

4 Throughout this work, a gap closing at the edge means that the states for which the gap closes are localized close to the edge with localization length much smaller than the system size.

FIG. 6. Diagram of the phases of the Wannier bands of the DMQI model, Eq. 1. On the yellow and blue lines a Wannier gap-closing transition occurs in the Wannier spectra W x and W y , respectively. The gap closing occurs either at ν i = 0 or ν i = 1 2 / as indicated. At the intersection of the two lines, marked with a red dot, the bulk spectrum is gapless. To each region, we can associate the Wannier band (nested) polarization p ν Q , with Q the associated maximal Wyckoff position. Black dots indicate specific configurations corresponding to the plots (a-e) in , where we show the evolution of W x and W y along the path 1, indicated with a dashed line. The phase (1 / , 2 1 / 2) 1 d is always separated from the remaining ones by a Wannier gap closing at the Wannier chemical potential µ = (1 2 1 / , / 2) , and hence is separated by a gap closing at the physical boundary for our choice of boundary termination. Choosing a different path from (a) to (e), for example dashed line 2, would still imply a boundary gap closing although at a different edge.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00011v2-figure-0502c8af-1ea8-444f-9feb-ef03ac828454.jpg)

closing the edge gap.

Using the preceding analysis, we can now show that for our choice of lattice termination, the state having p ν = (1 / 2 1 , / 2) is distinct from the other three Wannier polarization configurations. The latter three can all be continuously connected while preserving symmetry, but there is an edge obstruction when trying to go from one of (0 , 0) (1 , / 2 0) (0 1 , , , / 2) to (1 / , 2 1 / 2) . We can use Fig. 6 to illustrate this conclusion. This figure provides a diagram of all possible Wannier band phases as a function of the ratios γ /λ x x and γ /λ y y . The blue central square in the diagram corresponds to the Wannier polarization p ν = (1 / 2 1 , / 2) ; whereas other regions correspond to Wannier polarization equal to either p ν = (1 2 0) / , (orange regions), p ν = (0 1 2) , / (red regions), or p ν = (0 0) , (white regions). As we can see from the figure, going from the (1 / 2 1 , / 2) phase to any other phase involves at least one Wannier transition at ν x = 1 2 / or ν y = 1 2 / , whereas going between any of the three other phases can be achieved without a Wannier transition at ν i = 1 2 / . Wecan conclude that, in our convention, the (1 / 2 1 , / 2) phase is separated from the other

2 0)

ν y 1 20 40 0.5 1 ν x 1 20 40 0.5 1 (f) p ν Q = (0 0) , 1 a ν y = 0 It is also worth pointing out that the the equivalence of edge and Wannier spectra concerning the gaps only applies for a linear potential. In general, the two are continuously rather than smoothly connected such that a gap in one does not imply a gap in the other. This was pointed out in a recent work [47] which studied a modification of the DMQI with further neighbor hopping. This means that, in general, different choices of the function φ M ( x ) in (4) can move the 'phase boundaries' between different obstructed phases. On the other hand, the overall number of phases does not depend on such a choice since they can be diagnosed by discrete features such as filling anomalies or symmetry representations. We will discuss the detailed dependence on the boundary termination in Sec. III A.

p ν Q = (0 0) , 1 a p ν Q = (1 / , (a) (b) three phases by an edge gap-closing, whereas the three other phases are continuously connected, i.e., any pair of them can be smoothly connected without closing the edge gap or breaking symmetry. This is the resolution to the discrepancy between the topological characterization in terms of edge spectrum vs. Wannier spectrum. We note that if we changed our edge termination, e.g., choosing to terminate the unit cell at x = 0 , y = 1 2 / , the only effect would be to permute which of the four polarization configurations should be distinguished from the other three.

eigenstate 0 eigenstate 0 eigenstate eigenstate Having set up this background, we can now study the connection between Wannier polarization and the real space picture of the edge obstruction. Fig. 7 shows a detailed version of the deformation process in Fig. 4, and illustrates the different Wannier polarization configurations associated with each intermediate state. In Fig. 7(a), we have a phase characterized by p ν = (0 0) , . As the Wannier centers become separated horizontally away from 1 a , a first Wannier transition occurs between Fig. 7(a) and (b) where the y -Wannier gap closes at ν y ( π ) = 0 . If we interpret the edges of the geometry to be open instead of periodic then this transition, which changes p ν to (1 / 2 0) , , leaves the outermost Wannier centers (those on the left an right edges) unpaired [red open circles in Fig. 7(b,c)] and pinned to the position ν y = 0 . In contrast, the bulk Wannier centers are still free to move, and thus can adopt Wannier eigenvalues that come in pairs ν, -ν along both x or y.

We can verify this piece of the deformation through the explicit calculation of the Wannier spectra ν x ( ν y ) of the DMQI on a cylindrical geometry in which boundaries are closed along x ( y ) and open along y ( x ). In Fig. 7(f) we see that all of the Wannier eigenvalues come in ± ν i pairs, which matches our picture of having bulk Wannier functions with none isolated on the boundaries. In comparison we can look at Fig. 7(g), in which there is a pair of isolated Wannier eigenvalues (for every unit cell along the y -direction) exactly pinned at ν y = 0 . If we look at the eigenstates associated to these pinned Wannier eigenvalues, we find that one set is localized at one of the open, vertically oriented edges, and the other set is on the other vertical edge. Notice, however, that since their eigenvalues are pinned to ν y = 0 , these edge Wannier states lie at the center of the unit cell in the direction parallel to the edge, and thus amount to a trivial edge polarization. Also, notice that such edge isolated Wannier states do not exist for the

1

c

(

ν

x

x

x

FIG. 7. Correspondence between real-space configurations and Wannier bands in the DMQI model, Eq. 1. (a-e) Wannier center location and nested (Wannier band) polarizations along a path illustrated schematically in Fig.6 (path 1). A Wannier gap closing transition happens between panels (a) to (b) at ν y = 0 , and between (d) and (e) at ν x = 1 2 / . The latter happens at the Wannier chemical potential and corresponds to a boundary gap closing, and consequently to a topological phase transition. (f-h) Wannier spectra for cylindrical geometries for the three Wannier phases in (a-e). Wannier eigenvalues in red have eigenstates that are localized on the edges of the cylinder. These results exactly match the configurations of the real space pictures in (a-e). The panels are labelled by the Wannier band polarization p ν Q where the subscript is an indication of the maximal Wyckoff position Q that is associated to p ν . The arrows in (a-e) are drawn with their center at Q as visual guide.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00011v2-figure-f6c926e7-9f22-40fd-899e-5e930ed39a30.jpg)

ν x spectrum.

Now we can proceed with the deformation process along y. A second Wannier transition occurs between Fig. 7(d) and (e). At this transition, there is a gap closing in the x -Wannier bands at ν x ( π ) = 1 / 2 , so that in Fig. 7(e) p ν = (1 2 1 / , / 2) . Notice the appearance of unpaired Wannier centers at the horizontal edges with Wannier eigenvalues fixed to ν x = 1 2 / [red solid circles in Fig. 7(e) along horizontally oriented edges]. Since the vertical edge has reflection symmetry, the values ν y for states on this edge are quantized to 0 or 1/2. Being initially pinned at ν y = 0 after the first Wannier transition, they undergo an edge phase transition that changes their value to ν y = 1 2 / [red solid circles in Fig. 7(e) along vertical edges]. This is further confirmed by the explicit calculation of the Wannier spectrum on a cylinder shown in Fig. 7(h), which shows pairs of isolated Wannier eigenvalues (per unit cell along the periodic direction) at both ν x = 1 2 / and ν y = 1 2 / . If we look at the eigenstates corresponding to these isolated Wannier eigenvalues we find that they are localized at the horizontal or vertical edges respectively. We emphasize that this last transition is the only one during this process that is accompanied by an edge phase transition which closes the edge energy gap, and that this transition is unavoidable due to the fact that the isolated states (red) on the vertical edge constitute effective one-dimensional obstructed atomic limits protected by M y symmetry.

A final remark is due regarding Fig. 7(e). The overall configuration of Wannier centers is such that reflection symmetries are preserved in the bulk, but not at the edges. Indeed, it is not even possible to maintain a symmetric configuration of Wannier centers for this system in an open boundary given the number of occupied states we need to fill the lower energy bands in periodic boundary conditions. This is an example of a filling anomaly , i.e., a condition where the constraint of charge neutrality (integer band filling) is not compatible with the spatial symmetry [26]. We will provide a more detailed discussion of this concept in the next subsection.

## D. Corner charge and filling anomaly

As discussed briefly in Sec. II A, the p ν = (1 / , 2 1 / 2) phase exhibits zero energy corner modes, but only when particlehole or chiral symmetries are present. In the absence of these symmetries, these corner states are not pinned to zero energy and, as a result, can be moved up and down in energy, though they remain degenerate because of the spatial symmetry [33, 48]. Once they are pushed to the conduction or valence band they can hybridize with the delocalized states there, and do not necessarily remain exponentially localized at the corners. Thus, the edge obstructed phase is not associated with corner modes in general.

The absence of zero-energy corner modes, however, does not mean that other signatures such as corner charges (rather than states) are absent. To investigate this possibility, let us start with the particle-hole/chiral symmetric case and then add terms that break these symmetries. We note that the DMQI has 4 N states in total, where N is the number of unit cells. In the particle-hole/chiral symmetric case, we know there are 4 zero-energy corner states leaving 2 N -2 states in both the

conduction and valence bands. Once particle-hole/chiral symmetry is broken, the corner states can move in energy but they remain degenerate due to mirror symmetries, leading to two possibilities: they either remain inside the gap or move into the valence or conduction band. In the first case, they remain localized eigenstates of the Hamiltonian having degeneracy protected by the two mirror symmetries. As a result, the system can not be gapped 5 at half-filling since there are two electrons that need to occupy four degenerate states. 6 In the second case, these states hybridize with the states in the valence (conduction) band forming new and more extended eigenstates. As a result, the number of states in the valence (conduction) band is changed to 2 N + 2 , and the chemical potential will lie in the valence (conduction) band at half filling which implies that the system is again not gapped. This phenomenon, in which it is impossible for a system to be simultaneously gapped, symmetric, and charge-neutral at a certain filling has recently been referred to as 'filling anomaly" [26, 33, 48-51]. Crucially, if we impose the conditions that the system is a symmetric gapped insulator, then the filling anomaly will manifest in excess or deficient charges distributed symmetrically in the four corners. For the DMQI, such charge is equal to ( n + e/ 2) per corner for some integer n 7 .

Before considering the general case in the next section, let us build some intuition about filling anomalies in bulk and boundary obstructed phases by contrasting two important examples: the bulk obstructed SSH chain in one dimension protected by inversion symmetry; and the boundary obstructed DMQI. First, imagine we take an SSH chain having N unit cells (with two orbitals per unit cell) at half-filling so that the total number of electrons is N . Following [50, 51], we will now show that, in some cases, it is impossible for this system to be gapped, inversion-symmetric, and charge-neutral when considered with open boundaries. Since every 1D system is Wannierizable, we can assign real space position corresponding to the centers of the atomic orbitals in the gapped, periodic system, i.e., Wannier centers. We now open the boundaries and denote the inversion center of the open system by O . For a sufficiently large system, the filling at a small inversion-symmetric region surrounding O , which we denote by ν O , should not be affected by opening the boundary and thus should be the same in the open and periodic systems. The assumption of charge neutrality means that the number of electrons in the open system is N , the same as the periodic system. We now note that it is impossible to find an inversionsymmetric arrangement of real space positions away from the inversion center O if N -ν O is odd. Thus, one of the three

5 Agapped system is one with no charge excitations at zero or infinitesimally small energy.

6 We note that the usual paradigm of slightly breaking the symmetry so that two out of four modes could be filled is unavailable as we are strictly enforcing the symmetry. Indeed filling the modes this way essentially results in spontaneous symmetry breaking in the thermodynamic limit.

7 Corner charges associated with filling anomalies are only defined modulo an integer charge since we are free to add integer charges to the corners while keeping the system a symmetric insulator.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00011v2-figure-58198a04-a8ac-42b1-b298-6b5feae2952a.jpg)

/negationslash

(a) (b) (c) ∼ ∼ ¯ E ⊗ ¯ E ¯ E FIG. 8. Bulk obstruction and filling anomaly in an SSH chain protected by inversion symmetry. With open boundaries there is a single inversion center (red dot), while when the boundaries are identified an antipodal inversion center is created at the stitching point (blue dot). (a, left) Trivial state: the filling of the inversion center (red dot) is odd just as the total filling of the chain is odd, which implies the electrons can be filled while preserving inversion. (a, right) When the end points are identified only the red inversion center is filled. (b, left) Nontrivial state: at an odd filling, if the red inversion center is not filled there is no way to fill the chain while preserving inversion, and there is a filling anomaly associated with end charges that need to be fractional at neutrality to preserve the symmetry. (b, right) Identifying the end points cures the anomaly, i.e., we no longer require fractionalization at neutrality, and only the blue inversion center is filled. With open boundaries the two phases can be distinguished by the endpoint charges, while with closed boundaries a topological distinction remains since the charge at the blue site cannot be moved to the red site while preserving both the symmetry and the bulk gap. This is an example of a bulk obstruction.

assumptions - gap, inversion symmetry, charge neutrality must fail in the open system. For example, if N -ν O is odd and we require the gap and neutrality, then there will be a charge imbalance on the left and right halves of the system such that inversion symmetry is broken. If we instead balanced the charge to preserve inversion symmetry then the system would no longer be neutral.

Another way to intuitively understand the filling anomaly is to see that at a filling corresponding to neutrality an electron would be forced to fractionalize into two pieces to preserve the inversion symmetry, as illustrated in Fig. 8. If we require that the system be neutral then the filling anomaly in the SSH chain generates the bulk polarization, i.e., it implies that the inversion symmetry will be broken spontaneously and a dipole moment will appear that creates opposite fractional charges of ± e/ 2 bound at the two edges. If we require that the system remains symmetric then it will either be gapless since there will be degenerate low-energy modes, or we will need to add an extra electron (at the minimum), which will violate neutrality.

As an aside it is instructive to relate filling anomalies to the recently developed framework of layer construction of topological crystalline phases [52, 53] since we can use it for the description of the filling anomaly of the DMQI model. Within the layer construction, any D -dimensional topological crystalline phase can be built as follows. We start with some d -dimensional region Γ , with d &lt; D , that is invariant under the spatial symmetry group. We then decorate this region with some topological phases whose dimension is smaller than or equal to d . Finally, we build the D -dimensional system by symmetrically attaching or 'adjoining" lower dimensional units to Γ until the whole D -dimensional space is filled. The topological invariants of the D -dimensional topological phase will be given by those of the lower dimensional topological

phase modulo those related by the adjoining procedure. The layer construction was recently employed for the classification of higher-order topological phases in Refs. [25, 54]. It was also employed to study obstructed atomic insulators and fragile phases which can be built by layering 0D units, i.e., orbitals, in Ref. [29]. For example, in the SSH chain example considered above, Γ is the origin O . Its topological invariant is the filling ν O which is an integer invariant protected by U(1) charge conservation. The adjoining operation is the symmetric addition of zero-dimensional inversionrelated units (which are just a pair of atomic orbitals) which changes the filling ν O by an even number. Thus, after modding out by this operation, we end up with a Z 2 invariant diagnosing the filling anomaly in the system.

A similar analysis can be performed for the filling anomaly of the DMQI model. We start by considering N N x y unit cells at a filling of two electrons per unit cell. With open boundaries, the symmetry group of the model is the point group D π 2 relative to O , with irreducible representations given in Appendix A, Table IV. Similar to the discussion on the SSH model, we will now show that sometimes it is impossible for the system to be a symmetric charge-neutral atomic insulator 8 .

Our argument proceeds very similarly to the 1D case. We first select a symmetry invariant region Γ consisting of the union of the x = 0 and y = 0 lines. If the filling of this region, ν Γ , is the same in the open and periodic system (this assumption will be justified below), then there are 2 N N x y -ν Γ remaining electrons that should be symmetrically arranged throughout the rest of the open system (excluding the region Γ ) to obtain a symmetric charge neutral atomic insulator. However, this is only possible if 2 N N x y -ν Γ is a multiple of 4 since every point in the open system that does not lie inside the region Γ is mapped by M x , M y , and M M x y to three other distinct points, hence leading to a symmetry orbit of size 4. This can be used to establish the existence of a filling anomaly. Notice that ν Γ is always even. This can be seen by noting that it consists of (i) the filling at the origin, (ii) the filling at the line x = 0 excluding the origin, and (iii) the filling at the line y = 0 excluding the origin. The fillings (ii) and (iii) are manifestly even due to mirror symmetries, whereas (i) is even since the bulk electrons at maximal Wyckoff positions always transform as a 2D ¯ E irrep in the DMQI model, regardless of the phase. Thus, ν Γ is always even and the filling anomaly is determined by the Z 2 invariant N N x y -ν Γ / 2 mod 2 .

The previous argument has one missing ingredient: how can we show that ν Γ is the same in the periodic and open systems? In the 1D case, this was simple to show since opening the boundary did not affect the vicinity of the inversion center O (which corresponds to the region Γ in that case). For the 2D case, this can be established only if we make the extra assumption that the polarization per unit cell in both directions vanishes: p x = p y = 0 , which holds for the DMQI

8 Note that here we need the extra assumption that the system is an atomic insulator, i.e., it admits a Wannier representation, which was not needed in the 1D case since all 1D systems are Wannierizable.

model. This assumption means that if we consider a cylinder geometry where translation symmetry is retained along one of the directions, there is no filling anomaly and we can have a symmetric insulator whose filling per unit length along the periodic directions is the same in the open and periodic systems. Thus, for a cylinder with open boundaries in the x -direction and PBC in the y -direction, the filling of the x = 0 line is the same as in the periodic system. Obviously, this filling is unaffected by opening the boundary in the y -direction which takes place far away from the x = 0 line. Thus, by opening the boundary in two steps, first in x then in y , we can show that the filling at the x = 0 line is the same in the open and periodic systems. Similarly, we can open the boundary first in y then in x and show that the filling at y = 0 is the same in the open and periodic systems. This establishes that the filling at Γ , the union of x = 0 and y = 0 , is the same in the open and periodic systems.

We remark that for the point group D π 2 , the Z 2 filling anomaly explained above is the only filling anomaly possible. It is instructive to compare this to the point group D 2 of spinless electrons without the π flux. In this case, the filling at the origin ν O is not necessarily even, and can lead to a different filling anomaly diagnosed by the parity of the full filling at the origin ν O . This type of filling anomaly reflects a bulk obstruction, just like it does for the 1D SSH chain, and diagnoses the sum of the polarizations along the x and y directions, each of which must be quantized by the mirror symmetries. To identify the individual x and y polarizations, we would need to consider cylindrical geometries with an open boundary along one of the directions and a periodic boundary along the other.

Though we have now identified the filling anomaly for the DMQI mode, the previous discussion has not yet clarified the distinction between filling anomalies in boundary-obstructed models like the DMQI, compared to bulk-obstructed models such as the SSH chain. Our goal now is to elucidate this difference. The distinction can be understood by studying what happens when the anomaly is resolved, i.e., removed, by identifying some of the edges of the open system. What we mean by 'resolved" is that the change of geometry to a periodic system now allows for a symmetric, gapped, and neutral configuration at our filling of interest. For example, in the SSH case, we can resolve the anomaly by identifying the two endpoints, which generates a new inversion center O ′ . In this case, it is possible to have a symmetric, charge-neutral insulator whose filling N has different parity from ν O since there is now a new inversion center that can compensate for the parity mismatch. In the phase with (without) the filling anomaly, the new inversion center O ′ will have an odd (even) filling whose parity cannot be changed by symmetrically adding or removing electrons. Thus, the filling anomaly in the open system corresponds to a bulk topological distinction in the periodic system. That is, while the anomaly can be resolved by identifying the edges, the topological distinction between the two phases still remains with periodic boundaries; they are distinct atomic limits.

We can similarly resolve the filling anomaly of the DMQI by identifying either pair of opposite edges. When the boundary is stitched together, the symmetry of the corner and the

FIG. 9. Sewing the boundaries in a boundary obstructed atomic insulator (e.g., the DMQI quadrupole insulator). When the boundary is stitched along one direction, the filling anomaly is lifted since the two corners contribute with half a charge so that all sites will now have integer filling. Interestingly, in this case the bulk obstruction is also lifted: the seam line is a high-symmetry line where the two states can freely move towards each other and symmetrically reach all maximal Wyckoff positions, including the original symmetry origin.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00011v2-figure-9075d680-d84c-45e0-8e79-4e390ade98d5.jpg)

seam line is enhanced compared to the open geometry. The filling anomaly is resolved following this boundary identification since only two electrons are required to fill the edge states after the stitching, as illustrated in Fig. 9. Another way to say this is that there is a new high symmetry line (invariant under M x or M y depending on the identified edges) whose filling can compensate for the mismatch between ν Γ and 2 N N x y modulo 4. However, in contrast to the SSH chain, the resulting system is completely trivial since the two electrons at the new high symmetry line can now be symmetrically moved to the original center O . The reason is that, unlike the SSH case, the high symmetry positions at the center O , or at any mirror line, are all connected. Thus, the difference between filling anomalies associated with boundary and bulk obstructions is that when the anomaly is resolved by identifying some of the edges, the former leads to completely trivial bulk phases, whereas the latter leads to topologically non-trivial bulk phases (e.g., obstructed atomic limits).

↦

One final remark concerns the effect of adding some extrinsic degrees of freedom at the edge of the DMQI. For example, we can add an SSH chain with N SSH lattice sites to one edge (say the one parallel to the x -axis), and its mirror image to the opposite edge. This process changes the total number of states by 2 N SSH . It also introduces 2 ν O , SSH states at the y = 0 mirror line, where ν O , SSH is the filling at the center of each SSH chain. This changes ν Γ as ν Γ → ν Γ +2 ν O , SSH , and as a result, if the SSH chain has a filling anomaly indicated by parity mismatch between N SSH and ν O , SSH , then its addition to the DMQI in the anomalous (non-trivial) phase will cancel the filling anomaly of the DMQI. Such a process will, however, induce a filling anomaly in the non-anomalous (trivial) phase. Thus, the effect of symmetrically adding SSH chains at the edge essentially exchanges what we would identify as the trivial and non-trivial phase, but keeps the distinction between them . As we will see later, this is a general feature for boundary-obstructed phases.

## E. Symmetry representations

In the previous discussion we clarified that a topological obstruction in the DMQI model exists if and only if the system has open boundaries. This statement can directly be reiterated as a statement about how the filled electrons transform under the space group G with periodic boundaries, versus how it transforms under the point group F with open boundaries. Namely, the representations ρ of G for the filled electron states corresponding to different parameters of the DMQI model are isomorphic with periodic boundary conditions, while upon opening boundaries, the filled electron states may transform under distinct, i.e., non-isomorphic, representations of the point group F , depending on a choice of boundary. In the following we explicitly show this symmetry distinction by computing the G and F representations for this model as a means to build intuition for the general formulation that will be presented in Sec. III B.

## 1. Equivalence between band representations in the bulk periodic system

Now let us show explicitly using symmetry representations that with periodic boundaries there are no distinct topological phases in the DMQI model. As described in detail in Refs. [40, 45, 55], as well as in Sec. III B, all atomic bands transform under representations of a space group G, which are obtained from local representations of the site-symmetry group G q of a maximal Wyckoff position. In the DMQI case, these correspond to the four maximal Wyckoff positions 1 a, 1 b, 1 c, and 1 d of multiplicity one whose site-symmetry group is the point group F . Although they have isomorphic site symmetry groups, the four position centers correspond to distinct coset decompositions of the space group under translations F = G/T , which implies the groups F 1 a , F 1 b , F 1 c and F 1 d are isomorphic but not related by conjugation under any element of G . The representations π of F in the DMQI model with periodic boundaries, are given by the two-dimensional representation ¯ E , for all choices of parameters, guaranteeing the anti-commutation of M x and M y .

To prove that all Wannier configurations in the DMQI are topologically equivalent with periodic boundaries it suffices to find a symmetry-preserving unitary transformation that acts in the space of filled bands and transforms between the distinct representations ρ of G . In momentum space, the diagonal entries of these matrix representations form the band representations ρ k . We discuss how to find ρ k in Sec. III B following the standard references [32, 40, 41, 45]. The possible representations induced by the local representation ¯ E at the four maximal Wyckoff positions Q have distinct momentum dependencies ρ Q k , which we will now tabulate. First, we see that, because of the π flux,

$$\rho _ { k } ( C _ { 2 z } ^ { 2 } ) = - 1,$$

for all Wyckoff positions. Next, let us focus on position 1 a and consider the representations of ρ k ( M x ) and ρ k ( M y ) .

Since this position is chosen as the center of the unit cell, the representation matrices ρ k will be k -independent. One possible choice for ρ 1 a k ( M x,y ) is ρ 1 a k ( M x )= σ 1 and ρ 1 a k ( M y )= σ 3 which corresponds to the eigenbasis of M y . Alternatively, we can choose the eigenbasis of M x : ρ ′ 1 a k ( M x ) = σ 3 and ρ ′ 1 a k ( M x ) = σ 1 . Clearly, the two BRs are equivalent with the momentum-independent Hadamard transformation U k = ( σ x + σ z ) / √ 2 relating the two.

Next, we can compare the band representations induced from distinct maximal Wyckoff positions, located at parallel high symmetry lines. Inducing a representation from 1 b we find that ρ 1 b k ( M x ) = σ 1 and ρ 1 b k ( M y ) = exp( ik y ) σ 3 . The unitary transformation that relates ρ 1 a k ( M y ) and ρ 1 b k ( M y ) is

$$U _ { k } = \frac { 1 } { 2 } \begin{pmatrix} 1 + \exp ( i k _ { y } ) & 1 - \exp ( i k _ { y } ) \\ 1 - \exp ( i k _ { y } ) & 1 + \exp ( i k _ { y } ) \end{pmatrix}. \quad ( 7 ) \quad \text{reser} \quad \text{this} \, 1 \\ \text{reser}$$

Finding the other unitary transformations follows analogously. Since the band representations are used to distinguish the different atomic limit phases, we conclude from this argument that they are all equivalent to the trivial atomic limit when the system has periodic boundary conditions. We remark that since the unitary transformation, acting on the filled bands alone, changes the nested polarization p x,y → p x,y +1 2 / , it explicitly highlights that these values cannot be naively used as a bulk topological invariant.

## 2. Symmetry representations of the open system

Now let us consider the open system without periodic boundary conditions. Assuming the system has open boundaries implies that the symmetry group is restricted from G to F , i.e., translation symmetry is no longer a symmetry of the system, and the number of high symmetry lines in the lattice is substantially reduced to two perpendicular lines crossing a single point O which is the center of the lattice. Above, in the layer construction discussed in Sec. II D, this symmetric region was labelled by Γ .

We want to consider restricting the representation ρ of the full space group G of the periodic system to a (reducible) representation π of the point group F of the open system. This implies selecting only group elements that belong to F , which is symbolized by π = ρ ↓ F . Thus, starting from a representation ρ q of the site symmetry group G q , we first generate a representation ρ of the full space group G which is then restricted to a representation π of the point group F of the open system. In the language of representations a boundary obstruction is the statement that different site symmetry representations ρ q which generate the space group representation ρ can generate different representations π of the open system despite being equivalent in the bulk periodic system (as we just showed above). This provides a symmetry-based characterization of boundary obstructions which we will illustrate for the DMQI below and discuss in more general terms in Sec. III B. In fact, we will prove a stronger statement in what follows: not only are the representations π of the two DMQI 'phases' distinct, but only one of them is a local repre- sentation of the point group F . We say a representation of the point group F is local if it can be generated by the action of F on localized atomic orbitals in real space. This means that it is generated as a sum of site symmetry representations of the Wyckoff positions in the open system (the Wyckoff positions of the open system are defined in analogy to the periodic system, as a set of points forming a closed orbit under the action of the point group F ). The notion of a local representation thus provides a definition for a symmetric atomic insulator in the open system.

Let us now see how this explicitly applies to the DMQI. The Wyckoff positions consist of the center O with multiplicity 1, the x and y axes each with multiplicity 2, and the general position with multiplicity 4. We begin by assuming that the representations generated from both 1 a and 1 d are local and show this leads to a contradiction. First note that restricting a representation from G down to F cannot change the dimension of the representation. Thus, since π 1 a and π 1 d are obtained from restricting equivalent bulk representations ρ to the point group F , we have

$$\text{tr} \, \pi ^ { 1 a } ( 1 ) = \text{tr} \, \pi ^ { 1 d } ( 1 ),$$

which uses the fact that the dimension of the representation is the trace of the representation of the group identity element. Second, note that the dimension of a symmetry representation acting on Wyckoff position Q is necessarily a multiple of its multiplicity, thus the contribution from the general Wyckoff position to tr π (1) mod 4 drops out. If the representation π is local, then it is a sum of representations acting on the different Wyckoff positions of the open system, and furthermore it implies that tr π (1) mod 4 receives contributions only from the Wyckoff positions lying on high symmetry lines, i.e.,the region denoted by Γ . As a result, we find

$$\nu _ { \Gamma } = \text{tr} \, \pi ( 1 ) \mod 4.$$

On the other hand, we can read off ν Γ from the site filling of the periodic system assuming vanishing polarization in the x and y directions (see Sec. II D). For definiteness, let us assume both N x and N y are odd such that the boundary of the open system relative to the center of the unit cell coincides with the edge of the unit cell. Clearly for ρ 1 a , there are two orbitals at the center of the unit cell which lie inside Γ leading to ν Γ = 2( N x + N y -1) = 2 mod 4 . On the other hand, for position 1 d , there are clearly no orbitals lying on Γ leading to ν 1 d Γ = 0 . Thus,

$$\nu _ { \Gamma } ^ { 1 d } \neq \nu _ { \Gamma } ^ { 1 a } \quad \quad ( 1 0 )$$

glyph[negationslash]

We can see that Eq. 8, 9, and 10 lead to a contradiction which implies that the assumption that both π 1 a and π 1 d are local must be wrong. For the chosen boundary, it is easy to see that π 1 a is local since none of the filled sites lie at the boundaries, thus we conclude that π 1 d cannot be local, i.e., it cannot be described by a symmetric atomic insulator with the same filling as the bulk, which precisely defines a BOTP with a filling anomaly. Once we have established the existence of two distinct phases in the open system, we can classify represen-

tations obtained from any other site symmetry representation ρ q which will be equivalent to either π 1 a or π 1 d depending on whether ν Γ = 2 0 , . These results are summarized in Table I, from which we see that π 1 d is distinct from π 1 a , π 1 b , π 1 c , which are equivalent to each other.

## 3. Wannier band representations

We have seen in Sec. II C that a natural diagnosis tool for boundary topology is the Wannier spectrum. The Wannier spectrum, with eigenstates spanning the same space as the periodic Hamiltonian, transforms under representations related to the periodic and open boundary representations ρ and π , respectively. By constructing hybrid Wannier functions localized along the a direction dual to b ( a b · = 2 π ), each individual band in the Wannier spectrum W b ( k ⊥ ) preserves a subgroup G b of the space group G . Therefore, we can associate a lower dimension band representation to these bands, identified by the symmetry character of the Wannier states at high symmetry momenta in the (Wannier) Brillouin zone transverse to b spanned by the momenta k ⊥ . The representation of each Wannier band w b a ( k ⊥ ) is related to the representation ρ of the periodic boundary Hamiltonian by restricting this representation to the group elements of G b ; this process is denoted as w b = ρ ↓ G b . Since the Wannier spectrum can be decomposed into a set of disconnected bands, we can additionally write w b = ⊕ a w b a , where a is a band index. The explicit form of the representation is given in further detail in Sec. III B. If we look at the Wannier band representations in different directions b such that G = ⋃ b G b , then one can unambiguously determine from which original bulk representation they are restricted, as well as the specific coset decomposition G/T that identifies the maximal symmetry point in the unit cell that serves as its center.

In the DMQI model, the four coset decompositions, labelled by the unit cell centers 1 a to 1 d , can therefore be distinguished by the Wannier band representations. In this case, the Wannier band representations are given by ρ restricted to the subgroups G x = { T , M x x } and G y = { T , M y y } . We can now illustrate how to identify the bulk representation π from the set of Wannier band representations. Let us choose the bulk Wannier functions to be located at each of the possible four centers, where they transform under the ¯ E representation of F = D π 2 . When the center is 1 a , the hybrid Wannier functions are constructed out of linear combinations of localized Wannier functions at the center of the unit cell in the perpendicular plane. For example, if we choose to keep G y , the two states in the unit cell have opposite eigenvalues of M , y and these eigenvalues have no momentum dependence. This is analogous to the SSH chain where the representation of the (energy) bands do not have a momentum dependence when the Wannier centers occupy the inversion center in the unit cell. The Wannier bands transform under the band representations [ w x ± ] k y ( M y ) = ± 1 . Alternatively, when keeping G x the hybrid Wannier functions transform under the band representation [ w y ± ] k x ( M x ) = ± 1 . If instead the bulk Wannier centers occupied the 1 b position, which is located at y = 1 2 /

in the unit cell, then the restriction of ρ to G y would reflect this by acquiring a phase of exp { ik y } under the action of M y . That is [ w x ± ] k y ( M y ) = ± exp { ik y } . Indeed, all four coset decompositions G/T correspond to 1 a, 1 b, 1 c, 1 d generate distinct Wannier band representations, and are shown in Table I. These representations exhaust the possible restrictions of ρ , and do not depend on the possible symmetric displacement of the Wannier centers away from the high symmetry points. We see that all bulk representations of this model are isomorphic, and the different Wannier representations reflect a choice of basis.

As we have seen in the previous section, the difference between the basis choice becomes important with open boundaries. Namely, it is important to know if the center of the unit cell occupies a high symmetry region Γ in the open boundary system. This depends on the choice of boundary. In Table I we compare the representations of the Wannier spectrum, to the filling ν Γ of the open boundary system, where the corner position coincides with the high symmetry position 1 d , and the center O with 1 a with periodic boundaries. This is the conventional boundary termination, and the direct comparison of the Wannier spectrum representations obtained from the Hamiltonian for a given set of parameters can be compared with this table to predict if π is a local representation of the open boundary symmetry group F characterized by a filling anomaly.

| ρ             | k x = 0 w y ( M x )   | k x = π w y +   | k y = 0 k w x + ( M y )   | y = π w x +   |   ν Γ |
|---------------|-----------------------|-----------------|---------------------------|---------------|-------|
| (1 a, ¯ E )   | + +                   | +               | +                         | +             |     2 |
| (1 b, ¯ E )   | +                     | +               | +                         | -             |     2 |
| (1 c, ¯ E )   | +                     | -               | +                         | +             |     2 |
| ( 1 d , ¯ E ) | +                     | -               | +                         | -             |     0 |

TABLE I. Diagnosis of boundary obstructions through Wannier band representations [ w b a k ] ( g ) for the band a = + (selected to have a positive eigenvalue at k = 0 ) in the DMQI model evaluated for the mirror operator preserved by the perpendicular direction to b at different high symmetry momenta. The bold lines correspond to the phases with a boundary obstruction when the conventional Wannier chemical potential is chosen, i.e., µ = (1 2 1 2) / , / (obstructed high symmetry site 1 d ), and the center of the sample O = (0 0) , (high symmetry site 1 a ). In the last column we show the character ν Γ which differentiates between the anomalous and the non-anomalous configurations with this boundary choice.

## III. BOUNDARY OBSTRUCTIONS: GENERALITIES

## A. General definition

After discussing the properties of the DMQI model in detail, our goal now is to introduce a general definition to capture topological distinctions encoded in a boundary obstruction rather than a bulk obstruction. As we have seen for the DMQI, the distinguishing feature between the two phases of the model on an open boundary is a filling anomaly manifesting as a fractional corner charge of e/ 2 , which is present in

one of the phases, but absent in the other. However, unlike a standard SPT, e.g., a higher-order TI, the boundary of the DMQI is not anomalous since the corner charges can be removed by the addition of two SSH chains on opposite edges while preserving the mirror symmetries. This means that we should be more careful when defining an edge termination to distinguish surface features which arise from the bulk compared to those arising from the edge termination itself.

## 1. Boundary termination

Let us begin by recalling some relevant concepts. Given a Hamiltonian H and a chemical potential µ , we can define the projector on the filled bands via

$$P \coloneqq p ( \mathcal { H }, \mu ) = \sum _ { n } \theta ( \epsilon _ { n } - \mu ) \left | \psi _ { n } \right \rangle \left \langle \psi _ { n } \right |, \ \mathcal { H } \left | \psi _ { n } \right \rangle = \epsilon _ { n } \left | \psi _ { n } \right \rangle. \quad \text{ically i} \\ \intertext { T w - s u s t w - r \colon s u s t w - s u s t w - s u s t w - s u s t w - s u s t w - s u s t w - s u s t w - s u s t w - s u s t w - s u s t w - s u s t w - s u s t w - s u s t w - s u s t w - s u s t w - s u s t w - s u s t w - s u s t w - s u s t w - s u s t w - s u s t w - s u s t w - s u s t w - s u s t w - s u s t w - s u s t w - s u s t w - s u s t w - s u s t w - s u s t w - s u s t w - s u s t w - s u s t w - s u s t w - s s t w - s u s t w - s u s t w - s u s t w - s u s t w - s u s t w - s u s t w - s u s t w - s u s t w - s u s t w - s u s t w - s u s t w - s u s t w - s u s t w - s u s t w - s u s t w - s u s t w - s u s t w - s u s t w - s u s t w - s u s t w - s u s t w - s u s t w - s u s t w - s u s t w - s u s t w - s u s t w - s u s t w - s u s t w - s u s t w - s u s t w - s u s t w - s u s t w - s u s s t w - s u s t w - s u s t w - s u s t w - s u s t w - s u s t w - s u s t w - s u s t w - s u s t w - s u s t w - s u s t w - s u s t w - s u s t w - s u s t w - s u s t w - s u s t w - s u s t w - s u s t w - s u s t w - s u s t w - s u s t w - s u s t w - s u s t w - s u s t w - s u s t w - s u s t w - s u s t w - s u s t w - s u s t w - s u s t w - s u s t w - s u s t w - s u s t w - s u t w - s u s t w - s u s t w - s u s t w - s u s t w - s u s t w - s u s t w - s u s t w - s u s t w - s u s t w - s u s t w - s u s t w - s u s t w - s u s t w - s u s t w - s u s t w - s u s t w - s u s t w - s u s t w - s u s t w - s u s t w - s u s t w - s u s t w - s u s t w - s u s t w - s u s t w - s u s t w - s u s t w - s u s t w - s u s t w - s u s t w - s u s t w - s u s t w - s u s t u s t w - s u s t w - s u s t w - s u s t w - s u s t w - s u s t w - s u s t w - s u s t w - s u s t w - s u s t w - s u s t w - s u s t w - s u s t w - s u s t w - s u s t w - s u s t w - s u s t w - s u s t w - s u s t w - s u s t w - s u s t w - s u s t w - s u s t w - s u s t w - s u s t w - s u s t w - s u s t w - s u s t w - s u s t w - s u s t w - s u s t w - s u s t w - s u s t w - s u S t w - s u s t w - s u s t w - s u s t w - s u s t w - s u s t w - s u s t w - s u s t w - s u s t w - s u s t w - s u s t w - s u s t w - s u s t w - s u s t w - s u s t w - s u s t w - s u s t w - s u s t w - s u s t w - s u s t w - s u s t w - s u s t w - s u s t w - s u s t w - s u s t w - s u s t w - s u s t w - s u s t w - s u s t w - s u s t w - s u s t w - s u s t w - s u s t w - s u s t t w - s u s t w - s u s t w - s u s t w - s u s t w - s u s t w - s u s t w - s u s t w - s u s t w - s u s t w - s u s t w - s u s t w - s u s t w - s u s t w - s u s t w - s u s t w - s u s t w - s u s t w - s u s t w - s u s t w - s u s t w - s u s t w - s u s t w - s u s t w - s u s t w - s u s t w - s u s t w - s u s t w - s u s t w - s u s t w - s u s t w - s u s t w - s u s t w - s u s t W - s u s t w - s u s t w - s u s t w - s u s t w - s u s t w - s u s t w - s u s t w - s u s t w - s u s t w - s u s t w - s u s t w - s u s t w - s u s t w - s u s t w - s u s t w - s u s t w - s u s t w - s u s t w - s u s t w - s u s t w - s u s t w - s u s t w - s u s t w - s u s t w - s u s t w - s u s t w - s u s t w - s u s t w - s u s t w - s u s t w - s u s t w - s u s t w - s u s t w -s u s t w - s u s t w - s u s t w - s u s t w - s u s t w - s u s t w - s u s t w - s u s t w - s u s t w - s u s t w - s u s t w - s u s t w - s u s t w - s u s t w - s u s t w - s u s t w - s u s t w - s u s t w - s u s t w - s u s t w - s u s t w - s u s t w - s u s t w - s u s t w - s u s t w - s u s t w - s u s t w - s u s t w - s u s t w - s u s t w - s u s t w - s u s t w - s u s t w - s us t w - s u s t w - s u s t w - s u s t w - s u s t w - s u s t w - s u s t w - s u s t w - s u s t w - s u s t w - s u s t w - s u s t w - s u s t w - s u s t w - s u s t w - s u s t w - s u s t w - s u s t w - s u s t w - s u s t w - s u s t w - s u s t w - s u s t w - s u s t w - s u s t w - s u s t w - s u s t w - s u s t w - s u s t w - s u s t w - s u s t w - s u s t w - s u s t w - s u s t$$

The projector P is well-defined if and only if there are no eigenstates of H at the chemical potential, i.e., only if H is gapped at µ . Note that the definition above does not assume translation symmetry. Since we are interested in topological properties rather than dynamics, we will consider two pairs ( H , µ ) and ( H ′ , µ ′ ) equivalent if p ( H , µ ) = p ( H ′ , µ ′ ) . In the following, unless otherwise stated, we will choose µ = 0 such that p ( H ) = p ( H , 0) . Given an operator O ( H or P ), we can write it in the position basis as

$$O ( r, r ^ { \prime } ) \coloneqq \langle r | \hat { O } | r ^ { \prime } \rangle \quad \quad ( 1 2 ) \quad \ \ \ \ \.$$

For example, the position-dependent projector is given by

$$P ( r, r ^ { \prime } ) = \sum _ { n } \theta ( \epsilon _ { n } - \mu ) \psi _ { n } ( r ) \psi _ { n } ^ { \dagger } ( r ^ { \prime } ), \quad ( 1 3 ) \quad \text{men} \quad \text{pote} \quad \text{ }$$

which is a matrix-valued function of two spatial positions r and r ′ (for a model with N o orbitals, this will be an N o × N o matrix).

Westart by considering a periodic system with symmorphic space group G decomposed as a direct product of a group of translations T, and a point group F relative to a fixed point O . An F -symmetric boundary in d -dimensions is defined by specifying a ( d -1) -dimensional surface Σ dividing the space into an interior region Λ in and an exterior region Λ out with a 'thickness" parameter ξ such that for distances larger than ξ boundary effects are negligible. We can now define a surface termination as a continuous map σ which takes a translationally symmetric Hamiltonian H (with the corresponding projector P = p ( H ) ) to a Hamiltonian in the open system H obc (with the corresponding projector P obc = p ( H obc ) ), such that P obc ( r r , ′ ) = P ( r r , ′ ) deep inside the sample, and P obc ( r r , ′ ) = P 0 ( r r , ′ ) far outside, with P 0 denoting some reference "trivial" projector. For simplicity, we will take P 0 to denote projector where all states are empty. What we mean by deep inside and far outside the sample can be made precise by defining

$$d _ { r, \Sigma } = \min _ { r ^ { \prime } \in \Sigma } | r - r ^ { \prime } |, \quad \quad ( 1 4 ) \quad \text{ a}$$

which simply denotes the shortest distance from a point r to the surface Σ . Thus, we can define the most general boundary termination as

$$\text{we} \quad \text{$\mathcal{H}_{obc}=\sigma(\mathcal{H})$ such that for $d_{r,\Sigma}$, $d_{r^{\prime},\Sigma} \gg \xi$} \\ \text{m-} \quad \begin{cases} P _ { \text{obc}(r,r^{\prime})=P(r,r^{\prime})} & \colon r,r^{\prime}\in\Lambda_{in}$} \\ P _ { \text{obc}(r,r^{\prime})=P_{0}(r,r^{\prime})} & \colon r,r^{\prime}\in\Lambda_{out}$} \. \ \ (15 ) \\ P _ { \text{obc}(r,r^{\prime})=0 & \colon \text{otherwise} \end{cases}$$

Here, we also made the assumption that P ( r r , ′ ) vanishes when one of r r , ′ is deep inside and the other is far outside the sample. This means that there are no states which 'propagate' between the inside and outside of the sample.

We notice that in practice, boundary terminations are typically implemented by a sharp edge where the lattice terminates, or more generally by a potential well which separates the filled states inside the sample from the empty states outside. An example of the latter is ultracold atoms [56] where the lattice termination is usually implemented by a potential trap which confines the particles within a certain region. This motivates the consideration of a particularly simple boundary termination implemented by making the chemical potential spatially dependent such that it lies in the gap inside the sample, and below the smallest energy eigenvalue outside. This is equivalent to choosing the open system Hamiltonian considered in Ref. [46]

$$( 1 2 ) \quad \mathcal { H } _ { \text{obc} } = \sigma ( \mathcal { H } ) = p ( \mathcal { H } ) \phi ( r ) p ( \mathcal { H } ) + \phi _ { \text{out} } ( 1 - p ( \mathcal { H } ) ), \ ( 1 6 )$$

where φ ( r ) = φ in &lt; 0 deep inside the sample, and φ ( r ) = φ out &gt; 0 far outside. 9 This boundary termination, implemented by changing a single scalar parameter (the chemical potential), is the simplest way to interpolate between the bulk Hamiltonian and a fixed trivial reference without introducing any extra degrees of freedom at the boundary. It includes most surface terminations needed in practice, e.g., a linear potential, a sharp edge, etc. and has the advantage of connecting directly to the Wannier spectrum.

To make the notion of boundary obstructions introduced in the next section as general as possible, we will also find it useful to consider an augmented class of boundary terminations by allowing for the addition of gapped degrees of freedom at the boundary. Such additional degrees of freedom do not induce any additional boundary gap-closing, thereby preserving the notion of boundary obstructions. To define this more general class of boundary terminations rigorously, we note that since the open system Hamiltonian for any boundary termination is equivalent to (16) away from the boundary, we can write the Hamiltonian for the most general boundary termination we consider as

$$P _ { 0 } \quad \sigma ( \mathcal { H } ) = p ( \mathcal { H } ) \phi ( r ) p ( \mathcal { H } ) + \phi _ { \text{out} } ( 1 - p ( \mathcal { H } ) ) + \mathcal { H } _ { \Sigma }, \ \ ( 1 7 )$$

9 Wewill usually also impose the condition that φ ( r ) = 0 only at the boundary surface Σ , but this is not crucial for what follows.

where H Σ vanishes away from the boundary. H Σ is defined on a ( d -1) -dimensional shell which is topologically equivalent to S d -1 . We say a boundary termination σ is trivial if there exists φ ( r ) such that H Σ defined by (17) is symmetric and gapped for all H . A trivial termination allows us to compare distinct Hamiltonians with open boundaries by excluding terminations which introduce low energy (gapless) modes at the boundary. Such a restriction is necessary since BOTPs are extrinsic HOTIs, hence their surface states can always be removed through the addition of boundary degrees of freedom and it enables a direct correspondence between the bulk and boundary degrees of freedom.

## 2. Boundary obstructed phases

We are now ready to define the notion of boundary obstructions. Given a trivial termination σ , we say there is a boundary obstruction between two gapped translationally symmetric Hamiltonians H 1 , H 2 with periodic boundary conditions if:

- (i) there is a smooth trajectory in the space of symmetric gapped Hamiltonians H ( ) t for t ∈ [0 , 1] connecting H 1 and H 2 such that H (0) = H 1 and H (1) = H 2 ;
- (ii) the trajectory induced by the boundary termination σ ( H ( )) t necessarily involves closing the energy gap at a high-symmetry surface, i.e. there exists a t ∈ [0 , 1] such that σ ( H ( )) t has a zero energy (remember we use a convention that sets the chemical potential to 0) eigenstate localized in the direction perpendicular to a high symmetry surface.

A high symmetry surface (HSS) in d dimensions denotes any D -dimensional hyperplane ( D &lt; d ) on the boundary that is left invariant by at least one non-trivial point group symmetry operation in F . For example, in the simple rectangular termination considered for the DMQI, both edges are high symmetry surfaces since each one is invariant under one of the mirror symmetries. For more complicated boundaries, with arbitrary shape, we can still define high symmetry surfaces by first 'rounding the corners" to make the curvature finite everywhere, and then considering the tangent plane to the boundary at any given point. We then define the surface spectrum as the energy spectrum of the states which are localized in the direction perpendicular to the surface. The surface spectrum is a function of the momentum parallel to the surface which is a vector in the tangent plane at this point [23, 25, 33, 57]. See Appendix D for an illustration of the edge spectrum for a general boundary termination of the DMQI.

We can now define boundary-obstructed phases as equivalence classes of Hamiltonians under boundary obstructions defined earlier. That is, we consider a space of gapped Hamiltonians in the translationally-symmetric system Ξ pbc , which we take to be topologically trivial, i.e., all gapped phases can be continuously deformed to each other without closing the gap or breaking the symmetry. Given a trivial boundary termination σ , we define the group of boundary-obstructed topological phases to be the group of equivalence classes of boundary-terminated Hamiltonians σ (Ξ pbc ) under continuous deformations which do not close the gap at any high symmetry surface. Intuitively, this definition distinguishes different classes of Hamiltonians based on boundary obstructions, and enables a finer distinction between phases that are topologically identical from the perspective of bulk topology.

The above definition introduces a relative distinction between a pair of Hamiltonians given a certain trivial boundary termination. Such distinctions are generally termination dependent even for trivial terminations. For example, recall from the analysis of the DMQI in the previous section that the termination (which determines the Wannier chemical potential) determines which Hamiltonians will exhibit a filling anomaly. On the other hand, for any symmetric termination, there are always two distinct phases in the open system distinguished by the Z 2 filling anomaly. This suggests that although the boundary determines which BOTP class a given Hamiltonian ends up in, the topological classification (homotopy group) of BOTPs itself is independent of the boundary termination σ for sufficiently general topologically trivial boundaries. Here, by 'sufficiently general", we mean that the distinct topological sectors are stable to symmetry-preserving perturbations of the boundary which keep it gapped and Wannierizable. In particular, this includes the possibility of adding degrees of freedom on the surface which can trivialize the surface bands if they possess fragile topology [58]. While a general proof of this statement is beyond the scope of this work, we can directly verify it for all the examples of BOTPs considered in this work (the DMQI considered in Sec. II D and the 3D examples of Sec. V) by showing that the different BOTPs can be distinguished by a physical signature: either filling anomalies or gapless surface states, on any symmetric boundary.

In general, boundary obstructions may co-exist with bulk obstructions. For instance, a non-trivial bulk phase, such as a 3D higher-order topological insulator with chiral hinge states, can have several distinct patterns of hinge state. These patterns effectively differ by the addition of a BOTP, i.e., we can always move between such different surface configurations by adding a BOTP. Such a BOTP, which is manifestly trivial in the bulk, can be identified with the difference between the Hamiltonians corresponding to the two distinct patterns of surface states. Thus, when studying boundary obstructions, we can restrict ourselves to phases that are completely trivial in the bulk.

One important aspect in our definition is that only gapclosings at high-symmetry surfaces are relevant. Physically, this is a natural restriction since most simple realistic boundaries tend to be along high symmetry surfaces. One way to see this is by recalling the discussion of Sec. II D, in which the filling anomaly in the DMQI was diagnosed by the mismatch between the total filling and the number of states lying at the two mirror lines modulo 4. Thus, an edge gap-closing can change the filling anomaly by changing the number of states at the mirror lines if and only if this gap-closing occurs at one of the mirror invariant points at the edge. This means that our definition captures the same topological distinction described by the filling anomaly for the DMQI.

Another way to see the importance of distinguishing high

FIG. 10. Illustration of zero-energy states (open circles) for the chiral-symmetric DMQI model on a circle geometry. Zero energy states are localized on domain walls for an edge mass term that is invariant under the two mirror symmetries, but has a different sign at the x = 0 and y = 0 lines, thus changing sign four times along the edge. These zero-energy states can only be removed by annihilating pairwise at the high symmetry points denoted by the x symbol.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00011v2-figure-8c54929e-3cd7-46ce-bf67-a67c337cd0e7.jpg)

symmetry points/surfaces versus generic points on the surface is to consider the DMQI on a disk geometry whose boundary is a circle. If chiral/particle-hole symmetry is unbroken, the non-trivial phase of the model is characterized by four symmetry-related zero-energy states. These states can be associated with domain walls for a surface Dirac mass [22, 23, 25, 33, 35, 49, 59] that is invariant under both mirrors and changes sign four times as we go around the circle (cf. Fig. 10) as shown in Appendix D. In this case, we can see that the gap at any generic point on the edge can be closed by simply moving one of the domain walls through this point. Such a gap-closing cannot be associated with a topological distinction since that would imply the existence of uncountably-many distinct topological phases associated with all possible (continuously varying) positions of the domain walls. However, by restricting ourselves to gap-closing transitions occurring at the HSSs, i.e., at x = 0 or y = 0 , we can distinguish between the cases with and without corner states. This follows because the only way to get rid of the zero-energy domain wall states while preserving symmetry is by annihilating them pairwise at one of the mirror invariant points at the boundary.

## 3. Relationship to higher-order topological insulators

It is instructive to relate our understanding of the boundary obstruction in the DMQI in terms of surface domain walls to the corresponding understanding in higher-order topological insulators pioneered in Refs. [22, 23, 25, 49]. In these works, the surface states of higher-order topological insulators were understood in terms of symmetry-enforced surface domain walls where the transformation properties of the surface Dirac mass under spatial symmetry forced it to change sign. In the DMQI, the surface mass does not change sign under either mirror symmetry (see Appendix D). However, by requiring the edge spectrum to be gapped at HSSs, we can distinguish phases based on the sign of the mass term at these points which are associated with different patterns of surface domain walls. For instance, if the sign of the mass term at x = 0 differs from that at y = 0 , then there are four domain walls as we go around the circle (cf. Fig 10), which cannot be removed without changing the sign of one of these masses. Thus, we can understand the localized boundary modes of the DMQI as domain walls protected by the gap at HSSs. This is verified by an explicit calculation of the surface states of the DMQI in Appendix D. We note that the transformation properties of the surface mass term under the symmetries does not force it to change sign. Instead, two distinct, symmetryallowed configurations of mass signs at HSSs are possible, but they cannot be deformed to each other (without closing a gap at an HSS). This is in contrast to higher-order topological insulators where the sign change of the mass is enforced by the symmetry.

It is worth noting that our topological distinctions are similar to a real-space version of those used in diagnosing semimetals whose spectrum is gapped at high-symmetry momenta [31, 60]. For instance, a Weyl semimetal with two inversion-related Weyl points with opposite chirality can be smoothly deformed to a trivial insulator by bringing the two Weyl points together. However, this process will necessarily involve a gap-closing at one of the inversion-invariant momenta [61, 62].

Another related concept is that of the so-called extrinsic higher-order topological insulators discussed in Refs. [13, 24]. These are insulators hosting 'higher-order" surface states with co-dimension ( d -D ) ( D &gt; 1 ) that can be removed by symmetrically adding a lower-dimensional SPT on the boundary. Our boundary-obstructed phases fall under this definition since their surface states (or filling anomalies) can be removed by adding some SPTs on the boundary while preserving the symmetry. For example, the corner modes (and corner charge) in the DMQI can be removed by symmetrically adding a pair of SSH chains at two symmetry-related edges as discussed in Sec. II D. However, for a BOTP such a process only redefines what we identify as the trivial phase, and does not remove the distinction between different phases. In other words, the relative distinction between the phases is insensitive to what we add at the boundary as long as it is kept fixed when comparing the two phases, a requirement that is already encoded in our definition. Thus, we only consider boundary transitions driven by changes in the bulk, ruling out the cases in which a non-trivial SPT is glued to the surface of a trivial bulk.

In general, relative topological distinctions are captured by a mathematical structure called a torsor (rather than a group) [63]. A torsor can be thought of as a group without a clear notion of an identity element (trivial phase). We notice that this concept is not new as many topological distinctions for well-known bulk phases are also relative. This is the case, for example, for obstructed atomic phases where an arbitrary choice is made for the trivial atomic limit [31, 32]. A simple example is provided by the 1D SSH chain which has two different values of the polarization that can only be distinguished relative to each other. Such a relative distinction can be made more absolute by taking into account the background of positive ions or fixing a convention for the unit cell, either of

which can be used to distinguish the two polarization states in the SSH chain according to whether the charge centers lie on top of the positive ions (center of the unit cell) or not (edge of the unit cell) 10 . It is, however, important to keep in mind that the topological distinctions in the electronic Hamiltonian (which does not contain information about the positive ions) in these cases are strictly relative. A similar concept applies to the topological distinction provided by BOTPs.

## B. Band representations of BOTPs

In this section, we present an alternative formulation of BOTPs based on the formalism of band representations [55, 64] (see Appendix C for a brief review). A band representation corresponds to a set of bands (over the usual momentum space Brillouin zone) generated by decorating the real space Wyckoff positions by atomic orbitals that transform in some representations of the local site symmetry group. In short, a band representation is generated by the Fourier transform of a real space representation of the space group G generated by these local atomic orbitals. One of the recent developments of band theory [31, 55] was the classification of topological bands that identifies topological bands as any isolated set of bands that cannot be written as a band representation, i.e., cannot be generated from the symmetry action on real space orbitals. In the following, we will show these ideas can be extended to describe BOTPs.

As discussed in Sec. III A, we can restrict ourselves to BOTPs which admit a Wannier representation. This means that such BOTPs correspond to band representations of the space group G generated from real space atomic orbitals in the periodic system. When we discussed the notion of boundary obstructions in the case of the DMQI in Sec. II, we found that there were two distinct notions of obstruction on a given boundary. First, there is a relative notion which distinguishes phases which are equivalent in the bulk, but which are separated by edge-gap-closing transitions at the boundary. This notion does not specify which of the phases is trivial. Second, there is an absolute notion of obstruction identified by the presence of a filling anomaly which distinguishes the trivial phase (no filling anomaly) from the obstructed or non-trivial one (filling anomaly). In the following, we will use band representations to also define two notions of boundary obstructions.

Our main approach to define boundary obstructions is to identify band representations which do not describe atomic insulators in the open system. That is, a BOTP corresponds to an atomic insulator of the periodic system that cannot be written in terms of localized atomic orbitals which transform under the symmetry group F of the open system. This defines an absolute notion for boundary obstruction (once a boundary is fixed) similar to the notion of filling anomalies discussed

10 Wenote that there is still an integer ambiguity in defining polarization coming from the freedom to assign the electron position to any unit cell.

in Sec. II D. On the other hand, we can also define boundary obstructions in a relative way similar to Sec. II C such that two Hamiltonians are considered boundary obstructed if they belong to equivalent band representations of the space group in the bulk, but distinct representations of the point group in the open system. We will show below how these two notions - the absolute and relative - are related.

## 1. Symmetry obstruction with open boundaries

A band representation is specified by a set of Wyckoff positions and a set of irreps attached to each position. Any such band representation can be decomposed into a direct sum of so-called elementary band representations (EBRs). An EBR is a BR generated from a single site in a maximal Wyckoff position q [65], and will be the focus of what follows. Given a symmetry representation ρ q of the site symmetry group G q ⊂ G which leaves the site q invariant, we can generate an EBR ρ for the Wyckoff position invariant under the space group G by a procedure called induction denoted by ρ = ρ q ↑ G . This procedure essentially extends the symmetry action from the site symmetries of the Wyckoff positions in the unit cell to the space group [55, 64] (see appendix C for details). Note that band representations induced from different positions q and q ′ may turn out to be equivalent. For instance, we saw this to be the case in the DMQI model where the BRs generated by placing the two electrons in the ¯ E representation at each of the four maximal Wyckoff positions 1 a, . . . , 1 d turned out to be the same. In general, two BRs ρ and ρ ′ are considered equivalent if there exists a symmetry preserving interpolation matrix S t ( ) such that S (0) = ρ and S (1) = ρ ′ [55, 64]. The matrix function S t ( ) parameterizes the adiabatic paths between phases.

The crucial observation is that two representations ρ and ρ ′ which are equivalent in the periodic system may not be equivalent in the presence of a boundary. To see this, we note that the representation of the open system is obtained by restricting the representation of the space-group (i.e., the periodic system) to a representation of the point group F ⊂ G (i.e., the open system) relative to a global origin O for the point group. 11 Such a restriction, or 'subduction", which generates a representation π of the open system is denoted by π = ρ ↓ F . Thus, two open boundary representations π and π ′ which are generated from equivalent bulk representations may in fact not be the same depending on the original generating site q .

Let us now define the notion of boundary obstructions by employing the machinery of band representations reviewed in App. C. First, we say that a site symmetry representation generated from site q in the periodic system defines a (absolute) boundary obstructed phase in the open system with point

11 The explicit choice of boundary termination plays a role in the determination of the symmetry origin O . Hence the ambiguity in the definition of a BOTP due to the choice of boundary termination is implicitly embedded in the determination of the point-group relative to the origin O .

group F if π = ( ρ q ↑ G ) ↓ F is not a local representation of F . Here, a local representation of F means that it can be written as a sum of site symmetry representations of the Wyckoff positions of the open system. Alternatively, we can define a relative boundary obstruction between two phases described by induced representations of ρ q and ρ q ′ on a given open boundary with point group F if

$$\rho = ( \rho _ { q } \uparrow G ) \ \sim \ \rho ^ { \prime } \ = \ ( \rho _ { q ^ { \prime } } \uparrow G ) \quad \ \ ( 1 8 ) \quad \text{too\ } \text{a}$$

and

$$\pi = ( \rho _ { q } \uparrow G ) \downarrow F \ \rtimes \ \pi ^ { \prime } = ( \rho _ { q ^ { \prime } } \uparrow G ) \downarrow F \quad ( 1 9 ) \quad \text{trum.} \\ \text{ator} \,.$$

where the symbols ∼ and glyph[notsimilar] indicate equivalence or inequivalence of the representations as defined above.

The relative and absolute notions of boundary obstruction can be tied together by noting that a relative obstruction between ρ q and ρ q ′ implies that at least one of them is boundary obstructed in an absolute sense. In other words, there can never be a relative boundary obstruction between ρ q and ρ q ′ if both π = ( ρ q ↑ G ) ↓ F and π ′ = ( ρ q ′ ↑ G ) ↓ F are local and ρ = ρ q ↑ G and ρ ′ = ρ q ′ ↑ G are equivalent. To see this, we note that, by definition, we can decompose a local representation as a sum of irreps of the point group F as

$$\pi = \sum _ { l } n _ { l } ( \pi ) \pi _ { l } \quad \quad ( 2 0 ) \quad \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \
 \end{cases}$$

where l runs over the irreps of F, and n l ( π ) denotes the multiplicity of π l in the expansion of π . Here, we assumed that we moved all the electrons in the open system symmetrically to the center O which can always be done if π is local. The multiplicities can be obtained by employing the orthonormality of the irrep characters relative to the scalar product [66]:

$$\langle \pi _ { l } | \pi _ { m } \rangle = \frac { 1 } { | F | } \sum _ { g \in F } \text{tr} \, \pi _ { l } ( g ) \text{tr} \, \pi _ { m } ^ { * } ( g ) = \delta _ { l, m } \quad ( 2 1 ) \quad \text{plane} \quad \text{the s} _ { l } \quad \\ \text{tersec}$$

leading to

$$n _ { l } ( \pi ) = \frac { 1 } { | F | } \sum _ { g \in F } \text{tr} \, \pi ( g ) \text{ tr} \, \pi _ { l } ^ { * } ( g ). \quad \text{ (22)} \quad \text{inter} \, \\ \text{For} \, \text{defi}$$

Wenote now that for g ∈ F , tr π g ( ) = tr ρ g ( ) which by virtue of the equivalence ρ ∼ ρ ′ is equal to tr ρ g ( ′ ) = tr π ′ ( g ) leading to n l ( π ) = n l ( π ′ ) . That is, the two local representations π and π ′ are equivalent. This means that among the equivalence classes of representations of the open system induced from equivalent bulk representations, there is only one class containing the local representations, and they are all equivalent. This class can then be labelled as the 'trivial' phase. Any other phase with a boundary obstruction relative to it will also have an absolute obstruction.

## C. Diagnosing BOTPs using the Wannier spectrum

Since the BOTP classification is independent of the spe- cific termination (as long as it is topologically trivial), we are going to choose a particular termination for which the calculation of the boundary spectrum simplifies. This corresponds to choosing φ ( r ) in Eq. 16 to be a linearly increasing potential such that the boundary spectrum is the same as the Wannier spectrum, as discussed in Sec. II C [46]. We can then apply band representation theory to the Wannier spectrum to diagnose BOTPs in a periodic geometry without having to resort to an open system.

Let us begin by recalling some basic facts about the Wilson loop operator whose eigenvalues are the Wannier spectrum. As we briefly introduced before, the Wilson loop operator W b is a unitary operator defined by parallel transporting the occupied-band projection operator in a closed path in the Brillouin zone along the reciprocal lattice vector b . Its eigenvalues have the form exp 2 { πiν b ( k ) } , and its eigenstates are the hybrid Wannier functions that are localized in the direction along the lattice vector a which is dual to b , (i.e., a · b = 2 π ) for a fixed momentum along the perpendicular directions k ⊥ [36, 37]. The dispersion of the Wannier eigenvalues ν b ( k ⊥ ) with k ⊥ is known as the Wannier spectrum, and provides information about the k -resolved center of the charge along the a direction, for k ⊥ momenta perpendicular to a a , · k ⊥ = 0 . The Wannier spectrum is periodic, reflecting the translation symmetry of a unit cell along a . The lattice constant is here set to 1.

To see how boundary obstructions can be diagnosed using the Wannier spectrum, we note that since we only care about gaps at HSSs, we can restrict ourselves to boundary terminations consisting of intersections of ( d -1) -dimensional high-symmetry hyperplanes to study a d -dimensional BOTP (a high-symmetry surface hyperplane is one which is left invariant by at least one non-trivial spatial symmetry). We can then study the surface spectrum at different surface hyperplanes, or at the intersection of any number of such hyperplanes. For example, in three dimensions, we can consider the spectrum on 2D surfaces or on 1D hinges lying at the intersection of two such surfaces. In the following, we focus on boundary obstructions accompanied by a gap-closing at a surface hyperplane (i.e., surface obstructions) rather than the intersection of surface hyperplanes (i.e., hinge obstructions). For these types of geometries, a boundary obstruction in our definition is equivalent to the statement that any trajectory connecting the two Hamiltonians must close the surface gap for at least one high-symmetry surface hyperplane.

The key concept we apply is that the spectrum of a surface hyperplane perpendicular to a lattice vector a for a linear potential is the same as the Wannier spectrum ν b ( k ) [46]. This can be used to define a Wannier chemical potential µ b (cf. Sec. II C) at which the Wannier spectrum ν b ( k ) should remain gapped if the surface energy gap is to remain open. Thus, we can identify actual surface gap-closing transitions in a surface hyperplane perpendicular to the vector a with gapclosing transitions in the Wannier spectrum ν b ( k ) at the WCP µ b . Unless otherwise stated, we choose a convention that the boundary termination is consistent with (our choice of) the unit cell. The boundaries of the unit cell are chosen to be at 1 / 2 in units of the primitive lattice vectors, which corresponds

to a Wannier chemical potential µ b = 1 2 / for any primitive lattice vector a .

The discussion above provides a recipe for diagnosing boundary obstructions by using Wannier spectra as follows. Given two bulk Hamiltonians H 1 and H 2 , a set of high symmetry surface hyperplanes specified by a set of reciprocal lattice vectors { b } , and WCPs { µ b } describing the positions of these planes within the unit cell, then a boundary obstruction exists if every trajectory connecting H 1 and H 2 involves a gap closing in the Wannier spectra ν b ( k ) for some direction b at the Wannier chemical potential µ b .

## 1. Wannier band representations

A gap closing of the Wannier spectrum indicates a change in the symmetry representations or topology of the Wannier bands. Therefore, it is convenient to relate the symmetry representations of the Wannier spectrum to the periodic boundary and open boundary representations discussed in Sec. III B. The Wannier spectrum is obtained by the eigenstates of the Wilson loop operator W b ( k ⊥ ) , whose eigenstates are hybrid Wannier functions ∣ ∣ h b glyph[lscript]n ( k ⊥ ) 〉 (HWF), localized in the a -direction at a layer index glyph[lscript]. The HWFs satisfy W b ( k ⊥ ) ∣ ∣ h b glyph[lscript]n ( k ⊥ ) 〉 = exp 2 { πiν b n ( k ⊥ ) } ∣ ∣ h b glyph[lscript]n ( k ⊥ ) 〉 , and ∣ ∣ h b glyph[lscript]n ( k ⊥ ) 〉 = (1 2 / π ) ∫ 2 π 0 dk b exp {-iglyph[lscript]k b · a } U nm | ψ m ( k ) 〉 where U nm is fixed by the Wilson loop operator and reflects the gauge freedom in the choice of these functions.

Since the HWFs are just a linear combination of the occcupied states (i.e., a basis transformation of the occupied states), then the collection of the hybrid Wannier functions transforms under the same symmetry representation ρ as the collection of occupied Bloch states | ψ n ( k ) 〉 . The translation symmetry along a connects the HWFs in distinct layers glyph[lscript]. Thus we want to ask under what representation does the collection of HWFs transform at a given glyph[lscript] . Here this representation is denoted w b , and it is a representation of the subgroup G b ⊂ G that leaves glyph[lscript] invariant. Hence, w b is obtained by restricting ρ to sites in the layer glyph[lscript] , which can be expressed conveniently by w b = ρ ↓ G b . This restriction amounts to disregarding all entries of ρ acting on basis states outside of glyph[lscript] and group elements outside of G b . Alternatively, for atomic bands we can also obtain w b by induction from the site symmetry representation ρ q of a site q in this layer. These relations are symbolically written as

$$\rho = \rho _ { q } \uparrow G, \ \ w ^ { b } = \rho _ { q } \uparrow G _ { b }, \ \ w ^ { b } = \rho \downarrow G _ { b }. \quad ( 2 3 ) \quad \text{We a} \\ \text{throu} \quad \omega \ \omega$$

Note that due to translation symmetry, all layers transform under the same band representation, and changing the layer corresponds to changing the origin of G b . Since each glyph[lscript] preserves translations perpendicular to a , we can express the Wannier band representation (WBR) as a function of k ⊥ , and it is explicitly given by

$$w _ { k _ { \perp } } ^ { b } ( h ) = \sum _ { i j } \exp \{ i ( h _ { i } - h _ { j } ) q \cdot k _ { \perp } \} \rho _ { q } ( h _ { j } ^ { - 1 } h h _ { i } )$$

where h ∈ G , b and the summation is taken over the coset representatives h i ∈ G /G b q that generate all sites q i = h i q in the layer. This formula is analogous to the induction formula for the bulk representation

$$\rho _ { k } ( g ) = \sum _ { i j } \exp \{ i ( g _ { i } - g _ { j } ) q \cdot k \} \rho _ { q } ( g _ { j } ^ { - 1 } g g _ { i } ),$$

with g -1 j gg i ∈ G q , g ∈ G, and g i ∈ G/G q with a derivation presented in App. C to be self-contained.

With this structure in place, it is now crucial to notice that the Wannier band representation w b may correspond to disconnected Wannier bands,

$$w ^ { b } = \bigoplus _ { a } w _ { a } ^ { b },$$

labelled by an index a . In this case, the representations w b a may have their own topological characterization, and they are not required to be local even when the full bulk representation ρ is local. Indeed each w b a either corresponds to an elementary WBR, or to a topological band when its characters at high symmetry momenta are not compatible with a local atomic description. Reiterating, while w b in total forms an EBR, each individual component describing separate Wannier spectrum bands, w b a , does not need to be a an EBR. An example of topological Wannier bands is shown in Sec. V C. Such phases with topological Wannier bands result in anomalous boundary modes rather than the boundary charges that would appear if the Wannier bands were in obstructed atomic limit representations.

The consistency relations above allow us to infer the bulk representation ρ and the open boundary representation π from the different Wannier representations w b . Closing the Wannier spectrum implies a change in the Wannier representations w b a . The explicit calculation of the Wannier representations together with the bulk and open system representations are shown in Tables. I, II and III for the models discussed in this text.

## IV. BOTP AND FILLING ANOMALIES IN OTHER 2D SYSTEMS

Now that we have laid out the general framework for BOTPs we would now like to discuss the possibility of finding BOTP phases in other 2D systems besides the DMQI model. We are going to focus on BOTPs which can be identified through a filling anomaly in the open system similar to the DMQI (Sec. II D).

## A. Generalities on filling anomalies

Let us begin by distinguishing the spatial symmetry groups in a few different settings. First, we can consider an infinite translationally symmetric system with a space group G .

Among the possible wallpaper groups [67], we restrict ourselves to symmorphic ones which can be written as a direct product of a point group G relative to a point O and the group of translations in the infinite system. We can instead consider a finite translationally symmetric system with periodic boundary conditions such that for translation along a given lattice vector T a , there is an integer such that T n a = 1 , i.e., the group of translations is a product of cyclic groups. We can then consider the group of spatial symmetries of such a periodic system even if we allow translation symmetry to be broken. We will denote such a group by G which denotes the set of spatial symmetries on a torus without any translations. Finally, there is the symmetry group of the open system which is given by F, which can be non-trivial if the boundary is chosen symmetrically.

For example, in the SSH chain, the group of symmetries of the infinite system is generated by a reflection around a given point and translation by one unit cell. This includes mirror reflections around all the points n and n +1 2 / for any integer n . For a finite periodic system with length N , the symmetry group is obtained from the one for the periodic system by the condition T N x = 1 . For a system with periodic boundary conditions but not translations, the symmetry group is generated by two mirrors about the origin x = 0 and the point x = N/ 2 (which sends x to N -x which is equivalent under periodic boundary conditions to -x ). Finally, for the open system, the symmetry group is generated by a single mirror relative to a single fixed origin O which is determined by the choice of boundary termination.

In order to define filling anomalies for BOTP systems, we need a way to go back and forth between the open system and the periodic system. For this reason, we start by considering a special class of boundaries which can be folded back to a torus. We will call such boundaries foldable and define them as terminations which are obtained by cutting a torus into a d -dimensional region topologically equivalent to the disk such the the following conditions are satisfied: (i) the open system is invariant under the point group symmetry F relative to O , (ii) the d -dimensional torus can be obtained by identifying opposite edges/faces of the open region, and (iii) it includes an integer number of unit cells of the translationally symmetric system. This definition means that the open system has the shape of a (large) symmetric unit cell for the periodic system. This special class of boundaries enables us to compare the filling in the open and periodic systems and make the idea of 'resolving" the filling anomaly by identifying boundaries more precise. Once we have established the existence of a filling anomaly in the open system, we will show afterwards how we can define it for more general boundaries, hence establishing in the process that it captures a stable topological invariant of the open system which does not change under any boundary or bulk deformation which preserves the gap at high symmetry surfaces (as well as the bulk gap).

We now define a filling anomaly as follows. Given a periodic Hamiltonian with filling of ν per unit cell, we can define a corresponding open Hamiltonian (or projector) on a foldable boundary whose filling is νN with N denoting the number of unit cells. A Hamiltonian (or a projector) is characterized by a filling anomaly at νN if a gapped, symmetric, Wannierizable and charge-neutral ground state is possible in the periodic system, but impossible on the foldable open boundary. For a system with a filling anomaly we can find a gapped, symmetric and Wannierizable ground state in the open system if we remove the condition of charge-neutrality. Hence, the filling anomaly can be quantified by an integer ν fa denoting the number of electrons we need to add to obtain a gapped, symmetric and Wannierizable ground state in the open foldable system.

Now let us try to generalize this invariant to more generic types of boundaries. Our definition for boundary obstructions identifies Hamiltonians in an open system which are related by symmetric deformations which do not close the gap in the bulk or at any HSS on the boundary. This does allow for gap-closing transitions away from HSSs at the boundary which can change the total filling by pushing some eigenstates across the chemical potential. Such gap-closings have to occur at | F | symmetry related points/surfaces at the boundary which means that the total change in the filling induced by such a process has to be a multiple of | F | . As a result, the integer ν fa can change by a multiple of | F | but is invariant modulo | F | . This defines a topological invariant protected by the gap at HSS which is additive modulo | F | . As a result, the topological phases distinguished by filling anomaly invariants form a subgroup of the integers Z | F | . We can now lift the restriction of foldable boundaries by considering more general symmetric boundary terminations which can be deformed to a foldable boundary without closing the gap at a HSS. Since ν fa mod | F | is protected by the gap at HSSs, it will remain well-defined for such boundaries and can be used to define filling anomalies in a more general context.

To understand how a filling anomaly arises when we open the boundary, we notice that the little group of a point q , i.e. the subset of symmetries which leave q invariant, is generally larger on the torus compared to the open system. More specifically, the little group of q is defined in the open system and the torus as:

$$F _ { q } & = \{ g \in F, g q = q \}, \\ G _ { q } & = \{ g \in G, g q = q \}. \\. \quad. \quad. \quad. \quad. \quad. \quad. \quad. \quad. \quad. \quad. \quad. \quad.$$

It is clear from the definition that F q ⊆ G q since F ⊆ G . As a result, the size of the orbit of q under the action of the point group F is | F / F | | q | ( | F / G | | q | ) in the open (periodic) system. For example, in the DMQI a general point ( x, y ) whose little group contains only the identity element has an orbit with 4 elements in total given by ( x, y ) , ( -x, y ) , ( x, -y ) and ( -x, -y ) under the action of M x and M y . On the other hand, a point lying on one of the mirror lines, let's say ( x, 0) has an orbit with only two points: ( x, 0) and ( -x, 0) . When we open the boundary, a point q in the periodic system maps to n q points in the open system where n q is an integer larger than or equal one given by

$$n _ { q } = | G _ { q } | / | F _ { q } |.$$

Thus, a filling anomaly occurs if (and only if) there is a point

q whose filling in the periodic system is not divisible by n q . The anomaly invariant ν fa can be expressed in terms of n q , the size of the full group | G , | and the filling of q in the periodic system ν pbc q as

$$\nu ^ { \text{fa} } = \frac { \nu _ { q } ^ { \text{pbc} } | F | } { n _ { q } } \mod | F | = \frac { \nu _ { q } ^ { \text{pbc} } | F | | F _ { q } | } { | G _ { q } | } \mod | F |, \\ \text{which is a above on instance in a} \, E | \text{ is a divisible} \, \text{$k$} \, | C | + \, \text{$2D$ no}$$

which is always an integer since | F | is divisible by | G q | (which follows from G q ⊆ F ). Now we can state in particular that a filling anomaly describes a BOTP if there are two gapped bulk Hamiltonians which are symmetrically deformable in the periodic systems, but which correspond to distinct anomaly invariants ν fa in the open system.

## B. Filling anomalies in 2D systems

To investigate possible filling anomalies for BOTPs in general 2D systems, we consider the possible point group symmetries for which the anomaly condition derived above is satisfied. Crystallographic point groups in 2D belong to 2 categories: n -fold rotation groups C n , n = 1 2 3 4 6 , , , , or the dihedral groups D n generated by n -fold rotation and an inplane mirror reflection with n = 1 2 3 4 6 , , , , . We denote the filling of a point q in the open and periodic system by ν q and ν pbc q , respectively. We now note that a point q whose little group G q = C n will be characterized by n q = n/ G | q | . As a result, a filling anomaly is only possible if ν pbc q is not divisible by n . On the other hand, the value of ν pbc q modulo n cannot be symmetrically changed in the periodic system since the only possible symmetric operations involve the addition or removal of n electrons to q which will change its filling by multiples of n . As a result, all filling anomalies for points q with point group C n are associated with bulk phases not BOTPs.

glyph[negationslash]

Next, consider points q whose point group is G q = D n for which | G q | = 2 n such that a filling anomaly is realized whenever ν pbc q is not divisible by 2 n . In contrast to the C n case, there are symmetric deformations which can change the value of ν pbc q mod 2 n . Namely, by symmetrically moving n electrons in or out of q along the mirror invariant lines. As a result, a filling anomaly is associated with a bulk phase whenever ν pbc q = n mod 2 n, and it is associated with a BOTP when ν pbc q = n mod 2 n . For BOTPs this is equivalent to the conditions

$$\nu _ { q } ^ { \text{pbc} } = n, \quad | G _ { q } | = 2 n, \quad | F _ { q } | = 1, \quad ( 2 8 ) \quad \text{erate}$$

which gives a necessary (but not sufficient) condition for a BOTP filling anomaly. One possibility of realising (28) is at the corner of a rectangular sample which does not have any symmetry in the open system ( | F q | = 1 ) but lies at the intersection of two mirrors in the periodic system such that | G q | = 4 . Thus, the filling ν pbc q = 2 in the periodic system leads to a filling anomaly not associated with a bulk phase. This is what happens in the DMQI model. Other cases involve a point q lying at an edge in the open system with | F q | = 1 which becomes a mirror invariant line in the periodic system with G q = 2 . If the periodic filling ν pbc q = 1 ,

FIG. 11. (a) schematic illustration of a BOTP filling anomaly in the 2D point group D 4 associated with edge charges on a square geometry. (b) and (c) more complicated edge terminations which allows the localization of the zero modes so that they can only be removed by closing the gap at a high-symmetry edge.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00011v2-figure-764ddbd1-68bc-4e14-982a-a9cd07a5588a.jpg)

and this position can be symmetrically deformed to the origin, then it satisfies all the conditions of a boundary filling anomaly. An example is illustrated in Fig. 11 with the group D 4 with the points q lying at the edges of the square (with | G q | = 2 ) with ν pbc q = 1 . Such filling anomalies can be observed on generic D 4 symmetric boundaries likes the ones shown in Fig. 11b,c. However, on a simple square boundary with approximate translation symmetry along the edges, the localized charges (or alternatively the zero modes if we assume particle-hole or chiral symmetry) cannot exist at any arbitrary point along the edge and will be pushed towards the corners where they annihilate. Hence, for simple boundary terminations, e.g., a square or hexagon, we are not able to identify any interesting physical signatures of BOTPs other than what was already known for the DMQI. For the rest of this work, we will instead focus on BOTPs in 3D.

## V. 3D MODELS WITH BOUNDARY OBSTRUCTIONS

## A. Recipe for constructing 3D BOTPs

glyph[negationslash]

We now switch our attention to BOTPs in 3D. In the following we present a simple recipe that can be used to construct 3D models with different types of boundary obstructions. The recipe works in analogy to the construction of the 2D DMQI model. For that model, we took a pair of inversion symmetric SSH chains in each unit cell parallel to say the x -direction, and then coupled them with dimerized couplings, i.e., intra- and inter-cell couplings γ y and λ y along the y -direction (i.e., the stacking direction). Additionally, to generate the necessary symmetry structure, we threaded a π -flux per unit cell (Fig. 1). In the limit γ y = 0 , λ y = 0 , there is an isolated/unpaired SSH chain on each of the edges parallel to the x -axis that can be in the trivial ( | γ x | &gt; λ | x | ) or obstructed ( | γ x | &lt; λ | x | ) atomic limit, and thus may have an associated boundary obstruction.

This construction suggests a route to obtaining 3D BOTPs: we start with any 2D Hamiltonian H 2D ( k ) associated with boundary signatures such as corner/edge charges or 1D edge states. We then consider a pair of these Hamiltonians such that their sum is topologically trivial, and hence all the associated boundary signatures can be removed when they are coupled. Next, we couple these pairs via dimerized SSH-

glyph[negationslash]

glyph[negationslash]

like couplings γ , λ z z along the (stacking) z -direction. This arrangement guarantees that in the limit γ z = 0 , λ z = 0 , there is a single isolated copy of ±H 2D ( k ) at the upper/lower surface when the system is considered with open boundaries. Additionally, we need to to ensure that the bulk gap is open if either H 2D ( k ) is gapped or | γ z | = | λ z | . This can be done by choosing the 2D Hamiltonian to have alternating sign between layers. In most cases, this is realized by inserting π fluxes between the stacked layers (i.e., in the xz and yz plaquette types) possibly with the addition of some other modifications to the Hamiltonian parameters between layers as we will show later. The resulting 3D Hamiltonian is given by 1a 1f 1g 1c 1b 1d Bulk (a) 1 a 1 c 1 1 z unit cell (b) νz =0 νz =0 νz =0 νz =½ νz =½ νz =½ νx,y =(0,½) νx,y =(½,0) νx,y =(0,½)

1

surface obstructed

f

g

1

b

d

x

$$\mathcal { H } _ { 3 D } ( k ) = \mathcal { H } _ { 2 D } ( k ) \tau _ { 3 } + ( \lambda _ { z } \sin k _ { z } ) \tau _ { 2 } + ( \gamma _ { z } + \lambda _ { z } \cos k _ { z } ) \tau _ { 1 },$$

where τ 1 2 3 , , denote the Pauli matrices in the layer subspace along z . Due to the π fluxes, the different terms of the Hamiltonian anticommute which implies that it can only be gapless if each of the three terms separately vanish, i.e., only if H 2D ( k ) is gapless and | λ z | = | γ z | . As in the DMQI case, the existence of π fluxes means that spatial symmetries are generally combined with gauge transformations and transform projectively. In particular, since the π -fluxes here are only considered in the vertical xz and yz planes but not the horizontal xy plane, the resulting projective representation cannot be recast as a spinful representation as we will see later.

It is instructive to highlight the relationship between our construction and the layer construction used to build higherorder and topological crystalline phases from lower dimensional topological phases [52-54, 68]. In the conventional layer construction, one starts from a lower dimensional topological phase placed at a symmetry-invariant region and then symmetrically adds layers or copies of this lower dimensional system to realize a higher-dimensional system. A system constructed this way will always be topologically non-trivial in the bulk as long as the protecting spatial symmetry is unbroken. In contrast, we start with a lower dimensional system H 2D which is topologically trivial but consists of two components which are individually non-trivial. Upon repeating this lower-dimensional unit using the layering procedure, we can spatially separate these two non-trivial components such that they appear at opposite edges/surface of the system. The resulting system is topologically trivial in the bulk with lower dimensional non-trivial systems living at its boundary.

The type of boundary obstruction in the 3D Hamiltonian (29) depends on the type of obstruction for H 2D ( k ) , as well as the dimension of its surface states. For the top/bottom surfaces, H 3D ( k ) will have a surface (hinge) obstruction if H 2D ( k ) has a bulk (edge) obstruction. For side surfaces, H 3D ( k ) will have a surface (hinge) obstruction if H 2D ( k ) has edge states/charges (corner charges). The latter statement can be understood by investigating the dimension of the SPT needed to add to the side surfaces to cancel the surface states/charges (remember all BOTPs are also extrinsic higherorder topological insulators whose surface states can be removed by adding a lower dimensional SPT on the surface). If the 2D Hamiltonian has corner charges, then it is enough to

1

h

1

e

FIG. 12. Generic phase diagram for H 3D ( k ) obtained by the three dimensional construction in Sec. VA, where the two dimensional Hamiltonian H 2D ( k ) is parameterized by { ξ 1 , ξ 2 } . As an illustration we choose the 2D model to exhibit two topologically distinct phases depending on which of ξ 1 or ξ 2 is the largest. Hence, H 2D ( k ) has a bulk gap closing transition at ξ 1 = ξ 2 that will subsequently lead to a Wannier gap closing in H 3D ( k ) in the W z Wannier spectrum (blue plane), which can happen at ν z = 0 or ν z = 1 2 / as labeled in the figure. On the other hand, the gap in the W x and W y Wannier spectra closes when the stacks of H 2D ( k ) are uniformly coupled in the stacking direction, i.e., when λ z = γ z (yellow plane), again either at ν x,y = 0 or ν x,y = 1 2 / as labeled in the figure. The intersection between the two planes corresponds to a line where the bulk gap is closed in the phase diagram. This construction yields four phases A, B, C and D separated by Wannier gap closings. A choice of boundary may select one of the four to be topologically distinct from the other three (in the sense of a boundary obstruction) when the Wannier gap closing happens at the Wannier chemical potential. Here we highlighted phase D as a non-trivial boundary obstructed topological phase for a choice of boundary termination that coincides with the unit cell structure, µ i = 1 2 / .

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00011v2-figure-96e2edf8-7a61-4a08-a9b1-cb7c92efc330.jpg)

add 1D SSH chains to the vertical hinges to cancel them. On the other hand, a 2D Hamiltonian with edge states or charges requires the addition of a 2D SPT on the side surfaces to cancel these states/charges. It follows from the previous discussion that H 3D ( k ) has a surface (hinge) obstruction if and only if H 2D ( k ) has a bulk (edge) obstruction with 1D (0D) edge states/charges. We note that it is also possible to have mixed surface-hinge obstructions where connecting two BOTPs involves either a surface gap-closing on the upper/lower surfaces or a hinge gap-closing on the side surfaces. This would be the case if, for example, we take H 2D ( k ) to be a bulkobstructed atomic insulator with corner charge, e.g., the C 4 z quadrupole model. It is worth noting that the distinction between surface and hinge types of obstructions is only possible for boundaries consisting of several intersecting 2D planes. On a more smooth boundary like a sphere, such a distinction is ill-defined.

Nowlet us more explicitly describe the properties of our 3D parent model. A representative phase diagram of the above construction can be found in Fig. 12. There we consider a two dimensional Hamiltonian parametrized by generic parameters

y

glyph[negationslash]

ξ 1 and ξ 2 . For illustration we assume the 2D Hamiltonian admits a topologically nontrivial phase with boundary modes if ξ 1 &gt; ξ 2 , a trivial phase when ξ 1 &lt; ξ 2 and a critical point when ξ 1 = ξ 2 . When stacking into a three dimensional system, the bulk 3D Hamiltonian will be gapped even when ξ 1 = ξ 2 as long as there is a dimerization along z , λ z = γ z . As a function of ξ 1 , ξ 2 , and γ /λ z z there is a line in the phase diagram (red line) where the bulk is gapless. By construction, this model allows for four adiabatically connected bulk phases, separated by Wannier gap transitions where the Wannier spectra changes its topology. A choice of a symmetric boundary termination will single out one of the phases, where the Wannier gap closing happens at the WCP, (for our choice this is at ν x,y = 1 2 / and ν z = 1 2) / (c.f. Fig. 12). For this boundary termination, which coincides with the unit cell structure, phase D is a nontrivial boundary obstructed topological phase since λ z &gt; γ z leaving a dangling, unpaired H 2D ( k ) on both the top and bottom surfaces that is in its topologically nontrivial phase.

In the following two subsections, we will focus on models with surface obstructions while briefly discussing an example of a hinge obstruction. We will separately consider the cases where the surface obstruction is associated with a filling anomaly or gapless surface states.

## B. 3D BOTPs with filling anomalies

In this subsection, we will discuss 3D BOTPs associated with filling anomalies. We will propose a class of models with C 2 nh symmetry n = 1 2 3 , , exhibiting surface obstructions associated with a filling anomaly. After introducing the models, we analyze their properties using a real-space approach as well as through the Wannier spectrum. At the end, we will briefly discuss an example of a hinge-obstructed model with corner charge.

## 1. A new class of C 2 nh Hamiltonians

Now, let us introduce a new class of 3D models defined for the space groups 10, 83, and 175 which are symmorphic space groups whose point group F is C 2 nh with n = 1 2 3 , , , respectively. These symmetry groups are characterized by a 2 n -fold rotation C 2 nz accompanied by a mirror reflection perpendicular to the rotation axis M z .

We build the 3D Hamiltonians in each case following the recipe of Sec. VA using 2D models that realize a C 2 nz insulator at filling n. The full 3D model will be a 3D BOTP at filling of 2 n which has a surface obstruction associated with a filling anomaly. The 2D tightbinding Hamiltonians we use are shown in Fig. 13. Each of these 2D models has several distinct atomic limit phases separated by a bulk transition: one atomic limit where all the 2 n electrons are at the center of the unit cell (with site symmetry group C 2 nz ) and others where the 2 n electrons are distributed in pairs to n symmetry related edges of the unit cell (with site symmetry C 2 z ). The non-trivial obstructed atomic limit of these models is characterized by edge charges and bulk polarization so they can be understood as

FIG. 13. (a,b,c) Schematic illustration of the 2D hopping Hamiltonians used to build the 3D surface obstructed models with (a) C 2 z , (b) C 4 z , and (c) C 6 z symmetries. The hopping parameters λ 0 , λ x and λ y are denoted by red, green, and blue lines respectively. (d,e,f) Top view illustration of the deformation process explained in the text between position 1 a and position nc for (d) C 2 z , (e) C 4 z , and (f) C 6 z symmetries.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00011v2-figure-76cb5b32-3003-4207-b68f-9ff2712a43cf.jpg)

weak 2D phases if translation symmetry in the 2D plane is preserved.

Explicitly, the C 2 z model is given by

$$\text{for} \quad \mathcal { H } _ { 2 D } ^ { C _ { 2 x } } ( k ) = ( \lambda _ { x } & \sin k _ { x } + \lambda _ { y } \sin k _ { y } ) \sigma _ { 2 } \\ \text{ace} \quad \ \ \ + ( \lambda _ { 0 } + \lambda _ { x } \cos k _ { x } + \lambda _ { y } \cos k _ { y } ) \sigma _ { 1 }, \ \ ( 3 0 )$$

where σ 1 2 3 , , denote the Pauli matrices distinguishing the two atomic orbitals inside the unit cell. The hopping parameters λ 0 ,x,y are represented by the black, green, and blue links in Fig. 13a, respectively. The spectrum of Eq. (30) is given by

$$\overset { \i n z \text{ in} } { \text{TP at} } \quad \varepsilon _ { 2 D } ^ { 2 } = \lambda ^ { 2 } + 2 \lambda _ { 0 } ( \lambda _ { x } \cos k _ { x } + \lambda _ { y } \cos k _ { y } ) + 2 \lambda _ { x } \lambda _ { y } \cos ( k _ { x } - k _ { y } ), \\ \text{with a} \quad \dots \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \colon$$

with the shorthand λ 2 = λ 2 0 + λ 2 x + λ 2 y . The 2D model is gapped whenever the largest of λ 0 ,x,y exceeds the sum of the other two. In addition, the gapped phases at half filling correspond to distinct atomic insulators with a single electron localized either at the 1 a , 1 b, or 1 c positions whenever the largest λ is λ 0 , λ , x or λ y , respectively.

The 2D Hamiltonians for C 4 z and C 6 z symmetric cases can be constructed in a straightforward way. We simply stack two or three copies of H C 2 z 2D ( k ) related by C 4 z or C 6 z rotations,

respectively. They are shown schematically in Fig. 13b,c. Let us here explicitly show the resulting C 4 z model, by adding to the C 2 z model a π/ 2 rotated copy of itself,

$$\mathcal { H } _ { 2 D } ^ { C _ { 4 z } } ( k ) & = \left ( \begin{array} { c c } \mathcal { H } _ { 2 D } ^ { C _ { 2 z } } ( k ) & 0 \\ 0 & e ^ { i \frac { \pi } { 4 } \sigma _ { 1 } } \mathcal { H } _ { 2 D } ^ { C _ { 2 z } } ( C _ { 4 z } k ) e ^ { - i \frac { \pi } { 4 } \sigma _ { 1 } } \end{array} \right ), \quad \begin{array} { c } \Xi \\ \Xi \\ \end{array} \right ) \\ \text{share} \, \mathcal { H } _ { 2 } \, \mathcal { H } _ { 2 } \, = \, / \, \mathcal { H } _ { 2 } \, \mathcal { H } _ { 2 } \, \mathcal { H } _ { 2 } \, \mathcal { H } _ { 2 } \, \mathcal { H } _ { 2 } \, \mathcal { H } _ { 2 } \, \mathcal { H } _ { 2 } \, \mathcal { H } _ { 2 } \, \mathcal { H } _ { 2 } \, \mathcal { H } _ { 2 } \, \mathcal { H } _ { 2 } \, \mathcal { H } _ { 2 } \, \mathcal { H } _ { 2 } \, \mathcal { H } _ { 2 } \, \mathcal { H } _ { 2 } \, \mathcal { H } _ { 2 } \, \mathcal { H } _ { 2 } \, \mathcal { H H } _ { 2 } \, \mathcal { H H } _ { 2 } \, \mathcal { H H } _ { 2 } \, \mathcal { H H } _ { 2 } \, \mathcal { H H } _ { 2 } \, \mathcal { H H } _ { 2 } \, \mathcal { H H } _ { 2 } \, \mathcal { H H } _ { 2 } \, \mathcal { H H } _ { 2 } \, \mathcal { H H } _ { 2 } \, \mathcal { H H } _ { 2 } \, \mathcal { H H } _ { 2 } \, \mathcal { H H } _ { 2 } \, \mathcal { H H } _ { 2 } \, \mathcal { H H } _ { 2 } \, \mathcal { H H } _ { 1 } \, \mathcal { H H } _ { 2 } \, \mathcal { H H } _ { 2 } \, \mathcal { H H } _ { 2 } \, \mathcal { H H } _ { 2 } \, \mathcal { H H } _ { 2 } \, \mathcal { H H } _ { 2 } \, \mathcal { H H } _ { 2 } \, \mathcal { H H } _ { 2 } \, \mathcal { H H } _ { 2 } \, \mathcal { H H } _ { 2 } \, \mathcal { H H } _ { 2 } \, \mathcal { H H } _ { 2 } \, \mathcal { H H } _ { 2 } \, \mathcal { H H } _ { 2 } \, \mathcal { H H } _ { 2 } \, \cdots$$

where C 4 z k = ( k , y -k x ) is the momentum space action of C 4 z . This Hamiltonian describes a 2D four-band model.

The 3D models can be constructed following the recipe of Sec. V A stacking pairs of 2D models with dimerized couplings γ z and λ z along the z-direction. Additionally, we thread a π -flux per unit cell in the xz and yz plaquettes but not in the xy plaquette. This flux threading modifies the symmetry representation for the spatial symmetries to a projective representation by combining spatial symmetries with gauge transformations such that the gauge invariant condition

$$C _ { 2 z } M _ { z } C _ { 2 z } ^ { - 1 } M _ { z } ^ { - 1 } = - 1 \quad \quad \ ( 3 3 ) \quad \tilde { \text{clez} } \\ \quad \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \$$

is satisfied. This condition captures the π flux in the vertical ( xz and yz ) plaquettes. Notice that this condition implies that we cannot write this projective representation as a nonprojective spinless or spinful representation: the former has M 2 z = C 2 2 z = +1 whereas the latter has M 2 z = C 2 2 z = -1 with both having [ M ,C z 2 z ] = 0 [69]. As a result, the inversion operator which is given by I = C 2 z M z squares to +1 in any non-projective spinful or spinless representation whereas it squares to -1 when (33) is satisfied meaning that it cannot be consistent with any non-projective spinless or spinful representation. This arises from the special property of the models we consider where π -fluxes are attached only to the vertical plaquettes but not the horizontal ones.

For the 3D Hamiltonian defined by (29) with H 2D given in (30), the spatial symmetries C 2 z and M z are implemented as

$$M _ { z } \mathcal { H } ( k ) M _ { z } ^ { \dagger } & = \mathcal { H } ( m _ { z } k ), \quad \ \ M _ { z } = \sigma _ { 3 } \tau _ { 1 } \quad \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \
 \, \quad C _ { 2 z } \mathcal { H } ( k ) C _ { 2 z } ^ { \dagger } & = \mathcal { H } ( C _ { 2 z } k ), \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$$

where m z k = ( k , k x y , -k z ) and C 2 z k = ( -k , x -k , k y z ) . The eight-band C 4 h model has the symmetries

$$C _ { 4 z } & = e ^ { - i \frac { \pi } { 4 } } \begin{pmatrix} 0 & e ^ { i \frac { \pi } { 4 } \sigma _ { 1 } } \\ e ^ { i \frac { \pi } { 4 } \sigma _ { 1 } } & 0 \end{pmatrix} \tau _ { 0 }, \quad M _ { z } = \begin{pmatrix} \sigma _ { 3 } & 0 \\ 0 & \sigma _ { 3 } \end{pmatrix} \tau _ { 1 }. \quad \text{vertic} \\ \intertext { u s s } \sim \quad \text{and is the same as $e$ with the $u$ with the $u$} \, \infty \, \infty$$

Here σ 1 2 3 , , act in the space of orbitals within the same 2D layer, whereas τ 1 2 3 , , act in the space of orbitals in different layers. The C 6 z model has a similar structure.

## 2. Real space picture and filling anomaly

We now explain the boundary obstruction in the models described in the previous subsection. We can start with a similar picture to that of Sec. II B. In this picture, we consider two configurations of Wyckoff positions which are smoothly deformable in the bulk, but not in the presence of a surface. The precise meaning of this procedure follows from the machin-

FIG. 14. (a,b,c) local density of states at zero energy for the surface obstructed insulators constructed by stacking the tight binding models from Fig. 13 with (a) C 2 z , (b) C 4 z , or (c) C 6 z symmetry. We can clearly notice the existence of hinge-localized states associated with hinge charges if we require the system to be symmetric and gapped, or alternatively a partially filled hinge band if we require the system to be symmetric and neutral; both phenomena are a result of the filling anomaly. (d,e,f): illustration of the 0D surface states associated with the filling anomaly for a spherical surface termination. Such states can only be annihilated at the high symmetry regions at the equator or north/south poles leading to a gap-closing at these points.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00011v2-figure-eb8b92a2-25de-4c8a-a6cd-69d32ff979a9.jpg)

ery of Sec. III A wherein the deformation trajectory in the periodic system induces a corresponding trajectory in the open system via the boundary map σ . By tracing the movement of the charge centers, we can then show that such a trajectory cannot be both symmetric and smooth. We consider the real-space deformation between two distinct configurations of the electrons: one where all electrons are at the center of the unit cell 1 a (c.f. Fig. 15, red Q ) in the z = 0 plane, and the other with two electrons at the edge of the unit cell 1 c , 2 c or 3 c ( nc for shorthand) in the three different models (c.f. Fig. 15, purple Q ) in the z = 1 2 / plane. This can be done while preserving the symmetry by first moving the 2 n electrons at 1 a horizontally in the z = 0 plane into nc as shown in Fig. 13; and then moving the 2 electrons at this site vertically in opposite directions. We could also start by first moving n electrons vertically in each direction, then moving them horizontally to the position nc .

Now we can see that the first trajectory (which involves moving horizontally then vertically) will be obstructed in the presence of side surfaces at the unit cell edge. The reason is that at such side surfaces, the filling of the nc position is 1, and thus we cannot further move the electrons vertically while preserving the mirror symmetry in the z direction. Similarly, the second trajectory involving vertical then horizontal movement is not possible in the presence of the top/bottom surfaces at the z = 1 2 / plane. The reason is that the filling of the position 1 a at this surface is n (rather than 2 n which is its bulk filling). This does not allow us to symmetrically move the electrons in the plane to position nc without breaking the symmetry. This discussion parallels the one of Sec. II B for the DMQI, and it shows that in the presence of a surface ter-

mination consistent with the unit cell shown in Fig. 13, it is not possible to connect the atomic configuration where the 2 n electrons are at 1 a in the z = 0 plane to the one where 2 electrons are at each position nc at the z = 1 2 / plane, despite the fact that these two atomic configurations are symmetrically deformable to each other in the bulk (with periodic boundaries). We also note that, similar to the DMQI case, the possibility of interpolating the two atomic configurations in the bulk imposes some constraints on the commutation properties of the different symmetries, i.e.,

$$\{ C _ { 2 z }, M _ { z } \} = 0., \quad \quad \ \ ( 3 6 ) \quad \stackrel { \varkappa } { C _ { 2 } } \\ \text{$\text{$\text{$\text{$\text{$ioc}$} } }$$

which is shown to be a necessary condition for the smooth deformation between these different bulk trajectories in Appendix B.

Similar to the DMQI, the 3D models introduced in this section are chatacterized by filling anomalies which are not associated with bulk topology. They can be understood by employing a similar argument to that of Sec. II D. Let us start with the n = 1 model corresponding to 3D space group 10 whose point group is C 2 h . Taken on an open symmetric boundary, let us consider the filling of the low-dimensional symmetry-invariant region Γ given by the union of the xy plane and the z -axis. This region consists of the disjoint union of three regions: (i) the origin O , (ii) the z -axis excluding the origin (described by the line x = y = 0 ), and (iii) the xy -plane excluding the origin. The filling of (ii) and (iii) is necessarily even due to the action of M z and C 2 z , respectively. The filling of (i) is even due to the anticommutation of M z and C 2 z which enforces the (projective) symmetry representation to be two-dimensional. This two dimensional representation is neither spinless or spinfull, and we will refer to it in Table I as ¯ E is analogy to the DMQI model, however this is slight abuse of notation since this projective representation is not tabulated or equivalent to a purely spinless or spinfull representation [70]. As a result, the total filling of Γ is even. On the other hand, adding or removing electrons from outside this region can only change ν Γ by 4 corresponding to the four symmetry-related locations ( x, y, z ) , ( -x, -y, z ) , ( x, y, -z ) , and ( -x, -y, -z ) (that is to say that the general Wyckoff position in the open system has multiplicity 4). Thus, the parity of ν Γ / 2 describes a Z 2 invariant which cannot be changed without changing the filling of the high-symmetry region Γ . This is only possible by closing the gap at a high symmetry surface (or in the bulk). As in the DMQI case, the discrepancy between this number and the total filling indicates a filling anomaly described by the Z 2 invariant.

Such an anomaly can be resolved by identifying opposite side faces of the sample. This introduces new C 2 z -invariant two-fold rotation axes whose filling is even, but can take a value of 0 or 2 modulo 4, which makes it possible for it to compensate for the discrepancy between the total filling and ν Γ modulo 4. The orbitals in these new C 2 z -invariant lines can be symmetrically brought to Γ without closing the bulk gap or breaking the symmetry, thereby changing the value of ν Γ and removing the bulk signature of the filling anomaly observed in the open system. This indicates that such filling anomaly is associated with a BOTP rather than a bulk phase. Similar analysis can be carried for the models with point groups C 4 h and C 6 h .

Another approach to filling anomalies is the one introduced in Sec. IV. In this case, we choose a foldable boundary termination as shown in Fig. 14a,b,c which resembles a symmetric large cell for the periodic system. Recall the criterion derived in Sec. IV where a filling anomaly is realized whenever the filling of a position q on the periodic system ν pbc q is not divisible by n q , defined as the ratio of the order of its site symmetry group in the periodic and the open system (Eq. 27). For all the C 2 nh models discussed above, the bulk filling for the nc position at z = 1 2 / in the obstructed 'phase" is ν pbc q = 2 , and its symmetry group on the torus is G q = C 2 h (generated by C 2 z and M z ) which has four elements, | G q | = 4 . For the surface termination where the vertical surfaces are at the edges of the unit cell (shown in Fig. 13a,b,c), and the horizontal surfaces are at z = 1 2 / , the sites q lying at the horizontal hinges (see Fig. 14a,b,c) have no symmetry left. Thus | F q | = 1 for these points which implies that the condition of Eq. 28 is satisfied.

Finally, let us discuss the observable consequences of the filling anomaly. In the DMQI, the filling anomaly manifests itself in the existence of fractional corner charge. Here, we similarly expect fractional charges placed at 4 n symmetry related points at the corners of a generic boundary termination. These corner charges can be more easily undersood by first imposing chiral symmetry or particle-hole symmetry and noting that we can associate the anomaly with 2 n 0D surface states at generic symmetry related points on an arbitrary symmetric surface as shown in Fig. 14(d,e,f). Such 0D states can only be annihilated at the equator or the north/south poles which are HSSs. Following the discussion of Sec. II D, we can show that the filling anomaly will survive even after breaking chiral symmetry. In this case, it will appear as excess charge fractionally distributed among these symmetry related points for the symmetric gapped insulator. We notice that on highly symmetric foldable terminations such as the ones shown in Fig. 14, the edges have an approximate translation symmetry. As a result, the filling anomaly for a symmetric insulator implies the existence of fractional charges per unit cell along the hinge for a symmetric gapped insulator, i.e., a hinge charge. Alternatively, we can impose the condition of charge neutrality (and unbroken symmetry) to deduce that the filling anomaly gives rise to partially filled bands corresponding to gapless hinges in these models.

## 3. Wannier spectrum

We now turn our attention to the evolution of the Wannier bands and how they encode the boundary obstructions described in the previous subsection. We diagnose the boundary obstruction by analyzing the Wannier spectrum and choosing a lattice termination that coincides with the edges of the unit cell. That is, the WCP is given by µ = (1 / 2 1 , / 2 1 , / 2) , which means that only gap closings of the Wannier spectrum at these values correspond to actual boundary gap closings. A phase diagram of the models is shown for reference in Fig. 15.

glyph[negationslash]

C 2 h model -Agapped Wannier spectrum along ˆ z is possible provided H C 2 z 2D is in a gapped phase. This only happens when 2 λ M &gt; λ S (blue regions in Fig. 15), where we defined the parameters λ S = λ x + λ y + λ 0 and λ M = max( λ , λ x y , λ 0 ) . On the other hand, the Wannier spectrum along ˆ x or ˆ y is gapped as long as | γ z | = | λ z | (yellow plane in Fig. 15). The bulk Hamiltonian is gapped if either condition is satisfied, and is therefore only gapless in the intersections of the blue region and yellow plane in the phase diagram. The region of intersection is marked in red. It is evident that the Wannier gap may close while preserving the bulk gap, a necessary condition for a BOTP. We show the Wannier spectra ν z ( k , k x y ) and ν x ( k , k z y ) in Figs. 16a and 17a. In these figures, we can see six distinct regions with gapped Wannier spectra separated by Wannier gap closings, as will be explained in detail below.

Similar to the DMQI case, the different gapped phases of the C 2 h model are characterized by distinct Wannier band representations, which can be reflected in distinct values of the nested Wannier band polarizations. Let us first notice that the model has four quantized Wannier band polarizations: p z,x and p z,y are quantized due to C 2 z symmetry, and p x,z and p y,z are quantized by M z symmetry. In our model, p z,x and p z,y cannot be simultanuously non-zero, and p x,z and p y,z are equal, leading to six distinct gapped phases (see Fig. 15). Let us label the different phases by the maximal Wyckoff position Q around which the Wannier centers are located, as shown in Fig. 15b. While the electrons may be symmetrically displaced from this maximal Wyckoff position in some cases, they remain centered around Q , and both their band representations, and Wannier band representations, are topologically equivalent to those of the phase where the electrons lie exactly at Q . Therefore it is useful to label the phases by such maximal Wyckoff positions. We can use the notation ( p z,x , p z,y , p x,z , p y,z ) Q , to label the six distinct phases, which are given by: (0 0 0 0) , , , 1 a , (0 , 1 / 2 0 0) , , 1 b , (1 / 2 0 0 0) , , , 1 d , (0 , 0 1 , / 2 1 , / 2) 1 c , (0 1 , / 2 1 , / 2 1 , / 2) 1 g , (1 / 2 0 1 , , / 2 1 , / 2) 1 f . In our model we cannot tune parameters such that the Wannier centers lie at the 1 e or 1 h Wyckoff positions.

The Wannier transitions between the different phases can be studied numerically leading to the phase diagram in Fig. 15a. As we can see, all of the pairwise Wannier transitions between the phases 1 a , 1 b, and 1 d are associated with a gap-closing of the Wannier spectrum in the z -direction at ν z = 0 . On the other hand, the pairwise Wannier transitions between the phases 1 c , 1 g, and 1 f are associated with a gap-closing for the Wannier spectrum in the z -direction but at ν z = 1 2 / . This reflects the fact that transitions between 1 a , 1 b, and 1 d take place in the z = 0 plane, whereas transitions between 1 c , 1 f, and 1 g take place in the z = 1 2 / plane. The Wannier transitions 1 a ↔ 1 c , 1 b ↔ 1 f and 1 g ↔ 1 d are associated with a gap-closing in both Wannier spectra along the x and y directions occuring at ( ν x , ν y ) = (0 0) , , (0 , 1 2) / , and (1 / 2 0) , , respectively.

Following the previous discussion, we can see that the phases 1 f and 1 g have a surface obstruction, which can be diagnosed by p x,z = p y,z = 1 2 / and ( p z,x , p z,y ) = (0 , 1 / 2) or (1 / 2 0) , , respectively. They are characterized

FIG. 15. (a) Phase diagram for the C 2 h model introduced in Sec. VB1 with a fixed value of λ 0 = 1 2 / . The Wannier spectrum W z closes at the blue regions which span the γ /λ z z axis in the areas marked with W z ; while the Wannier spectrum W x closes at the yellow plane at γ /λ z z = 1 . Due to the constraints in the model, a gap closing in W x implies a gap closing in W y . In the intersection between the two, the red regions, the bulk is gapless The gapped phases with gapped Wannier spectra are marked by the representative maximal Q from which the Wannier band representations may be induced. (b) Spatial location of the maximal Q 's compatible with the C 2 h point group. With the conventional choice of Wannier chemical potential µ = (1 / , 2 1 / , 2 1 / 2) , phases induced from 1 f, 1 g , and 1 h may be surface obstructed; from 1 h they could also be hinge obstructed. In our model we realize 1 f and 1 g surface obstructions.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00011v2-figure-229662bb-d2fa-48fe-bd3c-063f551657da.jpg)

by a quantized hinge charge of 1 / 2 per unit cell, at the hinges along y or x respectively. A real space calculation of the low energy states is found in Fig. 14(a,b,c). As an example, in Figs. 16 and 17 we show how the two gapped phases, 1 a and 1 g , which correspond to the parameters ( λ , λ 0 x , λ y , γ z , λ z ) = (0 5 0 0 25 0 5 1) . , , . , . , and ( λ , λ 0 x , λ y , γ z , λ z ) = (0 5 1 0 25 1 5 1) . , , . , . , respectively, are separated by a Wannier gap closing transition at the WCP, and consequently a surface gap closing. We follow two distinct paths, either by first increasing λ x : 0 → 1 , and then increasing γ z : 0 5 . → 1 5 . , or vice versa. When increasing λ x first, the gap at the z = 1 2 / surface closes in the range of 0 25 . ≤ λ x ≤ 0 75 . , (c.f. gapless ν z ( k , k x y ) in Fig. 16 a2 and a3). Instead, if we increase γ z first we find the gap closing at the WCP at γ z = 1 , (c.f. gapless ν x ( k , k y z ) in Fig.17 b). Thus, a surface gap-closing at the z = 1 2 / or x = 1 2 / surface is unavoidable whenever we connect the 1 a and 1 g phases.

glyph[negationslash]

Finally, let us comment on the Wannier band representations. The Wannier spectrum of W z , corresponds to the upper and lower surfaces (which preserve C 2 ), while the Wannier spectrum of W x/y correspond to the side surfaces (which preserve M z ). The WBRs corresponding to w z and w x/y are split into two disconnected WBRs when 2 λ M &gt; λ S or | γ z | = | λ z | , respectively. The values that the WBRs take at high symmetry momenta may be used to determine the topological phases since in this case they can be uniquely associated with a representation π of the open boundary system. The symmetry eigenvalues for the WBRs are shown in Table. I. To determine these symmetry indicators we induce the bulk BR's from

FIG. 16. Evolution of the Wannier spectrum W z with eigenvalues ν z ( k , k x y ) for different values of λ x in the three-dimensional boundary obstructed models with filling anomalies and C 2 h , C 4 h and C 6 h symmetries. The other parameters are kept fixed at the values of λ 0 = 5 . and λ y = 25 . with λ z = 1 5 . and γ z = 1 . With gapped ν x (see text), we find that varying λ x allows for two gapped phases bordered by an extended region where ν z is gapless. This happens when . 25 &lt; λ x &lt; . 75 where the condition 2 max( λ , λ 0 x , λ y ) &lt; λ ( 0 + λ x + λ y ) is satisfied. Panels (a1) and (a4) correspond to the gapped phases 1 f and 1 g in the phase diagram of Fig.15, respectively.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00011v2-figure-801ad037-3a01-428a-be33-6b91c99b7b80.jpg)

TABLE II. Upper section of rows: Diagnosis of boundary obstructions through Wannier band representations in the C 2 h model introduced in Sec. V B 1. The left column indicates the maximal Wyckoff position Q that labels the phase with periodic boundaries. In the columns to the right we show the symmetry representations at high symmetry momenta of one of the two WBRs. They are evaluated for C 2 z and M z according to which symmetry is preserved in the respective direction. Lower section of rows: Other possible WBR configurations compatible with the C 2 h point group, but not realized by our model. Throughout the table, the bold lines correspond to the phases with a boundary obstruction when the conventional Wannier chemical potential is chosen, i.e., µ = (1 / 2 1 , / 2 1 , / 2) and the center of the sample is located at the position 1 a . The obstruction in these cases is of the filling anomaly type, where the multiplicity of the sites 1 g , 1 f and 1 h is larger in the open boundary system compared to the periodic system, see Fig. 15b. We see that for this model the obstructed phases are uniquely indicated by the symmetry eigenvalues of the Wannier bands. Here we use ¯ E in analogy to the DMQI model, however this projective representation is not tabulated in the usual symmetry representation of the three-dimensional point group C 2 h . The definition of this representation is discussed in the main text.

| k             | Γ             | X       | Y     | M     | Γ           | Z       |     |
|---------------|---------------|---------|-------|-------|-------------|---------|-----|
| ρ             | w z + ( C 2 z | ) w z + | w z + | w z + | w x + ( M z | ) w x + | ν Γ |
| (1 a, ¯ E )   | +             | +       | +     | +     | +           | +       | 2   |
| (1 b, ¯ E )   | +             | +       | -     | -     | +           | +       | 2   |
| (1 c, ¯ E )   | +             | +       | +     | +     | +           | -       | 2   |
| (1 d, ¯ E )   | +             | -       | +     | -     | +           | +       | 2   |
| ( 1 f , ¯ E ) | +             | +       | -     | -     | +           | -       | 0   |
| ( 1 g , ¯ E ) | +             | -       | +     | -     | +           | -       | 0   |
| (1 e, ¯ E )   | +             | -       | -     | +     | +           | +       | 2   |
| ( 1 h , ¯ E ) | +             | -       | -     | +     | +           | -       | 0   |

the maximal Q 's and subduce them to the subgroup maintained by the respective Wannier spectrum, as described in

Sec. III B. The subduction of the bulk BRs to the point group F for a fully open system (assuming the sample center to be at O = 1 a ) are also computed to directly associate the WBRs to the obstructed or trivial phases under this choice of boundary. The table shows exhaustively the WBRs compatible with this band representation, therefore we can observe that in this case the obstructed phases are symmetry indicated from the symmetry eigenvalues of the Wannier bands. As anticipated earlier, we can see that the C 2 z eigenvalues for 1 a/ c 1 , 1 b/ 1 f, and 1 d/ g 1 correspond to a 2D inversion symmetric insulator with a single electron localized at (0 0) , , (0 1 , / 2) , and (1 / 2 0) , , respectively. Similarly, The M z eigenvalues correspond to a 2D mirror-symmetric system where the electron is localized at the 0 or 1 / 2 mirror lines for 1 a/ b/ 1 1 d or 1 c/ 1 f/ g 1 , respectively.

C 4 h model -We now look at the 8 band model with C 4 h symmetry. The model differs from the C 2 h model in two main aspects. First, C 4 z symmetry dictates that the Wannier spectra W x and W y are identical, hence the nested polarizations p z,x and p z,y are equal. This means the model has only one nontrivial phase, and this phase has hinge charge along both the x and y directed hinges. In this phase, the occupied electrons are distributed around the 1 f and 1 g positions in Fig. 15b. Second, with four occupied electrons per unit cell, the Wannier bands allow for more complicated gap-closing patterns, as we will see below.

The analysis of the C 4 h model parallels the C 2 h model above. The model has four distinct regions of gapped Wannier spectra characterized by the values of the two distinct nested polarizations ( p z,x , p x,z ) equal to (0 0) , 1 a a 1 , (0 1 , / 2) 1 1 c c , (1 / 2 0) , 1 1 b d , and (1 / 2 1 , / 2) 1 f 1 g , and separated by Wannier gap-closings. Here the sequence of two Wyckoff position labels indicates one orbital in one position and another orbital in the second position. The transitions 1 a a 1 ↔ 1 1 c c and 1 1 b d ↔ 1 f 1 g are characterized by Wannier gap-closings at ν x = ν y = 0 and ν x = ν y = 1 2 / , respectively, whereas the

FIG. 17. Evolution of the Wannier spectrum W x with eigenvalues ν x ( k , k y z ) for different values of λ z and fixed λ 0 = 5 . , λ x = 1 , λ y = 25 . , and λ z = 1 . Across these that ν z is gapped, we find the gappless transition of ν x to happen when λ z = 1 . Panels (a) and (c) correspond to the gapped phases 1 g and 1 d in the phase diagram of Fig. 15, respectively.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00011v2-figure-26f2a159-2117-44a3-a754-c83ba575a4af.jpg)

transitions 1 a a 1 ↔ 1 1 b d and 1 1 c c ↔ 1 f 1 g are characterized by Wannier gap-closings at ν z = 0 and ν z = 1 2 / , respectively. As a result, the 1 1 b d phase is a BOTP that cannot be reached from any of the three other phases without closing a gap at the boundary for our chosen boundary termination (corresponding to a WCP µ x,y,z = 1 2 / ).

This can be seen more directly by studying the evolution of the Wannier bands. The Wannier bands in the z -direction calculated via W z are shown in Fig. 16b. The four bands are split into two groups of two bands. Akin to the C 2 h model, when λ M = λ 0 &gt; λ / S 2 , the two filled bands are located around 1 a with p z,x = p z,y = 0 . On the other hand, when λ M = λ x,y &gt; λ S / 2 , the two filled bands describe one electron at 1 f and one at 1 g , each having p z,x = p z,y = 1 2 / . In both cases, the Wannier bands form elementary bands and cannot be split further. In particular, this implies that the Wannier gap closing when λ M &lt; λ S / 2 involves all four bands.

Focusing on W x , we can similarly investigate the transition between the two distinct phases by tuning γ z . Here, we find a gap closing transition when | γ z | = | λ z | (c.f., Fig. 17). Unlike W z , the bands are split into four separate bands (four EBRs) since the symmetries of W x do not relate them. This follows from the fact that, with only mirror symmetry, all non-general Wyckoff positions have multiplicity 1. This implies that all EBRs that are not derived from the general position are onedimensional. As a result, the gap closing at the WCP takes place only between two of the four Wannier bands. This im- m = +0 5 . λ z = 1 5 . m = -0 5 . λ z = 1 5 . m = 0 λ z = 1 5 . (a) (b) (c) kx kx kx plies that we can identify the different phases by looking at the (nested) Wannier band polarization of a single band. The phase with hinge charge is characterized by p x,z = 1 2 / .

m = -0 5 . λ z = 1 0 . m = -0 5 . λ z = 1 5 . m = -0 5 . λ z = 0 5 . (d) (e) (f) ky kz ν z kz ν z kz ν z ky ky C 6 h model -We can define a model with C 6 h symmetry in a very similar fashion. The 2D Hamiltonian used to construct this model (via Eq. (29)) is schematically illustrated in the third panel of Fig. 13. The behavior of the model is qualitatively similar to the models above. Namely, there are two phases distinguished by a surface obstruction as shown in Figs. 16c and 17c. The non-trivial phase is characterized by a fractional hinge charge per unit cell (c.f. Fig. 14c).

x

ν

ν

x

ν

x

## 4. A hinge-obstructed 3D BOTP

kz kz kz One of the simplest 3D BOTPs we can construct is a trivial dimensional extension of the DMQI. Following the recipe of Sec. V A, we can create a 3D BOTP by stacking pairs DMQI's with π fluxes at each plaquette along the stacking direction. Since the DMQI is edge-obstructed and has corner charges, the resulting 3D BOTP has a hinge obstruction associated with a filling anomaly which translates to a corner charge of e/ 2 at each of the 8 corners.

The filling anomaly can be understood by noting the 3D model has mirror symmetries about three perpendicular planes, i.e., M ,M , x y and M z . When considered on a torus, i.e., with periodic boundary conditions, then there are 8 points q with maximal symmetry | G q | = 8 , where the three perpendicular mirror planes intersect. These points correspond to the 8 maximal Wyckoff positions of multiplicity 1. Whenever the filling of such points is a multiple of two, we can symmetrically move the electrons along the intersection of any two of the three mirror planes to reach any of the other symmetric positions. On the other hand, with open boundaries, the site q that lies on a corner of the open system has its site symmetry group reduced such that | F q | = 1 . At a bulk filling of ν pbc q = 4 (which is natural for a pair of DMQI's), the sites at the corners will have a filling of 1 / 2 signaling the existence of a BOTP filling anomaly. Such a filling anomaly is associated with a hinge rather than surface obstruction. Indeed, this model is the octupole insulator of Refs. 17 and 18, which was the first example of a 3D BOTP when protected by mirror symmetries (additionally, bulk-obstructions exist in this model when additional crystalline symmetries are considered). Since this model is hinge-obstructed, its boundary obstruction is encoded in the band representations of the nested Wannier bands rather than the WBR themselves.

## C. 3D BOTPs with gapless surface states

In the previous section, we discussed models with boundary obstructions associated with filling anomalies. Here, we construct two models with gapless hinge states due to a surface obstruction: the dimerized weak Chern insulator and the dimerized weak quantum spin Hall (QSH) insulator. Both of these models are constructed using the recipe of Sec. V A.

| k                          | Γ             | X       | Y     | M     |      |
|----------------------------|---------------|---------|-------|-------|------|
| ρ                          | w z + ( C 2 z | ) w z + | w z + | w z + | nu Γ |
| (1 a, { 1 ¯ E, 2 ¯ E } )   | +             | +       | +     | -     | 2    |
| ( 1 c , { 1 ¯ E, 2 ¯ E } ) | +             | +       | +     | -     | 0    |
| (1 a, { 1 ¯ E, 2 ¯ E } )   | +             | +       | +     | +     | 2    |
| (1 c, 1 ¯ E, 2 ¯ E )       | +             | +       | +     | +     | 2    |

{

}

TABLE III. Example of Wannier symmetry representations in the magnetic space group P 2 /m ′ generated by the irreducible representation of the 2 /m ′ point group { 1 ¯ E, 2 ¯ E } at different maximal Wyckoff positions allowed through different parameter choices in the model of Eq. (40). In this case M z is not a symmetry and does not label the phase. Top rows: Phases with These correspond to boundary obstructed phases with anomalous surface states, where the Wannier bands in W z are characterized by a Chern number. The bold line corresponds to the phase with a boundary obstruction when the conventional Wannier chemical potential is chosen, i.e., µ = (1 2 1 / , / , 2 1 / 2) . Neither the positions 1 a or 1 c admit a filling anomaly in this space group (See Fig. 15) however the bulk representation ρ admits a restriction the Wannier bands which are topological or in other words nonlocal. Such configuration makes the position 1 c obstructed. The Wannier configuration on the top rows are characterized by a Chern number and the obstructed cases will host anomalous chiral modes at the boundary. Note that this phase, unlike the models with filling anomalies presented above, is not uniquely indicated by the symmetry eigenvalues of the Wannier bands, but it is still diagnosed by the representation of the open boundary system.

Westart by considering a 2D Chern insulator whose Hamiltonian is given by

$$\mathcal { H } _ { \text{Cl}, \pm } ( k ) = \pm [ & \sin k _ { x } \sigma _ { 1 } + \sin k _ { y } \sigma _ { 2 } \\ & \quad + ( 2 + m - \cos k _ { x } - \cos k _ { y } ) \sigma _ { 3 } ]. \quad ( 3 7 ) \quad \mathcal { H } = \\ \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot$$

Let m be restricted to the interval -2 &lt; m &lt; 2 . In this case, the Hamiltonian H CI , ± has Chern number C = ∓ 1 for m&lt; 0 (with the sign determined by the sign in Eq. 37), Chern number 0 for m &gt; 0 , and it is gapless for m = 0 . The 3D dimerized Chern insulator is built by stacking copies of the 2D Chern insulator above with alternating sign which results in two layers per unit cell with alternating Chern number. Then we couple the layers with a dimerized coupling that alternates between two values λ z and γ z . Physically, the terms multiplying sin k x,y and cos k x,y correspond to hopping between unit cells, so flipping their sign is equivalent to threading a π -fluxes through all the vertical plaquettes. On the other hand, the term 2+ m denotes an onsite potential whose sign should be chosen to alternate between layers. The resulting 3D Hamiltonian is given in Eq. 29. This Hamiltonian is invariant under the product M z T = σ τ 1 1 K , and C 2 z = iσ 3 . This set of symmetries places our model in the magnetic space group P /m 2 ′ (10.45), which has the point group 2 /m ′ that admits a single spinful complex representation { 1 ¯ E, 2 ¯ E } [71].

The energy spectrum can be understood by noting that the 3D Hamiltonian has the form of a Dirac Hamiltonian

## Dimerized weak Chern insulator

Dimerized weak quantum spin-Hall insulator

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00011v2-figure-48077b90-fe39-4396-8209-cb1539665756.jpg)

FIG. 18. Three dimensional boundary obstructed phases without a Wannier representation: dimerized weak Chern insulator protected by M z T and C 2 z (top panel), and the dimerized weak quantum spin Hall insulator protected by M z , T , and C 2 z (bottom panel). Left: Scheme of the model where the arrows indicate chiral onedimensional modes that can be gapped by a coupling along z . In the two nontrivial cases, either a 2D Chern insulator or a 2D quantum spin Hall insulator is left unpaired at the surface. Center: Energy spectrum of a system with translation symmetry along y but open along both x and z . The ingap states are located at the top and bottom hinges. We have marked in black the band that is exponentially localized at one hinge. Right: illustration of the 1D surface states on a spherical surface termination. The equator and north/south pole are invariant under M z and C 2 z , respectively. The 1D surface states can only be symmetrically removed by closing the gaps at one of these high symmetry regions.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00011v2-figure-10d02c35-b34e-40ba-ab16-50a5c35f86d4.jpg)

$$\mathcal { H } = \sin k _ { x } \Gamma _ { 1 } + \sin k _ { y } \Gamma _ { 2 } + ( 2 + m - \cos k _ { x } - \cos k _ { y } ) \Gamma _ { 3 } \\ \text{this} \quad & + \lambda _ { z } \sin k _ { z } \Gamma _ { 4 } + ( \gamma _ { z } + \lambda _ { z } \cos k _ { z } ) \Gamma _ { 5 } \quad ( 3 8 )$$

with Γ 1 2 3 , , = σ 1 2 3 , , τ 3 , and Γ 4 5 , = τ 2 1 , such that C 2 z = Γ Γ 1 2 , and M z T = Γ Γ 1 4 K . It is then obvious that for this Hamiltonian to be gapped all 5 terms have to vanish simultaneously which is only possible for m = 0 and | γ /λ z z | = 1 . This set of points where the Hamiltonian is gapped forms a line in the three dimensional parameter space which means that it is always possible to connect any set of parameters ( m,γ ,λ z z ) for which the Hamiltonian is gapped without going through a bulk gap closing. That is, the bulk phase is always trivial. We note that we can also understand the triviality of the model from the fact that none of the symmetries C 2 z and M z T quantizes the θ angle in the axion term θ E B · [49, 61, 72].

In order to study this transition, we compute the Wannier spectra in Fig. 19. We can see in the top panel that changing the sign of m induces a gap-closing in the Wannier spectrum in the z -direction indicating gap-closing in the top/bottom surfaces. On the other hand, keeping the value of m fixed and negative, we find a gap-closing transition in the Wannier spectra along the x or y -directions (bottom panel). This

(a)

C

2

(d)

C

4

(g)

C

6

ν

means that for the chosen termination, the model exhibits a boundary obstruction separating the gapped phase with parameters ( | γ z | &gt; λ | z | , m &lt; 0) from ( | γ z | &lt; λ | z | , m ≶ 0) and ( | γ z | ≶ | λ z | , m &gt; 0) . ky x λ z = 0 5 . λ z = 0 5 . λ z = 0 5 . ky ν x ky ν x

ν

ν

An interesting aspect of this model is that the Wannier obstruction for W z has a different origin from W x and W y . In the Wannier spectrum of W z , the obstruction is associated with a non-zero Chern number, i.e., the Wannier bands do not separately form WBRs, as illustrated by the non-trivial spectral flow of the nested Wannier spectrum shown in Fig. 20a and e. On the other hand, the Wannier obstruction in W x or W y is associated with changes in the nested Wannier band polarization between values of 0 and 1 / 2 as shown in Fig. 20c, d, g, and h. Unlike the filling anomaly scenario, the obstruction in this model is not due to the change in multiplicity of a Wyckoff position when going between G and F, but rather the fact that generating the open boundary representation from the site 1 c is not compatible with a local representation, i.e., it cannot be represented as the symmetry action on localized atomic orbitals of the open system. . This is shown in Table III. When the Chern layers are in the topological phase (upper rows), the two Wannier bands of the Wannier spectrum W z form an irreducible representation of G z = { C 2 z , T x , T y } when grouped together. However, the separate Wannier bands are topological and their eigenvalues are not compatible with an atomic insulator. If the Wannier bands are both located around the center of the unit cell ( 1 a ), then this configuration is not obstructed, since for each unit cell the two Wannier bands of opposite Chern number combine to form a band representation whose local support lies on a site with site symmetry group F admitting the two dimensional irreducible representation. However, if the bulk band representation is generated from the position 1 c , with site symmetry group C 2 , which admits only one dimensional representations, it is necessarily obstructed for a boundary where 1 c lies on the boundary. Intuitively, the Chern bands of opposite Chern number centered around the edge of the unit cell are spatially separated leaving behind chiral modes exponentially localized to the hinges, which can either be removed by closing the gap of the Chern bands and inverting their mass (closing the gap on the poles, see Fig. 18), or changing the coset decomposition of G such that it can be generated from a Wyckoff position in the z = 0 plane, resulting in the surface gap closing at the equator. kz λ z = 1 λ z = 1 5 . λ z = 1 λ z = 1 5 . λ z = 1 λ z = 1 5 . kz kz kz ky x kz ky ν x kz ky ν x kz ky x kz ky ν x kz ky ν x (b) C 2 (e) C 4 (h) C 6 (c) C 2 (f) C 4 (i) C 6

It is important to highlight the role played by the C 2 z symmetry here. Although it is not needed as a traditional protecting symmetry for the Wannier topology, it is crucial to single out the upper and lower surfaces as high symmetry planes where Wannier gap-closings imply a surface phase transition. Similar to the discussion of Sec. III A, this can be illustrated by considering a spherical geometry as shown in Fig. 18c. Under M z T alone, only the equator is a high symmetry region, and we are generally allowed to smoothly deform the chiral hinge modes by shrinking them to a point without touching the equator. In the presence of C 2 z , however, the north and south poles become HSSs where surface gap-closings indicate a surface phase transition per our definition. Note that allowing for distinctions captured by gap-closings at generic, nonHSSs will generate the same spurious, uncountable boundary

FIG. 19. Evolution of the W z (a-c) and W x (d-f) Wannier spectra of the dimerized weak Chern insulator model for different values of m and λ z with fixed γ z = 1 . Entering the phase with m = -0 5 . and λ z = 1 5 . requires a Wannier gap closing at ν = 1 2 / both when tuning m or λ z . This phase is therefore separated by a boundary phase transition for a boundary with the conventional Wannier chemical potential µ x,y,z = 1 2 / .

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00011v2-figure-7c5543b4-7ea6-441e-bc1e-97fc9a7953d8.jpg)

FIG. 20. Nested Wannier polarization p z,x ( k y ) and p x,z ( k y ) for different values of m and λ z with fixed γ z = 1 in the dimerized weak Chern insulator (a-d) and weak quantum spin-Hall insulator (e-h). The phase with m = -0 5 . and λ z = 1 5 . (a,e) is characterized by a winding of the nested polarization which indicates that the Wannier bands in W z do not form band representations and are associated with hinge modes. On the other hand, p x,z ( k y ) is gapped with a different average value in the two phases, either 0 (c,g) or 1 2 / (d,h), reflecting the fact that W x decomposes into distinct Wannier band representations.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00011v2-figure-9296fc22-3a7a-4345-b79d-9f26fb1fe765.jpg)

distinctions as those discussed in Sec. III A.

## 2. Dimerized weak quantum spin Hall insulator

Similar to the dimerized Chern insulator, we can define a 3D BOTP by stacking layers of the 2D quantum spin Hall

(QSH) topological insulators with dimerized coupling. We start with the 2D Bernevig-Hughes-Zhang (BHZ) model [3] whose Hamiltonian is given by

$$\mathcal { H } _ { \text{BHZ} } ( k ) & = \sin k _ { x } \sigma _ { 1 } s _ { 3 } + \sin k _ { y } \sigma _ { 2 } \quad & \underset { \text{indica} } { \sim } \\ & \quad + [ 2 + m - \cos k _ { x } - \cos k _ { y } ] \sigma _ { 3 } \quad ( 3 9 ) \quad \text{eigen} \\. & \quad. & \quad. & \quad. & \quad. & \quad \text{end} \quad \text{$\Phi$}$$

where s 1 2 3 , , denote the Pauli matrices in spin space. We also take | m &lt; | 2 such that the Hamiltonian describes a topological insulator for m &lt; 0 and a trivial phase for m &gt; 0 . The Hamiltonian is invariant under spinful time-reversal symmetry T = is 2 K . It also has U(1) S z rotation symmetry, but that is not required for its topological protection [1, 2]. The Hamiltonian is also invariant under 2D inversion ( I 2D = σ 3 ) which, when combined with π S z rotation, yields the spinful two-fold rotation C 2 z = iσ s 3 3 .

We can then build a 3D BOTP by stacking layers of the 2D Hamiltonian (39) with alternating signs . As explained earlier, this amounts to threading π fluxes threaded in the vertical plaquettes and choosing the sign of the parameter 2+ m to alternate between layers. The resulting 3D Hamiltonian has the form (29) with H 2D = H QSH . Due to the Z 2 topology of the QSH, switching the sign of the Hamiltonian between layers does not alter the topological index, but it is needed to ensure that the bulk Hamiltonian describes a BOTP rather than a weak TI. This distinction will be encoded in the symmetry action as we will see now. In addition to time-reversal and C 2 z , the 3D Hamiltonian is also invariant under M z = iτ 1 σ s 1 1 . Similar to Secs. II and V B, we can think of the symmetry representations considered here as a projective representation which combines spatial symmetry transformations with internal and gauge transformations. We note that the 2D inversion symmetry I 2D anticommutes with M z as required by the presence of the π fluxes in the vertical plaquettes, whereas the spinful C 2 z rotation commutes with M z . Thus, the projective representation of the symmetry group generated by C 2 z and M z (together with translations) can be thought of as a nonprojective spinful (double) representation of the space group 10 (P2/m) (if we allow for perturbations which break S z spin rotation but not C 2 z or M z ).

As in the dimerized weak Chern insulator, the Hamiltonian has the form of a Dirac Hamiltonian

$$\mathcal { H } = \sin k _ { x } \Gamma _ { 1 } & + \sin k _ { y } \Gamma _ { 2 } + ( 2 + m - \cos k _ { x } - \cos k _ { y } ) \Gamma _ { 3 } \quad \text{try, ea} \\ & + \lambda _ { z } \sin k _ { z } \Gamma _ { 4 } + ( \gamma _ { z } + \lambda _ { z } \cos k _ { z } ) \Gamma _ { 5 } \quad ( 4 0 ) \quad \text{polar} \text{:} \\ & \quad \text{the K} \text{,}$$

with Γ 1 = σ s τ 1 3 3 , Γ 2 3 , = σ 2 3 , τ 3 , Γ 4 5 , = τ 2 1 , , and Γ 6 7 , = τ 3 σ s 1 2 1 , with C 2 z = Γ Γ 1 2 , T = Γ Γ 1 7 K , and M z = Γ Γ 4 7 . Similar to the dimerized weak Chern insulator, the Hamiltonian above is only gapless for m = 0 and | λ /γ z z | = 1 which means that all the parameters for which the Hamiltonian is gapped can be deformed to one another without closing the gap. One consistency check for the triviality of this Hamiltonian is to note that the Hamiltonian has the inversion symmetry I = C 2 z M z = τ 1 σ s 2 2 which squares to +1 . Thus, one can check its triviality by computing the symmetry indicators [23, 31] obtained from the inversion eigenvalues at the timereversal invariant momenta (TRIM). These correspond to a Z 4

index labelling the possible bulk phases (strong TI or secondorder TI), and three weak Z 2 invariants for the x , y , and z directions [23, 31]. There are four filled bands which come in pairs of Kramers' doublets at the TRIMs. The symmetry indicators are obtained by computing the product of inversion eigenvalues at each TRIM (counting only one state out of each Kramers' doublet) [23, 31, 73]. By an explicit calculation, we find that the product of inversion eigenvalues at every TRIM is equal to -1 . As a result, all symmetry indicators are trivial.

It is important to highlight the role of π -fluxes through the vertical plaquettes. Unlike the models considered in Sec. V B, it is possible here to choose a gauge such that M z and C 2 z commute. This arises from the fact that the total flux per plaquette in the vertical direction (including both spins) vanishes, although the there is a π -flux per spin. Thus, the the π -fluxes do not alter the commutation relations between M z and C 2

We can investigate the boundary obstruction numerically by computing the Wannier bands. Here again we consider a slightly modified model with the same topological properties but on a cubic lattice with the Hamiltonian given by

$$\text{the} \quad \mathcal { H } = ( 2 + m - \cos k _ { x } - \cos k _ { y } ) \Gamma _ { 1 } + \sin k _ { x } \Gamma _ { 2 } + \sin k _ { y } \Gamma _ { 7 } \\ \text{/ers} \quad \mathcal { \quad } + \lambda _ { z } \sin k _ { z } \Gamma _ { 4 } + ( \gamma _ { z } + \lambda _ { z } \cos k _ { z } ) \Gamma _ { 3 }. \quad ( 4 1 )$$

The Wannier spectra of the model are very similar to the weak Chern insulator since this model is just two time-reversed copies of the Chern insulator model. Thus, we see a Wannier transition at the WCP in W z when | λ z | &gt; γ | z | , and m = 0 ; as well as a Wannier transition at the WCP in W x and W y when m&lt; 0 but | γ z | = | λ z | . However, in contrast to the dimerized weak Chern insulator, this model has four filled bands. Each of the Wannier bands is two-fold degenerate at time-reversal invariant momenta, and the additional C 2 z symmetry causes this degeneracy to appear across the entire surface Brillouin zone.

The nested Wannier spectra of this model are shown in Fig. 20. We note that ν z,x ( k y ) resembles the spectral flow in a quantum spin Hall insulator. Two Wannier bands flow in opposite directions, crossing at k y = 0 and π where the degeneracy cannot be lifted due to Kramers' theorem. This nested Wannier winding captures the dangling QSH at the upper and lower surfaces. On the other hand, ν x,z ( k y ) in the nontrivial phase shows two bands related by time-reversal symmetry, each with average polarization of 1 / 2 . Although the total polarization vanishes, the polarization for each component of the Kramers' pair does not, and it is such a quantized (timereversal) polarization [74] that captures the Wannier transition in ν x ( k , k y z ) . In general one could use a description of the time-reversal polarization in terms of a Pfaffian invariant[74], but we will not pursue that further here. In this case the Wannier transition resembles the transition in a spinful SSH chain with time-reversal and inversion symmetry where the edge anomaly is associated with Kramers' pairs of electrons rather than an individual electron.

In summary, both the dimerized weak Chern model and the dimerized weak QSH model exhibit one-dimensional propagating hinge states that can be removed without closing the bulk gap. However, removing the hinge states requires clos-

ing the energy gap on at least one high-symmetry surface. These systems only admit a Wannier representation for periodic boundaries, but not for open boundaries. This serves to illustrate that boundary obstructions are not restricted to Wannier representable systems.

## VI. DISCUSSION

In this work, we introduced the notion of boundary topological obstructions that captures obstructions which do not exist on periodic boundaries, but exist on geometries with symmetric open boundaries. In more precise terms, we call two Hamiltonians H 1 and H 2 boundary-obstructed on a given symmetric open boundary if they can be symmetrically deformed to each other without closing the gap with periodic boundary conditions, but not for the open boundary. Another way to formulate this is by saying that any symmetric deformation of H 1 to H 2 involves a gap-closing at a highsymmetry surface rather than in the bulk , as in the case of conventional SPTs. Although the possibility of an obstruction between a given pair of Hamiltonians is termination dependent, the group of equivalence classes distinguished by boundary obstructions is termination-independent and is characterized by physical boundary signatures such as filling anomalies and gapless states.

It is worth stressing that boundary signatures associated with BOTPs are not anomalous as they can be removed by adding a lower-dimensional SPT on the boundary. As a result, we need to be careful when defining the boundary termination to distinguish surface features which arise from the bulk from those due to the termination. This is done by first considering a simple reference class of terminations achieved by 'interpolating" the bulk Hamiltonian to the vacuum Hamiltonian by means of a single real scalar function. This is always possible by changing the value of the chemical potential from the inside to the outside of the sample. In a more general boundary termination, we can always separate the contribution coming from the termination by comparing to this reference boundary (cf. Eq. 17). This enables us to isolate the boundary contribution to the open Hamiltonian making it possible to compare two open Hamiltonians with the same boundary .

We also offered an alternative definition of boundary obstructed phases using the theory of band representations. In this definition, BOTPs can be identified with equivalent band representations of the space group in the periodic system which become inequivalent upon restricting to a point group of the open system. This definition does not only bridge the idea of boundary obstruction to the bulk band representation, it also enables us to define an absolute notion of boundary obstruction, related to the existence of filling anomalies and gapless surface states, from the knowledge of a relative obstruction as defined above. This is achieved by observing that whenever two band representations which are equivalent in the periodic system are inequivalent in the open system, at least one of them has to be non-local and exhibit a filling anomaly or gapless surface states. This enables us to identify a trivial BOTP phase (the one corresponding to a local representation in the open system) such that all phases with a boundary obstruction relative to it are non-trivial.

The motivation for introducing the concept of BOTP is twofold. First, there is the advent of new ways to implement 2D topological systems, e.g., in synthetic materials [27, 7587] and cold atoms [88-95], which can be used to realize the DMQIandsimilar models not associated with a bulk topological invariant. The observation and robustness of corner modes in these systems motivate developing a concept that goes beyond standard SPTs to understand the features of these models responsible for these observables' stability. We note that these systems are characterized by a large degree of control over the boundary termination, which is usually done by sharply terminating the lattice or through a confining potential. Even for real crystals where there is less control regarding the termination, we expect our concept to be relevant to understanding surface states which may arise through tuning the bulk parameters by applying pressure, strain, magnetic field, etc. or by applying a spatially dependent gate voltage inside the sample to implement an edge termination as described in Eq. 17.

The second motivation is conceptual. Since the introduction of the concept of higher-order topological phases [17, 18, 20], there has been a large body of works that discussed models with corner or surface states (see Ref. [96] for a recent review). In many cases, the existence of surface states was established by numerically analyzing the Hamiltonian on a simple boundary without making the distinction whether these states arise from a bulk invariant or not. To distinguish these, the concept of intrinsic vs. extrinsic HOTIs was introduced in Ref. [13, 24] with the former corresponding to phases whose surface states are anomalous and cannot be removed with symmetric boundary addition, whereas the latter corresponding to systems whose surface states are removable through symmetric boundary additions. However, the concept of extrinsic HOTIs, while useful, is too broad since it lumps together models like the DMQI, which has robust boundary features when considered on most simple boundaries and which can be related to bulk quantities such as those described in this text, with models consisting of a completely trivial bulk attached to a lower dimensional SPT on the surface. The concept of BOTP helps make this distinction and explains why some models exhibit some topologically robust surface features without being associated with a bulk invariant.

In Sec. II, we used the double-mirror quadrupole insulator (DMQI) of Ref. [17, 18] as a prototypical example to illustrate the notion of BOTPs. Although this model has been intensively studied as one of the first examples of a higherorder topological insulator, one crucial aspect of it had so far been overlooked. Namely, the fact that in the absence of C 4 z symmetry, the two 'phases' of the model are not SPTs in the standard sense since they can be deformed to each without closing the bulk gap in periodic boundary conditions. Using the notion of BOTP, we provide a resolution to this problem by showing that such 'phases' are related by a boundary rather than bulk obstruction. This was investigated in detail from several different perspectives including real space orbital deformations, Wannier spectra, and symmetry representations.

The latter was particularly useful in establishing the connection between the boundary obstruction in the model and the existence of a filling anomaly. Indeed the existence of a filling anomaly allows for a topological boundary signature that does not rely on the existence of zero-energy corner states (which, for insulators, typically require fine-tuned symmetries such as particle-hole or chiral symmetry). Instead, it reflects the fact that the model, when filled with electrons to neutrality, cannot be simultanuously gapped, symmetric and charge-neutral. Thus, a symmetric, gapped DMQI at half-filling will necessarily have excess/deficit charge of 2 e equally distibuted among the four corners yielding a fractional corner charge of e/ 2 .

In Sec. V, we introduced several 3D models for BOTPs. The first family of such BOTPs have similar phenomenology to the DMQI and exhibit surface obstructions. The non-trivial BOTPs in this family are associated with filling anomalies similar to the one in the DMQI, and they manifest in fractional hinge charge (per unit cell) when the system is symmetric and gapped at half-filling. In addition, we introduced a second family of BOTPs characterized with chiral/helical propagating hinge modes. In both cases, the existence of a boundary obstruction is established explicitly by studying the Wannier spectra and the symmetry indicators of the corresponding WBRs.

Before closing, let us discuss the possible extension of the concepts discussed here beyond free fermions. We note that the existence of chiral/helical hinge states or filling anomalies are robust features that are expected to survive in the presence of interactions. Thus, we expect the boundary topological distinctions we found here to still be relevant for interacting systems. On the other hand, defining boundary obstructions is likely more tricky for interacting phases since it requires the implementation of a surface termination in a way which isolates the surface and bulk contributions. It is unclear whether this is generally possible which makes it difficult to decide what it means to 'keep the boundary fixed," or to compare two phases with the 'same boundary." We leave the investigation of such questions to future works.

## ACKNOWLEDGMENTS

We thank Barry Bradlyn, Jennifer Cano, Dominic Else, Adrian Po, Ryan Thorngren, and Ashvin Vishwanath for stimulating discussions. E. K. was supported by a Simons Investigator Fellowship, by NSF-DMR 1411343, and by the German National Academy of Sciences Leopoldina through grant LPDS 2018-02 Leopoldina fellowship. W. A. B. thanks the support of the Eberly Postdoctoral Fellowship at the Pennsylvania State University. T. L. H. thanks US National Science Foundation (NSF) grants EFMA-1627184 (EFRI), and DMR1351895 (CAREER) for support. R. Q. was funded by the Deutsche Forschungsgemeinschaft (DFG, German Research Foundation) - Projektnummer 277101999 - TRR 183 (project B03), the Israel Science Foundation, and the European Research Council (Project LEGOTOP).

Author Contributions: E. K. developed the general definition for boundary obstructions, the general correspondence between Wannier and boundary topology, proposed 3D model examples, aided the development of the real-space, band representation, and filling anomaly approaches to BOTPs, and was the primary manuscript writer. W. A. B. developed the connection between real-space Wannier center deformation and boundary obstructions, helped develop the Wannier band representation, real-space, and filling anomaly approaches to BOTPs, and helped in the writing of the manuscript. T. L. H. contributed to the refinement of all of the conceptual advances of the paper and helped in the writing of the manuscript. R. Q. developed the formalism of band representation theory of BOTPs, helped develop the real-space and Wannier spectrum approaches to BOTPs, developed the filling anomaly arguments in BOTPs, helped in the writing of the manuscript, and led the initiation of the collaboration.

## Appendix A: Character table of D π 2

TABLEIV. Character table of the point group D π 2 defined in the main text. The point group is isomorphic to D 4 in the absence of a π flux.

| Rep. / class   |   { 1 } { C |   2 2 z } { C |   2 z } { M y | } { M x }   |
|----------------|-------------|---------------|---------------|-------------|
| A 1            |           1 |             1 |             1 | 1 1         |
| A 2            |           1 |             1 |             1 | -1 -1       |
| B 1            |           1 |             1 |            -1 | 1 -1        |
| B 2            |           1 |             1 |            -1 | -1 1        |
| ¯ E            |           2 |            -2 |             0 | 0 0         |

## Appendix B: Real Space Proof of Anti-commuting Symmetries

## 1. Anticommuting Mirrors in the 2D DMQI

We now show that in order to carry out the deformation between the configurations at 1 a and 1 d in the DMQI, it is necessary that the two mirror operators anticommute. This is suggested by the fact that the two configurations actually correspond to different symmetry representations for commuting mirrors, but the same symmetry representation for anticommuting mirrors ( ¯ E ). We will see that the anticommuting mirrors can be inferred from the real space picture in Fig. 4. Let us start with the orbitals at position 1 a and then move them horizontally using position 2 e into position 1 c. The two orbitals in position 2 e are related by mirror symmetry M x and have the form | ( ± x, 0) 〉 , where | r 〉 denotes an orbital localized at point r . In this basis, M x is off-diagonal and acts generally as M x | ( ± x, 0) 〉 = e ± iφ | ( ∓ x, 0) 〉 (with some arbitrary phase φ ), leading to the eigenvectors | (+ x, 0) 〉± e iφ | ( -x, 0) 〉 whose eigenvalues are ± 1 respectively. Since these orbitals lie on the y = 0 mirror line, they are eigenvectors of M y . Next, to bring these orbitals to position 1 d , we move them vertically using position 2 h . In this position, the action of M y is off-diagonal, thus, following the same argument as for M x , we can deduce that it has two distinct eigenvalues +1 and -1 . Since the symmetry eigenvalues cannot change under

smooth, symmetry-preserving deformations, we deduce that M y = ± σ 3 in the 2 e basis. In the same basis, we showed that M x = σ 1 exp { iφσ 3 } , which in turn implies { M ,M x y } = 0 . Since anticommutation is basis-independent, we deduce that connecting the 1 a and 1 d positions at a filling of two electrons per cell requires M x and M y to anticommute. We note that the connection of mirror anticommutation in the DMQI model to important quadrupole properties was proposed in Refs. [17, 18] based on a study of the Wannier spectrum.

## 2. Anticommuting C 2 z and M z for 3D C 2 nh Models

To show that C 2 z and M z anticommute for the C 2 nh models that we discuss in Section V B 1, we note that in order to occupy a movable Wyckoff postion that interpolates between 1 a and nc in the M z invariant planes, the Wannier centers permute into each other under the action of C 2 nz . This implies the eigenvalues of the C 2 nz operator must span all the 2 n roots of unity e πil/n with l = 0 , . . . , 2 n -1 . Similarly, to occupy a vertically movable Wyckoff position at C 2 z invariant lines, the eigenvalues of M z should be ± 1 since its action is a permutation of the two sites in this position. Noting that the two electrons at position nc originate from electrons which have moved away from the center 1 a in diametrically opposite directions, and are thus related by C 2 z , we find that the requirement of deformability in the vertical direction forces any pair of C 2 z -related orbitals to have opposite mirror eigenvalues. Thus, in a basis where the mirror symmetry M z is diagonal, C 2 z is purely off-diagonal and it switches the positive and negative mirror sectors which implies the anticommutation condition

$$C _ { 2 z } M _ { z } C _ { 2 z } ^ { \dagger } = - M _ { z }. \quad \quad \ ( \mathbf B 1 ) \quad \text{du} ( \text{$\mathbf $the$} )$$

This discussion is very similar to the one for the DMQI above.

## Appendix C: Band representations of BOTPs

## 1. Induction of a band representation from a local representation

Restricting ourselves to atomic bands, it follows from the work of Zak [40] that one can write a representation ρ of a space group G from two ingredients. The first ingredient is the Wyckoff position Q , comprised of a site q with a site symmetry group G q and the sites q i obtained by the coset expansion of G with respect to G q . Explicitly, the full space group can be decomposed by G = ∑ i g G i q with g i ∈ G/G q . The lattice sites q i in the crystal are obtained by the left action of the coset representatives q i = g i q , with i = 1 , ..., N Q with N Q the multiplicity of Q , the number of times the site q is repeated over the lattice. The second ingredient needed to define ρ is the representation ρ q of G q under which the state | q 〉 transforms. That is, for any h ∈ G q , h | q 〉 = ρ q ( h ) | q 〉 . Finally, the collection of the representations ρ q i and the associate basis functions | q i 〉 , form a complete basis for the elec- tronic states in the system, and the complete representation ρ of the space group G , acts on this space. Naturally, ρ is not necessarily diagonal in this basis since it mixes orbitals at different lattice sites q i . Making use of translation symmetry in G , and the fact that the crystal momentum k diagonalizes lattice translations T , it is possible to express ρ in such basis, where the diagonal entries ρ k represent Bloch states labelled by k with little group G k , at a Brillouin zone point k . The collection of ρ k at high symmetry momenta determine the band representation [40, 41, 64, 65].

The construction of the band representation is straight forward. Let us start with the representation ρ q of G q for any element h ∈ G q , h | q 〉 = ρ q ( h ) | q 〉 , here | q 〉 denotes a basis state localized at site q . Other basis states | q i 〉 are obtained by the action coset elements g i as | g i q 〉 . Any group element g can be expressed as a product of an element of G q and elements of the coset as g = g hg j -1 i or alternatively g h j = gg i with h and g j uniquely determined by g and g i . Then, it follows that the action of any element g ∈ G on any basis state is given by

$$g | g _ { i } q \rangle = g _ { j } h | q \rangle = \rho _ { q } ( h ) | g _ { j } q \rangle. \quad \quad ( C 1 )$$

It follows that a symmetric state in the full Hilbert space transforms into itself under all elements of G under a (reducible) induced representation ρ of G , which in the real space basis | q i 〉 is given by

$$\Gamma _ { \beta } \quad [ \rho ( g ) ] _ { i j } = \rho _ { q } ( g _ { j } ^ { - 1 } g g _ { i } ) \text{ if } g _ { j } ^ { - 1 } g g _ { i } = h \in G _ { q }, \quad ( \mathbb { C } 2 )$$

and zero otherwise. This representation is nothing but the local representation ρ q induced to the entire Wyckoff position Q , and it is usually abbreviated as ρ = ρ q ↑ G . The induced representation is unique up to the choice ambiguity of the coset representatives g j .

While the statement above is completely general, we can take a shortcut in symmorphic space groups by inducing the representation from a site q with the point group G q = F as its site symmetry group. In this case, all the degrees of freedom in the unit cell collapse into the dimension of the local representation at q . The action of translations on the center of the unit cell O (fixed to coincide with the center of the lattice with open boundaries) defines the Bravais lattice which may or may not coincide with the lattice Q obtained by the action of translations on q . Given that G = F × T , the different lattices correspond to distinct coset decompositions G = ∑ i t F i . For each space group G there are usually more than one site q of maximal symmetry, where its site symmetry group is isomorphic but not the same. They can be distinguished by the symmetry action on the basis states and cannot be connected under conjugation by any element in G . Explicitly in the DMQI model, these points generate the Wyckoff positions 1 a to 1 d . The space group p mm 2 admits four distinct coset decompositions up to conjugation by an element in G , with four isomorphic factor groups F = G/T .

Directly applying the induction procedure we find that an element of G acts on the maximal lattice as g g | i q 〉 = [ ρ g ( )] ij | g j q 〉 with [ ρ g ( )] ij = ρ q ( g -1 j gg i ) if g -1 j gg i ∈ G , q

and zero otherwise. In order to describe the representation of Bloch states in the Brillouin zone it is convenient to change the basis | q i 〉 to the crystal momentum basis | k 〉 related by the basis transformation 〈 q i | k 〉 = exp( -i k · q i ) . After this Fourier transform we have g | k 〉 = [ ρ g ( )] kk ′ | k ′ 〉 with diagonal terms given explicitly by

$$\rho _ { k } ( g ) = \sum _ { i j } \rho _ { q } ( g _ { j } ^ { - 1 } g g _ { i } ) e ^ { i k \cdot ( g _ { j } - g _ { i } ) q }. \quad ( C 3 ) \quad \text{$Since$} \quad \text{$\logi$} \quad \text{$\text{$alyz$}$} \quad \text{$\text{$\text{$\text{$\text{$i$}$}$}$} \quad \text{$\text{$\text{$alyz$}$}$} \quad \text{$\text{$\text{$\text{$\text{$i$}$}$}$} \quad \text{$\text{$\text{$\text{$alz$}$}$}$} \quad \text{$\text{$\text{$alz$}$}$} \quad \text{$\text{$\text{$alz$}$}$} \quad \text{$\text{$\text{$alz$}$}$} \quad \text{$\text{$\text{$alz$}$}$} \quad \text{$\text{$\text{$alz$}$}$} \quad \text{$\text{$\text{$alz$}$}$} \quad \text{$\text{$\text{$alz$}$}$} \quad \text{$\text{$\text{$alz$}$}$} \quad \text{$\text{$alz$}$} \quad \text{$\text{$alz$}$} \quad \text{$\text{$alz$}$} \quad \text{$\text{$alz$}$} \quad \text{$\text{$alz$}$} \quad \text{$\text{$alz$}$} \quad \text{$\text{$alz$}$} \quad \text{$\text{$alz$}$} \quad \text{$\text{$alz$}$} \quad \text{$\text{$alz$}$} \quad \text{$\text{$alz$}$} \quad \text{$\text{$alz$}$} \quad \text{$\text{$alz$}$} \quad \end{$$

When ρ corresponds to an elementary band representation it is either connected in momentum space or disconnected into topological subsets, which separately are not compatible with any tabulated band representation [55].

## 2. Movable Wyckoff positions

Two representations ρ and ρ ′ are equivalent if and only if there is a unitary and smooth matrix-valued function S t, g ( ) respecting S (0 , g ) = ρ g ( ) and S (1 , g ) = ρ ′ ( g ) for all symmetry group elements g . In a crystalline system, the unitary matrix that generates these changes S t, g ( ) = U t ( -t , g ′ ) † S t , g ( ′ ) U t ( -t , g ′ ) , continuously implements a basis transformation of ρ . As discussed in Ref. [65], when ρ is an induced representation from a point q into G, and ρ ′ is an induced representation from a different site q ′ , the S matrices can be explicitly constructed by inducing a family of representations from intermediate sites p that continuously interpolate between q and q ′ . A Wyckoff position P with continuously tunable sites p is called movable or non-maximal when it can continuously preserve G p .

This structure is the mathematical background for the pictures presented in the main text: S t, g ( ) is the induced representation from the states located at p in a path that connects q and q ′ , and guarantees there exists a gapped, symmetrypreserving path between the different Wannier configurations located at q and q ′ . We can try to directly find the unitary transformation U (1 , g ) that transforms ρ g ( ) into ρ ′ ( g ) . If ρ and ρ ′ are band representations, labelled by k , this implies finding a unitary matrix that is periodic in k and satisfies

$$\rho _ { k } ( g ) = U _ { g k } ^ { \dagger } ( 1, g ) \rho _ { k } ^ { \prime } ( g ) U _ { k } ( 1, g ) \quad \ \ ( \mathbb { C } 4 ) \quad \log \frac { 1 } { 2 }$$

for all g and k . While we have argued that all representations along a movable path are equivalent, the Wannier states | R ; q 〉 and | R ; q ′ 〉 are generally different. Additionally, the Hamiltonian itself is not generally left invariant by these transformations, and the path is implemented by varying the Hamiltonian parameters. However, in some cases two different points along the path may simply correspond to a basis transformation . In these cases U k leaves the Hamiltonian invariant. This is relevant for the DMQI example, where the mirror operators anticommute, and the Wannier states can be either at the 2 e or the 2 g depending on which M x or M y we choose to act diagonally. Thus, whether we draw the Wannier functions in 2 e or 2 g is a choice of basis (at least when periodic boundary conditions are chosen).

## Appendix D: Surface states of the DMQI

In this appendix, we analyze the edge spectrum for the DMQI on a generic boundary termination. To do this, we start by noting that the DMQI can be obtained from the C 4 -symmetric quadrupole model by breaking the C 4 symmetry. Since the C 4 -symmetric model is a proper Higher-order topological phase, we can employ the same techniques used to analyze HOTI surface states as in Refs. [23, 25, 29, 33]. We note that recently, an analysis of the surface states of the C 4 quadrupole model was employed in Ref. [35] by adding a symmetry-breaking mass term to the surface of a 2D TI. Our analysis differs in the fact that we do not need to explicitly employ spinful representation of the symmetries and we can work directly in terms of the projective representations defined in the main text.

The C 4 -symmetric model is obtained from Eq. 1 by setting γ x = γ y = γ and λ x = λ y = λ . In this case, the Hamiltonian has the extra C 4 symmetry implemented as

$$\begin{smallmatrix} g \\ n \text{-} \\ i \text{i} \text{-} \\ - \\ i \text{is} \end{smallmatrix} \begin{array} { c c c c } C _ { 4 z } \mathcal { H } ( k ) C _ { 4 z } ^ { \dagger } = \mathcal { H } ( O _ { 4 } k ), \quad O _ { 4 } ( k _ { x }, k _ { y } ) = ( k _ { y }, - k _ { x } ) \\ C _ { 4 z } = \Gamma _ { 2 } e ^ { \frac { \pi } { 4 } ( \Gamma _ { 2 } \Gamma _ { 4 } - \Gamma _ { 1 } \Gamma _ { 3 } ) } = - \sigma _ { 2 } \tau _ { 2 } e ^ { i \frac { \pi } { 2 } \frac { \sigma _ { 0 } - \sigma _ { 3 } } { 2 } \tau _ { 2 } } \quad ( D 1 ) \\ i \text{is} \quad \dots \, \dots \, \dots \, \dots \, \dots \, \end{array} \begin{array} { c } C _ { 4 z } \mathcal { H } ( O _ { 4 } k ), \quad O _ { 4 } ( k _ { x }, k _ { y } ) = ( k _ { y }, - k _ { x } ) \\ \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \end{array}$$

We will find it useful to define the modified Γ matrices Γ ± = 1 √ 2 (Γ 2 ± Γ ) 4 such that

$$C _ { 4 z } \Gamma _ { \pm } C _ { 4 z } ^ { \dagger } = \pm \Gamma _ { \pm } \text{ \quad \ \ } \text{(D2)}$$

The Hamiltonian for the C 4 symmetric model can then be written as

$$\text{pic} ^ { - } _ { \text{pre} ^ { - } } \quad \mathcal { H } _ { C _ { 4 } } = \lambda [ \sin k _ { x } \Gamma _ { 3 } + \sin k _ { y } \Gamma _ { 1 } ] + m _ { + } ( k ) \Gamma _ { + } + m _ { - } ( k ) \Gamma _ { - } \\ \text{ects}$$

$$m _ { + } ( \boldsymbol k ) = \frac { 1 } { \sqrt { 2 } } [ 2 \gamma + \lambda ( \cos k _ { x } + \cos k _ { y } ) ], \quad ( \mathbf D 4 )$$

$$m _ { - } ( \vec { k } ) = \frac { 1 } { \sqrt { 2 } } ( \cos k _ { y } - \cos k _ { x } ) \text{ \quad \ \ (D 5 ) }$$

where m ± ( O 4 k ) = ± m ± ( k ) .

The C 4 symmetric model is characterized by two topologically distinct bulk phases depending on whether | λ/γ | is smaller or larger than 1. These can be distinguished by the eigenvalues of C 4 z symmetry at the Γ and M points as follows. For | λ | &lt; | γ | , the C 4 z eigenvalues are given by { e 3 iπ/ 4 , e -3 iπ/ 4 } at both the Γ and M whereas for | λ | &gt; | γ | , the C 4 z eigenvalues are given by { e 3 iπ/ 4 , e -3 iπ/ 4 } at Γ and by { e iπ/ 4 , e -iπ/ 4 } at M . Thus, the two phases of the C 4 z symmetric quadrupole model are distinguished by a bulk phase transition at M . Close to such transition, we can write a continuum model by expanding the Hamiltonian close to the M point and going to real space as

$$\det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det\det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \dotdet \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det  \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \dot \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \delta \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \dt \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \cot \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \det \end{c}$$

The action of the spatial symmetries on m ± ( r ) can be deduced from the conditions (D2) and M x,y Γ ± M † x,y = Γ ± lead-

FIG. 21. Illustration of the edge theory of the DMQI on a genetic surface. Starting from the C 4 -symmetric model, we can derive the edge theory by taking the mass term m + to change sign from the inside to the outside of the sample. The edge theory at any point can be expressed in terms of the normal and tangent vectors n and t . The variable r t parametrizes the 1D edge which can be alternatively expressed in terms of an angular variable φ .

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00011v2-figure-4596b108-5db6-40cf-bdb0-716a0faa0057.jpg)

ing to

$$m _ { \pm } ( O _ { 4 } r ) = \pm m _ { \pm } ( O _ { 4 } r ), \quad & ( \mathbf D 7 ) \\ r _ { \ } u \rangle = m _ { \ } t ( r _ { \ } u \rangle \quad & m _ { \ } t ( r _ { \ } - u \rangle = m _ { \ } t ( r _ { \ } u \rangle \quad & \mathbf w \mathbf.$$

$$m _ { \pm } ( - x, y ) = m _ { \pm } ( x, y ), \quad m _ { \pm } ( x, - y ) = m _ { \pm } ( x, y ) \quad \text{ where } \\ ( D 8 ) \quad \text{ proper}$$

The edge is implemented by choosing m + ( r ) to change sign from being negative inside the sample (corresponding to the non-trivial phase) to being positive outside (corresponding to the trivial phase). Let us denote the normal to the edge at point r by n r ( ) = (cos φ ( r ) , sin φ ( r )) with the corresponding tangent t r ( ) = ( -sin φ ( r ) , cos φ ( r )) (cf. Fig. 21). We note that, in general, the normal vector depends on the position on the edge. In the following, we will drop this r dependence to simplify the notation. We also introduce

$$\Gamma _ { n } = \cos \phi \Gamma _ { 3 } + \sin \phi \Gamma _ { 1 }, \quad \Gamma _ { t } = - \sin \phi \Gamma _ { 3 } + \cos \phi \Gamma _ { 1 } \quad \text{ trivial} \\ ( D 9 ) \quad \text{ gapp}$$

$$r _ { n } = r \cdot n, \ \ r _ { t } = r \cdot t \quad \quad ( D 1 0 ) \quad \text{all} \quad \text{at $t$}$$

which enables us to write

$$\mathcal { H } _ { C _ { 4 } } = i \lambda ( \Gamma _ { n } \partial _ { r _ { n } } + \Gamma _ { t } \partial _ { r _ { t } } ) + m _ { + } ( r ) \Gamma _ { + } + m _ { - } ( r ) \Gamma _ { - } \ ( D 1 1 ) \quad \text{states.}$$

- [1] Charles L Kane and Eugene J Mele, 'Quantum spin Hall effect in graphene,' Phys. Rev. Lett. 95 , 226801 (2005).
- [2] Charles L Kane and Eugene J Mele, ' Z 2 topological order and the quantum spin Hall effect,' Phys. Rev. Lett. 95 , 146802 (2005).
- [3] B. Andrei Bernevig, Taylor L. Hughes, and Shou-Cheng Zhang, 'Quantum spin hall effect and topological phase transition in hgte quantum wells,' Science 314 , 1757-1761 (2006), http://science.sciencemag.org/content/314/5806/1757.full.pdf.
- [4] Xiao-Liang Qi, Taylor L. Hughes, and Shou-Cheng Zhang, 'Topological field theory of time-reversal invariant insulators,' Phys. Rev. B 78 , 195424 (2008).

To obtain the edge theory, we will assume that the mass m -( r ) is much smaller than the value of m + ( r ) in the bulk. If we also consider wavefunctions which vary slowly along the tangent direction to the edge, then we find that the low energy states are spanned by eigenfunctions of the form

$$\psi ( r ) = \frac { 1 } { 2 } ( 1 - i \Gamma _ { n } \Gamma _ { + } ) e ^ { - \frac { 1 } { x } \int ^ { r } d r _ { n } ^ { \prime } m _ { + } ( r _ { n } ^ { \prime }, r _ { t } ) } \chi ( r _ { t } ) \pmod { 0 }$$

which are annihilated by iλ Γ n ∂ r n + m + ( r )Γ + for any function χ r ( t ) . Thus, the edge Hamiltonian is given by

$$\mathcal { H } _ { \text{edge} } = i \lambda \gamma _ { t } \partial _ { r _ { t } } + m _ { - } ( r _ { t } ) \gamma _ { - } \quad \text{(D13)}$$

where γ , γ t -are obtained by projecting Γ t , Γ -onto the space spanned by ψ ( r ) . The edge Hamiltonian is a one-dimensional Dirac Hamiltonian defined in terms of the variable r t along the direction parallel to the edge. Changing the variable from r t to the angular variable φ (cf. Fig. 21), and using (D7), we find that

$$m _ { - } ( \phi + \pi / 2 ) = - m _ { - } ( \phi ) \quad & \text{(D14)} \\ - \phi \partial = m _ { - } ( \phi \partial \quad m _ { - } ( \pi - \phi \partial = m _ { - } ( \phi \partial \quad \theta ) \quad \theta \text{b} \text{.}$$

$$m _ { - } ( - \phi ) = m _ { - } ( \phi ), \ \ m _ { - } ( \pi - \phi ) = m _ { - } ( \phi ) \quad ( \mathbf D 1 5 )$$

where the first condition follows from the transformation properties under C 4 rotation (D7) and the second follows from M x and M y (D8). Eq. D14 implies that the edge Dirac Hamiltonian hosts four domain walls where the mass term m -changes sign.

glyph[negationslash]

Notice now that if we break C 4 symmetry, these four domain walls can be freely moved to annihilate. However, Eq. D15 would dictate that they can only annihilate at φ = 0 or π , i.e. at a high symmetry surface. If we enforce the condition that m -(0) , m -( π ) = 0 , we can never get rid of these domain walls. We notice also that Eq. D15 is also satisfied by the trivial phase where m -( φ ) = const . leading to a completely gapped edge. Thus, For the C 4 -broken model, the symmetries allow for two possibilities with or without zero energy states at the edge. In contrast, for the C 4 symmetric model which is a HOTI, the symmetry forces the mass term to change sign enforcing the existence of four symmetry-related zero energy states.

- [5] Joel Moore, 'Topological insulators: The next generation,' Nat. Phys. 5 , 378-380 (2009).
- [6] MZahid Hasan and Charles L Kane, 'Colloquium: topological insulators,' Rev. Mod. Phys. 82 , 3045 (2010).
- [7] Xiao-Liang Qi and Shou-Cheng Zhang, 'Topological insulators and superconductors,' Rev. Mod. Phys. 83 , 1057 (2011).
- [8] Marcel Franz and Laurens Molenkamp, eds., Topological Insulators , Contemporary Concepts of Condensed Matter Science, Vol. 6 (Elsevier, 2013).
- [9] Alexei Kitaev, 'Periodic table for topological insulators and superconductors,' in AIP Conference Proceedings , Vol. 1134 (AIP, 2009) pp. 22-30.

- [10] Andreas P. Schnyder, Shinsei Ryu, Akira Furusaki, and Andreas W. W. Ludwig, 'Classification of topological insulators and superconductors,' in AIP Conference Proceedings , Vol. 1134 (AIP, 2009) pp. 10-21.
- [11] Shinsei Ryu, Andreas P. Schnyder, Akira Furusaki, and Andreas W. W. Ludwig, 'Topological insulators and superconductors: tenfold way and dimensional hierarchy,' New J. Phys. 12 , 065010 (2010).
- [12] A. Alexandradinata, Zhijun Wang, and B. Andrei Bernevig, 'Topological insulators from group cohomology,' Phys. Rev. X 6 , 021008 (2016).
- [13] Luka Trifunovic and Piet W. Brouwer, 'Higher-order bulkboundary correspondence for topological crystalline phases,' Phys. Rev. X 9 , 011012 (2019).
- [14] P Dziawa, BJ Kowalski, K Dybko, R Buczko, A Szczerbakow, MSzot, E Łusakowska, T Balasubramanian, Bastian M Wojek, MHBerntsen, et al. , 'Topological crystalline insulator states in Pb 1 -x Sn Se,' Nature materials x 11 , 1023-1027 (2012).
- [15] Y Tanaka, Zhi Ren, T Sato, K Nakayama, S Souma, T Takahashi, Kouji Segawa, and Yoichi Ando, 'Experimental realization of a topological crystalline insulator in SnTe,' Nature Physics 8 , 800-803 (2012).
- [16] Wladimir A. Benalcazar, Jeffrey C. Y. Teo, and Taylor L. Hughes, 'Classification of two-dimensional topological crystalline superconductors and majorana bound states at disclinations,' Phys. Rev. B 89 , 224503 (2014).
- [17] Wladimir A Benalcazar, B. Andrei Bernevig, and Taylor L Hughes, 'Quantized electric multipole insulators,' Science 66 , 61-66 (2017).
- [18] Wladimir A Benalcazar, B Andrei Bernevig, and Taylor L Hughes, 'Electric multipole moments, topological multipole moment pumping, and chiral hinge states in crystalline insulators,' Phys. Rev. B 96 , 245115 (2017).
- [19] Josias Langbehn, Yang Peng, Luka Trifunovic, Felix von Oppen, and Piet W. Brouwer, 'Reflection-symmetric second-order topological insulators and superconductors,' Phys. Rev. Lett. 119 , 246401 (2017).
- [20] Frank Schindler, Ashley M. Cook, Maia G. Vergniory, Zhijun Wang, Stuart S. P. Parkin, B. Andrei Bernevig, and Titus Neupert, 'Higher-order topological insulators,' Science Advances 4 , eaat0346 (2018).
- [21] Zhida Song, Zhong Fang, and Chen Fang, ' ( d -2) -dimensional edge states of rotation symmetry protected topological states,' Phys. Rev. Lett. 119 , 246402 (2017).
- [22] Chen Fang and Liang Fu, 'Rotation anomaly and topological crystalline insulators,' arXiv preprint arXiv:1709.01929 (2017).
- [23] Eslam Khalaf, Hoi Chun Po, Ashvin Vishwanath, and Haruki Watanabe, 'Symmetry indicators and anomalous surface states of topological crystalline insulators,' Phys. Rev. X 8 , 031070 (2018).
- [24] Max Geier, Luka Trifunovic, Max Hoskam, and Piet W. Brouwer, 'Second-order topological insulators and superconductors with an order-two crystalline symmetry,' Phys. Rev. B 97 , 205135 (2018).
- [25] Eslam Khalaf, 'Higher-order topological insulators and superconductors protected by inversion symmetry,' Phys. Rev. B 97 , 205136 (2018).
- [26] Wladimir A. Benalcazar, Tianhe Li, and Taylor L. Hughes, 'Quantization of fractional corner charge in C n -symmetric higher-order topological crystalline insulators,' Phys. Rev. B 99 , 245151 (2019).
- [27] Christopher W. Peterson, Tianhe Li, Wladimir A. Benalcazar, Taylor L. Hughes, and Gaurav Bahl, 'A fractional corner
- anomaly reveals higher-order topology,' Science 368 , 11141118 (2020).
- [28] Jeffrey C.Y. Teo and Taylor L Hughes, 'Topological defects in symmetry-protected topological phases,' Annual Review of Condensed Matter Physics 8 , 211-237 (2017).
- [29] Shang Liu, Ashvin Vishwanath, and Eslam Khalaf, 'Shift insulators: Rotation-protected two-dimensional topological crystalline insulators,' Phys. Rev. X 9 , 031003 (2019).
- [30] Tianhe Li, Penghao Zhu, Wladimir A. Benalcazar, and Taylor L. Hughes, 'Fractional disclination charge in twodimensional c n -symmetric topological crystalline insulators,' arxiv:1906.02752 (2019).
- [31] Hoi Chun Po, Ashvin Vishwanath, and Haruki Watanabe, 'Symmetry-based indicators of band topology in the 230 space groups,' Nature Communications 8 , 50 (2017).
- [32] Barry Bradlyn, L Elcoro, Jennifer Cano, MG Vergniory, Zhijun Wang, C Felser, MI Aroyo, and B Andrei Bernevig, 'Topological quantum chemistry,' Nature 547 , 298 (2017).
- [33] Yoonseok Hwang, Junyeong Ahn, and Bohm-Jung Yang, 'Fragile topology protected by inversion symmetry: Diagnosis, bulk-boundary correspondence, and wilson loop,' Phys. Rev. B 100 , 205126 (2019).
- [34] Oleg Dubinkin and Taylor L Hughes, 'Entanglement signatures of multipolar higher order topological phases,' arXiv preprint arXiv:2002.08385 (2020).
- [35] Benjamin J Wieder, Zhijun Wang, Jennifer Cano, Xi Dai, Leslie M Schoop, Barry Bradlyn, and B Andrei Bernevig, 'Strong and fragile topological Dirac semimetals with higherorder Fermi arcs,' Nature Communications 11 , 627 (2020).
- [36] Nicola Marzari and David Vanderbilt, 'Maximally localized generalized wannier functions for composite energy bands,' Phys. Rev. B 56 , 12847-12865 (1997).
- [37] Nicola Marzari, Arash A. Mostofi, Jonathan R. Yates, Ivo Souza, and David Vanderbilt, 'Maximally localized wannier functions: Theory and applications,' Rev. Mod. Phys. 84 , 1419-1475 (2012).
- [38] Rui Yu, Xiao Liang Qi, Andrei Bernevig, Zhong Fang, and Xi Dai, 'Equivalent expression of Ϝ 2 topological invariant for band insulators using the non-abelian berry connection,' Phys. Rev. B 84 , 075119 (2011).
- [39] Yizhi You, FJ Burnell, and Taylor L Hughes, 'Multipolar topological field theories: Bridging higher order topological insulators and fractons,' arXiv preprint arXiv:1909.05868 (2019).
- [40] J. Zak, 'Symmetry specification of bands in solids,' Phys. Rev. Lett. 45 , 1025-1028 (1980).
- [41] J Zak, 'Band representations of space groups,' Physical Review B 26 , 3010 (1982).
- [42] L Michel and J Zak, 'Elementary energy bands in crystalline solids,' EPL (Europhysics Letters) 50 , 519 (2000).
- [43] MGVergniory, L Elcoro, Zhijun Wang, Jennifer Cano, C Felser, MI Aroyo, B Andrei Bernevig, and Barry Bradlyn, 'Graph theory data for topological quantum chemistry,' Physical Review E 96 , 023310 (2017).
- [44] 'Wikipedia entry for wallpaper group,' https://en. wikipedia.org/wiki/Wallpaper\_group .
- [45] Jennifer Cano, Barry Bradlyn, Zhijun Wang, L Elcoro, MG Vergniory, C Felser, MI Aroyo, and B Andrei Bernevig, 'Building blocks of topological quantum chemistry: Elementary band representations,' Physical Review B 97 , 035139 (2018).
- [46] Lukasz Fidkowski, T. S. Jackson, and Israel Klich, 'Model characterization of gapless edge modes of topological insulators using intermediate Brillouin-zone functions,' Physical Review Letters 107 , 1-4 (2011).

- [47] Yan-Bin Yang, Kai Li, L.-M. Duan, and Yong Xu, 'Typeii quadrupole topological insulators,' Phys. Rev. Research 2 , 033029 (2020).
- [48] Jeffrey Goldstone and Frank Wilczek, 'Fractional quantum numbers on solitons,' Phys. Rev. Lett. 47 , 986-989 (1981).
- [49] Benjamin J Wieder and B Andrei Bernevig, 'The axion insulator as a pump of fragile topology,' arXiv preprint arXiv:1810.02373 (2018).
- [50] Guido van Miert and Carmine Ortix, 'Higher-order topological insulators protected by inversion and rotoinversion symmetries,' Phys. Rev. B 98 , 081110 (2018).
- [51] Jun-Won Rhim, Jan Behrends, and Jens H. Bardarson, 'Bulkboundary correspondence from the intercellular zak phase,' Phys. Rev. B 95 , 035421 (2017).
- [52] Hao Song, Sheng-Jie Huang, Liang Fu, and Michael Hermele, 'Topological phases protected by point group symmetry,' Physical Review X 7 , 011020 (2017).
- [53] Sheng-Jie Huang, Hao Song, Yi-Ping Huang, and Michael Hermele, 'Building crystalline topological phases from lowerdimensional states,' Phys. Rev. B 96 , 205106 (2017).
- [54] Zhida Song, Tiantian Zhang, Zhong Fang, and Chen Fang, 'Quantitative mappings between symmetry and topology in solids,' Nature Communications 9 , 3530 (2018).
- [55] Barry Bradlyn, L. Elcoro, Jennifer Cano, M. G. Vergniory, Zhijun Wang, C. Felser, M. I. Aroyo, and B. Andrei Bernevig, 'Topological quantum chemistry,' Nature 547 , 298 EP -(2017), article.
- [56] Christian Gross and Immanuel Bloch, 'Quantum simulations with ultracold atoms in optical lattices,' Science 357 , 995-1001 (2017), https://science.sciencemag.org/content/357/6355/995.full.pdf.
- [57] Jeffrey C. Y. Teo and C. L. Kane, 'Topological defects and gapless modes in insulators and superconductors,' Phys. Rev. B 82 , 115120 (2010).
- [58] Hoi Chun Po, Haruki Watanabe, and Ashvin Vishwanath, 'Fragile topology and wannier obstructions,' Phys. Rev. Lett. 121 , 126402 (2018).
- [59] Takahiro Fukui, 'Dirac fermion model associated with a second-order topological insulator,' Physical Review B 99 , 165129 (2019).
- [60] Zhida Song, Tiantian Zhang, and Chen Fang, 'Diagnosis for nonmagnetic topological semimetals in the absence of spinorbital coupling,' Phys. Rev. X 8 , 031069 (2018).
- [61] Taylor L Hughes, Emil Prodan, and B Andrei Bernevig, 'Inversion-symmetric topological insulators,' Phys. Rev. B 83 , 245132 (2011).
- [62] Ari M Turner, Yi Zhang, Roger SK Mong, and Ashvin Vishwanath, 'Quantized response and topology of magnetic insulators with inversion symmetry,' Phys. Rev. B 85 , 165120 (2012).
- [63] John Baez, 'Torsors made easy,' (2009).
- [64] Barry Bradlyn, L. Elcoro, Jennifer Cano, M. G. Vergniory, Zhijun Wang, C. Felser, Mois I. Aroyo, and B. Andrei Bernevig, 'Topological quantum chemistry,' Nature 547 , 298305 (2017), arXiv:NIHMS150003.
- [65] Jennifer Cano, Barry Bradlyn, Zhijun Wang, L. Elcoro, M. G. Vergniory, C. Felser, Mois I. Aroyo, and B. Andrei Bernevig, 'Building blocks of topological quantum chemistry: Elementary band representations,' Physical Review B 97 , 035139 (2018), arXiv:1709.01935.
- [66] Robert A Evarestov, Quantum chemistry of solids: the LCAO first principles treatment of crystals , Vol. 153 (Springer Science &amp;Business Media, 2007).
- [67] Benjamin J. Wieder, Barry Bradlyn, Zhijun Wang, Jennifer Cano, Youngkuk Kim, Hyeong-Seok D. Kim,
- Andrew M. Rappe, C. L. Kane, and B. Andrei Bernevig, 'Wallpaper fermions and the nonsymmorphic dirac insulator,' Science 361 , 246-251 (2018), https://science.sciencemag.org/content/361/6399/246.full.pdf.
- [68] Heejae Kim, Ken Shiozaki, and Shuichi Murakami, 'Glidesymmetric magnetic topological crystalline insulators with inversion symmetry,' Phys. Rev. B 100 , 165202 (2019).
- [69] Benjamin J. Wieder and C. L. Kane, 'Spin-orbit semimetals in the layer groups,' Phys. Rev. B 94 , 155108 (2016).
- [70] Mois Ilia Aroyo, Juan Manuel Perez-Mato, Cesar Capillas, Eli Kroumova, Svetoslav Ivantchev, Gotzon Madariaga, Asen Kirov, and Hans Wondratschek, 'Bilbao crystallographic server: I. databases and crystallographic computing programs,' Zeitschrift für Kristallographie-Crystalline Materials 221 , 1527 (2006).
- [71] A. P. Crackwell and K. C. Wong, 'Double-valued corepresentations of magnetic point groups,' Australian Journal of Physics 20 , 173 (1967).
- [72] Nicodemos Varnava and David Vanderbilt, 'Surfaces of axion insulators,' Physical Review B 98 , 245117 (2018).
- [73] Liang Fu and C L Kane, 'Topological insulators with inversion symmetry,' 10.1103/PhysRevB.76.045302.
- [74] Liang Fu and Charles L Kane, 'Time reversal polarization and a z 2 adiabatic spin pump,' Physical Review B 74 , 195312 (2006).
- [75] Jiho Noh, Wladimir A. Benalcazar, Sheng Huang, Matthew J. Collins, Kevin P. Chen, Taylor L. Hughes, and Mikael C. Rechtsman, 'Topological protection of photonic mid-gap defect modes,' Nature Photonics (2018), 10.1038/s41566-018-01793.
- [76] Christopher W. Peterson, Wladimir A. Benalcazar, Taylor L. Hughes, and Gaurav Bahl, 'A quantized microwave quadrupole insulator with topologically protected corner states,' Nature 555 , 346 EP - (2018).
- [77] Marc Serra-Garcia, Roman Süsstrunk, and Sebastian D. Huber, 'Observation of quadrupole transitions and edge mode topology in an lc circuit network,' Phys. Rev. B 99 , 020304 (2019).
- [78] Stefan Imhof, Christian Berger, Florian Bayer, Johannes Brehm, Laurens W. Molenkamp, Tobias Kiessling, Frank Schindler, Ching Hua Lee, Martin Greiter, Titus Neupert, and Ronny Thomale, 'Topolectrical-circuit realization of topological corner modes,' Nature Physics 14 , 925-929 (2018).
- [79] Bi Ye Xie, Hong Fei Wang, Hai-Xiao Wang, Xue Yi Zhu, Jian-Hua Jiang, Ming Hui Lu, and Yan Feng Chen, 'Secondorder photonic topological insulator with corner states,' arxiv:1805.07555 (2018).
- [80] Xiang Ni, Matthew Weiner, Andrea Alu, and Alexander B. Khanikaev, 'Observation of bulk polarization transitions and higher-order embedded topological eigenstates for sound,' arxiv:1807.00896 (2018).
- [81] Haoran Xue, Yahui Yang, Fei Gao, Yidong Chong, and Baile Zhang, 'Acoustic higher-order topological insulator on a kagome lattice,' arXiv:1806.09418 (2018).
- [82] S. N. Kempkes, M. R. Slot, J. J. van den Broeke, P. Capiod, W. A. Benalcazar, D. Vanmaekelbergh, D. Bercioux, I. Swart, and C. Morais Smith, 'Robust zero-energy modes in an electronic higher-order topological insulator,' Nature Materials 18 , 1292-1297 (2019).
- [83] Haoran Xue, Yahui Yang, Guigeng Liu, Fei Gao, Yidong Chong, and Baile Zhang, 'Realization of an acoustic thirdorder topological insulator,' Phys. Rev. Lett. 122 , 244301 (2019).
- [84] Sunil Mittal, Venkata Vikram Orre, Guanyu Zhu, Maxim A. Gorlach, Alexander Poddubny, and Mohammad Hafezi, 'Pho-

- tonic quadrupole topological phases,' Nature Photonics 13 , 692-696 (2019).
- [85] Jiacheng Bao, Deyuan Zou, Weixuan Zhang, Wenjing He, Houjun Sun, and Xiangdong Zhang, 'Topoelectrical circuit octupole insulator with topologically protected corner states,' Phys. Rev. B 100 , 201406 (2019).
- [86] Haoran Xue, Yong Ge, Hong-Xiang Sun, Qiang Wang, Ding Jia, Yi-Jun Guan, Shou-Qi Yuan, Yidong Chong, and Baile Zhang, 'Quantized octupole acoustic topological insulator,' (2019), arXiv:1911.06068 [cond-mat.mes-hall].
- [87] Xiang Ni, Mengyao Li, Matthew Weiner, Andrea Alù, and Alexander B. Khanikaev, 'Demonstration of a quantized acoustic octupole topological insulator,' (2019), arXiv:1911.06469 [cond-mat.mes-hall].
- [88] M. Aidelsburger, M. Atala, S. Nascimbène, S. Trotzky, Y.-A. Chen, and I. Bloch, 'Experimental realization of strong effective magnetic fields in an optical lattice,' Phys. Rev. Lett. 107 , 255301 (2011).
- [89] Xiaopeng Li, Erhai Zhao, and W. Vincent Liu, 'Topological states in a ladder-like optical lattice containing ultracold atoms in higher orbital bands,' Nature Communications 4 , 1523 (2013).
- [90] Marcos Atala, Monika Aidelsburger, Julio T. Barreiro, Dmitry Abanin, Takuya Kitagawa, Eugene Demler, and Immanuel Bloch, 'Direct measurement of the zak phase in topological bloch bands,' Nature Physics 9 , 795-800 (2013).
- [91] Marcos Atala, Monika Aidelsburger, Michael Lohse, Julio T. Barreiro, Belén Paredes, and Immanuel Bloch, 'Observation of chiral currents with ultracold atoms in bosonic ladders,' Nature Physics 10 , 588-593 (2014).
- [92] M. Aidelsburger, M. Lohse, C. Schweizer, M. Atala, J. T. Barreiro, S. Nascimbène, N. R. Cooper, I. Bloch, and N. Goldman, 'Measuring the chern number of hofstadter bands with ultracold bosonic atoms,' Nature Physics 11 , 162-166 (2015).
- [93] M. Mancini, G. Pagano, G. Cappellini, L. Livi, M. Rider, J. Catani, C. Sias, P. Zoller, M. Inguscio, M. Dalmonte, and L. Fallani, 'Observation of chiral edge states with neutral fermions in synthetic hall ribbons,' Science 349 , 1510-1513 (2015), https://science.sciencemag.org/content/349/6255/1510.full.pdf.
- [94] Tracy Li, Lucia Duca, Martin Reitter, Fabian Grusdt, Eugene Demler, Manuel Endres, Monika Schleier-Smith, Immanuel Bloch, and Ulrich Schneider, 'Bloch state tomography using wilson lines,' Science 352 , 1094-1097 (2016), https://science.sciencemag.org/content/352/6289/1094.full.pdf.
- [95] Martin Leder, Christopher Grossert, Lukas Sitta, Maximilian Genske, Achim Rosch, and Martin Weitz, 'Real-space imaging of a topologically protected edge state with ultracold atoms in an amplitude-chirped optical lattice,' Nature Communications 7 , 13112 (2016).
- [96] Luka Trifunovic and Piet W Brouwer, Higher-order topological band structures , Tech. Rep., arXiv:2003.01144v1.
Eslam Khalaf Wladimir A. Benalcazar Taylor L. Hughes Raquel Queiroz
Boundary-obstructed topological phases
Symmetry protected topological (SPT) phases are gapped phases of matter that
cannot be deformed to a trivial phase without breaking the symmetry or closing
the bulk gap. Here, we introduce a new notion of a topological obstruction that
is not captured by bulk energy gap closings in periodic boundary conditions.
More specifically, we say two bulk Hamiltonians belong to distinct boundary
obstructed topological `phases' (BOTPs) if they can be deformed to each other
on a system with periodic boundaries, but cannot be deformed to each other for
symmetric open boundaries without closing the gap at at least one high symmetry
region on the surface. BOTPs are not topological phases of matter in the
standard sense since they are adiabatically deformable to each other on a torus
but, similar to SPTs, they are associated with boundary signatures in open
systems such as surface states or fractional corner or hinge charges. We show
that the double-mirror quadrupole model of [Science, 357(6346), 2018] is a
prototypical example of such phases, and present a detailed analysis of several
aspects of boundary obstructions in this model. In addition, we introduce
several three-dimensional models having boundary obstructions, which are
characterized either by surface states or fractional corner or hinge charges.
We also provide a general framework to study boundary obstructions in
free-fermion systems in terms of Wannier band representations (WBR), an
extension of the recently-developed band representation formalism to Wannier
bands. WBRs capture the notion of topological obstructions in the Wannier bands
which can then be used to study topological obstructions in the boundary
spectrum by means of the correspondence between the Wannier and boundary
spectra. This establishes a form of bulk-boundary correspondence for BOTPs by
relating the bulk band representation to the boundary topology.