Introduction

The gapless spin-polarized states of topological insulators (TIs) at the edge (in two-dimensional (2D) TIs) or at the surface (in three-dimensional (3D) TIs) open up exciting possibilities for applications of these materials in spintronics1,2,3. The classification of TIs is based on topological invariants (v0; v1 v2 v3) and allows to attribute 3D TIs to one of two classes. The v0 = 1 identifies strong TIs, which are characterized by metallic surface states forming an odd number of Dirac cones. These states are robust against perturbations that do not break the time-reversal symmetry. The non-zero v1, v2 and v3 indices define (at v0 = 0) so-called weak TIs that have an even number of Dirac cones. Spin-locked topological surface states (TSSs) in weak TIs should exist on any crystal surface which is not parallel to the plane with Miller indices (v1v2v3)4. The surface states of a weak TI exhibit weaker topological protection than those of a strong TI and can be gapped without breaking the time-reversal symmetry5.

In contrast to the extensive studies on strong 3D TIs, only a limited number of accounts on electronic properties of weak TIs is available at present, e.g. honeycomb compounds XYZ (X = K, Na, Li; Y = Hg, Cd/Au, Ag; Z = Sb, As, P/Te, Se)6, a theoretically modelled octahedron-decorated cubic lattice7, PbTe/SnTe superlattices8, Bi14Rh3I99,10,11 and Bi2TeI12. With an exception of refs 9,11, the above-mentioned papers do not address the (v1v2v3) crystal surface.

Classes of materials with non-trivial band structures are not restricted to TIs. So-called topological crystalline insulators (TCI) represent another type of topological insulating materials with a band gap inverted by the strong spin-orbit coupling (SOC). However, the topological phase therein is protected by a symmetry differing from the time-reversal symmetry, namely by the crystal mirror symmetry13,14. Consequently, these materials possess topologically protected surface states on the surfaces that are perpendicular to the mirror plane.

It was established that the majority of known strong TIs, like Bi2Te3 with the trigonal crystal structure, also accommodate a TCI phase and, thus, the Dirac state residing on the (111) surface is simultaneously protected by both the time-reversal and the crystal mirror symmetries. A magnetic field applied perpendicular to the mirror plane in Bi2Te3 destroys the topological phase and, hence, unravels the dual topological nature of the Dirac state since the TCI phase is preserved and the Dirac surface state remains gapless15.

From this viewpoint, weak TIs are perspective for realization of a dual topological phase where the Dirac surface state of the topological phase and the topological states protected by the crystal symmetry can be manipulated separately without applying a magnetic field, owing to the fact that they appear at different crystal surfaces. In the following, we show that the weak 3D TI Bi2TeI is a promising candidate.

This chemically stable crystalline compound is built by a periodic stack of two types of 2D fragments: BiTeI trilayers and Bi-bilayer16. Previously it was theoretically predicted that Bi2TeI is a weak TI characterized by the (0; 0, 0, 1) invariant12. Accordingly, calculations of the electronic structure of the (010) surface, which is not a natural cleavage surface of the compound, revealed TSS inside the bulk gap12.

Herewith we present an ab-initio density-functional-theory study of the electronic structure of the Bi2TeI natural cleavage (001) surface with different types of possible terminations. By means of these simulations we demonstrate for the first time that the (v1v2v3) surface of a weak 3D TI can accommodate a gapless spin-split surface state in the band gap. The preferable Te termination holds topological surface state, which opposite-spin branches cross at non-symmetric points of the 2D Brilloin zone (BZ) lying in the directions (≈0.3 ) where they are protected by mirror symmetry of the system of the system. Away from the mirror plane a tiny gap of ≈8 meV opens up since the crossing is avoided by a symmetry constraint. On the contrary, the iodine- and [Bi2]-terminated surfaces of Bi2TeI exhibit Rashba-like spin-split bands in the BZ center in addition to an even number of gapless TSSs.

Results

The crystal structure of Bi2TeI was elucidated from a single-crystal X-ray diffraction experiment16. The compound crystallizes in a centrosymmetric monoclinic unit cell (space group C12/m1, a = 7.586 Å, b = 4.380 Å, c = 17.741 Å, β = 98.20°) with 16 atoms which can be represented as a Niggli-reduced cell (a = b = 4.380 Å, c = 17.741 Å, α = β = 82.9047°, γ = 60.0025°) with 8 atoms [Fig. 1(a)]. We performed crystal growth experiments based on the synthetic protocol from16 and determined that the stacking sequence of layers and the lattice parameters are in agreement with the earlier reported results.

