Tight-binding models tight-binding models. Phys540.nb Example 1: a one-band model

Transcription

1 Phys540.nb 77 6 Tight-binding models 6.. tight-binding models Tight-binding models are effective tools to describe the motion of electrons in solids. Here, we assume that the system is a discrete lattice and electrons can only stay on the lattice site. The kinetic energy is included by allowing electrons to hop from one site to another Example : a one-band model Lets consider a D lattice with one atom per unit cell. For each atom, we consider only one quantum state. The creation (annihilation) operator c i (c i ) creates and annihilates one particle on site i. The non-interacting Hamiltonian can be written as H =- ij t ij c i c j + t ji c j c i + i V i c i c i (6.) The first term ij t ij c i c j describes hoppings from site j to i. The second term ij t ji c j c i describes the hopping from i to j. The last term is the potential energy, which tells us how much energy we need to put an electron on each site V i. Because the lattice contains only one type of * atoms, V i = constant (i.e. translational symmetry). Because H is Hermitian, we find that t ij = t ji and V i is a real number. Proof: H = H - ij t ij c i c j + t ji c j c i + i V i c i c i =- ij t ij c i c j + t ji c j c i + i V i c i c i (6.) (6.3) - ij t ij * c j c i + t ji c i c j + i V i * c i c i =- ij t ij c i c j + t ji c j c i + i V i c i c i (6.4) By comparing the two sides, we find that t ij * = t ji and V i = V i *. Therefore, we can simplify the Hamiltonian H =- ij t ij c i c j + t ij * c j c i + V i c i c i =- ij t ij c i c j + t ij * c j c i + V N (6.5) In the last term, N = i c i c i is the total number of electrons in the system. Because VN is a constant, this term just shifts the total energy by a constant, and thus has no other physical contribution (can be ignored if we are not interested in the total energy). For tight-binding models, a typically approximation is to assume that electrons can only hop to its nearest-neighbor sites. In reality, long-range hopping is allowed, but their amplitudes are small (decay exponentially as distance increases). Therefore, in many cases, we just need the nearest-neighbor hopping terms to describe the systems. Due to the translational symmetry, all nearest-neighbor hoppings shall have the same hopping strength. If we assume that this hopping strength is a real number t, the Hamiltonian is H =-t ij c i c j + c j c i =-t i c i c i+ + h.c. (6.6) Here, ij implies that i and j must be neighbors. On the right hand side, h.c. means Hermitian conjugate. Notice that the second term is in fact the Hermitian conjugate of the first term, so we just use h.c. to represent it.

2 78 Phys540.nb Fourier series: typically, Fourier series is used to describe a periodic function in the real space, which will have a describe set of wave vectors in the k space. Here, it is the opposite. We have a continuous k-space and it is periodic (the Brillouin zone), but the real space is discrete. c k = a p - kx c i i (6.7) c i = a p kc k  kx (6.8) Because c i, c j =d ij, it is easy to check that c k, c k' = c i - kxi, i c j  k' x j = j pa pa c i, c j i, j pa - kxi  k' x j = d i, j i, j pa - k-k' xi = i pa - k-k' ia = a dk - k ' a =dk - k ' (6.9) If we know c k, c k' =dk - k ', we can also show that c i, c j =d ij. c i, c j = kc k  kxi, p a p a k ' c k' - k' x j = pa k k ' c k, c k'  kxi - k' x j = pa k k ' dk - k '  kxi - k' x j = pa k  k xi-x j =d ij Hamiltonian in k-space We can transfer the Hamiltonian into the k-space. In real space, it looks like H =-tc i c i+ + h.c. i (6.0) (6.) The first term is j c j c j+ = j p a kc k - ka j p a k ' c k'  k' j+ a = j p a k k ' c k c k'  k' a - k'-k aj = k k ' c k c k'  k' a dk - k ' = kc k c k  ka The second term is the Hermitian conjugate of the first term, so h.c. = c j c j+ = j kc k c k  ka = kc - ka k c k (6.) (6.3) H =-t i c i c i+ + h.c. =-t kc k c k  ka - t kc k c k - ka =- t kc k c k cos ka= k - t cos ka c k c k (6.4) For a solid, we know that the total energy of electrons (ignore interactions) is E = n k e n k n n k (6.5) where n sums over all bands, e n k is the dispersion relation for band n and n n k is the occupation number for the Bloch wave state in band n with momentum k. In second quantization, this formula implies that the Hamiltonian is H = n k e n k g n,k g n,k (6.6) Here, g n,k is the creation operator for a Bloch wave y nk r = u nk r  kr. If we compare this formula with the tight-binding Hamiltonian (in k- space), we find immediately that the tight-binding model we considered here has only one energy band. And our c k operator is in fact the creation operator for Bloch waves. And the dispersion relation for this band is e k =-tcos ka. Notice that e k is a periodic function of k with periodicity pa, which is exactly what we expect for Bloch waves Discrete Fourier transformation and Fourier series

3 Phys540.nb 79 In the formula above, we treat k as a continuous variable. In many cases, it is useful to treat k as discrete values (where k turns into k ). Here, we briefly discuss about the discrete Fourier transform. Consider a discrete function f i, where i =,, 3 N marks different lattice site. In addition, we assume the periodic boundary condition f N+i = f i. We can define f è k = N - ka j f j j (6.7) where a is the lattice constant and j marks different lattice sites. And it is easy to prove that f j = N f è ka j k k  (6.8) Same as j, the wave vector k here is also a discrete variable. This is because f N+i = f N (our system has a finite size) f j+n = N k f è k  k a j+ an = N k f è k  ka j  kan (6.9) f j = N f è ka j k k  (6.0) If we compare the two equations, we find that  kan =, which implies k = p m N a = p m L Here L = Na is the size of the system. For infinite systems, L, the discrete sum k turns into an integral k. (6.) Same as j, k also have a periodicity, and the periodicity is also N. This periodicity is just the Brillouin zone. Define k m = p m L and by definition, f è k m = N j f j - k m aj = N j f j - p Na m aj = N f j - p m j N j (6.) For k m+n, by definition f è k m+n = N j f j - k m+n aj = N j f j - p Na m+n aj = N f j - p m j N - p j j (6.3) Because j is an integer, the factor - p j =. Therefore f è k m+n = N f j - p m j N = f è j k m (6.4) As a result, we can limit the value of m to be -N m < N. If m is not in this range, the value of f è k m can be obtained using the periodic condition f è k m+n = f è k m m =- N, - N +, N -, (6.5) For the wave vector k, this means that f è k+ pa = fè k (6.6) So, we can confine the value of k into the range of -pa k <pa, which is the first Brillouin zone. For k outside the first Brillouin zone, we can find the corresponding f è k using the periodicity fè k+ pa = fè k k =- p a, - p a + μ p L, - p a + μ p L, p a - p L, Two useful Identities: (6.7)

