Robust spin-polarized midgap states at step edges of topological crystalline insulators

See allHide authors and affiliations

Science  09 Dec 2016:
Vol. 354, Issue 6317, pp. 1269-1273
DOI: 10.1126/science.aah6233

An edge that is hard to get rid of

A distinguishing characteristic of topological insulators (TIs) is that they have conducting states on their boundary—a surface for a three-dimensional (3D) TI or a line edge for a 2D TI. Sessi et al. used scanning tunneling spectroscopy to discover unusual 1D states in a 3D crystalline TI. The states appeared on the edge of a particular kind of step in the crystal and survived large magnetic fields and increased temperatures. This robustness bodes well for the potential use of these states in practical applications.

Science, this issue p. 1269


Topological crystalline insulators are materials in which the crystalline symmetry leads to topologically protected surface states with a chiral spin texture, rendering them potential candidates for spintronics applications. Using scanning tunneling spectroscopy, we uncover the existence of one-dimensional (1D) midgap states at odd-atomic surface step edges of the three-dimensional topological crystalline insulator (Pb,Sn)Se. A minimal toy model and realistic tight-binding calculations identify them as spin-polarized flat bands connecting two Dirac points. This nontrivial origin provides the 1D midgap states with inherent stability and protects them from backscattering. We experimentally show that this stability results in a striking robustness to defects, strong magnetic fields, and elevated temperature.

The recent theoretical prediction and experimental realization of topological insulators (TIs) have considerably extended the notion of a phase of matter. Within this framework, it has been shown that—based on some topological invariants—the electronic properties of materials can be classified into distinct topological classes (1, 2). In topologically nontrivial materials, unconventional boundary modes have been experimentally detected by several different techniques (39). In two-dimensional (2D) TIs, counter-propagating spin-momentum–locked 1D edge modes develop along the sample boundary; in contrast, 3D TIs (4) have boundary modes that are linearly dispersing chiral surface states. Although a large variety of 3D TIs have been reported, only very few 2D TIs are known [HgTe (3), InAs (10) quantum wells, and Bi bilayers (11)]. These 2D TIs are delicate and difficult to realize experimentally because they all require the fabrication of precisely controlled thin film heterostructures. Properties such as small band gaps (3, 10), strong substrate-induced hybridization effects (11), or the existence of residual trivial states (10, 11) make helical edge states not only challenging to study but also of limited appeal for applications. Furthermore, their topological properties are protected only as long as time-reversal symmetry is preserved.

Here, we report that 2D topological surfaces, in turn, can be the mother state for nontrivial 1D midgap states (12), suggesting a dimensional hierarchy of boundary states in topological insulators. Specifically, we report on the discovery of 1D topological spin-filtered channels that naturally develop at step edges of 3D topological crystalline insulators (TCIs)—i.e., materials where the existence of surface Dirac states is guaranteed by crystal symmetries.

Figure 1A displays the rock-salt structure of Pb1−xSnxSe (x ≤ 0.4). Depending on Sn-content x, these compounds have been reported to belong to two topologically distinct phases (13) that can be stoichiometrically controlled. Starting from PbSe, which is topologically trivial, the substitutional solid solution Pb1−xSnxSe turns into a topologically nontrivial phase as the Sn concentration exceeds x ≈ 0.24 (14, 15). Its bulk inverted band gap cannot be adiabatically connected to a trivial state as long as some crystal symmetries are preserved. The electronic properties of high-symmetry surfaces of these TCIs are characterized by topologically protected linearly dispersing Dirac states (13, 14, 1619). Figure 1B illustrates this scenario for the nonpolar (001) surface, which is commonly exposed after cleaving bulk crystals. It hosts four Dirac cones centered in close proximity to the Embedded Imageand Embedded Image points of the Brillouin zone. Figure 1C shows a typical image of the (001) surface acquired by scanning tunneling microscopy (STM) on a freshly cleaved Pb0.67Sn0.33Se bulk crystal (20)—i.e., a material safely within the topological regime at x ≥ 0.24. An atomically resolved image showing the Se sublattice (21) is displayed as an inset. The profile taken along the gray line shows that several atomically flat terraces exist, separated by step edges of different heights.

Fig. 1 Electronic properties of Pb0.67Sn0.33Se terraces and step edges probed by STS.

