Exotic topology on the surface
Analyzing the spatial symmetries of three-dimensional (3D) crystal structures has led to the discovery of exotic types of quasiparticles and topologically nontrivial materials. Wieder et al. focus on the symmetry groups of 2D surfaces of 3D materials—the so-called wallpaper groups—and find that some of them allow for an additional topological class. This class hosts a single fourfold-degenerate Dirac fermion on the surface of the material and, on the basis of the authors' calculations, is expected to occur in the compound Sr2Pb3.
Science, this issue p. 246
Abstract
Materials whose gapless surface states are protected by crystal symmetries include mirror topological crystalline insulators and nonsymmorphic hourglass insulators. There exists only a very limited set of possible surface crystal symmetries, captured by the 17 “wallpaper groups.” Here we show that a consideration of symmetry-allowed band degeneracies in the wallpaper groups can be used to understand previously described topological crystalline insulators and to predict phenomenologically distinct examples. In particular, the two wallpaper groups with multiple glide lines, pgg and p4g, allow for a topological insulating phase whose surface spectrum consists of only a single, fourfold-degenerate, true Dirac fermion, representing an exception to a symmetry-enhanced fermion-doubling theorem. We theoretically predict the presence of this phase in Sr2Pb3 in space group 127 (P4/mbm).
Topological phases stabilized by crystal symmetries have a number of distinct realizations. The first class that was theoretically proposed is composed of rotation and mirror topological crystalline insulators (TCIs) that host surface fermions protected by the projection of a bulk mirror plane or rotation axis onto a particular surface (1–3). These fermions have been observed in SnTe (4, 5) and related compounds (6, 7). Recent efforts to expand this analysis to nonsymmorphic systems with surface glide mirrors—symmetries composed of a mirror and a half-lattice translation—have yielded additional exotic free-fermion topological phases, which can exhibit so-called surface gapless “hourglass fermions” and the glide spin Hall effect (8–12). The theoretical proposal of (9, 10) has recently also received incipient experimental support (13).
In addition, topological insulators (TIs), crystalline or otherwise, provide exceptions to fermion-doubling theorems. These theorems impose fundamental bounds on phenomena in condensed matter physics. For example, in two dimensions, a single Kramers degeneracy in momentum space must always have another, partner Kramers crossing elsewhere in the Brillouin zone (BZ); otherwise, the Berry phase of a loop enclosing the degeneracy suffers from ambiguity (14). The discovery of the TI provided the first exception to this theorem: In these systems, Kramers pairs are allowed to be isolated on a single two-dimensional (2D) surface because they are connected to a 3D topological insulating bulk and have their partners on the opposite surface (15, 16).
Bulk fermions with higher-fold degeneracy, such as Dirac points, which are stabilized by crystal symmetry, may come with their own fermion-doubling theorems (17–20). As noted in (21), a single fourfold-degenerate Dirac fermion cannot be stabilized by 2D crystal symmetries as the only nodal feature at a given energy; it must always have at least one partner or accompanying hourglass nodal points. This is because a single Dirac point in two dimensions represents the quantum critical point (QCP) separating a trivial insulator (NI) from a TI. As shown in more detail in section 1B of (22), stabilizing just one of these Dirac points with crystal symmetries would therefore force the broken-symmetry NI and TI phases to be related by just a unitary transformation, violating their
The crucial requirement is that the surface preserves multiple nonsymmorphic crystal symmetries. Until now, most attention has been paid to crystal systems with surfaces that preserve only a single glide mirror. However, of the 17 2D surface symmetry groups, called wallpaper groups, there are two that host two intersecting glide lines (23). As we show in section 2 of (22), the algebra of the two glides requires that bands appear with fourfold degeneracy at a single time-reversal–invariant momentum (TRIM) at the edge of the BZ.
We study the noninteracting topological phases allowed in bulk crystals with surfaces that are invariant under the symmetries of the two wallpaper groups mentioned above, pgg and p4g. We show that, in addition to generalizations of the hourglass fermions introduced in (9), they host a topological phase characterized by a single, symmetry-enforced fourfold Dirac surface fermion—i.e., twice the degeneracy of a traditional TI surface state. This Dirac fermion is pinned by nonsymmorphic symmetry to the QCP between a TI and a NI, allowing for controllable topological phase transitions of the 2D surface under spin-independent glide-breaking strain.
Wallpaper groups pgg and p4g
The surface of a nonmagnetic crystal is itself a lower-dimensional crystal that preserves a subset of the bulk crystal symmetries; all 2D nonmagnetic surfaces are geometrically constrained to be invariant under the action of the 17 wallpaper groups. The set of spatial wallpaper group symmetries is restricted to those 3D space group symmetries that preserve the surface normal vector: rotations about that vector and in-plane lattice translations, mirror reflections, and glide reflections. Surface band features will therefore be constrained by the irreducible co-representations of these symmetries and their momentum-space compatibility relations (24, 25). Focusing on the 17 wallpaper groups that describe the surfaces of 3D crystals with strong spin-orbit coupling (SOC) and time-reversal symmetry [section 2 of (22)], we designate all topological 2D surface nodes formed from symmetry-enforced band degeneracies “wallpaper fermions.” In the language of group theory, wallpaper fermions are therefore fully captured by the irreducible co-representations of the wallpaper groups, so that, for example, the fourfold rotation–protected quadratic topological node introduced in (2) is a “spinless wallpaper fermion”; in contrast, the magnetic twofold surface degeneracy in (26), for which both bands have the same co-representations and are instead prevented from gapping by an additional rotation- and time-reversal–enforced
Although other works have focused on the mathematical classification of topological phases protected by wallpaper group symmetries (27, 28), here we detail a further topological classification, as well as a search for topological insulating phases with phenomenologically distinct surface states. For the symmorphic wallpaper groups, topological rearrangements of the minimal surface band connectivities, recently enumerated for 3D bulk systems in (25), allow for only quantum spin Hall (QSH) phases and rotational variants of the mirror TCI phases, which exhibit twofold-degenerate free fermions along high-symmetry lines in the surface BZ (4–7).
For the four wallpaper groups with nonsymmorphic glide lines (pg, pmg, pgg, and p4g), this picture is enriched. Even in two dimensions, glide symmetries require that groups of four or more bands be connected, an effect that frequently manifests itself in hourglass-like band structures (21, 29, 30). For the wallpaper groups with only a single glide line, pg and pmg, surface bands can be connected in topologically nontrivial patterns of interlocking hourglasses (9, 10) or in “Möbius strip” connectivities if time-reversal symmetry is relaxed (26, 31), both of which exhibit twofold-degenerate surface fermions along some of the momentum-space glide lines. In the remaining two wallpaper groups with multiple glide symmetries, pgg and p4g, higher-degenerate wallpaper fermions are uniquely allowed.
We consider a z-normal surface with glides
Shown are the unit cells and Brillouin zone (BZ) for their two-site realizations. The A and B sites are characterized by T-symmetric internal degrees of freedom (blue arrows) that are transformed under crystal symmetry operations. These correspond physically to nonmagnetic properties that transform as vectors, such as atom displacements or local electric dipole moments. Glide lines (green) exchange the sublattices through fractional lattice translations. In p4g, there is an extra C4 symmetry (⊗) about the surface normal with axes located on the sites. The combination of this C4 symmetry and the glides produces additional diagonal symmorphic mirror lines (red) in p4g. In the BZ of p4g, C4 relates
For bulk insulators, the glide-preserving bulk topological phase and, consequently, z-normal surface states can be determined without imposing a surface by classifying the allowed connectivities of the z-projection Wilson loop holonomy matrix (10, 32–38), a bulk quantity defined by
By generalizing the
To enumerate the allowed topological phases shown in Fig. 2, we consider possible restrictions on (χx, χy). Although χx and χy can individually take on values 0, 1, 2, and 3, only pairs that satisfy χx + χy mod 2 = 0 are permitted in bulk insulators. This can be understood as follows: If χx + χy is odd, the 2D surface consisting of the four partial planes defined by 0 ≤ kx ≤ π ky = 0, π and kx = 0, π 0 ≤ ky ≤ π possesses an overall Chern number, which implies the existence of a bulk gapless point (40, 41), contradicting our original assumption that the system is insulating. We present a rigorous proof in section 5B of (22) and show that the remaining collection of eight insulating phases is indexed by the group
Shown are the eight topologically distinct Wilson band connectivities for bulk insulators with crystal surfaces that preserve 2D glide reflections on the projections of two orthogonal bulk glide planes. Each band structure is labeled by its two
For χx,y = 1, 3, the system is a strong topological insulator (STI). The four “double-glide spin Hall” STI phases possess the usual twofold-degenerate Kramers pairs at
When χx,y = 0, 2, the system is in a topological crystalline phase. For (χx, χy) = (0, 2) or (2, 0), which is only permitted in a C4z-broken surface pgg, a variant of the hourglass insulating phase (9) is present on the surface. For example, when (χx, χy) = (0, 2), either time-reversed partners of twofold-degenerate free fermions live along
Lastly, but most interestingly, for χx = χy = 2, we find that the system exists in a previously uncharacterized “nonsymmorphic Dirac insulating” phase, capable of hosting just a single fourfold-degenerate Dirac surface fermion at
Materials realizations
We applied density functional theory (DFT), including the effects of SOC, to predict topological phases stabilized by wallpaper groups pgg and p4g in previously synthesized, stable materials. The details of these large-scale calculations are provided in section 6 of (22). We find double-glide topological phases on the (001) surface (wallpaper group p4g) of three members of the space group 127 (P4/mbm) A2B3 family of materials: Sr2Pb3 (42, 43), Au2Y3 (44), and Hg2Sr3 (45, 46). We find that Sr2Pb3 has a gap at the Fermi energy at each crystal momentum, in spite of the presence of electron and hole pockets (Fig. 3). A Wilson loop calculation of the bands up to this gap (Fig. 3D) indicates that this material possesses the bulk topology of a (2, 2) nonsymmorphic Dirac insulator [section 6B of (22)]. Calculating the surface spectrum through surface Green’s functions, as described in section 6B of (22) (Fig. 4), we find that the (001) surface of Sr2Pb3, although displaying an overall metallic character, develops a gap of 45 meV at the Fermi energy at
This compound is in space group 127 (P4/mbm). (A) The unit cell of Sr2Pb3. Mz is a symmorphic mirror reflection through the z axis. (B) The bulk BZ and the (001)-surface BZ (wallpaper group p4g). (C) The electronic bands obtained using DFT. The Fermi level is set to 0 eV. At each point in the BZ, there is a band gap near the Fermi energy, indicated by the dashed red line; insets show magnified images of the boxed regions. (D) The (001)-directed Wilson bands; red (blue) points indicate Wilson bands with positive (negative) surface glide eigenvalues,
Shown is the (001)-surface band structure of Sr2Pb3 (wallpaper group p4g) calculated using surface Green’s functions [section 6B of (22)]. The color bar indicates the relative degree of surface localization. The Fermi level is set to zero. The fourfold surface Dirac fermion appears at
Unlike Sr2Pb3, Au2Y3 and Hg2Sr3 in space group 127 (P4/mbm) are gapless, with bulk C4z-protected Dirac nodes (47) present near the Fermi energy. In section 6B of (22), we show that under weak (100) strain, these Dirac nodes can be gapped to induce the (0, 2) topological hourglass phase in these two materials.
We additionally find that the (001) surface (wallpaper group pgg) of the narrow-gap insulator Ba5In2Sb6 in space group 55 (Pbam) (48) hosts a double-glide topological hourglass fermion. We find that Ba5In2Sb6 develops an indirect band gap of 5 meV (direct band gap, 17 meV) (Fig. 5). The Wilson loop spectrum obtained from the occupied bands (Fig. 5D) demonstrates that this material is a (2, 0) double-glide topological hourglass insulator [section 6A of (22)]. We find that the (001) surface of Ba5In2Sb6 has a projected insulating bulk gap that is spanned along
This compound is in space group 55 (Pbam). (A) The unit cell of Ba5In2Sb6. (B) The bulk BZ and the (001)-surface BZ (wallpaper group pgg). (C) The electronic bands obtained using DFT. The Fermi level is set to 0 eV. There is an insulating gap at the Fermi energy, indicated by the dashed black line. (D) The (001)-surface Wilson bands; red (blue) points indicate Wilson bands with positive (negative) surface glide eigenvalues,
Shown is the (001)-surface band structure of Ba5In2Sb6 (wallpaper group pgg) calculated using surface Green’s functions [section 6A of (22)]. The Fermi level is set to zero. Surface bands from the top of one hourglass fermion and the bottom of another connect the valence and conduction manifolds along
Discussion
We have demonstrated theoretically the existence of a nonsymmorphic Dirac insulator—a topological crystalline material with a single fourfold-degenerate surface Dirac point stabilized by two perpendicular glides. After an exhaustive study of the 17 time-reversal–symmetric, strong-SOC wallpaper groups, only pgg and p4g are shown to be capable of supporting this fourfold fermion. This phase is one of eight topologically distinct phases that can exist in insulating orthorhombic crystals with surfaces that preserve two perpendicular glides; we have classified all eight phases by topological indices (χx, χy) that characterize the connectivity of the z-projection Wilson loop spectrum. We report the prediction of the nonsymmorphic Dirac insulating phase in Sr2Pb3 and of related double-glide topological hourglass phases in Ba5In2Sb6, as well as in (100)-strained Au2Y3 and Hg2Sr3. We also report the theoretical prediction of a set of previously unknown double-glide spin Hall phases. Although their surface Kramers pairs and fourfold Dirac fermions should be distinctive in ARPES experiments, characterization of transport in the double-glide spin Hall phases remains an open question.
We also find that a simple intuition exists for the topological crystalline phases χx,y = 0, 2. In section 7A of (22), we present an eight-band tight-binding model that, when half-filled, can be tuned to realize all
Lastly, because the (2, 2) topological surface Dirac point is symmetry-pinned to the QCP between a 2D TI and a NI, we examine its potential for hosting strain-engineered topological physics. Consider the two-site surface unit cell in wallpaper group pgg from Fig. 1. In the (2, 2) nonsymmorphic Dirac insulating phase, the surface Dirac fermion can be captured by the k · p Hamiltonian near
Supplementary Materials
This is an article distributed under the terms of the Science Journals Default License.