Abstract
We present a scheme to explicitly construct and classify general topological states jointly protected by an onsite symmetry group and a spatial symmetry group. We show that all these symmetry-protected topological states can be adiabatically deformed into a special class of states we call topological crystals. A topological crystal in, for example, three dimensions is a real-space assembly of finite-sized pieces of topological states in one and two dimensions protected by the local symmetry group alone, arranged in a configuration invariant under the spatial group and glued together such that there is no open edge or end. As a demonstration of principle, we explicitly enumerate all inequivalent topological crystals for noninteracting time-reversal symmetric electronic insulators with spin-orbit coupling and any one of the 230 space groups. This enumeration gives topological crystalline insulators a full classification.
INTRODUCTION
Symmetry-protected topological (SPT) phases are gapped many-body ground states that can only be adiabatically deformed into product states of local orbitals by breaking a given symmetry group or by closing the energy gap (1). Typical examples are topological insulators, topological superconductors, and the Haldane spin chain [hereinafter, if not explicitly stated otherwise, the terms “topological insulator” and “strong topological insulator” (STI) will refer to the three-dimensional (3D) topological insulator protected by time-reversal symmetry]. The best-understood SPT phases are those of noninteracting fermions; not long after the discovery of topological band insulators, free-fermion topological phases were completely classified for systems with internal (i.e., nonspatial) symmetries (2–5). Of course, crystalline symmetries play a central role in solid-state physics, so attention naturally began to turn to topological crystalline insulators (TCIs), which are electronic insulators whose topologically nontrivial nature is protected, in part, by point group or space group symmetry (5–8). Sparked by the prediction and observation of TCIs in SnTe (7, 9–11), remarkable theoretical (12–24) and experimental progress (25–27) has followed over the last few years.
Despite these developments, somewhat unexpectedly, a unified picture of the classification of noninteracting electron TCIs has yet to emerge. The primary tool for the classification of free-fermion topological phases with spatial symmetries has been K-theory (3) and equivariant K-theory (28). A number of concrete classification results have been obtained (12, 14, 16, 17, 22), but reflecting the complexity of K-theory, there is a paucity of concrete results for 3D (d = 3) insulators with general space group symmetry and time-reversal symmetry. Moreover, it is not understood how or whether electron interactions can be included within K-theory. Therefore, there is a need to develop alternate means to classify TCIs and other crystalline SPT (cSPT) phases jointly protected by internal and spatial symmetries. Ideally, to provide a useful complement to K-theory, these methods will be real space based and physically transparent and allow for interactions to be included.
Here, we propose a general method for classifying and constructing cSPT phases, which is then applied to the case of electronic TCIs in all 230 space groups, with time-reversal symmetry and spin-orbit coupling. We also extend these results to include STIs. Our approach is based on recent developments in the seemingly harder problem of classifying interacting cSPT phases (29–33). The key idea is to first argue that any cSPT phase is adiabatically connected to a real-space crystalline pattern of lower-dimensional topological states, which we refer to as a topological crystal (29, 32). One then develops a classification of phases of matter in terms of topological crystals. For bosonic cSPT phases with only space group symmetry, the resulting classification (32) agrees with that obtained in complementary approaches based on tensor network states and gauging crystalline symmetry (30, 31). Recently, Shiozaki et al. (34) have discussed how to formulate the topological crystal approach within K-theory via the Atiyah-Hirzebruch spectral sequence.
Our approach is related to, but goes beyond, layer constructions of TCIs. Any construction of a TCI in terms of decoupled layers, including the archetype of weak topological insulators as stacks of d = 2 topological insulators, is a topological crystal. However, by comparison to the recent systematic study of layer constructions in (35), we show that in certain nonsymmorphic space groups, there are TCIs that do not have a layer construction but can still be realized as topological crystals. We emphasize that all the topological states implied by symmetry eigenvalues proposed in (19, 20) are contained in our classification. The work in (35) shows that only five symmetry-eigenvalue-implied TCIs—five weak topological insulators—are not layer constructable. Here, we explicitly construct these weak topological insulators as topological crystals.
The results we obtain are related to the recent work of Khalaf et al. (36), who considered TCIs with anomalous surface states (dubbed sTCIs) and proposed a classification for sTCIs with point group and space group symmetry via the surface states of doubled STIs. The TCI classifications produced by our approach, which focuses on the bulk and does not assume anomalous surface states or a description in terms of Dirac fermions, agree with the sTCI classifications in (36). This agreement shows that all the TCIs we classify have anomalous surface states for some surface termination.
RESULTS
Topological crystals
We begin by considering a d = 3 system with symmetry G = Gint × Gc, where Gint is some internal symmetry and Gc is either a crystalline site symmetry (i.e., point group) or space group. (The assumption that G is a direct product of Gint and Gc is not necessary and is only made for simplicity of discussion and because it holds for the electronic TCIs to be later discussed.) We assume the system is in an SPT phase (which could be the trivial phase); that is, below an energy gap, the ground state ∣ψ> is unique and symmetry preserving, and moreover, ∣ψ> is adiabatically connected to a trivial product state if the symmetry G is broken explicitly. Moreover, we restrict to those SPT phases that only remain nontrivial in the presence of crystalline symmetry; that is, ∣ψ> is adiabatically connected to a product state if Gc is explicitly broken, even if Gint is preserved. To avoid complications associated with gapless boundary states, we consider periodic boundary conditions.
To proceed, we identify an asymmetric unit (AU), which is the interior of a region of space that is as large as possible, subject to the condition that no two points in the region are related by a crystalline symmetry. The AU is then copied throughout space using the crystalline symmetry, and we denote the resulting union of nonoverlapping AUs as 𝒜. This construction gives 3D space a cell complex structure (see Methods), where the 3-cells are the individual (nonoverlapping) copies of the AU in 𝒜. The 2-cells lie on faces where two 3-cells meet, with the property that no two distinct points in the same 2-cells are related by symmetry. Similarly, 1-cells are edges where two or more faces meet, and 0-cells are points where edges meet. The 3-cells are in one-to-one correspondence with elements of Gc: arbitrarily choosing one 3-cell to correspond to the identity element, each other 3-cell is the unique image of this one upon acting with g ∈ Gc. The 2-skeleton X2 is the complement of 𝒜. An example of this cell structure for the space group
(A) The AU 0 < x, y < 1 and
The work in (32) argued that ∣ψ> is adiabatically connected to a product of a trivial state on 𝒜, with a possibly nontrivial state on X2 [see section VI of (32)]. More precisely, one considers a thickened version of X2, with characteristic thickness w, and its complement. For the argument to go through, it is important that w ≫ ξ, where ξ is any characteristic correlation or entanglement length of the short-range entangled state ∣ψ>. If Gc is a point group symmetry, then this requires no assumptions because w can be taken sufficiently large. For space group symmetry, w is limited by the unit cell size, and one must make the assumption that, by adding a fine mesh of trivial degrees of freedom, it is possible to make the correlation length ξ as small as desired. This assumption not only allows the reduction to a topological crystal state but also implies that the correlation length of the topological crystal state itself is much smaller than the unit cell size; this is important because it allows us to associate a well-defined lower-dimensional state with each cell of X2. While we believe that this assumption is likely to hold, it is not proven, and strictly speaking, it should be treated as a conjecture. If this conjecture is false, then our approach is simply restricted to those cSPT phases whose correlation length is not bounded below upon adding trivial degrees of freedom. We note that a preliminary version of the idea of reduction to X2 was discussed in (32); there, unlike in the present work, the idea was not developed into a tool to obtain classifications.
The result of this reduction procedure is a topological crystal state. The state on X2 can be understood by associating a db-dimensional topological phase with each db-cell of X2, where db = 0, 1, and 2. These lower-dimensional states are referred to as the “building blocks” of the topological crystal, and db is the block dimension. The building blocks must be glued together so as to eliminate any gapless modes in the bulk while preserving symmetry; for instance, db = 2 blocks will generally have gapless edge modes, which must gap out at the 1-cells and 0-cells where the building blocks meet. Whereas crystals are formed by periodically arranged atoms, i.e., zero-dimensional objects, topological crystals are “stacked” from building blocks which are themselves topological states in lower dimensions.
TCI classification
First, we note that a number of topological invariants are already known that distinguish different phases and should be a part of any classification of TCIs. In particular, these invariants include the weak ℤ2 invariants (37, 38), mirror Chern numbers (ℤ invariant) (39), and ℤ2 invariants associated with rotation (40), glide reflection (17, 18), inversion (21, 41, 42), rotoreflection (35, 36), and screw rotation (35, 36) symmetries. All possible combinations of these invariants that can be realized in TCIs with a layer construction were enumerated in (35). While one could not prove a priori that these seven quantities exhaust all independent topological invariants, in this work, we show they play a special role as a complete list of invariants for TCIs. That is, we find that any two inequivalent TCIs differ by at least one of these invariants.
To further apply the tool of topological crystals to the case of electronic TCIs, we consider a system of noninteracting electrons with spin-orbit coupling, with internal symmetries of charge conservation and time reversal; that is,
The next step is to understand what kind of topological crystal states can be placed on X2. First, we consider topological crystals built out of d = 2 topological states. There are two kinds of 2-cells, those that coincide with a mirror plane and those that do not. 2-Cells coinciding with a mirror plane can host a d = 2 mirror Chern insulator (MCI) state, which is characterized by a ℤ invariant (39). The MCI state can be understood by diagonalizing the mirror operation σ : z → −z, where z is the coordinate along the normal direction of the mirror plane. Because σ2 = (−1)F, one-electron wave functions can be divided into two sectors with mirror eigenvalue ±i. Because σT = Tσ, time reversal exchanges these two sectors, which therefore have equal and opposite Chern numbers, leading to a ℤ invariant. Each sector can be understood as a d = 2 fermion system in class A, which has a ℤ classification. We see that the relevant symmetry class is, thus, effectively modified from AII to A on a mirror plane; this modification of symmetry class is familiar from classifications of reflection-symmetric free-fermion topological phases in momentum space (12). For 2-cells not coinciding with a mirror plane, the symmetry class remains AII, and these cells can host a d = 2 STI (2dTI), which is characterized by a ℤ2 invariant.
We also have to consider the possibility of topological crystals built from the d = 1 and d = 0 states. 1-Cells only host trivial states: The effective symmetry class on a 1-cell can be either AII or A (see the discussion on MCIs above), and in either case, the classification in d = 1 is trivial. On the other hand, there are nontrivial topological crystals built from db = 0 building blocks, which are atomic insulators formed from patterns of localized filled orbitals. Although there are distinct atomic insulators constituting different quantum phases of matter and although these distinctions may be a source of interesting physics, all these states, being product states of localized orbitals, are, in a sense, topologically trivial. Therefore, we ignore distinctions among atomic insulators in our classification. Formally, this is accomplished by taking a certain quotient (Methods). In the future, our results could be extended to include db = 0 building blocks, which would facilitate a more direct comparison with K-theory–based approaches, which do include such states.
We, thus, see that there are two kinds of TCIs, both built from db = 2 blocks. We refer to TCIs built from MCI blocks as mirror TCIs (MTCIs), while TCIs built from 2dTI blocks are dubbed ℤ2 TCIs. Of course, a general TCI can have mixed MTCI and ℤ2 TCI character, and the classification is a product of MTCI and ℤ2 TCI classifications. To proceed, we consider the requirement that the building blocks must be glued together to eliminate any gapless modes in the bulk. As shown in Methods, this requirement implies that MTCIs can always be decomposed into decoupled planar MCI layers.
ℤ2 TCIs are not quite as simple. If we consider placing a 2dTI on some subset of the 2-cells of X2, then it can be shown these building blocks can be glued together into a topological crystal if and only if every 1-cell is the edge of an even number of 2dTI blocks (Methods). While, sometimes, these states can be decomposed into decoupled 2dTI layers, this is not always true. For example, in space group P42212 (#94), for which the possible topological crystals are described below, we find a topological crystal that is beyond the scope of layer construction, as shown in Fig. 2B. In this state, the 2-cells decorated with 2dTIs form a complicated yet connected structure. Intuitively, one may lower the two yellow facets at
(A) The symmetry-inequivalent 2-cells
Having described TCIs in terms of topological crystals, we next use these states to classify TCIs. First, we discuss equivalence relations among topological crystal states and argue that two distinct topological crystals on X2 are in different phases. Following (29, 32), we need to consider an additional equivalence relation, beyond those for the d = 2 phases of matter on the 2-cells. Within an AU and all its copies under symmetry, we create a small bubble of 2dTI and expand the bubble until it joins with the AU boundary; this process can be achieved by adiabatic evolution, preserving symmetry, within a finite time, so any two states related in this way belong to the same phase. (Equivalently, this process can be achieved by acting with a finite depth symmetry-preserving quantum circuit.). The reason we consider a bubble of 2dTI and not something else is that this is the only nontrivial d = 2 state that can exist within the AU, where the only symmetries are charge conservation and time reversal. The source of this equivalence relation is the arbitrary width w of the thickened X2 space in the dimensional reduction procedure; making w larger corresponds to bringing in additional degrees of freedom from the “bulk.” However, in the present case, this equivalence operation has a trivial effect because every 2-cell is joined with two layers of 2dTI, one on each side of the 2-cell. In Methods, we give an example (without time-reversal symmetry) where this equivalence relation has a nontrivial effect.
Therefore, any two distinct topological crystal states are in different quantum phases of matter. So, to obtain a classification of TCIs, we need to enumerate possible topological crystals. First, we observe that topological crystals form an Abelian group 𝒞 under stacking, i.e., upon superposing two different states in the same space. Because MCIs (2dTIs) are characterized by ℤ (ℤ2) invariants, 𝒞 is a product of ℤ and ℤ2 factors, with the ℤ factors generated by MTCIs and the ℤ2 factors generated by ℤ2 TCIs. Because the MTCIs can be decomposed into decoupled planar layers, there is one ℤ factor for each symmetry-inequivalent set of mirror planes. For any particular crystalline symmetry of interest, the classification 𝒞 is easily worked out by considering possible colorings of the faces of the AU with the MCI and 2dTI states.
To provide a concrete illustration, we here explicitly work out the topological crystals for space group P42212. P42212 has a tetragonal lattice and is generated from three translations {1∣100}, {1∣010}, and {1∣001}, a fourfold screw
Solving Eq. 1, we get three independent states that generate all possible topological crystals under stacking: (i) n1 = n4 = 1, n2 = n3 = 0 (Fig. 2B); (ii) n1 = n2 = 1, n3 = n4 = 0 (Fig. 2C); and (iii) n3 = n4 = 1, n1 = n2 = 0 (Fig. 2D). While states (ii) and (iii) are obviously layer constructions, state (i) is not layer constructible, as discussed above.
Topological invariants
Now, we turn to the topological invariants characterizing topological crystals. First, all MTCIs are characterized by real-space Chern numbers associated with certain mirror planes, and the mirror Chern numbers in momentum space for each of them are listed in (35). Therefore, we focus on ℤ2 TCIs. Given a ℤ2 TCI and its corresponding topological crystal, for each symmetry operation g ∈ Gc, we assign a ℤ2 number δ(g). We arbitrarily choose one AU and let r be a point inside, then we set δ(g) = 1 [δ(g) = 0] if a path connecting r to gr crosses through an odd (even) number of 2dTI 2-cells. It is shown in Methods that (1) δ(g) is well defined, independent of the arbitrary choices of AU, r, and the path connecting r to gr, and (2) δ(g1g2) = δ(g1) + δ(g2). The latter property implies that δ is a homomorphism from Gc to ℤ2 [or, equivalently, an element of H1(Gc, ℤ2)], which means that it is enough to specify δ(g) for the generators of Gc. δ(g) encodes all the ℤ2 invariants for TCIs listed earlier, by choosing different operations g. For instance, if g is a translation, then δ(g) is the corresponding ℤ2 weak invariant, if g is inversion, then δ(g) is the ℤ2 inversion invariant, and so on.
As an example, the topological crystal shown in Fig. 2B has δ({1∣100}) = 0, δ({1∣010}) = 0, δ({1∣001}) = 1,
As suggested by the above discussion of invariants, the classification 𝒞 of TCIs has a simple relationship with H1(Gc, ℤ2), which allows us to efficiently compute 𝒞 and obtain the full TCI classification 𝒞 for all space groups (Methods); the results are given in Table 2.
Moreover, we find that there are 12 groups hosting topological crystals beyond layer construction (Fig. 3); for these non–layer-constructable states, we tabulated their invariants and symmetry-based indicators (19, 20) in the Supplementary Materials, completing the mapping from indicators to TCI invariants (35). Table 1 gives the classification of TCIs protected by point group symmetry for the 32 crystallographic point groups in three dimensions.
(A) to (L) show these states for space groups Pnn2 (#34), Pnnn (#48), P42 (#77), P42/n (#86), P4222 (#93), P42212 (#94), P42cm (#102),
N/A denotes a trivial classification. Full classifications including strong topological insulators can be easily obtained from this table by using the group extension rules given in Results and are also provided as a table in the Supplementary Materials.
Last, given the classification of TCIs, we obtained simple rules that extend this classification to include STIs. The key fact is that upon stacking two identical STIs together, one can either obtain a trivial state or a nontrivial TCI. Upon identifying the state thus obtained, we obtain a full classification of all topologically nontrivial insulators of noninteracting electrons with time-reversal symmetry and spin-orbit coupling.
Unified classification of STIs and TCIs
We have focused thus far on classifying TCIs, where crystalline symmetry is required to protect a nontrivial cSPT phase. Here, we address a more general problem, namely, the classification of all d = 3 free-electron insulators with time-reversal symmetry, spin-orbit coupling, and arbitrary crystalline point group or space group symmetry. We still ignore distinctions among atomic insulators, so the one new state that must be added as a generator of the classification is the STI, which is, of course, robust even upon breaking crystalline symmetry.
First, we assume that the STI is compatible with an arbitrary crystalline symmetry. We expect that this is true, but to our knowledge, it has not been proved rigorously. One argument in favor of this expectation is to note that the STI can be described by a continuum theory of a massive Dirac fermion, which is invariant under arbitrary rigid motions of 3D space. This symmetry can be broken down to an arbitrary space group or point group symmetry, for instance, by adding a periodic potential, which produces a model of an STI with arbitrary space group symmetry. This is not quite a rigorous argument because one has to show that it is possible to regularize the continuum theory in a manner compatible with an arbitrary lattice symmetry. Another argument is to note that any centrosymmetric space group has a ℤ4 indicator (20), and according to the Fu-Kane formula (37), the root state with z4 = 1 is an STI. While it was argued that any symmetry indicator can be realized by a band structure, there is no guarantee that the resulting band structure is an insulator (20). Assuming an STI can indeed be found for each centrosymmetric space group, then we need only note that every noncentrosymmetric space group is a subgroup of some centrosymmetric space group, so an STI compatible with the latter is compatible with the former. To show this expectation holds rigorously, a straightforward approach would be to find a small number of space groups that contain all space groups as subgroups and exhibit a model realizing an STI for each of these symmetry groups.
Next, we would like to compute the classification 𝒞full including both TCIs and STIs. The topological crystal picture tells us that TCIs are a subgroup (i.e., 𝒞 ⊂ 𝒞full) because stacking two TCIs produces another TCI or a trivial state. It is also true that 𝒞full/𝒞 ≃ ℤ2, because this quotient corresponds to ignoring the distinctions among TCIs, which leaves only aℤ2 generated by the STI. It follows that ∣Cfull∣ = 2 ∣𝒞∣ when the TCI classification is finite. One might expect 𝒞full = 𝒞 × ℤ2, with the ℤ2 factor generated by the STI, but this is not true in general because stacking two identical STIs can result in a nontrivial TCI. Put another way, 𝒞full can be a nontrivial group extension of 𝒞 by ℤ2, and we need to solve this group extension problem.
We proceed by choosing a particular STI state and stacking this state with itself to get a state we call (STI)2. We know that (STI)2 has trivial strong index and is, thus, either a nontrivial TCI or a trivial state, that is, (STI)2 ∈ 𝒞. We need to determine the element of 𝒞 given by (STI)2. First, we observe that our choice of STI is arbitrary under stacking with a TCI because this stacking does not change the strong invariant. It is obvious that stacking STI with a ℤ2 TCI does not affect (STI)2, but stacking STI with an MTCI can change the ℤ invariants of (STI)2 by arbitrary even integers, depending on the choice of MTCI. We, thus, see that the information in (STI)2 that is independent of the arbitrary choice of STI is precisely the information preserved under the map
Denoting the homomorphism given by π((STI)2) by δ : Gc → ℤ2, it is natural to conjecture that δ(g) = 0 when g is a rigid motion preserving the orientation of space (e.g., translations and rotations), and δ(g) = 1 when g reverses orientation (e.g., inversion, reflections, and glide reflections). This conjecture is natural because the map δ should depend only on the crystalline symmetry Gc, and there does not seem to be any other nontrivial map that can be defined in a uniform way for all Gc.
We can establish this conjecture using results of Khalaf et al. (36), where the authors studied surface theories obtained by stacking two STIs. For a crystalline symmetry
The argument is not quite complete because the crystalline symmetry Gc cannot generally be preserved by a surface termination. However, we found that specifying the seven types of invariants listed in the Results section uniquely determines a TCI phase (element of 𝒞). Therefore, we can take
This result determines the group structure of 𝒞full. There are three cases: (i) If Gc contains only orientation-preserving operations, then 𝒞full = 𝒞 × ℤ2. (ii) If Gc contains orientation-reversing operations but no mirror reflections, and hence 𝒞 has no ℤ factors, then one of the ℤ2 factors in 𝒞 is replaced in 𝒞full by a ℤ4 factor. (iii) If Gc contains mirror reflections, then (STI)2 generates a ℤ factor in 𝒞. In this case, 𝒞full and 𝒞 have the same group structure, but in 𝒞full, the generator of one of the ℤ factors is an STI state. These rules easily allow one to obtain 𝒞full for all the crystallographic point groups and space groups, and the results are given as tables in the Supplementary Materials. Because π((STI)2) is known, it is also straightforward to explicitly construct the topological crystal corresponding to (STI)2, up to the arbitrariness in defining STI. We note that recently, (43) used K-theory to obtain classifications for crystallographic point groups in d = 3. To compare these results with ours for all point groups would, for instance, require extending our results to include topological crystals with db = 0 bulding blocks. However, for the nine nontrivial point groups where all operations share a common fixed line or fixed plane (Cn, Cs, and Cnv), states with db = 0 building blocks are not relevant, and for these point groups, our results agree with those of (43).
DISCUSSION
Our method provides a unified, real-space perspective for TCIs, complementary to the momentum-space perspective usually taken for free-fermion systems. Real-space constructions have an advantage when electron interactions are considered. For example, because all ℤ2 TCIs are made from 2dTIs, from the stability of the 2dTI phase against interactions, we immediately know that the classification of ℤ2 TCIs is not affected by interactions. On the other hand, the ℤ classification of d = 2 MCI states collapses to ℤ8 in the presence of interactions (44). This implies that the classification of MTCIs does collapse, but the ℤ invariants characterizing MTCIs are robust to interactions modulo 8.
One can easily use the idea of topological crystals to classify free-electron d = 3 insulators with time-reversal but without spin-orbit coupling, that is, with SU(2) spin rotation symmetry. For these systems, time reversal and crystalline symmetry can be taken to act trivially on electron operators, with no factors of fermion parity, and in particular T2 = 1. It can then be seen that all 2-cells have symmetry class AI, while 1-cells can be in class AI or A. In all these cases, only trivial states are possible, and no free-electron TCIs can occur. This means that any nonzero symmetry-based indicator implies some topological nodes in the bulk, proved by exhaustion in (45).
The topological crystal approach developed here can be applied in many other physical settings. For instance, one can classify topological crystalline superconductors, described at the free-fermion level by Bogoliubov–de Gennes Hamiltonians. Moreover, as other works have begun to explore, it is also possible to use topological crystals to classify interacting fermionic cSPT phases.
METHODS
Cell complex structure
A cell complex is a topological space constructed by gluing together points (0-cells) and n-dimensional balls (n-cells). In more detail, to construct a cell complex X, one starts with a set of discrete points X0, referred to as the 0-skeleton. Next, one forms the 1-skeleton X1 by attaching a set of 1-cells to X0. To attach a 1-cell, one starts with a closed interval on the real line (whose interior is the 1-cell), and the two endpoints are identified with points in X0. The resulting space X1 is given by attaching the 1-cells to X0. The process continues in the natural way; for instance, to attach a 2-cell to X1, we start with a 2D disc D with boundary (whose interior is the 2-cell) and identify ∂D with a subset of X1 using a continuous map from ∂D to X1. The n-skeleton Xn is given by attaching the n-cells to Xn−1. A more detailed discussion can be found in the book by Hatcher (46).
Here, we describe in more detail how 3D space ℝ3 can be given a cell complex structure upon choosing an AU. An AU is an open subset of ℝ3 that is as large as possible, subject to the condition that no two points in the AU are related by the action of Gc. The choice of AU is not unique. While not strictly necessary, we can always choose an AU such that the boundary of the closure of the AU consists of segments of flat planes; that is, in the case of space group symmetry, the AU can be chosen as the interior of a polyhedron. Once we choose an AU, it and its copies under the action of Gc form the 3-cells, which are in one-to-one correspondence with the elements of Gc. The union of all the 3-cells is denoted by 𝒜, and its complement X2 = ℝ3 − 𝒜 is the 2-skeleton of the cell complex.
We choose 2-cells of X2 satisfying three properties: (1) Each 2-cell is a subset of a face where two 3-cells meet. To be precise, we say that two 3-cells meet at a face when the intersection of their closures is homeomorphic to a 2-manifold (possibly with boundary), and we define the face where they meet to be this intersection. (2) No two distinct points in the same 2-cell are related under the action of Gc. Note that a 2-cell may be a subset of a mirror plane, in which case the mirror symmetry will take every point in the 2-cell to itself. This property ensures that each 2-cell has no spatial symmetries; if there is a mirror symmetry, it acts on the 2-cell effectively as an internal symmetry. (3) The 2-cell structure on X2 respects the Gc symmetry. Precisely, given a 2-cell e2 and a symmetry operation g ∈ Gc, the image g(e2) is also a 2-cell.
It is always possible to choose a set of 2-cells satisfying these properties: We started with the set of faces where pairs of 3-cells meet and took their interiors as 2-manifolds. This gives a set of 2-cells for X2 satisfying properties (1) and (3), but property (2) need not be satisfied. This can be rectified by dividing up the 2-cells until property (2) is satisfied.
Letting 𝒜2 be the union of all the 2-cells, the 1-skeleton X1 is X2 − 𝒜2. We choose 1-cells to satisfy three properties very similar to those for 2-cells. A difference from the 2-cell case is that different numbers of 2-cells can meet at a 1-cell; we would like to ensure that the same set of n 2-cells meets everywhere along the extent of a given 1-cell. We therefore modified property (1) as follows: Each 1-cell is a subset of an edge where exactly n 2-cells meet. More precisely, we say that n 2-cells meet at an edge when the intersection of their closures is homeomorphic to a 1-manifold (possibly with boundary), and the edge where they meet is defined to be this intersection. Apart from these n 2-cells, we require the edge to have empty intersection with the closure of any other 2-cell. Properties (2) and (3) are required to hold with the obvious modifications.
The 0-cells are just the points where two or more 1-cells meet. Formally, letting 𝒜1 be the union of all 1-cells, the 0-cells are the points of X1 − 𝒜1.
We illustrate this rather abstract discussion with some examples. First, we considered Gc = Ci, the point group generated by inversion symmetry. We took the AU to be the half space z > 0, and the 3-cells are then the two half-spaces z > 0 and z < 0. There are two 2-cells, which are z = 0 half planes with y > 0 and y < 0, and two 1-cells, which are z = y = 0 half lines with x > 0 and x < 0. Last, the single 0-cell is the point at the origin.
As a second example, we took Gc to be space group #1, which consists only of translation symmetry. We set the lattice constant to unity and take the three Bravais lattice basis vectors to be (1,0,0), (0,1,0), and (0,0,1). A natural choice for an AU is simply the interior of a unit cell, i.e., the region 0 < x, y, z < 1. The 3-cells are then the copies of the AU under translation. There are three kinds of 2-cells. One type consists of the xy plane (i.e., z = 0 plane) region with 0 < x, y < 1 and its images under translation, and the other two types are similar but lie in the xz and yz planes. Similarly, there are three kinds of 1-cells, with one type consisting of the x = y = 0 region with 0 < z < 1 and its images under translation. The other two types are similar regions oriented along the x and y axes. Last, the 0-cells are points (nx, ny, nz), with nx, ny, and nz integers.
Formal structure of classification resolved by block dimension
Here, we describe the Abelian group structure of TCIs (32) and how taking a certain quotient allows us to ignore distinctions among atomic insulators. We let 𝒟db be the Abelian group classifying insulators whose nontrivial building blocks are dimension db and below. That is, states classified by 𝒟db can be reduced to a state on X2 where all n-cells with n > db host a trivial state. Clearly, db = 0, 1, and 2, and we have the sequence of subgroups 𝒟0 ⊂ 𝒟1 ⊂ 𝒟2. 𝒟2 as the classification of all TCIs or, at least, all those that can be classified in terms of topological crystals. The observation that all 1-cells are trivial implies that 𝒟0 = 𝒟1. Phases in 𝒟0 are atomic insulators, which we wish to exclude from consideration. Although there are distinct atomic insulators constituting different quantum phases of matter, all atomic insulators are, in some sense, topologically trivial. We can eliminate these states by taking the quotient 𝒞 = 𝒟2/𝒟0, which gives the desired classification of TCIs.
Gluing MCI building blocks: Planar decomposition of MTCIs
Here, we address the consequences of the gluing conditions for MTCIs, i.e., topological crystals built from MCI building blocks placed on the 2-cells of X2. In particular, we show that MTCIs can always be decomposed into decoupled planar MCI states placed on mirror planes. We consider a mirror plane P, and note that we must have P ⊂ X2. Therefore, up to a set of measure zero, P is a union of 2-cells of X2. We consider a 2-cell
To proceed, we consider a number of cases. In case (1), there exists an element g ∈ Gc that maps
The relativistic Dirac Hamiltonian as discussed in the Supplementary Materials allows us to determine the presence or absence of the (−1)F factor. We choose coordinates so that the mirror plane is the z = 0 plane, and the action of σ on the Dirac field Ψ(r) is given in the Supplementary Materials. We consider a symmetry operation g that takes the mirror plane into itself, with action on the Dirac field
Now we use this claim to show that the MCI states on
To complete the discussion of case (1), we need to address gluing of the two MCI states at e1. It is enough to consider only symmetries that take the set of cells
Now, we move on to case (2), where there is no element g ∈ Gc that maps
In case (2b),
Last, in case (2b.ii), e1 is contained in a C3v axis, where the C3v symmetry is generated by σ and a 3-fold rotation C3. Here, there are six 2-cells that coincide with the mirror planes meeting at e1. These 2-cells come in three pairs, with each pair contained in one of the three mirror planes that intersect e1.
The images of these fermions under C3 rotation are
Using these relations, we find
Now, focusing only on the mirror reflection symmetry, we can change variables to diagonalize σ and find that the cR2, cR3, cL2, and cL3 fermion modes can be gapped out. This leaves the cR1 and cL1 edge of the MCI state on
Gluing condition for ℤ2 TCIs
Here, we consider the gluing condition for ℤ2 TCIs, i.e., the requirement that there are no gapless modes within the bulk. If we consider placing 2dTI states on a subset of the 2-cells of X2, then we will show that the gluing condition can be satisfied if and only if an even number of 2dTI edges meet at each 1-cell. This is illustrated in fig. S4 for the space groups P3 (no. 143) and P4 (no. 75). One direction is trivial to show: If the gluing condition is satisfied, then an even number of 2dTI edges must meet at each 1-cell e1 because, otherwise, time reversal would forbid e1 from being gapped.
Now, we suppose that we place 2dTI states on 2-cells of X2 such that an even number of 2dTI edges meets at each 1-cell. We would like to show that each 1-cell can be gapped; thus, the gluing condition is satisfied. We do this by considering the different possible point group symmetries of a 1-cell e1, which may impose constraints on gluing of 2dTI edge modes along e1. If e1 has trivial point group symmetry, then the only symmetries are charge conservation and time reversal, and an even number of 2dTI edge modes can always be gapped. If e1 is contained in a mirror plane, then 2dTI edge modes come in mirror-symmetric pairs above and below the mirror plane, and we have already shown in the above that each such pair can be gapped.
Next, we consider the case where e1 is contained in a Cn axis. When n = 2, 2dTI edge modes come in pairs related by C2 symmetry, and the action of C2 symmetry on such a pair is identical to that of the σ′ symmetry discussed above (see Eq. 7 and the surrounding discussion). Therefore, these pairs of edge modes can be gapped. When n > 2 is even, edge modes come in groups of n related by Cn symmetry, and these can be grouped into n/2 pairs related by C2 symmetry. We can focus on one such pair and gap it out, then use the Cn symmetry to “copy” its mass term to the other n/2 − 1 pairs, which gaps out all the modes respecting the Cn symmetry. Last, for C3 symmetry, edge modes must come in groups of six, with two sets of three edge modes related by symmetry. We can take a pair of edge modes unrelated by symmetry and gap these out and then use the C3 symmetry to copy the resulting mass term to the other two pairs of edge modes.
Similar arguments can be applied when e1 is contained in a Cnv axis. Here, 2-cells that can host 2dTIs lie away from the mirror planes, so that edge modes come in groups of 2n related by Cnv symmetry. Viewing Cnv as generated by a mirror reflection and Cn rotation, we can start by gapping out a pair of neighboring edge modes related by the mirror symmetry and then copying its mass term using the rotational symmetry.
Gluing at 0-cells
In the analysis of the gluing condition in the above, we started with a set of topological states on 2-cells and considered gluing these states together at 1-cells. In principle, this may not be the end of the story; we needed to consider gluing at 0-cells. That is, we needed to ensure that there are no localized protected gapless states at 0-cells, which would violate the gluing condition. However, there is a simple reason this cannot occur: consider the set of gapped 1-cells that meet at a 0-cell. We can lump these 1-cells together and view them as a single gapped d = 1 system, with the 0-cell as its endpoint. Upon decomposing the spectrum into irreducible representations of any site symmetry at the 0-cell, this d = 1 system divides into sectors that are either in class A or AII. To have a protected gapless state on the 0-cell, the d = 1 system would need to be topologically nontrivial, but class A and AII have a trivial classification in d = 1.
Nontrivial example of equivalence operation
In the Results section, we discussed an equivalence operation involving creating bubbles of 2dTI within each AU, which turns out to have a trivial effect on the classification of time-reversal symmetric electronic TCIs. To clarify the nature and role of this equivalence operation, here we discuss an example without time-reversal symmetry where it does modify the classification. We consider d = 3 insulators without time-reversal symmetry with point group
(A) The 2-cells:
First, we consider the gluing condition on X2. X2 decompose into two 2-cells:
We now consider the equivalence operation. We create a bubble of Chern insulator with Chern number 1 in the z > 0 region. Because of the inversion symmetry, another bubble will be created in the z < 0 region as shown in Fig. 4B. We enlarge the two bubbles respecting the inversion symmetry, until the bottom surface of the z > 0 bubble joins with X2 and the rest of the bubble moves to infinity. At the same time, the top surface of the z < 0 bubble also joins with X2. Therefore, after this operation, the total Chern number on X2 is increased by 2, which implies that the classification is reduced from ℤ to ℤ2.
Topological crystals, topological invariants, and H1(G, ℤ2)
In this section, we give a more detailed discussion of the ℤ2-valued function of invariants δ(g) characterizing a topological crystal. On the basis of this discussion, we established a connection between the classification 𝒞 of TCIs and H1(Gc, ℤ2), which allows 𝒞 to be computed with the aid of standard computer algebra tools such as GAP (49). In particular, we show that 𝒞 and H1(Gc, ℤ2) have the same number of generators.
In the text, we defined δ(g) only for a ℤ2 TCI. It will be useful for our present purposes to give a definition valid for an arbitrary topological crystal. First, we introduce the notion of a ℤ2 coloring of X2, which is given by associating a ℤ2 number to each 2-cell of X2 so that each 2-cell is either colored and assigned 1, or empty and assigned 0. These ℤ2 numbers must be assigned to respect the crystalline symmetry and satisfy a gluing condition, namely, for each 1-cell, an even number of the 2-cells meeting there must be colored. ℤ2 colorings can be added using the ℤ2 addition law, and this makes ℤ2 colorings into a group that we denote by
Next, we observe that there is a map from topological crystals (elements of 𝒞) to
Last, given a ℤ2 coloring of X2, we define δ(g) as in the Results section. That is, we arbitrarily choose one AU, let r be a point inside, and define δ(g) = 1 if a path connecting r to gr crosses an odd number of colored 2-cells, while δ(g) = 0 if the path crosses an even number of colored 2-cells. The path should be chosen to avoid 0-cells and 1-cells but is otherwise an arbitrary continuous path.
We now establish some properties of δ(g) quoted in the Results section. First, we show that δ(g) is independent of the path chosen to connect r to gr and is, thus, a well-defined function mapping Gc to ℤ2. Any two such paths are related by a finite number of moves, where a segment of the path is passed through a 1-cell. The gluing condition says that an even number of colored 2-cells meet at every 1-cell, so these moves do not affect δ(g). At this stage, we have not yet shown that δ is independent of the arbitrary choice of AU.
Next, we show that δ is a homomorphism from Gc to ℤ2, that is, δ(g1g2) = δ(g1) + δ(g2). We compute δ(g1g2) by considering a path from r to g1g2r that first goes from r to g1r and then goes from g1r to g1g2r. The number of colored 2-cells modulo 2 crossed by the first segment is δ(g1) by definition. The second segment is related by symmetry to a path joining r to g2r, so δ(g2) is the number of colored 2-cells (modulo 2) crossed by the second segment. Therefore, δ(g1g2) = δ(g1) + δ(g2).
Last, we show that δ does not depend on the arbitrary choice of AU. Let δ(g) be the function defined by choosing an AU with a point r inside, and let δ′(g) be the function defined by choosing a different AU, which contains a point g0r for some g0 ∈ Gc. Then, δ′(g) is the number of colored 2-cells (modulo 2) crossed by a path connecting g0r to gg0r. By symmetry, this number is the same as for a path joining r to
Our construction of δ gives a map
Now that we have shown that
Using the above discussion, we showed that the map
These results make it a simple matter to compute the TCI classification 𝒞. First, we compute H1(Gc, ℤ2); this can be done using GAP (49). Then, we know that the number of ℤ factors in 𝒞 is nM, the number of symmetry-inequivalent sets of mirror planes. We then obtain 𝒞 from H1(Gc, ℤ2) by replacing nM of the ℤ2 factors with ℤ factors. For a Gc crystallographic point group or space group, nM can be obtained immediately from information tabulated in the International Tables for Crystallography (50). The results of this procedure are presented in Table 1 for crystalline point groups, and in Table 2 for space groups. As discussed in the Introduction, we note that Khalaf et al. (36) have also obtained the results in these tables via a mathematically equivalent procedure.
N/A denotes a trivial classification. Full classifications including strong topological insulators can be easily obtained from this table by using the group extension rules given in Results and are also provided as a table in the Supplementary Materials.
While useful, we emphasize that this largely automated procedure is not a substitute for explicit real-space construction of topological crystals for a given symmetry group of interest, as described in Results. The latter procedure not only results in the same group structure but also provides additional physical insight and a starting point for further analysis, by giving an explicit real-space construction of each of the TCI phases classified. We also obtained the results in Table 2 by automating the explicit real-space constructions.
SUPPLEMENTARY MATERIALS
Supplementary material for this article is available at http://advances.sciencemag.org/cgi/content/full/5/12/eaax2007/DC1
Fig. S1. Illustration of the 2-cells
Fig. S2. The 3-cells
Fig. S3. Cross section through e1 and the 2-cells that meet at e1, in case (2b.ii), where e1 is contained in a C3v axis.
Fig. S4. Illustration of the effects of the gluing condition for ℤ2 TCIs for two representative space groups.
Table S1. Topological crystals beyond layer construction.
Table S2. Classification 𝒞full of STIs and TCIs with time-reversal symmetry and spin-orbital coupling for all point groups.
Table S3. Classification 𝒞full of STIs and TCIs with time-reversal symmetry and spin-orbital coupling for all space groups.
This is an open-access article distributed under the terms of the Creative Commons Attribution-NonCommercial license, which permits use, distribution, and reproduction in any medium, so long as the resultant use is not for commercial advantage and provided the original work is properly cited.
REFERENCES AND NOTES
- Copyright © 2019 The Authors, some rights reserved; exclusive licensee American Association for the Advancement of Science. No claim to original U.S. Government Works. Distributed under a Creative Commons Attribution NonCommercial License 4.0 (CC BY-NC).