Figure 1
figure 1

Crystal structure and bulk electronic spectrum of the Bi2TeI.

(a) A view of the Bi2TeI crystal structure. (b) The 3D and 2D Brillouin zones of a Niggli-reduced cell. (c) Experimental diffraction patterns of the [100]* (left) and the [010]* (centre) zones for a 50 nm thick lamella cut out from a Bi2TeI crystal. Right: a HRTEM image of the [010]* zone with an inserted simulated HRTEM (d = 28 nm, t = 10 nm) outlined by a white dashed line. According to simulations, the bright spots correspond to atoms. The unit cell is outlined by a solid white line. Bulk band spectra of Bi2TeI calculated without (d) and with (e) SOC included. The atomic contributions are color-coded. (f) The same spectra magnified in the vicinity of the Z point. Deep-pink and blue colors denote the states localized in the BiTeI trilayer and the Bi-bilayer, respectively. A schematic view of the inversion of the bulk-band-gap edges for Topological Insulator (g) and Bi2TeI (h). The bands with opposite parities are outlined by dashed and dotted lines; the violet ovals I–III indicate areas of the band inversion.

The layered Bi2TeI structure comprises two 2D building blocks relevant for TI studies: the Bi-bilayer, [Bi2], which is theoretically predicted to be a 2D TI17,18 and BiTeI trilayer, which is a structural unit of the BiTeI compound with giant Rashba-like spin splitting19,20,21,22,23,24. The BiTeI compound can be tuned into a TI phase by external pressure25. Moreover, a single BiTeI trilayer holds the giant Rashba state itself26. These two types of fragments alternate in the stack in such a manner that the Bi-bilayers are always inserted between iodine atomic layers forming [Te-Bi-I] · [Bi2] · [I-Bi-Te] sandwiches. In contrast to covalent bonding within the [Bi2] and the [I-Bi-Te] fragments, the sandwiches are held together by significantly weaker van-der-Waals interactions. As a result, the natural cleavage surface for Bi2TeI is perpendicular to the (001) direction and this ensures absence of trivial surface states within the bulk band gap. The material can be easily cleaved with scotch tape into thin flakes with three possible terminations: tellurium or iodine planes of the [Te-Bi-I] block, or bismuth bilayers. Our experimental results on cleavage show that the Te termination is predominant (more than 80% of instances). Furthermore, we estimated energies for the Te-Te and the Bi2-I cleavages and found that the cleavage between adjacent [Te-Bi-I] · [Bi2] · [I-Bi-Te] sandwiches, is ≈50 times more favorable than that between the [Bi2] and the [I-Bi-Te] fragments, which is consistent with the experimentally found preference for the Te termination.

Stacked layers in the periodic Bi2TeI structure are slightly shifted with respect to each other in the ab plane (this shift amounts to ca. 0.001 Å only within the unit cell). This forces the reduction of the trigonal point symmetry of the [Bi2] and [I-Bi-Te] building blocks, respectively, down to the monoclinic symmetry of the entire Bi2TeI crystal lattice. As a result, the BZ of the (001) surface is a slightly distorted hexagon with the base angles equal to 120.005 and 119.995 degrees, respectively. For this reason, the crystal structure can be regarded as a pseudo-hexagonal one in which only one mirror plane, , is retained, whereas the other two mirror planes of the hexagonal lattice are transformed into pseudo-mirror planes ( and ), see Fig. 1(b). Below the influence of this small distortion from the trigonal symmetry on the electronic structure is discussed.

