Lecture 19

04/27/23, Th.

Today

1. Kane-Mele model of 2D TI

Consider two copies the Haldane models, which are time-reversed images of each other. As the Berry curvature or flux changes sign under time reversal,

(1)time reversal symmetry: Ω(k)=Ω(k)inversion symmetry: Ω(k)=Ω(k)

So the combined system has no berry flux through the Brillouin zone. But there will be a pair of degenerate edge states that are Kramer's pair of each other.

image-20230426151447774

Kramers' degeneracy comes from time-reversal symmetry. For spin-1/2, the time-reversal operation is

(2)Θ=iσyK

where K takes the complex conjugate of what ever number follows it. So

(3)Θψk=eikrΘukψk

which is a Bloch function transforms as the irreducible representation (不可约表示) of k, so that

(4)ψk|Θ|ψk=0

If the Hamiltonian has time-reversal symmetry, this means

(5)[H,Θ]=0.

Therefore, ψk and Θψk are degenerate and called a Kramer's pair. This degeneracy is called the Kramer's degeneracy.

This is precisely what the Kane-Mele model does, except the staggered magnetic flux is replaced by spin-orbit coupling. To write down the Kane-Mele model, we need three sets of Pauli matrices

σ : sublattice A,B τ : valley K,K ' s : spin ,

So the low-energy effective Hamiltonian is

(6)h0τ(k)=vF(τkxσx+kyσy)+mKMτszσz=[mKMvF(τqxiqy)00vF(τqx+iqy)mKM0000mKMvF(τqxiqy)00vF(τqx+iqy)mKM]

So the bulk spectrum is two-fold degenerate within each valley τ

(7)ετq=±mKM2+(vFq)2

Initially, Kane and Mele thought there will be a quantized spin current, so they give the name quantum spin Hall effect, describable by a spin Chern number

(8)cs=cc

But the problem is if spin is not good quantum number (with symmetry being broken ever so slightly), so it cannot be a conserved quantity and the quantization gets destroyed.

It turns out, what is important is the Kramer's degeneracy protected by TR symmetry. As edge states ψk and Θψk have exactly opposite spins, if the impurity has time reversal symmetry there will be no back scattering. This pair of counter-propagating edge states related by TR, and are called helical edge states.

It is argued that even with strong disorder, the reflected amplitude must be odd under time-reversal symmetry, so it should be zero. The consequence is that the helical edge states are immune to strong disorder and Anderson localization. From the Landauer formula, we have

(9)σ=2e2h per edge.

The helical edge states protected by the time-reversal symmetry. But what is the bulk topological invariant? The exact formulation is algebraically involved. We will get a first idea with a pictorial argument.

Consider the edge states. We note that if k=k, then obviously the bands at k are two-fold degenerate. The quasimomenta that satisfy this relation

(10)Θk=k+g

where g is a reciprocal lattice vector, are called time-reversal invariant (quasi)momenta (well, apart from an ). Then at these TRIM, the edge states are doubly degenerate due to the Kramer degeneracy. If there spin-orbit interaction, the degeneracy is generally lifted, and each band is spin polarized.

In the case of a 2D insulator, its edge Brillouin zone has two TRIMS: k¯=0,π. The key is, how the doubly degenerate bands at k¯=0 and k¯=π are connected by spin polarized bands. In fact there are only two scenarios, as depicted in the cartoon below.

image-20230426192407905

The case on the left is considered trivial, because it can be continuously deformed into an insulator, by changing the chemical potential for example. The case on the right is the topological. So from the edge state point of view, the topology protected by time reversal symmetry has two distinct classes, entailing a Z2 classification.

In bulk Z2 invariant can be defined with the knowledge of occupied wavefunctions for systems with time-reversal symmetry. It can be simplified for systems with inversion symmetry, for which

(11)(1)ν=k0δ(k0)

where k0 are TRIM of the bulk Brillouin zone, and δ(k0) is the parity of occupied states at k0

(12)δ(k0)=m occ λm

image-20230426193544568

If ν=1, the system is then a topological insulator, with topologically protected edge/surface states. If ν=0, the system is then called a trivial (平庸) insulator.

 

2. Weyl semimetal

A minimal model for a Weyl semimetal can be written again in a two-band Hamiltonian.

(13)H=χ=±1,kΨkhχ(k)Ψk,Ψk=[akbk]

where

(14)hχ(k)=χvF(kχk0)σ

Here χ=±1 corresponding to the chirality of a Weyl node. For example with k0=(0,0,k0),

image-20230427084157160

The 2D Brillouin zone with kz[k0,k0] gives a Chern insulator (quantum anomalous hall insulator), with a chiral edge mode on the surface. Out side this interval, the 2D Brillouin zone for each kz is trivial. So the surface Fermi line terminate at k¯z=±k0, forming a finite surface Fermi arc.

I will stop here our discussion on these topics in view of time, and switch gear to the more traditional topic of another topological system.

 

3. Integer quantum Hall effect

3.1 Phenomenology

The quantum Hall effect is typically induced in a 2-dimensional electron gas (2DEG) by a perpendicular magnetic field.

image-20230426200141802

It is called quantum Hall effect because the Hall conductivity is quantized

(15)σxy=νe2hρyx=1νhe2

where the filling factor ν is an integer.

In a typical Hall measurement, we have current flowing in the longitudinal direction, with the transverse current to be zero (open circuit).

image-20230426223406146

In Drude theory

