Home | About | Journals | Submit | Contact Us | Français |

**|**Scientific Reports**|**PMC5333630

Formats

Article sections

Authors

Related links

Sci Rep. 2017; 7: 43666.

Published online 2017 March 2. doi: 10.1038/srep43666

PMCID: PMC5333630

I. A. Nechaev,^{a,}^{1,}^{2,}^{3} S. V. Eremeev,^{2,}^{3,}^{4,}^{5} E. E. Krasovskii,^{5,}^{6,}^{7} P. M. Echenique,^{5,}^{6} and E. V. Chulkov^{1,}^{2,}^{3,}^{5,}^{6}

Received 2016 October 27; Accepted 2017 January 26.

Copyright © 2017, The Author(s)

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/

The quantum spin Hall insulators predicted ten years ago and now experimentally observed are instrumental for a break- through in nanoelectronics due to non-dissipative spin-polarized electron transport through their edges. For this transport to persist at normal conditions, the insulators should possess a sufficiently large band gap in a stable topological phase. Here, we theoretically show that quantum spin Hall insulators can be realized in ultra-thin films constructed from a trivial band insulator with strong spin-orbit coupling. The thinnest film with an inverted gap large enough for practical applications is a centrosymmetric sextuple layer built out of two inversely stacked non-centrosymmetric BiTeI trilayers. This nontrivial sextuple layer turns out to be the structure element of an artificially designed strong three-dimensional topological insulator Bi_{2}Te_{2}I_{2}. We reveal general principles of how a topological insulator can be composed from the structure elements of the BiTe*X* family (*X*=I, Br, Cl), which opens new perspectives towards engineering of topological phases.

Two-dimensional (2D) topological insulators (TIs)—a new electronic phase also referred to as a quantum spin Hall (QSH) insulator—are characterized by an absolute band gap induced by spin-orbit coupling (SOC) and helical gapless edge states inside the gap^{1}. These states protected by time-reversal symmetry provide perfectly conducting spin-filtered channels, meeting the demands of low-power nanoelectronics and spintronics. The existence of such states as the fingerprint of a topologically non-trivial 2D insulator was first predicted in Refs ^{2}^{,}^{3}. It was also suggested that the QSH effect can be observed in graphene, where SOC opens a gap at the two inequivalent Dirac points. This gap in graphene appears to be too small for practical use, so heavy-elements based analogs of graphene must be sought. Actually, the 2D materials with low-buckled honeycomb-lattice structures^{4}^{,5}^{,6}^{,7}^{,8}^{,9}^{,10}—silicene, germanene, and stanene—possess a significantly larger SOC-induced gap at the Dirac points (up to ~0.1eV in stanene), and the spin-polarized edge channels could be detected at easily accessible temperatures. However, the QSH effect in such systems has not been experimentally observed so far.

Further effective enhancement of the SOC to make the gap larger can be realized by chemical functionalization of the above 2D materials^{11}^{,12}. Such a functionalization substantially enlarges the gap, in fact “destroying” the Dirac cones, and it may lead to a SOC-induced band inversion at the time reversal invariant momentum (TRIM) (normally at **k**=0) with an absolute gap of several hundred meV at this momentum. If the band inversion occurs, the resulting 2D system is a 2D TI that should support the QSH effect^{12} similar to the inverted HgTe quantum wells predicted to be QSH insulators in ref. ^{13}. It is important that this prediction has found experimental confirmation: In the inverted HgTe/CdTe and InAs/GaSb quantum wells^{14}^{,15}^{,16} the QSH effect was observed despite the very small gaps in these quantum wells, less than 10meV. It has spurred a rising tide of theoretical propositions of different 2D TIs with honeycomb- or square-lattice structures and a large inverted gap enabling room-temperature operating^{17}.

*Ab initio* approaches to electronic structure, especially those based on the density functional theory (DFT), have become a powerful tool to search for new materials with unique properties. At the same time, the effective models that proved indispensable in predicting the QSH effect in graphene-like systems and quantum wells are currently widely used to analyze the effect of strain, quantum confinement, and external fields in 2D TIs, i.e., to solve the problems that presently are not accessible with *ab initio* methods. Thus, to efficiently model the nanoelectronics and spintronics devices, the microscopic methodology must be bridged with the effective Hamiltonian approach based on symmetry considerations and on the **k·p** perturbation theory.

