Chengwang
Niu
a,
Ying
Dai
*a,
Meng
Guo
a,
Yandong
Ma
a,
Baibiao
Huang
a and
Myung-Hwan
Whangbo
b
aSchool of Physics, State Key Laboratory of Crystal Materials, Shandong University, Jinan 250100, People's Republic of China. E-mail: daiy60@sina.com
bDepartment of Chemistry, North Carolina State University, Raleigh, North Carolina 27695-8204, USA
First published on 22nd October 2012
The ternary topological insulator Bi2Te2Se is exceptional as it has a larger bulk resistivity than previously studied topological insulators. Here we explored the possibility of modifying the topological surface states and achieving the insulating massive Dirac fermion state in Bi2Te2Se on the basis of density functional calculations. Substitution of O for the outmost-layer Te leads to tunable surface states with an ideal Dirac cone. The co-substitution of Cr and O, as well as that of Mn and F, for the surface Bi and Te places the Dirac point inside the bulk band gap and opens a band gap at the Dirac point, hence creating the insulating massive Dirac fermion state. The co-substitution of magnetic and non-magnetic elements is a promising way of tuning the properties of topological insulators.
Recently, Chen et al. investigated the evolution of the band structure of Bi2Se3 during photon-assisted O surface doping, and demonstrated the raising of the Dirac point by O doping.6 Song et al.29 examined the possibility of an interfacial Dirac cone from a periodic array of alternating band and topological insulators, which showed that the surface Se atoms in every third van der Waals gap replaced by O in Bi2Se3 induce an interfacial and ideal Dirac cone. The robustness of topological order in Bi2Se3 against oxidation was revealed by Wang et al. based on first-principles calculations.30 These studies suggest that a promising way of tuning electronic structure is to replace the chalcogen atoms (Se, Te) on the surface of a TI with O atoms. However, the O-substitution alone cannot break the TRS and introduce a band gap. A magnetic element, Cr, which has the same stable oxidation state (+3) as does Bi, is known to induce ferromagnetism in Bi-based TI.31 In the present work we show on the basis of density functional calculations that an ideal Dirac cone can be realized in Bi2Te2Se by introducing O atoms on the outermost layers, and that the IMDF state can be created by co-substituting Cr and O (also Mn and F) atoms for the Bi and Te atoms, respectively, without introducing free carriers.
For our density functional calculations for the pristine and the substituted slabs, we employed the projector augmented wave method implemented in the Vienna ab initio simulation package33 using the generalized gradient approximation of Perdew–Burke–Ernzerhof for the exchange-correlation functional,34 the plane-wave cutoff energy of 400 eV, the energy convergence within 10−6 eV, and a set of 9 × 9 × 1 k points to sample the surface Brillouin zone. The effective on-site parameter (Ueff = U − J) proposed by Dudarev et al.35 with typical values36 of U = 4.5 eV and J = 0.5 eV on transition metal elements is adopted for Cr 3d and Mn 3d states. Spin–orbit coupling (SOC) is the key ingredient for creating the topological insulator phase and thus SOC is included in our calculations with the spin orientation of the magnetic ions taken along the c-axis unless mentioned otherwise. In our calculations the lattice constants (both the in-plane and out of plane) and the internal atomic coordinates were fully optimized. In the following discussion, our results refer to the 6QL-slab calculations, unless stated otherwise.
Fig. 1 The band structures of Bi2Te2Se calculated for (a) the pristine 6QL-slab, (b) the O-substituted 6QL-slab, (c) the O-substituted 6QL-slab with no vacuum layer, (d) the O-substituted 4QL-slab, (e) the Cr-substituted 6QL-slab, and (f) the Mn-substituted 6QL-slab. The inset in panel (a) shows the 2D Brillouin zone, and those in (b)–(f) show zoomed-in views of the bands of the highlighted area. The Fermi level is indicated by the dashed line at 0 eV. |
Fig. 2 Comparison of the local structures and the 3D band structures around the Γ point calculated for the pristine and O-substituted (111) surfaces of Bi2Te2Se: (a) the optimized local structure, (b) the conduction band and (c) the valence band of the pristine (111) surface. (d) The optimized local structure, (e) the conduction band and (f) the valence band of the O-substituted (111) surface. The numbers in (a) and (d) are bond lengths in Å. The constant-energy contour plots are given at the bottom in (b), (c), (e) and (f). |
Fig. 3 (a) The total and projected DOS plots calculated for the pristine and the O-substituted 6QL-slabs. To make the impurity states more visible, the values of the projected DOS values for the Te 5p and O 2p states are multiplied by a factor of 10. The vertical dashed lines denote the Fermi levels. All of the DOS plots are aligned by using the Te 5p states of the deepest QL layer. The inset shows the local structure of the O–Bi layer with 25% O substitution. (b) The relative distance between the Dirac point and the Fermi level (EF) as a function of the SOC strength (λ0 is the actual SOC strength). The inset shows the evolution of the Dirac point with different SOC strength. (c) The band structures calculated for the 6QL-slab with 0%, 25%, 50%, 75%, and 100% O-substitution in the outermost Te layers. The highest occupied and the lowest unoccupied bands are plotted to clearly show the evolution of the TSS states. The Fermi level is indicated by the dashed line at 0 eV. |
We now examine how the position of the Dirac point depends on the extent of the O substitution for Te by using the 6QL-slab of Bi2Te2Se. The band structures calculated for the 0%, 25%, 50%, 75%, and 100% O-substitutions of the outmost Te layers are summarized in Fig. 3(c). The pristine (111) surface (0% substitution) shows that the Dirac point is buried in the bulk valence band. For 25% substitution, a ‘V’-shaped Dirac cone forms at the top of the valence band with the Fermi level crossing the upper Dirac cone. For 50% substitution, the Fermi level no longer crosses the upper Dirac cone and the concave lower Dirac cone becomes convex. With greater than 50% substitution, the ideal Dirac cone, retaining the Kramers degeneracy at the Dirac point and the massless character, are obtained, and the position of the Dirac point at the Γ point moves upward with increasing the O coverage. As shown in Fig. 3(c), the position of the Dirac point does not depend monotonically on the O contents; initially it shifts downward and then upward. This is attributable to complex changes in O–Bi layers. For the 25% substitution a stronger O–Bi bonding occurs leading to a shorter O–Bi bond (2.342 Å), which lowers the O 2p states much more than the other cases (Fig. 3(a)). For the band structure of 4QL-slab that is the critical slab thickness to maintain the Dirac cone26 with 100% O-for-Te substitution, the major features of the Dirac cone are obtained but a tiny finite gap at the Dirac point appears as shown in Fig. 1(d). This is a direct result of the shorter Bi–O bond, and the coupling between the upper and the lower surface states becomes strong enough to open up this whole insulating gap.32,38,39
Now that the ideal Dirac cone with the in-gap Dirac point is obtained for the 6QL-slab with the O-substitution, its electronic structure can be modified further to create the IMDF state. For this purpose, we consider an additional Cr-substitution to break the TRS. First, we investigate the situation of Cr-substitution. The band structure of the (111) surface with the Cr-substitution, presented in Fig. 1(e), shows a band gap in the surface state near the original Dirac point but has one energy band crossing the Fermi level. A magnetic moment of 3.0 μB with out-of-plane spin orientation is introduced by Cr-substitution of Bi. A gap at the Dirac point appears because the surface electrons carry real spins and there are exchange interactions between the introduced magnetic impurities and the surface electrons. To study the magnetic order, two dopant Cr atoms are placed at two Bi sites in one QL with the smallest separation (about 2.83 Å). The total energies are calculated with SOC for the ferromagnetic (FM) and antiferromagnetic (AFM) spin arrangements (EFM and EAFM, respectively) to find that ΔE = EFM − EAFM is −159 meV (the calculations with different AFM configurations are checked and show the same result), which means that the FM state is the ground state. This is essential for many striking phenomena including the quantized anomalous Hall effect.9,10 To see the effect of co-substitution on the Dirac point, the band structure of the 6QL-slab with the Cr/O co-substitution is calculated and presented in Fig. 4(a). It shows a gap at the Dirac point, separating the upper and lower Dirac cones with the Fermi level located in both the surface and bulk band gaps, namely, the IMDF state is produced. It is interesting to note that a similar Dirac point at DP′ appears (Fig. 1(e) and Fig. 4(a)) although the magnetic impurity has exchange interactions with the surface electrons. This is very similar to the chemisorption of nonmagnetic impurity on the topological surface of Bi2Se3.30 The Cr-substitution only changes the electronic state of the topmost surface, and opens a gap for the corresponding Dirac point, while the electronic state at the other surface is unchanged and its Dirac point remains. As shown in Fig. 5(a) and (b), the charge density associated with the bands forming the Dirac point shows that the Dirac point originates mainly from the O and Se atoms on the outermost QL for the O substitution, but mainly from the Te and Bi atoms on the second QL for the Cr/O co-substitution. But for the bulk case with the Cr/O co-substitution, as shown in Fig. 4(a), a band (C1) formed above the Fermi level. The band decomposed charge density analysis show that the band originates mainly from the doped Cr atom and its first-shell neighbors (Fig. 5(c)).
Fig. 4 The electronic structures of Bi2Te2Se with the Cr/O and Mn/F co-substitutions. The band structures calculated for the (a) Cr/O co-substituted and (b) Mn/F co-substituted 6QL-slab with SOC, and their corresponding bulk band projections are represented by the shaded areas. The total and projected DOS plots for the 6QL-slab (c) with the Cr/O co-substitution and (d) with the Mn/F co-substitution without SOC, the shaded areas indicate the occupied states. The Fermi level is indicated by the dashed line at 0 eV. The insets in panels (c) and (d) show the side view of the outermost QLs with the Cr/O and Mn/F co-substitutions, respectively. |
Fig. 5 The charge density associated with the bands forming the Dirac point of 6QL-slab Bi2Te2Se with (a) the O-substitution and (b) the Cr/O co-substitution, and that associated with the band C1 of bulk Bi2Te2Se with (c) Cr/O co-substitution and (d) Mn/F co-substitution. |
To examine the magnetic property, the total DOS plot and the projected DOS plots calculated without SOC for the Cr 3d and O 2p states for the 6QL-slab with the Cr/O co-substitution are plotted in Fig. 4(c). The Cr 3d states are strongly spin-polarized, resulting in a high-spin d3 configuration with 3.0 μB per Cr. The up-spin 3d states of each Cr3+ ion located at an octahedral site are split into t2g↑ < eg↑, and the energy gap between the t2g↑ and eg↑ levels is large enough to accommodate the upper and lower Dirac cones inside this gap. Similar to the Fe adatoms on the surface of the TI Bi2Se3,40 calculations of the energy of the Cr/O co-substituted system as a function of the Cr spin-moment direction show that the ‖c spin-orientation is more stable than the ⊥c spin-orientation by 0.32 meV per Cr. Then the different degrees of Cr substitution for Bi with 100% O-for-Te substitution are investigated, and the IMDF states are formed for all cases (see Fig. 6).
Fig. 6 The band structures calculated for the Cr/O co-substituted 6QL-slab of Bi2Te2Se with (a) 25%, (b) 50%, and (c) 75% Cr-substitution and 100% O-substitution. The Fermi level is indicated by the dashed line at 0 eV. |
In the following, we show the realization of IMDF states without the ideal Dirac cone. As found in previous studies,41,42 the introduction of a Mn atom in Bi2Te2Se leads to a high-spin Mn3+ ion with moment 4.0 μB indicating the (t2g↑)3(eg↑)1 configuration for the Mn3+ ion, which has one hole in the up-spin 3d-states of Mn3+. For the case when two Mn ions have the smallest separation (about 2.76 Å), the FM spin arrangement is more stable than the AFM spin arrangement by 264 meV. Fig. 1(f) shows the band structure of the 6QL-slab with the Mn-substitution. As in the case of the 6QL-slab with the Cr-substitution, a band gap appears in the surface state near the Dirac point of the pristine surface. A new Dirac point occurs and the dispersion of the new surface states still keeps the nearly perfect ‘V’ shape. The hole of the up-spin 3d-states of Mn3+ is filled by the F/Mn co-substitution, and the magnetic moment increases to 5.0 μB. As shown in Fig. 4(b), this co-substitution introduces a band gap between the upper and lower Dirac cones of the pristine surface due to the loss of the TRS, and lowers the Fermi level so that it is now positioned inside the bulk and surface band gaps, with the similar Dirac point at DP′ and band (C1) appearing. In the substitution of F for Te, an F− ion replaces a Te2− ion. Each F atom requires one less electron than does Te from the Bi2Te2Se lattice, so the extra electron is transferred to a Mn3+ ion converting it to Mn2+ with electron configuration (t2g↑)3(eg↑)2. Thus, a gap between the spin-up and spin-down 3d-states of Mn2+ accommodates the upper and lower Dirac cones, as can be seen from Fig. 4(d) (collinear calculations without including SOC). As a result, the Mn/F co-substitution also induces an IMDF state into Bi2Te2Se. In addition, for the Mn/F co-substituted system, the ‖c spin-orientation is calculated to be more stable than the ⊥c spin-orientation (by 0.53 meV per Mn). We note that substitution of Fe for Bi, which induces the same magnetic moment as Mn/F co-substitution, cannot introduce the IMDF state. Therefore, the formation of the IMDF state is likely set by the TRS breaking disorder potential introduced on the surface and the interaction of magnetic atoms with the surface atoms.
Finally, we examine the feasibility of the Cr/O and Mn/F co-substitutions by calculating the defect formation energy Eform using the following equations:43
Eform(M) = E(M-substituted) − E(pristine) + nBi(μBi − μM), |
Eform(M/L) = E(M/L co-substituted) − E(pristine) + nBi(μBi − μM) + nTe(μTe − μL), |
Substitution | E f (eV) | |||
---|---|---|---|---|
Te-poor | Te-rich | M tot (μB) | ||
a The bulk calculations employed the 60-atom supercell, and the surface calculations the 6QL-slab. | ||||
Bulk | Cr | 0.823 | 0.082 | 3.0 |
Cr/O | −0.352 | −0.352 | 3.0 | |
Mn | 1.782 | 1.041 | 4.0 | |
Mn/F | −2.721 | −2.721 | 5.0 | |
Surface | Cr | 1.887 | 1.146 | 3.0 |
Cr/O | 0.039 | 0.78 | 3.0 | |
Mn | 1.904 | 1.163 | 4.0 | |
Mn/F | −2.385 | −2.385 | 5.0 |
Fig. 7 The difference of formation energies between M-substituted and M/L co-substitution, ΔEform = Eform(M/L) − Eform(M), in the bulk Bi2Te2Se and on the (111) surface of Bi2Te2Se. The bulk calculations employed the 60-atom supercell, and the surface calculations the 6QL-slab. |
(1) Bi2Te2Se possesses a highly insulating bulk band, but the Dirac point is buried in the bulk valence band. The O-substitution for the surface layer Te can raise the Dirac point to have an ideal Dirac cone with the in-gap Dirac point.
(2) The substitution Cr for Bi can induce a magnetic moment of 3.0 μB with out-of-plane spin orientation, break the time-reversal symmetry and open a gap at the Dirac point. Together with the O-substitution induced ideal Dirac cone, the IMDF state is formed by the Cr/O co-substitution.
(3) The introduction of F atom in Mn-doped Bi2Te2Se leads to an increase in the out-of-plane magnetic moment and a change in the interaction between the Mn and surface atoms. As a result, the Mn/F co-substitution also induces an IMDF state into Bi2Te2Se.
The co-substitution of magnetic and non-magnetic elements into topological insulators is a promising way of tailoring their electronic structures and properties. Recent experimental progresses have shown that various impurity elements (such as Cr, Mn, and O) can be substituted into the Bi-based topological insulators.6,41,44,45 The predictions of our proposed co-substitution, which can be prepared by the current technical capability of materials engineering, can be experimentally verifiable.
This journal is © The Royal Society of Chemistry 2013 |