Interplay between the magnetic structures and the surface states in MnBi2Te4from
first-principles studies
Zujian Dai,1, 2 Gan Jin,1, 2 and Lixin He1, 2, ∗
1CAS Key Laboratory of Quantum Information, University of Science and Technology of China, Hefei 230026, Anhui, China
2Synergetic Innovation Center of Quantum Information and Quantum Physics,
University of Science and Technology of China, Hefei, 230026, China
The antiferromagnetic (AFM) topological insulator MnBi2Te4was believed to have a topologi-
cal surface state (TSS) with large band gap due to the ferromagnetic (FM) order on the surface,
and be able to host the long-sought axion states. However, recent angle-resolved photoemission
spectroscopy (ARPES) experiments indicate that the TSS is gapless, contradicting the theoretical
predictions. Meanwhile, several experiments have suggested that there is robust out-of-plane FM
order on the surface of MnBi2Te4. To understand these seemingly contradictory results, we carry
out comprehensive first-principles calculations to investigate the interplay between the surface mag-
netism and the TSS. Our calculations provide direct evidence that in a wide range of parameters,
the (nearly) gapless TSS can coexist with the surface FM order, therefore solving the paradox of the
surface magnetism and the gapless TSS. We further show that proximity effects can be a promising
route to open the gap in the TSS of MnBi2Te4. Our research deepens the understanding of the
relationship between surface magnetism and TSS.
The interplay between magnetism and nontrivial band
topology may lead to rich physical phenomena that not
only are interesting for fundamental physics but also
have important potential applications in spintronic de-
vices [1, 2]. Therefore, MnBi2Te4, the first synthesized
intrinsic magnetic topological insulator (MTI)[3–5], has
attracted great attention since its appearance [6–10].
MnBi2Te4has an A-type antiferromagnetic (AFM) struc-
ture, which enables unique thickness-dependent topolog-
ical properties: for a thin film with an odd number of
layers, it is the QAH insulator [11–15], whereas for a thin
film with an even-number of layers, it is the long-sought
axion insulator[9, 13, 14, 16–24].
However, despite intensive research, the nature of the
topological surface states (TSSs) of this material is still
very controversial[6, 18, 25–38]. First-principles calcu-
lations predicted that MnBi2Te4has a TSS with con-
siderable energy gap (larger than 60 meV) due to the
ferromagnetic (FM) spin order on the surface[6, 8, 9, 25–
27], which is promising for achieving QAH at rather high
temperatures. Nevertheless, the zero-field QAHE was
observed only at rather low temperatures (below 1.6 K)
in this system [12, 13]. The gapped TSS was reported
in early ARPES experiments [6, 25–27]. However, more
recent ARPES measurements have observed a nearly per-
fect Dirac cone or a strong reduction of the gap at the
Dirac point on the MnBi2Te4(0001) surface [18, 28–
36, 38].
To understand the origin of the gapless TSS [37], three
scenarios have been proposed, including surface mag-
netic reconstruction [18, 31, 39], geometric reconfigu-
ration [36, 40–43] and hybridization of the surface and
bulk bands [44, 45]. For surface magnetic reconstruc-
tion, it has been shown that three types of surface spin
reorientation may led to gapless TSSs, including para-
magnetism (PM), in-plane A-type AFM, and G-type
AFM[29]. However, both time-resolved angle-resolved
photoemission spectroscopy (ARPES)[46], and magnetic
force microscopy [47] suggested that there is robust out-
of-plane FM order on the surface of MnBi2Te4. A key
perplexity here is weather the gapless TSSs can coexist
with surface FM order[46].
In this Letter, we carry out comprehensive first-
principles calculations to investigate the interplay be-
tween the surface magnetism and the TSS. Our calcu-
lations provide solid evidence that in a wide range of
parameters, the (nearly) gapless TSSs can coexist with
the surface FM order, therefore solving the paradox of
the surface magnetism and the gapless TSS. We further
show that proximity effects can be a promising route to
open the gap in the TSSs of MnBi2Te4.
The first-principle calculations are carried out with
the Atomic orbtial Based Ab-initio Computation at
UStc (ABACUS) code[48, 49]. The Perdew-Burke-
Ernzerhof[50] (PBE) exchange-correlation functional is
adopted and the DFT-D3 correction is used to account
the van der Waals (vdW) interactions[51]. A Hubbard-
like Uvalue of 4.0 eV is used for the half-filled, strongly
localized Mn 3d orbitals [52, 53]. The ABACUS code
is developed to perform large-scale density functional
theory calculations based on numerical atomic orbitals
(NAO)[48]. The optimized norm-conserving Vanderbilt
(ONCV) [54] fully relativistic pseudopotentials [55] from
the PseudoDojo library[56] are used. The valence elec-
trons for Mn, Bi and Te are 3s23p63d54s2, 5d106s26p3
and 4d105s25p4, respectively, and the NAO bases for Mn,
Bi and Te are 4s2p2d1f, 2s2p2d and 2s2p2d, respectively.
In the self-consistent and band structure calculations, the
energy cutoff for the wave functions is set to 120 Ry. Ex-
perimental lattice parameters [57] have been used. The
atomic positions are fully optimized until all forces are
less than 0.01 eV/˚
A.
arXiv:2210.14009v1 [cond-mat.mtrl-sci] 25 Oct 2022