4 80 Phys540.nb N k  ka j - ka j' =d j, j' N j  ka j - k' aj =d k,k' (6.8) (6.9) Example : a two-band model in D Now, let us consider a slightly more complicated situation, i.e. a d chain formed by two different types of atoms a and b). H =-t i a i b i + b i a i+ + h.c. + V a i a i a i + V b i b i b i (6.30) Here, the position of a sites in the jth unit cells is r = aμ j + r a (6.3) and the position of b sites in the jth unit cells is r = aμ j + r b (6.3) where, a is the lattice constant (the size of a unit cell) Discrete Fourier transform (as discussed in the previous section) a k = N i a i  kx (6.33) a i = N - kx a k k (6.34) b k = N i a i  kx (6.35) b i = N - kx a k k (6.36) Therefore, we know that j a j b j = j N k k' a k b k' k a k - k a j+ra N b k' k'  k' a j+rb = N a k k' k b k' - k a j+ra  k' a j+rb = j N j - k-k' aj - kra  k' rb = k k' a k b k' d k,k' - kra  k' rb = k a k b k  k rb-ra This conclusion is generic. In k-space, hoppings from b to a result in a term in Hamiltonian k a k b k k rbra. The phase factor here is determined by how far the electron hop, i.e. r b - r a. If we apply the same conclusion to the hopping from a to b, we find immediately that j b j a j+ = k b k a k  kr a-r b If r b - r a = a, (6.37) (6.38) j a j b j = k a k b k  ka j b j a j+ = k b k a k  ka (6.39) (6.40) For the potential term, we also know that i a i a i = k a k a k (6.4) and

5 Phys540.nb 8 i b i b i = k b k b k (6.4) Therefore, in the k-space, the Hamiltonian looks like H =-t a k k b k  ka + b k a k  ka + b k a k - ka + a k b k - ka + V a a k k a k + V b b k k b k = - t a k k b k cos ka + b k a k cos ka + V a a k a k + V b b k b k We can write it in a matrix form: (6.43) H = k a k, b k V a - t cos ka - t cos ka V b a k b k For any tight-binding models with two quantum states per unit cell, the Hamiltonian can be written in terms of a two-by-two matrix in the k space: (6.44) H = k a k, b k Hk a k b k (6.45) where H is a Hermitian matrix as a function of k. It is called the kernel of the Hamiltonian. Because Hk contains all the information of the Hamiltonian H, it is often called the Hamiltonian in literature. However, it is important to keep in mind that Hk is only part of the Hamiltonian. The Hamiltonian H must be gauge invariant, but Hk is not. For example, if we change a a Âf, the Hamiltonian H is invariant, but the kernel Hk is NOT. For the model considered here, Hk = V a - t cos ka - t cos ka V b (6.46) For more generic cases, if one have m quantum states per unit cell, H(k) will be a m m Hermitian matrix. When we have a matrix, we know what to do. We need to find the eigenvalue and eigenfunctions of the Hamiltonian, i.e. diagonalizing the matrix. The next two sections discuss the physical meanings of the eigenvalues and eigenvectors of Hk Eigenvalues of Hk The eigenvalues of Hk give us the dispersion relations. At each k point, one can define a unitary transformation c k d k = U k - a k b k c k d k = a k b k U k (6.47) (6.48) where U k - = U k H = k c k, d k U k + H U k c k d k (6.49) If we choose U k such that U k H U k is a diagonal matrix, U k H U k - = e c 0 0 e d (6.50) Then, the Hamiltonian becomes H = k c k, d k e c 0 0 e d c k d k = k e c k c k c k + k e d k d k d k (6.5) We know that using Bloch waves (which are eigenstates of the Hamiltonian), the total energy is E = k,m e m k n m k = k e n k + k e n k + k e 3 n 3 k + (6.5) where e m k is the dispersion relation for the mth band and n m k is the occupation number for the quantum state in the mth band with momentum k. Therefore, in terms of Bloch waves, the Hamiltonian should take the form:

6 8 Phys540.nb H = k e k g,k g,k + k e k g,k g,k + k e 3 k g 3,k g 3,k + (6.53) where g n,k is the creation operator for a Bloch wave in band n at momentum k. If we compare this formula with Eq. 6.5, we find that in this model we get two energy bands. c k creates a Bloch wave with momentum k in one of the bands. d k creates a Bloch wave with momentum k in the other band. For the bands created by c k, its dispersion relation is e c k, while the other hand has a dispersion relation e d k. For the problem we considered here, the dispersions are e c = V a + V b - t cos ka + V a - V b (6.54) e d = V a + V b + t cos ka + V a - V b which are the eigenvalues of the matrix Hk. (6.55) If we have a model with m orbitals per unit cell, Hk would be a m m matrix. This matrix will have m eigenvalues, which results in a m-band model and each eigenvalue e n k is the dispersion relation for one of the bands. Bottom line: eigenvalues of H(k) tell us the dispersion relation Eigenvectors of Hk The eigenvectors Hk gives us the Bloch waves. Let s start from the two band model discussed above as an example. Here, the two-by-two matrix Hk has two eigenvalues, e c k and e d k. For each eigenvalue, the corresponding eignevector is a -component vector: H v- ak v - b k =e c v - ak v - b k (6.56) H v+ ak v + b k =e d v + ak v + b k (6.57) It can be easily checked that c k = v - ak a k + v - bk b k d k = v + ak a k + v + bk b k (6.58) (6.59) We know that a and b are the creation operators for plane waves (because they come from Fourier transformations) and we also know that c k and d k create Bloch waves. The relation between plane waves and Bloch waves is y n,k x = u n, k x  kx N (6.60) So, the coefficients v a and v bk here are in fact the discrete version of u n, k x. Bottom line: eigenvectors of H(k) tell us u n, k. We have learned that the Berry connection A n k is defined as A n =-Âun,k u n,k k (6.6) For continuous models, this formula means A n =- x u n, k x * u n, k x k (6.6) For lattice models (tight-binding models)