Spin-orbit interaction (SOC) plays a crucial role in this material. The electronic spectrum of Bi2TeI without SOC included (Fig. 1(d)) has semimetallic character with a zero gap at the Z point. Switching on SOC transforms the electronic spectrum of Bi2TeI from a semimetallic to an insulating one (Fig. 1(e)). Our calculation of the , determined from the product of the parity eigenvalues of the occupied states at the time-reversal invariant momenta (TRIM), confirms the (0; 0, 0, 1) indices obtained by Tang et al.12. However, the SOC-induced band inversion is more complicated in this case than in conventional TIs, like Bi2Te3. In the latter, SOC lifts up the Kramers degeneracy in a valence-band-edge state (formed by Te p-orbitals) and a conduction-band-edge state (formed by Bi p-orbitals) at the TRIM. Since the strength of SOC is larger than the gap width, it leads to an inversion of the edge bands and formation of a new insulating gap owing to their hybridization (see a schematic picture in Fig. 1(g)). As can be seen from Fig. 1(d,e), the SOC-induced band inversion for the Bi2TeI cannot be easily described in terms of atomic orbitals. Nevertheless, if one combines the orbitals that belong to different blocks (Fig. 1(f)) it can be deduced, firstly, that the bulk-band inversion in Bi2TeI occurs between the states of the [Bi2] and the [BiTeI] structural blocks, and, secondly, that two pairs of bands contribute to the SOC-induced band inversion in the resulting spectrum (shown schematically in Fig. 1(h)), in contrast to conventional TIs. The consequence of such complicated band inversion is that the band-gap-inversion area is not centered at the TRIM (see the indicated area I in Fig. 1(h)), as in the case of conventional TI. Furthermore, it triggers formation of additional hybridization gaps between the inverted bands in the valence and conduction bands (areas II and III, respectively). Based on these premises, emergence of a band-gap topological surface state at non-symmetric points as well as the valence- and conduction-band TSSs can be expected.

Simulations of the preferable Te-terminated surface were performed on a slab composed of four [Te-Bi-I] · [Bi2] · [I-Bi-Te] sandwiches, i. e. of 32 atomic layers. Figure 2(a) shows the surface band-structure calculated without inclusion of SOC. This spectrum exhibits semimetallic character, so that the valence and conduction bands touch only at the point and no surface states exist, as can be expected for a surface formed by cleavage through a van-der-Waals gap. The surface spectrum with included SOC (Fig. 2(b)) features a spin-polarized TSS that crosses the band gap along the and directions. Besides, additional spin-split TSSs arise in the valence and conduction bands. They are marked by violet circles II and III, respectively, in Fig. 2(b). Henceforward we focus solely on the band-gap TSS. This state is gapless with a crossing point (CP) lying at ≈0.3 near the conduction band, whereas in the direction it has a tiny gap of 8 meV (Fig. 2(b)). Noteworthy is that the TSS comprises two degenerate states, one for each slab surface. For the considered slab thickness these states are slightly coupled, introducing an artificial minigap (of 2 meV) in the crossing point along the direction. The doubling of the slab in the z direction halves the size of this minigap, whereas the gap in TSS along the remains unchanged. Moreover, identical minigaps are also found in the and directions, thus leading to a conclusion that slight deviations from the trigonal structure in the title compound affect the energy spectrum even weaker than artifacts of the slab model. From these finding we can infer that the crystal-symmetry protection is tolerant (at least for the slab geometry) toward small structural distortions.

Figure 2
figure 2

Electronic structure of the Te-terminated surface.

The band structure of a Te-terminated Bi2TeI (001) slab (black lines) and the projected bulk band structure (gray background) without (a) and with (b) spin-orbit coupling. The dots in the panel (b) represent weights of the states in the outer layers of the slab multiplied by a value of in-plane spin components Sx and Sy (red and blue colors denote the positive and negative projections of the spin vector on Cartesian axes, respectively). The violet ovals indicate: the gapless band-gap (I), the valence-band (II) and the conduction-band (III) TSSs. The inset shows a magnified view of the tiny gap in the band-gap TSS. (c) A 3D view of the band-gap topological surface state. (d) Spatial distribution of the charge density integrated over (x, y) planes, , for both branches of the band-gap TSS. Yellow and green colors accord with those in the panel (c). (e) Spin-resolved constant-energy contours (CECs) taken in the middle of the bulk gap (see the orange dash-doted line in (b)). (fi) CECs for the cuts 1–4 defined in the panel (b).

Energy dependence of the TSS in the full 2D BZ shows that the spin branches with opposite spins form two continuous surfaces with different spatial localization around the center of the BZ (see Fig. 2(c)). One surface (highlighted in yellow) is strongly localized in the outer [Te-Bi-I] trilayer, while the other one (in green) is localized in both [Te-Bi-I] and adjacent [Bi2] blocks (see Fig. 2(d)). There is a tiny gap at the intersection of these two branches at all away of the mirror-planes, where the crossing is avoided by a symmetry constraint27. The TSS remains gapless only in six CPs lying along the directions (marked by black dots in Fig. 2(c)), where it is protected by the crystal mirror symmetry.

