## Pursuing High-Temperature Quantum Anomalous Hall Effect in MnBi Te /Sb Te 2 4 2 3 Heterostructures

Shifei Qi, 1, 2 Ruiling Gao, 3, 2 Maozhi Chang, 3, 2 Yulei Han, 1 and Zhenhua Qiao 1, ∗ 1 ICQD, Hefei National Laboratory for Physical Sciences at Microscale, CAS Key Laboratory of Strongly-Coupled Quantum Matter Physics, and Department of Physics, University of Science and Technology of China, Hefei, Anhui 230026, China 2 Department of Physics, Hebei Normal University, Shijiazhuang, Hebei 050024, China 3 Institute of Materials Science, Shanxi Normal University, Linfen, Shanxi 041004, China

(Dated: December 2, 2019)

Quantum anomalous Hall effect (QAHE) has been experimentally realized in magnetically-doped topological insulators or intrinsic magnetic topological insulator MnBi 2 Te 4 by applying an external magnetic field. However, either the low observation temperature or the unexpected external magnetic field (tuning all MnBi 2 Te 4 layers to be ferromagnetic) still hinders further application of QAHE. Here, we theoretically demonstrate that proper stacking of MnBi 2 Te 4 and Sb Te 2 3 layers is able to produce intrinsically ferromagnetic van der Waals heterostructures to realize the high-temperature QAHE. We find that interlayer ferromagnetic transition can happen at T C = 42 K when a fivequintuple-layer Sb 2 Te 3 topological insulator is inserted into two septuple-layer MnBi 2 Te 4 with interlayer antiferromagnetic coupling. Band structure and topological property calculations show that MnBi Te /Sb Te /MnBi Te 2 4 2 3 2 4 heterostructure exhibits a topologically nontrivial band gap around 26 meV, that hosts a QAHE with a Chern number of C = 1. In addition, our proposed materials system should be considered as an ideal platform to explore high-temperature QAHE due to the fact of natural charge-compensation, originating from the intrinsic n-type defects in MnBi 2 Te 4 and p-type defects in Sb 2 Te . 3

Introduction-. As a typical topological state, quantum anomalous Hall effect (QAHE) has been firstly observed in magnetic-element doped topological insulator (TI) thin films [1-5]. Due to its robustness against local perturbations, the corresponding chirally propagating edge mode has immense potential in future lowpower or dissipationless electronic device applications [68]. QAHE was initially predicated in a honeycomb lattice toy model by F. M. Haldane in 1988 [9]. After then, much effort has been made to exploring new systems for realizing the QAHE in these related systems [10-21]. Because of the inherently strong spin-orbit coupling in Z 2 topological insulators [22, 23], it was predicted that ferromagnetic TI thin films can give rise to QAHE. One convenient approach to induce magnetism in TIs is to dope magnetic elements [13, 24-29]. This has indeed led to the first experimental observation of QAHE [1]. So far, all the experimentally observed QAHEs in doped TIs were achieved at extremely low temperatures, typically ∼ 30 mK, making drastically increasing the QAHE observation temperature a daunting challenge both fundamentally and for potential applications.

control of interfaces, van der Waals (vdW) layered magnetic materials exhibit great advantage in designing highquality FMI/TI heterostructures.

Aside from the magnetic-element doped TIs, sandwiched TI heterostructures between two ferromagnetic insulator (FMI) layers were considered as another practical platform, because of the more 'clean' ferromagnetism with higher Curie temperature than that in the magnetic-element doped TIs [30, 31]. Unfortunately, due to the poor interfaces, the observed anomalous Hall resistance cannot reach the quantized value, though several FMI/TI heterostructures have been successfully fabricated in experiments [32-38]. Considering the better

Recently, MnBi Te 2 4 (vdW compounds) was successfully synthesized and confirmed as an intrinsic magnetic TI in experiments [40, 41]. It possesses layered structure with a triangle lattice, as displayed in Fig. 1(a) and 1(b). So far, rich topological states were proposed in MnBi 2 Te 4 due to its characteristics of intralayer ferromagnetic and interlayer antiferromagnetic exchange interactions [3942]. In particular, QAHE in MnBi Te 2 4 system has been observed at 4.5 K by applying an external magnetic field to convert the initial antiferromagnetism between layers to be ferromagnetism [43]. This is an undoubtedly obvious enhancement of QAHE observation temperature. To realize QAHE without applying any external field in MnBi Te 2 4 systems, interlayer antiferromagnetic coupling is a big obstacle. It is known that single layer MnBi Te 2 4 is a trivial FMI. Thus, one can construct the vdWMnBi Te /TI heterostructures with atomically per2 4 fect interface and intrinsic magnetism. We noted that MnBi Te /Bi Te 2 4 2 3 (not a TI) heterostructures have been successfully fabricated in experiment [44]. Although similar MnBi Te /Bi Te 2 4 2 3 heterostructure was theoretically predicted to be good candidate for realizing QAHE [45], the naturally existing n-type defects in both MnBi Te 2 4 and Bi Te 2 3 make their heterostructures metallic. On the contrary, considering the intrinsic p-type defects in Sb Te 2 3 (a 3D TI), MnBi Te /Sb Te 2 4 2 3 heterostructure should be charge compensated, which resolves the obstacle for realizing QAHE in MnBi Te 2 4 system.