7 Phys540.nb 83 A - =-Â v - a k * v - bk * k A + =-Â v + a k * v + bk * k v - ak v - b k v + ak v + b k (6.63) (6.64) In general, for a m band model, where Hk is a mμm matrix, each of the eigenvector gives us the Bloch wavefunction for one of the m bands. We can use each eigenvector to compute the Berry connection, the Berry curvature and the Chern number for that band Bottom line: the eigenvectors are Bloch wave functions Summary The key properties of a tight-binding model is coded in H k, which is a mμm Hermitian matrix. Each component is a function of k. (In other words, H k is a matrix function of k.) The eigenvalues if H (as a function of k) gives the band structure (the dispersion relation) for each bands e n k with n =,, m The eigenvectors as a function of k gives the Bloch wave: u n,k a, where n=,,, m is the band index and a=, m marks different site in a unit cell. Using these u n,k a, we can define the Berry connection, Berry curvature and the Chern number. Just replace the integral in the real space by summing over a= m Before 80s, people stop at step # (computing the eigenvalues), without calculating the eigenvectors. This is because eigenvectors are the wavefunction, which cannot be measured directly in experiments. After the discovery of topological insulators, physicists realized that although the wavefunction can not be measured directly, it contains the topological information, which is a physical observable (e.g. the Hall conductivity) One band model In a one-band model, the eigenvector is trivial, which is just identity. Therefore, the Berry connection A n =-Â x u n, k x * u n, k x must k be trivial A n = 0. This tells us that a one-band model can never show quantum Hall effect. We need at least two bands. In fact, one can prove that the total Chern number summing over all bands must be zero. 6.. an example of a topologically nontrivial insulator Let s consider a -band model, whose kernel is Hk = Dsin k x +Âsin k y t cos k x + t cos k y +m - t cos k x - t cos k y -m Dsin k x -Âsin k y (6.65) Here, we assume D>0 and t>0 and -4 t <m<4 t. In real space, this Hamiltonian corresponds to a square lattice and on each site there are two quantum states. Please notice that here I use m to refer to a control parameter, which is NOT the chemical potential Two-band models and Pauli matrices For any two-band models, H is a two-by-two Hermitian matrix. For a two-by-two Hermitian matrix, one can always separate it into the identity and Pauli matrices I, and s i H = H 0 k I + H x k s x + H y k s y + H z k s z where H 0 k, H x k, H y k and H z k are real functions of k. If we define Hk = H x k, H y k, H z k, we find that (6.66) H = H 0 k I + Hkÿs (6.67) This Hamiltonian is rather similar to a spin S=/ under magnetic fields. H = constant +mb ÿs (6.68) For the case studied here, H 0 k = 0 (6.69)

8 84 Phys540.nb H x k =Dsin k x H y k =Dsin k y H z k =- t cos k x - t cos k y - m (6.70) (6.7) (6.7) The eigenvalues of H are E = H 0 k Hk = H 0 k H x k + H y k + H z k (6.73) This is pretty much the same as the spin system, where the eigen-energies of the are E = constant m B (6.74) For a two band model, the top band has energy E + H 0, while the lower band has E - H 0. The energy gap between the two bands is: Dk = E + k - E - k = Hk (6.75) As long as H k 0, the two band will not cross with each other. Let s now focus on the lower band E - ), its eigenvector is u I -k = N I H z k - Hk H x k +ÂH y k (6.76) This wavefunction is singular if H x k = H y k = 0 and H z > 0, u I -k = N I H z - Hk H x k +ÂH y k = N I H z - H z k 0 = N I H z - H z k 0 = N I 0 0 (6.77) This is indeed a problem for the Hamiltonian shown above. At k x = k y =p. H 0 k = 0 H x k =Dsin k x = 0 H y k =Dsin k y = 0 H z k =-tcos k x - t cos k y -m=4t-m>0 (6.78) (6.79) (6.80) (6.8) Remember that we assumed m < 4 t. In fact, there is another way to write down the same eigenvector (i.e. a phase shift) u II -k = N I H z - Hk H x k +ÂH y k μ H z+ Hk H x k+â H y k H z+ Hk H xk+â H yk = N II -H xk -H yk H xk+â H yk H z + Hk = N II -H x k +ÂH y k H z + Hk (6.8) The two wavefunctions u I -k and u II -k differ by a phase u II -k = u I -k Âfk (6.83) where Âfk = H z + Hk H x k+â H y k H z+ Hk H x k+â H y k (6.84) This new wavefunction is well defined at H x k = H y k = 0 and H z k > 0. However, it is NOT well defined at H x k = H y k = 0 and H z k < 0. There, we find

9 Phys540.nb 85 u II -k = N II -H x k +ÂH y k H z k + Hk = N II 0 = 0 H z k + H z k N II H z k + H z k = N II 0 H z k - H z k = 0 (6.85) For the Hamiltonian shown above, the origin is such a point. k x = k y = 0 H 0 k = 0 H x k =Dsin k x = 0 H y k =Dsin k y = 0 H z k =-tcos k x - t cos k y -m=-4t-m<0 (6.86) (6.87) (6.88) (6.89) Therefore, we need to cut the into two areas and use two different wave functions to describe the Bloch waves. They are connected by a gauge transformation u II -k = u I -k Âfk (6.90) A - II k = A - I k + k fk (6.9) This is in strong analogy to the magnetic monopole case, where we need to different gauge fields. Some comments: These singularities are NOT physical. If one measures any physical observables, there is no singularity anywhere in the momentum space. However. for the wavefunction and the Berry connection, which are not physical observables, there is always some singularity for this Hamiltonian. The location of these singularity points depends on the gauge (phase) choice. In other words, the location of the singularities has no physical meaning either. Only one thing about these singularities is physical, its existence. There must be some singularity points. This statement is independent of gauge choice and it tells us that the topological index is nonzero. Let s compute the topological index for this model. For these model, there are four special points which satisfy H x = H y = 0. They are k = 0, 0, k = p, p, k = 0, p and k = p, 0. At these four points, the values of Hz are: H z =-4t - m, 4 t - m, -m and -m respectively case I: m<-4 t If m<-4 t, all the four special points has H z > 0. So we can use u II -k for the whole Brillouin zone, and there is no singularity points. u II -k = N I H z k - Hk H x k +ÂH y k (6.9) Then we can get Berry connection A - II =-Âu II -k k u II - k (6.93) The Berry curvature for the lower band is W - = μa - (6.94) Thus. the total Berry curvature (p times the Chern number) is k W - = k μa - = k A - = 0 (6.95) So, the Chern number is zero C = p k W - = 0. (6.96) Marginal case: m=-4 t H x k = H y k = H z k = 0 at k = p, p Because the gap between the two bands is