(A) Rock-salt crystal structure and (B) schematic band structure of (Pb,Sn)Se. (C) Topographic STM image of a cleaved Pb0.67Sn0.33Se surface (scan parameters: U = −75 mV; I = 50 pA). (Inset) Atomic resolution image of the Se sublattice. Two steps are visible in the main panel. The line section (bottom panel) measured along the gray line shows that their heights correspond to a single- (right) and a double-atomic step (left), respectively. Whereas the periodicity is maintained for even step edges (D), odd step edges lead to a structural show π shift (E). (F) dI/dU map (top) measured at the same location as (C). The line section (bottom) reveals an enhanced conductance at the position of the single-atomic step edge. (G) Local tunneling spectra measured with the STM tip positioned at the locations indicated in (F). The spectra measured on atomically flat terraces (1, 3, and 5) and at even step edges (2) display the typical V shape with a minimum at the Dirac energy (ED = −75 meV) surrounded by two maxima indicating van Hove singularities (L = −110 meV and L+= −30 meV); the spectrum measured at the position of the odd step (4) exhibits a strong peak at the Dirac energy.

Given the equal probability of breaking the crystal bonds at every atomic layer, all steps can be mapped onto two different classes: those corresponding to an integer number n of the conventional unit cell heights a—i.e., na, from here on even steps—and those corresponding to half-integer multiples—i.e., (1/2 + n)a, called odd steps hereafter. Whereas even steps maintain the translational symmetry of the surface lattice, odd steps introduce a structural π shift (Fig. 1, D and E, shows the two situations from the top). As will be described below, this has far-reaching implications for the electronic structure of the step edges.

Figure 1F shows a differential conductance (dI/dU) map at an energy of EEF = −75 meV, which was measured simultaneously with the topographic image presented in Fig. 1C. Whereas the local density of states (DOS) is similar for all terraces and also remains essentially unchanged for the even step edge, a strong enhancement can be recognized along the odd step. Its intensity is symmetrically distributed on both sides of the step and has a width of ~10 nm (line profile in Fig. 1F). Remarkably, this DOS is very homogeneously distributed along the step; no evidence for scattering-induced quasiparticle interference pattern was found [for the Fourier-transformed dI/dU intensity along the step edge, see (20) and fig. S1].

To unravel the origin of this prominent electronic feature, we measured local scanning tunneling spectroscopy (STS) data on the three atomically flat terraces and the two step edges. The precise locations where the STS curves displayed in Fig. 1G have been obtained are indicated by numbers 1 to 5 in Fig. 1F. As expected, all spectra obtained on top of Pb0.67Sn0.33Se(001) terraces (panels 1, 3, and 5) show the same features. In particular, similar to earlier observations on TCI surfaces (17, 22), the minimum visible about −75 meV marks the position of the Dirac point; the peaks L and L+ identify van Hove singularities, which signal the energy position of the two saddle points below and above the Dirac point, respectively.

Whereas this local DOS profile remains essentially unperturbed at the even step edge (position 2), completely different STS curves are observed at odd step edges—e.g., at position 4. Here, the spectrum is characterized by a strong peak at the Dirac point. More generally, its emergence is associated with a strong redistribution of the spectral weight over a relatively large energy range, as evidenced by the disappearance of the features associated with the van Hove singularities. This scenario is consistently also found for even and odd steps of higher order—i.e., triple and quadruple step edges (see fig. S2). Furthermore, the absence of scattering [similar to recent experiments on weak topological insulators (23)], the energy position locked at the Dirac point, and the association with translational symmetry breaking all point to a topological origin of these 1D channels.

To unequivocally prove that the emergence of these 1D edge states is linked to the existence of a nontrivial bulk band structure, we have performed measurements on crystals where the Sn concentration has been systematically changed, thereby spanning the entire range from trivial to topological surfaces (14) (Fig. 2). For all three concentrations—i.e., (A,B) x = 0, (C,D) x = 0.24, and (E,F) x = 0.33—the topographic images displayed in Fig. 2, A, C, and E, show one even and one odd step edge (see respective line profiles at the bottom of each panel). Irrespective of the step’s even- or oddness, no particular edge feature is visible in the dI/dU map of the topologically trivial material—i.e., of pure PbSe without any Sn (Fig. 2B). An edge state with a slightly enhanced dI/dU intensity develops once x = 0.24 (Fig. 2D)—i.e., just at the critical Sn concentration required for the transition between a normal and an inverted band gap, thereby leading to a weakly protected topological nontrivial state (14). By further increasing the amount of Sn to x = 0.33, the Pb1−xSnxSe compound is driven deeply into the topologically nontrivial regime. Correspondingly, we observe a strong enhancement of the local dI/dU signal at the position of the odd step edge (Fig. 2F).