In this Letter, we theoretically show that high-

FIG. 1. (a) Side view and (b) top view of crystal structures of ABC-ABC stacking heterostructure constructed by 1SL MnBi Te 2 4 and 5QL Sb Te ; (c) The corresponding top view 2 3 and (d) side view of crystal structure of ABC-ACB stacking case. Arrows (P) indicate the direction of spontaneous electric polarization. Red-solid frames in panels (a) and (d) label different stacking of 1SL MnBi Te 2 4 and 5QL Sb Te . 2 3 (e) Schematics of six different heterostructures with MBT and ST representing MnBi Te 2 4 and Sb Te , respectively. 2 3

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00498v2-figure-caa34e12-713d-4c95-933c-df5a0b6e5a34.jpg)

temperature QAHE can be realized in MnBi Te /Sb Te 2 4 2 3 heterostructure without additional magnetic field. We first demonstrate that an interlayer ferromagnetic transition can occur at T C = 42 K when a five quintuplelayer (QL) Sb Te 2 3 is inserted between two septuple-layer (SL) MnBi Te 2 4 with interlayer antiferromagnetic coupling. Further calculation indicates that the sandwiched heterostructure is an intrinsic QAHE with a bulk gap of 26 meV and a nonzero Chern number of C = 1. Our study not only reveals the possibility of producing ferromagnetic MnBi Te 2 4 to remove the unexpected external magnetic field in observing QAHE, but also provides an experiment-friendly sandwiched system to realize hightemperature QAHE.

Methods and Systems-. Our first-principles calculations were performed by using the projected augmentedwave method [46] as implemented in the Vienna Ab-initio Simulation Package (VASP) [47]. The generalized gradient approximation (GGA) of Perdew-Burke-Ernzerhof (PBE) type was used to treat the exchange-correlation interaction [48]. A 1 × 1 MnBi Te /Sb Te 2 4 2 3 supercell is chosen in our study. For thin film calculations, the Sb Te 2 3 film thickness was chosen to be 5QL. A vacuum buffer space of 20 Awas used to prevent the ˚ coupling between adjacent slabs. The kinetic energy cutoff was set to be 350 eV. During structural relaxation, all atoms were allowed to relax until the Hellmann-Feynman force on each atom is smaller than 0.01eV/ ˚ . The BrillouinA zone integration was carried out by using 9 × 9 × 1 Monkhorst-Pack grids for different heterostructure systems. Unless mentioned otherwise, spin-orbit coupling and the GGA+U method were used with the on-site repulsion parameter U=3.00 and 5.34 eV [49] and the exchange parameter J=0 eV, where U is for the more localized 3 d orbitals of Mn atoms in all calculations. Different U values provide similar results in our calculations. In addition, the vdW interactions were considered in our numerical calculations [50]. The Curie temperature T C was estimated within the mean-field approximation k B T C = 2 3 / Jx [51],where k B is the Boltzmann constant, x is the dopant concentration, and J is the exchange parameter obtained from the total energy difference between ferromagnetic (FM) and antiferromagnetic (AFM) configurations in different heterostructures.

MnBi Te 2 4 is composed of SL blocks stacked one by one along the [0001] direction via vdW interaction [40, 41] [see Fig. 1(a)]. The resulting magnetic order appears to be interlayer antiferromagnetic, where the ferromagnetic Mn layers of neighboring blocks are coupled in an antiparallel manner [40, 41]. Among various 3D TIs, Bi Te 2 3 should be the most optimal candidate to construct MnBi Te /TI heterostructures due to the 2 4 small lattice mismatch (1.3%). However, as observed in related experiments, MnBi Te 2 4 and Bi Te 2 3 are respectively electron conducting due to their intrinsic n-type defects. Thus, their heterostructures are metallic (i.e., the Fermi level is out of the band gap). Considering the intrinsic ptype conducting character of Sb 2 Te 3 (lattice mismatch is about 2.3%), the MnBi Te /Sb Te 2 4 2 3 heterostructure can be tuned to be charge compensated, which is an important prerequisite condition for realizing QAHE. Therefore, in our work, we focused on different heterostructures [See Fig. 1(e)] constructed by MnBi Te 2 4 and Sb Te . 2 3

