Topological Insulators with Inversion Symmetry

TL;DR

The work develops a parity-based framework for diagnosing $Z_2$ topological insulators in systems with inversion symmetry, reducing the problem to parity eigenvalues at time-reversal-invariant momenta. It connects bulk topological invariants to robust surface states and provides concrete predictions for strong topological insulators among Bi$_{1-x}$Sb$_x$, strained α-Sn, HgTe, and Pb$_{1-x}$Sn$_x$Te. The paper also analyzes tight-binding models (Graphene, Diamond, BHZ) to illustrate the method and discusses experimental implications, including surface transport, ARPES signatures, and potential half-quantized Hall responses. These results offer a practical route to identifying and studying strong topological insulators in real materials, while highlighting the role of inversion symmetry and the limitations imposed by disorder and interactions.

Abstract

Topological insulators are materials with a bulk excitation gap generated by the spin orbit interaction, and which are different from conventional insulators. This distinction is characterized by Z_2 topological invariants, which characterize the groundstate. In two dimensions there is a single Z_2 invariant which distinguishes the ordinary insulator from the quantum spin Hall phase. In three dimensions there are four Z_2 invariants, which distinguish the ordinary insulator from "weak" and "strong" topological insulators. These phases are characterized by the presence of gapless surface (or edge) states. In the 2D quantum spin Hall phase and the 3D strong topological insulator these states are robust and are insensitive to weak disorder and interactions. In this paper we show that the presence of inversion symmetry greatly simplifies the problem of evaluating the Z_2 invariants. We show that the invariants can be determined from the knowledge of the parity of the occupied Bloch wavefunctions at the time reversal invariant points in the Brillouin zone. Using this approach, we predict a number of specific materials are strong topological insulators, including the semiconducting alloy Bi_{1-x} Sb_x as well as α-Sn and HgTe under uniaxial strain. This paper also includes an expanded discussion of our formulation of the topological insulators in both two and three dimensions, as well as implications for experiments.

Figures (7)

• Figure 1: (a) A two dimensional cylinder threaded by magnetic flux $\Phi$. When the cylinder has a circumference of a single lattice constant $\Phi$ plays the role of the edge crystal momentum $k_x$ in band theory. (b) The time reversal invariant fluxes $\Phi=0$ and $h/2e$ correspond to edge time reversal invariant momenta $\Lambda_1$ and $\Lambda_2$. $\Lambda_a$ are projections of pairs of the four bulk time reversal momenta $\Gamma_{i=(a\mu)}$, which reside in the two dimensional Brillouin zone indicated by the shaded region. (c) In 3D the generalized cylinder can be visualized as a "corbino donut", with two fluxes, which correspond to the two components of the surface crystal momentum. (d) The four time reversal invariant fluxes $\Phi_1$, $\Phi_2 = 0$, $h/2e$ correspond to the four two dimensional surface momenta $\Lambda_a$. These are projections of pairs of the eight $\Gamma_{i=(a\mu)}$ that reside in the bulk 3D Brillouin zone.
• Figure 2: Schematic representations of the surface energy levels of a crystal in either two or three dimensions as a function of surface crystal momentum on a path connecting $\Lambda_a$ and $\Lambda_b$. The shaded region shows the bulk continuum states, and the lines show discrete surface (or edge) bands localized near one of the surfaces. The Kramers degenerate surface states at $\Lambda_a$ and $\Lambda_b$ can be connected to each other in two possible ways, shown in (a) and (b), which reflect the change in time reversal polarization $\pi_a\pi_b$ of the cylinder between those points. Case (a) occurs in topological insulators, and guarantees the surface bands cross any Fermi energy inside the bulk gap.
• Figure 3: Diagrams depicting four different phases indexed by $\nu_0; (\nu_1\nu_2\nu_3)$. The top panel depicts the signs of $\delta_i$ at the points $\Gamma_i$ on the vertices of a cube. The bottom panel characterizes the band structure of a 001 surface for each phase. The solid and open circles depict the time reversal polarization $\pi_a$ at the surface momenta $\Lambda_a$, which are projections of pairs of $\Gamma_i$ which differ only in their $z$ component. The thick lines indicate possible Fermi arcs which enclose specific $\Lambda_a$.
• Figure 4: (a) Honeycomb lattice of graphene, with a unit cell indicated by the dashed lines. (b) Brillouin zone, with the values of $\delta_i$ associated with the time reversal invariant momenta labeled. $\tau_{1/2}$ describes the loop enclosing half the zone used in Eq. \ref{['spinless']}.
• Figure 5: Schematic representation of band energy evolution of ${\rm Bi}_{1-x} {\rm Sb}_x$ as a function of $x$. Adapted from Ref. lenoir.