With a few exceptions, none of the theoretically proposed 2D materials has been hitherto fabricated^{17}. Thus, the intensive search of robust and easily fabricated materials remains to be actual. In particular, it was suggested that 2D TIs can be produced from a thin film of layered 3D (three-dimensional) TIs of the Bi_{2}Se_{3} family, where the hybridization between the opposite surfaces of the film opens a gap at the Dirac point (DP). Depending on film thickness, the 2D system may “oscillate” between band insulator and QSH insulator as was predicted by the 4-band effective **k·p** model (see refs ^{18}^{,}^{19}). The thinnest known topologically nontrivial film consists of at least two structural elements—quintuple or septuple layers^{20}.

Besides the studies on the thin films of 3D TIs, recently it was heuristically suggested that a 3D TI can be constructed artificially via stacking 2D bilayers that are topologically trivial^{21}. It encourages our search for 2D TIs built out of trivial band-insulator constituents. These constituents should have strong spin-orbit coupling (SOC), and prospective candidates are bismuth tellurohalides BiTe*X* with *X*=I, Br, and Cl, among which the polar semiconductor BiTeI demonstrates the strongest spin-orbit coupling providing the biggest known Rashba spin-splitting of bulk and surface states^{22}^{,23}. The structure element of BiTeI is a trilayer (TL) with the I-Bi-Te stacking. A single TL that possesses the Rashba spin-split band structure^{24} can be grown epitaxially on a suitable substrate or be easily exfoliated from the bulk BiTeI, where the adjacent TLs couple through a weak van-der-Waals (vdW) interaction. The samples of BiTeI always contain a large number of randomly distributed bulk stacking faults, which leads to a mixture of terminations at the surface, as experimentally observed in refs ^{25}^{, }^{26}^{, }^{27}^{, }^{28}^{, }^{29}^{, }^{30}. This implies that adjacent TLs may have different sequence order along the hexagonal **z** axis.

Here, based on DFT calculations we demonstrate that a centrosymmetric sextuple layer (SL) constructed from two BiTeI TLs with facing Te-layer sides and a typical vdW spacing is a 2D TI with the gap of 70meV at . The vdW interaction between these TLs is crucial to realize such a QSH insulator phase: The SL becomes topologically trivial with increasing the vdW spacing by 5% only. We consider the *nontrivial* SL as a structure element, a repetition of which along the **z** axis results in thin films that are found to “oscillate” between trivial and nontrivial phases with the number of SLs. The corresponding bulk system composed of SLs turns out to be a *strong* 3D TI (hereafter referred to as Bi_{2}Te_{2}I_{2}). It is energetically unfavourable by only 0.5meV compared with the non-centrosymmetric BiTeI. This makes it plausible to suppose that crystals of BiTeI grown by the Bridgman method already contain the desired SLs, and that the alternative stacking can be experimentally observed and controllably manufactured. To describe the low-energy properties of Bi_{2}Te_{2}I_{2} and its films, we derive four-band **k·p** Hamiltonians from the *ab initio* wave functions. They are similar to the Hamiltonians constructed for Bi_{2}Se_{3} -family 3D TIs and their thin films^{18}^{,31}^{,32}. For a more accurate description of the SL, we derive an eight-band Hamiltonian that involves Rashba-split valence and conduction bands of the stand-alone TLs. We thus demonstrate that due to the bonding-antibonding splitting the inversion occurs between one of the Te-related valence bands and one of the conduction bands formed by Bi orbitals. The proposed materials illustrate the effectiveness of the new way to design 2D TIs from trivial band insulators with giant-Rashba spit bands for room-temperature operating.

Figure 1(a–c) show the band structure of 1 and 5 SLs and the bulk crystal of Bi_{2}Te_{2}I_{2} obtained with the extended linearized augmented plane wave (ELAPW) method^{33} within the local density approximation (LDA) for the exchange-correlation functional and with the use of the full potential scheme of ref. ^{34}. (Details on the equilibrium bulk atomic structure, the bulk-truncated slab geometry of the related thin films, and the calculations performed can be found in Supplementary Note 1.)