(16)p˙=e(E+pm×B)pτ.

In steady state p˙=0 so,

(17)[1/τωcωc1/τ][pxpy]=e[ExEy]

where the cyclotron frequency is

(18)ωc=eBm.

Setting py=0, we have

(19)Ey=mωcne2jx=Bnejx.

That is

(20)ρyx=Bne

Note the Hall resistivity is independent of scattering time.

In the quantum Hall measurement, the Hall resistance RH is quantized

(21)RH=ρyx=1νhe2=Bne=Φ/(h/e)Nhe2=NΦNhe2

so the filling factor ν is

(22)ν=NNΦ

the number of electrons per flux quantum threading the system.

image-20230426225843100

Von Klitzing et al. Phys. Rev. Lett. 45, 494 (1983).

image-20230426230739398

At the Hall plateau

The question is: why the precise quantization? By now, you probably have guessed it is the same thing as the quantum anomalous Hall effect. And you are right. But let's see how this can be understood.

 

3.2 Disorder is necessary

A curious thing about the integer quantum Hall effect is that disorder is important! Let's consider a perfect 2DEG that is translationally and rotationally invariant. In the lab frame

(23)E=0B=Bz^

There is no current.

Now in a moving frame with a constant velocity v, there appears to be current

(24)j=nev;

and the electric and magnetic fields of the Lorentz boost is

(25)E=v×BB=B

such that

(26)E=Bnej×z^ρ=Bne(0+110)

There can be no quantization! As the argument only involves density, it applies whether the theory is quantum or classical.

To have quantization, we must break the Lorentz covariance. We need inhomogeneity (disorder) to break the translational or rotational invariance.

 

3.3 Landau levels

Choosing the Landau gauge for B=Bz^

(27)a=Bxy^

the single-electron Hamiltonian reads

(28)H=12m[px2+(pyeBx)2]

which is translationally invariant in y-direction. So the wavefunction can be written as

(29)Ψ(x,y)=eikyψk(x)
(30)HΨ(x,y)=eiky12m[px2+(keBx)2]ψk(x)Hk=12m[px2+12mωc2(xxk)2]

where

(31)cyclotron frequency:ωc=eBmmagnetic length:lB=eBguiding center:xk=lB2k

So Hk looks exactly like a Hamiltonian of a 1DHO centered at xk with a characteristic frequency ωc. So the energies of these discrete Landau levels are

(32)εnk=(n+12)ωc.

which is highly degenerate as it is k independent.

image-20230427085416641

We use n as the quantum number for the energy levels of the 1DHO, and the wavefunction is

(33)ψnk=1Aexp[(xxk)22lB2]Hn[(xxk)lB]

where Hn(z) is the Hermite polynomial.

image-20230427085454227

What exactly is the degeneracy of a Landau level? We need xk[0,L],

(34)0<k<LxlB2

States on the left and right edges can have vastly different k values. Noting k is ky, the number of states in a LL is

(35)N=Ly0Lx/lB2dk2π=LxLy2πlB2=LxLyh/eB=NΦ

Thus, there there is one state per flux quantum in an LL. So when ν takes integer values, the Landau levels are either fully filled or empty in the ground state. Thus, the quantum Hall states are insulating, incompressible, like a band insulator.

Consider a wavepacket and its time evolution is clearly

(36)Ψ(r,t)=Ly2πndk an(k)ψnk(r)ei(n+12)ωct

Now consider

(37)Ψ(r,t+2π/ωc)=Ψ(r,t)

so the motion is clearly periodic, like the classical cyclotron orbit motion.

3.4 Current response

Now we want to calculate the current response of a quantum Hall system. We will restrict ourselves to the lowest Landan level (LLL)

(38)ψk(r)=1π1/2LylBeikye122(xklB2)2

The current density for a given k is

(39)jk(r)=emψk(r)[i+ea(r)]ψk(r)

So the current in y-direction is

(40)Iy(k)=emeπ1/2LylBdx e12lB2(xxk)2(keBx)e12lB2(xxk)2=eωcπ1/2lB1Lydx e1lB2(xxk)2(xxk)=0

Here, the canonical momentum in y-direction k is exactly cancelled by the vector potential, so the mechanical momentum vanishes.

There is no current in the ground state, for sure. This is expected since εnk is independent of k, so the electron group velocity is zero, in any LL.

But, this independence of energy on k comes from the fact that energy is independent of xk. This can be immediately changed if an electric field is applied

(41)V(r)=eEx,

which also preserves translation invariance in y-direction.

image-20230427093232949

The Hamiltonian can be re-written

(42)Hk=px22m+12mωc2(xxk)2+eEx=px22m+12mωc2(xxk)2+eExk+12mvd2

where we have new guiding center

(43)xk=xkeEmωc

and introduced a drift velocity

(44)vd=EB

So the LL energies become

(45)εnk=(n+12)ωc+eExk+12mvd2

So the current becomes

(46)Iy(k)=eωcπ1/2lB1Lydx e1lB2(xxk)2(xxk)=eωcπ1/2lB1LyπlB(xkxk)=1Lyeωc(eEmωc2)=evdLy

So the current density in y direction is

(47)jy=LyLx0Lx/lB2dk2πIy(k)=evd1LxLx2πlB2=e2hE

So

(48)σyx=e2h

So the filled LLL has a quantized Hall conductivity.

We will then try and understand the topological nature of the quantized Hall effect, by figuring out the bulk topological invariant and the corresponding edge states.