The 2D Fermi Surface (FS) is composed of two Γ-centered contours (Fig. 2(e)) for various positions of the chemical potential within the entire bulk band gap. While the inner contour is almost circular, the outer one is subject to strong hexagonal warping. The inner branch of the TSS exhibits clockwise spin-rotation with a negligible out-of-plane spin component, whereas the outer, hexagonally warped branch demonstrates generally counter-clockwise helicity with a more complex spin-texture.

Let us consider the evolution of the 2D FS topology of the TSS at the energies above the conduction band minimum. Right above E = 0, where the TSS has a tiny gap in the directions (see cut 1 in Fig. 2(b,f)), two -centered contours transform into six pockets. At higher energies (cut 2 in Fig. 2(b,g)), the tiny TSS gap shifts away from the point on the high-symmetry direction which initiates further transformations in the FS from six to twelve pockets. Each pocket is centered at the points residing in high-symmetry directions. At E = ECP (cut 3 in Fig. 2(b,h)) the FS again evolves into two distinct -centered contours so that the inner, camomile-like contour touches the outer one at the points lying along the directions. Further shift of the chemical potential towards higher energies (cut 4 in Fig. 2(b,i)) leads to contraction(expansion) of the inner(outer) FS contour.

Let us now regard the less frequent iodine-terminated surface of Bi2TeI. This case was approximated as an [I-Bi-Te] overlayer on top of the Te-terminated slab, so we consider mainly the changes that this add-on introduces in the electronic structure of the Te-terminated surface. The surface spectrum calculated without taking spin-orbit coupling into account (see Fig. 3(a)) is generally similar to the Te-terminated surface spectrum (see Fig. 2(a)) with an exception of a surface state residing at the point at ≈0.2 eV in a local gap of the conduction band (marked by a deep pink curve). This state originates from the splitting of the upper edge of the first conduction band which is caused by positive band bending that is provided by the [I-Bi-Te] overlayer. This case bears similarity to the effects introduced by the iodine-terminated surface in the BiTeI compound21.

Figure 3
figure 3

Electronic structure of the iodine-terminated surface.

An energy spectrum of the iodine-terminated surface with spin-orbit coupling switched off (a) and on (b). (c) Spatial distribution of the charge density for the TSS branches corresponding to the positive (orange) and negative (green) spin projections (note the small colored squares in (b)) and for the trivial surface state at with (black shading) and without (deeppink curve) including SOC. (d) Constant energy contours at the center of the bulk band gap.

The switched-on SOC provokes significant modification of the surface electronic structure as compared to the case of the Te-terminated surface. First of all, SOC induces the bulk-band inversion, so that the bulk conduction band, from which the surface state is split off, dives into the valence band and, consequently, a trivial surface state emerges in the band gap (Fig. 3(b)). In spite of this modification the trivial surface state maintains its localization within the surface trilayers with a maximum in the outer trilayer (Fig. 3(c), deep pink and black curves in the outmost right panel) regardless of whether SOC is taken into account or not. Another consequence of the activated SOC is emergence of a Rashba-type spin splitting at small for the trivial surface state at . At larger the splitting acquires a more complicated character owing to hybridization of the trivial surface state with the topological one. This hybridization also substantially modifies the gapless spin-helical topological surface state which survives at , although the energy of the CP noticeably lowers (with respect to its position at the Te-terminated surface) and approaches the bulk valence band (Fig. 3(b)). As far as localization of the TSS is concerned, the spatial distribution profiles for the opposite spin branches (Fig. 3(c), yellow and green curves) resemble those at the Te-terminated surface (Fig. 2(e)) with the difference that now they penetrate considerably into the [I-Bi-Te] overlayer.

Strong alternation of the electronic spectrum of the Bi2TeI surface induced by the [I-Bi-Te] overlayer causes substantial changes in the 2D Fermi surface. In contrast to two -centered contours provided by the topological surface state (Fig. 2(f)), it is now formed by six isolated egg-shaped pockets enclosing the TSS degeneracies plus two concentric contours from the Rashba-like trivial surface state (RS) in the vicinity of (see Fig. 3(d)).