The 1SL film is constructed from two BiTeI trilayers with facing Te-layer sides, Fig. 1(d). It is noteworthy that the band structure of this film with the gap of 56meV (70meV in the relaxed geometry, see Supplementary Fig. 1) differs substantially from that of its constituents (cf. Supplementary Fig. 2): there is no trace of Rashba-type split bands. The band structure of the 5SL film exhibits a gapless Dirac state residing in the band gap of 151meV, see Fig. 1(b) and Supplementary Fig. 3. This is a signature of the topological character of the respective bulk band structure (Fig. 1(c)), which has an inverted gap of 234meV at Γ and the fundamental gap of 169meV in the Γ–A line close to the A point. As seen in Fig. 2(a), the Dirac surface state almost completely resides within the outer SL. Moreover, this state is localized stronger than the Dirac state of TIs like Bi_{2}Te_{3}, since 70% of its weight falls in the outermost half of the SL, i.e., in the surface BiTeI TL (see also Fig. 1(e)).

The spin texture of the Dirac state is illustrated in Fig. 2(b) and (c), which show spin-resolved constant energy contours for the lower and upper cones of the Dirac surface state. Apart from the in-plane polarization – clockwise above the DP and counterclockwise below it, both contours also have an out-of-plane spin component, which is an intrinsic feature of the hexagonal surface. However, in this case *S*_{z} is extremely small and varies in the range of ≈±0.01.

As has been shown for the Dirac state in Bi_{2}Se_{3} both experimentally and theoretically^{35}^{,36}^{,37}^{,38}^{,39}, the spin textures of *p*_{x}, *p*_{y}, and *p*_{z} orbitals are remarkably different, which leads to the dependence of the spin polarization of photoelectrons on the polarization of light. The spin texture provided by *p*_{z} orbitals has clockwise (counterclockwise) chirality for the upper (lower) cone, while the projections of the total spin on *p*_{x} and *p*_{y} orbitals are not chiral, and their spins are opposite to each other. Similar spin-orbital texture we find in the Bi_{2}Te_{2}I_{2}, see Fig. 2(d–f) for the upper Dirac cone (for the lower Dirac cone, the coupling of spin and orbital textures is opposite, not shown). As can be seen, the spin orientations for *p*_{x} and *p*_{y} projections are antiparallel at each *k*_{||} point, whereas the spin orientation of *p*_{z} projection coincides with the total spin.

It is noteworthy that similar spin-orbital texture has been observed for the spin-polarized Rashba state in BiTeI^{40}. For the Te-terminated surface, it was found that the outer Rashba branch demonstrates the same spin orientations for *p*_{x}, *p*_{y}, and *p*_{z} projections as those for the upper Dirac cone in Bi_{2}Te_{2}I_{2}, and the inner Rashba branch has opposite texture, i.e., the same as in the lower cone. Because the surface of the Bi_{2}Te_{2}I_{2} slab has iodine termination and its spin-texture is reversed due to the opposite orientation of the **z** axis, the spin-orbit texture of the upper (lower) Dirac cone in Bi_{2}Te_{2}I_{2} is the same as the texture of the inner (outer) Rashba branch in BiTeI.

To construct a simple effective **k·p** model for the centrosymmetric Bi_{2}Te_{2}I_{2} we derive a model Hamiltonian of a desired dimension and accurate up to the second order in **k** from the LDA spinor wave functions **Ψ**_{n↑(↓)} of the doubly degenerate bands *E*_{n} found at **k**=0 (see Supplementary Note 2 and ref. ^{41} for details). The subscripts ↑ or ↓ in **Ψ**_{n↑(↓)} refer to the *z*-component of the total angular momentum **J**=**L**+**S** in the atomic sphere that has the largest weight in the *n*-th band, see Fig. 1(e). The Hamiltonian is constructed in terms of the matrix elements^{42}
of the velocity operator , where *n* and *m* run over the relativistic bands (from semi-core levels up to high-lying unoccupied bands). Here, * σ* is the vector of the Pauli matrices, and