Fig. 2 Sn concentration-dependent electronic properties of (Pb,Sn)Se.

Topography (top), dI/dU maps (bottom), and their corresponding profiles taken along the indicated line (bottom of each panel), measured on Pb1−xSnxSe crystals with different Sn content—i.e., (A and B) x = 0, (C and D) x = 0.24, and (E and F) x = 0.33—thereby spanning the range from trivial to topological surfaces. Step edges on the trivial compound (x = 0) carry no particular edge feature, irrespective of their even- or oddness. In contrast, a weak and strong enhancement of the local DOS is indicated by the high dI/dU signal measured at odd step edges for x = 0.24 and x = 0.33, respectively. Scan parameters: U = −310 mV, I = 30 pA (x = 0); U = −115 mV, I = 50 pA (x = 0.24); U = −70 mV, I = 100 pA (x = 0.33). T= 4.8 K.

To understand the accumulation of midgap states at odd step edges, it is best to think of them as electronic domain walls in the surface states (24). This effective 2D electronic surface state features four Dirac cones, and domain walls are created by interchanging the two atoms in the unit cell. The simplest 2D toy model displaying an even number of Dirac cones and a two-sublattice structure is the staggered flux model—i.e., a square lattice model restricted to nearest-neighbor hopping terms, in which neighboring plaquets are threaded by opposite magnetic fluxes +ϕ and −ϕ, with 0 < ϕ < π. We model a domain wall by interchanging the relative position of the two sublattices along the [11] direction (Fig. 3A). On the domain wall, the 2D Brillouin zone (BZ) of the staggered flux lattice is projected onto a 1D BZ, as indicated in green in Fig. 3B. Thereby, the two Dirac cones are not projected onto each other. We solve this system by Fourier transformation in the Embedded Image direction only, while keeping a real space index for the direction perpendicular to the wall. We obtain a chain model with alternating, real hopping amplitudes Embedded Image, with Embedded Image, and a domain wall (Fig. 3C). Except for the positions of the two Dirac cones in the 1D BZ, Embedded Image and Embedded Image, where Embedded Image, the 1D model is gapped. Because closing the chain in the presence of the domain wall requires an odd number of sites, whereas particle-hole symmetry demands that all states with Embedded Image are arranged in pairs with energies ±E, we necessarily obtain a midgap state at E = 0 localized around the domain wall. We obtain one such midgap state for each value of Embedded Image —i.e., one for each real space unit cell. The spectrum will hence schematically look as depicted in Fig. 3D.

Fig. 3 Theoretical analysis of the TCI step edge electronic structure.

(A) Staggered flux lattice with a domain wall in the [11] direction. (B) Brillouin zone including the projection onto the 1D domain wall. (C) Effective 1D hopping model describing the staggered flux lattice. (D) Schematic view of the resulting spectrum. (E) Tight-binding results of the electronic structure (momentum Embedded Image parallel to the edge) for a 3D TCI with [110]-oriented step edges of different height. The continuum of bulk states and surface contributions are shown as gray lines and red dots, respectively. Whereas strongly localized narrow bands connecting the two Dirac cones emerge at odd step edges, no such spectral weight is observed for even step edges. The two insets show the spin polarization in the presence of a single-atomic step edge: in-plane component perpendicular to the step edge (x, left) and the out-of-plane component (z, right) with the color (red and blue) indicating opposite spin directions.

To elevate this minimal model result to an accurate description of the experimental scenario observed in Fig. 1, we have performed tight-binding calculations adapted to reproduce the band structure of a topologically nontrivial crystalline insulator (13). The model is solved for a geometry with two step edges on the topmost surface; the bottom surface is left unperturbed. Figure 3E shows the calculated band structure of an infinite terrace (left panel), a single-, a double-, and a triple-step edge (right) as a function of Embedded Image, the momentum along the step edge. In close analogy to the staggered flux model described above, we consider the case of step edges parallel to the [110] direction. A discussion of the step edges along the [100] direction (as those shown in the sketch of Fig. 1, D and E), as well as further details of the model and the calculation, can be found in (20). In all cases, the electronic structure obtained for the step-edge-free surface corresponds to the results expected for a TCI surface—i.e., with two Dirac cones symmetrically shifted with respect to Embedded Image (gray lines). The size of the red dots is proportional to the degree of localization at the step edge.

