Robustness of topological order and formation of quantum well states in topological insulators exposed to ambient environment
See allHide authors and affiliations
Edited by* E. Ward Plummer, Louisiana State University, Baton Rouge, LA, and approved January 18, 2012 (received for review September 24, 2011)

Abstract
The physical property investigation (like transport measurements) and ultimate application of the topological insulators usually involve surfaces that are exposed to ambient environment (1 atm and room temperature). One critical issue is how the topological surface state will behave under such ambient conditions. We report high resolution angle-resolved photoemission measurements to directly probe the surface state of the prototypical topological insulators, Bi2Se3 and Bi2Te3, upon exposing to various environments. We find that the topological order is robust even when the surface is exposed to air at room temperature. However, the surface state is strongly modified after such an exposure. Particularly, we have observed the formation of two-dimensional quantum well states near the exposed surface of the topological insulators. These findings provide key information in understanding the surface properties of the topological insulators under ambient environment and in engineering the topological surface state for applications.
The topological insulators represent a novel state of matter where the bulk is insulating but the surface is metallic, which is expected to be robust due to topological protection (1⇓⇓⇓–5). The topological surface state exhibits unique electronic structure and spin texture that provide a venue not only to explore novel quantum phenomena in fundamental physics (6⇓⇓⇓–10) but also to show potential applications in spintronics and quantum computing (2,5,11). The angle-resolved photoemission spectroscopy (ARPES) is a powerful experimental tool to directly identify and characterize topological insulators (12). A number of three-dimensional topological insulators have been theoretically predicted and experimentally identified by ARPES (13⇓⇓⇓⇓⇓⇓⇓–21); some of their peculiar properties have been revealed by scanning tunneling microscopy (STM) (22⇓⇓⇓–26). The application of the topological surface states depends on the surface engineering that can be manipulated by incorporation of nonmagnetic (27⇓⇓⇓–31) or magnetic (27, 28, 31⇓–33) impurities or gas adsorptions (27, 33⇓–35). While the ARPES and STM measurements usually involve the fresh surface obtained by cleaving samples in situ under ultrahigh vacuum, for the transport and optical techniques, which are widely used to investigate the intrinsic quantum behaviors of the topological surface state (36⇓⇓⇓–40), and particularly the ultimate applications of the topological insulators, the surface is usually exposed to ambient conditions (1 atm air and room temperature) or some gas protection environment. It is therefore crucial to investigate whether the topological order can survive under the ambient conditions and, furthermore, whether and how the surface state may be modified after such exposures.
Results and Discussion
We start by first looking at the electronic structure of the prototypical topological insulators Bi2(Se,Te)3 under ultrahigh vacuum. The Fermi surface and the band structure of the Bi2(Se3-xTex) topological insulators depend sensitively on the composition, x, as shown in Fig. 1. The single crystal samples here were all cleaved in situ and measured at 30 K in an ultrahigh vacuum (UHV) chamber with a base pressure better than 5 × 10-11 torr. For Bi2Se3, a clear Dirac cone appears near -0.36 eV (Fig. 1 D and E); the corresponding Fermi surface (Fig. 1A) is nearly circular but with a clear hexagon shape in the measured data (41). It is apparently of n type because the Fermi level intersects with the bulk conduction band. On the other hand, the Dirac cone of the Bi2Te3 sample lies near -0.08 eV (Fig. 1 H and I), much closer to the Fermi level than that reported before (-0.34 eV in ref. 16). The corresponding Fermi surface (Fig. 1C) becomes rather small, accompanied by the appearance of six petal-like bulk Fermi surface sheets. These results indicate that our Bi2Te3 sample is of p type because the Fermi level intersects the bulk valence band along the direction. This is also consistent with the positive Hall coefficient measured on the same Bi2Te3 sample (42). This difference of the Fermi surface topology and the location of the Dirac cone from others (16) may be attributed to the different carrier concentration in Bi2Te3 due to different sample preparation conditions. In our Bi2(Se3-xTex) samples, we have seen a crossover from n-type Bi2Se3 to p-type Bi2Te3. In order to eliminate the interference of the bulk bands on the surface state near the Fermi level, we fine tuned the composition x in Bi2(Se3-xTex) and found that, for x = 2.6, nearly no spectral weight can be discerned from the bulk conduction band, as seen from both the Fermi surface (Fig. 1B) and the band structure (Fig. 1 F and G). A slight substitution of Te by Se in Bi2(Se0.4Te2.6) causes a dramatic drop of the Dirac point to -0.31 eV (Fig. 1 F and G) and an obvious hexagon-shaped Fermi surface (Fig. 1B). It is interesting to note that the hexagon shape of Bi2(Se0.4Te2.6) (Fig. 1B) is rather pronounced, although its Fermi surface size is smaller than that of Bi2Se3 (Fig. 1A). The hexagonally shaped Fermi surface observed in the topological surface states reflects the hybridization of surface electronic states with the bulk states and can be theoretically explained by considering the higher order terms in the k·p Hamiltonian (43).
Fermi surface and band structure of Bi2(Se3-xTex) (x = 0, 2.6, 3) topological insulators cleaved in situ and measured at 30 K in ultrahigh vacuum. (A–C) Fermi surface of Bi2Se3, Bi2(Se0.4Te2.6), and Bi2Te3, respectively. The Fermi surface here, and in other figures below, are original data without involving artificial symmetrization. The band structures along two high symmetry lines and
are shown in D and E for Bi2Se3, in F and G for Bi2(Se0.4Te2.6), and H and I for Bi2Te3.
In order to directly examine how the topological surface state behaves under ambient conditions in the topological insulators, we carried out our ARPES measurements in different ways. (1). We first cleaved the sample in situ and performed ARPES measurement in the ultrahigh vacuum (UHV) chamber. The sample was then pulled out to another chamber filled with 1 atm N2 gas, exposed for about 5 min, before transferring back to the UHV chamber to do ARPES measurements (2). We cleaved and measured the sample in the UHV chamber and then pulled the sample out to air for 5 min before transferring back to the UHV chamber for the ARPES measurements; (3). We cleaved the sample in air and then transferred it to the UHV chamber to do the ARPES measurement. Our measurements show that the above procedures of exposure to air or N2 produce similar and reproducible results for a given sample.
The surface exposure of the topological insulators to air or N2 gives rise to a dramatic alteration of the surface state, as shown in Figs. 2⇓–4, for Bi2Se3, Bi2(Se0.4Te2.6), and Bi2Te3, respectively, when compared with those for the fresh surface (Fig. 1). The first obvious change is the shifting of the Dirac cone position relative to the Fermi level. For Bi2Se3, Bi2(Se0.4Te2.6), and Bi2Te3, it shifts from the original -0.36 eV (Fig. 1 D and E), -0.31 eV (Fig. 1 F and G), -0.08 eV (Fig. 1 H and I) for the fresh surface to -0.48 eV (Fig. 2B), -0.40 eV (Fig. 3 A and B), and -0.28 eV (Fig. 4 C–F) at 30 K for the exposed surface, respectively. In all these cases, the shift of the Dirac cone to a larger binding energy indicates an additional doping of electrons into the surface state. The exposure also gives rise to a dramatic change of the surface-state Fermi surface. For Bi2Se3, in addition to a slight Fermi surface size increase, an obvious change occurs in the Fermi surface shape that the hexagon shape becomes much more pronounced in the exposed surface (Fig. 2D) than that in the fresh sample (Fig. 1A). For Bi2(Se0.4Te2.6), one clearly observes the much-enhanced warping effect in the exposed surface (Fig. 3C) when compared with the nearly standard hexagon in the fresh surface (Fig. 1B). The most dramatic change occurs for Bi2Te3 where not only the Fermi surface size increases significantly but also the warping effect in the exposed surface (Fig. 4I) becomes much stronger. Overall, the exposure causes the lowering of the Dirac cone position, an increase of the surface Fermi surface size, and an obvious enhancement of the Fermi surface warping effect in the Bi2(Se3-xTex) system. A careful comparison of the energy bands and Fermi surface between the fresh and exposed surfaces indicates that, if we take the Dirac cone energy as a common reference energy, the Fermi surface and bands nearly overlap with each other (see Fig. S1). These indicate that the Fermi surface change in the exposed samples is mainly due to chemical potential shift, not from the Fermi surface deformation.
Fermi surface and band structure of Bi2Se3 cleaved in air and measured in the ultrahigh vacuum (UHV) chamber. (A) Band structure of the fresh Bi2Se3 cleaved and measured in the UHV chamber at 30 K along direction. (B) Band structure of Bi2Se3 cleaved in air and measured in UHV at 30 K along
direction. (C) Band structure of Bi2Se3 cleaved in air and measured in UHV at 300 K along
direction. (D and E) Fermi surface of Bi2Se3 cleaved in air and measured in UHV at 30 K and 300 K, respectively. Black dashed lines in B and C mark the parabolic bands above the Dirac point from the two-dimensional electron gas.
Emergence of quantum well states in Bi2Te2.6Se0.4 after exposing to N2. (A and B) Band structure measured at 30 K along and
, respectively. Black dashed lines in B mark the quantum well states formed in the bulk conduction band (BCB) above the Dirac point. (C) The corresponding Fermi surface. It shows threefold symmetry where three corners of M points are strong while the other three are weak. This is also in agreement with the asymmetric band structure in Fig. 3B. (D) Schematic band structure showing the possible formation of the quantum well states near the sample surface in the bulk conduction band. The blue dotted lines between the bulk valence band (BVB) and bulk conduction band (BCB) represent the topological surface states while the blue solid lines represent quantum well states.
Persistence of topological surface state and formation of quantum well states in Bi2Te3 after exposure to N2 or air. The sample was first cleaved and measured in UHV at 30 K. (A and B) The corresponding band structure along the and
directions. The sample was then pulled out from the UHV chamber and exposed to N2 at 1 atm for 5 min before transferring back into UHV chamber for the ARPES measurement. (C and D) The band structure of the N2-exposed sample along the
and
directions. The black dashed lines in C illustrate the quantum well states formed in the bulk valence band below the Dirac point. The sample was then pulled out again and exposed to air for 5 min before putting back in vacuum for ARPES measurement. (E and F) The band structure of the air-exposed sample at 30 K along the
and
directions. (G and H) The measurements at 300 K, and I and J show their corresponding second-derivative images in order to highlight the bands. (K) Fermi surface of N2-exposed sample. (L) First principle calculation of the band structure of Bi2Te3 slab with seven quintuple layers.
The topological order in the Bi2(Se3-xTex) topological insulators is robust even when the surface is exposed to ambient conditions, in spite of all the alterations mentioned above. One clearly observes the persistence of the Dirac cone in the exposed surface as in Bi2Se3 (Fig. 2 B and C), in Bi2(Se0.4Te2.6) (Fig. 3 A and B), and Bi2Te3 (Fig. 4 C–H). This is particularly the case for the surface exposed to air and measured at room temperature (Fig. 2C for Bi2Se3 and Fig. 4 G and H for Bi2Te3). On the other hand, after the exposure, although the signal of the surface state gets weaker for Bi2Se3 (Fig. 2), it remains rather strong for Bi2(Se0.4Te2.6) (Fig. 3) and Bi2Te3 (Fig. 4). This is in stark contrast to the conventional trivial surface state where minor surface contamination will cause the extinction of the surface state (44). The robustness of the topological order to Coulomb, magnetic, and disorder perturbations has been reported before (33, 34). Our present observations directly demonstrate the robustness of the topological order against absorption and thermal process under ambient conditions, presumably due to the protection of the time-reversal symmetry (3, 4).
The surface exposure to air or N2 in the Bi2(Se3-xTex) topological insulators produces two-dimensional electronic states near the surface. In Bi2Se3, the exposure gives rise to additional parabolic bands, as schematically marked by the dashed line in Fig. 2 B and C. Correspondingly, this leads to additional Fermi surface sheet(s) inside the regular topological surface state (Fig. 2 D and E). In Bi2(Se0.4Te2.6), this effect gets more pronounced and the newly emerged bulk conduction band splits into several discrete bands, as marked by the dashed lines in Fig. 3B. While the band quantization effect occurs in the bulk conduction band in Bi2(Se0.4Te2.6), it shows up in the valence band in the exposed Bi2Te3 surface, as shown in Fig. 4 C–H, where one can see a couple of discrete M-shaped bands. The quantized bands are obvious at low temperature and get slightly smeared out when the temperature rises to room temperature (Fig. 4 G and H). One may wonder whether these two-dimensional electronic states might come from the formation of a different phase on the surface due to the exposure. We believe this is unlikely because, as shown in Figs. 3 and 4, the two-dimensional electronic states are rather different although the composition of Bi2(Se0.4Te2.6) is only slightly different from Bi2Te3.
It is interesting to note that the Dirac structure shows an obvious change with temperature. As shown in Fig. 2 B and C, as the measurement temperature increases from 30 K to 300 K, the Dirac cone location of Bi2Se3 shifts upward from -0.48 eV to -0.38 eV. For Bi2Te3, the Dirac cone also shifts to lower binding energy upon increasing temperature from -0.28 eV at 30 K (Fig. 4E) to -0.25 eV at 300 K (Fig. 4G). On the other hand, the two-dimensional quantum well states show little change with temperature, as shown in Fig. 2 B and C for Bi2Se3. The upward energy shift of the Dirac cone with increasing temperature is possibly due to desorption process on the sample surface which results in electron removal from the surface.
The formation of the split bands in the exposed surface of the topological insulators is reminiscent of the quantum well states observed in the quantum confined systems (45) and in some topological insulators (29⇓–31, 35, 46, 47). There are a couple of possibilities that the quantum well states may be formed. One usual way is due to the band-bending effect. The surface exposure to air or N2 causes an electron transfer to the surface of the topological insulators. The accumulation of these additional electrons near the surface would lead to a downward bending of the bulk bands near the surface region, as schematically shown in Fig. 3D, resulting in a V-shaped potential well where the bulk conduction band of electrons can be confined. This picture, as proposed before (29, 46), becomes questionable to explain the present observation.
Generally speaking, this is because the valence band top has a hole-like component, so the downward band-bending that acts as a quantum well potential for electrons will no longer be a potential for holes. Specifically, Bianchi et al. (29) argued that, when the total bandwidth of the valence band is narrower than the band-bending depth, it is possible to form quantum well states for the valence bands. This scenario does not work for our Bi2Te3 case for a couple of reasons. First, although the M-shaped valence band has a V-shape in the middle that acts like electrons and could get quantum well states, the two wings of the bands remain hole-like and they should not generate quantum well states. This is not consistent with the experimental observations that the entire M-shaped band exhibits quantum well states. Second, as pointed out in ref. 29, in order for the band-bending picture to work for the valence band, a necessary condition is that the total bandwidth of the valence band must be smaller than the band-bending depth. This condition is not satisfied in Bi2Te3. As seen in Fig. S2, our measured band width of Bi2Te3 along the Γ-K direction is 300 meV (Fig. S2A). The band structure calculations give a band width along the Γ-K direction of 350 meV (Fig. S2B) (16). Using a similar procedure as in ref. 29, the band-bending depth is 227 meV as determined from the position difference of the Dirac point between the freshly cleaved sample (Fig. 4A) and the exposed sample (Fig. 4B). Therefore, the total bandwidth of Bi2Te3 valence band is obviously larger than that of the band-bending depth, which does not satisfy the necessary requirement proposed in ref. 29. The observation of quantum well states in the valence band of Bi2Te3 cannot be explained by the picture proposed in ref. 29. Therefore, the band-bending is not a general picture to explain the formation of the quantum well states in all the samples on the same footing.
An alternative scenario is the expansion of van der Waals spacings in between the quintuple layers (QLs) caused by the intercalation of gases (48). The observation of multiple split bands with different spacings would ask for multiple van der Waals gaps with different expansions. Whether and how these can be realized in the exposed surface remains to be investigated. We note that our observations of multiple split bands are similar to those seen in the ultrathin films of Bi2Se3 (49) and Bi2Te3 (50). From our first principle band structure calculations on Bi2Te3 with different number of quintuple layers, we also find that a detached slab with a thickness of seven quintuple layers can give a rather consistent description (Fig. 4L) of our observed results in terms of the quantitative spacings between the three resolved bands (VB0, VB1, and VB2 bands as marked in Fig. 4 C and L). In addition, the distance between the conduction band bottom (CB0 band in Fig. 4 I and L) and the first valence subband bottom (VB0 band in Fig. 4 I and L) is rather consistent between the measured and calculated results. These seem to suggest that a “confined surface slab” with nearly seven quintuple layers may be formed after the exposure that acts more or less independently from the bulk. More work needs to be done to further investigate whether such a confined surface slab can be thermodynamically stable. Overall, the formation of the two-dimensional quantum well states is a general phenomenon for the exposed surface of the Bi2(Se3-xTex) topological insulators; the effect depends sensitively on the composition x of the samples, which may facilitate manipulation of these quantum well states.
The present work has significant implications on the fundamental study and ultimate applications of the topological insulators. Many experimental measurements, such as some transport measurements, involve samples exposed to ambient conditions. The practical applications may involve sample surface either exposed to ambient condition or in contact with other magnetic or superconducting materials. On the one hand, the robustness of the topological order under ambient conditions sends a good signal for these experimental characterization and practical utilizations. The formation of the quantum well states may give rise to new phenomena to be studied and utilized. The sensitivity of the surface state to the Bi2(Se3-xTex) composition provides a handle to manipulate these quantum states. On the other hand, the strong modification of the electronic structure and the formation of additional quantum well states in the exposed surface have to be considered seriously in interpreting experimental data and in surface engineering. The observed change of resistivity and Hall coefficient with time can be understood as a result of electron doping on the air-exposed surface (51). It is critical to realize beforehand that the surface under study or to be utilized may exhibit totally different behaviors as those from the fresh surface cleaved in ultrahigh vacuum. In addition to the alteration of electronic states upon exposure, the transport properties of the topological surface state may be further complicated by the formation of quantum well states. In this sense, the transport measurements need to be checked because no considerations were made before on the formation of quantum well states that may affect transport analysis (36⇓⇓–39).
Methods
Crystal Growth Methods.
Single crystals of Bi2(Se3-xTex) (x = 0, 2.6, and 3) were grown by the self-flux method. Bismuth, selenium, and tellurium powders were weighed according to the stoichiometric Bi2(Se3-xTex) (x = 0, 2.6, and 3) composition. After mixing thoroughly, the powder was placed in alumina crucibles and sealed in a quartz tube under vacuum. The materials were heated to 1,000 °C, held for 12 h to obtain a high degree of mixing, and then slowly cooled down to 500 °C over 100 h before cooling to room temperature. Single crystals of several millimeters in size were obtained. The crystal structure of the resulting crystals was examined by use of a rotating anode X-ray diffractometer with Cu Kα radiation (λ = 1.5418 Å). The chemical composition of the crystals was analyzed by the energy dispersive X-ray spectroscopy (EDAX) and the induction-coupled plasma atomic emission spectroscopy (ICP-AES). The resistivity of the crystals was measured by the standard four-probe method.
Laser-ARPES Methods.
The angle-resolved photoemission measurements were carried out on our vacuum ultraviolet (VUV) laser-based angle-resolved photoemission system (52). The photon energy of the laser is 6.994 eV with a bandwidth of 0.26 meV. The energy resolution of the electron energy analyzer (Scienta R4000) is set at 1 meV, giving rise to an overall energy resolution of approximately 1 meV, which is significantly improved from 10 ∼ 15 meV from regular synchrotron radiation systems (15, 16). The angular resolution is approximately 0.3°, corresponding to a momentum resolution of approximately 0.004 Å-1 at the photon energy of 6.994 eV, more than twice improved from 0.009 Å-1 at a regular photon energy of 21.2 eV for the same angular resolution. Our superior instrumental resolution of laser ARPES has made the measured features of topological insulators in this work much sharper. The Fermi level is referenced by measuring on a clean polycrystalline gold that is electrically connected to the sample. The samples were all measured in vacuum with a base pressure better than 5 × 10-11 torr.
Acknowledgments
We thank Prof. Xianhui Chen for providing us samples at the initial stage of the project and Prof. Liling Sun and Prof. Zhong-xian Zhao for their help in the characterization of the samples. This work is supported by the National Natural Science Foundation of China (91021006) and the Ministry of Science and Technology of China (973 program 2011CB921703).
Footnotes
- ↵1To whom correspondence should be addressed. E-mail: XJZhou{at}aphy.iphy.ac.cn.
Author contributions: Chaoyu Chen and X.J.Z. designed research; Chaoyu Chen and H.W. performed research; S.H., W.Z., L.Z., H.L., X.J., D.M., S.L., J.H., Y.P., Y.F., Z. Xie, G.L., X. Dong, J.Z., X.W., Q.P., Z.W., S.Z., F.Y., Chuangtian Chen, and Z. Xu contributed new reagents/analytic tools; Chaoyu Chen, X. Dai, Z.F., and X.J.Z. analyzed data; and Chaoyu Chen and X.J.Z. wrote the paper.
The authors declare no conflict of interest.
*This Direct Submission article had a prearranged editor.
This article contains supporting information online at www.pnas.org/lookup/suppl/doi:10.1073/pnas.1115555109/-/DCSupplemental.
Freely available online through the PNAS open access option.
References
- ↵
- ↵
- Qi XL,
- Zhang SC
- ↵
- ↵
- ↵
- ↵
- Qi X-L,
- Hughes TL,
- Zhang S-C
- ↵
- ↵
- Qi X-L,
- Li R,
- Zang J,
- Zhang S-C
- ↵
- ↵
- Yu R,
- et al.
- ↵
- ↵
- Hasan MZ
- ↵
- ↵
- ↵
- ↵
- Chen YL,
- et al.
- ↵
- ↵
- ↵
- ↵
- ↵
- ↵
- ↵
- ↵
- ↵
- ↵
- Hanaguri T,
- et al.
- ↵
- Wray LA,
- et al.
- ↵
- ↵
- ↵
- ↵
- Pan Z-H,
- et al.
- ↵
- Chen YL,
- et al.
- ↵
- ↵
- ↵
- ↵
- Qu D-X,
- Hor YS,
- Xiong J,
- Cava RJ,
- Ong NP
- ↵
- ↵
- ↵
- ↵
- ↵
- ↵
- ↵
- ↵
- ↵
- Speer NJ,
- Tang S-J,
- Miller T,
- Chiang T-C
- ↵
- Bianchi M,
- et al.
- ↵
- ↵
- Eremeev SV,
- et al.
- ↵
- ↵
- ↵
- ↵
Citation Manager Formats
Article Classifications
- Physical Sciences
- Applied Physical Sciences