10 86 Phys540.nb Dk = E + - E - = Hk (6.97) Dk = 0 at k = p, p. The gap closes at k = p, p, i.e. the two energy bands touch with each other at k = p, p. This gives us a Dirac point. Here, the system is not an insulator, but a metal (a semi-metal). So we cannot define a topological index. Remember that we showed early on, to change the topology of an insulator, one need to close the gap. And here is the gap closing case II: -4 t <m<0 For -4 t <m<0, H z < 0 at k = 0, 0, and H z > 0 for the other three points k = 0, p, p, 0 and p, p. Therefore, we need two wavefunctions. First, we draw a small circle around the origin. Inside this small circle, which we will call region D I, we use u I -k = N I H z k - Hk H x k +ÂH y k Outside the circle, which we will call region D II, we choose (6.98) u II -k = N I H z k - Hk H x k +ÂH y k (6.99) The wavefunctions for these two regions are connected by a gauge transformation u II -k = u I -k Âfk (6.00) And the Berry connection are related also by the same gauge transformation. A - II k = A - I k + k fk (6.0) The total Berry curvature (p times the Chern number) is k W - = DI k μa - I + DII k μa - II = DI k A - I + DII k A - II = DI k A - I - DI ka - II = DI ka - I - A - II =- DI k k fk =- 0 p q q f=-fq= p -fq=0 =fq=0 -fq= p (6.0) We know that f q=0 -fq= p is quantized fq=0 -fq= p = p n (6.03) Therefore, the Chern number is quantized. For the model we considered here: H 0 k = 0 H x k =Dsin k x H y k =Dsin k y H z k =-tcos k x - t cos k y - m (6.04) (6.05) (6.06) (6.07) If we choose D I to be a very small circle with radius k ~0, then around the circle, we can expand everything as a power series of small k. H x k =Dsin k x ºDk x + Ok H y k =Dsin k y ºDk y + Ok H z k =- t cos k x - t cos k y -m=- t -m+ok (6.08) (6.09) (6.0) Âfk = H z + Hk H xk+â H yk H z + Hk H x k+â H y k = H xk+â H yk H x k+â H y k = H x k +ÂH y k H x k +ÂH y k = D k x +ÂDk y D k x +ÂDk y = k x +Âk y k x +Â k y = k Âq k Âq = q = -Â Âq (6.)

11 Phys540.nb 87 Here, we transfer k x +Âk y into polar coordinate, which is k  q. This result tells us that fk =-q (6.) Therefore, k W - =fq=0 -fq= p = p = p (6.3) So, C = p k W - = This is a topological insulator with C =. (6.4) Marginal case: m=0 The gap closes at k = 0, p and p, 0. The system is not an insulator, so we cannot define the Chern number case III: 0 <m<4 t For 0 < m < 4 t, H z > 0 at k=(p,p), and H z < 0 for the other three points. Here,we can draw a small circle centered at (p,p) to cut the system into two parts. Inside the circle, near (p,p), we use u II and outside the circle, near (0,0) we use u I. Same as case II, we get a topological insulator with C = Marginal case: m=4 t The gap closes at k = p, p. Not an insulator case IV: m>4 t All the four special points has H z < 0. So we just use u I for the whole. No singularity, and thus C = the top band? The top band has the opposite Chern number C + =-C -. This is because for any tight-binding models, the total Chern number for all the bands is always 0. So here, we have C + + C - = Summary This model has four phases. Two topological phases with C=+ and two trivial insulator phase with C=0. A topological phase and a trivial insulator phase are always separated by a phase transition, which is known as a topological phase transition. Across the topological phase transition, the topological index changes its value Across a topological transition, the insulating gap closes and then reopens (generically true) Gap closing is the necessary condition for a topological transition, but it is not sufficient. One may close and reopen the gap without changing the topological index. e.g. the m=0 point here The model of Haldane on a honeycomb lattice the Honeycomb lattice (graphene)

12 88 Phys540.nb Fig.. a honeycomb lattice. Two sites per unit cell (two sublattices) a and b (red and black respectively). The green arrow indicates the lattice vectors v and v. If we shift the lattices by m v + nv with m and n being integers, the lattice is invariant. Vr = Vr + mv + nv (6.5) Define a to the length of the nearest-neighbor (NN) bonds. Then we can show that v = 3 a, 0 and v = - 3 a, 3 a. If we only consider hoppings between the nearest-neighbor (NN) sites, the Hamiltonian is H =-t <i, j> a i b j - t <i, j> b i a j (6.6) where i,j means nearest neighbors Band structures There are three type of NN bonds: () along the y axis with q=p/, () along q=p/+p/3=7p/6, (3) along q=p/+4p/3=p/6, and we need to write them out separately. H =-t a i i r b ri +e - t a i i r b ri +e - t a i i r b ri +e + h.c. 3 with e = 0, a e = - 3 a, - a (6.7) (6.8) (6.9) e3 = 3 a, - a (6.0) Here the sum i sums over all unit cells (or say all red sites). If we go to the momentum space, by performing the D Fourier transformation, a k = N i a i  k ÿr (6.) a i = N k a k - k ÿr (6.) b k = N i a i  k ÿr (6.3) b i = N k a k - k ÿr (6.4)

13 Phys540.nb 89 the first term in the Hamiltonian becomes -t i a ri b ri +e =- t N i k k' -t k k' a k b k' - k' ÿe N i  k a k  k ÿr i bk' -k' - k' ÿr i+e = ÿr i =-t k a k' k b k' - k' dk,k' =-t k a k b k ÿe - k' ÿe (6.5) Notice that, as we proved early on, we get a phase factor expâ k ÿr a - r b, which only depends on the separation between a and b sites. If we repeat the same procedure for the other two bonds, we get H =-t i a ri b ri+e - t i a ri b ri+e - t i a ri b ri+e + h.c. =-t a k k b k - k ÿe + - k = k a k b k + - k ÿe 0 H k H k 0 ÿe 3 a k b k - t k b k a k  k ÿe +  k +  k ÿe ÿe 3 (6.6) H k =-texp-â k ÿe + exp-â k ÿe + exp-â k ÿe3 (6.7) H k = H k * =-texpâ k ÿe + expâ k ÿe + expâ k ÿe3 (6.8) The kernel of the Hamiltonian: Hk = 0 H k H k 0 (6.9) The eigenvalues of H(k) gives the dispersion relation. e k = H k = t 3 + cos 3 k x a + 4 cos 3 k x a Cos 3 k y a (6.30)

