Magnetizing topological surface states of Bi2Se3 with a CrI3 monolayer

See allHide authors and affiliations

Science Advances  31 May 2019:
Vol. 5, no. 5, eaaw1874
DOI: 10.1126/sciadv.aaw1874


To magnetize surfaces of topological insulators without damaging their topological feature is a crucial step for the realization of the quantum anomalous Hall effect (QAHE) and remains as a challenging task. Through density functional calculations, we found that adsorption of a semiconducting two-dimensional van der Waals (2D-vdW) ferromagnetic CrI3 monolayer can create a sizable spin splitting at the Dirac point of the topological surface states of Bi2Se3 films. Furthermore, general rules that connect different quantum and topological parameters are established through model analyses. This work provides a useful guideline for the realization of QAHE at high temperatures in heterostructures of 2D-vdW magnetic monolayers and topological insulators.


Topological insulators (TIs) are emergent quantum materials that have nontrivial bandgaps in their bulks along with topologically protected surface or edge states at boundaries (1). For the design of the next-generation spintronic devices (2, 3), especially those based on the quantum anomalous Hall effect (QAHE), it is crucial to find efficient ways to magnetize their topological surface states (TSSs) while maintaining their topological features (46). The conventional approach is doping magnetic ions into TIs (711), and the QAHE has been successfully realized in Cr- or V-doped (Bi, Sb)2Te3 thin films (911). However, it is still very challenging to control the distribution and magnetic order of dopants in TIs, and the critical temperature (Tc) for the observation of QAHE is extremely low (30 mK) (9). A promising alternative way to magnetize TSSs is through the interfacial magnetic proximity effect by putting three-dimensional (3D) magnetic insulators on TIs (1218). Unfortunately, the interfacial hybridization in most of these heterostructures is too strong, and TSSs are either damaged or shifted away from the Fermi level (12, 14, 19). To this end, recently discovered 2D-vdW ferromagnetic monolayers (MLs) such as CrI3 (20) appear to offer an optimal way to magnetize the TSSs of TIs. First, the Curie temperature of the CrI3 ML is as high as 45 K (20), and it is hence conceivable that the Tc of the QAHE in CrI3/TI heterostructures can be far higher than that in the Cr-doped (Bi, Sb)2Te3 thin films (9). Furthermore, the CrI3 ML is a semiconductor, and thus, the transport properties of TIs should be preserved. Nevertheless, magnetic ions in most 2D-vdW magnetic materials such as Cr3+ ions in CrI3 are typically covered by nonmagnetic layers in both sides; it is thus questionable whether their spin polarization can be sensed by TSSs, although the latter have a fairly large spatial extension.

Here, we report results of systematic computational studies based on the density functional theory (DFT) that suggest the possibility of using a CrI3 ML to magnetize the TSSs of a prototypical 3D TI: Bi2Se3 (BS). We build up inversion-symmetric CrI3/BS/CrI3 heterostructures with a varying thickness for the BS films, from three to seven quintuple layers (QLs). We find that the TSSs of BS can be effectively magnetized by the CrI3 ML in all cases. However, we show that CrI3/BS/CrI3 becomes a Chern insulator only when the BS film is thicker than five QLs. This results from the competition between the exchange field from CrI3 and the remaining interaction between two surfaces of the BS film. This work reveals the subtleness of designing topological spintronic materials and provides a general guidance for the realization of the QAHE at high temperatures by combining 3D TIs with 2D-vdW magnetic MLs (20, 21).


Electronic properties of CrI3/BS/CrI3 heterostructures

CrI3 ML is a robust 2D semiconducting ferromagnet. It has a perpendicular magnetic anisotropy and a reasonably high Curie temperature (45 K) (20). Structurally, CrI6 octahedrons form a honeycomb lattice. The optimized lattice constant of the pristine CrI3 ML is aCrI3 = 7.04 Å, consistent with the previous theoretical result (22). This value is only 2.5% smaller than the size of the 3×3 BS supercell (7.22 Å), and hence, we stretch the CrI3 ML so as to use a manageable unit cell for the simulation of CrI3/BS/CrI3. We find that all main properties of the CrI3 ML are not substantially affected by the small lattice stretch (see fig. S1 and table S1).