Structural Stability-. Let us first study the structural stability of heterostructures. It is noted that the stacking of MnBi Te 2 4 and Sb Te 2 3 is different from the similar system MnBi Te /Bi Te 2 4 2 3 [45], where the ABC-ABC stacking is most favorable. Two different stacking heterostructures of MnBi Te 2 4 and Sb Te 2 3 are provided. One is in the ABC-ABC stacking [Fig. 1(a)], and the other is in the ABC-ACB stacking [Fig. 1(d)]. Their top views are respectively displayed in Figs. 1(b) and 1(c). Our numerical calculations show that the ABC-ACB stacking is more preferred based on the calculated result that the total energy of ABC-ACB stacking structure is -30.8 meV/unit cell lower than that of ABC-ABC stacking structure. The underlying reason is closely related to the dipole-dipole interaction between MnBi Te 2 4 and Sb Te . As highlighted in Figs. 1(a)and 1(d), the dipole2 3 moment directions of MnBi Te 2 4 and Sb Te 2 3 are opposite in the ABC-ABC stacking case, while they are the same in the ABC-ACB stacking one. All these indicate that

TABLE I. Total magnetic moments, magnetic order, energy difference between FM and AFM configurations ∆ E FM -AFM , and Curie temperature T C in 2SLs MnBi Te 2 4 and different MnBi Te 2 4 and Sb Te 2 3 heterostructures.

| Systems           | Moments/ µ B   | Magnetic Order   | ∆E/meV   | T C /K   |
|-------------------|----------------|------------------|----------|----------|
| 2SLs MBT          | 0 . 0          | AFM              | 9 . 3    | -        |
| MBT/ST            | 5 . 0          | FM               | -        | 14       |
| MBT/ST/MBT        | 10 . 0         | FM               | - 27 . 9 | 42       |
| ST/MBT/ST         | 5 . 0          | FM               | -        | 14       |
| MBT/ST/MBT/ST     | 0 . 0          | AFM              | 0 . 3    | -        |
| MBT/ST/MBT/ST/MBT | 5 . 0          | AFM              | 2 . 4    | -        |
| ST/MBT/ST/MBT/ST  | 0 . 0          | AFM              | 11 . 4   | -        |

the dipole-dipole interactions are respectively repulsive and attractive in ABC-ABC and ABC-ACB stacking. This finding can also be reflected by the larger interfacial distance in the ABC-ABC stacking than that in the ABC-ACB stacking [see Figs. 1(a) and 1(d)]. According to these finding, hereinbelow, we adopted the ABC-ACB stacking system to further investigate the magnetic properties, band structures, and topological properties.

Magnetic Properties-. As displayed in Tab. I, we first confirm the magnetic properties of 2SLs MnBi Te . 2 4 From our calculations, the 2SLs film indeed exhibits interlayer antiferromagnetism (i.e., 9.3 meV lower than that in ferromagnetic order). The local magnetic moment in 1SL is 5 µ B , agreeing with previous report [41], even if it is inserted into different heterostructures. Among all six studied heterostructures, only three of them are in the ferromagnetic ordering state. Especially for the MBT/ST/MBT sandwiched system, we found that a ferromagnetic phase transition occurs when a 5QL Sb Te 2 3 is inserted between two 2SL MnBi Te 2 4 thin films with interlayer antiferromagnetic coupling. This finding agrees well with a recent experimental result [44], i.e., the interlayer antiferromagnetic exchange interaction can be gradually weakened with increasing the separation between ferromagnetic layers in the MnBi Te /Bi Te 2 4 2 3 heterostructure and a ferromagnetic state could be established below 5K [44]. The total energy difference between FM and AFM states is -27.9 meV in the MBT/ST/MBT system and the estimated Curie temperature T C is around 42K. The other favorable FM state can also be obtained in MBT/ST, ST/MBT/ST, and odd-SL MBT/ST/MBT/ST/MBT systems, which are highly desirable candidates for realizing high-temperature QAHE. In fact, the magnetic ground states of MBT/ST, MBT/ST/MBT, and ST/MBT/ST heterostructure have been calculated with the inclusion of spin-orbit coupling, and the results indicate that out-of-plane magnetic structures are more preferred, which are respectively about 0.01, 0.06, and 0.03 meV(per Mn atom) lower than those of in-plane magnetic structures of MBT/ST, MBT/ST/MBT, and ST/MBT/ST heterostructures.

