Topology of electronic bands and Topological Order

Transcription

1 Topology of electronic bands and Topological Order R. Shankar The Institute of Mathematical Sciences, Chennai TIFR, 26 th April, 2011

2 Outline IQHE and the Chern Invariant Topological insulators and the Z 2 invariant Topological order Exactly solved models with topological order A model with everything!

3 Outline IQHE and the Chern Invariant Topological insulators and the Z 2 invariant Topological order Exactly solved models with topological order A model with everything!

4 Outline IQHE and the Chern Invariant Topological insulators and the Z 2 invariant Topological order Exactly solved models with topological order A model with everything!

5 Outline IQHE and the Chern Invariant Topological insulators and the Z 2 invariant Topological order Exactly solved models with topological order A model with everything!

6 Outline IQHE and the Chern Invariant Topological insulators and the Z 2 invariant Topological order Exactly solved models with topological order A model with everything!

7 Outline IQHE and the Chern Invariant Topological insulators and the Z 2 invariant Topological order Exactly solved models with topological order A model with everything!

8 Quantum Hall Plateaus R = (18) n Why is R so accurately measurable? Ω

9 An answer

10 Tight binding models H = iα jβ c iα h(r i R j ) αβ c jβ R i = a i a e a, α, β = 1,..., N B Fourier transform, Reciprocal lattice basis k a c iα = k e i k R i c kα G a e b = δ ab k a Ga, k Ri = a i a k a Hence k a k a + 2π, the Brillouin zone is always a topologically a torus.

11 Momentum space wave functions H = k c α(k) h(k) αβ c β (k) αβ The single particle hamiltonian is an N B N B matrix at every k. h(k) αβ uβ n (k) = ɛ n(k)uα(k) n ɛ n (k), n = 1,..., N B form the energy bands. N B c n (k) (u n (k)) αc α (k) GS = α=1 ɛ n(k)<ɛ F c n(k) 0

12 The phase that launched a thousand scripts The wave functions, u n (k) and e iωn (k) u n (k) represent the same physical state. The phase picked up by the wave-function when the state is Adiabatically transported defines a "Pancharathnam-Berry" connection: u n (k 1 ) e i k 2 k A n 1 i (k)dk i u n (k 2 ) A n i (k) = 1 ( ) (u n (k)) i u n (k) h.c 2i F n ij (k) = ia n j (k) ja n i (k) The Pancharathnam-Berry curvature field, Fij n (k) is independent of the phase convention used to define the eigenstates.

13 The Chern invariant For 2-d systems, the integral of F over the Brillouin zone is an integer invariant, the Chern number. ν n = 1 d 2 k ɛ ij Fij n 8π (k) e2 h Chern Number=Quantized Hall Conductance Thouless, Kohmoto, Nightingale and den Nijs, (1982) Chern Number=Number of chiral edge channels Hatsugai, (1993)

14 The Chern number as a winding number At every k we have a N B level system. eg. N B = 2 we have a 2 level system. h(k) = ɛ(k) + 1 ɛ(k)ˆn(k) τ 2 The physical states at every k form a CP NB 1 manifold. eg at N B = 2, we have the Bloch sphere (CP 1 = S 2 ). ( ) cos θ(k) u(k) = 2 sin θ(k) 2 eiφ(k) The wavefunctions for every band define a map from the Brillouin zone to the space of physical states, CP NB 1.

15 The Chern number as a winding number When d = 2, we have maps from T 2 CP NB 1. These are topologically equivalent to maps from S 2 CP NB 1 since all loops in CP NB 1 can be shrunk to points. The winding number of the map is exactly the Chern number. ν ± = ± 1 d 2 k ˆn 1ˆn 2ˆn 4π = ± 1 d 2 k F 12 (k) 4π

16 Physical Effects: Anomalous velocity Semiclassical equations of motion for wave packets: ẋ i = k i ɛ(k) + F ij (k) k j k i = e x i V (x) + eb ij (x)ẋ j Ganesh Sundaram and Qian Niu, Phys. Rev. B23, (1999); F. D. M. Haldane, PRL 93, (2004) Leads to anomalous hall conductivity: σ H = ɛ(k)<ɛf e 2 2h ɛ ijf ij (k) If the Fermi level is in a gap and all the bands are either filled or empty, then the hall conductivity is quantised.

17 Explicit models on the honeycomb lattice

18 Graphene h(k) = α x p x (k) + α y p y (k) + βm ( ) M p = x ip y p x + ip y M p x (k) = t(1 + cos k 1 + cos k 2 ) p y (k) = t(sin k 1 sin k 2 )