There are three possible highly symmetric alignments between the 3×3 BS supercell and the CrI3 ML, i.e., with Cr ions taking the Se, hollow, or Bi sites on BS (Fig. 1A), respectively. The calculated binding energies suggest that Cr3+ cations prefer to sit on the top of Se2− anions (see fig. S2). The optimized interlayer distances between the CrI3 ML and the BS surface are in a range from 3.07 to 3.12 Å, depending on the thickness of the BS film (see fig. S3). Note that the energy change is very small as we shift the CrI3 ML in the lateral plan, so it is crucial to impose the inversion symmetry during the structural optimization. Here, we only discuss electronic and magnetic properties of the most stable configuration.

Fig. 1 Atomic structure and charge difference in CrI3/BS/CrI3.

(A) Side view of the most stable CrI3/BS/CrI3. The green, red, pink, and blue balls are for I, Cr, Se, and Bi atoms, respectively. Dashed lines show the alignments between atoms in CrI3 and BS layers; d0 denotes the optimized vdW gap. (B) Real-space distribution of the charge difference Δρ = ρtotal − ρBS − ρCrI3 and (C) planar-averaged spin density Δσ = ρ − ρ in the interfacial region of the CrI3/6QL-BS/CrI3 heterostructure.

To shed some light on the interaction across the CrI3/BS interface, we plot the charge density difference Δρ of CrI3/BS/CrI3, using six-QL BS as an example in Fig. 1B. The charge redistribution only occurs within a couple of atomic layers in BS, and the magnitude of Δρ is small. In particular, there is no observable net charge transfer between CrI3 and BS. This is understandable since the bandgap of the freestanding CrI3 ML is rather wide (larger than 1 eV). The interfacial Se atoms acquire a small but meaningful magnetic moment of −0.003 μB, which aligns antiparallelly with the magnetic moments of Cr3+ ions. The planar-averaged spin density Δσ (Fig. 1C) shows that the negative spin polarization penetrates through the top half of the first QL of BS, following the direction of spin polarization of the iodine atoms.

The calculated band structures of CrI3/BS/CrI3 with different thicknesses of BS (four, five, and six QLs) are shown in Fig. 2 (A to C, respectively). One important feature is that bands near the Fermi level are predominantly from BS, and CrI3 states lie either 0.7 eV above the Fermi level or at least 0.3 eV below the Fermi level. As a result, the Dirac cone feature of TSSs are well maintained, or more explicitly, (i) all three cases have bandgaps of several millielectron volts at the Γ point and (ii) the spin degeneracies of TSSs are lifted, indicating that the TSSs of BS are magnetized by the CrI3 ML. These results suggest the suitability of using vdW magnetic MLs to realize QAHE, instead of using conventional ferromagnetic or antiferromagnetic films that damage TSSs (14, 15, 19, 23).

Fig. 2 Band structures of CrI3/BS/CrI3.

DFT-calculated band structures of (A) CrI3/4QL-BS/CrI3, (B) CrI3/5QL-BS/CrI3, and (C) CrI3/6QL-BS/CrI3. Insets in (A), (B), and (C) are the fine band structures around the Fermi level. Colors in the main panels indicate the weights of bands from BS (blue) and CrI3 (yellow). Colors in insets indicate the spin projections.

Topological properties of CrI3/BS/CrI3 heterostructures

To examine whether the magnetized TSSs of BS by CrI3 MLs are topologically nontrivial, we take CrI3/4QL-BS/CrI3 and CrI3/5QL-BS/CrI3 as examples and calculate their Chern numbers CN by integrating the Berry curvature in the first Brillouin zone with the Wannier90 package (see figs. S4 and S5). The Chern numbers of CrI3/4QL-BS/CrI3 and CrI3/5QL-BS/CrI3 are CN4QL=0 and CN5QL=1, respectively. This indicates that the former is a normal insulator but the latter is a Chern insulator. As shown in Fig. 3A, the QAHE state of CrI3/5QL-BS/CrI3 is further confirmed by the presence of one chiral edge state that connects the valance and conduction bands.

Fig. 3 Topological properties of CrI3/BS/CrI3.