For the bulk Bi_{2}Te_{2}I_{2}, in the basis of the two valence bands and two conduction bands , our *ab initio* four-band Hamiltonian reads:

where , , , and the direct matrix product of the Pauli matrices * τ* and

The parameters in Eq. (1) obtained within the LDA are: C_{0}=0.03 eV, *C*_{z}=0.13 a.u., *C*_{||}=4.19 a.u., *M*_{0}=−0.12eV, *M*_{z}=1.35 a.u., *M*_{||}=5.88 a.u., *V*_{||}=0.52 a.u., and *V*_{z}=0.13 a.u. (we use Rydberg atomic units: ). Since the basis functions explicitly refer to the valence and conduction bands rather than to atomic orbitals, the parameter *M*_{0} that defines the band gap at **k**=0 is negative and does not change sign upon moving from the topologically non-trivial insulator to the trivial one. The eigenvalues *E*(**k**) of the Hamiltonian (1) with the above parameters are shown in Fig. 1(c) by red lines, nicely reproducing the LDA curves over a quite large **k**-region and providing an absolute gap in the **k·p** spectrum. Moreover, these parameters reflect the band inversion and meet the condition of the existence of topological surface states (see, e.g., ref. ^{43}) in accord with the _{2} topological invariant ν_{3D}=1 obtained from the parities of the bulk LDA wave functions at the TRIM points^{44}. Actually, the diagonal dispersion term *M*_{z(||)} is positive, and it is larger than the electron-hole asymmetry: |*C*_{z(||)}|<*M*_{z(||)}.

For the Bi_{2}Te_{2}I_{2} thin films, we derive the Hamiltonian in the basis as

where , , and * τ* refers now to the two decoupled sets of massive Dirac fermions. The Hamiltonian (2) is similar to the one obtained for 3D TI thin films within the effective continuous model based on the substitution in the Hamiltonian of ref.

All the considered Bi_{2}Te_{2}I_{2} films are characterized by the velocity *V*_{||}=0.45±0.01 a.u. and the electron-hole asymmetry *C*_{||}=4.15±0.10 a.u., which are weakly sensitive to the number of SLs, where the±ranges indicate the variations of *V*_{||} and *C*_{||} in moving from 1 to 5 SLs. On the contrary, as seen in Fig. 1(g) and (h) the parameters *M*_{0} and *M*_{||} depend strongly on the film thickness, approaching monotonically zero.

In order to explicitly indicate whether a given film is a QSH insulator, in Fig. 1(g) we also plot the gap parameter with *ν*_{2D} being the _{2} invariant obtained from the parities of the wave functions at the TRIM points of the 2D Brillouin zone. This parameter is negative for a topologically non-trivial film and positive for a trivial one. As follows from the figure, **Δ** “oscillates” with the period of 2 SLs within the examined thickness interval. (The parity of and is (+) and (−), respectively, for **Δ**<0, and it is (−) and (+) for **Δ**>0). As in 3D TI films^{45}, the thickness dependence of **Δ** may be sensitive to the quasi-particle approximation employed, and it may change if many-body corrections beyond DFT are introduced. However, even the simplest quasi-particle method, the *GW* approximation for the self-energy, is methodologically challenging and computationally too demanding to study a large series of complex systems. Thus, DFT remains the method of choice, and its good performance for a wide range of TIs justifies the use of the Kohn-Sham band structure as a reasonable starting point.

The diagonalization of the Hamiltonian (2) then leads to *E*(**k**) shown by red lines in Fig. 1(a) and (b). The absence of the absolute gap in the resulting **k·p** spectrum is the general feature of all the films studied. It is caused by the rather big electron-hole asymmetry *C*_{||} compared with the diagonal dispersion parameter *M*_{||}, Fig. 1(h). It should be noted that the conclusion on whether the edge states exist in a TI film is often made based on the signs and relative values of the parameters *M*_{0}, *M*_{||}, and *C*_{||}. On the contrary, we find that the asymmetry |*C*_{||}| is larger than |*M*_{||}| everywhere, breaking one of the conditions for the film to be a QSH insulator, see, e.g., refs ^{18}^{,}^{46}. Focusing on the behaviour of the diagonal dispersion *M*_{||} (as, e.g., in the topology analysis of ref. ^{15}), we note that it is positive for all the thicknesses, Fig. 1(h). Along with the negative *M*_{0}, this should signify an inverted band gap for the respective films. However, it does not correlate with the oscillating **Δ**, Fig. 1(g).