19 Graphene k = K 1 + q, k = K 2 + q, h(k) = t 3 2 (α xq x + α y q y ) + βm F ij (k) = 1 M 4π (q 2 + M 2 ) 3/2 h(k) = t 3 2 (α xq x α y q y ) + βm F ij (k) = 1 4π M (q 2 + M 2 ) 3/2

20 The Haldane model h(k) = α x p x (k) + α y p y (k) + βm(k) p x (k) = t(1 + cos k 1 + cos k 2 ) p y (k) = t(sin k 1 sin k 2 ) M(k) = M + (sin k 1 + sin k 2 + sin( k 1 k 2 )) M >>, ν = 0, M <<, ν = ±1. Topological transition from ν = 0 to ν = ±1 phase at M = Magnetic field not necessary for non-zero Chern number. Other time reversal symmetry breaking terms can also induce it. Bands can exchange units of Pancharathnam-Berry flux when they touch at Dirac points.

21 The Haldane model

22 Outline IQHE and the Chern Invariant Topological insulators and the Z 2 invariant Topological order Exactly solved models with topological order A model with everything!

23 Topological Insulators: Theory

24 Topological Insulators: Real materials

25 Time reversal in Quantum mechanics ψ 2 = e ih(t 2 t 1 ) ψ 1 If there is an operator T such that, T ψ 1 = e ih(t 2 t 1 ) T ψ 2 Then the system is time reversal invariant. Time reversal for Bloch hamiltonians, T u(k) α = (σ y u (k)) α if σ y h (k)σ y = h( k) Then system is time reversal invariant

26 Band Pairs T u n (k) = u n ( k) T u n (k) = u n ( k) ɛ n (k) = ɛ n ( k) Fij n (k) = n Fij

27 The 2-d Z 2 invariant Consider a system with 2 occupied time reversed bands: The total Berry flux carried by them will always be zero. If the hamiltonian is smoothly perturbed, maintaining time reversal invariance, and they touch (h(k) becomes degenerate), it will always happen at pairs of points (k, k). Time reversal invariance ensures that the flux exchanged by bands at these two points is always equal. The change of flux in each band is always even and hence the Chern index of each band modulo 2 is invariant under smooth time reversal symmetric perturbations. in general, is invariant. δ = 1 N B ν n modulo 2 2 n=1

28 The 2-d Z 2 invariant

29 Consequences at the edge Even number of edge pairs (per edge) for δ = 0 Odd number of edge pairs (per edge) for δ = 1 No backscattering for a single pair due to time reversal symmetry.

30 The 3-d Z 2 invariants Parameterise the 3-d torus by π k i π, i = x, y, z. There are many time reversal invariant planes. eg k i = 0, π. The 2-d Z 2 invariants of these planes are all topological invariants. Of these, there are 4 independent invariants which can be chosen to be, δ x δ kx =π, δ y δ ky =π, δ z δ kz=π δ 0 δ kz=0δ kz=π If δ 0 = 1, then the system is called a Strong topological insulator".

31 (δ 0, δ x, δ y, δ z )

32 (δ 0, δ x, δ y, δ z )

33 Consequences at the surface Even number of surface Dirac cones for δ = 0 Odd number of surface Dirac cones for δ = 1 Gaplessness for odd number of Dirac points protected by time reversal symmetry.

34 Outline IQHE and the Chern Invariant Topological insulators and the Z 2 invariant Topological order Exactly solved models with topological order A model with everything!

35 FQHE: Partially filled Landau Bands ν = p 2mp+1 Jainendra Jain

36 Topological order in FQHE For ν = 1/3: unique ground state on the sphere but 3-fold degenerate ground state on the torus. In general FQHE states have a genus dependent degeneracy States characterised by gapped, fractionally charged quasi-particles obeying anyonic statistics. Gapless chiral excitations at the edge. The longwavelength physics described by topological field theories (Chern-Simons theories).

37 Strongly correlated systems

38 Topological entanglement entropy

39 Topological order Topological order realised in quantum Hall liquids and quantum spin liquids. Characterised by: Genus dependent degeneracy of all states Emergent gauge fields Quasi-particles with fractional quantum numbers and statistics Related to pattern of quantum entanglement"?

40 Outline IQHE and the Chern Invariant Topological insulators and the Z 2 invariant Topological order Exactly solved models with topological order A model with everything!

41 Kitaev s honeycomb model Alexei Kitaev, Anyons in an exactly solved model and beyond", cond-mat/ , Annals of Physics. H = J x σi x σj x <ij> + J y σ y i σ y j <ij> + J z σi z σj z <ij>