The last possible cleavage surface of Bi2TeI is terminated by the Bi-bilayer and can be approximated as a [Bi2] · [I-Bi-Te] overlayer on top of the Te-terminated surface or, alternatively, as a [Bi2] overlayer on top of the iodine-terminated surface.

First, the electronic structure of a freestanding [Bi2] · [I-Bi-Te] overlayer is addressed. As can be seen from Fig. 4(a), its spectrum has a gap at the point with a valence-band Rashba-split state at the Fermi level that is mostly localized in the [Bi2]-block in the vicinity of the point, while away from this point the [I-Bi-Te] trilayer contributes predominantly to this state.

Figure 4
figure 4

Electronic structure of a free-standing [Bi2] · [I-Bi-Te] film and the Bi2-terminated surface.

(a) Electronic spectrum of a [Bi2] · [I-Bi-Te] film. (b) Band structure of the Bi2-terminated surface. (c) A magnified view of topological surface state in the direction (note a light blue rectangle in the panel (b)). (d) Charge density distribution of the Rashba-like and topological surface states. The color scheme accords with the colors of the small squares in the panels (b,c). (e) Spin-resolved energy contour for topological surface state below the CP (note a pink line in the panel (c)).

The Rashba-like state formed by the [Bi2] · [I-Bi-Te] overlayer at the Fermi level is retrieved as the most prominent feature in the surface spectrum of the Bi2-terminated Bi2TeI. Comparison of the spectra of a free-standing [Bi2] · [I-Bi-Te] overlayer and the Bi2-terminated surface shows that the dispersion and spatial localization of the discussed state do not change (Fig. 4(d), left) as it remains localized in the [Bi2] block. The major difference is, however, that the spectrum is gapless in the case of the Bi2-terminated surface (see a light blue rectangle in Fig. 4(b)). This state with the CP at ≈0.25 (Fig. 4(c)) below the conduction band can be regarded as the survived topological state which penetrated deeply into the adjacent sublayers down to a second trilayer (Fig. 4(d), right). The spin texture of this state is highly unusual because of the strong hybridization with the trivial Rashba state (Fig. 4(e)), namely, the spin, being mostly in-plane, always has a positive Sx component around the contour. The counterpart gapless state located in the direction has an opposite spin direction and, thus, the net spin equals zero over the Brillouin zone.

Methods

Crystal growth

A stoichiometric mixture of Bi, Te and BiI3 (sublimated in vacuum prior to use) was placed into a silica ampoule that was then evacuated and sealed. The ampoule was heated to 823 K with a rate of 10 K/h, tempered for 2 hours and subsequently cooled down to ambient temperature with a rate 2 K/h. Largely overgrown crystalline platelets were found on the batch along with some single crystals that grew on the ampoule’s walls. The products were characterized by semi-quantitative energy dispersive X-ray analysis (SU8020 (Hitachi) SEM, Silicon Drift Detector (SDD) X − MaxN (Oxford)), X-ray diffractometry (X’Pert Pro MPD diffractometer (PANalytical), Ge(111) monochromator, Cu-Kα radiation) and TEM methods (FEI Titan F20 microscope with CS-correction operating at 80 kV). A typical SAED and HRTEM images for Bi2TeI are given in Fig. 1(c) of the main text. The observed stacking sequence of layers and lattice parameters are in consistent with the those determined from the structure elucidation in ref. 16.

DFT calculations

Electronic structure calculations were carried out within the density functional theory using the projector augmented-wave method28 as implemented in the VASP code29,30 and ABINIT code31. The PAW data sets in ABINIT code were taken from ref. 32. The exchange-correlation energy was treated using the generalized gradient approximation33. The Hamiltonian contained the scalar relativistic corrections and the spin-orbit coupling was taken into account. The calculation of invariant was carried by using the parity of the wave functions obtained in the framework of the Full Potential Linearized Augmented Plane wave (FLAPW) method implemented in FLEUR code34. The bulk and surface spectra obtained using different codes are in full agreement.

Additional Information

How to cite this article: Rusinov, I. P. et al. Mirror-symmetry protected non-TRIM surface state in the weak topological insulator Bi2TeI. Sci. Rep. 6, 20734; doi: 10.1038/srep20734 (2016).