Whereas even step edges give rise to two Dirac cones and a relatively wide open mouth-shaped feature connecting them, the odd-edge spectrum is markedly different. The states in between the Dirac cones have a narrow dispersion and therefore account for the overwhelming part of the step edge spectral weight (20). In contrast to the result of the staggered flux model presented above, these narrow states are not completely flat because the tight-binding model is not perfectly particle-hole symmetric. Furthermore, the midgap states do not extend over the whole Brillouin zone, as the Dirac cones are shifted away from high-symmetry points. For odd step edges, these states are laterally (i.e., perpendicular to the step) confined to a few nm, whereas in the z direction, they extend up to four lattice spacings. As shown in the left inset of Fig. 3E, our results for the single-atomic step edge correctly reproduce the spin-momentum locking of electronic states in the Dirac cone with an in-plane spin polarization. The right inset reveals that the two pairs of narrow, counter-propagating in-gap bands exhibit instead a high degree of out-of-plane spin polarization (maximum expectation value ≈ 0.7 for the operator 〈σz〉).

In comparison to 2D TIs, the Pb1−xSnxSe topological edge state reported here is expected to offer superior properties because it is protected by a bulk band gap that is up to 200 meV wide (14). We experimentally tested the robustness of the edge state to perturbations, such as hybridization with adjacent edge states, external magnetic field, or enhanced temperatures. For example, Fig. 4A shows a topographic STM image of a sample location with four odd step edges. Whereas the two right step edges run almost in parallel, thereby maintaining their separation, the two left step edges are inclined with respect to one another and eventually merge close to the bottom of the image. As qualitatively evidenced by the dI/dU map of Fig. 4B and quantitatively supported by the line sections plotted at the bottom of this panel, the edge state disappears as soon as the step-step separation decreases below the spatial extent of the edge state (25)—i.e., about 10 nm. We further analyze the response of these edge states to high magnetic fields B. Figure 4C reports a dI/dU map (top) and STS data acquired on an odd step edge (bottom) at B = 11 T; contrary to the quantum spin Hall state found in HgTe, the 1D TCI state investigated here is robust against time-reversal symmetry breaking perturbations. Finally, Fig. 4D shows that the edge state also persists at elevated temperatures (T = 80 K). Despite the reduced intensity evidenced by the STS spectrum, a well-defined 1D channel is still clearly present.

Fig. 4 Robustness of the edge state against perturbations.

(A) Topographic STM image showing a Pb0.67Sn0.33Se surface area exhibiting a total of four odd step edges. The separation between the two leftmost step edges continuously reduces from the top to the bottom of the image. (B) The dI/dU map reveals that the intensity of the edge state fades away as the edges approach below about 10 nm. Line sections along the lines indicated by numbers are plotted in the bottom panels. dI/dU maps taken around odd steps at (C) B = 11 T and (D) T = 80 K show the robustness of the 1D edge state against perturbations.

The observation of a distinct type of one-dimensional states at odd step edges of topological crystalline insulators with relatively wide bulk band gaps opens up opportunities for the use of topological materials for sensing and information processing purposes well beyond existing materials (4, 10, 11). Furthermore, the absence of scattering and the high degree of spin polarization observed in tight-binding calculations indicate that the 1D midgap state might be useful for spintronics applications. By patterning the step-and-terrace structure of TCI surfaces, this may allow for the creation of well-separated conductive channels with a width of only about 10 nm. This may lead to interconnections between functional units at ultrahigh packing densities. To fully explore whether the 1D midgap state found at odd TCI step edges display quantum conductance effects, further investigations by, for example, four-probe transport measurements, will be needed.

Supplementary Materials

Materials and Methods


Figs. S1 to S7

References (26, 27)

References and Notes

  1. See supplementary materials on Science Online.
  2. Acknowledgments: This research was supported by DFG (through SFB 1170 “ToCoTronics”; projects A02, B04, and C05) and by the Polish National Science Centre NCN grants 2014/15/B/ST3/03833 and 2012/07/B/ST3/03607. We further acknowledge support by the European Research Council (ERC) through ERC-StG-TOPOLECTRICS-Thomale-336012. D.D.S. and G.S. gratefully acknowledge the Gauss Centre for Supercomputing e.V. ( for funding this project by providing computing time on the GCS Supercomputer SuperMUC at Leibniz Supercomputing Centre (LRZ,

Stay Connected to Science

Navigate This Article