The magnetic Weyl semimetallic state can lead to intriguing magnetotransport, such as chiral anomaly and the layered quantum Hall effect. Mn3X (X = Sn, Ge) is a noncollinear antiferromagnetic semimetal where a Weyl semimetallic state is stabilized by time-reversal symmetry breaking. Compared to the well-studied Mn3Sn, the Weyl fermion-induced magnetotransport in Mn3Ge has been merely studied. Here, we report an in-depth study on the magnetotransport in a microfabricated Mn3Ge single crystal from room temperature to 10 K. We reveal an anomalous anisotropic magnetoresistance with fourfold symmetry and a positive high-field longitudinal magnetoresistance below the critical temperature (160–170 K). The possible origin is the temperature-induced tilting of the Weyl nodes. Our study helps to understand the magnetotransport properties in the Weyl fermion system.

Three-dimensional gapless topological excitation in Weyl and Luttinger semimetals can manifest novel properties in condensed-matter systems.1–5 In the Weyl semimetallic state, the nondegenerate energy bands linearly touch near the Fermi level and form the Weyl nodes due to the underlying symmetry breaking.6 A pair of Weyl nodes with opposite chirality act as the source and sink of the large Berry curvature, leading to nontrivial magnetotransport.7,8 The fingerprint of Weyl fermions is the chiral anomaly.9 Under the electric field applied parallel to the magnetic field, the charge pumping between two Weyl nodes with opposite chirality results in a negative longitudinal magnetoresistance (LMR). Generally, symmetry constraints involving either the inversion or time-reversal symmetry (TRS) topologically protect the Weyl fermion state.10 Therefore, the Weyl nodes cannot annihilate but tilt in the momentum space via any perturbation, respecting the translational symmetry.2 Moreover, the tilting of Weyl nodes may result in additional magnetotransport responses where the conduction electrons acquire the shifted distribution and anomalous velocity.11–14 

Among the Weyl fermion systems, the noncollinear Weyl semimetals Mn3X (X = Sn, Ge) have been intensively studied because of their giant spontaneous anomalous Hall effect (AHE) even with vanishingly small net magnetization.15–18 These compounds host a hexagonal crystal structure with an antichiral spin texture.19 The time-reversal symmetry breaking stabilizes the magnetic Weyl fermion state. The correlated electronic structure has been theoretically calculated20 and experimentally verified by angle-resolved photoemission spectroscopy.3 Despite the similar crystal structure, Mn3Ge exhibits many more pairs of Weyl nodes than Mn3Sn, connected by the Fermi arcs on the Fermi surface.20 It may result in a more complicated magnetotransport in Mn3Ge. Pioneering work on the chiral anomaly and anomalous Nernst effect has been reported for Mn3Ge.21,22 Still, the study of magnetotransport driven by Weyl fermions in Mn3Ge is far less than the well-studied Mn3Sn.3,21,23 Here, we report the magnetotransport in a Mn3Ge single crystal device microfabricated by focus ion beam (FIB) from room temperature down to 10 K. In addition to the twofold magnetotransport consistent with the chiral anomaly, it is discovered that a fourfold symmetric AMR behavior occurs below 160 K. It is in good agreement with the theoretical results, considering the tilting Weyl node.13 Meanwhile, the LMR under a high magnetic field regime (>2 T) changes the sign from negative to positive below 170 K.

The bulk Mn3Ge single crystal is grown by the bismuth flux method shown in our previous report.24 The chemical composition is determined to be 3.03:0.97 by energy dispersive x-ray (EDX) measurement. We use the FIB to microfabricate the device, with the configuration shown in Fig. 1(a) for AMR measurement. The device is about 50 µm in length, 8 µm in width, and 2 µm in thickness. The applied magnetic field is rotated in the Kagome plane with an angle φ. Before the AMR measurement, we confirmed the FIB fabricated sample quality by the temperature dependence of the longitudinal resistivity, as shown in Fig. 1(b). It shows the metallic–semimetallic transition behavior at ∼230 K, the same as the bulk. Figure 1(c) shows the AMR signals as a function of φ measured from 340 K down to 10 K under a magnetic field of 1.5 T. We further plot the AMR ratio = (ρρ)/ρ × 100% and Δρxx = ρρ as a function of the temperature in Fig. 1(d), where ‖ and ⊥ refer to the parallel and perpendicular directions of the magnetic field relative to the applied current, respectively. Mn3Ge shows a negative AMR, varying from −0.1% to −0.5%. The amplitude of the AMR ratio increases and then decreases with decreasing temperature, reaching the maximum at 100 K.