(A) Chiral edge state of the Chern insulator CrI3/5QL-BS/CrI3 ribbon. (B) Dependence of Chern numbers and gaps of CrI3/BS/CrI3 on the number (N) of QLs of BS. The inset shows the dependence of bandgaps of CrI3/6QL-BS/CrI3 on the CrI3-BS interlayer distance d. d0 is the optimized interfacial distance. (C) and (D) show the fitting parameters Δ and M in different CrI3/BS/CrI3 heterostructures, respectively. The horizontal dashed line in (B) shows the position of 3.5 meV.

We use a low-energy effective four-band Hamiltonian to thoroughly study the topological properties of all CrI3/BS/CrI3 films. Note that states of CrI3 MLs are far away from the Fermi level and that the effect of CrI3 ML on TSSs of BS can be represented by an exchange field (HZeeman) and an interfacial potential (HInterface). Hence, the low-energy effective four-band Hamiltonian (2, 5, 6) within the basis set of {∣t, ↑ 〉, ∣t, ↓ 〉, ∣b, ↑ 〉, and ∣b, ↓ 〉} isH(kx,ky)=Hsurf(kx,ky)+HZeeman(kx,ky)+HInterface(kx,ky)=A(kx2+ky2)+[0ivFkMk0ivFk+00MkMk00ivFk0MkivFk+0]+[Δ0000Δ0000Δ0000Δ]+[Vp0000Vp0000Vp0000Vp](1)Here, t and b denote the top and bottom surface states, respectively; ↑ and ↓ represent the spin up and down states, respectively; vF and k± = kx ± iky are the Fermi velocity and wave vectors, respectively; and Mk = MB(kx2 + ky2) describes the coupling between the top and bottom TSSs. Note that M represents the gap produced by this coupling rather than magnetization as conventionally used, Δ is the exchange field from CrI3, and Vp represents the magnitude of the asymmetric interfacial potential. Although Vp vanishes in the models with inversion symmetry, the situation in experiments can be changed by the misalignment of two CrI3 MLs and by the presence of substrates and cover layers. Therefore, we keep the Vp term to hold the generality of our discussions. As shown in fig. S6 and table S2, the DFT bands of CrI3/BS/CrI3 around the Γ point can be fitted well by this Hamiltonian, indicating its applicability to these systems. The Berry curvature Ω(k) and Chern number CN are calculated on the basis of the formulas in (24).

There is a clear phase boundary in Fig. 3B for the topological feature of CrI3/BS/CrI3, i.e., CN = 0 (CN = 1), when BS is thinner (thicker) than five QLs, consistent with the Chern numbers obtained by Wannier functions in CrI3/BS/CrI3 with four- and five-QL BS. Therefore, CrI3/BS/CrI3 heterostructures with less than five QLs of BS are normal insulators, and five QLs of BS are required to realize the QAHE in experiments. The calculated bandgaps of CrI3/5QL-BS/CrI3, CrI3/6QL-BS/CrI3, and CrI3/7QL-BS/CrI3 are 2.3 meV (27 K), 3.6 meV (42 K), and 5.9 meV (68 K), respectively. Equation 1 shows that its bandgap is 2∣M − Δ∣ (here, M > 0 and Δ> 0). On the basis of the data in Fig. 3 (C and D), we see that the increase of bandgap with the BS thickness mainly results from the monotonic decrease of M. Moreover, the nontrivial bandgap of CrI3/BS/CrI3 noticeably increases with the reduction of d (see the inset in Fig. 3B for CrI3/6QL-BS/CrI3), indicating that it is beneficial to make d smaller through an external pressure for the realization of QAHE. Considering that the ML CrI3 has a Tc of 45 K (20) and no other complex factors, such as uncontrollable distribution of dopants, local magnetic ordering, and charge transfer, are involved in these systems, we believe that QAHE should be observable in CrI3/BS/CrI3 heterostructures up to a few tens of kelvins, much higher than the temperature achieved with the doping approach (911).