For the TI/magnetic insulators heterostructures, sur- face carrier-mediated interaction is found to be responsible for long-range ferromagnetism. In especially, Kou et al. found that the penetration depth of the surface states is a typical value of 2-3 QLs, and then they observe that fewer surface Dirac fermions will participate into the RKKY mediating process with larger TI layer depth, therefore weaking the surface-related RKKY coupling strength [52]. For our MBT/ST heterostructures, surface carrier-mediated interaction is responsible for the long-range ferromagnetism. In the MBT/ST/MBT system, as shown in the Figs. 2(a), the FM spin polarizations mainly appear at the interfaces of MBT and ST topological insulator and long-range FM coupling exists in this system. According to Kou et al's finding, the penetration depth of the surface states is 2-3 QLs. For 5QLs ST topological insulator, the penetration depths of surface states of the top and bottom layers are all 2-3 QLs. The Figs. 2(b) implies exchange interactions between different QLs. Hence, long-range FM coupling can be realized by surface carriers in the MBT/ST/MBT system.

For the MBT/ST/MBT/ST/MBT system, there should be three possible magnetic states including AFM1( ↑↑↓ ), AFM2 ( ↑↓↑ ), and FM ( ↑↑↑ ) where only magnetic moments of MBT are shown. From our calculated results, AFM1 is the most stable one, then the relative energies of FM and AFM2 states are respectively 2.4 and 3.6 meV higher than that of AFM1 state. Hence, the magnetic property and band structure of MBT/ST/MBT/ST/MBT system are corresponding to the AFM1 state. In the AFM1, the MBT/ST/MBT part of MBT/ST/MBT/ST/MBT system is still FM coupling, only loner range FM coupling do not realized (FM state). This is because that longer range FM coupling (FM state) is not only dependent on surface carriers in one ST unit mediation, but also surface carriers from the first and the second ST units of the MBT/ST/MBT/ST/MBT system. The total penetration depth of surface states of the top layer of the first ST unit and the bottom layer of the second layer is only 4-6 QLs, which is smaller than distance between the first MBT and the third MBT. Hence FM state is not the most stable one for the MBT/ST/MBT/ST/MBT system.

Band Structures and QAHE-. To confirm that such

/s40/s97/s41

FIG. 2. (a) The space distribution of spin density of MBT/ST/MBT heterostructure. The yellow and blue isosurfaces in the image respectively represent majority and minority spin.(b) The layer averaged spin density along the z direction of MBT/ST/MBT heterostructure.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00498v2-figure-f0d87d07-9e40-4ad7-866d-3fd3fbe1da13.jpg)

FIG. 3. Band structures of four different ferromagnetic heterostructures along high-symmetry lines. (a)-(d) are respectively for MBT/ST, MBT/ST/MBT, ST/MBT/ST, and MBT/ST/MBT/ST/MBT. (a1)-(d1) are the corresponding Zooming-out near Fermi level. The dashed line denotes the Fermi level.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00498v2-figure-628fed9a-63ef-43d2-85c2-1a32e059c77d.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00498v2-figure-9838e5ce-5bdc-4cb9-b03e-cb3141ceb7af.jpg)

systems can indeed realize QAHE, band structures of these heterostructures are necessary. First, we presented that bulk band gaps can be opened up around the Dirac point Γ, which in principle guarantees the insulating state [23]. Figure 3 displays the band structures along high-symmetry lines for four different heterostructures. One can see that for MBT/ST system, the opened band gap is only 3 meV [see panel (a)]; while for MBT/ST/MBT [see panel (b)] and odd-SL MBT/ST/MBT/ST/MBT [see panel (d)] systems, the same sizeable band gaps of 26 meV are opened; and for ST/MBT/ST [see panel (c)], no band gap opens but a Dirac point is formed.

Why is the Dirac point formed in the ST/MBT/ST heterostructure? Based above discussion, the surface carrier mediated FM can penetrate into 2-3 QLs of the first

FIG. 4. (a) The crystal structure and label numbers of some important atoms in the ST/MBT/ST heterostructure. (b) The characters band near the Fermi level for ST/MBT/ST heterostructure by projecting Kohn-Sham states to the local orbitals of a single atom. The contribution of each atom is charactered by the size of circles in different colors.