14 90 Phys540.nb Fig.. The contour plot of the energy dispersion for the lower band (e - as a function of k x and k y ). The red dash lines mark the first, which is a hexagon. The dispersion e is a periodic function of k space (the hexagon repeats itself in the figure shown above). For the lower band, its energy minimum is located at k = 0. The maximum of e - is reached at the corners of the. The first is a hexagon. It has six corners. However, due to the periodic structure in k-space, the three corner points with q=0, q=p/3 and q=4p/3 are the same point (their momenta differ by a reciprocal vector). Similarly, the three corners with q=p, q=p+p/3 and q=p+4p/3 are the same point. Therefore, there are only two different corner points and they are known as the K and K points, where 4 p K =,0 (6.3) 3 3 a and K ' = - 4 p,0 3 3 a (6.3) For most of the momentum points, one of the two bands has positive energy e + > 0 and the other one has negative energy e - < 0. However at the corner of the, K and K, the two bands are degenerate e + =e - = 0. e K = H K = t 3 + cos 3 k x a + 4 cos 3 k x a Cos 3 k y a = t 3 + cos 3 4 p a + 4 cos 3 3 a 3 4 p a Cos a μ0 = (6.33) t 3 + cos 4 p cos p 3 = t 3 + μ μ - = 3-3 = 0 Two bands have the same energy at K and K, and they are in fact band crossing points.

15 Phys540.nb 9 Fig. 3. The energy dispersion for both bands (e as a function of k x and k y ) in the first (within the red dash lines mark in the figure above). Near K or K points, the energies of the two bands are linear functions of momentum e k - K. To see this, we expand e k near k ~ K and k ~ K ' e k =e K + q = t 3 + cos 3 4 p + q x a + 4 cos 3 3 a 3 4 p + q x a Cos a q y a = (6.34) 3 t a q x + q y + Oq = 3 t aq+ Oq Fig. 4. A zoom in to the Dirac point. Q: Fermions with linear dispersion e k. What are they? A: A Massless Dirac fermion. Near the K point, we have Hk = HK + q = 3 ta 0 q x -Âq y q x +Âq y 0 Near the K point, we have Hk = HK ' + q =- 3 ta 0 q x +Âq y q x -Âq y 0 + Oq º 3 taq x s x + q y s y = cq ÿs + Oq º- 3 taq x s x - q y s y =-cq ÿs x s s x (6.35) (6.36) Here c is the speed of light. Each of them is a Weyl fermion. These two Weyl fermions with opposite chirality forms a Dirac fermion. Dirac fermions we learned in quantum mechanics: Â t y= 0 cq ÿs cq ÿs 0 y (6.37) If one make an unitary transformation (changing the basis, which doesn t change any physics)

16 9 Phys540.nb y' = The equation turns into y  t y' = -cq ÿs x s s x 0 0 cq ÿs y' (6.39) The Dirac point (band crossing points) in a honeycomb lattice is stable as long as the space inversion symmetry r -r and the time reversal symmetry t -t are preserved. We will come back to this point latter. No matter how one perturb the systems (e.g. adding longer-range hoppings), the Dirac point is always there as long as the two symmetries mentioned above are preserved. For graphene, the lower band is filled and the upper band is empty, which is known as half-filling. The half here means that the number of electrons N e over the number of sites N s is N e N S =. However, it is worthwhile to notice that the ratio between N e and the number of unit cells N is, because there are two sites in each unit cell. So one of the two bands are totally filled. By gating, one can tune the Fermi energy slightly away from the Dirac point, which is known as a dopped graphene Aharonov Bohm effect and complex hopping Q: Can hopping strength t be complex? A: Yes, and the phase can come from the Aharonov Bohm effect in the presence of a magnetic field or spin-orbital couplings. The Aharonov Bohm effect: if one moves a charged particle around a closed contour. The phase difference between the final and initial states is proportional to the magnetic flux enclosed by the contour f=eñbÿ S = eñ A. l Now, let s consider a discrete lattice. Consider three sites a, b and c. The hopping strength between these three sites are t ab t bc and t ca respectively. If a particle hops from a to b and then to c, the hopping strength around this loop is: t ab t bc t ca = t ab Âfab * t bc Âfbc * t ca Âfca = t ab t bc t ca  fab+fbc+fca (6.40) The phase picked up by the electron is: f ab +f bc +f ca = e Ñ Bÿ S (6.4) If B is nonzero inside the triangle formed by these three sites, the phase for these hoppings are nonzero. Please notice that () t ab and t ba has opposite phase, due to the Hermitian condition (One can also prove this using the Fermi s golden rule, which says that the tunneling amplitude from the state f> to i> is the complex conjugate of the tuning amplitude from i> to f>). * t ab = t ba (6.4) () The individual phases for t ab, t bc and t ca have no physical meaning and their phases depend on the gauge choice. However, the total phase around a loop is independent of gauge, and in fact it is a physical observable, i.e. the magnetic flux. Proof: We know that the phase of the hopping term is f ab = eñ a baÿ l (6.43) Underage transformation: A A + c f ab f ab ' = eñ a b A + cÿ l = eñ a b Aÿ l + eñ a b cÿ l =f ab + c b -c a e Ñ (6.44) Obviously, f ab is not a physical observables, since it relies on the gauge choice. However, the total phase around a loop is different. It is a loop integral of A, which is gauge independent and the physics meaning of this Integral is the magnetic flux.