FIG. 1.

(a) Device configuration of an FIB microfabricated Mn3Ge single crystal for AMR measurement. The magnetic field is applied in the Kagome plane with an angle φ. (b) Temperature dependence of longitudinal resistivity in microfabricated Mn3Ge. (c) Angular dependence of the AMR signal from 340 to 10 K. (d) The AMR ratio and Δρxx as a function of temperature. The AMR ratio is defined as Δρxx/ρ × 100%, where Δρxx = ρρ.

FIG. 1.

(a) Device configuration of an FIB microfabricated Mn3Ge single crystal for AMR measurement. The magnetic field is applied in the Kagome plane with an angle φ. (b) Temperature dependence of longitudinal resistivity in microfabricated Mn3Ge. (c) Angular dependence of the AMR signal from 340 to 10 K. (d) The AMR ratio and Δρxx as a function of temperature. The AMR ratio is defined as Δρxx/ρ × 100%, where Δρxx = ρρ.

Close modal

In general, the AMR signal follows ρxxφ=ρ0+Δρxxcos2φ, which exhibits a twofold symmetry. However, ρxxφ deviates from the twofold function at low temperatures, as shown in Fig. 1(c). To understand the origin, we performed AMR measurements from 340 K down to 10 K under magnetic fields of 1.5, 4, 7, and 9 T. Figure 2(a) shows the angular dependence of −Δρxx under various magnetic fields at 280 K. The AMR signals have similar twofold symmetric angular dependence, but the amplitude increases when increasing the magnetic field from 1.5 to 9 T. It can be well-fitted with the cos2φ function, as shown in Fig. 2(b). We also performed the angular dependence measurement at 10 K, as shown in Fig. 2(c). Surprisingly, none of them follow the cos2φ function from 1.5 to 9 T. The AMR begins to have an additional fourfold symmetric contribution (cos4φ), besides the cos2φ component at 10 K. Therefore, we decompose the angular dependence curve under 9 T into cos2φ and cos4φ contributions, as shown in Fig. 2(d). Noticeably, the magnitude of the cos4φ component becomes as large as the cos2φ component at 10 K, a clear contrast to the 280 K case. We summarize both cos2φ and cos4φ contributions extracted from the angular variation in −Δρxx as a function of the magnetic field at 280 and 10 K cases in Fig. 2(e). It is clear that the magnitude of the cos2φ contribution is proportional to the magnetic field at 280 K but is not at 10 K. The magnitude of the cos4φ contribution is nearly zero at 280 K while it is finite and independent of the magnetic field at 10 K. It indicates that the fourfold symmetric AMR signal is caused by the temperature rather than the magnetic field. To verify the temperature-caused fourfold symmetric AMR effect, we plot the cos2φ and cos4φ contributions of −Δρxx as a function of temperature under magnetic fields of 1.5 and 9 T, as shown in Fig. 2(f). The cos4φ contribution becomes observable as the temperature decreases below 160 K under magnetic fields of 1.5 and 9 T. It is clear that the temperature induces the additional fourfold symmetric AMR. In addition, the cos2φ component changes the sign when the temperature is lower than 130 K.

FIG. 2.

(a) Angular dependence of AMR signals at 280 K under various magnetic fields. (b) Fitting of the AMR signal with cos2φ at 280 K under 9 T. (c) Angular dependence of AMR signals at 10 K under various magnetic fields. (d) Fitting of the AMR signal with cos2φ and cos4φ contributions at 10 K under 9 T. The magenta and blue curves show the cos2φ and cos4φ contributions, respectively. (e) Field dependence of AMR signals for both cos2φ and cos4φ contributions at 280 and 10 K. (f) Temperature dependence of AMR signals for both cos2φ and cos4φ contributions at 1.5 and 9 T. The error bars indicate fitting uncertainties.

FIG. 2.