and second ST TIs adjent to MBT. But the distance between top layer of the first ST (5QLs) and bottom layer of the second ST(5QLs) is larger than penetration depths of the bottom layer of the first ST and top layer of the second layer. Hence, long-range FM coupling between top layer of the first ST and bottom layer of the second ST can not be realized. Based on above analysis, the ST/MBT/ST system should be still a TI one. Thus, the Dirac point is roboust for the ST/MBT/ST heterostructure. We also provide the detailed atomic contribution to the Dirac point of the ST/MBT/ST system. As shown in the Figs. 4, the bands around the Dirac point are mainly contributed by the atoms from the bottom QL of the second ST, which is in agreement with above analysis based on FM coupling of the ST/MBT/ST system. In addi-

FIG. 5. (a)-(d)Band structures, (e) total magnetic moments and magnetocrystalline anisotropy energy of MBT/ST heterostructure under zero electric field and electric field (Direction of electric-field strength E is from Sb 2 Te 3 to MnBi Te , 2 4 which is same to that of the built-in electric field of the system due to MnBi Te 2 4 with intrinsic n-type defects and Sb 2 Te 3 with intrinsic p-type defects).

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00498v2-figure-8806e072-6366-47f2-9421-82681185f5cc.jpg)

tion, we also calculated the Z2 value of the ST/MBT/ST heterostructure and obtained Z2=1. It confirms that the ST/MBT/ST is a topological nontrivial system.

In experiments, the fabricated MnBi Te 2 4 has n-type defects and Sb Te 2 3 has p-type defects, which might induce a built-in electric field in MnBi Te /Sb Te 2 4 2 3 heterostructure. We have performed some calculations to estimate the effect of built-in electric field on the band structure and magnetic property of MnBi Te /Sb Te 2 4 2 3 heterostructure. As shown in the Figs. 5(a)-(d), in comparsion with the band structure without built-in field, the band gap changes quite slightly when a small electric field is applied, but it will lead to gap closed when electric field reaches to 0.1 V/ ˚ . From the Figs. 5(e), we can A obtain that electric field almost has not effect on magnetic moment of the MBT/ST system, but it will obviously make the out-of-plane magnetic structure more stable than that of in-plane magnetic structure of MBT/ST system. Hence, low concentration of intrinsic defects in Sb Te 2 3 and MnBi Te 2 4 will have neglected influence on the band gaps, but may enhance magnetocrystalline anisotropy energy of MBT/ST heterostructures.

We now verify that the opened band gap in MBT/ST/MBT system hosts a topologically nontrivial insulator, i.e., QAHE, by calculating the Chern number, which can be obtained by integrating the Berry curvature of the occupied valence bands [53, 54]. Figure 6(a) displays the Berry curvature distribution along the high symmetry lines, where a large negative peak appears near

FIG. 6. Band structures along high-symmetry lines (a) and near the Γ point (b) of MBT/ST/MBT heterostructure. The corresponding Berry curvatures are indicated in the (a). (c) Hall conductivity as a function of energy. (d) Band structures of semi-infinite systems including the chiral edge state.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00498v2-figure-11e460d6-d067-433c-8ba4-10562d4e6094.jpg)

the Γ point and zero elsewhere. As a consequence, the total integration or the Hall conductance with the Fermi level lying inside the band gap must be nonzero. The nonzero Berry curvature means that the electrons flow in the curl field and contribute to the quantized Hall conductivity. Figure 6(b) displays the Hall conductance σ xy as a function of the energy. One can clearly see that σ xy = + e /h 2 in the energy interval near Fermi level. Hence, the nontrivial band gap for MBT/ST/MBT is 26 meV, which is beneficial for realizing high-temperature QAHE. The chiral edge states are also provided by using the surface Green function method as implemented in the WannierTools [55]. For MBT/ST/MBT system, one edge state [see Fig. 6(c)] emerges inside the band gap, connecting the valence and conduction bands, and corresponding to the Chern number of C = 1. How about the QAHE observation temperature in MBT/ST/MBT system? Based on the previously obtained Curie temperature of T C = 42 K and the gap size of 26 meV, one can conclude that the QAHE observation temperature can reach up to 42 K.

Summary-. In conclusion, our findings report that a natural magnetic van der Waals heterostructures composed of MnBi Te 2 4 and Sb Te 2 3 can be constructed to realize high-temperature QAHE. We provide proof-ofprinciple numerical demonstration that: (1) an interlayer ferromagnetic phase transition can occur at T C =42K when a 5QL Sb Te 2 3 is inserted into two SL MnBi Te 2 4 with interlayer antiferromagnetic coupling; (2) Band structures and topological property calculations verify that MnBi Te /Sb Te /MnBi Te 2 4 2 3 2 4 heterostructure system is an intrinsic QAHE with a surface gap 26 meV and the corresponding Chern number is C = 1. In addition, considering the charge compensation fact due to intrinsic n -type defects in MnBi Te 2 4 and p -type defects in Sb Te 2 3 in experiments, MnBi Te /Sb Te /MnBi Te 2 4 2 3 2 4