17 Phys540.nb 93 f ab +f bc +f ca = eñ Aÿ l f ab ' +f bc ' +f ca ' = eñ A + cÿ l = eñ Aÿ l + eñ cÿ l = eñ Aÿ l =f ab +f bc +f ca (6.45) The complex hopping strength induced by B fields breaks the time-reversal symmetry because B -B under time reversal. (in other words, under time-reversal one needs to flip the sign of all these phases) Complex next-nearest-neighbor (NNN) hoppings (breaking T-symmetry using B fields) Ref: F. D. M. Haldane, Model for a Quantum Hall Effect without Landau Levels: Condensed-Matter Realization of the Parity Anomaly, Phys. Rev. Lett. 6, (988). Now let us add some NNN hoppings and assume their hopping amplitudes are complex. For simplicity, we choose the amplitude and the complex phase to be the same for all NNN bonds. If the hopping is along the arrows marked in the figure, the phase of the hopping strength is f. If one hops in the opposite direction, the phase is -f. Fig. 5. The model of Haldane. This complex phases can be realized (in theory) by applying a staggered B field, which is positive near the center of each hexagon and negative near the edges. The NNN hoppings are from an a-site to another a-site (and from a b-site to another b-site). For a-to-a hoppings, there are three different types of NNN bonds, along q=0, p/3 and 4p/3. Same is true for b-to-b hoppings. So the Hamiltonian is: H NNN =-t ' Âf i a ri a ri +v - t ' Âf i a ri +v a ri -v 3 - t ' Âf i a ri -v 3 a ri + h.c.+a b and f-f (6.46) Here, v and v are marked on the first figure of this section with v = 3 a, 0 and v = - 3 a, 3 a. v = - 3 a, -3 a H NNN =-t ' Âf i a ri a ri+v - t ' Âf i a ri+v a ri-v 3 - t ' Âf i a ri-v 3 a ri + h.c.+a b and f-f =-t ' Âf k a k a k -Â kÿv + -Â kÿv + -Â kÿv 3 + h.c. +a b and f-f =- t ' k a k a k cosk ÿv -f + cosk ÿv -f + cosk ÿv 3 -f - (6.47) t ' k b k b k cosk ÿv +f + cosk ÿv +f + cosk ÿv 3 +f H = H NN + H NNN = k a k b k H k H k H k H k a k b k (6.48) H k =-t -Â k' + -Â k' ÿe + -Â k' ÿe ÿe 3 (6.49)

18 94 Phys540.nb H k = H k * =-t  k' ÿe +  k' +  k' H k =- t ' cosk ÿv -f + cosk ÿv -f + cosk ÿv 3 -f H k =- t ' cosk ÿv +f + cosk ÿv +f + cosk ÿv 3 +f ÿe ÿe 3 (6.50) (6.5) (6.5) If one computes the eigenvalues of H(k), one find that the two bands never cross with each other for any k (if t is non-zero and f is NOT an integer multiplied by p). Fig. 6. The energy dispersion for both bands (e as a function of k x and k y ) in the first (within the red dash lines mark in the figure above). Using Pauli matrices: H = H 0 k I + H x k s x + H y k s y + H z k s z H 0 = H k + H k =-t ' cos f cosk ÿv + cosk ÿv + cosk ÿv 3 H z = H k - H k =-t ' sin f sink ÿv + sink ÿv + sink ÿv 3 H x = ReH k =-tcosk ÿe + cosk ÿe + cosk ÿe3 (6.53) (6.54) (6.55) (6.56) H y = ImH k =-tsink ÿe + sink ÿe + sink ÿe3 (6.57) The energy dispersions: e k = H 0 k H x k + H y k + H z k (6.58) The gap between the two bands: e + k -e - k = H x k + H y k + H z k (6.59) At K or K, H x k = H y k = 0 and the gap is: e + k -e - k = H x k + H y k + H z k = H z k = K = 6 3 t ' sin f (6.60) In fact, at the K point, H z =-3 3 t ' sin f at the K point (6.6) H z = 3 3 t ' sin f (6.6) They have opposite signs (as long as f is not n p). Based on what we learned last time, this means that one cannot define the wavefunction in

19 Phys540.nb 95 the whole. We need to cut the systems into two regions. The region I contains the K point, the region II contains the K points. And we need to use different eigenvectors for these two regions. Using the same method we learned in the last lecture, one finds that the Chern number here is ±. Here, we can draw a small circle around the K point. Inside the circle, we use wavefunction u I. Outside it, we use u II. Âfk = H z+ Hk H x k+â H y k H z+ Hk = H x k +ÂH y k H x k +ÂH y k = q x +Âq y q x +Â q y = q = -Â Âq (6.63) H x k+â H y k Potential energy (breaking inversion symmetry) Let s keep NNN hoppings to be zero for now and add some potential energy to the Hamiltonian. H Potential = V + M i a i a i + V - M i b i b i = VN+ M i a i a i - M i b i b i (6.64) The V part (average potential between a and b sites) just adds a constant term to the energy, since the total particle number N is conserved. So we can drop the V term and only consider the difference between the potential energies at a and b sites (M). In k-space H Potential = M i a i a i - M i b i b i = M k a k a k - M k b k b k = k a k b k M s z a k b k (6.65) It adds a s z component to the Hamiltonian. Same as the NNN complex hopping, this term also opens a gap at the Dirac points. At K and K, H now has the same sign (H = M for any k). So, we can use one wave-function for the whole and thus the system is a trivial insulator. Q: What will happen if we have both H Potential and H NNN? A: We just need to look at the signs of H z at K and K points. If they have the same sign, we can use one wavefunction to cover the whole, so C = 0 (trivial insulator). If they have opposite signs, the system is a topological insulator with C=± At the K point, H z = M - 3 at the K point H z = M t ' sin f 3 t ' sin f (6.66) (6.67) Therefore, as long as M < 3 3 t ' sin f, the system is an topological insulator (H z flips sign). If M > 3 3 t ' sin f, H z is always positive (or negative) and thus the system is topologically trivial. The marginal case M = 3 3 t ' sin f is a topological transition, where H z = 0 at either the K point or the K point. Because H x = H y = 0 at these two points, the gap must be zero at one of the two points. Remarks: The model of Haldane is the first example of a topological insulator beyond quantum Hall effect. It demonstrates that topological insulator is a generic concept, which may appear in any insulating systems (NOT just quantum Hall). It also demonstrates that as long as the topological index is nonzero, one will observe all the topological phenomena expected for a quantum Hall state, including the quantized Hall conductivity and the existence of the edge states. The key differences between the model of Haldane and the quantum Hall effects are () the B field is on average zero in the model of Haldane while the QHE has a uniform B field and () there is a very strong lattice background in the model of Haldane while the QHE requires weak lattice potential. Systems similar to the Haldane s model are known as topological Chern insulators or Chern insulators (average B is 0 and have a strong lattice potential). But sometimes, Chern insulators are also used to refer to the quantum Hall effect. The model of Haldane is also the foundation to explore more complicated and exotic topological states. For example, the time-reversal invariant topological insulators was first proposed using a modified Haldane s models, which we will study latter in the semester.