(a) Angular dependence of AMR signals at 280 K under various magnetic fields. (b) Fitting of the AMR signal with cos2φ at 280 K under 9 T. (c) Angular dependence of AMR signals at 10 K under various magnetic fields. (d) Fitting of the AMR signal with cos2φ and cos4φ contributions at 10 K under 9 T. The magenta and blue curves show the cos2φ and cos4φ contributions, respectively. (e) Field dependence of AMR signals for both cos2φ and cos4φ contributions at 280 and 10 K. (f) Temperature dependence of AMR signals for both cos2φ and cos4φ contributions at 1.5 and 9 T. The error bars indicate fitting uncertainties.

Close modal

Here, we try to analyze the possible origin of the temperature-induced anomalous magnetotransport. As is shown in the supplementary material, there is no distinct difference in the magnetization curves with temperature. We thus exclude the effect of magnetism on the magnetotransport. Generally, a slight change in the magnetic structure can affect magnetotransport. Therefore, we mainly focus on the temperature-driven magnetic structure variation. While both E1g and B2g symmetries with opposite chirality are energetically equivalent in Mn3Ge,25 the neutron diffraction has demonstrated that the ground-state spin structure with E1g symmetry in Mn3Ge is stable due to the Dzyaloshinskii–Moriya interaction over the whole temperature range.19,26,27 Thus, the temperature should not change the chirality of the magnetic structure. The previous theoretical report shows that the tilting Weyl nodes under E1g symmetry lead to a higher fold symmetric AMR and planar Hall contribution, consistent with our experimental observation.13 As shown in Figs. 3(b) and 3(a), the tilting vector t = (cos θ, sin θ, 0) is assumed relative to the reference position under E1g symmetry. Meanwhile, the spins are rotated by −φ when the magnetic field B = B (cos φ, sin φ, 0) is applied as shown in Fig. 3(a). As a result, the Weyl nodes acquire a new tilting angle θ – φ, as shown in Fig. 3(b). The ultimate longitudinal conductivity is expressed as σxx=n=06āncos(nφ), where the sign and magnitude of coefficient ān are determined by chirality and θ. For the chirality of E1g symmetry, only even n components are nonzero. Thus, the fourfold symmetric component can be expected. As for the role of temperature, it may slightly change the magnetic structure, for example, by modifying the in-plane lattice constant parameter27 and may affect the tilting angle θ. Thus, the coefficient ān for the fourfold component may vanish at a high-temperature regime.

FIG. 3.

(a) Rotation of the spin structure under a magnetic field B = B (cos φ, sin φ, 0). (b) The distribution of Weyl node pairs with a chirality of χ± in the momentum space. The dashed pattern shows the original location of the Wely point with a tilting angle θ. The solid patterns are the movement of nodes due to the spin rotation, and the Weyl nodes obtain an ultimate tilting angle θφ.

FIG. 3.

(a) Rotation of the spin structure under a magnetic field B = B (cos φ, sin φ, 0). (b) The distribution of Weyl node pairs with a chirality of χ± in the momentum space. The dashed pattern shows the original location of the Wely point with a tilting angle θ. The solid patterns are the movement of nodes due to the spin rotation, and the Weyl nodes obtain an ultimate tilting angle θφ.

Close modal

The AMR effect is attributable to the chiral anomaly in the Weyl fermion system. Therefore, we expect the corresponding temperature-driven anomalous chiral anomaly behavior in Mn3Ge. A negative LMR usually evidences the chiral anomaly. The Landau levels shown in Fig. 4(a) characterize the Weyl semimetallic states. The zeroth Landau levels form the linear dispersions with opposite chirality (χ = ±1) in the Weyl fermion system. Once the electric field E is applied parallel to the magnetic field B, the electrons will be pumped from one Landau level with χ = +1 to the other with χ = −1.28 Consequently, the population difference of electrons will cause chemical potential imbalance and facilitate conduction. Meanwhile, the imbalance relaxes via the scattering between two nodes and finally reaches an equilibrium state.29, Figure 4(b) shows the LMR for the magnetic field under B // E // x and BE configurations at 280 K. Clearly, the negative LMR occurs under the B // E configuration but not under the BE configuration. This negative LMR was also reported in Mn3Sn in the whole temperature.3,23 However, when the temperature decreases down to 10 K, Mn3Ge exhibits a negative LMR under the low magnetic field regime (<2 T) while it exhibits a positive LMR at the high magnetic field regime (>2 T), as shown in Fig. 4(c). Interestingly, Chen et al.21 measured the chiral anomaly of Mn3Ge bulk single crystals under the B // E // z configuration at 0.3 K, but there is no sign change. It indicates that the temperature may only affect the in-plane magnetic structure. To study the temperature-induced anomalous LMR, we plot /dB as a function of temperature, as shown in Fig. 4(d). We can clearly see that the LMR starts to change the sign from negative to positive when the temperature is lowered below 170 K. This fact corresponds to the appearance of the fourfold symmetric AMR effect due to the tilting Weyl nodes. Noticeably, the positive LMR under a high magnetic field may also originate from the Coulomb interaction among the electrons.30 As temperature decreases, the relaxation rate between two nodes slows down. The increasing population of conduction electrons within one node enhances the Coulomb interaction under the high magnetic field regime, which may contribute to the positive LMR. In addition, a recent theory31 suggested that the fourfold symmetric AMR in ferromagnetic metals may be caused by the relaxation time anisotropy brought on by variations in the Fermi-level density of states. It may also play a similar role in antiferromagnetic Weyl semimetals. However, such a contribution should be small in a hexagonal system, such as Mn3Ge.