system should be an very appropriate candidate for exploring high-temperature QAHE.

This work was financially supported by the National Key R&amp;D Program (2017YFB0405703), NNSFC (11974098, 61434002 and 11474265), Natural Science Foundation of Hebei Province (A2019205037), Science Foundation of Hebei Normal University (2019B08), Anhui Initiative in Quantum Information Technologies. We are grateful to AMHPC and Supercomputing Center of USTC for providing the high-performance computing resources.

- ∗ Correspondence author: qiao@ustc.edu.cn
- [1] C.-Z. Chang, J. S. Zhang, X. Feng, J. Shen, Z. C. Zhang, M. Guo, K. Li, Y. Ou, P. Wei, L.-L. Wang, Z.-Q. Ji, Y. Feng, S. H. Ji, X. Chen, J. F. Jia, X. Dai, Z. Fang, S.-C. Zhang, K. He, Y. Y. Wang, L. Lu, X.-C. Ma, and Q.-K. Xue, Science 340 , 167 (2013).
- [2] J. G. Checkelsky, R. Yoshimi, A. Tsukazaki, K. S. Takahashi, Y. Kozuka, J. Falson, M. Kawasaki, and Y. Tokura, Nat. Phys. 10 , 731 (2014).
- [3] X. Kou, S.-T. Guo, Y. Fan, L. Pan, M. Lang, Y. Jiang, Q. Shao, T. Nie, K. Murata, J. Tang, Y. Wang, L. He, T.-K. Lee, W.-L. Lee, and K. L. Wang, Phys. Rev. Lett. 113 , 137201 (2014).
- [4] C. Z. Chang, W. Zhao, D. Y. Kim, H. Zhang, B. A. Assaf, D. Heiman, S.-C. Zhang,C. Liu, M. H. W. Chan, and J. S. Moodera, Nat. Mater. 14 , 473 (2015).
- [5] M. Mogi, R. Yoshimi, A. Tsukazaki, K. Yasuda, Y. Kozuka, K. S. Takahashi, and Y. Tokura, Appl. Phys. Lett. 107 , 182401 (2015).
- [6] H. Weng, R. Yu, X. Hu, X. Dai, and Z. Fang, Adv. Phys. 64 , 227 (2015).
- [7] Y. Ren, Z. Qiao, and Q. Niu, Rep. Prog. Phys. 79 , 066501 (2016).
- [8] K. He, Y. Wang, and Q.-K. Xue, Annu. Rev. Condens. Matter Phys. 9 , 329 (2018).
- [9] F. D. M. Haldane, Phys. Rev. Lett. 61, 2015 (1988).
- [10] M. Onoda and N. Nagaosa, Phys. Rev. Lett. 90, 206601 (2003).
- [11] C. X. Liu, X. L. Qi, X. Dai, Z. Fang, and S. C. Zhang, Phys. Rev. Lett. 101, 146802 (2008).
- [12] C. Wu, Phys. Rev. Lett. 101, 186807 (2008).
- [13] R. Yu, W. Zhang, H. J. Zhang, S. C. Zhang, X. Dai, and Z. Fang, Science 329 , 61-64 (2010).
- [14] Z. H. Qiao, S. A. Yang, W. X. Feng, W.-K. Tse, J. Ding, Y. G. Yao, J.Wang, and Q. Niu, Phys. Rev. B 82 , 161414(R) (2010).
- [15] Z. F. Wang, Z. Liu, and F. Liu, Phys. Rev. Lett. 110, 196801 (2013).
- [16] K. F. Garrity and D. Vanderbilt, Phys. Rev. Lett. 110, 116802 (2013).
- [17] J. Hu, Z. Zhu, and R. Wu, Nano Lett. 15, 2074 (2015).
- [18] C. Fang, M. J. Gilbert, and B. A. Bernevig, Phys. Rev. Lett. 112, 046801 (2014).
- [19] J. Wang, B. Lian, H. Zhang, Y. Xu, and S. C. Zhang, Phys. Rev. Lett. 111, 136801 (2013).
- [20] G. F. Zhang, Y. Li, and C. Wu, Phys. Rev. B 90, 075114 (2014).
- [21] H. Z. Lu, A. Zhao, and S. Q. Shen, Phys. Rev. Lett. 111, 146802 (2013).
- [22] M. Z. Hasan and C. L. Kane, Rev. Mod. Phys. 82 , 30453067 (2010).
- [23] X. L. Qi and S. C. Zhang, Rev. Mod. Phys. 83 , 10571110 (2011).
- [24] T. Jungwirth, J. Sinova, J. Masek, J. Kucera and A. H. MacDonald, Rev. Mod. Phys. 78 , 809-863 (2006).
- [25] Y. S. Hor, P. Roushan, H. Beidenkopf, J. Seo, D. Qu, J. G. Checkelsky, L. A. Wray, D. Hsieh, Y. Xia, S. Y. Xu, D. Qian, M. Z. Hasan, N. P. Ong, A. Yazdani, and R. J. Cava, Phys. Rev. B 81 , 195203 (2010).
- [26] C. Niu, Y. Dai, M. Guo, W. Wei, Y. Ma, and B. Huang, Appl. Phys. Lett. 98 , 252502 (2011).
- [27] P. P. J. Haazen, J. B. Laloe, T. J. Nummy, H. J. M. Swagten, P. Jarillo-Herrero, D. Heiman, and J. S. Moodera, Appl. Phys. Lett. 100 , 082404 (2012).
- [28] S. Qi, Z. Qiao, X. Deng, E. D. Cubuk, H. Chen, W. Zhu, E. Kaxiras, S. Zhang, X. Xu, and Z. Zhang, Phys. Rev. Lett. 117 , 056804 (2016).
- [29] Y. Ou, C. Liu, G. Y. Jiang, Y. Feng, D. Y. Zhao, W. X. Wu, X. X. Wang, W. Li, C. L. Song, L. L. Wang, W. B. Wang, W. D. Wu, Y. Y. Wang, K. He, X. C. Ma, and Q. K. Xue, Adv. Mater. 30 , 1703062 (2018).
- [30] K. He and Q.-K. Xue, Natl. Sci. Rev. 6 , 202 (2019).
- [31] X. L. Qi, T. L. Hughes, and S. C. Zhang, Phys. Rev. B 78 , 195424 (2008).
- [32] P. Wei, F. Katmis, B. A. Assaf, H. Steinberg, P. JarilloHerrero, D. Heiman, and J. S. Moodera, Phys. Rev. Lett. 110, 186807 (2013).
- [33] C. Lee, F. Katmis, P. Jarillo-Herrero, J. S. Moodera, and N. Gedik, Nat. Commun. 7, 12014 (2016).
- [34] F. Katmis, V. Lauter, F. S. Nogueira, B. A. Assaf, M. E. Jamer, P. Wei, B. Satpati, J. W. Freeland, I. Eremin, D. Heiman, P. Jarillo-Herrero, and J. S. Moodera, Nature 533, 513 (2016).
- [35] L. D. Alegria, H. Ji, N. Yao, J. J. Clarke, R. J. Cava, and J. R. Petta, Appl. Phys. Lett. 105, 053512 (2014).
- [36] M. Lang, M. Montazeri, M. C. Onbasli, X. Kou, Y. Fan, P. Upadhyaya, K. Yao, F. Liu, Y. Jiang, W. Jiang, K. L. Wong, G. Yu, J. Tang, T. Nie, L. He, R. N. Schwartz, Y. Wang, C. A. Ross, and K. L. Wang, Nano Lett. 14, 3459 (2014).
- [37] Z. Jiang, C.-Z. Chang, C. Tang, P. Wei, J. S. Moodera, and J. Shi, Nano Lett. 15, 5835 (2015).
- [38] C. Tang, C.-Z. Chang, G. Zhao, Y. Liu, Z. Jiang, C.X. Liu, M. R. McCartney, D. J. Smith, T. Chen, J. S. Moodera, and J. Shi, Sci. Adv. 3, e1700307 (2017).
- [39] D. Zhang, M. Shi, T. Zhu, D. Xing, H. Zhang, and J. Wang Phys. Rev. Lett. 122, 206401 (2019).
- [40] Y. Gong, J. Guo, J. Li, K. Zhu, M. Liao, X. Liu, Q. Zhang, L. Gu, L. Tang, X. Feng, D. Zhang, W. Li, C. Song, L. Wang, P. Yu, X. Chen, Y. Wang, H. Yao, W. Duan, Y. Xu, S.-C. Zhang, X. Ma, Q.-K. Xue, K. He, Chin. Phys. Lett. 36, 076801 (2019).
- [41] J. Li, Y. Li, S. Du, Z. Wang, B. L. Gu, S. -C. Zhang, K. He, W. Duan, Y. Xu, Sci. Adv. 5, eaaw5685 (2019).
- [42] H. Wang, D. Wang, Z. Yang, M. Shi, J. Ruan, D. Xing, J. Wang, and H. Zhang, arXiv 1907, 03380 (2019).
- [43] Y. Deng, Y. Yu, M. Z. Shi, J. Wang, X. H. Chen, and Y. Zhang, arXiv 1904, 11468 (2019).
- [44] J. Wu, F. Liu, M. Sasase, K. Ienaga, Y. Obata, R. Yukawa, K. Horiba, H. Kumigashira, S. OKUMA, T. Inoshita, and H. Hosono, arXiv 1905, 02385 (2019).

