1 Introduction
In the previous article, we have reviewed our recent works on quantum transport phenomena in topological Weyl and Dirac semimetals under magnetic fields [
1], which covers weak (anti-)localization [
2,
3], negative magnetoresistance in Weyl/Dirac semimetals [
4], and magnetotransport in the quantum limit [
5,
6]. In this review, we summarize our recent works that were not addressed in Ref. [
1], specifically, on the 3D quantum Hall effect [
7], quantum oscillations in Weyl/Dirac [
8] and nodal-line semimetals [
9], and negative magnetoresistance in topological insulators [
10], Weyl-node annihilation, [
11] and vanishing backscattering [
12] in the quantum limit. Following the structure in the previous review, these phenomena are organized depending on the strength of the magnetic field. There are several recent review papers on topological semimetals [
13-
16].
The structure of this review is as the following. In Section 2, we introduce the models we used for describing Weyl, Dirac, nodel-line semimetals, and topological insulators. In Section 3, we survey the negative magnetoresistance without chiral anomaly in topological insulators [
10]. In Section 4, we discuss quantum oscillations with the anomalous phase shift in topological semimetals [
8]. In Section 5, we present the rules for phase shifts of quantum oscillations in topological nodal-line semimetals [
9]. In Section 6, we predict a 3D quantum Hall effect of Fermi arcs in topological semimetals [
7]. In Section 7, we present the theory for the Weyl fermion annihilation in the quantum limit. The signature is recently observed as a sharp reversal in the Hall resistance of the topological Weyl semimetal TaP at extremely strong magnetic fields [
11]. In Section 8, we propose that resistance dips in the quantum limit can serve as a signature for topological insulators because of forbidden backscattering [
12]. Finally, we remark on future works in Section 9.
2 Models for topological semimetals and insulators
In this section, we introduce the models we used for describing Weyl, Dirac, nodel-line semimetals, and topological insulators.
2.1 Topological Weyl semimetal
The topological Weyl and Dirac semimetal are 3D topological states of matter [
17-
28]. Their energy bands touch at the Weyl nodes which host monopoles. The Weyl nodes have been verified in the Dirac semimetals Na
3Bi [
24,
29,
30] and Cd
3As
2 [
4,
26,
31-
36], and the Weyl semimetal TaAs family [
37-
43,
43-
52], TaIrTe
4 [
53] and YbMnBi
2 [
54]. Also, Weyl semimetals can be induced from half-Heusler compounds by applying magnetic fields [
55-
57].
According to Section 2.1 of Ref. [
1], a Weyl semimetal can be described by a two-node Hamiltonian [
58-
60]
where
are the Pauli matrices,
is the wave vector,
A,
M, and
kw are model parameters. The eigen energies of the Hamiltonian are
, with + for the conduction band and - for the valence band. The two Weyl nodes are at
, and it has been demonstrated that the model is of all the topological semimetal properties [
60], in particular, the Fermi arcs [
6], different from the
model [
2,
61,
62]. This model has the topological properties because of the
term [
63]. For the 3D quantum Hall effect in Section 6, the above model is added with two trivial
D terms [
58-
60]
2.2 Topological Dirac semimetal
A Dirac semimetal can be regarded as a Weyl semimetal and its time-reversal partner. Dirac semimetals can be studied by using the Hamiltonian [
24,
26,
64]
where
gs is the
g factor for the
s band,
gp is the
g factor for the
p band [
64],
,
, and
. The
x,
y, and
z axes in the Hamiltonian are defined along the [
100], [010], and [001] crystallographic directions, respectively.
2.3 Topological nodal-line semimetal
Nodal-line semimetals [
65-
68], in which the cross sections of conduction and valence bands are closed rings [Fig. 3(a)] in momentum space [
69-
72], have the drumhead surface states, of which the direct evidence is still missing [
73]. Recently, a new type of surface state called the floating band was discovered in the nodal-line semimetal ZrSiSe [
74]. Besides, when the symmetries are broken, the nodal-line semimetals may develop into Dirac semimetals, topological insulators, and surface Chern insulators [
75]. The nodal lines are predicted in many materials, such as HgCr
2Se
4 [
20], graphene networks [
76], Cu
3(Pd/Zn)N [
77,
78], SrIrO
3 [
79,
80], TlTaSe
2 [
81], Ca
3P
2 [
82,
83], CaTe [
84], compressed black phosphorus [
85], CaAg(P/As) [
86], CaP
3 family [
87], PdS monolayer [
88], Zintl compounds [
89], Ba
MX3 (
, Nb and Ta,
, Se) [
90], rare earth monopnictides [
91], alkaline-earth compounds [
92-
94], other carbon-based materials [
69,
95], and metallic rutile oxides
XO
2 (
, Os, Rd) [
96]. So far, the nodal lines have been verified in ZrSiS [
97-
99], PbTaSe
2 [
100,
101], InBi [
102] and PtSn
4 [
103] by ARPES.
The Hamiltonian of nodal-line semimetals can be described as [
81,
100]
where
is the wave vector,
τ,
σ are the Pauli matrices,
λ,
m, and
u are model parameters [
81,
100]. The eigen energies of the Hamiltonian are
. When
u is positive, two bands intersect at zero energy as
, which describes the nodal-line ring. The radius of the nodal ring is
[Fig. 3(a)]. When
, the dispersions result in a torus Fermi surface [Fig. 3(a)]. When
, the Fermi surface evolves into a drum-like structure [Fig. 3(a)]. Because of the low carrier density of the samples in experiments [
97,
98], we focus on the case that
. Moreover, the model in Eq. (4) is of the mirror reflection symmetry [
81,
100]. Nodal lines can also be protected by other symmetries [
67,
70,
78,
104], for example, two-fold screw rotation [
79,
80], four-fold rotation, inversion [
67] and non-symmorphic symmetry through a glide plane [
97,
98], etc.
2.4 Topological insulator
3D topological insulators can be described by the
Hamiltonian [
58,
105,
106]
where
,
,
Ci,
Mi and
Vi are model parameters. The model depicts a 3D strong topological insulator as
and
[
58]. There are four energy bands
near the
Γ point, two conduction bands and two valence bands (see Fig. 2 in Ref. [
10]). The model has been shown effective to give proper descriptions for the topological surface states [
63,
107,
108] and explain the negative magnetoresistance in topological insulators [
10,
109-
112].
In the presence of the magnetic field, the Zeeman Hamiltonian reads
where are Landé g-factors for valence/conduction bands along the z direction and in the x–y plane and is the Bohr magneton.
3 Weak field: Negative magnetoresistance in topological insulators
Recently discovered topological semimetals can host the chiral anomaly, namely, the violation of the conservation of chiral current [
113-
115], which is widely believed to be the cause of the negative magnetoresistance [
4,
35,
36,
45,
46,
116-
123]. Nevertheless, in topological insulators, where the chiral anomaly is not well defined in the momentum space, a negative magnetoresistance can also be observed. This results in great confusion [
109-
112,
124-
126] on the explanation of the negative magnetoresistance. Lately, it is found that the chiral anomaly in real space can be defined in a quantum spin Hall insulator [
127]. In Ref. [
10], we use the semiclassical Boltzmann formalism with the Berry curvature and orbital moment, to explain the negative magnetoresistance in topological insulators, and show a quantitative agreement with the experiments (see Fig. 1).
3.1 Anomalous velocity
The Berry curvature and orbital moment will induce the anomalous velocity, which may lead to the negative magnetoresistance. In experiments, the negative magnetoresistance exists above
K [
111], thus quantum interference mechanisms can be excluded. In addition, due to the poor mobility in the topological insulators Bi
2Te
3 and Bi
2Se
3 [
129], when the magnetic field is up to 6 Tesla, the Landau levels cannot be well-formed. In the semiclassical regime, the electronic transport can be described by the equations of motion [
130]
where both the position
r and wave vector
k appear simultaneously,
and
are their time derivatives,
is the electron charge,
E and
B are external electric and magnetic fields, respectively.
,
is the band dispersion,
m is the orbital moment induced by the semiclassical self-rotation of the Bloch wave packet, and
Ωk is the Berry curvature [
131]. In the linear-response limit (
E=0), Eq. (7) yields an effective velocity
where is the correction to the density of states, and
Due to the Berry curvature, the velocity develops an anomalous term, which is proportional to
B. Note that the conductivity is the current–current (velocity–velocity) correlation [
132], thus the presence of the anomalous velocity is expected to generate an extra conductivity, which grows with the magnetic field, namely, a negative magnetoresistance. It has been implied that the negative magnetoresistance in topological semimetals is related to the Berry curvature [
1,
133,
134]. In our previous review [
1], we have shown that the Berry curvature [
130] can lead to a conductivity correction that grows with magnetic field
B,
where
kF is the Fermi wave vector and the carrier density
n is directly proportional to
. This is an alternative understanding to the negative magnetoresistance induced by the chiral anomaly. Ref. [
10] shows that this mechanism is large enough in topological insulators as those observed in the experiments, where the relative magnetoresistance can exceed
in a parallel magnetic field of several Tesla [
109-
112,
124-
126].
3.2 Magnetoresistance formula
In our calculation, the relative magnetoresistance is defined as
. In the semiclassical Boltzmann formalism, the longitudinal conductivity
is contributed by all the bands crossing the Fermi energy, and for band
n [
134]
where
n is suppressed for simplicity,
Dk and
are given by Eq. (9),
is the Fermi distribution in equilibrium, the transport time
τ is assumed to be a constant in the semiclassical limit [
135]. For the
n-th band of the Hamiltonian
H, the
ξ component of the Berry curvature vector can be found as
, where
stand for
, and
is the Levi-Civita anti-symmetric tensor, and
where . The orbital moment m can be found as
The Zeeman energy can induce a finite distribution of
Ω and
m [
10].
3.3 Comparison with negative magnetoresistance in experiments
Figure 1 shows that the numerically calculated relative magnetoresistance in parallel magnetic fields are negative and decrease monotonically with the magnetic field. They can be fitted by
at small magnetic fields and conform to the Onsager reciprocity MR
=MR
. We also use a tight-binding model [
136] to justify the calculation. Figure 1 shows a good agreement on the negative magnetoresistance between the experiments and our numerical calculations. The current direction and temperature are from the experiments and the model parameters are from the
calculations [
106] and experiments [
128]. In the experiment by Wang
et al.[
109], the temperature is 1.8 K, so the original data (orange triangles) has a positive magnetoresistance near zero field due to the weak anti-localization [
137-
140]. The weak anti-localization induces a positive magnetoresistance [
2], which is subtracted before the comparison. In the experiments by Wiedmann
et al. [
111] and He
et al.[
110], the temperatures are 29 K and 300 K, there is no weak anti-localization effect. The negative magnetoresistance does not change much with temperature [
10], consistent with the experiments and showing the semi-classical nature of the negative magnetoresistance. The negative magnetoresistance becomes enhanced as the Fermi level approaches the band bottom [
10], indicating the role of the Berry curvature. Ref. [
10] also shows that the signs of
g-factors in the Zeeman coupling determine the signs of magnetoresistance qualitatively. In the experiment, the techniques used for topological insulators, for example, electron spin resonance and quantum oscillations, cannot determine the signs of
g-factors but only their absolute values [
128,
141]. Transport measurements can determine the sign of the
g-factor only in specific setups [
142]. The orbital moment has been neglected in most of the literature studying the magneto-transport using the semiclassical formalism [
133,
134]. As shown in Refs. [
143,
144], the orbital moment is essential for the magnetoresistance anisotropy in a Weyl semimetal. Moreover, the correction
m·
B to
ε(
k) can enhance the band separation and the negative magnetoresistance. The orbital moment effectively enhances the MR
x a few times larger. MR
z can be even positive without
m.
3.4 Discussions
The conventional equations of motion in the low-field semiclassical regime (
) are only accurate to the linear order in the external fields (
E and
B), thus for studying magnetoconductivity which is an intrinsically nonlinear coefficient, the obtained results would not be complete. Based on a recently developed semiclassical theory with second-order accuracy [
145,
146], a complete theory of magnetoconductivity for general 3D nonmagnetic metals was formulated within the Boltzmann framework with the relaxation time approximation [
144]. The work shows several surprising results. First, there is an important previously unknown Fermi surface contribution
to the magnetoconductivity, termed as the intrinsic magnetoconductivity, because the ratio
is independent of the relaxation time. Here
is the conductivity at
. Second, a pronounced
term can lead to the violation of Kohler’s rule. Previously, any deviation from Kohler’s rule is usually interpreted as from factors beyond the semiclassical description or from the presence of multiple types of carriers or multiple scattering times. The result here reveals a new mechanism for the breakdown of Kohler’s rule. Third,
can lead to a positive longitudinal magnetoconductivity (or negative longitudinal magneto-resistivity). The effect is independent of chiral anomaly effect for the Weyl/Dirac fermions, and it can occur for a generic doped semiconductor without any band crossings. This indicates that positive longitudinal magnetoconductivity measured in the semiclassical regime alone cannot be regarded as smoking-gun evidence for the existence of topological band crossings. The intrinsic contribution generally exists in 3D metals with nontrivial Berry curvatures, and should be taken into account when interpreting experimental results. It may already play an important role behind the puzzling magnetotransport signals observed in recent experiments on TaAs
2 and related materials.
In the quantum limit where only the lowest Landau band is occupied, magnetoresistance depends subtly on scattering mechanisms [
5,
6,
147], rather than the Berry curvature and orbital moment. The current-jetting effect is usually induced by inhomogeneous currents when attaching point contact electrodes to a large bulk crystal and may also hamper the interpretation of the negative magnetoresistance data [
148]. A recent work also has pointed out that the negative magnetoresistance may exist without the chiral anomaly [
149]. In (Bi
1-xIn
x)
2Se
3, it is proposed that the in-plane negative magnetoresistance is due to the topological phase transition enhanced intersurface coupling near the topological critical point [
126]. In addition, it is also found that the magnetoresistance is robust against the deviation from the ideal Weyl Hamiltonian, such as the shifted Fermi energy, nonlinear dispersions, and the Weyl node annihilation [
150].
4 Strong field: Quantum oscillation in Weyl and Dirac semimetals
When applying a magnetic field in the
z-direction, the energy spectrum evolves into a series of 1D Landau bands [
6,
60] (see Figs. 1(b) and (c) of [
8]), which result in the Shubnikov-de Haas (SdH) oscillation of resistance. The oscillation of the resistivity
ρ can be demonstrated by the Lifshitz–Kosevich formula [
151]
where φ is the phase shift, F is the oscillation frequency and B is the magnitude of magnetic field. F and φ can reveal valuable details on the Fermi surface of the material.
The phase shift of each frequency component can be argued as the following:
where
is a correction, which emerges only in 3D, and
is the Berry phase [
131,
152]. The curvature of the Fermi surface along the direction of the magnetic field determines the sign of
[
153-
156]. When the cross section is maximum, for electron carriers
and for hole carriers
; when the cross section is minimum, for electron carriers
and for hole carriers
. For the sphere Fermi surface, there is only a maximum, thus for electron carriers
and for hole carriers
. For a parabolic energy band, it does not have Berry phase, thus the phase shift is
and
in 2D and 3D, respectively. However, a linear energy band (for example, Weyl and Dirac semimetals [
17-
20,
157]) has an extra
π Berry phase [
152], thus the phase shift is 0 [
158] and
[
159] in 2D and 3D, respectively. In a nodal-line semimetal [
65-
68], an electron can collect a nontrivial
π Berry phase around the loop encircling the nodal line [
79]. The phase shifts of 2D and 3D bands with linear and parabolic dispersions are summarized in Table 1.
Topological Weyl/Dirac semimetals and nodal-line semimetals provide a novel platform to explore the nontrivial Berry phase [
3,
119,
122,
159-
175]. In Ref. [
8], we show that, near the Lifshitz point, the phase shift of the quantum oscillation can go beyond recognized values of
or
and nonmonotonically move toward a wide range between
and
. However, these values in experiments may be misunderstood as
. For Dirac semimetals, the total phase shift adopts the discrete values of
or
. In recent experiments of electron carriers, the positive phase shifts are observed and may be explained by our findings. In addition, a new beating pattern, resulting from the topological band inversion, is found. Up to now, quantum oscillations have been inspected experimentally in HfSiS [
176], ZrSiS [
177-
182], ZrSi(Se/Te) [
183] and ZrGe(S/Se/Te) [
184], but the phase shifts have been concluded different. In Ref. [
9], the phase shifts and frequencies (Table 5) of nodal-line semimetals are extracted from analytic results of the calculated resistivity. We also summarize the generic rules for phase shifts in random cases (Table 4). The generic rules assist us to handle several materials, for example, ZrSiS and Cu
3PdN [
177-
182].
4.1 Quantum oscillation in linear and parabolic limits of a Weyl semimetal
The resistivity is calculated in two direction configurations according to linear response theory [
185-
188]. For the longitudinal configuration, the resistivity
is examined along
z direction. For the transverse configuration, the resistivity
is examined along
x direction. The magnetoresistivity in the linear dispersion limit and parabolic dispersion limit takes the form of Eq. (14). In Table 2, we list the analytic expressions for the phase shift
φ and frequency
F in these two limits.
4.2 Resistivity peaks and integer Landau indices
In the experiment, the peak positions or valley positions, that are on the
B axis, are given the integer Landau indices
n, then the phase shift
φ and frequency
F can be extracted from a plot of
n and
[see inset of Fig. 1(d) of [
8]]. Nevertheless, it is still in debate that whether the peaks [
122,
159,
160,
162-
165] or valleys [
167,
168,
171,
189] should be given the Landau indices. Our results explicitly reveal that the resistivity peaks of
and
emerge near the Landau band edges and are in correspondence with the integer Landau indices. We evaluate the resistivity components theoretically from the conductivity components [
190,
191]. For the longitudinal configuration, the resistivity
, where the conductivity
is along
z direction. Near band edges, the conductivity
shows valleys due to vanishing velocities, thus
shows peaks. For the transverse configuration,
, and the longitudinal Hall conductivity and field-induced Hall conductivity are as follows:
where
is the oscillation part and
represents the zero-field conductivity. In
, the
δ term, which comes from the disorder scattering, was rarely considered before. A consequence of the
δ term is that
, up to the first order of
δ. As
and
are both proportional to
, their peaks are lined up for the random ratio of
to
(but when
, the oscillation is so weak that it is hard to be observed). This new finding comes from the disorder scattering term
δ in the Hall conductance. At the same time, the
valleys are lined up with the
peaks, since
, which comes from diffusion, is in proportion to the scattering times. However,
, which arises from hopping is inversely in proportion to the scattering times [
60,
186,
192]. To sum up, the peak positions follow the relation
, thus
and
present peaks near Landau band edges and their phase shifts are the same.
4.3 Anomalous phase shift near the Lifshitz point of Weyl and Dirac semimetals
Figure 2 shows the results of the frequency
F and the phase shift
φ calculated numerically for the model in Eq. (1). In Fig. 2(c), the numerical results are in agreement with the analytical prediction
for the linear limit (
) and
for the parabolic limit (
). We define
and
. For
, the
φ–
EF curves break as the beating patterns emerge. In Fig. 2(c), for
, the phase shift drops below
rather than shift monotonically from
to
around the Lifshitz transition point (namely,
). This is because there is no simple
dependence [
8]. At the Lifshitz point, we can analytically show that the phase shift is
, which is in agreement with that in Fig. 2(c). It is equivalent to
, which is usually considered to originate from the
π Berry phase as the Fermi sphere encircles the single Weyl node. Nevertheless, the Fermi sphere encircles two Weyl nodes when the Fermi energy is at the Lifshitz point. For
, there is no nonmonotonicity in
.
A Weyl semimetal combine with its time-reversal counterpart can compose a Dirac semimetal. The model of a Dirac semimetal can be made by
H(
k) in Eq. (1), integrated with its time-reversal counterpart
H*(-
k), where the asterisk indicates complex conjugate. This model can be treated as a building block for Weyl semimetals, which respect time-reversal symmetry and meanwhile break inversion symmetry [
37-
39,
42-
48]. Here, for the Dirac semimetal, the total phase shift may take two values,
for
and
or
for
. Around the Lifshitz point, the total phase shift may change between the two values.
The electron carriers is supposed to yield negative phase shifts and the hole carriers is supposed to yield positive phase shifts [
159]. Nevertheless, in experiments of the Dirac semimetal Cd
3As
2, the phase shift for electron carriers takes positive values [
160,
162,
167]. One explanation may be that, for the phase shift
to
in the experiments, their actual values are around
to
due to the
periodicity. Our numerical results show that the total phase shift adopts these values from around the Lifshitz point to higher Fermi energies, which is also consistent with the carrier density in the experiments. In Table 3, we propose the counterparts for the experimental values of the phase shift. In Ref. [
8], we also demonstrate that beating patterns will emerge because of the band inversion, that is different from orbital quantum interference [
193], the Zeeman splitting [
165,
166,
194] and nested Fermi surfaces [
162].
5 Strong field: Quantum oscillation in nodal-line semimetals
5.1 Phase shifts of nodal-line semimetals
The frequencies
F and phase shifts
φ can be analyzed from the Fermi surface described by Eq. (4). From the Onsager relation, we have
, where
A is the area of the extremal cross section on the Fermi surface perpendicular to the magnetic field [
195]. When the nodal-line plane is perpendicular to a magnetic field, specifically
Bz here, there are two extremal cross sections at the
plane [Fig. 3(b)]. Combining this dispersion with the Onsager relation, we find that the high frequency is
for the outside circle and the low frequency is
for the inside circle. These two frequencies may lead to a beating pattern [
196].
Phase shifts are more complicated in the nodal-line semimetals. First of all, whether the magnetic fields are in-plane or out-of-plane, the torus Fermi surface has both maximum cross section and minimum cross section. Secondly, the Berry phase is 0 along the circle parallel to the nodal line and is π along the circle enclosing the nodal line. Therefore, depending on direction of the magnetic field, the quantum oscillation allows different phase shifts, which is summarized in Table 1. When the nodal-line plane is perpendicular to a magnetic field, there are two cross sections α and β, shown in Fig. 3. Along loops of the cross sections α and β, the Berry phase is 0, thus phase shifts take values or . For the electron carriers, (see Table 4), the phase shifts of the maximum cross sectionα are and . The phase shifts of the minimum cross section β are and , that is equivalent to since the oscillation has periodicity.
5.2 Magnetoresistivity of nodal-line semimetals
To testify the above conclusion, we calculate the resistivities along the z (, out of the nodel-line plane) and x (, in the nodel-line plane) directions, respectively, according to and . The calculations show that for both and , there are two terms in the magnetoresistivities,
where is the phase shifts and are the oscillation frequencies. Their analytic expressions are listed in Table 5. We show that the resistivity calculations and the Fermi surface analysis are equivalent for the phase shifts of the quantum oscillation in the nodal-line semimetal. The phase shifts, in a magnetic field parallel to the nodal-line plane [Fig. 3(c)], can also be found in a similar way, as listed in Table 4.
5.3 Discussions
For nodal-line semimetals, most of the quantum oscillation experiments have been done for the ZrSiS family materials [
176-
184], in which there are both electron and hole pockets at the Fermi energy [
97,
98,
177,
182]. ZrSiS (see Fig. 3 in Ref. [
9]) has the diamond-shaped electron pockets encircling the nodal line and the quasi-2D tubular-shaped hole pockets at the
X points [
182]. When the magnetic field is normal to the diamond-shaped Fermi surface, there are three extremal cross sections, the outer (
α) and inner (
β) cross sections of the diamond-shaped electron pocket, and the tubular-shaped hole pockets (
γ). The
α and
β cross sections of the diamond-shaped Fermi surface take phase shifts of
and
, respectively. However, the frequencies of the
α and
β pockets are so large (about
T) that it is hard to be extracted from the experiments. In contrast, the
γ nodal-line hole pockets have been identified as the origin of the high frequency (about 210 T) component [
182]. This pocket encircles a nodal line, thus it has a
π Berry phase. In addition, it has
due to its quasi-2D nature. From Eq. (15), we find that the
γ pocket has a phase shift of
, which is in agreement with the results obtained by Ali
et al. [
178]. In Ref. [
9], we also analyzes the phase shifts for Cu
3 PdN [
77,
78].
When the symmetry protecting the nodal line is broken, there is a finite gap Δ that separates the conduction and valence bands. The Berry phase becomes
The nodal-line semimetals may also be distinguished from their weak localization behaviors that cross between 2D weak anti-localization and 3D weak localization [
197].
6 Strong field: 3D quantum Hall effect
The discovery of the quantum Hall effect in 2D opens the door to the field of topological phases of matter [
198,
199]. In 3D electron gases, the extra dimension along the magnetic field direction prevents the quantization of the Hall conductance. Thus, the quantum Hall effect is normally observed in 2D systems [
158,
198,
200-
202]. In Ref. [
7], we show a 3D quantum Hall effect in a topological semimetal. The topological semimetal can be considered as a 2D topological insulator for momenta (
kz here) between the Weyl nodes, resulting in the topologically protected surface states [in Fig. 4(c)] at the surfaces parallel to the Weyl node separation direction. The topologically protected states form the Fermi arcs [
24,
26,
26,
29-
34,
37-
39,
42,
44,
203,
204] on the Fermi surface [red curves in Figs. 4 (a) and (b)]. The transport signature of the Fermi arcs is an intriguing topic [
205-
209].
6.1 Wormhole tunneling via Fermi arcs and Weyl nodes
The topological nature requires that only a region between the Weyl nodes can be occupied by the states of Fermi arcs [
6] [Fig. 4(b)]. At one surface, a closed Fermi loop, which is essential to the quantum Hall effect, cannot be formed by the Fermi arcs. However, in a topological semimetal slab, the Fermi arcs from opposite surfaces [Fig. 4(c)] can form the required closed Fermi loop [Fig. 4(d)]. Thus electrons can tunnel between the Fermi arcs at opposite surfaces via the Weyl nodes [Figs. 4(e)–(g)]. The Fermi loop formed by the Fermi arcs at opposite surfaces via the Weyl nodes can support a 3D quantum Hall effect. To be specific, the Weyl nodes act like “wormholes” that connect the top and bottom surfaces, and an electron can complete the cyclotron motion. Since the Weyl nodes are singularities in momentum, the wormhole tunneling can be infinite in real space, according to the uncertainty principle. The time scale of the “wormhole” tunneling is 0 as the Weyl node on the top surface and the Weyl node on the bottom surface are the same one and are described by the one coherence wavefunction. In experimental materials, the tunneling distance is limited by the mean free path, which can be comparable to or longer than 100 nm in high-mobility topological semimetals [
43,
45,
46,
119,
160,
162,
167,
210,
211], even up to 1 μm [
212], thus the thickness in the calculation is chosen to be 100 nm. The wormhole effect has been addressed in different situations in topological insulators [
213]. The quantum Hall effect solely from the Fermi arcs requires the bulk carriers to be depleted by tuning the Fermi energy to the Weyl nodes [
214]. Compared to the novel quantum oscillations [
212,
215], the quantum Hall effect of the Fermi arcs contributes a quantized complement to the Fermi arc dominant electronic transports. The Weyl semimetals TaAs family [
37-
39,
42-
46,
49-
52,
210] and the Dirac semimetals Cd
3As
2 and Na
3Bi have extremely high mobilities [
119,
160,
162,
167,
211] required by the quantum Hall effect. Low carrier densities [
4,
35,
36] and gating [
36] have also been achieved. The 3D quantum Hall effect of the Fermi arcs is expected in slabs of the TaAs family [
37-
39,
42-
46,
210,
214], [110] or
Cd
3As
2 [
35,
216-
219], and [100] or [010] Na
3Bi.
6.2 Quantized Hall conductance
We can calculate the Hall conductivity from the Kubo formula [
7,
199,
220-
224]. Figure 5(b) presents the sheet Hall conductivity
for the topological semimetal slab. When the Fermi energy is far away from the Weyl nodes, the sheet Hall conductivity obeys the usual
dependence. The closer the Fermi energy moves towards the Weyl nodes, the smaller the slope becomes, which indicates that the carrier density is decreasing. Moreover, when the Fermi energy moves towards the Weyl nodes, the quantized plateaus of
begin to arise. Note that a 100-nm slab is still a 3D object.
6.3 3D distribution of the edge states
Figures 4(h) and (i) show that the edge states of the Fermi arcs have a unique 3D spatial distribution. Specifically, the top edge states propagate to the left (green arrow) and the bottom edge states to the right (orange arrow). This unique 3D distribution of the edge states of the Fermi arcs can be probed by scanning tunneling microscopy [
225] or microwave impedance microscopy [
226]. Different from topological insulators [
201,
202], the Fermi-arc quantum Hall effect requires the collaboration of the two surfaces.
6.4 Topological Dirac semimetals
A single surface of the Dirac semimetal can form a complete Fermi loop needed by the quantum Hall effect due to the time-reversal symmetry. However, the single surface Fermi arc loop is not that robust and may be distorted [
227]. Therefore, it may present different characteristics compared to the two-surface Fermi arc loop. For the [
112] and [
110] Cd
3As
2 and [010] Na
3Bi [
30] slabs, the parameters from Ref. [
228] yield the 3D quantum Hall effect, which may exhibit a fourfold degeneracy.
7 Extremely strong field: Weyl fermion annihilation
The topological properties of the Weyl nodes can be revealed by studying the high-field transport properties of a Weyl semimetal. The lowest Landau bands of the Weyl cones remain at zero energy unless a strong magnetic field couples the Weyl fermions of opposite chirality. The coupled Weyl fermions lost their chiralities and acquire masses, two of the most characteristic features of the Weyl fermion. In this sense, the Weyl fermions are annihilated. In the Weyl semimetal TaP, we achieve such a magnetic coupling [
11]. Their lowest Landau bands move above chemical potential, leading to a sharp sign reversal in the Hall resistivity at a specific magnetic field corresponding to the W
1 Weyl node separation. In the following, we use a model calculation to show the physics.
7.1 Calculation of Landau bands with magnetic fields in the x–z plane
In Section 2.1 of Ref. [
1] and Eq. (1), we have given a minimal model for a Weyl semimetal
where
σ are the Pauli matrices,
,
is the wave vector, and
A,
are model parameters. When
, the intersections of the two bands are at
where
(see Fig. 1 of [
1]), leading to the topological semimetal phase.
In Section 2.6 of Ref. [
1], we have given the Landau bands in a magnetic field along the
z direction. Now we generalize the case to an arbitrary field applied normal to the
y direction
, where
φ is the angle between the
z and field directions. The Landau gauge can be chosen as
. Under the Pierls replacement
the Hamiltonian becomes
where and . Define the guiding center
and the ladder operators [
229]
the Hamiltonian becomes
where
We have defined , which is the summation of the projections of kx and kz along the direction of the magnetic field and can serve as a good quantum number.
7.2 Landau bands in the z-direction magnetic field
At , i.e., the magnetic field is applied along the z direction, the Hamiltonian reduces to
where , , . With the trial wave functions for (later denoted as ) and for , where ν indexes the Hermite polynomials, the eigen energies E can be found from the secular equation
for , and for , where . The eigen energies are found as
They represent a set of Landau energy bands (ν as band index) dispersing with kz. The eigen states for are
and for is
where , and the wave functions are found as
where , is the area of the sample, the guiding center , are the Hermite polynomials. As the dispersions are not explicit functions of kx, the number of different kx represents the Landau degeneracy in a unit area in the x–y plane.
7.3 Landau bands in the x-direction magnetic field
If , the Hamiltonian can be solved numerically. We can write Eq. (24) as
where , and . can be readily diagonalized to give the band spectrum
where . The wave functions for the conduction and valence bands are
with , . Denote
we can calculate the off-diagonal matrix elements in the basis of as
Then the energy spectrum along arbitrary directions in the x–z plane can be solved numerically.
Figure 6 shows the Landau bands of the Weyl semimetal. When the magnetic field is along the
z direction (
), the lowest Landau band (red) crosses the Fermi energy (dashed line). As the magnetic field is rotated from the
z direction to the
x direction (
), the lowest Landau band is shifted and evolving. When the magnetic field is along the
x direction, the spectrum of the Landau bands is particle-hole symmetric and there is a gap due to the coupling between the Weyl fermions at the opposite Weyl nodes. This gap is why there is a sharp sign reversal in the Hall resistivity in the strong-field quantum limit of the Weyl semimetal TaP [
11]. Because of the gap, the Weyl fermions acquire masses and lose their chiralities. Since having chirality and no mass are two features of the Weyl fermion, the Hall signal therefore indicates that the Weyl fermions are annihilated.
8 Extremely strong field: Forbidden backscattering and resistance dip in the quantum limit
The discovery of 3D topological insulators [
58,
230,
231], whose characteristics is of topologically protected 2D surface states, shines light on the exploring of exotic topological phases [
232-
237]. Therefore, distinguishing the bulk-state transport to identify topological insulators is an intriguring topic.
In strong magnetic fields, 1D Landau bands are formed from the quantization of the bulk states of a 3D topological insulator. In 2D, lowest Landau levels cross each other, which servers as a signature for the quantum spin Hall phase [
238,
239]. In 2D, other approaches can be deployed to probe the quantum spin Hall phase, for example, interference effects [
240]. However, in 3D, it is seldom addressed that whether the lowest Landau band could be used to identify a topological insulator. In Ref. [
12], we study the resistance of a 3D topological insulator in the strong-field quantum limit, namely, only the lowest Landau band is occupied [Fig. 7(b)]. We find that the backscattering can be totally suppressed in the quantum limit at a critical magnetic field, which can be used to identify the topological insulator phases. Besides, this forbidden backscattering is absent in topological semimetals [
5,
6,
147]. This theory is consistent with the recent experiment [Figs. 7(c) and (d)]. Moreover, this mechanism will be practical for those materials with small gap, for example, the ZrTe
5 [
241,
242] families and the Ag
2Te [
243].
8.1 Forbidden backscattering in the quantum limit
In a strong magnetic field B along the z direction, the energy spectrum quantizes into a series of 1D Landau bands [Figs. 7(a) and (b)]. The energies of the lowest two Landau bands, denoted as and , are , where the magnetic length , the electron charge , and the mass term
We can determine the gap between the two lowest Landau bands by m with .
Next, we will concentrate on an electron-doped quantum limit, namely, the Fermi energy intersects only with the Landau band, of which the eigenstate is
where we have defined
and
is the state of a usual zeroth Landau level multiplying a plane wave function along the
z direction [
5].
In solids, the electronic transport is relatively influenced by the backscattering, which plays a dictating role in the presence of the 1D Landau band, since there are only two states at the Fermi energy, as denoted by kF and in Fig. 7. The backscattering between these two states is characterized by the scattering matrix element between them. From the spinor eigenstate in Eq. (37), we find that the modular square of the scattering matrix element between the kF and states is in proportion to the form factor
IS vanishes when
, that is, the backscattering between state
kF and state
is forbidden. According to Eq. (36),
m vanishes at a critical magnetic field
Bc evaluated by
, where
kz is equal to the Fermi wave vector
kF at the Fermi energy. For a topological insulator,
and
, thus
Bc holds finite solutions, at which the backscattering is totally suppressed. A rich phase diagram can be found in Fig. 8. This forbidden backscattering will result in a dip in the resistance as a function of the magnetic field, which can be probed in experiments and can serve as a signature for topological insulator phases. This forbidden backscattering is an eigenstate property and thus is new and different from the mechanism of Landau level crossing [
238,
239], which is a spectrum property.
8.2 Conductivity in the quantum limit
Along the direction of the magnetic field, there is no Hall effect, thus the resistivity is the inverse of the conductivity, namely,
. In the quantum limit, only band
contributes to the conductivity. Following the methods in Refs. [
5,
6,
244], the resitivitity can be expressed as
where
is the conductivity independent of the spinor inner product part and
IS is the form factor. For different types of scattering potential,
takes different forms. However, the form factor
IS in Eq. (40) indicates that, for a topological insulator, the resistance always has a dip, regardless of
. Figure 4 in Ref. [
12] reveals the resistivity when the Gaussian and screened Coulomb potentials are present. They both show clear dips in the resistivity for some weak topological phases (
and
in row 1, columns 1 and 2) and strong topological insulator phases (
and
in row 1, column 3). Furthermore, the positions of the minima on the
B axis do not change for various potentials. The spin-orbit scattering can improve the above picture and result in a better fitting to the experiment. Specifically, in Fig. 7(d), the spin-orbit coupling is included in the screened Coulomb scattering potential. By choosing proper fitting parameters, we obtain a -40% change at the dip of the resistance, which is in agreement with the experiment, as shown in Figs. 7(c) and (d).
9 Remarks and perspective
The theories of the quantum oscillations need to be improved to match the semiclassical argument and full quantum mechanics calculations. Recently, a quantum theory of intrinsic magnetoresistance for three-dimensional Dirac fermions in a uniform magnetic field is proposed, which shows that the relative magneto-resistance is inversely quartic of the Fermi wave vector and only determined by carrier density, and a formula for the phase shift in SdH oscillation is present as a function of the mobility and the magnetic field [
245]. Furthermore, new discoveries on quantization rules in oscillations have been found for graphene, 2D materials, topological metals, topological crystalline insulator, and Dirac and Weyl semimetals [
246,
247]. Topological contributions have also been found in Bloch oscillations [
248]. Generalizations of these classical notions to nodal-line systems will be topics of fundamental interests in the future. Whether the rules for phase shift can be generalized to extremal orbits shared by electron and hole pockets in type-II Weyl semimetals [
246,
249] will be an outstanding problem. Quantum oscillations in type-II Dirac semimetals PdTe
2 [
250] and nodal-line systems [
251,
252] have also been addressed.
Lately, the quantized Hall resistance plateaus have been experimentally observed in the topological semimetal Cd
3As
2 [
217,
218,
253], with thickness ranging from 10 to 80 nm. They cannot be regarded as 2D. Nevertheless, several questions still hold. First, Cd
3As
2 is a Dirac semimetal, composed of two time-reversed Weyl semimetals. At a single surface, there is a complete 2D electron gas, formed by two time-reversed half 2D electron gases of the Fermi-arc surface states. There may be also the trivial quantum Hall effect on a single surface. Second, the 3D bulk states quantize into 2D subbands for those thicknesses. If the 3D bulk states cannot be depleted entirely, they also have the trivial quantum Hall effect. The two issues may explain the 2-fold and 4-fold degenerate Hall resistance plateaus observed in the experiments. To deplete the 3D bulk states, the Fermi energy has to be placed exactly at the Weyl nodes. How to distinguish these trivial mechanisms from the 3D quantum Hall effect will be an interesting direction. Previously, when studying the geometric phase, the parameter space is usually either in real space or momentum space [
131]. The Weyl orbit formed by the Fermi arcs and Weyl nodes is a new physics, because part of the geometric phase is accumulated in real space and part in momentum space, quite different from the parameter spaces studied before. In particular, the geometric phase has a thickness dependence when accumulated along the path as electrons tunnel between the top and bottom surfaces [
215,
254]. Recently, a new experiment uses this thickness-dependent phase shift to demonstrate the contribution of the Weyl orbit in the observed quantized Hall resistance [
219]. The 3D quantum Hall effect can also be supported by the CDW mechanism [
255], which has been observed recently in ZrTe
5 [
256]. More works will be inspired to verify the mechanism and realize the 3D quantum Hall effect in the future [
257].
In the quantum limit, our theory have shown that up to two resistance dips may appear if the system is a 3D strong topological insulator [
12]. Surprisingly, recent experiments report up to five oscillations in the quantum limit [
258,
259]. The oscillation as a function of the magnetic field follows a logithimic scale invariance law, much like those in the Efimov bound states of cold atoms. The Efimov bound state is a three-body bound state arising from the two-body interactions between atoms. Efimov-like bound states have been used to understand the unexpected oscillations [
260,
261]. A direct fitting of the resistance in the experiment have also been demonstrated [
258,
259,
261]. Nevertheless, other mechanisms that may lead to the scale invariance oscillations will be topics of broad interest.
In solids, the space group can protect energy nodes with other degeneracies, such as three-, six- and eight-fold one, which may lead to massless fermions that have no counterpart particles in high-energy physics [
262-
265]. Triply-degenerate nodal-point semimetals have been proposed with symmorphic space group symmetry of WC type crystal structure, including TaN, ZrTe and MoP, and observed by angle-resolved photoemission spectroscopy [
266,
267]. The unconventional three component fermions in them are formed by crossing of nondegenerate and double degenerate bands, protected by both rotational and mirror symmetries. As an intermediate fermion between Dirac and Weyl fermion, the host semimetal is expected to have different magnetoresistance [
268-
270]. More topics will also be of broad interest, including the quantum transport in magnetic Weyl semimetals [
56,
271-
279], the kagome ferromagnet Fe
3Sn
2 [
280], double-Weyl semimetals [
20,
281], type-II Weyl semimetals [
51,
282-
291], hopf-link nodal-line semimetals [
292-
295], Aharonov-Bohm effect [
296], quasiparticle interference on the surfaces of Weyl semimetals [
297-
300] and Fano effect [
301], exploring axial-gravitational anomaly through thermoelectrical transport [
302].
Higher Education Press and Springer-Verlag GmbH Germany, part of Springer Nature 2019