FIG. 4.

(a) Schematic diagram of the chiral anomaly in the Weyl fermion system. The red and blue curves represent the Landau levels. When the electric field E is applied parallel to the magnetic field B, the charge carriers with chirality χ = +1 will be pumped into the Landau levels with χ = −1, resulting in negative LMR. (b) The field dependence of LMR with B // E (red) and BE (blue) at 280 K. Δρ is defined as ρBρ(0). (c) The field dependence of LMR with B // E (red) and BE (blue) at 10 K. (d) The slopes of LMR with B // E as a function of temperature. The slope defined as dΔρ/dB is obtained from linear fitting in the high magnetic field regime (from 4 to 9 T).

FIG. 4.

(a) Schematic diagram of the chiral anomaly in the Weyl fermion system. The red and blue curves represent the Landau levels. When the electric field E is applied parallel to the magnetic field B, the charge carriers with chirality χ = +1 will be pumped into the Landau levels with χ = −1, resulting in negative LMR. (b) The field dependence of LMR with B // E (red) and BE (blue) at 280 K. Δρ is defined as ρBρ(0). (c) The field dependence of LMR with B // E (red) and BE (blue) at 10 K. (d) The slopes of LMR with B // E as a function of temperature. The slope defined as dΔρ/dB is obtained from linear fitting in the high magnetic field regime (from 4 to 9 T).

Close modal

In conclusion, we have reported temperature-induced anomalous magnetotransport in the Weyl semimetal Mn3Ge. An additional fourfold symmetric AMR and a positive LMR occur below ∼160 K. The possible origin is the temperature-induced tilting of Weyl nodes. The fourfold symmetric AMR has been experimentally observed in various material systems, but the origin is still poorly understood.31 Our study may promote exploration of the microscopic mechanism underlying fourfold symmetric AMR.

Note added in proof. While preparing the article, we became aware of a similar work, Ref. 27. They proposed another possible origin: the high-field positive LMR is driven by the metallic–semimetallic transition near 200 K. Their neutron diffraction demonstrates that there is no observable magnetic structure change with temperature but there is a small vibration of the in-plane lattice constant. We think that this lattice constant vibration may lead to Weyl node tilting with temperature. This can be the reason that the sign change in the chiral anomaly occurs under the B // E // x configuration but does not occur under the B // E // z configuration with decreasing temperature.

See the supplementary material for the magnetization characterization of Mn3Ge at room and low temperatures.

We thank S. Kurosawa for the helpful comments. This work was financially supported by JST-CREST (Grant No. JPMJCR18T3), the JST-Mirai Program (Grant No. JPMJMI20A1), the Japan Science and Technology Agency. M.W. would like to acknowledge the support from the JSPS Research Program for Young Scientists (Grant No. 21J21461). The work at the Institute for Quantum Matter, an energy frontier research center, was funded by the DOE, Office of Science, Basic Energy Sciences, under Award No. DE-SC0019331.

The authors have no conflicts to disclose.

Mingxing Wu: Data curation (equal); Formal analysis (equal); Investigation (equal); Methodology (equal); Writing – original draft (equal). Kouta Kondou: Methodology (supporting). Taishi Chen: Investigation (supporting); Methodology (supporting). Satoru Nakatsuji: Funding acquisition (equal); Writing – review & editing (supporting). Yoshichika Otani: Conceptualization (equal); Funding acquisition (equal); Project administration (equal); Supervision (equal); Writing – review & editing (equal).