In general, we find rules that govern topological properties of CrI3/BS/CrI3, i.e., (i) CN = 1 when Δ2 > M2 and (ii) CN = 0 when Δ2 < M2. This is consistent with the phase diagram in (5) and gives simple ways to manipulate the topological properties of systems with a 2D magnetic ML on TI. As shown in the Fig. 3C, the exchange field Δ slightly depends on the thickness of the BS film and fluctuates around 3.5 meV. This result is understandable since exchange field Δ mainly depends on the interfacial interaction between CrI3 and BS. As expected, M has a strong dependence on the thickness of the BS film and exponentially decreases as the BS film becomes thicker (Fig. 3D). Note that the very small M = 0.5 meV in CrI3/7QL-BS/CrI3 is consistent with the experimentally observed gap closing in BS films thicker than six QLs (25). Since Vp is zero in CrI3/BS/CrI3 films with the inversion symmetry, the Hamiltonian equation 1 can be rewritten in terms of the basis of {∣ + , ↑ 〉, ∣ − , ↓ 〉, ∣ + , ↓ 〉, and ∣ − , ↑ 〉} with ±,=(t,±b,)/2 and ±,=(t,±b,)/2 as (2, 5)H˜(kx,ky)=A(kx2+ky2)+[H+(kx,ky)00H(kx,ky)](2)In Eq. 2, H±(kx, ky) = vFkyσxvFkxσy + [M ± Δ − B(kx2 + ky2)]σz and σx, y, z are Pauli matrices. From this Hamiltonian, we can get that the condition of Δ2 > M2 leads to only one of the band inversions, either in H+(kx, ky) or in H(kx, ky), and a Chern number CN = Δ/∣Δ∣ (see table S3). This is understandable because the exchange splitting Δ can overcome the coupling-induced gap M in this condition and cause the band inversion of a pair of spin subbands (26). In CrI3/BS/CrI3 with four-QL or thinner BS, we have Δ2 < M2 and hence CN = 0. In contrast, thicker CrI3/BS/CrI3 films satisfy the condition of Δ2 > M2 and hence become topologically nontrivial with CN = 1. Even in the case that the top and bottom TSSs are completely decoupled, i.e., both M = 0 and B = 0, the exchange field Δ from the CrI3 ML still leads to the QAHE in CrI3/BS/CrI3 as a result of the unique half-integer quantum Hall of the 2D massive Dirac Fermion in the two surfaces of the 3D TI (5, 6, 27).

To explore the effect of inversion symmetry breakdown, we investigate the evolution of band structure of CrI3/6QL-BS/CrI3 as the asymmetric interface potential Vp increases. The original topologically nontrivial bands (Fig. 4A) gradually change to M- and W-shaped bands (Fig. 4C) for Vp > 2.8 meV. Similar M- and W-shaped bands induced by Vp are also reported in a previous work (28). The bottom panel of Fig. 4D shows the Berry curvature of the occupied bands as Vp = 5.0 meV. One may see regions with both negative and positive Berry curvatures around the Γ point. This is different from the Berry curvature for the Vp = 0 case shown in the top panel of Fig. 4D, where only a huge positive Berry curvature appears. With the presence of Vp, the criterion for the topological phase transition can be extended as follows (5): (i) When Δ2>M2+Vp2, CN = 1, and (ii) when Δ2<M2+Vp2, CN = 0. Note that it is counterintuitive that M- and W-shaped bands induced by the asymmetric interface potential Vp is topologically trivial. This shows the subtleness of identifying a topological state of CrI3/BS/CrI3 heterostructures or probably all analogous material systems. Note that the presence of Vp quickly reduces the topologically nontrivial gap in the Chern insulator region (Fig. 4E). This shows the importance of maintaining the inversion symmetry of CrI3/BS/CrI3 heterostructures for the realization of QAHE and other exotic quantum properties (26).

Fig. 4 Effect of Vp on the topological properties of CrI3/6QL-BS/CrI3.

(A) to (C) show band evolutions with Vp increasing from 0 to 5.0 meV. (D) Berry curvatures of occupied bands around the Γ point for Vp = 0 (A, top) and Vp = 5.0 meV (C, bottom). The large positive Berry curvatures in the bottom panel are highlighted by black circles. (E) Dependence of Chern numbers and gaps on Vp Δ=3.1 meV and M = 1.4 meV are adopted in these calculations.


Magnetizing TSSs of TIs is of great importance as it gives rise to QAHE and many other novel phenomena, especially the hallmark of 3D TIs—the topological magnetoelectrical effect (27). We noted that large spin splitting might be produced in TSSs through δ codoping of Mn and Se(Te) in the topmost QL of TIs (29), but this may involve issues regarding the distribution and magnetic ordering of dopants. As magnetic insulators such as MnSe and EuS damage TSSs or shift TSSs away from the Fermi level (12, 14, 19), the emergent 2D-vdW magnetic MLs (20, 21) offer unique opportunities to magnetize TSSs through the proximity effect as demonstrated in the present work.