Let us now analyze the behaviour of the diagonal dispersion parameter *M*_{||} together with the topological invariant under a continuously varying geometry. We choose the 1SL film–the thinnest film, for which the **k·p** prediction of the band inversion does not contradict the actual topological property–and gradually expand the van-der-Waals spacing *d*_{vdW}. The evolution of the band structure with increasing *d*_{vdW} is shown in Fig. 3. According to the gap parameter **Δ**, see Fig. 4(a), a topological phase transition occurs at *d*_{vdW} that is just around the mentioned 5% larger than its bulk value, and the 1SL film becomes topologically trivial. Further expansion leads to a larger band gap at , which is not inverted anymore. It is noteworthy that such a behaviour of **Δ** as a function of *d*_{vdW} with the topological phase transition around 5% is stable with respect both to the choice of the approximation to the DFT exchange-correlation functional (LDA, GGA, dispersion corrected GGA) and to the SL geometry (bulk truncated or relaxed). In the limit of very large *d*_{vdW}, when the BiTeI trilayers composing the 1SL film are too far from each other, the band structure is identical to that of a free-standing BiTeI trilayer (see Supplementary Fig. 2). Similarly, artificial reduction of the spin-orbit interaction strength *λ* relative to its actual value *λ*_{0} in the equilibrium SL leads to a decrease in the gap, which closes at *λ*/*λ*_{0}=0.95. A further decrease in *λ* causes a widening of the already uninverted gap of the trivial phase. In general, the dependence of the relative gap-width on the spin-orbit interaction strength is almost linear and can be approximated as .

The 1SL parameters of the 4-band **k·p** Hamiltonian (2) strongly depend on *d*_{vdW} (the respective eigenvalues *E*(**k**) of this Hamiltonian are shown by red lines in Fig. 3(a–c)). With the *d*_{vdW} expansion (given in percents of the bulk value ) up to 50%, the velocity *V*_{||} decreases monotonically from 0.470 a.u. to 0.342 a.u., and the electron-hole asymmetry *C*_{||} becomes smaller as well, Fig. 4(c). At , *C*_{||} is already smaller than *M*_{||}, ensuring an absolute gap in the 4-band **k·p** spectrum, see Fig. 3(c). With further increasing it even becomes negative, but it remains |*C*_{||}|<*M*_{||}. A stepwise behaviour of the parameter *M*_{||} that changes sign at the small indicates that *M*_{||} keeps following the actual *ν*_{2D} and, thus, predicts a gap without inversion. With increasing *d*_{vdW} this parameter again goes through zero around , telling us that the band gap becomes inverted again, and at ~35% with the given *C*_{||} and *M*_{0} meets the conditions of the existence of the edge states^{18}^{,46}. However, as seen in Fig. 4(a), the 1SL film is too far from a topological phase transition at such . With this example we illustrate the strong limitations of the predictive capabilities of the effective continuous model.

Let us now analyze the formation of the SL band structure with the inverted band gap. Starting from well-separated layers, Fig. 3(f), and going back to the bulk value of the van-der-Waals spacing, Fig. 3(a), we retrace the valence bands (*ν*_{1} with the energy and *ν*_{2} with ) with the predominant contribution coming from the *p*_{z} orbitals of Te and the conduction bands (*c*_{1} with and *c*_{2} with ) mainly formed by Bi *p*_{z} orbitals, see Supplementary Fig. 2. We derive an 8-band Hamiltonian which is presented in Supplementary Note 3. Its eigenvalues are shown in Fig. 3 by blue lines, and the corresponding parameters as a function of the *d*_{vdW} expansion are depicted in Fig. 4(b) and (d), see also Supplementary Fig. 4.

As seen in Fig. 4(b), at in the large-*d*_{vdW} limit there are two doubly degenerate energy levels, and . Upon decreasing *d*_{vdW}, the TLs start to interact primarily by their Te-layer sides to cause the bonding-antibonding splitting of the two degenerate levels: The Te-related energies as a function of *d*_{vdW} disperse stronger than those of Bi. Near the bulk value , the splitting is large enough to invert the order of the and levels, ensuring the topological phase transition. Thus, the stacking procedure that leads to the 3D TI is based on SL building blocks principally different from the Rashba bilayers used in ref. ^{21}. It is essential that in our case the two Rashba constituents of the block (the stand-alone TLs) bring not only the Rashba-split conduction band but also the valence band, see Supplementary Note 3. Then the gap in the SL (which may be inverted or not) is quite naturally the gap between the valence and conduction bands, in contrast to the scenario of ref. ^{21}, where the band gap in the bilayer block is achieved by a dispersive “finite quantum tunneling” between the two Rashba constituents – the 2D electron gases of the adjacent layers.

Figure 4(d) shows the behaviour of the inverse effective masses of the chosen bands over the *d*_{vdW} interval considered. We find that the conduction-band inverse masses, which are equal in the large-*d*_{vdW} limit, , change smoothly with decreasing *d*_{vdW}: At the parameter becomes twice as large, while falls below zero. On the contrary, the valence-band inverse masses ( in the large-*d*_{vdW} limit) “diverge” because the band *ν*_{1} moves down and “goes through” the I-orbital dominated bands, and *ν*_{2} moves up and hybridizes with Te *p*_{x,y} bands, see Fig. 3 and Supplementary Fig. 2. Finally, at the parameter reaches its large *d*_{vdW} limit, while becomes negative. Thus, in the topologically non-trivial 1 SL we have and , where *C*_{||} and *M*_{||} are the 1SL parameters of the Hamiltonian (2). At that, the interband coupling of the bands *ν*_{2} and *c*_{1} is equal to *V*_{||} of the 4-band **k·p** description. This reveals a close relation between the 4-band and 8-band Hamiltonians. However, already with 8 bands there is an absolute gap (see Fig. 3(a)), which is reasonably accurate and quite suitable for the theoretical research on linear response, Hall conductance, and motion of Dirac fermions in external fields.

**How to cite this article:** Nechaev, I. A. *et al*. Quantum spin Hall insulators in centrosymmetric thin films composed from topologically trivial BiTeI trilayers. *Sci. Rep.*
**7**, 43666; doi: 10.1038/srep43666 (2017).

**Publisher's note:** Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

This work was supported by the Spanish Ministry of Economy and Competitiveness MINECO (Project Nos FIS2013-48286-C2-1-P, FIS2013-48286-C2-2-P, and FIS2016-76617-P), the Basque Country Government, Departamento de Educación, Universidades e Investigación (Grant No. IT-756-13) and Saint Petersburg State University (Grant No. 15.61.202.2015).

The authors declare no competing financial interests.

**Author Contributions** I.A.N. and S.V.E. conceived the idea and designed the research within the projects coordinated by P.M.E. and E.V.C., S.V.E. performed atomic structure optimization and GGA band structure calculations. I.A.N. and E.E.K. derived the small-size **k·p** Hamiltonians from *ab initio* wave functions and performed the corresponding LDA-based calculations. I.A.N., S.V.E., and E.E.K. wrote the manuscript. All the authors discussed the results and commented on the manuscript.

- Qi X.-L. & Zhang S.-C. Topological insulators and superconductors. Rev. Mod. Phys. 83, 1057 (2011).
- Kane C. L. & Mele E. J.
_{2}Topological Order and the Quantum Spin Hall Effect. Phys. Rev. Lett. 95, 146802 (2005). [PubMed] - Kane C. L. & Mele E. J. Quantum Spin Hall Effect in Graphene. Phys. Rev. Lett. 95, 226801 (2005). [PubMed]
- Cahangirov S., Topsakal M., Aktürk E., Şahin H. & Ciraci S. Two- and One-Dimensional Honeycomb Structures of Silicon and Germanium. Phys. Rev. Lett. 102, 236804 (2009). [PubMed]
- Lalmi B. et al. . Epitaxial growth of a silicene sheet. Appl. Phys. Lett. 97, 223109 (2010).
- Liu C.-C., Feng W. & Yao Y. Quantum Spin Hall Effect in Silicene and Two-Dimensional Germanium. Phys. Rev. Lett. 107, 076802 (2011). [PubMed]
- Liu C.-C., Jiang Hua & Yao Y. Low-energy effective Hamiltonian involving spin-orbit coupling in silicene and two-dimensional germanium and tin. Phys. Rev. B 84, 195430 (2011).
- Ezawa M. A topological insulator and helical zero mode in silicene under an inhomogeneous electric field. New Journal of Physics 14 033003 (2012).
- Dávila M. E., Xian L., Cahangirov S., Rubio A. & Lay G. Le. Germanene: a novel two-dimensional germanium allotrope akin to graphene and silicene. New Journal of Physics 16 095002 (2014).
- Zhu Feng-feng et al. . Epitaxial growth of two-dimensional stanene. Nature Materials 14, 1020 (2015). [PubMed]
- Ma Y., Dai Y., Niu C. & Huang B. Halogenated two-dimensional germanium: candidate materials for being of Quantum Spin Hall state. J. Mater. Chem. 22, 12587 (2012).
- Xu Y. et al. . Large-Gap Quantum Spin Hall Insulators in Tin Films. Phys. Rev. Lett. 111, 136804 (2013). [PubMed]
- Bernevig B. A., Hughes T. L. & Zhang S.-C. Quantum Spin Hall Effect and Topological Phase Transition in HgTe Quantum Wells. Science 314, 1757 (2006). [PubMed]
- König M. et al. . Quantum Spin Hall Insulator State in HgTe Quantum Wells. Science 318, 766 (2007). [PubMed]
- König M. et al. . The Quantum Spin Hall Effect: Theory and Experiment. Journal of the Physical Society of Japan 77, 031007 (2008).
- Knez I., Du R.-R. & Sullivan G. Evidence for Helical Edge Modes in Inverted InAs/GaSb Quantum Wells. Phys. Rev. Lett. 107, 136603 (2011). [PubMed]
- Ren Y., Qiao Z. & Niu Q. Topological phases in two-dimensional materials: a review. Rep. Prog. Phys. 79, 066501 (2016). [PubMed]
- Lu H.-Z., Shan W.-Y., Yao W., Niu Q. & Shen S.-Q. Massive Dirac fermions and spin physics in an ultrathin film of topological insulator. Phys. Rev. B 81, 115407 (2010).
- Liu C. X. et al. . Oscillatory crossover from two-dimensional to three-dimensional topological insulators. Phys. Rev. B 81, 041307 (2010).
- Bihlmayer G., Koroteev Yu. M., Menshchikova T. V., Chulkov E. V. & Blügel S. Ab Initio Calculations of Two-Dimensional Topological Insulators In Topological Insulators, Fundamentals and Perspectives (Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim, 2015).
- Das T. & Balatsky A. V. Engineering three-dimensional topological insulators in Rashba-type spin-orbit coupled heterostructures. Nature Commun. 4, 1972 (2013). [PMC free article] [PubMed]
- Ishizaka K. et al. . Giant Rashba-type spin splitting in bulk BiTeI. Nature Mater. 10, 521 (2011). [PubMed]
- Eremeev S. V., Nechaev I. A., Koroteev Yu. M., Echenique P. M. & Chulkov E. V. Ideal two-Dimensional electron systems with a giant Rashba-type spin splitting in real materials: surfaces of bismuth tellurohalides. Phys. Rev. Lett. 108, 246802 (2012). [PubMed]
- Eremeev S. V., Tsirkin S. S., Nechaev I. A., Echenique P. M. & Chulkov E. V. New generation of two-dimensional spintronic systems realized by coupling of Rashba and Dirac fermions. Sci. Rep. 5, 12819 (2015). [PMC free article] [PubMed]
- Landolt G. et al. . Disentanglement of Surface and Bulk Rashba Spin Splittings in Noncentrosymmetric BiTeI. Phys. Rev. Lett. 109, 116403 (2012). [PubMed]
- Crepaldi A. et al. . Giant Ambipolar Rashba Effect in the Semiconductor BiTeI. Phys. Rev. Lett. 109, 096803 (2012). [PubMed]
- Tournier-Colletta C.
et al. . Atomic and electronic structure of a Rashba
*p-n*junction at the BiTeI surface. Phys. Rev. B 89, 085402 (2014). - Fiedler S. et al. . Defect and structural imperfection effects on the electronic properties of BiTeI surfaces. New J. Phys. 16, 075013 (2014).
- Butler C. J. et al. . Mapping polarization induced surface band bending on the Rashba semiconductor BiTeI. Nature Commun. 5, 4066 (2014). [PMC free article] [PubMed]
- Fiedler S.
et al. . Termination-dependent surface properties in the giant-Rashba semiconductors BiTe
*X (X*=Cl, Br, I). Phys. Rev. B 92, 235430 (2015). - Zhang H.
et al. . Topological insulators in Bi
_{2}Se_{3}, Bi_{2}Te_{3}, and Sb_{2}Te_{3}with a single Dirac cone on the surface. Nat. Phys. 5, 438 (2009). - Zhang F., Kane C. L. & Mele E. J. Surface states of topological insulators. Phys. Rev. B 86, 081303(R) (2012).
- Krasovskii E. E. Accuracy and convergence properties of the extended linear augmented-plane-wave method. Phys. Rev. B 56, 12866 (1997).
- Krasovskii E. E., Starrost F. & Schattke W. Augmented Fourier components method for constructing the crystal potential in self-consistent band-structure calculations. Phys. Rev. B 59, 10504 (1999).
- Zhang H., Liu C.-X. & Zhang S.-C. Spin-Orbital Texture in Topological Insulators. Phys. Rev. Lett. 111, 066801 (2013). [PubMed]
- Jozwiak C. et al. . Photoelectron spin-flipping and texture manipulation in a topological insulator. Nature Phys. 9, 293 (2013).
- Cao Yue et al. . Mapping the orbital wavefunction of the surface states in three-dimensional topological insulators. Nature Phys. 9, 499 (2013).
- Zhu Z.-H.
et al. . Photoelectron Spin-Polarization Control in the Topological Insulator Bi
_{2}Se_{3}. Phys. Rev. Lett. 112, 076802 (2014). [PubMed] - Xie Z. et al. . Orbital-selective spin texture and its manipulation in a topological insulator. Nature Commun. 5, 3382 (2014). [PubMed]
- Maaß H. et al. . Spin-texture inversion in the giant Rashba semiconductor BiTeI. Nature Commun. 7, 11621 (2016). [PMC free article] [PubMed]
- Nechaev I. A. & Krasovskii E. E.
Relativistic k·p Hamiltonians for centrosymmetric topological insulators from
*ab initio*wave functions. Phys. Rev. B 94 201410(R) (2016). - Krasovskii E. E. Microscopic origin of the relativistic splitting of surface states. Phys. Rev. B 90, 115434 (2014).
- Shan W.-Y., Lu H.-Z. & Shen S.-Q. Effective continuous model for surface states and thin films of three-dimensional topological insulators. New J. Phys. 12, 043048 (2010).
- Fu L. & Kane C. L. Topological insulators with inversion symmetry. Phys. Rev. B 76, 045302 (2007).
- Förster T., Krüger P. & Rohlfing M.
*GW*calculations for Bi_{2}Te_{3}and Sb_{2}Te_{3}thin films: Electronic and topological properties. Phys. Rev. B 93, 205442 (2016). - Zhou B., Lu H.-Z., Chu R.-L., Shen S.-Q. & Niu Q. Finite Size Effect on Helical Edge States in a Quantum Spin-Hall System. Phys. Rev. Lett. 101, 246807 (2008). [PubMed]

Articles from Scientific Reports are provided here courtesy of **Nature Publishing Group**

PubMed Central Canada is a service of the Canadian Institutes of Health Research (CIHR) working in partnership with the National Research Council's national science library in cooperation with the National Center for Biotechnology Information at the U.S. National Library of Medicine(NCBI/NLM). It includes content provided to the PubMed Central International archive by participating publishers. |