The data that support the findings of this study are available from the corresponding author upon reasonable request.

1.
S.-Y.
Xu
,
I.
Belopolski
,
N.
Alidoust
,
M.
Neupane
,
G.
Bian
,
C.
Zhang
,
R.
Sankar
,
G.
Chang
,
Z.
Yuan
,
C.-C.
Lee
,
S.-M.
Huang
,
H.
Zheng
,
J.
Ma
,
D. S.
Sanchez
,
B.
Wang
,
A.
Bansil
,
F.
Chou
,
P. P.
Shibayev
,
H.
Lin
,
S.
Jia
, and
M. Z.
Hasan
,
Science
349
,
613
(
2015
).
2.
B. Q.
Lv
,
H. M.
Weng
,
B. B.
Fu
,
X. P.
Wang
,
H.
Miao
,
J.
Ma
,
P.
Richard
,
X. C.
Huang
,
L. X.
Zhao
,
G. F.
Chen
,
Z.
Fang
,
X.
Dai
,
T.
Qian
, and
H.
Ding
,
Phys. Rev. X
5
,
031013
(
2015
).
3.
K.
Kuroda
,
T.
Tomita
,
M.-T.
Suzuki
,
C.
Bareille
,
A. A.
Nugroho
,
P.
Goswami
,
M.
Ochi
,
M.
Ikhlas
,
M.
Nakayama
,
S.
Akebi
,
R.
Noguchi
,
R.
Ishii
,
N.
Inami
,
K.
Ono
,
H.
Kumigashira
,
A.
Varykhalov
,
T.
Muro
,
T.
Koretsune
,
R.
Arita
,
S.
Shin
,
T.
Kondo
, and
S.
Nakatsuji
,
Nat. Mater.
16
,
1090
(
2017
).
4.
Y.
Machida
,
S.
Nakatsuji
,
S.
Onoda
,
T.
Tayama
, and
T.
Sakakibara
,
Nature
463
,
210
(
2010
).
5.
T.
Kondo
,
M.
Nakayama
,
R.
Chen
,
J. J.
Ishikawa
,
E. G.
Moon
,
T.
Yamamoto
,
Y.
Ota
,
W.
Malaeb
,
H.
Kanai
,
Y.
Nakashima
,
Y.
Ishida
,
R.
Yoshida
,
H.
Yamamoto
,
M.
Matsunami
,
S.
Kimura
,
N.
Inami
,
K.
Ono
,
H.
Kumigashira
,
S.
Nakatsuji
,
L.
Balents
, and
S.
Shin
,
Nat. Commun.
6
,
10042
(
2015
).
6.
X.
Wan
,
A. M.
Turner
,
A.
Vishwanath
, and
S. Y.
Savrasov
,
Phys. Rev. B
83
,
205101
(
2011
).
7.
8.
S.
Jia
,
S.-Y.
Xu
, and
M. Z.
Hasan
,
Nat. Mater.
15
,
1140
(
2016
).
9.
H. B.
Nielsen
and
M.
Ninomiya
,
Phys. Lett. B
130
,
389
(
1983
).
11.
Y.
Ferreiros
,
A. A.
Zyuzin
, and
J. H.
Bardarson
,
Phys. Rev. B
96
,
115202
(
2017
).
12.
D.
Ma
,
H.
Jiang
,
H.
Liu
, and
X. C.
Xie
,
Phys. Rev. B
99
,
115121
(
2019
).
13.
G.
Chaudhary
,
A. A.
Burkov
, and
O. G.
Heinonen
,
Phys. Rev. B
105
,
085108
(
2022
).
14.
A.
Sakai
,
Y. P.
Mizuta
,
A. A.
Nugroho
,
R.
Sihombing
,
T.
Koretsune
,
M.-T.
Suzuki
,
N.
Takemori
,
R.
Ishii
,
D.
Nishio-Hamane
,
R.
Arita
,
P.
Goswami
, and
S.
Nakatsuji
,
Nat. Phys.
14
,
1119
(
2018
).
15.
S.
Nakatsuji
,
N.
Kiyohara
, and
T.
Higo
,
Nature
527
,
212
(
2015
).
16.
N.
Kiyohara
,
T.
Tomita
, and
S.
Nakatsuji
,
Phys. Rev. Appl.
5
,
064009
(
2016
).
17.
A. K.
Nayak
,
J. E.
Fischer
,
Y.
Sun
,
B.
Yan
,
J.
Karel
,
A. C.
Komarek
,
C.
Shekhar
,
N.
Kumar
,
W.
Schnelle
,
J.
Kübler
,
C.
Felser
, and
S. S. P.
Parkin
,
Sci. Adv.
2
,
e1501870
(
2016
).
18.
S.
Nakatsuji
and
R.
Arita
,
Annu. Rev. Condens. Matter Phys.
13
,
119
(
2022
).
19.
Y.
Chen
,
J.
Gaudet
,
S.
Dasgupta
,
G. G.
Marcus
,
J.
Lin
,
T.
Chen
,
T.
Tomita
,
M.
Ikhlas
,
Y.
Zhao
,
W. C.
Chen
,
M. B.
Stone
,
O.
Tchernyshyov
,
S.
Nakatsuji
, and
C.
Broholm
,
Phys. Rev. B
102
,
54403
(
2020
).
20.
H.
Yang
,
Y.
Sun
,
Y.
Zhang
,
W.-J.
Shi
,
S. S. P.
Parkin
, and
B.
Yan
,
New J. Phys.
19
,
015008
(
2017
).
21.
T.
Chen
,
T.
Tomita
,
S.
Minami
,
M.
Fu
,
T.
Koretsune
,
M.
Kitatani
,
I.
Muhammad
,
D.
Nishio-Hamane
,
R.
Ishii
,
F.
Ishii
,
R.
Arita
, and
S.
Nakatsuji
,
Nat. Commun.
12
(
1
),
572
(
2021
).
22.
C.
Wuttke
,
F.
Caglieris
,
S.
Sykora
,
F.
Scaravaggi
,
A. U. B.
Wolter
,
K.
Manna
,
V.
Süss
,
C.
Shekhar
,
C.
Felser
,
B.
Büchner
, and
C.
Hess
,
Phys. Rev. B
100
,
085111
(
2019
).
23.
S.
Kurosawa
,
T.
Tomita
,
I.
Muhammad
,
M.
Fu
,
A.
Sakai
, and
S.
Nakatsuji
, arXiv:2204.00882 (
2022
).
24.
M.
Wu
,
H.
Isshiki
,
T.
Chen
,
T.
Higo
,
S.
Nakatsuji
, and
Y.
Otani
,
Appl. Phys. Lett.
116
,
132408
(
2020
).
25.
T.
Nagamiya
,
S.
Tomiyoshi
, and
Y.
Yamaguchi
,
Solid State Commun.
42
,
385
(
1982
).
26.
J.-R.
Soh
,
F.
De Juan
,
N.
Qureshi
,
H.
Jacobsen
,
H.-Y.
Wang
,
Y.-F.
Guo
, and
A. T.
Boothroyd
,
Phys. Rev. B
101
,
140411
(
2020
).
27.
V.
Rai
,
S.
Jana
,
M.
Meven
,
R.
Dutta
,
J.
Perßon
, and
S.
Nandi
,
Phys. Rev. B
106
,
195114
(
2022
).
28.
J.
Xiong
,
S. K.
Kushwaha
,
T.
Liang
,
J. W.
Krizan
,
M.
Hirschberger
,
W.
Wang
,
R. J.
Cava
, and
N. P.
Ong
,
Science
350
,
413
(
2015
).
29.
B.
Cheng
,
T.
Schumann
,
S.
Stemmer
, and
N. P.
Armitage
,
Sci. Adv.
7
,
eabg0914
(
2021
).
30.
X.
Huang
,
L.
Zhao
,
Y.
Long
,
P.
Wang
,
D.
Chen
,
Z.
Yang
,
H.
Liang
,
M.
Xue
,
H.
Weng
,
Z.
Fang
,
X.
Dai
, and
G.
Chen
,
Phys. Rev. X
5
,
031023
(
2015
).
31.
Y.
Dai
,
Y. W.
Zhao
,
L.
Ma
,
M.
Tang
,
X. P.
Qiu
,
Y.
Liu
,
Z.
Yuan
, and
S. M.
Zhou
,
Phys. Rev. Lett.
128
,
247202
(
2022
).

Supplementary Material