20 96 Phys540.nb 6.4. Some symmetry properties of the Berry curvature F and the Chern number time-reversal transformation T and time-reversal symmetry In this section, we consider systems with time-reversal symmetry and shows that in the presence of the time-reversal symmetry, the Chern number must be zero. Therefore, to get a topological insulator with nonzero C, one must break the time-reversal symmetry. Dispersion relation Dispersion relation under time-reversal transformation: T e n k =e n -k (6.68) This is because T turns k into -k. There are two ways to see that T turns k into -k. k is the momentum. Under time-reversal, velocity changes sign, so does the momentum. k =-Â r. Under T, r remains invariant. However, because T is an anti-unitary transformation (which changes a complex number to its complex conjugate), it flips the sign of Â. Therefore, T k = T -Â r =Â r =-k. If the system is time-reversally invariant, the time-reversal symmetry implies that e n k is invariant under time-reversal transformation T e n k =e n k (6.69) Compare the two equations, we have: e n k =e n -k (6.70) Bottom line: for systems with time-reversal symmetry, the dispersion is an even function of the momentum k. Bloch wave function y n,k r = u n,k r kr Bloch wave under time-reversal transformation: T u n, k r = u n,-k r * (6.7) Here, T changes k to -k. In addition, because u n,k r is in general a complex function, T changes the function to its complex conjugate. If the system is time-reversally invariant, T u n, k r = Âfk u n, k r (6.7) By comparing the two equations above, we found that for a system with time-reversal symmetry, u n,-k r * = Âfk u n, k r (6.73) This arbitrary phase is here because wavefunctions are defined up to an arbitrary phase Berry connection A n k =- u n,k k u n,k The Berry connection is defined as: A n k =-Âu n,k k u n,k =-Â r u n,k r * k u n,k r (6.74) The under time-reversal, the new A n k is defined as TA n k =-ÂT u n,k k T u n,k = -Â r u n,-k r k u n,-k r * =Â r k u n,-k r u n,-k r * =-Â r u n,-k r * -k u n,-k r =-Âu n,-k -k u n,-k = A n -k (6.75) If the system has time-reversal symmetry TA n k = A n k + k ck (6.76) So we have A n -k = A n k + k ck (6.77) Bottom line: For systems with time-reversal symmetry, A n k and A n k differs by a gauge transformation k ck.

21 Phys540.nb 97 Berry curvature F n k =- e ij i u n,k j u n,k The Berry curvature is defined as: F n k =-Âe ij ki u n,k kj u n,k =-Âe ij r ki u n,k r * kj u n,k r (6.78) The under time-reversal, the new F n k is defined as TF n k =-Âe ij ki Tu n,k kj Tu n,k = -Â e ij r ki Tu n,k r * kj Tu n,k r =-Âe ij r ki u n,-k r kj u n,-k r * =-Âe ij r kj u n,-k r * ki u n,-k r (6.79) If we swap i and j here TF n k =-Âe ji r ki u n,-k r * kj u n,-k r =Âe ij r ki u n,-k r * kj u n,-k r =-F n -k (6.80) If the system has time-reversal symmetry TF n k = F n k (6.8) So we have F n -k =-F n k (6.8) Bottom line: For systems with time-reversal symmetry, F n k is an odd function of k Because the integral of an odd function over the whole Brillouin zone must be zero, (the contribution from k and -k cancels each other), the Chern number for a time-reversally invariant system must be C=0. The Chern number C = kf n k p TC = T p k x k y F n k =- p k x k y F n -k =- p k x k y F n k =-C (6.83) If the system has time-reversal symmetry, TC = C, (6.84) So we have C =-C, which means C = 0. Bottom line: to have a nontrivial Chern number, the system must break the time-reversal symmetry (using external B field or some other method) Space-inversion transformation I and space-inversion symmetry Space inversion transformation I turns a vector r into -r. It also changes k into -k, which is similar to T. However, I is an unitary transformation, while T is anti-unitary. Dispersion relation Dispersion relation under time-reversal transformation: Ie n k =e n -k (6.85) If the system is invariant under space inversion I e n k =e n k (6.86) Compare the two equations, we have: e n k =e n -k (6.87) Bottom line: For systems invariant under space inversion, the dispersion is an even function of the momentum k. If one want make a system where the dispersion is NOT an even function, one need to break both space inversion symmetry and timereversal symmetry.

22 98 Phys540.nb Bloch wave function y n,k r = u n,k r kr I u n, k r = u n,-k r (6.88) If the system is invariant under space inversion I u n, k r = Âfk u n, k r (6.89) By comparing the two equations above, we found that for a system with time-reversal symmetry, u n,-k r = Âfk u n, k r (6.90) Berry connection A n k =- u n,k k u n,k IA n k = A n -k (6.9) If the system is invariant under space inversion IA n k = A n k + k ck (6.9) So we have A n -k = A n k + k ck (6.93) Bottom line: For systems invariant under space inversion, A n k and A n k differs by a gauge transformation k ck. Berry curvature F n k =- e ij i u n,k IF n k = F n -k j u n,k (6.94) If the system is invariant under space inversion IF n k = F n k (6.95) So we have F n -k = F n k (6.96) Bottom line: For systems invariant under space inversion, F n k is an even function of k Time-reversal symmetry tell us F n k is an odd function of k. If system have both T and I symmetry, F must be both an even function and an odd function. So F=0 at any k point. The Chern number C = kf n k p IC = C (6.97) I T transformation and IT symmetry Under IT, ITF n k = I-F n -k =-F n k (6.98) If system is invariant under IT, ITF n k =-F n k (6.99) So, we have F n k =-F n k (6.00) Bottom line: For systems invariant under IT, F n k = 0 at every k. For T symmetry, the integral of F is zero (C=0), but F can be nonzero at each k point. For IT symmetry, F is zero at every single k point. This is a stronger statement! Please notice that all these conclusions assumes that the system is gapped. For gapless systems with IT symmetry, F can actually be nonzero and quantized, which will be discussed in the next section.