42 Conserved quantities W p σ1 x σy 2 σz 3 σx 4 σy 5 σz 6 [W p, H] = 0 = [W p, W p ] W 2 p = 1 Conserved quantities have the properties of magnetic fluxes of a Z 2 gauge theory

43 µ il = m<l Jordan-Wigner transformation for S = 1 2 The disorder variables: ( σ z im ) Jordan-Wigner fermions: ξ in = σi x n µ in η in = σ y i n µ in {ξ in, ξ im } = 2δ nm = {η in, η im } {ξ in, η im } = 0

44 Majoranisation H = J x iξ i ξ j + J y iξ i ξ j + J z iξ i u ij ξ j ij ij u ij = iη i η j Hamiltionian of Majorana fermions, ξ i, interacting with static Z 2 gauge fields, u ij, in the gauge, u ij = u ij = 1 ij

45 Generalisations

46 Topological order in the Kitaev Model Fermionic spin-1/2" excitations. Non-abelian anyonic excitations in a particular phase: Unpaired Majorana fermions bound to a vortex (flux defect). 4-fold degeneracy on a torus corresponding to Z 2 fluxes passing through the holes". Saptarshi Mandal PhD thesis

47 Our work Spin liquid: Very short ranged spin-spin correlations. G. Baskaran, Saptarshi Mandal and R. Shankar, PRL 98, (2007) True at all S. Huge classical spin degeneracy: Very much like frustrated quantum antiferromagnets: G. Baskaran, Diptiman Sen and R. Shankar, Phys. Rev. B (2008); D. Dhar, Kabir Ramola and Samarth Chandra Phys. Rev. B First order transitions to magnetically ordered phases (confined phases) when an Ising interaction is added. Criterion that perturbation does not induce long range correlations: Saptarshi Mandal, Subhro Bhattacharjee, Krishnendu Sengupta, R. Shankar and G. Baskaran Simple chains where Majorana Fermions can be manipulated: Abhinav Saket and S. R. Hassan and R. Shankar, Phys. Rev. B 82, (2010) Creating unpaired Majorana fermions by coupling to an impurity spin: Kusum Dhochak and Vikram Tripathi and R. Shankar, Phys. Rev. Lett. 105, (2010)

48 Outline IQHE and the Chern Invariant Topological insulators and the Z 2 invariant Topological order Exactly solved models with topological order A model with everything!

49 Proposals for engineered realisations

50 The Demler et. al. Model H = ij t C i C j + h.c+ ) t (C i σx C j + h.c ij x ++ ( ) t (C i σy C j + h.c )+ ij z t C i σz C j + h.c ij y +U i n i n i S.R. Hassan, S. Goyal and R. Shankar (IMSc.) David Senechal and A -M. Tremblay (Sherbrooke)

51 U = 0 ( ) 0 Σ(k) h(k) = Σ (k) 0 ( ) ( ) Σ(k) = t 1 + e ik 2 + e ik 1 +t σ z + σ x e ik 2 + σ y e ik 1 = Σ ( k) At t = 0, At t = 0, h(k) = σ y h ( k)σ y h(k) = βσ y h ( k)σ y β ( ) I 0 β 0 I At t, t 0, time reversal symmetry is broken.

52 Eigenvalues: t = 0.5

53 Eigenvalues: t = 1.0 The bands stop overlapping at t = 6 3 = 0.717

54 Chern Numbers

55 Pancharatnam-Berry field: t = 0.1

56 Pancharatnam-Berry field: t = 1

57 Chiral edge states

58 U > 0: Variational Cluster Perturbation Theory Choose a cluster hamiltonian H with variational parameters approximating the effect of the environment. (We choose a 6 site cluster and the hopping parameters as the variational parameters). Solve the cluster hamiltonian exactly and the compute the self energy functional, Ω t [G]. Find the saddle point of Ω t to determine G.

59 The Greens Function

60 The Greens Function

61 Conclusions The model proposed by Duan, Demler and Lukin to realise the Kitaev hamiltonian in cold atom systems has time reversal non-invariant spin-orbit couplings except at t = 0 (graphene) and t = 0 For certain parameter regimes, this model has four non-overlapping bands that carry non-zero Chern numbers. This implies a non-zero angular momentum of the ground state at quarter filling. Thus we will have a rotating condensate in a static optical lattice. The gap at quarter filling presists in a region in the U t space. There seem to be two phases, A quantum Hall state and a Chiral Metal. At t = 0, half filling we have a topological insulator. FQHE states at large U, partially filled bands?? (Next problem)

62 THANK YOU!