Abstract
The surface of a three-dimensional topological insulator (TI) hosts two-dimensional massless Dirac fermions (DFs), the gapless and spin-helical nature of which leads to their high transmission through surface defects or potential barriers. Here, we report the behaviors of topological surface states (TSS) in a triangular quantum corral (TQC) which, unlike a circular corral, is supposed to be totally transparent for DFs. By real-space mapping of the electronic structure of TQCs, both the trapping and detrapping behaviors of the TSS are observed. The selection rules are found to be governed by the geometry and spin texture of the constant energy contour of TSS upon the strong hexagonal warping in Bi2Te3. Our work indicates the extended nature of TSS and elucidates the selection rules of the trapping of TSS in the presence of a complicated surface state structure, giving insights into the effective engineering of DFs in TIs.
INTRODUCTION
The existence of topological surface states (TSS), with their low-energy behavior described by the Dirac equation rather than the Schrödinger equation, is ensured by the twisted bulk band topology resulting from the strong spin-orbital coupling of the material (1, 2). Unlike the massive electron systems having a forbidden gap, such as any trivial band insulator and even massive Dirac systems, the gapless Dirac fermion (DF) nature of TSS makes the potential barriers leaky (3). Furthermore, as long as time-reversal symmetry is respected, the backscattering of TSS is prohibited because of its helical nature (3–5), leading to the full transmission of normal incident electrons through potential barriers. This situation is more prominent for a family of three-dimensional (3D) topological insulators (TIs) (Bi2Se3, Bi2Te3, and Sb2Te3) (6–9), which have a single Dirac cone in the surface state structure and only one Fermi surface on the constant energy contour (CEC) within the bulk gap, further eliminating possible backscattering channels. Thus, the two fundamental properties of TSS, gaplessness and helical spin texture, may render any conventional quantum corral (10) highly leaky and complete trapping of TSS impossible, similar to the case of DFs in graphene (11).
For a 1D barrier on the TI surface, although reflection for normal incident electrons is forbidden, the reflection at oblique angles is possible, and even dominates the tunneling process at large incident angles, approaching one at grazing angles (12), which is governed by the spin texture of the CEC (13, 14). This makes the quasi-bounding of TI surface states possible via whispering-gallery modes as observed in a circular graphene pn junction (15–17) or, in the case of antimony, by constructing a Fabry-Pérot resonator (3, 12), the “infinity” of which along the 1D barriers helps to trap the TSS. However, antimony has a more complicated surface state structure (18), which facilitates backscattering and makes the trapping behavior of surface states more conventional. The abovementioned approaches both use the high reflection probability at near-grazing angles of DFs. Their intriguing attributes pique interest to construct different resonators, in which electrons can no longer keep bouncing at grazing angles, on TIs but with simple CECs.
Here, we successfully constructed such resonators as triangular quantum corrals (TQCs) on the surface of a 3D TI, Bi2Te3, which has a single Dirac cone and only one Fermi surface on the CEC in the surface state structure within the gap (6). Although the trapping of DFs by TQCs seems impossible because of the above argument, our experiments show well-defined interference patterns in some energy ranges of the TSS band. The trapping of Dirac electrons is found to be assisted by the hexagonal warping effect in the TSS, which in Bi2Te3 (19) is much stronger than that in Bi2Se3 (20). As shown below, the strong warping effect of TSS in Bi2Te3 yields energy-dependent trapping behaviors and selection rules of TSS in TQCs. Careful analysis of the quantum interference patterns of quasi-bound states yields the corresponding wave vectors of trapped TSS, through which two trapping mechanisms favoring momenta in different directions are uncovered.
RESULTS
MBE fabrication of TQCs and the observation of trapped TSS
Our TQCs are constructed by heteroepitaxially growing a Bi bilayer (BL) on the surface of the Bi2Te3 thin film (see Materials and Methods). The previous work (21) shows that at appropriate Bi coverages (typically 0.9 BL), Bi2Te3 vacancy islands surrounded by the Bi BL film can be formed. Taking an island with the side length L of ~35 nm, a close scrutiny uncovers quasi-equilateral triangular vacancy islands (Fig. 1A), bounded by
(A) Typical equilateral triangular Bi2Te3 vacancy island surrounded by Bi BLs (brighter regions) (42 nm × 42 nm, tunneling condition: 0.5 V, 0.1 nA). The Bi BL was grown on the Bi2Te3 film by MBE. The thickness of the Bi2Te3 film is 15 quintuple layers. The atomic resolution image in the white square shows the defect-free Te-terminated Bi2Te3 (16 nm × 16 nm, −0.4 V, 0.4 nA).
The Dirac energy of the pristine Bi2Te3 film locates at ~−170 meV (fig. S1) and is buried below the bulk valence band. Compared with pristine films, the dI/dV spectrum taken at the TQC center shows a downward shift of ~130 meV (blue curve in Fig. 1B), indicating the electron-doping effect of the Bi BL on the Bi2Te3 film (21). The Dirac energy of TSS in the uncovered region should be around −300 meV. The presence of a Bi BL changes the boundary condition of the Bi2Te3 film, and the Dirac energy is lifted to ~−210 meV. Therefore, there is a mismatch between the Dirac energies across the interface along
We note that the peaks shown in the spectrum taken at the TQC center (Fig. 1B) reveal the occurrence of quantized states, which is against our previous analysis that a TQC cannot trap massless DFs. The spatial patterns of these states were imaged by taking dI/dV maps in the energy range from −400 to 500 meV with a 25-meV interval (see fig. S2). Figure 2 shows the dI/dV maps at some representative energies that capture the main feature of trapping behaviors of TSS. The regular patterns can be seen below −250 meV or above −175 meV, and no pattern can be seen in between (except a bulk state–related pattern at −240 meV, which can be seen in Fig. 1B and more clearly in the Discussion). The patterns in the energy window from −175 to −100 meV are relatively weak and get much stronger above −75 meV.
Spatial dI/dV maps of TQC in Fig. 1A (42 nm × 42 nm) at different energies (−300, −275, −250, −225, −200, −175, −150, −125, −100, −75, 75, and 100 meV) (see fig. S2 for more images at other energies). The dashed triangles are guides for the eye, indicating the position of triangular barriers. The inset in each image is the simulated pattern of a certain index, which matches the interference pattern best. In this simulation, the spinor part of the wave function is omitted, which may only affect the overall intensity and does not contribute to the spatial variation of the pattern.
Modeling of the interference patterns
A simple particle-in-a-box model (23) is used and diagrammed in Fig. 3 (A and B) to simulate the above patterns and to elucidate the exact nature of these trapped states as well as the TSS selection rules in such a TQC. Provided the reflections are appreciable, an electron in a TQC with a side length of L can be possibly trapped by successive reflections at the triangular barrier walls, generating six wave vectors—k1, k2, k3, k4, k5, and k6—which have the c3v symmetry and are related to each other by in-plane 120° rotations and by reflections with respect to
(A) Schematic diagram of the scattering process of surface states inside the TQC. The thick (red) arrow indicates where the spatially dependent dI/dV spectra in Fig. 4 were taken. (B and C) Scattering process shown on the energy-dependent CECs of TSS with a strong warping effect. The dashed lines in (B) indicate the possible scattering vectors between wave vectors k1 to k6. q1, q2, and q3 denote three inequivalent reflections. kti and kvi are the wave vectors close to the tip and valley positions on the warped CECs, respectively. At high energies, the CECs are labeled with a spin texture exclusively for TSS with a strong warping term. Note that because of the warping effect, the spin-momentum helicity deviates from those at low energies, acquiring an out-of-plane spin component in the valley positions, as indicated by the crosses and dots in (C). The arrows with different lengths show the in-plane spin components with different magnitudes.
The simulated eigenmodes are shown in the insets of Fig. 2, which fit the interference patterns best at different energies. The patterns at and below −250 meV are found (see Discussion) to be caused by the band bending of the bulk states, except those at −275 and −300 meV, which are attributed to the (2, 1) and (3, 1) states, respectively. The patterns from −175 to −100 meV all have higher m indexes (m ≥ 2) but are relatively weak, much weaker than the spectral intensity of the Bi BL, especially the one at −175 meV, which is nearly indiscernible. At energies above −75 meV, the interference patterns get stronger and switch to low m indexes (m = 1).
In the diagram of Fig. 3C, the vector series corresponding to low or high m–indexed states are indicated on the energy-dependent CECs of TSS in Bi2Te3, as well as the complicated spin texture due to the strong hexagonal warping. For the (n, 1) states, the vectors kvi are close to valley positions in
DISCUSSION
The identification of trapped states below the highest valence subband
Figure 4A provides the spatially resolved energy levels of trapped states by taking the dI/dV spectra along the line across the TQC (thick arrow in Fig. 3A). Ec (−75 meV) denotes the conduction band minimum as identified by the dI/dV spectra (Fig. 1B). By using the pattern in Fig. 4A, the correction of the energies of the trapped states can be made. In addition, two trapped surface states (2, 1) and (3, 1) are identified. The spots at or below −240 meV, which correspond to the peaks in the dI/dV spectra of different sized TQCs (35 and 17 nm; see fig. S3), reflect the bulk valence subbands. The highest valence subband (denoted as Ev ~−240 meV) also shows as a peak (~−110 meV) in the dI/dV spectrum on the pristine Bi2Te3 film (Fig. 1B) but with an energy shift caused by different doping. These bulk-related patterns are identified to be the density variation of valence band states caused by band bending in the lateral direction near the steps; as noted above, the Bi BL has a strong electron-doping effect on Bi2Te3. The (2, 1) state centered at −280 meV in Fig. 4A, corresponding to the pattern at −275 meV, obeys a sinusoidally spatial intensity distribution, which is in contrast to other patterns below Ev. The (3, 1) pattern at ~−300 meV (~–320 meV after correction) does not have a node in the TQC center because of its weak intensity and the relatively stronger density variation of bulk states caused by band bending.
(A) Bias and spatially dependent dI/dV spectra along the line in Fig. 3A (tunneling condition: 0.2 V, 0.2 nA). The dashed line indicates the center position of the TQC. The sidebar indicates the energy regions as described in the main text having different trapping behaviors. (B) Blow-up maps from −300 to 100 meV. The red arrows indicate the localized states near the barriers. (C) Schematic of the band structure of Bi2Te3 thin films. The two short lines indicate the states (2, 1) and (3, 1) (see fig. S6 for details of the green dashed region). (D) Spectrum (Enm,
Seven different energy regions in the band structure of Bi2Te3
Thus, by a careful investigation of Fig. 4 (A and B), as well as the interference patterns at different energies (Fig. 2 and fig. S2), seven energy regions with distinct trapping behaviors or interference patterns are identified. The first region is below −320 meV, in which, except for the patterns caused by band bending, no trapped state can be observed. The second region is from −320 meV to Ev. Besides the patterns coming from the band-bending effect, two trapped states, (3, 1) at −320 meV and (2, 1) at −280 meV, are found. The (3, 1) state has a clover shape that cannot be explained by the band-bending effect. The pronounced peak of the (2, 1) state occurs in the dI/dV spectra of the TQC (L = 35 nm) but does not appear in that of a smaller TQC (L = 17 nm; see fig. S3) and that of the pristine film, indicating its trapped nature. The third region is from Ev to Ew (175 meV), in which no pattern appears (Figs. 2 and 4B). Here, Ew denotes approximately the energy above which the CEC of TSS becomes concave. The fourth region is from Ew to Ec (−75 meV), in which there are very weak patterns (Fig. 2). These weak patterns all match with the states of (n, m ≥ 2). The fifth region is from Ec to 75 meV. In this region, the patterns are stronger but seem mixed such that there are no clear spots or nodes in the spatial map of energy levels (Fig. 4B), although the dominating contribution comes from the (n, 1) modes. The pattern mixing is possible because of the level broadening caused by the finite lifetime of the trapped states and the unequal spaced energy levels of (n, m) states (some states with different n and different m may be very close). This means that the trapping probabilities of different channels are competing in such an energy range. The sixth region is from 75 to 400 meV, in which strong (n ≥ 10, 1) patterns are obvious at a glance. The indexing of trapped states can be double-checked by the symmetry of the (n, m) states. This can be done by noting that (3i + 2, 3j + 1) and (3i + 1, 3j + 2) states (i and j are nonnegative integers) have an antinode at the center of the TQC (dashed line in Fig. 4A). The last region is above 400 meV, where there are no more clear patterns. Note that the absence of trapped states in regions 1 and 7 is a result of the finite energy range of TSS.
Detailed nature of the trapped TSS in four regions (6, 5, 4, and 2)
On the basis of the above observation and discussion, the schematic band structure of our Bi2Te3 thin film is shown in Fig. 4C. Figure 4D gives a summary of the indexed states trapped by the TQC (L = 35 nm). The energies of the trapped states are all corrected by Fig. 4A, except those above 250 meV (half-filled squares). Here, we may focus on the details of the abovementioned trapping behaviors (in regions 6, 5, 4, and 2). For the (n ≥ 6, 1) states (regions 6 and 5; filled squares in Fig. 4D), the vectors kvi are close to valley positions in
For the higher m–indexed (n, m ≥ 2) states (region 4, hollow squares from −175 to −100 meV in Fig. 4D), the vectors kti are close to tip positions in
The (2, 1) and (3, 1) states below Ev (region 2) are also from TSS by using the argument for the states above Ec and also by noting that there is only one (2, 1) state despite the multiple valence subbands. However, the appearance of the (2, 1) state, especially its strong intensity, seems like a puzzle based on the arguments above, contradicting TRS.
The mechanism of three kinds of trapping behaviors
Our interpretation of these three kinds of trapping (in regions 2, 4, 5, and 6) behaviors is as follows. First, the Dirac energy is buried well below Ev. This makes the trapping of TSS below Ev (region 2) possible because there may exist two Fermi pockets with inverted spin textures similar to the antimony case (6, 18), especially for the (2, 1) state, as shown in fig. S6, where the backscattering is no longer forbidden. The deviation of
The other two trapping behaviors (in regions 4, 5, and 6) are found to be caused by the strong hexagonal warping in Bi2Te3, which leads to heavily deformed CECs and modified spin textures of TSS band (19). The energy-dependent CECs evolve from a circle near the Dirac point to a hexagon and then to a snowflake (6). In the case of a concave CEC above a certain energy, there may exist spin-polarized local modes along a line defect extended in the
Furthermore, the complexity of the spin texture (diagrammed in Fig. 3C) due to the strong hexagonal warping (19) may account for the appearance and the growing intensity of the (n ≥ 6, 1) states with energy. In the modified spin texture, the spin polarization near the
Despite all the abovementioned reflectivity-improving mechanisms that help to trap the DFs in TIs, the lifetime of the trapped states is still limited. The typical lifetime estimated from the peak widths (from 20 to 40 meV) ranges from 20 to 40 fs, which corresponds roughly to the reflectivity of TSS at the TQC edges from ~0.12 to ~0.35 (see fig. S8 for details). Here, only the states (n, 1) are discussed because the high-m states are too weak. The relatively low reflectivity indicates the high transmission of TSS at the TQC edges.
CONCLUSIONS
To conclude, despite the promised transparency of TQCs for TSS, three situations to trap the TSS are directly visualized in our work. One is the presence of double Fermi surfaces having opposite spin textures. For Bi2Te3, this region lies below the bulk valence band maximum. The other two take advantage of the strong hexagonal warping of TSS in Bi2Te3, which leads to the existence of possibly magnetic bound states and a more complicated spin texture. While the former opens the forbidden backscattering channels that facilitate the trapping of the high m–indexed (n, m ≥ 2) states of TQCs, the latter favors the trapping of the (n, 1) states, selecting the states near the tips (kti) and those close to the valleys (kvi) on the CEC, respectively. In the situation where all the aforementioned conditions are absent, the TQCs are completely transparent for the TSS, in agreement with the argument presented at the beginning. This region in Bi2Te3 is from Ev to Ev + 85 meV, where the CECs of TSS are supposed to be convex, constituting a truly forbidden region to trap the TSS. Note that the (3, 1) and (4, 1) states (Fig. 4D) above the Dirac point are missing. Similar experiments can be conducted on TIs without strong hexagonal warping (Bi2Se3 or Sb2Te3), which may have different trapping behaviors and may further validate our interpretation of current results. Thus, our findings elucidate the trapping and the detrapping behaviors of TSS in the presence of complicated CEC geometry and spin textures, giving insight into the manipulation of TSS, which is especially intriguing owing to their massless, relativistic and spin 1/2 DF nature.
MATERIALS AND METHODS
MBE growth and scanning tunneling microscopic characterizations
The main experiments were performed in ultrahigh vacuum (base pressure, ≤1.0 × 10−10 torr) using a commercial molecular beam epitaxy (MBE)–low temperature (LT) scanning tunneling microscopic (STM) combined system (Unisoku). Bi(111) BL films were heteroepitaxially grown on pristine Bi2Te3 films (thickness, ~15 quintuple layers) supported by graphitized SiC substrates. After growth, the samples were immediately transferred into in situ LT-STM so that the vacancy island shape could be frozen at 4.36 K to avoid Bi migration. The mechanically sharpened PtIr tips were used for STM after calibration on standard Pb/Si(111). Point scanning tunneling spectrum (STS or dI/dV) was performed with feedback loop open (lock-in bias modulation: 1 to 3 mVrms at 987.5 Hz), and dI/dV maps were obtained simultaneously with constant current images (8 to 10 mVrms at 2997 Hz). All the energies are shown with respect to the Fermi level (EF). The STM images were processed using WSxM software.
SUPPLEMENTARY MATERIALS
Supplementary material for this article is available at http://advances.sciencemag.org/cgi/content/full/5/5/eaaw3988/DC1
Fig. S1. Band structure of Bi/Bi2Te3 on (111) surface.
Fig. S2. Successive dI/dV maps of the TQC (35 nm) at energies with an interval of 25 meV.
Fig. S3. dI/dV spectra on TQCs of different sizes.
Fig. S4. dI/dV maps on a smaller TQC (17 nm).
Fig. S5. dI/dV spectra of Fig. 4A in the line plot form.
Fig. S6. A zoom-in view of the schematic band structure in Fig. 4C.
Fig. S7. Section profiles of patterns at 50, 75, and 100 meV along the lines indicated in fig. S2.
Fig. S8. A particular situation adopted to give a rough estimation of the trapped states’ lifetime in which the path of the electron in a TQC roughly follows a regular pattern.
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).