23 Phys540.nb p flux and Dirac points magnetic field, Aharonov Bohm effect and magnetic p- flux Q: Why B changes sign under T? A: Many ways to see it. Here, I use the A-B effect. In the presence of B field, Aharonov and Bohm tells us that if we move a charged-particle around a closed loop, the electron will pick up a phase, and the phase is the total magnetic field enclosed by the loop (times e/ñ) f=eñ Bÿ S = eñ A. l (6.0) Under time-reversal Âf turns into -Â f, because time-reversal is anti-unitary. So the A-B phase changes sign Tf=-f (6.0) In other words, the integral of B changes sign T D Bÿ S =- D Bÿ S (6.03) Because this equation is true for any region D, we have B =-B. If a system has time-reversal symmetry, everything must remain the same before and after we flip the arrow of time, including the A-B phase Âf = -Â f (6.04) Therefore, we have f=-f. So we have f=0 for any region D, which implies that B = 0. Q: Do we really need to have B = 0 to preserve the time-reversal symmetry. A: NO! This is because phase is only well-define up to mod p. The time-reversal symmetry requires Âf = -Â f, which does NOT implies f=-f. In fact, as long as f=-f+np, the A-B phase will be the same before and after the time-reversal transformation. This means that to preserve the time-reversal symmetry we just need f=np while n don t need to be 0. So we need to have D Bÿ S = n p for any D. This means that B can contains some delta functions. B = Ñ e i n i pdr - r i (6.05) For this B field, the integral of B in region D is just f=eñ Bÿ S = i ' ni p (6.06) Here the sum i ' is over all r i inside D. A delta function in B is a magnetic flux, which is known as a n p-flux. These fluxes don t break the time-reversal symmetry. In addition, it worthwhile to point out that for the A-B phase is only well defined up to mod p. So all the np fluxes give the same A-B phase, which is 0. All the (n+)p fluxes give the same A-B phase, which is p Berry flux We can do the same thing for the Berry curvature F(k). If we goes around a contour in the k-space, the phase change of the Bloch wave function around this contour is f= D F k (6.07) This is known as the Berry phase. Under IT, D is invariant so IT f=it D F k k =- D F k k =-f (6.08) If the system is IT invariant, we must have f=-f (but up to mod p).

24 00 Phys540.nb f=-f+ n p (6.09) so f=n p (6.0) As a result, if we have IT symmetry, D F k =np for any region D. So F must be either 0 or some delta functions. F k = i n i pdk - k i (6.) These delta functions are known as Berry fluxes. Because the Berry phase is only well defined up to mod p, we have in general two types of Berry fluxes: 0 and p. And a p flux is a Dirac point Dirac point Near a Dirac point, the kernel of the Hamiltonian is: 0 kx-âky H = v F kx +Âky 0 (6.) Using polar coordinates H = v F k 0 -Â f Âf 0 (6.3) The eigenvalues for this matrix is e = v F k. The eigenvectors are: u + = -Â f and u - = - -Â f (6.4) The Berry connection in the polar coordinates has two components, radius and angle. A +k =-Âu + k u + = 0 (6.5) A +f =- Â k u + f u + =- Â k - Âf f - -Â f =- Â k Âf -Â -Â f =- k (6.6) Let s choose a small circle around the Dirac point, the flux through this circle is: p p c=f k = Aÿ k = A f k f=- 0 0 k k f=- 0 p f=-p (6.7) Why is a Dirac point stable? This is because the Berry fluxes are very stable. If one has TI symmetry, the Berry flux is quantized to integer values. Once its value is fixed (say p), no perturbation can change it. The only thing that can change continuously is the location of the delta functions. In other words, without break the TI symmetry, one can only move the Dirac points around in the k-space, but they cannot just disappear. To get ride of Dirac points, one need to move two Dirac points together and let them annihilate each other (let a p flux meet a -p flux, so that we get zero flux). Another way to get ride of Dirac points is to break either the T or the I symmetry Edge states (numerical calculations) Why energy is a function of k x and k y Consider one electron moving in an infinite D lattice (for simplicity, we consider a square lattice here). Because the lattice is invariant under lattice translations x x + a and y y + a, lattice momentum is a good quantum number (a conserved quantity). Momentum conservation means that k x, H = k y, H = 0 (6.8) In addition, we know k x, k y = 0. With three operators which commute with each other, quantum mechanics tells us that we can find common eigenstates for all these three operators, and we can use these eigenstates as the basis for the Hilbert space. This set of basis looks like:

25 Phys540.nb 0 y en,k x,k y (6.9) This state is the eigenstate of H, k x and k y with eigenvalues e n, k x and k y respectively. Here, one need a sub-index (n) to distinguish different bands. So we have the band structure: e n k x, k y (6.0) What if we have an finite system? To study edge states, we must have an edge. So we cannot use an infinite D plan. Let s consider the simplest case here, a infinitely long stripe. We assumes that the system is infinitely long along x, but it has a finite width along y. Now, the translational symmetry along x is still preserved. But there is no translational symmetry along y (due to the existence of the edge). So k x is still a good quantum number but k y is NOT. k x, H = 0 but k y, H 0 (6.) So we cannot define common eigenstates for H, k x and k y, but we can still find eigenstates for H and k x, because they commute with each other. y em,k x (6.) Here, again, we need an extra index distinguish different states with the same k x. Here I emphasize that m is NOT just the band index n. It includes information for both n and k y. So now, energy is a bunch of function of k x e m k x (6.3) where m =,, 3 N S. We have N S functions of k x. We can draw these N S functions on the k x -e plane, which gives us N S curves, which is our D dispersion Examples We consider the model of Haldane for a system with infinite size along x, by finite width along y. The Hamiltonian is invariant under translations along the x axis x x + 3 a, so the x-component of the momentum is a conserved quantity k x, H = 0. So we can find common eigenstates of k x and H. The system has two edges: one on the top, one at the bottom. We can consider this lattice as a bunch of D (horizontal) lines coupled together, so the width are described by the number of chains one have in this infinite stripe N. For e m k x, we have m=,, N. This factor of comes from the fact that we have two sites in each unit cell for a honeycomb lattice. First, lets consider the case without NNN hopping (graphene, with two Dirac cones). Here, for the first panel, we show e m k x for a system with width N=0 so there are N S = N = 0 lines. The second panel is the same calculation for a system with width N=00 N S = N = 00. The last panel shows the bulk band structure (side view) Fig. 7. the energy spectrum for a infinite-long ribbon (the honeycomb lattice) The D figure of e m k x is in some sense the 3D band structure e n k x, k y projected onto a D plane. We can see here an upper band and a lower band and Dirac cones. The system is gapless. Now, let s consider the topological case, with complex NNN hoppings. The first panel shows the energy dispersion for an infinite stripe and the second one shows the bulk band structure for an infinite D system.

26 0 Phys540.nb Fig. 8. the energy spectrum for a infinite-long ribbon (the model of Haldane) Here, I choose the parameter f=p/, t NN = and t NNN = 0.. As one can see, for the infinite stripe, there are two bands, whose energy spectrum coincide with the bulk band structure for infinite systems. But in addition, there are two extra lines shown inside the gap, which are the metallic edge states. By looking at the wavefunctions (see the interactive figures on the course website), one find that one of the in gap state is localized near the top edge, while the other is localized near the bottom edge. If we consider the momentum region (0,p), instead of (-p,p), we find that one of the two edge states has positive slop and the other one has negative slop (left moving). Because the slop is the velocity of electrons (v = e k), one of the edge states has positive velocity (right moving) while the other one is right moving (negative velocity) Fig. 9. the velocity of the edge states