Abstract
New threedimensional (3D) topological phases can emerge in superlattices containing constituents of known twodimensional topologies. Here we demonstrate that stoichiometric Bi_{1}Te_{1}, which is a natural superlattice of alternating two Bi_{2}Te_{3} quintuple layers and one Bi bilayer, is a dual 3D topological insulator where a weak topological insulator phase and topological crystalline insulator phase appear simultaneously. By density functional theory, we find indices (0;001) and a nonzero mirror Chern number. We have synthesized Bi_{1}Te_{1} by molecular beam epitaxy and found evidence for its topological crystalline and weak topological character by spin and angleresolved photoemission spectroscopy. The dual topology opens the possibility to gap the differently protected metallic surface states on different surfaces independently by breaking the respective symmetries, for example, by magnetic field on one surface and by strain on another surface.
Introduction
Topological insulators (TIs) are bulk insulating materials that exhibit metallic conductivity on their boundaries via electronic edge (in twodimensional (2D) TIs) or surface states (in threedimensional (3D) TIs), which are guaranteed by the topology of the bulk band structure^{1,2}. Electrons in these boundary states are spin polarized. Their spin and momentum are locked to each other by spin–orbit coupling, creating helical spin textures, which make TIs highly attractive for spintronic applications^{3}. One of the most favourable aspects of 3D TIs is the fact that their surface inevitably hosts these metallic surface states as long as the symmetry defining the topological index is not broken^{4,5}. In a strong TI (STI), timereversal symmetry protects these states on all surfaces. Weak TIs (WTIs), on the other hand, display protected metallicity only at surfaces with a certain orientation, while other surfaces do not contain topologically protected surface states. The latter can be understood in a simple picture, where a stack of 2D TIs forms a WTI with metallic surface states inherited from the edge states of the 2D TI but with an insulating surface plane (the dark side) normal to the stacking direction. Finally, in topological crystalline insulators (TCIs), where the symmetry with respect to a mirror plane defines the topology, metallic surface states can be found on surfaces perpendicular to these mirror planes^{6,7}.
Bi_{2}Te_{3} was the first material predicted to be both a STI and a TCI. Since it exhibits two topological properties, it was termed a dual TI^{8}. Such a combination opens the possibility that controlled symmetry breaking would destroy certain surface states while keeping others intact. For example, one could imagine a material that is both a WTI and a TCI and has all surfaces covered by metallic surface states, that is, the mirror plane of the TCI is normal to the dark side of the WTI. Then, a magnetic field would destroy the topological protection of the states caused by the WTI character, while the mirrorsymmetryprotected states remain intact such that the dark surface is metallic. Likewise, small structural distortions can break the mirror symmetry without affecting the surface states arising from timereversal symmetry, thus rendering only the dark side insulating.
In the search for such a material, we start from Bi_{2}Te_{3} with the properties mentioned above and from the Bi bilayer (BL) which is known to be a 2D TI^{9,10,11}. It is possible to produce natural superlattices [Bi_{2}]_{x}[Bi_{2}Te_{3}]_{y} from hexagonal, metallic BLs and semiconducting Bi_{2}Te_{3} quintuple layers (QLs) in a wide range of x and y^{12,13,14,15}. Within this series Bi_{2}Te_{3} consists of QL building blocks only, while the unit cell of Bi_{1}Te_{1} exhibits a stacking sequence of a single BL interleaved with two subsequent QLs. The size of the unit cell along the stacking direction, that is, the c lattice constant, varies quite severely among the different stable compounds, which makes them easily distinguishable in a diffraction experiment. Recently, a similar superlattice Bi_{4}Se_{3} (that is, x=y=1) was investigated in some detail and characterized as a topological semimetal^{16}. For practical applications, however, an insulating bulk material is preferable.
In this article, we identify the stoichiometric natural superlattice Bi_{1}Te_{1} (that is, x=1, y=2) as a semiconductor with a small bandgap of about 0.1 eV with the desired properties: Bi_{1}Te_{1} is both a WTI and a TCI, and hence a novel type of dual TI with the favourable property of suitable surfaces as described above. Our density functional theory (DFT) calculations predict a class of (0;001) and a mirror Chern number . We find two characteristic surface states of the TCI on the (0001) surface, regardless of the surface termination. A similar situation has been reported recently for Bi_{2}TeI, theoretically^{17}. We demonstrate that Bi_{1}Te_{1} can be grown in the form of highquality thin films on Si(111) by molecular beam epitaxy (MBE). Its layered structure is confirmed by scanning transmission electron microscopy (STEM) and Xray diffraction (XRD), depicting a repeated stacking sequence of 2QL of Bi_{2}Te_{3} and a single Bi BL. We investigate the electronic structure of Bi_{1}Te_{1} by means of spin and angleresolved photoemission (spinARPES) on the (0001) surface. Our spectra taken along nonhighsymmetry lines reveal band crossings away from timereversal invariant momenta (TRIM) points that can be associated with surface states protected by mirror symmetry and the TCI character of Bi_{1}Te_{1}. SpinARPES reveals that the surface states surrounding Γ and being close to E_{F} exhibit a nearly vanishing spin polarization in contrast to the directly compared timereversal symmetrydriven TSS in Bi_{2}Te_{3} (x=0, y=3, y is 3 here because the unit cell consists of three QLs), which is in line with WTI character with the dark surface being (0001).
Results
Ab initio calculations
Figure 1a depicts a schematic model of the crystal structure of Bi_{1}Te_{1}, indicating Bi (green) and Te (orange) atoms as well as one unit cell defined by the lattice constant c along the stacking direction. The separation of the layered structure into QLs and Bi BLs is marked. The bulk band structure of Bi_{1}Te_{1} in the relaxed structural geometry is presented in Fig. 1b. Spin–orbit coupling is included in this calculation and the colour represents the localization of the electronic states at the BL (green) or at the QL (red). As one can see, there are no states at the Fermi level, E_{F}, reflecting the insulating character with an energy gap of 73 meV. The states around the Fermi level alternate between BL and QLrelated, where the highest occupied levels at the timereversal invariant momenta Γ and A stem from QLs (red), while the lowest unoccupied states originate from the BL (green), and vice versa for the M and L points. Since the crystal possesses spatial inversion symmetry, the parity of the states can be calculated and the topological index can be deduced according to ref. 18 based on the product of the parities of all occupied bands at the eight TRIMs, that is, one at Γ, one at A, three at M and three at L. The result is shown in Fig. 1c for the corresponding TRIMs, Γ and A have parity products of +1 (red ‘+’) while M and L have −1 (blue ‘−’), leading to a topological invariant . Therefore, Bi_{1}Te_{1} is a weak TI with the (0001) surface, which is perpendicular to the stacking direction, being the dark surface and being free of timereversal symmetryprotected surface states.
It is tempting to relate the WTI property to the fact that both the BL and the 2QLs Bi_{2}Te_{3} are 2D TIs such that the WTI results from a simple stacking of 2D TIs in the cdirection. However, our band structure calculations in Fig. 1d, which introduce artificially expanded distances between the BL and the QLs, show a more complex scenario: if the BL is sufficiently separated from the 2QLs, the states can be decomposed into contributions from the two components (green=Bi BL and red=2QLs, respectively). However, due to charge transfer, the inverted gap of the BL is shifted above the Fermi level and, accordingly, some of the 2QL Bi_{2}Te_{3} conduction band states are below E_{F}. Only the hybridization of the BL states with the QL states opens up the gap that leads to the insulating bulk structure in Fig. 1b, as can be deduced from the changing colour of the bands along the kdirections. Nevertheless, the topological character of the stacked film remains nontrivial. A similar complexity is also found for the first confirmed, stacked weak TI Bi_{14}Rh_{3}I_{9} (refs 19, 20, 21).
Next, we examined the surface band structure of the (0001) surface, which heavily depends on the precise surface termination. Due to the layered crystal structure and the weak van der Waals bonds between the subsequent building blocks, there exist three natural cleavage planes and thus surface terminations, that is, 1Bi BL, 1QL and 2QL. Figure 2 depicts the spinresolved surface band structure for the latter two terminations. In all cases an even number of Fermi level crossings is found in direction (see Supplementary Notes II and III for all terminations and convergence with film thickness). We observed that, although there are surface states, the band structure is compatible with the fact that we look at the dark surface of a WTI phase. Nevertheless, it is remarkable that the bands along show a band crossing for all possible terminations, reminiscent of Diraclike cones observed in TCIs. Additional evidence that this crossing is protected by a mirror symmetry in the crystal comes from the observation that the crossing is lifted when the surface atoms are displaced in direction, breaking this symmetry in one of the directions (Fig. 2c).
To finally confirm that Bi_{1}Te_{1} is a TCI, we determined the mirror Chern number of the bulk phase. In the (k_{x},k_{z}) plane in reciprocal space (Fig. 1c) all Bloch states can be distinguished by their eigenvalues with respect to a mirror operation in the plane. To calculate their corresponding Berry phases as well as the Chern numbers, we construct a tightbinding Hamiltonian based on the maximally localized Wannier functions^{22}. The Chern numbers of all occupied bands for the opposite mirror eigenvalues +i and −i are n_{+i}=2 and n_{−i}=−2, respectively, and therefore the mirror Chern number (ref. 6), given as , is , confirming the fact that Bi_{1}Te_{1} is a TCI. The modulus of shows that we have to expect two linear crossings along the highsymmetry line formed by the mirror and the surfaceplane^{23}, the sign determines the spin orientation^{6} shown in the Supplementary Note II. The mirror Chern number is incompatible with STI phase since it produces even number of crossings. In this context it might be interesting to note that the dual TI Bi_{2}Te_{3} has , and consequently a single Fermi level crossing^{8}.
Let us discuss the individual features of the different terminations shown in Fig. 2 in more detail. The QLterminated surface (Fig. 2b) shows near the Fermi level only the linear crossing of bands, which is protected by mirror symmetry. In the 2QL case, we observe, in addition, an upwarddispersing, Rashbatype spinsplit parabolic surface band with 0.15 eV binding energy at . It touches the bulk valence band along and hybridizes with the protected surface state around −0.1 eV in direction. The BLterminated surface shows characteristic downwards dispersing states, very similar to features observed for a single BL on Bi_{2}Te_{3} (refs 24, 25), that are not seen in the experimental spectra shown below (calculations are presented in the Supplementary Notes II and III). All terminations show strongly spinpolarized surface states around 1.0 eV (Fig. 3a), which are similar to the Rashbatype surface states that also characterize the surface of Bi_{2}Te_{3} (ref. 26) or Sb_{2}Te_{3} (ref. 27).
Crystallographic structure
Figure 4 shows the experimental characterization of the bulk crystal structure of our Bi_{1}Te_{1} thin films via XRD (Fig. 4a–c) and STEM (Fig. 4d). From the ω/2θ scans in Fig. 4a, the crystal phase was determined by comparing the peak positions with the calculated Bragg reflections, both for Bi_{1}Te_{1} and Bi_{2}Te_{3}. The inplane and outofplane lattice constants a and c were determined precisely from the reciprocal space maps around the (1,0,−1,20) reflection for Bi_{2}Te_{3} (Fig. 4b) and the (1,0, −1,16) reflection for Bi_{1}Te_{1} (Fig. 4c). The peaks were fitted with Gaussians for determination of the experimental error. We find a=4.37±0.12 Å and c=30.51±0.46 Å for Bi_{2}Te_{3} and a=4.45±0.02 Å and c=24.0±0.1 Å for Bi_{1}Te_{1}.
In addition, the stoichiometries of the samples were also checked by Rutherford backscattering spectroscopy which confirmed the 50:50 ratio of Bi:Te (see Supplementary Fig. 1).
Figure 4d depicts a highangular annular dark field image of a representative section of a 39nmthick Bi_{1}Te_{1} film recorded by STEM. The observed clear contrast is related to the difference between individual atomic columns of Bi and Te (Bi atomic columns appear brighter than Te columns). Distinct van der Waals gaps, separating QLs from BLs, are visible and the arrangement of BL and QL matches the expected 1:2 composition ratio. Furthermore, by extracting a line profile (yellow frame, which also defines the scale of the image) and fitting Gaussians to the peaks (green=Te; red=Bi) the atomic positions can be determined precisely. Using this method, the size of one bulk unit cell was confirmed to be c=23.99±0.02 Å, which is in good agreement with the results from XRD.
As we have seen, due to the superlattice character of Bi_{1}Te_{1}, there is more than one possible surface termination, but neither XRD nor STEM probe the surface. Supplementary Notes II present a spectroscopic study of the chemical composition of the surface of Bi_{1}Te_{1} and the influence of noble gas sputtering on the surface termination. It turns out that our growth conditions result in Bipoor surfaces (that is, a higher amount of QLterminated surface regions), while ion sputtering leads to Birich surfaces (that is, a higher amount of BLterminated surface areas). This is the reason why, in the following, we will mostly focus on 1QL and 2QLterminated Bi_{1}Te_{1}, since the samples that have been investigated by highresolution ARPES (HR ARPES) were vacuumtransferred after growth and have never been exposed to ambient conditions.
However, since we expect differently terminated surface terraces to be in the order of few micrometer in size^{16,28}, and we employ beam spot sizes of 400 μm (HR ARPES) or even 1 mm (spinARPES) in the ARPES experiments, we have to assume that our electronic structure investigations always probe a superposition of different terminations. Therefore, due to the rich variety of surfacerelated states (Fig. 2) and the fact that our ARPES technique averages over the beam spot size, a detailed distinction of the surface electronic features is challenging.
Surface electronic structure by ARPES
The comparative results of our ARPES investigations on vacuumtransferred, asgrown Bi_{2}Te_{3} and Bi_{1}Te_{1} thin films are summarized in Fig. 3. In the case of the prototypical STI Bi_{2}Te_{3} our results reproduce earlier findings^{26}. In general, the spectra exhibit sharp features and a good signaltonoise ratio revealing the high crystalline quality of the thin films. Figure 3a,b depict wide range binding energy E_{B} versus wavevector k_{,x} maps of Bi_{2}Te_{3} (a) and Bi_{1}Te_{1} (b), respectively, along trajectories in the direction which traverse the point of the surface BZ, recorded with hν=21.2 eV. Both samples are of ntype nature with the conduction band minimum being cut by the Fermi level.
On the first glance the spectra of the two samples show a lot of similarities but a closer analysis reveals some differences, as it can be best seen in energy distribution curves (Fig. 3c,d), which were plotted for k_{,x}=0 Å^{−1} (Fig. 3c) (normal emission) and for k_{,x}=0.39 Å^{−1} (Fig. 3d). The Bi_{1}Te_{1} spectra globally seem to be downshifted which we attribute to a possible electron donation of the BLs to the QLs^{16}.
The spinpolarized surface electronic structure slabcalculations are superimposed onto the ARPES maps (for 1QL surface termination in the case of Bi_{1}Te_{1}). The Fermi level in the calculation needed to be shifted upwards by 250 meV (100 meV) to fit better to the experimental data of Bi_{2}Te_{3} (Bi_{1}Te_{1}). The prominent and intense Rashbatype surface state located between E_{B}=0.7–1.05 eV (0.95–1.3 eV) has been used as a gauge to match the calculation to the ARPES data. As one can see, the agreement between the data and DFT simulation is reasonable and most of the features can be matched.
Figure 3e–g, h–k show the comparison of the magnified electronic structure close to the Fermi level for Bi_{2}Te_{3} and Bi_{1}Te_{1}, respectively. In Bi_{1}Te_{1} (Fig. 3h) the predicted gap opening along direction for 1QL termination is not reproduced in the experiment. However, the superposition (Fig. 3k) of the calculated spectra of 1QL (Fig. 3i) and 2QLterminated Bi_{1}Te_{1} (Fig. 3j) agrees well with the experimental spectrum where 2QL closes the gap. The gap opening in the surface state along is expected since there is no mirror symmetry protecting the bands along this direction. Note, that the Fermi level in the 2QLterminated calculation was shifted by roughly 100 meV more than in the 1QL case pointing to different charge transfer. Supplementary Note II depicts wider energy range simulations where one can see that the Fermi level of the 2QLterminated case indeed needs to be shifted further to match the deeper lying intense Rashbatype features. Additionally, in Supplementary Note IV we analyse the photon energy dependence of the states closest to E_{F} in both Bi_{2}Te_{3} and Bi_{1}Te_{1}. We find in both cases a negligible outofplane dispersion that reveals the surface state character of these features.
In summary, we believe that we do not find a gapopening in the ARPES data of the surface band close to E_{F}, due to the lack of lateral resolution of our measurement technique, probing different terminations simultaneously.
A strong experimental evidence of the topological nature of a state is the verification of the helical spin polarization of a single surface state^{2}. This becomes apparent for the bands marked by red arrows close to E_{F} (Fig. 3f,j), where two bands of opposing chirality nearly overlap in the case of Bi_{1}Te_{1}. Figure 5 summarizes our findings on the spin polarization of the interesting states close to E_{F} of sputtered, that is, Birich, Bi_{1}Te_{1} (Fig. 5a–c). Again the data from Bi_{1}Te_{1} is compared to measurements on Bi_{2}Te_{3} and the spin polarization of the prototypical TSS (Fig. 5d–f).
Figure 5a,d show the wide range ARPES maps of Bi_{1}Te_{1} and Bi_{2}Te_{3} from Fig. 3a,b, respectively, which illustrate along which opposing kpoints, marked by the red dashed area, the spin polarization is measured. Figure 5b,c as well as Fig. 5e,f depict the wide range and nearFermi level (inplane) spinresolved partial intensities I_{left} and I_{right} along the indicated kpoints. The spectra were corrected by the asymmetry function, sometimes called the Sherman function, of S=0.27, and the net spin polarization is shown underneath. Both samples show quite similar and rather high inplane spin polarization of 40–50% in the bands at higher binding energies, around E_{B}≈3.2, ≈2.1 and ≈0.9–1.1 eV in Fig. 5b,e. The full reversal of the spin polarization between the two opposing kpoints confirms the helical nature of these states in both samples. Further, the TSS of Bi_{2}Te_{3} shows a helical spin polarization of up to 40% in Fig. 5f, which confirms its topological nature and is in agreement with what was reported earlier^{26}. On the contrary, Fig. 5c reveals that the most interesting states in Bi_{1}Te_{1} at the Fermi level exhibit only very small (although nonvanishing) inplane spin polarization of max. 10% without a clear reversal at the opposing kpoints. Such weak spin polarization is expected due to spin–orbit coupling (Rashba effect) in the topologically trivial surface states (see states in the calculation in Fig. 3b,c). Therefore, this measurement reveals clear difference to the prototypical TSS and thus gives a strong experimental indication but no final proof for the nontopological character of these states at E_{F} in Bi_{1}Te_{1}.
Finally, we describe our experimental evidence for the mirrorsymmetryprotected band crossings at nonTRIM points, which are a consequence of the topological crystalline character of Bi_{1}Te_{1}. Figure 6 depicts experimental and calculated spectra along nonhighsymmetry lines which reveal a region in kspace, where the TCIinduced states can be identified and no other states interfere. Figure 6a shows a constant energy contour, that is, k_{,x} versus k_{,y} map, at E_{B}=E_{F} with red dashed lines marking ARPES spectra taken along different cut directions at k_{,y}=0.087 Å^{−1} (Fig. 6b), 0.120 Å^{−1} (Fig. 6c), 0.152 Å^{−1} (Fig. 6d), 0.184 Å^{−1} (Fig. 6e), 0.216 Å^{−1} (Fig. 6f) and 0.258 Å^{−1} (Fig. 6g). Figures 6h–m depict the corresponding spinpolarized surface electronic structure calculations, which are performed for a 1QLterminated sample. Besides some additional features in the experimental spectra at 0.4–0.6 eV, the agreement is good and one can identify most bands close to the Fermi level. Importantly, we find the mirrorsymmetryprotected crossing point around k_{,y}=0.184 Å^{−1} and E≈150 meV, which exhibits a gap opening after minor symmetry breaking (see Fig. 2c). Additionally, Fig. 6n depicts the energy distribution curve along the red dashed line in spectrum (Fig. 6g), which shows the gap opening away from the Dirac point in agreement with the DFT results. We attribute these crossing bands at a nonTRIM point to be the consequence of the TCI character of Bi_{1}Te_{1}.
For the k_{} range presented in Fig. 6, the calculation for 1QLterminated Bi_{1}Te_{1} predicts a region without any states at higher binding energies between 400 and 600 meV. In the spectra at Xe 8.4 eV excitation, the experimental spectral weight in that binding energy range is indeed decreased; however, one can still identify dispersive features crossing this region. This is related to a small portion of the surface exhibiting termination other than 1QL or 2QL, which we describe in Supplementary Note V and the related discussion there. There, we show that for the probed k_{} range the spectral weight of the states between 400 and 600 meV strongly depends on the photon energydependent photoemission matrix element and nearly vanishes at He I 21.2 eV excitation. In the related theoretical calculations, this region without bands is still present for the 2QL termination, but not for the Bi BL termination.
Discussion
In summary, we predicted by DFT and demonstrated by ARPES the dual TI character of the stoichiometric natural superlattice Bi_{1}Te_{1}.
Our study theoretically predicts by ab initio DFT calculations that Bi_{1}Te_{1} exhibits a dark surface perpendicular to the stacking direction which is free of timereversal symmetryprotected surface states at the TRIM points, due to weak topological indices (0;001). Moreover, we identify an additional protection of topological states with crossings at nonTRIM points in the mirror plane direction governed by the crystal mirror symmetry with nonzero mirror Chern number . This dual WTI and TCI character leads to the existence of topological states on every surface of the crystal, which are protected either by time reversal or by mirror symmetries, respectively.
Confronting the theoretical predictions with the experiment, we successfully realized thin films of Bi_{1}Te_{1} on Si(111) by MBE growth, characterized the bulk crystal structure as well as the surface chemistry, and thoroughly investigated the (spin) electronic structure. We identified significant differences to the prototypical STI Bi_{2}Te_{3}, that is, the spin polarization of the surfacerelated features at the Fermi level is negligible for Bi_{1}Te_{1}, which points to a nontopological character of this state, and, hence, is an indication of the WTI nature of Bi_{1}Te_{1}. Furthermore, we reveal mirrorsymmetryprotected band crossings at nonTRIM points in excellent agreement to DFT calculations, which we attribute to the TCI character of Bi_{1}Te_{1}.
In future work, the weak topological nature of Bi_{1}Te_{1} could be confirmed by probing the topologically protected onedimensional electron edge channels at step edges of the dark surface, for example, by STS study (similar to what was reported in refs 19, 29, 30).
The dual topological character opens up new vistas because topological states are protected by different symmetries and can be potentially switched on and off separately by breaking one of the symmetries.
Methods
Sample growth
All samples for this study were grown as thin films via MBE. First, 10 × 10 mm^{2} Si(111) samples were prepared by a standard set of wafer cleaning steps (RCAHF) to remove organic contaminations and the native oxide. A consecutive HF dip passivates the Si surfaces with hydrogen for the transfer into the MBE chamber (base pressure (5 × 10^{−10} mbar). To desorb the hydrogen from the surface, the substrates were heated up to 700 °C for 10 min and finally cooled down to 275 °C. For the evaporation of Te and Bi, standard effusion cells were heated to T_{Te}=260 °C and T_{Bi}=460°C, resulting in a growth velocity of Bi_{1}Te_{1} of v=2.5 nm h^{−1}. The tellurium shutter was opened several seconds in advance to terminate the silicon surface by Te, which saturates the dangling bonds. While Bi_{2}Te_{3} is grown in a tellurium overpressure regime^{31}, Bi_{1}Te_{1} requires equal vapour pressures of tellurium and bismuth. The 1:1 ratio between bismuth and tellurium changes the structure from solely QLs in Bi_{2}Te_{3} into a layered structure with additional BLs between every two QLs in Bi_{1}Te_{1}.
After growth, the samples were transferred from the MBE chamber into the ARPES apparatus (<1 × 10^{−10} mbar) without breaking the vacuum, by a ultrahigh vacuum shuttle with a base pressure below 1 × 10^{−9} mbar. The surface of such asgrown samples is, due to the growth mode, expected to be Bipoor, that is, mostly QLterminated. Nevertheless, the surface exhibits all three different terminations (see Supplementary Note II).
Structural characterization
For characterizing the bulk crystal structure, XRD measurements were carried out, employing a highresolution Bruker D8 diffractometer. Additionally, crosssectional specimens were measured in an aberrationcorrected STEM with an electron beam of 0.8 Å (FEI Titan 80–200) for structural investigations on the atomic scale. For this, selected specimens are prepared by focused ion beam etching with firstly 30 keV and subsequently 5 keV Ga ions. Later Ar ion milling using the Fishione NanoMill was performed to reduce the FIBinduced damage. Highresolution STEM images made in highangular annular dark field contain chemical information, since the contrast scales approximate with the square of the atomic number Z^{2}, allowing to distinguish between Bi and Te atoms.
Spectroscopy
The labbased highresolution ARPES investigation was performed at T=25 K with an MBS A1 electron spectrometer, using either a nonmonochromatized He Iα radiation of hν=21.2 eV from a focused helium lamp (Focus HIS 13) or light from a microwavedriven xenon discharge lamp (MBS) producing hν=8.4 eV photons. The beam spot size is about 400 μm in the former and 1 mm in the latter case and the light is unpolarized. The analyser measures E_{B} versus k_{,x} dispersion maps at once. Fermi surface mapping is achieved by rotating the sample with respect to the entrance slit of the spectrometer. The overall energy resolution is estimated to be 10 meV and the angular resolution is <0.02 Å^{−1}.
For spinresolved ARPES measurements we used photons of hν=22 eV, a Scienta SES2002 spectrometer, and a Focus SPLEED polarimeter at beamline BL5 of the DELTA synchrotron in Dortmund with the sample kept at room temperature, resulting in an energy resolution of ≈100 meV^{32}. Here, clean sample surfaces were prepared by sputtering and annealing after sample transfer through air, which resulted in Birich sample surfaces (see Supplementary Note II).
Electronic structure calculations
The DFT calculations are performed for the bulk phase and thin films with three different surface terminations, namely a single and a double QL (1QL and 2QL), and a Bi BL. The bulk unit cell consists of two QLs and one Bi BL, and the hexagonal atomic planes are all assumed to have a facecentred cubiclike (A–B–C) stacking. To simulate a Bi BLterminated surface, a symmetric 26 layer film with BL–QL–QL–BL–QL–QL–BL stacking was used. For the 1QL and 2QL termination, symmetric 24 and 34 layer films were set up. We further checked that slabs of 72 and 84 layers did not change the topological character. We employ the fullpotential linearized augmented plane wave method as implemented in the Fleur code (see http://flapw.de) with the relaxed lattice parameters from the Vienna ab initio simulation package^{33,34}. The generalized gradient approximation of Perdew–Burke–Ernzerhof form^{35} is used for the exchange correlation potential. Spin–orbit coupling is included selfconsistently in the calculations.
From the DFT calculations, we obtain structural parameters that are in good agreement with the experimental data. The size of the bulk unit cell in cdirection is 25.0 Å. It consists of two QLs of 7.48 Å thickness each and a Bi BL of 1.68 Å. The BLQL separation is 2.66 Å and the distance between the QLs is 3.04 Å. At the surfaces, these distances contract slightly, for example, the QL–QL distance decreases by 0.06 Å at the 2QLterminated surface, while the QL–BL distance is reduced only by 0.04 Å for the 1QL termination. For the BL termination, the interlayer distance changes even less. The step height between a BLterminated and a 2QLterminated surface is thus 1.68+2.66=4.34 Å.
It is well known that generalized gradient approximation overestimates the bond length in van der Waalsbonded systems. Therefore, we relaxed the structure also using the DFTD2 method^{36} and obtained a slight contraction of the unit cell in cdirection. The resulting lattice constant of 24.0 Å is in good agreement with the experimental value while the topolocical features are unchanged.
Data availability
The data that support the findings of this study are available from the corresponding author on request.
Additional information
How to cite this article: Eschbach, M. et al. Bi_{1}Te_{1} is a dual topological insulator. Nat. Commun. 8, 14976 doi: 10.1038/ncomms14976 (2017).
Publisher’s note: Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
References
 1
Hasan, M. Z. & Kane, C. L. Colloquium: topological insulators. Rev. Mod. Phys. 82, 3045–3067 (2010).
 2
Ando, Y. Topological insulator materials. J. Phys. Soc. Jpn 82, 102001 (2013).
 3
Hsieh, D. et al. Observation of unconventional quantum spin textures in topological insulators. Science 323, 919–922 (2009).
 4
Fu, L., Kane, C. & Mele, E. Topological insulators in three dimensions. Phys. Rev. Lett. 98, 106803 (2007).
 5
Murakami, S. Phase transition between the quantum spin Hall and insulator phases in 3D: emergence of a topological gapless phase. New J. Phys. 9, 356 (2007).
 6
Teo, J. C. Y., Fu, L. & Kane, C. L. Surface states and topological invariants in threedimensional topological insulators: application to Bi1−xSbx . Phys. Rev. B 78, 045426 (2008).
 7
Fu, L. Topological crystalline insulators. Phys. Rev. Lett. 106, 106802 (2011).
 8
Rauch, T., Flieger, M., Henk, J., Mertig, I. & Ernst, A. Dual topological character of chalcogenides: theory for Bi2Te3 . Phys. Rev. Lett. 112, 016802 (2014).
 9
Murakami, S. Quantum spin Hall effect and enhanced magnetic response by spin–orbit coupling. Phys. Rev. Lett. 97, 236805 (2006).
 10
Wada, M., Murakami, S., Freimuth, F. & Bihlmayer, G. Localized edge states in twodimensional topological insulators: ultrathin Bi films. Phys. Rev. B 83, 121310 (2011).
 11
Yang, F. et al. Spatial and energy distribution of topological edge states in single Bi(111) bilayer. Phys. Rev. Lett. 109, 016801 (2012).
 12
Bos, J., Zandbergen, H., Lee, M.H., Ong, N. & Cava, R. Structures and thermoelectric properties of the infinitely adaptive series (Bi2)m(Bi2Te3)n . Phys. Rev. B 75, 195203 (2007).
 13
Bos, J.W., Faucheux, F., Downie, R. & Marcinkova, A. Phase stability, structures and properties of the (Bi2)m (Bi2Te3)n natural superlattices. J. Solid State Chem. 193, 13–18 (2012).
 14
Cava, R. J., Ji, H., Fuccillo, M. K., Gibson, Q. D. & Hor, Y. S. Crystal structure and chemistry of topological insulators. J. Mater. Chem. C 1, 3176 (2013).
 15
Isaeva, A., Rasche, B. & Ruck, M. Bismuthbased candidates for topological insulators: chemistry beyond Bi2Te3 . Phys. Status Solidi RRL 7, 39–49 (2013).
 16
Valla, T. et al. Topological semimetal in a BiBi2Se3 infinitely adaptive superlattice phase. Phys. Rev. B 86, 241101 (2012).
 17
Rusinov, I. et al. Mirrorsymmetry protected nonTRIM surface state in the weak topological insulator Bi2TeI. Scientific Rep. 6, 20734 (2016).
 18
Fu, L. & Kane, C. Topological insulators with inversion symmetry. Phys. Rev. B 76, 045302 (2007).
 19
Pauly, C. et al. Subnanometrewide electron channels protected by topology. Nat. Phys. 11, 338–343 (2015).
 20
Rasche, B. et al. Stacked topological insulator built from bismuthbased graphene sheet analogues. Nat. Mater. 12, 422–425 (2013).
 21
Pauly, C. et al. Electronic structure of the dark surface of the weak topological insulator Bi14Rh3I9. ACS Nano 10, 3995–4003 (2016).
 22
Mostofi, A. A. et al. wannier90: a tool for obtaining maximallylocalised Wannier functions. Comput. Phys. Commun. 178, 685–699 (2008).
 23
Hsieh, T. H. et al. Topological crystalline insulators in the SnTe material class. Nat. Commun. 3, 982 (2012).
 24
Hirahara, T. et al. Interfacing 2D and 3D topological insulators: Bi(111) bilayer on Bi2Te3 . Phys. Rev. Lett. 107, 166801 (2011).
 25
Miao, L. et al. Quasiparticle dynamics in reshaped helical Dirac cone of topological insulators. Proc. Natl Acad. Sci. USA 110, 2758–2762 (2013).
 26
Herdt, A. et al. Spinpolarization limit in Bi2Te3 Dirac cone studied by angle and spinresolved photoemission experiments and ab initio calculations. Phys. Rev. B 87, 035127 (2013).
 27
Pauly, C. et al. Probing two topological surface bands of Sb2Te3 by spinpolarized photoemission spectroscopy. Phys. Rev. B 86, 235106 (2012).
 28
Gibson, Q. D. et al. Terminationdependent topological surface states of the natural superlattice phase Bi4Se3 . Phys. Rev. B 88, 081108 (2013).
 29
Li, X.B. et al. Experimental observation of topological edge states at the surface step edge of the topological insulator ZrTe5 . Phys. Rev. Lett. 116, 176803 (2016).
 30
Wu, R. et al. Evidence for topological edge states in a large energy gap near the step edges on the surface of ZrTe5 . Phys. Rev. X 6, 021017 (2016).
 31
Borisova, S., Krumrain, J., Luysberg, M., Mussler, G. & Grützmacher, D. Mode of growth of ultrathin topological insulator Bi2Te3 films on Si(111) substrates. Cryst. Growth Des. 12, 6098–6103 (2012).
 32
Plucinski, L., Oelsner, A., Matthes, F. & Schneider, C. M. A hemispherical photoelectron spectrometer with 2dimensional delayline detector and integrated spinpolarization analysis. J. Electron Spectroscop Relat Phenom 181, 215 (2010).
 33
Kresse, G. & Furthmüller, J. Efficient iterative schemes for ab initio totalenergy calculations using a planewave basis set. Comput. Mater. Sci. 6, 15 (1996).
 34
Kresse, G. & Furthmüller, J. Efficient iterative schemes for ab initio totalenergy calculations using a planewave basis set. Phys. Rev. B 54, 11169–11186 (1996).
 35
Perdew, J., Burke, K. & Ernzerhof, M. Generalized gradient approximation made simple. Phys. Rev. Lett. 77, 3865 (1996).
 36
Grimme, S. Semiempirical GGAtype density functional constructed with a longrange dispersion correction. J. Comput. Chem. 27, 1787 (2006).
Acknowledgements
The authors acknowledge financial support from the priority program SPP1666 of the DFG (projects MU3187/31, BI823/2, and MO858/132) and the Virtual Institute for Topological Insulators (VITI) of the HGF. We thank Volkmar Hess, Samuel Königshofen, Frank Matthes and Daniel Bürgler for additional characterization of the surface chemistry of our samples in their laboratorybased XPS chamber. We also thank B. Holländer for the precise determination of the samples stoichiometry via RBS. Further, the authors acknowledge the technical support by B. Küpper and A. Bremen. The authors gratefully acknowledge the computing time granted on the supercomputer JURECA at Jülich Supercomputing Centre (JSC).
Author information
Affiliations
Contributions
M.E., Ma.La. and C.N. contributed equally to this work. Ma.La., P.S. and G.M. grew the samples via MBE. M.E., E.M., P.G., J.K. and L.P. carried out the labbased HR ARPES experiments. C.N., G.B. and S.B. provided the ab initio DFT band structure calculations. M.E., M.G. and S.D. performed spinresolved ARPES measurements at BL5 of DELTA in Dortmund. Ma.La., E.N. and Ma.Lu. prepared the FIB lamellas and performed HR STEM measurements. G.M. characterized the samples by XRD. B.H. measured RBS for stoichiometry analysis. M.E., L.P., C.N. and G.B. wrote the paper with contributions from all coauthors. The project was supervised by M.M., S.B., D.G. and C.M.S.
Corresponding author
Ethics declarations
Competing interests
The authors declare no competing financial interests.
Supplementary information
Supplementary Information
Supplementary Figures, Supplementary Notes and Supplementary References (PDF 3462 kb)
Rights and permissions
This work is licensed under a Creative Commons Attribution 4.0 International License. The images or other third party material in this article are included in the article’s Creative Commons license, unless indicated otherwise in the credit line; if the material is not included under the Creative Commons license, users will need to obtain permission from the license holder to reproduce the material. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/
About this article
Cite this article
Eschbach, M., Lanius, M., Niu, C. et al. Bi_{1}Te_{1} is a dual topological insulator. Nat Commun 8, 14976 (2017). https://doi.org/10.1038/ncomms14976
Received:
Accepted:
Published:
Further reading

Observation and control of the weak topological insulator state in ZrTe5
Nature Communications (2021)

Visualizing coexisting surface states in the weak and crystalline topological insulator Bi2TeI
Nature Materials (2020)
Comments
By submitting a comment you agree to abide by our Terms and Community Guidelines. If you find something abusive or that does not comply with our terms or guidelines please flag it as inappropriate.