1 Introduction
The past decade has been the rapid development of topological materials due to their novel properties and huge potential applications in quantum computing and manufacturing quantum devices. Topological materials contain a variety of confirmed members with more candidates waiting to join. The topological phases were firstly introduced aiming at explaining the properties of the quantum Hall effect (QHE) [
1,
2]. The strong magnetic field in the two-dimensional (2D) electron gas can give rise to Landau-level quantization, and in high quality samples with high mobility the system exhibits the QHE [
2,
3]. The most striking feature in QHE is that chiral edge states existing between Landau-levels lead to the quantized Hall conductance for a range of magnetic field, which can be understood in terms of topological invariants known as Chern or Thouless–Kohmoto–Nightingale–Nijs (TKNN) numbers [
4,
5]. Due to the macroscopic separation of protected chiral edge states, the backscattering is forbidden. Thus, the Hall resistance plateaus survives in macroscopic systems, insensitive to disorder and dephasing effect. However, QHE requires the strong magnetic field, and unfortunately, it is almost impossible to generate such a strong magnetic field outside the laboratory. Alternatively, the quantum anomalous Hall effect (QAHE) was then proposed to harbor gapless chiral edge states with ferromagnetism. Haldane theoretically proposed a 2D honeycomb-lattice QAHE model resulting from breaking of time-reversal symmetry (TRS) without magnetic field [
6]. However, since the 2D honeycomb-lattice model was not realized until 2004 and the manipulation of magnetic flux has been challenging in the condensed matter physics, the first observation of QAHE was reported in Cr-doped (Bi, Sb)
2Te
3 thin films in 2013 [
7].
The 2D topological insulator was firstly predicted in graphene, lately known as quantum spin Hall effects (QSHE) by Kane and Mele in 2005 [
8,
9]. The intrinsic spin–orbit coupling (SOC) converts graphene from an ideal semi-metallic state to a QSH one, which has metallic helical edge states protected by a bulk gap. As opposed to QHE, the electrons in QSH edge channel are spin-momentum locked. They counter-propagate with opposite spins locked perpendicular to their momentums. Thus, the QSHE can be regarded as two copies of QAHE with opposite Chern numbers, which leads to the quantized spin Hall conductance. Furthermore, TRS was found to play a key role in QSHE [
9]. A topological nontrivial state described by TKNN numbers can only occur when TRS is broken, while a QSH state is found to be characterized by a topological
Z2 invariant with TRS [
9]. Two counter-propagating edge states form a Kramers pair in QSHE which are immune to non-magnetic impurity backscattering, regardless of spin-rotation symmetry. Due to TRS and helical edge states, there is no Hall resistance in QSHE, whereas the longitudinal resistance exhibits the quantized plateaus. Though Kane and Mele suggested a structure to realize a QSH phase, the intrinsic SOC in graphene is extremely weak, about
meV, allows the QSHE to occur only at an unrealistical low temperature experimentally [
10-
14]. Fortunately, Bernevig
et al. theoretically proposed that QSHE can be realized in zinc-blende semiconductors such as GaAs [
15] and HgTe quantum wells (QWs) [
16]. Soon after the theoretical proposal, the QSHE is experimentally verified in HgTe QWs [
17], and then observed in InAs/GaSb heterojunctions [
18]. Nowadays, there is a boom in QSHE researches by theoretical predictions [
19-
34], experimental realizations and characterizations [
35-
42] and demonstrations of exotic properties [
43-
45]. And of all the experimental characterizations, the quantized longitudinal resistance was regarded as the smoking gun for the QSHE. However, up to now, the quantized longitudinal resistance can only be observed in mesoscopic samples, which deviates from the quantized value in macroscopic systems [
17,
38,
39].
A topological nontrivial state characterized by Chern numbers can only exist in even dimensions. By contrast, the state featured by a
Z2 topological invariant can be generalized into a 3D version, known as a topological insulator (TI). Fu
et al. firstly predicted that alloy Bi
1-xSb
x would be a 3D TI [
46] and is then observed [
47]. A 3D TI phase has gapped bulk states and 2D helical surface states. However, the surface states in Bi
1-xSb
x turn out to be very complex. The tetradymite semiconductors Bi
2Te
3, Bi
2Se
3, and Sb
2Te
3 are expected to be 3D TIs [
48]. Bi
2Se
3 is stoichiometric, and hence can be prepared as a sample with higher purity than an alloy Bi
1-xSb
x. Hasan’s group observed the single Dirac-cone surface state of Bi
2Se
3 samples [
49,
50]. Then, a pure topological insulator phase without bulk carriers was first manifested in Bi
2Te
3 [
51]. The elegant electronic structure and mature manufacturing techniques greatly promote researches on 3D TIs [
52-
57]. Analogy with 2D TIs, spin-momentum locked helical surface states are also the hallmark of 3D TIs. More intriguingly, the TI surface exhibits a helical spin texture which leads to the exotic transport properties and spintronics applications [
49,
58]. The TRS guarantees the robustness of helical surface states against weak disorder. In common, the helical surface states are stable except the bulk gaps are closed and reopen. Besides, the topological states can be extended into other discrete symmetries [
59,
60].
Dephasing is the process of carriers’ phase coherence loss by definition. In a realistic sample, dephasing process is generally classified into two types: normal dephasing and spin dephasing. Normal dephasing, referring to electron–phonon interactions or electron–electron interactions, etc. [
61,
62], contributes to phase relaxation but does not flip spin of electrons. Conversely, spin dephasing process affects not only phase relaxation but also spin memory. Jiang
et al. investigated the longitudinal resistance in QSHE under two dephasing effects. By phenomenologically numerical study, they found that the quantized longitudinal plateaus is robust against normal dephasing but fragile with spin dephasing [
63]. The theoretical results are in agreement with the first QSHE experiments observation that the quantized longitudinal plateaus can only be observed in mesoscopic samples [
17,
38,
39]. Further transport measurements in InAs/GaSb QWs strongly confirmed that the spin dephasing is the dominant edge scattering process in this system [
44]. Therefore, investigating spin dephasing effects in QSHE is a continuing concern. In addition to a single mechanism, e.g., magnetic impurities or spin fluctuation individually [
64,
65], the latest developments of microscopic mechanisms of spin dephasing rely on the joint effects below: (i) electron–electron interaction and the momentum linear
k-order Rashba SOC for a extended QSH state [
66,
67]; (ii) the
k3-order Rashba SOC term and electron–phonon interaction for a extended QSH state [
68]; (iii)
k-order Rashba SOC and electron–phonon interaction for a bound QSH state [
69]. Next, the dephasing effects on transport properties in 3D TIs are introduced. The 2D spin-momentum locked surface states exhibit a chiral spin texture. Thus, it is meaningless to distinguish normal dephasing from spin dephasing for 3D TIs. Liu
et al. found that the combination of normal dephasing and impurity scattering also leads to dramatic backscattering [
70] and causes surface anomalous gap without TRS broken. Besides, Liao
et al. showed that variable-range-hopping (VRH) mechanism in 3D TI thin films provides a new dephasing mechanism [
71].
The outline of the article is as follows. Section 2 introduces dephasing effects in 2D TIs. The works fall under two categories: phenomenological investigation of dephasing effects on transport properties and latest development of mircosopic mechanisms of spin dephasing. In Section 3, we present how the dephasing process affects the transport properties in 3D TIs. Finally, a brief conclusion and an outlook for future dephasing research are given in Section 4.
2 Dephasing effects in 2D topologicalinsulators
In this section, we focus on the dephasing effects in 2D TIs. First, we introduce dephasing effects on transport properties by phenomenological analysis. Jiang
et al. showed that spin dephasing is the dominant edge scattering process in QSHE [
63]. Next, we present the latest development of microscopic spin dephasing caused by the joint effects of puddles and normal dephasing mainly based on Refs. [
67,
69].
2.1 Phenomenological investigation of dephasing effects on transport properties
A QSH state is featured by helical edge states. Previous theories predicted the quantized longitudinal resistance plateaus [
16]. However, the milestone QSHE experiment shows that the quantized longitudinal resistance plateaus can be only observed for short samples with a size of about 1 μm [
17], which gives close to the value
. When the sample size reaches about 20 μm, the longitudinal resistance deviates from the value
and increases strongly. Afterwards, these phenomena are confirmed by two independent groups in InAs/GaSb QWs [
38] and WTe
2 [
39]. Furthermore, the latter also manifested that the longitudinal resistance increases linearly with the sample length regardless of the sample width. Previous theoretical studies are short of understanding of dephasing effects in QSHE. In the following, we describe that how the dephasing effects influence on the transport properties of a QSH state.
2.1.1 Longitudinal resistance behaviours
First, we review the dephasing effects on the longitudinal resistance in the helical edge states. The work starts from considering a six-terminal device as shown in Fig. 1(a). For simplicity, QSHE is regarded as two copies of QHE with opposite Chern numbers. The Hamiltonian of such a QSHE device in the tight binding model can be described as [
72]
where the first term describes the central region and six terminals.
(
) is the creation (annihilation) operator of an electron on the lattice site
i in the central region,
is the nearest hopping matrix element with the lattice constant
a,
is the spin-dependent phase with
for the real spin up (
) and down (
) due to SOC [
73]. The summation of four
ϕij along each unit satisfies
, here
is an effective magnetic field caused by SOC. The second term represents the Hamiltonian of the ideal virtual leads and the coupling with central regions.
(
) is the creation (annihilation) operator of an electron in the ideal virtual leads connected with site
i. The size of the central region is
as shown in Fig. 1(a). The longitudinal resistance is computed by applying the Landauer–Büttiker formalism combined with the non-equilibrium Green function method [
74-
76]. Dephasing effects only take place in the central region. In the numerical study, the dephasing effect is introduced by Büttiker’s virtual probes [
76,
77] with the dephasing strength Γ, which is directly related to the phase coherence length
[
78]. The total current under dephasing effects contains the contributions of phase-coherent part and phase-incoherent part. The ratio of the phase-incoherent part increases along with the increase in dephasing strength Γ. At a certain value Γ, both parts could have equal percentage, and the system size
L is recognized to be equal to the phase coherent length
. A small bias
between lead 1 and lead 4 is added to drive a current
flowing along the longitudinal direction.
The longitudinal resistance versus for the normal and spin dephasing with several isolated dephasing strength Γ is shown in Figs. 1(b, c). The solid lines in both figures represent the perfect quantized plateaus of a QSH state at when . Here, the index is responsible for multi-channels. When increasing the dephasing strength Γ, the normal dephasing has a slight impact on the longitudinal resistance . The quantized plateaus is almost preserved in Fig. 1(b). That’s a big contrast when introducing the spin dephasing in the central region. The longitudinal resistance in Fig. 1(c) deviates from the quantized plateaus seriously even with small Γ. In conclusion, the quantized plateaus in QSHE is robust against the normal dehasing but fragile with the spin dephasing. The two distinct behaviours of can be related to the concept of backscattering. Since normal dephasing doesn’t flip spin, the carriers only lose phase memories. The backscattering is only allowed between two different edges of the sample. Thus, the backscattering is very weak except when the Fermi energy is near a Landau level center, and the quantized plateaus still sustains even with very large Γ. However, spin dephasing can flip a carrier’s spin. A carrier can be scattered in one single edge. The backscattering is much stronger and the longitudinal resistance is much larger. In that case, the quantized plateaus is broken due to strong backscattering.
In addition, the phase coherence length
is calculated to visualize the destruction of phase relevant to the dephasing strength Γ shown in Figs. 2(d) and (e). Specifically, spin dephasing is common due to magnetic impurities and nuclear spin fluctuations in a realistic system [
64,
65], etc., and hence, the quantized plateaus of longitudinal resistance of QSHE only exists in systems with sample size smaller than spin dephasing length
. This is the reason why the quantized plateau cannot be observed in samples with large lengths [
17,
38,
39]. The plateau of
is preserved under normal dephasing effects regardless of the sample length
L or the width
W [see Fig. 2(b)]. However, for spin dephasing,
is linearly increasing with the sample size
L and almost keeps stable with the disorder
W [see Fig. 2(a)]. It turns out the quantized plateaus cannot sustain for large samples.
approaches the plateau characteristics with the dephasing strength Γ altered, although the plateau may deviate from the idealized quantized one of
without spin dephasing [see Fig. 2(c)]. The feature is consistent with the experimental results. To sum up, the experimental findings from the perspective of the behavior of
and its dependence on temperature and sample size can be explained by spin dephasing effects [
17].
2.1.2 Quantized spin Hall resistance
Since the quantized plateaus of longitudinal resistance
cannot be observed for macroscopic samples due to spin dephasing, it is required to find an observable quantized quantity to reflect the topological nature of the QSHE. Here, spin Hall resistance
is defined as a new quantized quantity and shown in Figs. 3(c) and (d). Notice though, that the quantized plateaus of
Rs is preserved regardless of strong normal or spin dephasing. Specifically, the sample length at
exceeds
by one order of magnitude,
Rs is well quantized. Thus, the plateaus of spin Hall resistance is visible in macroscropic QSHE systems. The robustness of
Rs against both dephasing mechanisms traces an analogy with the Hall resistance plateaus in QHE. The origin is explained below. In Figs. 3(a) and (b), the chemical potential along the boundary is present. The chemical potentials are always spin-independent with
in the left(right) leads. However, the chemical potential
is spin-dependent in the central region. The left leads inject both spin-up and spin-down electrons. Since normal dephasing doesn’t flip spin, the spin-up electron flows into the right lead and the spin-down one goes back to the left lead [see Fig. 3(a)]. When the spin dephasing is introduced, the flow directions of electrons are random due to spin flip. The chemical potential
descends along the longitudinal direction as shown in Fig. 3(b). Nevertheless, the current
necessitates a constant, and hence stable
at any position along the sample and the spin dephasing strength Γ, since the current
is carried by the edge states between
and
. Therefore, the spin Hall resistance
is quantized in macroscopic samples. The voltage difference of
is observable experimentally [
79-
84].
2.2 Latest development of microscopic mechanisms of spin dephasing
In above subsection, we have introduced the significant influences on transport properties in QSHE caused by spin dephasing effects. Here, we describe several joint effects resulting in spin dephasing mechanisms. In the past, it is known that spin dephasing can come from a single mechanism, e.g. magnetic impurities or spin fluctuation [
64,
65]. Apart from the single mechanism, several joint effects which lead to spin dephasing mechanisms in the extended edge states in QSHE generate much interest recently: (1) electron–electron interaction and the momentum linear
k-order Rashba SOC for a extended QSH state [
66,
67]; (2) the
k3-order Rashba SOC term and electron–phonon interaction for a extended QSH state [
68]. Nevertheless, latest experiments study QSHE on a local scale, since puddles widely exist in doped semiconductor hetero-structures due to the constraints of inhomogeneity of donors and acceptors inevitably introduced by sample growth [
17,
43,
85]. These researches identify the enhancement of spin dephasing in the helical edge states by puddles. Thus, we turn to discuss the latest developments of puddles enhanced microscopic mechanisms of spin dephasing which rely on the joint effects.
2.2.1 Puddles combined with electron–phonon interaction enhanced spin dephasing mechanism
In confined QSH QDs, the combination of k-order Rashba SOC and electron–phonon interaction provides an effective spin dephasing mechanism. Hereinafter, Rashba SOC term is treated as perturbation. The whole Hamiltonian of our model, a 4×4 matrix, is written as
A metallic puddle is simulated in QSHE as a QSH QD, described by the massive Dirac Hamiltonian as shown in Fig. 4(a). The above matrix element denote the Hamiltonian for the spin-up part, here is the time-reversal counterpart. The Hamiltonian can be represented in the cylindrical coordinates as follows:
with the mass term
where
m1 and
m2 has positive values and
v is the Fermi velocity. When the sign of the mass
m(
r) is different inside and outside the boundary
, there always exists a bound solution at
[
86,
87]. Thus, one call the mass sign change at
a topological boundary. The bound solution is a pair of helical states circulating around the disk due to the TRS, as shown in Fig. 4(a). Since the
z-component of the total angular momentum
can be a good quantum number [
88], one would like to adopt
to label the helical state in the QSH QD.
is the angular momentum in
z direction.
A Rashba SOC term
hR may come from the axial symmetry breaking [
89] and
is the time-reversal counterpart of
hR. Specifically, the electric potential at the boundary of puddles causes a perpendicular electric field affecting the Rashba SOC, which is present in the cylindrical coordinate with the strength
α below:
The eigenvector of Eq. (2) is taken as . and are spinors associated with the spin-up and spin-down states, respectively. The bound solutions of are solved when and . By introducing the spherical spinor waves into the Hamiltonian in Eq. (3), the Hamiltonian can be separated into two parts. The radial part of becomes
Only bound states near the junction are considered for here. Assume a trial solution . should be positive and real to meet the conditions of . The pure imaginary solutions, when , indicate that the corresponding wave functions spread all over the space. By substituting into the radial Hamiltonian hr, they can obtain the modified Bessel equation of fn and the corresponding spinor , which requires . is the modified Bessel function of second kind. Next, a similar procedure should be adopted in the region . The solution is and with (A, B are the normalization factors). is the modified Bessel function of first kind. By requiring the continuity at , the solution is a transcendental equation of as (for simplicity, they take the mass and thus )
The angular momentum is labeled as in the spin-down case. For a system with TRS, the Kramers pair and have the opposite angular momentum, i.e., , so the transcendental equation for energies of is written as
The energy levels of topological boundary of a QSH QD verse radius R0 are shown in Figs. 4(b, c) numerically. The smaller R0 is, the more discrete the energy levels become. Specifically, they can deduce when . For an arbitrary angular momentum n, the two states have opposite velocities forming a pair of helical bound states (see Fig. 4). Furthermore, the energy spectrum has the Dirac-like form of by defining a new wavevector . Thus, the asymptotic solution can recall the extended helical edge states.
Here, we review how the Rashba SOC and electron–phonon interaction affect backscattering of bound helical edge states in a TI QD. In this model, Rashba SOC is treated as perturbation. Thus, the scattering matrix
can be calculated. The whole scattering matrix contains two parts: the inside region
and the outside region
. To get the nonzero value, the angular part of the integration must meet
and thus
. In other words, the scattering only occurs between two states when their angular momentum differ by one. Though
is nonzero, backscattering is prohibited for a physical system due to the energy conservation when only Rashba SOC is considered. The energy difference of the initial and final states can be overcome by considering phonons. Thus, both the SOC and phonon-electron interaction are essential for inelastic backscattering in our model. The temperature-dependent spin dephasing rate can be used to describe the inelastic backscattering [
90] which is written as
.
is the spin-flip Eliashberg function [
91,
90], here
denotes the energy difference,
ξ is the width of optical phonon energy,
and
represents the optical phonon energy. The local electron–phonon scattering matrix is taken as
with the optical deformation potential
, the atomic mass density
and the lattice constant
a [
92].
So far, the effective spin dephasing due to the joint effect of linear
k-order Rashba SOC and the electron–phonon interaction in local helical edge states with TRS preserved has been demonstrated. Figure 5 shows the spin dephasing relaxation time
verse the radius
R0 in two typical QWs: HgTe/CdTe QWs and InAs/GaSb QWs. The puddles in this model has the size of tens-of-nanometer. Notice that, for an arbitrary temperature
T,
varies with the puddle sizes and has a peak when the level spacing energy approaches the phonon frequency
. The peak shows that the spin dephasing is much stronger and the spin dephasing time is much shorter. Pay attention to the non-monotonic behaviour of the dephasing rate in Fig. 5. The scattering is very weak at both ends of the curves, but has different origins at different end. First, the large energy spacing between neighboring local edge states is scarcely overcome by the optical phonon energy for small puddles. On the other hand, for puddles with large size, the edge states tend to be extended, and this is beyond the scope of our model. In such case, it has been proved that the combination of linear
k-order Rashba SOC and normal dephasing cannot contribute to the spin dephasing in extended edge states [
68]. In a realistic system, there always exist puddles due to inhomogeneity caused by long-range disorder or the external gate. The extended helical edge states at the boundary of sample can couple with the local helical edge states in the puddle several times. Thus, the electrons in extended helical edge states can also undergo the inelastic backscattering through the tunneling.
2.2.2 Puddles combined with electron–electron interaction enhanced spin dephasing mechanism
Except for the puddles combined with electron–phonon enhanced spin dephasing mechanism, the robustness of the quantized conductivity for free-electrons could be destroyed by electron–electron scattering concerning with the puddles as well [
67]. The puddles are stimulated by the QDs (see Fig. 6). TRS guarantees the quantization of conductivity
G0 with value
each edge. Each edge has a pair of the right-moving (R) and left-moving (L) carrier. Whereas there is no backscattering in the extended edge states, electron–electron interactions in the QDs lead to inelastic backscattering. They consider the long-range Coulomb interaction only in QDs which they assume respects the TRS. The long-range Coulomb interaction may cause two backscattering processes: one is for one electron (RR-RL) and the other one is for two electrons (RR-LL), which gives the conductivity correction
and
. Specifically,
and
in low temperature limit [
67]. The helical edge states can go through the inelastic backscattering via the coupling with the QDs via a point contact.
3 Dephasing effects in 3D topological insulators
In this section, we discuss how the dephasing effects influence the transport properties in 3D TIs. One of the hallmark of 3D TIs is the particular chiral spin texture in the 2D surface. Thus, the microscopic mechanisms of spin dephasing in 3D TIs are quite different from the case in QSHE. Rashba SOC is not essential for spin flip due to the intrinsic spin texture [
49,
58]. However, electrons in the helical edge states can only propagate forward and backward along the edge with spin perpendicular to the momenta in QSHE. In a realistic sample, impurities always exist to some degree. Scattering by impurities can change a carrier’s momenta, hence its spin orientation due to spin-momentum locking mechanism. The joint effect of impurities and normal dephasing commonly exist in 3D TIs, and analogy to the 2D case, can cause spin dephasing. Such spin dephasing effects have been observed in recent experiments [
96-
99]. On the other hand, the variable-range-hopping (VRH) transport due to the wide existence of puddles results in a new spin dephasing mechanism. In the following, we introduce the recent developments of dephasing effects in 3D TIs mainly based on Refs. [
70,
71].
3.1 Impurities induced dephasing effects on the helical surface states
3D TIs are characterized by the gapless surface states which are protected by the bulk gap. Recent researchers have focused on the experiments of the alloy BiTl(S
1-xSe
x)
2 [
96-
99]. A topological quantum phase transition (QPT) would take place at
from a trivial state (
) to a topological state (
) when increasing the S/Se ratio
x. It is found that the surface states in BiTl(S
1-xSe
x)
2 acquires an unexpected mass for
which is so-called anomalous Dirac gap at the Γ point. It is thought that a Dirac gap in helical surface states occurs when TRS is broken. However, zero
z-spin component in BiTl(S
1-xSe
x)
2 by spin-resolved ARPES indicates that no magnetic order exists in the samples and thus TRS is preserved. The existing accounts fail to resolve the contradiction between the theory and experimental result. Besides, it reminds us the importance of the impurities in view of the fact that this gap grows when increasing the S concentration [
97]. The anomalous “gap-like” feature arouses the investigation of the microscopic spin dephasing in the helical surface states of 3D TIs with TRS. Liu
et al. found that the phenomenon can be explained by the spin dephasing effects in the helical surface states. Moreover, they intensively studied the spin dephasing effects on transport properties in 3D TIs.
3.1.1 Origin of the normal dephasing
An incident particle is reflected along a random path “1” (
). In the phase coherent regime, there is a time-reversed path “1R” (
). The schematic plot of the two paths is shown in Fig. 7. In a realistic system, an electron is scattered not only by the impurities but also by the phonons or other electrons which are two typical dephasing origins. Once the dephasing mechanism is considered, the quasiparticles which go through path “1” and time-reversed path “1R” gain a random phase difference
. Traditionally, the random phase uncertainty at the finite temperature satisfies
for parabolic-like massive carriers, e.g.,
for electron–electron interaction [
62,
100] and
for electron–phonon interaction [
61,
101]. However, the situation is very different in Dirac electron systems. They show the temperature dependence of
for both dephasing mechanisms below.
Firstly, the electron–phonon interaction induced dephasing effects is considered. The random phase uncertainty
is associated with the temperature-dependent dephasing rate
by relation
.
is the time between two adjacent scatters. The dephasing rate due to electron phonon coupling is determined by the Eliashberg function (introduced in Section 2.2.2) [
101]:
where
T is temperature. For simplicity, only non-Umklapp scattering process is included. The low temperature regime with
and the high temperature regime with
are considered in view of high Debye frequency
. Here,
vL denotes the longitudinal phonon velocity and
q0 represents the longitudinal phonon self-energy correction. In the low temperature regime, the dephasing rate
shows the
T3 dependence as
[
70]. Compared to this, the dephasing rate in the high temperature regime, displays the
dependence as
[
70]. Thus, the random phase
shows the same temperature dependence with
due to Eq. (9).
Next, the electron–electron interaction induced dephasing mechanism is introduced. In particular, the random phase uncertainty
between two counterpart paths is shown in Fig. 7(a). Assume the electrons travel in the environment with a time-dependent electric potential
, so the phase uncertainty can be expressed as
. Electron-electron induced potential has the form
. By applying the RPA approximation [
100], the random phase has the form
[
70]. Here,
σ denotes the conductivity,
vf is the Fermi velocity and
d0 is the skin depth of helical surface states. It is found that the random phase still linearly depends on temperature.
3.1.2 Backscattering cross section
In this section, we discuss the dephasing effects in 3D TIs. The main focus is the metallic surface states which are the hallmark of the 3D TIs here. The helical surface states of the 3D TIs are described by the two-component Dirac Hamiltonian as [
70]
where vf is the Fermi velocity, are Pauli matrix for real spin, and represents the up and bottom surface, respectively. By defining a new wavevector , the Hamiltonian is rewritten as to keep it compact and replace as k for generality. In the following part, only the up surface states with is considered. The eigenstates of the top surface can be expressed in the plane wave basis as , here ± denotes the particle and hole carrier respectively and represents the direction angle. The dephasing effects in 3D TIs by scattering theory are explained below.
Considering the ultrarelativistic Dirac case, the asymptotic scattering spinon wave function is
, where
r=|
r|,
k=|
kin|,
kin denotes the incident wave vector, the
is the scattering amplitude, and the cross section
[
102]. Since
differs
by a
π phase due to the perfect coherence, there is no backscattering theoretically (see detailed derivation below). The quasiparticles which go through path “1” and time-reversed path “1R” gain a random phase difference
concerning with the dephasing effets. By analogy to the case in QSHE, the existence of the random phase uncertainty
cannot give rise to the backscattering individually. In fact, only the combination of impurities and normal dephasing effects contributes to backscattering. In the following part, the above conclusion by calculating the scattering amplitude and backscattering cross-section are confirmed.
The dilute impurities with the potential
V(
r-
Ri) is considered, here
Ri denotes the
i-th impurity location. If the effect of the scatter is not very strong, Born approximation can be used to compute the the scattering amplitude and backscattering cross-section below. The first-order scattering wave function [
70] can be deduced far away from the finite potential range,
where
. The factor
in the first-order Born amplitude is zero due to the
π phase. Even if a random phase
is added to the wave function
, the first-order backscattering is still zero. Therein, no backscattering process takes place in the first-order Born approximation. Next, take the higher-order Born approximation into account. The second-order scattering process is related to the Boltzmannian trajectories in Fig. 7(a). Assume that the average distance
l between two scatters
R1 and
R2 are much larger than the finite potential range, and set the orientation angle of
R2-
R1 as
φ. The second-order scattering wavefunction contains contributions from two inteference paths [
70]:
where the phase factor between two paths is
, and the scattering strength is written as
. The first part in Eq. (11) represents the contributions of path “1” and the
part is related with path “1R” which is associated with the first part by replacing
φ with
and changes
R1 for
R2. The second-order Born amplitude is deduced as [
70]
The phase factor indicates that the backscattering is also eliminated in second-order Born approximation. To simulate the realistic environment, the dephasing mechanism is introduced with a random phase difference by phenomenological analysis. The second-order scattering cross-section after the angle average of φ is
In conclusion, two complementary roles are essential in the backscattering process of the helical surface states: impurities and normal dephasing process. The backscattering process cannot take place with either impurities or dephasing process alone. Notice though, that backscattering only exists in high-order Born approximation.
3.1.3 Anomalous “gap-like” feature
The anomalous gap in the helical surface states is manifested by the co-existence of impurities and normal dephasing process [
70]. In order to understand the experimental results deeply, the effects of the impurities should be investigated in more details. The long-range Coulomb potential
is taken into account, here Δ represents the dimensionless potential intensity. The backscattering cross-section for long-range potential
is inversely proportional to
k2, which indicates that the long-range Coulomb potential scattering occurs around the Dirac point. On the basis of Coulomb scattering, The quasiparticle lifetime is obtained from the transport cross section:
, here
ni is the impurity concentration and
is the total cross section [
70]. Furthermore, the quasiparticle lifetime
τ is included in the imaginary part of self energy, which is related to the Coulomb scattering and dephasing effect by [
70]
where
, the Fermi velocity
m·s
-1, the coupling constant
eV·m, and a moderate dephasing strength with
s is considered [
103]. In Eq. (14), the first part is the scattering lifetime without dephasing, and the dephasing effect is present in the second part. The spectrum function is calculated as
with
. The quasi-particle dispersion of TI surface states by altering the impurity density
ni, is shown in Fig. 7. The quasiparticle dispersion behaves in a linear fashion except the adjacent Dirac point. When the impurity density
ni increases, it seem to form a gap near Dirac point [especially in the Fig. 7(c)]. The band broadening effect is caused by the strong backscattering due to the combination of the normal dephasing and the long-range Coulomb potential. Thus, the anomalous Dirac gap in BiTl(S
1-xSe
x)
2 can be clearly explained. Due to the enhanced backscattering around the Dirac point, weak (anti-) localization in such systems can be largely affected.
3.2 Puddles induced dephasing effects on the helical surface states
The dominant source of phase relaxation is usually electron–electron interactions at low temperatures. However, the situation may be more complicated in 3D TIs. Liao
et al. have reported a new dephasing mechanism in highly tunable (Bi,Sb)
2Te
3 thin films [
71]. The phase coherent transport in such a 3D TI thin films can be tuned continuously from the bulk-conducting regime to the decoupled surface-transport regime. In the bulk-conducting regime, both the bulk and surface states are conducting well. As opposed to this case, the latter refers to the insulating bulk states and decoupled metallic surface states. In the bulk-conducting regime, experimental results have consistently shown the electron–electron induced dephasing rate
. This is so-called Nyquist dephasing mechanism [
104], showing a
T-linear dependence (we also introduced it in Section 3.1.1). Whereas in the decoupled surface-transport regime, the measured dephasing rate has a sublinear power-law temperature dependence (
), deviating seriously from the electron–electron scattering (see Fig. 8). The enhancement of dephasing effects at low temperatures are attributed to the joint effect of the helical surface states and localized charge puddles in the 3D bulk states. Mott’s theory of VRH transport can be used to explain the abnormal behavior of dephasing rate [
105], in which electrons hop between localized states at different locations by thermal activations. VRH transport is a phonon-mediated hopping process, and hence affects the dephasing rate which is present as
with
. Thus,
and
stand for the VRH transport in 2D and 3D electron systems [
71]. Therefore, given the hybrid of 2D surface states and 3D bulk puddles in 3D TI thin films, the theoretical value is consistent with the experimental result in view of dephasing rate.
4 Conclusions and outlook
Thus far, we have provided a brief summary of the latest developments of dephasing effects in 2D and 3D TIs based on reviewing our own dephasing studies and other related investigations. In 2D TIs, the quantized longitudinal resistance is seriously damaged by spin dephasing. Thus, the quantized longitudinal resistance plateaus can only be observed in a sample smaller than spin dephasing length. In particular, the mechanism of spin dephasing in QSHE can originate from the joint effect of puddles and normal dephasing. In 3D TIs, it is meaningless to distinguish normal dephasing and spin dephasing due to the helical spin texture in the 2D surfaces, since any scattering process that changes the momentum can also change the spin direction. Thus, for 2D surface states in 3D TIs, the normal dephasing and spin dephasding are of the same origin. And backscattering takes place within TRS when considering the joint effect of impurities and normal dephasing. In spite of lots of significant progresses, there is abundant room for further studies in dephasing effects in TIs and other topological systems. We discuss the opportunities for dephasing studies in the future before the end of the review.
One of the remarkable theme is the intersect of topological phases and non-Hermitian systems which are open and out of equilibrium. A number of researchers make efforts to seek out topological non-Hermitian phases [
106-
114]. In particular, the system which go through dephasing process is indeed a non-Hermitian one. When applying spin dephasing into the 2D TIs, the quantized longitudinal resistance plateaus would be damaged. Instead, the spin Hall resistance is found to be a macroscopic observable quantity, which enlightens us on the finding of observable quantities in topological non-Hermitian phases.
On the other hand, searching for a more robust QSH state attracts much attention. Recently, Wang’s group realized a “QSH-like” state by a bilayer of QAH insulators with opposite Chern numbers and the insulating layer sanwiched [
115]. The two QAH insulators are applied by different coercive fields, hence have opposite magnetization. As a result, the two edge states counter-propagate with each other harboring opposite spin directions, which are very similar to the helical edge states in QSHE. More intriguingly, the quantized longitudinal resistance plateaus survives over macroscopic length scale even in strong disordered environment due to the spatial separation of two chiral edge states. Besides, it is convenient to tune the potentials of the top and bottom layers, respectively. Therefore, it enables the direct measurements of spin Hall resistance. Furthermore, the systems experience normal dephasing process when a lead is coupled with the single-layer QAH. Since the carrier propagates only in the single QAH layer, it will not change its chirality. The situation is such different when the lead is coupled with the whole sample. In this case, the chiral carriers enter into the lead, randomly lose their spin memories and turn into the opposite QAH layer, which suffers from the spin dephasing process. In other words, voltages can easily tune the types of dephasing effects in such a “QSH-like” state, which makes it a good platform for dephasing effect investigations.
Moreover, metal-insulator transition (MIT) is one of the most classical and important questions in physics. It is known that disorder can lead to the MIT [
116]. When dephasing is applied in such systems, the universal class of phase transitions can be changed [
117]. Questions have been raised about the dephasing effects on the MIT in topological materials [
118-
122]. In recent years, there has been an increasing interest in other topological materials, such as topological crystalline insulators [
123], Weyl semimetals (WSMs) [
124], and higher-order topological insulators (HOTIs) [
125], etc. The dephasing effects on their transport properties, e.g., magnetoresistance, await a similar degree of development.
Higher Education Press and Springer-Verlag GmbH Germany, part of Springer Nature 2019