In summary, using CrI3/BS/CrI3 heterostructures as modeling systems, we demonstrated that the emergent 2D-vdW magnetic MLs can effectively magnetize TSSs of TIs and also maintain their topological characteristics around the Fermi level. Furthermore, our analyses with an effective four-band Hamiltonian verify that CrI3/BS/CrI3 are Chern insulators when the BS film is five QLs or thicker. Their bandgaps are a few millielectron volts, and we believe that the QAHE is observable in these heterostructures at a temperature of a few tens of kelvins. As the family of 2D-vdW magnetic MLs steadily expands, it is foreseeable that even a higher-temperature QAHE can be achieved. According to discussions above, promising 2D-vdW magnetic MLs should have strong magnetization, large spatial extension of spin density, and shorter interfacial distance to TI surfaces.


We used the Vienna Ab initio Simulation Package at the level of the generalized gradient approximation (30, 31) in this work. The projector-augmented wave pseudopotentials were adopted to describe the core-valence interaction (32, 33), and the energy cutoff for the plane-wave expansion was set to 500 eV (31). As shown in Fig. 1A, we constructed inversion-symmetric CrI3/BS/CrI3 slab models with a3×3 BS supercell in the lateral plane to match the CrI3 ML. The vacuum space between adjacent slabs was set to 15 Å. Atomic structures were fully optimized with a criterion that requires the force on each atom being less than 0.01 eV/Å. To correctly describe the weak interaction across CrI3 and BS layers and the strong relativistic effect in Bi and I atoms, we included the nonlocal vdW functional (optB86b-vdW) (34) and the spin-orbit coupling (SOC) term in the self-consistent iterations. Furthermore, the local-spin-density approximations plus U (LSDA+U) method (35), with the on-site Coulomb interaction U = 3.0 eV and the exchange interaction J = 0.9 eV, was adopted to take the strong correlation effect of Cr 3d electrons into account. Topological properties of CrI3/BS/CrI3 slabs were studied by using effective Hamiltonians and the Wannier90 package (36). Test calculations indicate that the topological properties of CrI3/BS/CrI3 are very robust against the choices of U (fig. S7).


Supplementary material for this article is available at

Table S1. Magnetic moment, bandgap, magnetocrystalline anisotropy energy, and Heisenberg exchange interactions J1 and J2.

Table S2. Parameters of the effective four-band model (see Eq. 1 in the main text) are used to fit the DFT-calculated band of CrI3/BS/CrI3.

Table S3. Analysis of band inversions and the Chern number CN with respect to parameters Δ, M, and B.

Fig. S1. DFT + SOC + U–calculated band structure of the CrI3 ML with pristine and stretched lattice constants.

Fig. S2. Top views of the three highly symmetric alignments between CrI3 and BS in CrI3/BS/CrI3 heterostructures.

Fig. S3. Dependence of the binding energies and the vdW gaps of CrI3/BS/CrI3 on the number (N) of QL of BS.

Fig. S4. Topological properties of CrI3/4QL-BS/CrI3.

Fig. S5. Topological properties of CrI3/5QL-BS/CrI3.

Fig. S6. DFT-calculated (black lines) and fitted band structures (red dashed lines) of CrI3/BS/CrI3 based on the effective four-band model (see Eq. 1 in the main text).

Fig. S7. Effect of U on the topological properties of CrI3/BS/CrI3.

This is an open-access article distributed under the terms of the Creative Commons Attribution-NonCommercial license, which permits use, distribution, and reproduction in any medium, so long as the resultant use is not for commercial advantage and provided the original work is properly cited.


Acknowledgments: DFT calculations were performed on parallel computers at NERSC. Funding: Work was supported by DOE-BES (grant no. DE-FG02-05ER46237). Author contributions: R.W. conceived the idea of this study. Y.H. performed DFT calculations. J.K. calculated the edge states. Y.H. and R.W. wrote the manuscript. All authors commented on the manuscript. Competing interests: The authors declare that they have no competing interests. Data and materials availability: All data needed to evaluate the conclusions in the paper are present in the paper and/or the Supplementary Materials. Additional data related to this paper may be requested from the authors.

Stay Connected to Science Advances

Navigate This Article