- [45] M. M. Otrokov, T. V. Menshchikova, M. G. Vergniory, I. P. Rusinov, A. Y. Vyazovskaya, Y. M. Koroteev, G. Bihlmayer, A. Ernst, P. M. Echenique, A. Arnau, and E. V. Chulkov, 2D Mater. 4, 025082 (2017).
- [46] P. E. Bl¨chl, Phys. Rev. B o 50 , 17953 (1994).
- [47] G. Kresse and J. Furthmuller, Phys. Rev. B 54 , 11169 (1996).
- [48] J. P. Perdew, J. A. Chevary, S. H. Vosko, K. A. Jackson, M. R. Pederson, D. J. Singh, and C. Fiolhais, Phys. Rev. B 46 , 6671 (1992).
- [49] V. I. Anisimov, J. Zaanen, and O. K. Anderson, Phys. Rev. B 44 , 943 (1991).
- [50] S. Grimme, J. Comput. Chem. 27 , 1787 (2006).
- [51] L. Bergqvist, O. Eriksson, J. Kudrnovsky, V. Drchal, A.
- Bergman, L. Nordstrom, and I. Turek, Phys. Rev. B 72 , 195210 (2005).
- [52] X. F. Kou, L. He, M. Lang, Y. Fan,K. Wong, Y. Jiang, T. Nie, W. Jiang, P. Upadhyaya, Z. K. Xing, Y. Wang, F. X. Xiu, R. N. Schwartz, and K. L. Wang, Nano Lett. 13, 4587 (2013).
- [53] Y. G. Yao, L. Kleinman, A. H. MacDonald, J. Sinova, T. Jungwirth, D.-S. Wang, E. Wang, and Q. Niu, Phys. Rev. Lett. 92 , 037204 (2004).
- [54] D. Xiao, M. Chang, and Q. Niu, Rev. Mod. Phys. 82 , 1959-2007 (2010).
- [55] Q. S. Wu, S. N. Zhang, H. -F. Song, M. Troyer, and A. A. Soluyanov, Comput. Phys. Commun. 224 , 405-416 (2018).
Shifei Qi Ruiling Gao Maozhi Chang Yulei Han Zhenhua Qiao
Pursuing High-Temperature Quantum Anomalous Hall Effect in MnBi$_2$Te$_4$/Sb$_2$Te$_3$ Heterostructures
Quantum anomalous Hall effect (QAHE) has been experimentally realized in
magnetically-doped topological insulators or intrinsic magnetic topological
insulator MnBi$_2$Te$_4$ by applying an external magnetic field. However,
either the low observation temperature or the unexpected external magnetic
field (tuning all MnBi$_2$Te$_4$ layers to be ferromagnetic) still hinders
further application of QAHE. Here, we theoretically demonstrate that proper
stacking of MnBi$_2$Te$_4$ and Sb$_2$Te$_3$ layers is able to produce
intrinsically ferromagnetic van der Waals heterostructures to realize the
high-temperature QAHE. We find that interlayer ferromagnetic transition can
happen at $T_{\rm C}=42~\rm K$ when a five-quintuple-layer Sb$_2$Te$_3$
topological insulator is inserted into two septuple-layer MnBi$_2$Te$_4$ with
interlayer antiferromagnetic coupling. Band structure and topological property
calculations show that MnBi$_2$Te$_4$/Sb$_2$Te$_3$/MnBi$_2$Te$_4$
heterostructure exhibits a topologically nontrivial band gap around 26 meV,
that hosts a QAHE with a Chern number of $\mathcal{C}=1$. In addition, our
proposed materials system should be considered as an ideal platform to explore
high-temperature QAHE due to the fact of natural charge-compensation,
originating from the intrinsic n-type defects in MnBi$_2$Te$_4$ and p-type
defects in Sb$_2$Te$_3$.