Wannier functions, Modern theory of polarization 1 / 51
|
|
- Sharleen Knight
- 5 years ago
- Views:
Transcription
1 Wannier functions, Modern theory of polarization 1 / 51
2 Literature: 1 R. D. King-Smith and David Vanderbilt, Phys. Rev. B 47, Nicola Marzari and David Vanderbilt, Phys. Rev. B 56, Raffaele Resta, Phys. Rev. Lett. 80, Raffaele Resta and David Vanderbilt, Theory of Polarization: A Modern Approach, in Physics of Ferroelectrics: a Modern Perspective, ed. by K.M. Rabe, C.H. Ahn, and J.-M. Triscone (Springer-Verlag, 2007, Berlin), pp (local preprint). 5 Jenő Sólyom. Fundamentals of the Physics of Solids: Volume II, Electronic properties (Springer Verlag, Berlin;Heiderberg, 2009). 6 Rui Yu, Xiao Liang Qi, Andrei Bernevig, Zhong Fang, and Xi Dai, Phys. Rev. B 84, (2011). 7 J. K. Asbóth, L. Oroszlány, and A. Pályi, arxiv: (A Short Course on Topological Insulators, Lecture Notes in Physics, Vol. 919, 2016, Springer Verlag). 2 / 51
3 Further reading: 1 E. Blount, Formalisms of band theory, Solid State Physics, 13, page New York, Academic Press, (1962). 2 S. Kivelson, Phys. Rev. B 26, 4269 (1982). 3 Raffaele Resta, What makes an insulator different from a metal?, arxiv preprint cond-mat/ , J. Zak, Phys. Rev. Lett. 62, 2747 (1989). 3 / 51
4 Overview Reminder: Bloch functions, periodic boundary condition Wannier functions: introduction Maximally localized Wannier functions Modern theory of polarization Polarization and Wannier centers Calculation of Berry phase in discrete k space Position operator and periodic boundary conditions Calculation of Wannier centers in 1D systems 4 / 51
5 Bloch functions, periodic boundary condition We want to describe the bulk properties of crystalline materials. In non-interacting approximation: N e ( ˆp 2 ) Ĥ = i +U(r i ) 2m e i=1 U(r) = U(r+R n ) is lattice-periodic potential, R n lattice vector. 5 / 51
6 Bloch functions, periodic boundary condition We want to describe the bulk properties of crystalline materials. In non-interacting approximation: N e ( ˆp 2 ) Ĥ = i +U(r i ) 2m e i=1 U(r) = U(r+R n ) is lattice-periodic potential, R n lattice vector. Hamiltonian is a sum of single-particle Hamiltonians, the solution of this N e electron problem can be obtained using the solutions of the Schrödinger equation ] [ 2 2 +U(r) Ψ i (r) = E i Ψ i (r) 2m e 5 / 51
7 Bloch functions, periodic boundary condition In order to capture the translational invariance of the bulk, we use use periodic boundary condition the solutions satisfy (Bloch theorem) Ψ k (r+r m ) = e ik Rm Ψ k (r) wavevector k can take on discrete values in the Brillouin zone 6 / 51
8 Bloch functions, periodic boundary condition In order to capture the translational invariance of the bulk, we use use periodic boundary condition the solutions satisfy (Bloch theorem) Ψ k (r+r m ) = e ik Rm Ψ k (r) wavevector k can take on discrete values in the Brillouin zone Equivalent formulation: Ψ k (r) = e ik r u k (r) where u k (r+r m ) = u k (r) lattice periodic and [ ] (ˆp + k) 2 +U(r) u n,k (r) = E n,k u n,k (r) 2m e 6 / 51
9 Bloch functions, periodic boundary condition For equivalent wavenumber vectors k = k+g (G lattice vector of the reciprocal lattice) Ψ n,k = e ik r u nk (r) = e i(k G) r u nk (r) = e ik r u nk (r) where u nk (r) = u nk+g (r) = e ig r u nk (r). One can also show that E nk+g = E nk. 7 / 51
10 Overview Reminder: Bloch functions, periodic boundary condition Wannier functions: introduction Maximally localized Wannier functions Modern theory of polarization Polarization and Wannier centers Calculation of Berry phase in discrete k space Position operator and periodic boundary conditions Calculation of Wannier centers 8 / 51
11 Wannier functions Since Ψ n,k+g (r) = e i(k+g) r u n,k+g = e ik r u n,k = Ψ n,k (r) for fixed r the Bloch function Ψ n,k (r) is periodic in the reciprocal space. 9 / 51
12 Wannier functions Since Ψ n,k+g (r) = e i(k+g) r u n,k+g = e ik r u n,k = Ψ n,k (r) for fixed r the Bloch function Ψ n,k (r) is periodic in the reciprocal space. it can be expanded into a Fourier series: Ψ n,k (r) = 1 N R j : Lattice vectors in real space R j Φ n (r,r j )e ik Rj 9 / 51
13 Wannier functions Since Ψ n,k+g (r) = e i(k+g) r u n,k+g = e ik r u n,k = Ψ n,k (r) for fixed r the Bloch function Ψ n,k (r) is periodic in the reciprocal space. it can be expanded into a Fourier series: Ψ n,k (r) = 1 N R j : Lattice vectors in real space Reverse transformation: Φ n (r,r j ) = 1 N R j Φ n (r,r j )e ik Rj k BZ e ik R j Ψ n,k (r) 9 / 51
14 Wannier functions Φ n (r,r j ) is a function of r R j : Φ n (r+r n,r j +R n ) = = = 1 N k Bz 1 N k BZ 1 N k BZ = Φ n (r,r j ) e ik (R j+r n) Ψ nk (r+r n ) e ik (R j+r n) e ik Rn Ψ nk (r) e ik R j Ψ nk (r) 10 / 51
15 Wannier functions Φ n (r,r j ) is a function of r R j : Φ n (r+r n,r j +R n ) = = = 1 N k Bz 1 N k BZ 1 N k BZ = Φ n (r,r j ) e ik (R j+r n) Ψ nk (r+r n ) e ik (R j+r n) e ik Rn Ψ nk (r) e ik R j Ψ nk (r) For R n = R j : Φ n (r,r j ) = Φ n (r R j,0) it depends only on r R j 10 / 51
16 Wannier functions Φ n (r,r j ) is a function of r R j : Φ n (r+r n,r j +R n ) = = = 1 N k Bz 1 N k BZ 1 N k BZ = Φ n (r,r j ) e ik (R j+r n) Ψ nk (r+r n ) e ik (R j+r n) e ik Rn Ψ nk (r) e ik R j Ψ nk (r) For R n = R j : Φ n (r,r j ) = Φ n (r R j,0) it depends only on r R j (See, e.g., Ref[5]). Φ n (r R j,0) = W n (r R j ) are Wannier functions 10 / 51
17 Wannier functions Orthonormality and completeness: Bloch functions: Ψ n,k (r)ψ n,k (r)dr = δ n,n δ k,k Ψ n,k (r)ψ n,k(r ) = δ(r r ) n,k Wannnier functions W n(r R j )W n (r R j )dr = δ Rj,R j δ n,n (See, e.g., Ref[5]). n,r j W n (r R j)w n (r R j ) = δ(r r ) 11 / 51
18 Wannier centers Define the Wannier center of W n (r R j ) as r (j) n = W n (r R j ) r W n (r R j ) In general, a non-trivial question if the the above expectation value is finite in extended systems 12 / 51
19 Wannier centers Define the Wannier center of W n (r R j ) as r (j) n = W n (r R j ) r W n (r R j ) In general, a non-trivial question if the the above expectation value is finite in extended systems localization properties of Wannier functions 12 / 51
20 Wannier centers Define the Wannier center of W n (r R j ) as r (j) n = W n (r R j ) r W n (r R j ) In general, a non-trivial question if the the above expectation value is finite in extended systems localization properties of Wannier functions For explicit calculation of the Wannier functions in the SSH model, see Phys. Rev. B 26, 4269 (2016). 12 / 51
21 Wannier centers in the Rice-Mele model Consider the Rice-Mele model (N site long, periodic boundary condition) Bloch functions: Ψ n (k) = k u n (k) 13 / 51
22 Wannier centers in the Rice-Mele model Consider the Rice-Mele model (N site long, periodic boundary condition) Bloch functions: Ψ n (k) = k u n (k) n = 1,2 (conduction and valence band), u n (k): eigenstates of the momentum-space Hamiltonian H(k) 13 / 51
23 Wannier centers in the Rice-Mele model Consider the Rice-Mele model (N site long, periodic boundary condition) Bloch functions: Ψ n (k) = k u n (k) n = 1,2 (conduction and valence band), u n (k): eigenstates of the momentum-space Hamiltonian H(k) k = 1 N e ikxl, k x { 2π N N,22π N,...,N2π N } l=1 13 / 51
24 Wannier centers in the Rice-Mele model Consider the Rice-Mele model (N site long, periodic boundary condition) Bloch functions: Ψ n (k) = k u n (k) n = 1,2 (conduction and valence band), u n (k): eigenstates of the momentum-space Hamiltonian H(k) Position operator: k = 1 N e ikxl, k x { 2π N N,22π N,...,N2π N } ˆx = l=1 N m( m,a m,a + m,b m,b ) m=1 m: unit cell index; A, B, site index in a unit cell 13 / 51
25 Wannier centers in the Rice-Mele model ˆxW n (j) = ˆx 1 N k x e ijkx Ψ n (k x ) 1 = e ijkxˆx 1 N e ilkx l u n (k x ) N k x N l = 1 me imkx m u n (k x ) N k x e ijkx m 14 / 51
26 Wannier centers in the Rice-Mele model For N k x continous, use partial integration: ˆxW n (j) = 1 2π 2π 0 = 1 2π ( i) + 1 2π + i 2π 2π 0 2π 0 dk x e ijkx me imkx m u n (k x ) 2π 0 m dk x e i(m j)kx m u n (k x ) k x m dk x j e ijkx m u n (k x ) dk x m m m k x u n (k x ) 15 / 51
27 Wannier centers in the Rice-Mele model This can be simplified: ˆxW n (j) = ( i 2π + j 2π + i 2π e i(m j)kx m u n (k x ) ) m 2π 0 2π 0 dk x e ijkx m u n (k x ) m dk x m kx u n (k x ) m 2π 0 16 / 51
28 Wannier centers in the Rice-Mele model ˆxW n (j) = ( i 2π + j 2π + i 2π e i(m j)kx m u n (k x ) ) m 2π 0 2π 0 dk x e ijkx m u n (k x ) = j W n (j) m dk x m kx u n (k x ) m 2π 0 0 (periodic function) 17 / 51
29 Wannier centers in the Rice-Mele model ˆxW n (j) = One finds: ( i 2π + j 2π + i 2π e i(m j)kx m u n (k x ) ) m 2π 0 2π 0 dk x e ijkx m u n (k x ) = j W n (j) m dk x m kx u n (k x ) m 2π 0 0 (periodic function) W n (j) ˆx W n (j) = i 2π dk x u n (k x ) kx u n (k x ) +j 2π 0 17 / 51
30 Wannier centers in the Rice-Mele model ˆxW n (j) = One finds: ( i 2π + j 2π + i 2π e i(m j)kx m u n (k x ) ) m 2π 0 2π 0 dk x e ijkx m u n (k x ) = j W n (j) m dk x m kx u n (k x ) m 2π 0 0 (periodic function) W n (j) ˆx W n (j) = i 2π dk x u n (k x ) kx u n (k x ) +j 2π 0 can be expressed in terms of Berry-phase (Zak s phase) 17 / 51
31 Wannier centers in the Rice-Mele model ˆxW n (j) = One finds: ( i 2π + j 2π + i 2π e i(m j)kx m u n (k x ) ) m 2π 0 2π 0 dk x e ijkx m u n (k x ) = j W n (j) m dk x m kx u n (k x ) m 2π 0 0 (periodic function) W n (j) ˆx W n (j) = i 2π dk x u n (k x ) kx u n (k x ) +j 2π 0 can be expressed in terms of Berry-phase (Zak s phase) Wannier centers are equally spaced 17 / 51
32 Overview Reminder: Bloch functions, periodic boundary condition Wannier functions: introduction Maximally localized Wannier functions Modern theory of polarization Polarization and Wannier centers Calculation of Berry phase in discrete k space Position operator and periodic boundary conditions Calculation of Wannier centers in 1D systems 18 / 51
33 Maximally localized Warnier functions Figure: Phys Rev B 56, (1997). 19 / 51
34 Maximally localized Warnier functions Gauge freedom in choosing the Bloch orbitals: Ψ nk (r) e iφn(k) Ψ nk (r) [u nk (r) e iφn(k) u nk (r)] describes the same electron density 20 / 51
35 Maximally localized Warnier functions Gauge freedom in choosing the Bloch orbitals: Ψ nk (r) e iφn(k) Ψ nk (r) [u nk (r) e iφn(k) u nk (r)] describes the same electron density Wannier functions are not unique! 20 / 51
36 Maximally localized Warnier functions Gauge freedom in choosing the Bloch orbitals: Ψ nk (r) e iφn(k) Ψ nk (r) [u nk (r) e iφn(k) u nk (r)] describes the same electron density Wannier functions are not unique! If there are band degeneracies for some k (symmetries, e.g., time reversal): not sufficient to consider isolated bands consider a set of J bands (composite bands) 20 / 51
37 Maximally localized Warnier functions Gauge freedom in choosing the Bloch orbitals: Ψ nk (r) e iφn(k) Ψ nk (r) [u nk (r) e iφn(k) u nk (r)] describes the same electron density Wannier functions are not unique! If there are band degeneracies for some k (symmetries, e.g., time reversal): not sufficient to consider isolated bands consider a set of J bands (composite bands) A more general U(J) gauge freedom: u nk (r) p U (k) pn u pk (r) 20 / 51
38 Maximally localized Warnier functions (MLWF) Define for the Wannier functions at the origin R j = 0 the following: 1 Wannier center r n = W n (r) r W n (r) 2 second moment r 2 n = W n (r) r 2 W n (r) 21 / 51
39 Maximally localized Warnier functions (MLWF) Define for the Wannier functions at the origin R j = 0 the following: 1 Wannier center r n = W n (r) r W n (r) 2 second moment r 2 n = W n (r) r 2 W n (r) A measure of the spread (delocalization) of the Wannier functions: Ω = n [ r 2 n r 2 n] One can then try to minimize Ω with respect to the unitary transformations U (k) pn Maximally localized Wannier functions (MLWF) 21 / 51
40 Maximally localized Warnier functions (MLWF) Define for the Wannier functions at the origin R j = 0 the following: 1 Wannier center r n = W n (r) r W n (r) 2 second moment r 2 n = W n (r) r 2 W n (r) A measure of the spread (delocalization) of the Wannier functions: Ω = n [ r 2 n r 2 n] One can then try to minimize Ω with respect to the unitary transformations U (k) pn Maximally localized Wannier functions (MLWF) This approach can be used to obtain MLWFs from DFT calculations in plane-wave basis. 21 / 51
41 Maximally localized Warnier functions (MLWF) One can decompose Ω = Ω 1 +Ω 2, where Ω 1 = r 2 n Rm r 0n 2 n R,m Ω 2 = n Rm 0n Rm r 0n 2 22 / 51
42 Maximally localized Warnier functions (MLWF) One can decompose Ω = Ω 1 +Ω 2, where Ω 1 = r 2 n Rm r 0n 2 n R,m One can show that Ω 2 = n Rm 0n Rm r 0n 2 Ω 1 is gauge invariant, i.e., does not depend on U (k) pn transformation Ω 2 is positive definite 22 / 51
43 Maximally localized Warnier functions (MLWF) One can decompose Ω = Ω 1 +Ω 2, where Ω 1 = r 2 n Rm r 0n 2 n R,m One can show that Ω 2 = n Rm 0n Rm r 0n 2 Ω 1 is gauge invariant, i.e., does not depend on U (k) pn transformation Ω 2 is positive definite minimalization of Ω means minimalization of Ω 2 22 / 51
44 Special case: MLWF in 1D Notation: W n (x R x ) = R x,n 23 / 51
45 Special case: MLWF in 1D Notation: W n (x R x ) = R x,n Consider a group of bands, e.g., the valence band(s): Projectors: P = R x,n R x,n = Ψ n,kx Ψ n,kx R x,n n,k x 23 / 51
46 Special case: MLWF in 1D Notation: W n (x R x ) = R x,n Consider a group of bands, e.g., the valence band(s): Projectors: P = R x,n R x,n = Ψ n,kx Ψ n,kx R x,n n,k x In 1D the eigenfunctions of the projected position operator P ˆx P are MLWF. 23 / 51
47 Special case: MLWF in 1D Notation: W n (x R x ) = R x,n Consider a group of bands, e.g., the valence band(s): Projectors: P = R x,n R x,n = Ψ n,kx Ψ n,kx R x,n n,k x In 1D the eigenfunctions of the projected position operator P ˆx P are MLWF. Let 0m be an eigenfunction of P ˆx P with eigenvalue x 0m. R n ˆx 0m = R n P ˆx P 0m = x 0m δ R,0 δ m,n 23 / 51
48 Special case: MLWF in 1D Notation: W n (x R x ) = R x,n Consider a group of bands, e.g., the valence band(s): Projectors: P = R x,n R x,n = Ψ n,kx Ψ n,kx R x,n n,k x In 1D the eigenfunctions of the projected position operator P ˆx P are MLWF. Let 0m be an eigenfunction of P ˆx P with eigenvalue x 0m. R n ˆx 0m = R n P ˆx P 0m = x 0m δ R,0 δ m,n Ω 2 vanishes, Ω 1 gauge invariant 23 / 51
49 Special case: MLWF in 1D Notation: W n (x R x ) = R x,n Consider a group of bands, e.g., the valence band(s): Projectors: P = R x,n R x,n = Ψ n,kx Ψ n,kx R x,n n,k x In 1D the eigenfunctions of the projected position operator P ˆx P are MLWF. Let 0m be an eigenfunction of P ˆx P with eigenvalue x 0m. R n ˆx 0m = R n P ˆx P 0m = x 0m δ R,0 δ m,n Ω 2 vanishes, Ω 1 gauge invariant 0m is a MLWF 23 / 51
50 Special case: MLWF in 1D Notation: W n (x R x ) = R x,n Consider a group of bands, e.g., the valence band(s): Projectors: P = R x,n R x,n = Ψ n,kx Ψ n,kx R x,n n,k x In 1D the eigenfunctions of the projected position operator P ˆx P are MLWF. Let 0m be an eigenfunction of P ˆx P with eigenvalue x 0m. R n ˆx 0m = R n P ˆx P 0m = x 0m δ R,0 δ m,n Ω 2 vanishes, Ω 1 gauge invariant 0m is a MLWF Argument does not work in 3D, because P ˆx P, Pŷ P, PẑP do not commmute 23 / 51
51 Summary I bulk properties of materials periodic boundary conditions Bloch functions alternative description: Wannier functions Wannier-centers (expectation value of the position operator): in Rice-Mele model Berry-phase Maximally localized Wannier functions: in 1D eigenfunctions of the projected position operator P ˆxP 24 / 51
52 Overview Reminder: Bloch functions, periodic boundary condition Wannier functions: introduction Maximally localized Wannier functions Modern theory of polarization Polarization and Wannier centers Calculation of Berry phase in discrete k space Position operator and periodic boundary conditions Calculation of Wannier centers in 1D systems 25 / 51
53 Modern theory of polarization Figure: Theory of Polarization: A Modern Approach, in Physics of Ferroelectrics: a Modern Perspective, ed. by K.M. Rabe, C.H. Ahn, and J.-M. Triscone (Springer-Verlag, 2007, Berlin), pp (local preprint) 26 / 51
54 Modern theory of polarization Macroscopic polarization P: fundamental concept in the phenomenological description of dielectrics 27 / 51
55 Modern theory of polarization Macroscopic polarization P: fundamental concept in the phenomenological description of dielectrics Some materials possess polarization without external field (ferroelectricity) or become polarized upon applying strain (pieozoelectricity) 27 / 51
56 Modern theory of polarization Macroscopic polarization P: fundamental concept in the phenomenological description of dielectrics Some materials possess polarization without external field (ferroelectricity) or become polarized upon applying strain (pieozoelectricity) How to measure the polarization? Figure: Pieozoelectricity: surface or bulk effect? (Figure from Ref[4]). 27 / 51
57 Modern theory of polarization Macroscopic polarization P: fundamental concept in the phenomenological description of dielectrics Some materials possess polarization without external field (ferroelectricity) or become polarized upon applying strain (pieozoelectricity) How to measure the polarization? Figure: Pieozoelectricity: surface or bulk effect? (Figure from Ref[4]). While the crystal is strained, a transient electrical current flows through the sample 27 / 51
58 Modern theory of polarization Fundamental relation: the change in polarization P is accompanied with a transient current j(t) flowing through the sample: dp dt = j(t) 28 / 51
59 Modern theory of polarization Fundamental relation: the change in polarization P is accompanied with a transient current j(t) flowing through the sample: Change in polarization: dp dt = j(t) P = P( t) P(0) = If the change is slow enough adiabatic limit t 0 j(t) 28 / 51
60 Modern theory of polarization Fundamental relation: the change in polarization P is accompanied with a transient current j(t) flowing through the sample: Change in polarization: dp dt = j(t) P = P( t) P(0) = t If the change is slow enough adiabatic limit adiabatic perturbation theory to calculate j(t) 0 j(t) 28 / 51
61 Formal description of the theory: polarization P Introduce the parameter λ [0, 1]: dimensionless adiabatic time 29 / 51
62 Formal description of the theory: polarization P Introduce the parameter λ [0, 1]: dimensionless adiabatic time P = λ = 0:initial state; λ = 1:final state 1 0 dλ dp dλ 29 / 51
63 Formal description of the theory: polarization P Introduce the parameter λ [0, 1]: dimensionless adiabatic time P = λ = 0:initial state; λ = 1:final state Assumptions: 1 0 dλ dp dλ system remains insulating for all values of λ system bulk retains crystalline periodicity for all λ 29 / 51
64 Formal description of the theory: current j (λ) 30 / 51
65 Formal description of the theory: current j (λ) Because of crystalline periodicity the solutions of the Schrödinger equation are Bloch-functions: Ψ (λ) k (r) = eik r u λ k (r) 30 / 51
66 Formal description of the theory: current j (λ) Because of crystalline periodicity the solutions of the Schrödinger equation are Bloch-functions: Ψ (λ) k (r) = eik r uk λ(r) Schrödinger equation for the lattice periodic part u (λ) k (r) [ ] (ˆp + k) 2 +U (λ) (r) u (λ) 2m n,k (r) = E(λ) n,k u(λ) n,k (r) e 30 / 51
67 Formal description of the theory: current j (λ) Because of crystalline periodicity the solutions of the Schrödinger equation are Bloch-functions: Ψ (λ) k (r) = eik r uk λ(r) Schrödinger equation for the lattice periodic part u (λ) k (r) [ ] (ˆp + k) 2 +U (λ) (r) u (λ) 2m n,k (r) = E(λ) n,k u(λ) n,k (r) e Adiabatic change: the current j (λ) can be calculated using adiabatic perturbation theory 30 / 51
68 Formal description of the theory: current j (λ) Because of crystalline periodicity the solutions of the Schrödinger equation are Bloch-functions: Ψ (λ) k (r) = eik r uk λ(r) Schrödinger equation for the lattice periodic part u (λ) k (r) [ ] (ˆp + k) 2 +U (λ) (r) u (λ) 2m n,k (r) = E(λ) n,k u(λ) n,k (r) e Adiabatic change: the current j (λ) can be calculated using adiabatic perturbation theory For a adiabatically changing time periodic 1D lattice [H(k,t) = H(k,t +T)] this is done in Chapter 5 of the Lecture Notes (Ref[7]) The same steps can be done here: t λ, λ does not need be periodic 30 / 51
69 Formal description of the theory: current j (λ) Current from a single filled band n: dp n (λ) dλ = j(λ) n = ie (2π) 3 dk[ k u (λ) n,k λu (λ) n,k λu (λ) n,k ku (λ) n,k ] BZ 31 / 51
70 Formal description of the theory: current j (λ) Current from a single filled band n: dp n (λ) dλ = j(λ) n = ie (2π) 3 dk[ k u (λ) n,k λu (λ) n,k λu (λ) n,k ku (λ) n,k ] BZ Total change in polarization ( number of pumped particles): P = ie (2π) 3 N n=1 BZ N: number of filled bands dk λ 0 dλ[ k u (λ) n,k λu (λ) n,k λu (λ) n,k ku (λ) n,k ] 31 / 51
71 Polarization as Berry phase For simplicity, consider a 1D system: P = ie (2π) N n=1 BZ dk λ 0 dλ[ k u (λ) n,k λu (λ) n,k λu (λ) n,k ku (λ) n,k ] 32 / 51
72 Polarization as Berry phase For simplicity, consider a 1D system: P = ie (2π) N n=1 BZ dk λ 0 dλ[ k u (λ) n,k λu (λ) n,k λu (λ) n,k ku (λ) n,k ] Make use of the following to re-write the integral for P partial integration with respect to λ k u (λ) n,k λu (λ) n,k = ku (λ) n,k λu (λ) n,k + u(λ) n,k k λ u (λ) n,k 32 / 51
73 Polarization as Berry phase For simplicity, consider a 1D system: P = ie (2π) N n=1 BZ dk λ 0 dλ[ k u (λ) n,k λu (λ) n,k λu (λ) n,k ku (λ) n,k ] Make use of the following to re-write the integral for P partial integration with respect to λ k u (λ) n,k λu (λ) n,k = ku (λ) n,k λu (λ) n,k + u(λ) n,k k λ u (λ) n,k P = ie (2π) N n=1 BZ dk [ u (λ) n,k ku (λ) n,k dλ ] k u(λ) n,k λu (λ) n,k 32 / 51
74 Polarization as Berry phase Note, that 1 0 dλ dk BZ k u(λ) n,k λu (λ) n,k = 0 33 / 51
75 Polarization as Berry phase Note, that 1 dλ dk 0 BZ k u(λ) n,k λu (λ) n,k = 0 because the Bloch functions are periodic in reciprocal space: Ψ (λ) k,n (x) = Ψ(λ) k+g,n (x) 33 / 51
76 Polarization as Berry phase Note, that 1 dλ dk 0 BZ k u(λ) n,k λu (λ) n,k = 0 because the Bloch functions are periodic in reciprocal space: Ψ (λ) k,n (x) = Ψ(λ) k+g,n (x) u(λ) n,k (x) = eigx u (λ) n,k+g (x) 33 / 51
77 Polarization as Berry phase Note, that 1 dλ dk 0 BZ k u(λ) n,k λu (λ) n,k = 0 because the Bloch functions are periodic in reciprocal space: Ψ (λ) k,n (x) = Ψ(λ) k+g,n (x) u(λ) n,k (x) = eigx u (λ) n,k+g (x) u (λ) n,k λu (λ) n,k is periodic in k its integral over the BZ is zero. 33 / 51
78 Polarization as Berry phase Note, that 1 dλ dk 0 BZ k u(λ) n,k λu (λ) n,k = 0 because the Bloch functions are periodic in reciprocal space: Ψ (λ) k,n (x) = Ψ(λ) k+g,n (x) u(λ) n,k (x) = eigx u (λ) n,k+g (x) u (λ) n,k λu (λ) n,k is periodic in k its integral over the BZ is zero. One finds: where P (λ) = ie (2π) P = P (λ=1) P (λ=0) N n=1 BZ dk u (λ) n,k ku (λ) n,k 33 / 51
79 Polarization as Berry phase P (λ) = ie (2π) N n=1 BZ dk u (λ) n,k k u (λ) n,k the right hand side is the integral of the Berry-connection over the BZ Berry phase (Zak s phase) 34 / 51
80 Polarization as Berry phase P (λ) = ie (2π) N n=1 BZ dk u (λ) n,k k u (λ) n,k the right hand side is the integral of the Berry-connection over the BZ Berry phase (Zak s phase) in case the path in λ space takes the crystal from its centrosymmetric reference state (P (λ=0) = 0) to its equilibrium polarized state spontaneous polarization 34 / 51
81 Polarization as Berry phase P (λ) = ie (2π) N n=1 BZ dk u (λ) n,k k u (λ) n,k the right hand side is the integral of the Berry-connection over the BZ Berry phase (Zak s phase) in case the path in λ space takes the crystal from its centrosymmetric reference state (P (λ=0) = 0) to its equilibrium polarized state spontaneous polarization Note: the total electrical polarization of any material has an ionic contribution as well (but this will not be important in the discussion of topological properties) 34 / 51
82 Polarization as Berry phase: gauge dependence There is a gauge freedom in defining u n,k : ũ n,k = e iβ(k) u n,k gives the same electron density. 35 / 51
83 Polarization as Berry phase: gauge dependence There is a gauge freedom in defining u n,k : ũ n,k = e iβ(k) u n,k gives the same electron density. β(k) is not arbitrary: Bloch functions obey Ψ n,k+g (x) = Ψ n,k (x) β(2π/a) β(0) = 2πj, j integer 35 / 51
84 Polarization as Berry phase: gauge dependence There is a gauge freedom in defining u n,k : ũ n,k = e iβ(k) u n,k gives the same electron density. β(k) is not arbitrary: Bloch functions obey Ψ n,k+g (x) = Ψ n,k (x) β(2π/a) β(0) = 2πj, j integer the polarization is well defined only modulo 2π: dk ũ (λ) n,k k ũ (λ) n,k = dk u (λ) n,k k u (λ) n,k +2iπj BZ BZ 35 / 51
85 Polarization as Berry phase: gauge dependence There is a gauge freedom in defining u n,k : ũ n,k = e iβ(k) u n,k gives the same electron density. β(k) is not arbitrary: Bloch functions obey Ψ n,k+g (x) = Ψ n,k (x) β(2π/a) β(0) = 2πj, j integer the polarization is well defined only modulo 2π: dk ũ (λ) n,k k ũ (λ) n,k = dk u (λ) n,k k u (λ) n,k +2iπj BZ BZ In 3D: P n = P n + er V cell, R a lattice vector. 35 / 51
86 Overview Reminder: Bloch functions, periodic boundary condition Wannier functions: introduction Maximally localized Wannier functions Modern theory of polarization Polarization and Wannier centers Calculation of Berry phase in discrete k space Position operator and periodic boundary conditions Calculation of Wannier centers in 1D systems 36 / 51
87 Polarization and Wannier centers Remember the relation between the Bloch functions and Wannier functions M: number of unit cells Ψ n,k (x) = 1 M M m=1 W n (x R m )e ik Rm 37 / 51
88 Polarization and Wannier centers Remember the relation between the Bloch functions and Wannier functions M: number of unit cells Ψ n,k (x) = 1 M M m=1 u n,k (x) = 1 M M m=1 W n (x R m )e ik Rm W n (x R m )e ik (R m x) 37 / 51
89 Polarization and Wannier centers Remember the relation between the Bloch functions and Wannier functions M: number of unit cells Ψ n,k (x) = 1 M M m=1 u n,k (x) = 1 M M m=1 W n (x R m )e ik Rm W n (x R m )e ik (R m x) Using the Berry-phase formula, one can express the polarization in terms of the Wannier functions: P (λ) = e N W n (λ) (x) x W n (λ) (x) n=1 37 / 51
90 Polarization and Wannier centers Remember the relation between the Bloch functions and Wannier functions M: number of unit cells Ψ n,k (x) = 1 M M m=1 u n,k (x) = 1 M M m=1 W n (x R m )e ik Rm W n (x R m )e ik (R m x) Using the Berry-phase formula, one can express the polarization in terms of the Wannier functions: P (λ) = e N W n (λ) (x) x W n (λ) (x) n=1 Polarization sum of the Wannier centers of the occupied bands (for one given R m ) 37 / 51
91 Tracking the Wannier centers for cyclic λ Until now, λ is not periodic 38 / 51
92 Tracking the Wannier centers for cyclic λ Until now, λ is not periodic The formalism remains valid for cyclic changes: H λ=1 = H λ=0 P cyc = 1 0 dλ dp dλ 38 / 51
93 Tracking the Wannier centers for cyclic λ Until now, λ is not periodic The formalism remains valid for cyclic changes: H λ=1 = H λ=0 P cyc = 1 0 dλ dp dλ The Wannier center must return to their initial location at the end of the cyclic evolution. But this is possible in two different ways: Figure: Fig.10 of Theory of Polarization: A Modern Approach 38 / 51
94 Tracking the Wannier centers in the Rice-Mele model In the case of Rice-Mele model, λ = t time, the adiabatic charge pumping can be visualized by tracking the Wannier centers: 39 / 51
95 Tracking the Wannier centers in the Rice-Mele model In the case of Rice-Mele model, λ = t time, the adiabatic charge pumping can be visualized by tracking the Wannier centers: W n (t) (j) ˆx W n (t) (j) = i 2π dk x u n (t) (k x ) kx u n (t) (k x ) +j 2π 0 Figure: Figure 4.5(a) of the Lecture Notes. Time evolution of the Wannnier centers of the bands. Solid line: valence band, dashed line: conduction band. The Chern number is / 51
96 Overview Reminder: Bloch functions, periodic boundary condition Wannier functions: introduction Maximally localized Wannier functions Modern theory of polarization Polarization and Wannier centers Calculation of Berry phase in discrete k space Position operator and periodic boundary conditions Calculation of Wannier centers in 1D systems 40 / 51
97 Polarization as Berry phase: calculation in discrete k-space We only consider the 1D case P n (λ) = e 2π Φ(λ) n ; L length of sample Φ n (λ) = Im BZ dk u (λ) n,k k u (λ) n,k 41 / 51
98 Polarization as Berry phase: calculation in discrete k-space We only consider the 1D case P n (λ) = e 2π Φ(λ) n ; L length of sample Φ n (λ) = Im BZ dk u (λ) n,k k u (λ) n,k Typically, we need to calculate this on a discrete grid of k-points dk u n,k k u n,k dk k j dk u n,k k u n,k k=kj 41 / 51
99 Polarization as Berry phase: calculation in discrete k-space We only consider the 1D case P n (λ) = e 2π Φ(λ) n ; L length of sample Φ n (λ) = Im BZ dk u (λ) n,k k u (λ) n,k Typically, we need to calculate this on a discrete grid of k-points dk u n,k k u n,k dk k j dk u n,k k u n,k k=kj Note the following: u n,k+dk u n,k + k u n,k dk 41 / 51
100 Polarization as Berry phase: calculation in discrete k-space We only consider the 1D case P n (λ) = e 2π Φ(λ) n ; L length of sample Φ n (λ) = Im BZ dk u (λ) n,k k u (λ) n,k Typically, we need to calculate this on a discrete grid of k-points dk u n,k k u n,k dk k j dk u n,k k u n,k k=kj Note the following: u n,k+dk u n,k + k u n,k dk u n,k u n,k+dk 1+ u n,k k u n,k dk 41 / 51
101 Polarization as Berry phase: calculation in discrete k-space We only consider the 1D case P n (λ) = e 2π Φ(λ) n ; L length of sample Φ n (λ) = Im BZ dk u (λ) n,k k u (λ) n,k Typically, we need to calculate this on a discrete grid of k-points dk u n,k k u n,k dk k j dk u n,k k u n,k k=kj Note the following: u n,k+dk u n,k + k u n,k dk u n,k u n,k+dk 1+ u n,k k u n,k dk ln[ u n,k u n,k+dk ] = ln[1+ u n,k k u n,k dk] u n,k k u n,k dk 41 / 51
102 Polarization as Berry phase: calculation in discrete k-space dk u n,k k u n,k dk k j ln[ u n,k u n,k+dk ] k=kj 42 / 51
103 Polarization as Berry phase: calculation in discrete k-space dk u n,k k u n,k dk k j ln[ u n,k u n,k+dk ] k=kj Take k j = 2πj J a, a lattice constant, J, j = 0,...J 1 integer, dk = 2π/(J a) 42 / 51
104 Polarization as Berry phase: calculation in discrete k-space dk u n,k k u n,k dk k j ln[ u n,k u n,k+dk ] k=kj Take k j = 2πj J a, a lattice constant, J, j = 0,...J 1 integer, dk = 2π/(J a) Φ (λ) n = Im j ln[ u (λ) n,k j u (λ) J 1 n,k j+1 ] = Imln u (λ) n,k j u (λ) n,k j+1 j=0 42 / 51
105 Polarization as Berry phase: calculation in discrete k-space dk u n,k k u n,k dk k j ln[ u n,k u n,k+dk ] k=kj Take k j = 2πj J a, a lattice constant, J, j = 0,...J 1 integer, dk = 2π/(J a) Φ (λ) n = Im j ln[ u (λ) n,k j u (λ) J 1 n,k j+1 ] = Imln u (λ) n,k j u (λ) n,k j+1 j=0 Note, the product u (λ) n,k 0 u (λ) n,k 1 u (λ) n,k 1 u (λ) n,k 2... does not depend on the phase of u (λ) n,k j s gauge independent 42 / 51
106 Polarization as Berry phase: calculation in discrete k-space dk u n,k k u n,k dk k j ln[ u n,k u n,k+dk ] k=kj Take k j = 2πj J a, a lattice constant, J, j = 0,...J 1 integer, dk = 2π/(J a) Φ (λ) n = Im j ln[ u (λ) n,k j u (λ) J 1 n,k j+1 ] = Imln u (λ) n,k j u (λ) n,k j+1 j=0 Note, the product u (λ) n,k 0 u (λ) n,k 1 u (λ) n,k 1 u (λ) n,k 2... does not depend on the phase of u (λ) n,k j s gauge independent Can be considered a 1D Wilson loop for a single non-degenerate band 42 / 51
107 Overview Reminder: Bloch functions, periodic boundary condition Wannier functions: introduction Maximally localized Wannier functions Modern theory of polarization Polarization and Wannier centers Calculation of Berry phase in discrete k space Position operator and periodic boundary conditions Calculation of Wannier centers in 1D systems 43 / 51
108 A subtle issue: position operator and periodic boundary conditions In 1D the eigenfunctions of the operator PˆxP, where P is a projector onto a set of bands, yields a localized set of Wannier functions. The eigenvalues of P ˆxP are the Wannier centers. 44 / 51
109 A subtle issue: position operator and periodic boundary conditions In 1D the eigenfunctions of the operator PˆxP, where P is a projector onto a set of bands, yields a localized set of Wannier functions. The eigenvalues of P ˆxP are the Wannier centers. The projectors are most easily calculated using Bloch functions: P = n,k x Ψ n,kx Ψ n,kx = n,k x u n,kx u n,kx 44 / 51
110 A subtle issue: position operator and periodic boundary conditions In 1D the eigenfunctions of the operator PˆxP, where P is a projector onto a set of bands, yields a localized set of Wannier functions. The eigenvalues of P ˆxP are the Wannier centers. The projectors are most easily calculated using Bloch functions: P = n,k x Ψ n,kx Ψ n,kx = n,k x u n,kx u n,kx The projectors u n,kx u n,kx are periodic in x because u n,kx (x) is periodic. However, the operator ˆx is not periodic. 44 / 51
111 A subtle issue: position operator and periodic boundary conditions In 1D the eigenfunctions of the operator PˆxP, where P is a projector onto a set of bands, yields a localized set of Wannier functions. The eigenvalues of P ˆxP are the Wannier centers. The projectors are most easily calculated using Bloch functions: P = n,k x Ψ n,kx Ψ n,kx = n,k x u n,kx u n,kx The projectors u n,kx u n,kx are periodic in x because u n,kx (x) is periodic. However, the operator ˆx is not periodic. We want to work in the Hilbert space spanned by u n,kx (x). 44 / 51
112 A subtle issue: position operator and periodic boundary conditions In 1D the eigenfunctions of the operator PˆxP, where P is a projector onto a set of bands, yields a localized set of Wannier functions. The eigenvalues of P ˆxP are the Wannier centers. The projectors are most easily calculated using Bloch functions: P = n,k x Ψ n,kx Ψ n,kx = n,k x u n,kx u n,kx The projectors u n,kx u n,kx are periodic in x because u n,kx (x) is periodic. However, the operator ˆx is not periodic. We want to work in the Hilbert space spanned by u n,kx (x). We need to find a periodic form of PˆxP 44 / 51
113 A subtle issue: position operator and periodic boundary conditions Raffaele Resta [3]: if periodic boundary conditions are used, expectation values that involve the operator ˆx should be calculated using the unitary operator ˆX = e i 2π L ˆx L: imposed periodicity in the system 45 / 51
114 A subtle issue: position operator and periodic boundary conditions Raffaele Resta [3]: if periodic boundary conditions are used, expectation values that involve the operator ˆx should be calculated using the unitary operator ˆX = e i 2π L ˆx L: imposed periodicity in the system to obtain well-localized Wannier functions while using periodic boundary conditions: find the eigenfunction(s) of ˆX P = P ˆXP P: projects onto the occupied bands (see Phys Rev B 84, ) 45 / 51
115 A subtle issue: position operator and periodic boundary conditions Raffaele Resta [3]: if periodic boundary conditions are used, expectation values that involve the operator ˆx should be calculated using the unitary operator ˆX = e i 2π L ˆx L: imposed periodicity in the system to obtain well-localized Wannier functions while using periodic boundary conditions: find the eigenfunction(s) of ˆX P = P ˆXP P: projects onto the occupied bands (see Phys Rev B 84, ) Note: in general ˆXP is not a Hermitian operator, only for L 45 / 51
116 Overview Reminder: Bloch functions, periodic boundary condition Wannier functions: introduction Maximally localized Wannier functions Modern theory of polarization Polarization and Wannier centers Calculation of Berry phase in discrete k space Position operator and periodic boundary conditions Calculation of Wannier centers in 1D systems 46 / 51
117 P ˆXP operator in 1D lattice system For a 1D system consisting of M unit cells: M ˆx = R m m,α m,α m=1 α α band index, R m labels the mth unit cell 47 / 51
118 P ˆXP operator in 1D lattice system For a 1D system consisting of M unit cells: M ˆx = R m m,α m,α m=1 α α band index, R m labels the mth unit cell M ˆX = e iδ kr m m,α m,α m=1 α δ k = 2π Ma a: lattice constant 47 / 51
119 P ˆXP operator in 1D lattice system For simplicity, consider a single occupied band and the corresponding Bloch functions Ψ k 48 / 51
120 P ˆXP operator in 1D lattice system For simplicity, consider a single occupied band and the corresponding Bloch functions Ψ k First, the matrix elements of ˆX: Ψ k ˆX Ψ k = M m =1 e i k R m = 1 M u(k ) u(k) M m u(k ) M m=1 = u(k ) u(k) δ k+δk k,0 M e iδ kr m e i k Rm m u(k) m=1 e i(k+δ k k )R m 48 / 51
121 P ˆXP operator in 1D lattice system For simplicity, consider a single occupied band and the corresponding Bloch functions Ψ k First, the matrix elements of ˆX: Ψ k ˆX Ψ k = M m =1 e i k R m = 1 M u(k ) u(k) M m u(k ) M m=1 = u(k ) u(k) δ k+δk k,0 M e iδ kr m e i k Rm m u(k) m=1 e i(k+δ k k )R m ˆX P = k,k Ψ k Ψ k ˆX Ψ k Ψ k = k u(k +δ k ) u(k) Ψ k+δk Ψ k 48 / 51
122 P ˆXP operator in 1D lattice system For a system consisting of M unit cells, the k values are discrete and there are M x Bloch-wavefunctions Ψ km : k m = 2π Mam, m = 1...M 49 / 51
123 P ˆXP operator in 1D lattice system For a system consisting of M unit cells, the k values are discrete and there are M x Bloch-wavefunctions Ψ km : k m = 2π Mam, m = 1...M Taking matrix elements of ˆX P in the basis of Bloch wavefunctions: X P = u(k 1 ) u(k M ) u(k 2 ) u(k 1 ) u(k 3 ) u(k 2 ) u(k 4 ) u(k 3 ) X P is a M M matrix 49 / 51
124 P ˆXP operator in 1D lattice system For a system consisting of M unit cells, the k values are discrete and there are M x Bloch-wavefunctions Ψ km : k m = 2π Mam, m = 1...M Taking matrix elements of ˆX P in the basis of Bloch wavefunctions: X P = u(k 1 ) u(k M ) u(k 2 ) u(k 1 ) u(k 3 ) u(k 2 ) u(k 4 ) u(k 3 ) X P is a M M matrix One can easily check, that (X P ) M = w1, where w C and 1 is the unit matrix 49 / 51
125 P ˆXP operator in 1D lattice system For a system consisting of M unit cells, the k values are discrete and there are M x Bloch-wavefunctions Ψ km : k m = 2π Mam, m = 1...M Taking matrix elements of ˆX P in the basis of Bloch wavefunctions: X P = u(k 1 ) u(k M ) u(k 2 ) u(k 1 ) u(k 3 ) u(k 2 ) u(k 4 ) u(k 3 ) X P is a M M matrix One can easily check, that (X P ) M = w1, where w C and 1 is the unit matrix Note, we used that u(k 2 ) u(k 1 ), u(k 3 ) u(k 2 ) etc are complex numbers they commmute 49 / 51
126 P ˆXP operator in 1D lattice system (X P ) M = w1 where w = u(k 1 ) u(k M ) u(k 2 ) u(k 1 ) u(k 3 ) u(k 2 )... u(k M ) u(k M 1 ) w is a 1D Wilson loop, gauge independent 50 / 51
127 P ˆXP operator in 1D lattice system (X P ) M = w1 where w = u(k 1 ) u(k M ) u(k 2 ) u(k 1 ) u(k 3 ) u(k 2 )... u(k M ) u(k M 1 ) w is a 1D Wilson loop, gauge independent Similar to the discrete Berry-phase 50 / 51
128 P ˆXP operator in 1D lattice system (X P ) M = w1 where w = u(k 1 ) u(k M ) u(k 2 ) u(k 1 ) u(k 3 ) u(k 2 )... u(k M ) u(k M 1 ) w is a 1D Wilson loop, gauge independent Similar to the discrete Berry-phase The eigenvalues λ m of X P are the M roots of w: w = w e iθ θ = Imln[w] 50 / 51
129 P ˆXP operator in 1D lattice system (X P ) M = w1 where w = u(k 1 ) u(k M ) u(k 2 ) u(k 1 ) u(k 3 ) u(k 2 )... u(k M ) u(k M 1 ) w is a 1D Wilson loop, gauge independent Similar to the discrete Berry-phase The eigenvalues λ m of X P are the M roots of w: w = w e iθ θ = Imln[w] λ m = M w e i(2πm+θ)/m, m = 1,...M 50 / 51
130 P ˆXP operator in 1D lattice system (X P ) M = w1 where w = u(k 1 ) u(k M ) u(k 2 ) u(k 1 ) u(k 3 ) u(k 2 )... u(k M ) u(k M 1 ) w is a 1D Wilson loop, gauge independent Similar to the discrete Berry-phase The eigenvalues λ m of X P are the M roots of w: w = w e iθ θ = Imln[w] λ m = M w e i(2πm+θ)/m, m = 1,...M Note: since ˆX P is not Hermitian, the eigenvectors are not orthogonal in general, only in the M limit (we will actually not need the eigenstates) 50 / 51
131 Wannier centers and the eigenvalues of P ˆXP Remember from earlier: 1) polarization of a single filled band in a 1D lattice x j Wannier centers, x j = W(j) x W(j), j is unit cell index 51 / 51
132 Wannier centers and the eigenvalues of P ˆXP Remember from earlier: 1) polarization of a single filled band in a 1D lattice x j Wannier centers, x j = W(j) x W(j), j is unit cell index 2) In discrete k space J 1 x Imln u n,kj u n,kj+1 j=0 j different discrete value of k in the 1D Brillouin zone 51 / 51
133 Wannier centers and the eigenvalues of P ˆXP Remember from earlier: 1) polarization of a single filled band in a 1D lattice x j Wannier centers, x j = W(j) x W(j), j is unit cell index 2) In discrete k space J 1 x Imln u n,kj u n,kj+1 j=0 j different discrete value of k in the 1D Brillouin zone If there are M unit cells, k is discretized as k m = 2π Ma m The Wannier centers x m can be obtained from λ m as x m = M 2π Imln[λ m] = θ Imln[w] +m = +m 2π 2π 51 / 51
Adiabatic particle pumping and anomalous velocity
Adiabatic particle pumping and anomalous velocity November 17, 2015 November 17, 2015 1 / 31 Literature: 1 J. K. Asbóth, L. Oroszlány, and A. Pályi, arxiv:1509.02295 2 D. Xiao, M-Ch Chang, and Q. Niu,
More informationTopological Physics in Band Insulators IV
Topological Physics in Band Insulators IV Gene Mele University of Pennsylvania Wannier representation and band projectors Modern view: Gapped electronic states are equivalent Kohn (1964): insulator is
More informationBerry phase, Chern number
Berry phase, Chern number November 17, 2015 November 17, 2015 1 / 22 Literature: 1 J. K. Asbóth, L. Oroszlány, and A. Pályi, arxiv:1509.02295 2 D. Xiao, M-Ch Chang, and Q. Niu, Rev. Mod. Phys. 82, 1959.
More informationElectrons in a periodic potential
Chapter 3 Electrons in a periodic potential 3.1 Bloch s theorem. We consider in this chapter electrons under the influence of a static, periodic potential V (x), i.e. such that it fulfills V (x) = V (x
More informationarxiv:cond-mat/ v1 [cond-mat.mtrl-sci] 17 Jan 1998
Electronic polarization in the ultrasoft pseudopotential formalism arxiv:cond-mat/9801177v1 [cond-mat.mtrl-sci] 17 Jan 1998 David anderbilt Department of Physics and Astronomy, Rutgers University, Piscataway,
More informationLecture 4: Basic elements of band theory
Phys 769 Selected Topics in Condensed Matter Physics Summer 010 Lecture 4: Basic elements of band theory Lecturer: Anthony J. Leggett TA: Bill Coish 1 Introduction Most matter, in particular most insulating
More informationCHEM6085: Density Functional Theory
Lecture 11 CHEM6085: Density Functional Theory DFT for periodic crystalline solids C.-K. Skylaris 1 Electron in a one-dimensional periodic box (in atomic units) Schrödinger equation Energy eigenvalues
More informationTwo Dimensional Chern Insulators, the Qi-Wu-Zhang and Haldane Models
Two Dimensional Chern Insulators, the Qi-Wu-Zhang and Haldane Models Matthew Brooks, Introduction to Topological Insulators Seminar, Universität Konstanz Contents QWZ Model of Chern Insulators Haldane
More informationJ. Phys. Chem. Solids 61, 147 (2000). Berry-phase theory of proper piezoelectric response David Vanderbilt Department of Physics and Astronomy, Rutgers University, Piscataway, New Jersey 08855-0849 (July
More informationTight binding models from band representations
Tight binding models from band representations Jennifer Cano Stony Brook University and Flatiron Institute for Computational Quantum Physics Part 0: Review of band representations BANDREP http://www.cryst.ehu.es/cgi-bin/cryst/programs/bandrep.pl
More informationFloquet Topological Insulator:
Floquet Topological Insulator: Understanding Floquet topological insulator in semiconductor quantum wells by Lindner et al. Condensed Matter Journal Club Caltech February 12 2014 Motivation Motivation
More informationNumerical construction of Wannier functions
July 12, 2017 Internship Tutor: A. Levitt School Tutor: É. Cancès 1/27 Introduction Context: Describe electrical properties of crystals (insulator, conductor, semi-conductor). Applications in electronics
More informationElectron levels in a periodic potential: general remarks. Daniele Toffoli January 11, / 46
Electron levels in a periodic potential: general remarks Daniele Toffoli January 11, 2017 1 / 46 Outline 1 Mathematical tools 2 The periodic potential 3 Bloch s theorem and Born-von Karman boundary conditions
More informationIntroduction to topological insulators. Jennifer Cano
Introduction to topological insulators Jennifer Cano Adapted from Charlie Kane s Windsor Lectures: http://www.physics.upenn.edu/~kane/ Review article: Hasan & Kane Rev. Mod. Phys. 2010 What is an insulator?
More informationFirst-Principles Calculation of Topological Invariants (Wannier Functions Approach) Alexey A. Soluyanov
First-Principles Calculation of Topological Invariants (Wannier Functions Approach) Alexey A. Soluyanov ES'12, WFU, June 8, 212 The present work was done in collaboration with David Vanderbilt Outline:
More informationTime-Reversal Symmetric Two-Dimensional Topological Insulators: The Bernevig-Hughes-Zhang Model
Time-Reversal Symmetric Two-Dimensional Topological Insulators: The Bernevig-Hughes-Zhang Model Alexander Pearce Intro to Topological Insulators: Week 5 November 26, 2015 1 / 22 This notes are based on
More informationThe Quantum Heisenberg Ferromagnet
The Quantum Heisenberg Ferromagnet Soon after Schrödinger discovered the wave equation of quantum mechanics, Heisenberg and Dirac developed the first successful quantum theory of ferromagnetism W. Heisenberg,
More informationBerry s phase in Hall Effects and Topological Insulators
Lecture 6 Berry s phase in Hall Effects and Topological Insulators Given the analogs between Berry s phase and vector potentials, it is not surprising that Berry s phase can be important in the Hall effect.
More informationSymmetry, Topology and Phases of Matter
Symmetry, Topology and Phases of Matter E E k=λ a k=λ b k=λ a k=λ b Topological Phases of Matter Many examples of topological band phenomena States adiabatically connected to independent electrons: - Quantum
More informationarxiv: v1 [cond-mat.other] 20 Apr 2010
Characterization of 3d topological insulators by 2d invariants Rahul Roy Rudolf Peierls Centre for Theoretical Physics, 1 Keble Road, Oxford, OX1 3NP, UK arxiv:1004.3507v1 [cond-mat.other] 20 Apr 2010
More informationIntroduction to Condensed Matter Physics
Introduction to Condensed Matter Physics The Reciprocal Lattice M.P. Vaughan Overview Overview of the reciprocal lattice Periodic functions Reciprocal lattice vectors Bloch functions k-space Dispersion
More informationProof that the maximally localized Wannier functions are real
1 Proof that the maximally localized Wannier functions are real Sangryol Ri, Suil Ri Institute of Physics, Academy of Sciences, Unjong District, Pyongyang, DPR Korea Abstract The maximally localized Wannier
More informationTheory of Polarization: A Modern Approach
Theory of Polarization: A Modern Approach Raffaele Resta 1,2 and David Vanderbilt 3 1 INFM DEMOCRITOS National Simulation Center, Via Beirut 4, I-34014, Trieste, Italy resta@democritos.it 2 Dipartimento
More informationElectronic Structure of Crystalline Solids
Electronic Structure of Crystalline Solids Computing the electronic structure of electrons in solid materials (insulators, conductors, semiconductors, superconductors) is in general a very difficult problem
More informationWannier functions. Macroscopic polarization (Berry phase) and related properties. Effective band structure of alloys
Wannier functions Macroscopic polarization (Berry phase) and related properties Effective band structure of alloys P.Blaha (from Oleg Rubel, McMaster Univ, Canada) Wannier functions + + Wannier90: A Tool
More informationElectrons in Crystals. Chris J. Pickard
Electrons in Crystals Chris J. Pickard Electrons in Crystals The electrons in a crystal experience a potential with the periodicity of the Bravais lattice: U(r + R) = U(r) The scale of the periodicity
More informationPOEM: Physics of Emergent Materials
POEM: Physics of Emergent Materials Nandini Trivedi L1: Spin Orbit Coupling L2: Topology and Topological Insulators Reference: Bernevig Topological Insulators and Topological Superconductors Tutorials:
More informationTopological insulators. Pavel Buividovich (Regensburg)
Topological insulators Pavel Buividovich (Regensburg) Hall effect Classical treatment Dissipative motion for point-like particles (Drude theory) Steady motion Classical Hall effect Cyclotron frequency
More informationHomogeneous Electric and Magnetic Fields in Periodic Systems
Electric and Magnetic Fields in Periodic Systems Josef W. idepartment of Chemistry and Institute for Research in Materials Dalhousie University Halifax, Nova Scotia June 2012 1/24 Acknowledgments NSERC,
More informationTopological Properties of Quantum States of Condensed Matter: some recent surprises.
Topological Properties of Quantum States of Condensed Matter: some recent surprises. F. D. M. Haldane Princeton University and Instituut Lorentz 1. Berry phases, zero-field Hall effect, and one-way light
More informationMD simulation: output
Properties MD simulation: output Trajectory of atoms positions: e. g. diffusion, mass transport velocities: e. g. v-v autocorrelation spectrum Energies temperature displacement fluctuations Mean square
More informationElectronic Structure Methodology 1
Electronic Structure Methodology 1 Chris J. Pickard Lecture Two Working with Density Functional Theory In the last lecture we learnt how to write the total energy as a functional of the density n(r): E
More informationLecture 6 Photons, electrons and other quanta. EECS Winter 2006 Nanophotonics and Nano-scale Fabrication P.C.Ku
Lecture 6 Photons, electrons and other quanta EECS 598-002 Winter 2006 Nanophotonics and Nano-scale Fabrication P.C.Ku From classical to quantum theory In the beginning of the 20 th century, experiments
More information4 Electric polarization
4 Electric polarization Clearly the electric polarization P of a crystalline insulator should be defined so as to carry the meaning of dipole moment per unit volume, but a precise formulation is surprisingly
More informationOrbital magnetization in insulators with broken time-reversal symmetry. David Vanderbilt Rutgers University
Orbital magnetization in insulators with broken time-reversal symmetry David Vanderbilt Rutgers University Collaboration Collaboration Timo Thonhauser (Rutgers) David Vanderbilt (Rutgers) Davide Ceresoli
More informationBerry-phase Approach to Electric Polarization and Charge Fractionalization. Dennis P. Clougherty Department of Physics University of Vermont
Berry-phase Approach to Electric Polarization and Charge Fractionalization Dennis P. Clougherty Department of Physics University of Vermont Outline Quick Review Berry phase in quantum systems adiabatic
More informationwoptic: Transport Properties with Wannier Functions and Adaptive k-integration
woptic: Transport Properties with Wannier Functions and Adaptive k-integration Elias Assmann Institute of Solid State Physics, Vienna University of Technology Split, 2013-09-29 About this presentation
More informationTopological Insulators and Superconductors
Topological Insulators and Superconductors Lecture #1: Topology and Band Theory Lecture #: Topological Insulators in and 3 dimensions Lecture #3: Topological Superconductors, Majorana Fermions an Topological
More informationOrganizing Principles for Understanding Matter
Organizing Principles for Understanding Matter Symmetry Conceptual simplification Conservation laws Distinguish phases of matter by pattern of broken symmetries Topology Properties insensitive to smooth
More informationFirst-Principles Wannier Functions of Silicon and Gallium. Arsenide arxiv:cond-mat/ v1 [cond-mat.mtrl-sci] 22 Nov 1996.
First-Principles Wannier Functions of Silicon and Gallium Arsenide arxiv:cond-mat/9611176v1 [cond-mat.mtrl-sci] 22 Nov 1996 Pablo Fernández 1, Andrea Dal Corso 1, Francesco Mauri 2, and Alfonso Baldereschi
More informationC2: Band structure. Carl-Olof Almbladh, Rikard Nelander, and Jonas Nyvold Pedersen Department of Physics, Lund University.
C2: Band structure Carl-Olof Almbladh, Rikard Nelander, and Jonas Nyvold Pedersen Department of Physics, Lund University December 2005 1 Introduction When you buy a diamond for your girl/boy friend, what
More informationLecture 10 Waves in Periodic Potentials Today: Questions you should be able to address after today s lecture:
Lecture 10 Waves in Periodic Potentials Today: 1. Direct lattice and periodic potential as a convolution of a lattice and a basis. 2. The discrete translation operator: eigenvalues and eigenfunctions.
More information1 Mathematical preliminaries
1 Mathematical preliminaries The mathematical language of quantum mechanics is that of vector spaces and linear algebra. In this preliminary section, we will collect the various definitions and mathematical
More informationwien2wannier and woptic: From Wannier Functions to Optical Conductivity
wien2wannier and woptic: From Wannier Functions to Optical Conductivity Elias Assmann Institute of Solid State Physics, Vienna University of Technology AToMS-2014, Bariloche, Aug 4 Outline brief introduction
More informationInterband effects and orbital suceptibility of multiband tight-binding models
Interband effects and orbital suceptibility of multiband tight-binding models Frédéric Piéchon LPS (Orsay) with A. Raoux, J-N. Fuchs and G. Montambaux Orbital suceptibility Berry curvature ky? kx GDR Modmat,
More information1 Commutators (10 pts)
Final Exam Solutions 37A Fall 0 I. Siddiqi / E. Dodds Commutators 0 pts) ) Consider the operator  = Ĵx Ĵ y + ĴyĴx where J i represents the total angular momentum in the ith direction. a) Express both
More informationElectrons in periodic potential
Chapter 9 Electrons in periodic potential The computation of electronic states in a solid is a nontrivial problem. A great simplification can be achieved if the solid is a crystal, i.e. if it can be described
More informationQuantum Physics III (8.06) Spring 2007 FINAL EXAMINATION Monday May 21, 9:00 am You have 3 hours.
Quantum Physics III (8.06) Spring 2007 FINAL EXAMINATION Monday May 21, 9:00 am You have 3 hours. There are 10 problems, totalling 180 points. Do all problems. Answer all problems in the white books provided.
More informationTopological phases of matter give rise to quantized physical quantities
Quantized electric multipole insulators Benalcazar, W. A., Bernevig, B. A., & Hughes, T. L. (2017). Quantized electric multipole insulators. Science, 357(6346), 61 66. Presented by Mark Hirsbrunner, Weizhan
More informationLes états de bord d un. isolant de Hall atomique
Les états de bord d un isolant de Hall atomique séminaire Atomes Froids 2/9/22 Nathan Goldman (ULB), Jérôme Beugnon and Fabrice Gerbier Outline Quantum Hall effect : bulk Landau levels and edge states
More informationBerry Phase Effects on Charge and Spin Transport
Berry Phase Effects on Charge and Spin Transport Qian Niu 牛谦 University of Texas at Austin 北京大学 Collaborators: Shengyuan Yang, C.P. Chuu, D. Xiao, W. Yao, D. Culcer, J.R.Shi, Y.G. Yao, G. Sundaram, M.C.
More informationEnergy bands in two limits
M. A. Gusmão IF-UFRGS 1 FIP10601 Text 4 Energy bands in two limits After presenting a general view of the independent-electron approximation, highlighting the importance and consequences of Bloch s theorem,
More informationproperties Michele Catti Dipartimento di Scienza dei Materiali Università di Milano Bicocca, Italy
Elastic and piezoelectric tensorial properties Michele Catti Dipartimento di Scienza dei Materiali Università di Milano Bicocca, Italy (catti@mater.unimib.it) 1 Tensorial physical properties of crystals
More informationMaximally localized Wannier functions: Theory and applications
Maximally localized Wannier functions: Theory and applications Nicola Marzari Theory and Simulation of Materials (THEOS), École Polytechnique Fédérale de Lausanne, Station 12, 1015 Lausanne, Switzerland
More informationQuantum Mechanics Solutions. λ i λ j v j v j v i v i.
Quantum Mechanics Solutions 1. (a) If H has an orthonormal basis consisting of the eigenvectors { v i } of A with eigenvalues λ i C, then A can be written in terms of its spectral decomposition as A =
More informationDensity Functional Theory: from theory to Applications
Density Functional Theory: from theory to Applications Uni Mainz November 29, 2010 The self interaction error and its correction Perdew-Zunger SIC Average-density approximation Weighted density approximation
More informationPhysics 221A Fall 1996 Notes 16 Bloch s Theorem and Band Structure in One Dimension
Physics 221A Fall 1996 Notes 16 Bloch s Theorem and Band Structure in One Dimension In these notes we examine Bloch s theorem and band structure in problems with periodic potentials, as a part of our survey
More informationMP463 QUANTUM MECHANICS
MP463 QUANTUM MECHANICS Introduction Quantum theory of angular momentum Quantum theory of a particle in a central potential - Hydrogen atom - Three-dimensional isotropic harmonic oscillator (a model of
More informationExploring topological states with cold atoms and photons
Exploring topological states with cold atoms and photons Theory: Takuya Kitagawa, Dima Abanin, Erez Berg, Mark Rudner, Liang Fu, Takashi Oka, Immanuel Bloch, Eugene Demler Experiments: I. Bloch s group
More informationĝ r = {R v} r = R r + v.
SUPPLEMENTARY INFORMATION DOI: 1.138/NPHYS134 Topological semimetal in a fermionic optical lattice Kai Sun, 1 W. Vincent Liu,, 3, 4 Andreas Hemmerich, 5 and S. Das Sarma 1 1 Condensed Matter Theory Center
More informationHIGHER INVARIANTS: TOPOLOGICAL INSULATORS
HIGHER INVARIANTS: TOPOLOGICAL INSULATORS Sponsoring This material is based upon work supported by the National Science Foundation Grant No. DMS-1160962 Jean BELLISSARD Georgia Institute of Technology,
More informationReciprocal Space Magnetic Field: Physical Implications
Reciprocal Space Magnetic Field: Physical Implications Junren Shi ddd Institute of Physics Chinese Academy of Sciences November 30, 2005 Outline Introduction Implications Conclusion 1 Introduction 2 Physical
More informationThe 3 dimensional Schrödinger Equation
Chapter 6 The 3 dimensional Schrödinger Equation 6.1 Angular Momentum To study how angular momentum is represented in quantum mechanics we start by reviewing the classical vector of orbital angular momentum
More informationSSH Model. Alessandro David. November 3, 2016
SSH Model Alessandro David November 3, 2016 Adapted from Lecture Notes at: https://arxiv.org/abs/1509.02295 and from article: Nature Physics 9, 795 (2013) Motivations SSH = Su-Schrieffer-Heeger Polyacetylene
More informationPage 404. Lecture 22: Simple Harmonic Oscillator: Energy Basis Date Given: 2008/11/19 Date Revised: 2008/11/19
Page 404 Lecture : Simple Harmonic Oscillator: Energy Basis Date Given: 008/11/19 Date Revised: 008/11/19 Coordinate Basis Section 6. The One-Dimensional Simple Harmonic Oscillator: Coordinate Basis Page
More information3.15. Some symmetry properties of the Berry curvature and the Chern number.
50 Phys620.nb z M 3 at the K point z M 3 3 t ' sin 3 t ' sin (3.36) (3.362) Therefore, as long as M 3 3 t ' sin, the system is an topological insulator ( z flips sign). If M 3 3 t ' sin, z is always positive
More informationDirac semimetal in three dimensions
Dirac semimetal in three dimensions Steve M. Young, Saad Zaheer, Jeffrey C. Y. Teo, Charles L. Kane, Eugene J. Mele, and Andrew M. Rappe University of Pennsylvania 6/7/12 1 Dirac points in Graphene Without
More informationTopological Phases in Floquet Systems
Rahul Roy University of California, Los Angeles June 2, 2016 arxiv:1602.08089, arxiv:1603.06944 Post-doc: Fenner Harper Outline 1 Introduction 2 Free Fermionic Systems 3 Interacting Systems in Class D
More informationLecture. Ref. Ihn Ch. 3, Yu&Cardona Ch. 2
Lecture Review of quantum mechanics, statistical physics, and solid state Band structure of materials Semiconductor band structure Semiconductor nanostructures Ref. Ihn Ch. 3, Yu&Cardona Ch. 2 Reminder
More informationTopological Phases of Matter Out of Equilibrium
Topological Phases of Matter Out of Equilibrium Nigel Cooper T.C.M. Group, Cavendish Laboratory, University of Cambridge Solvay Workshop on Quantum Simulation ULB, Brussels, 18 February 2019 Max McGinley
More informationCoordinate and Momentum Representation. Commuting Observables and Simultaneous Measurements. January 30, 2012
Coordinate and Momentum Representation. Commuting Observables and Simultaneous Measurements. January 30, 2012 1 Coordinate and Momentum Representations Let us consider an eigenvalue problem for a Hermitian
More informationPOEM: Physics of Emergent Materials
POEM: Physics of Emergent Materials Nandini Trivedi L1: Spin Orbit Coupling L2: Topology and Topological Insulators Tutorials: May 24, 25 (2017) Scope of Lectures and Anchor Points: 1.Spin-Orbit Interaction
More information-state problems and an application to the free particle
-state problems and an application to the free particle Sourendu Gupta TIFR, Mumbai, India Quantum Mechanics 1 2013 3 September, 2013 Outline 1 Outline 2 The Hilbert space 3 A free particle 4 Keywords
More informationv. Tε n k =ε n k T r T = r, T v T = r, I v I = I r I = v. Iε n k =ε n k Berry curvature: Symmetry Consideration n k = n k
Berry curvature: Symmetry Consideration Time reversal (i.e. motion reversal) 1 1 T r T = r, T v T = v. Tε n k =ε n k n k = n k Inversion Symmetry: 1 1 I r I = r, I v I = v. Iε n k =ε n k n k = n k θ
More informationIntroduction to Quantum Mechanics PVK - Solutions. Nicolas Lanzetti
Introduction to Quantum Mechanics PVK - Solutions Nicolas Lanzetti lnicolas@student.ethz.ch 1 Contents 1 The Wave Function and the Schrödinger Equation 3 1.1 Quick Checks......................................
More informationLinear Algebra in Hilbert Space
Physics 342 Lecture 16 Linear Algebra in Hilbert Space Lecture 16 Physics 342 Quantum Mechanics I Monday, March 1st, 2010 We have seen the importance of the plane wave solutions to the potentialfree Schrödinger
More informationarxiv: v1 [cond-mat.mes-hall] 26 Sep 2013
Berry phase and the unconventional quantum Hall effect in graphene Jiamin Xue Microelectronic Research Center, The University arxiv:1309.6714v1 [cond-mat.mes-hall] 26 Sep 2013 of Texas at Austin, Austin,
More informationTopological insulator with time-reversal symmetry
Phys620.nb 101 7 Topological insulator with time-reversal symmetry Q: Can we get a topological insulator that preserves the time-reversal symmetry? A: Yes, with the help of the spin degree of freedom.
More informationMean-Field Theory: HF and BCS
Mean-Field Theory: HF and BCS Erik Koch Institute for Advanced Simulation, Jülich N (x) = p N! ' (x ) ' 2 (x ) ' N (x ) ' (x 2 ) ' 2 (x 2 ) ' N (x 2 )........ ' (x N ) ' 2 (x N ) ' N (x N ) Slater determinant
More informationNoncollinear spins in QMC: spiral Spin Density Waves in the HEG
Noncollinear spins in QMC: spiral Spin Density Waves in the HEG Zoltán Radnai and Richard J. Needs Workshop at The Towler Institute July 2006 Overview What are noncollinear spin systems and why are they
More informationOn the K-theory classification of topological states of matter
On the K-theory classification of topological states of matter (1,2) (1) Department of Mathematics Mathematical Sciences Institute (2) Department of Theoretical Physics Research School of Physics and Engineering
More informationPLEASE LET ME KNOW IF YOU FIND TYPOS (send to
Teoretisk Fysik KTH Advanced QM (SI2380), Lecture 2 (Summary of concepts) 1 PLEASE LET ME KNOW IF YOU FIND TYPOS (send email to langmann@kth.se) The laws of QM 1. I now discuss the laws of QM and their
More informationStatistical Interpretation
Physics 342 Lecture 15 Statistical Interpretation Lecture 15 Physics 342 Quantum Mechanics I Friday, February 29th, 2008 Quantum mechanics is a theory of probability densities given that we now have an
More informationMaximally localized Wannier functions for periodic Schrödinger operators: the MV functional and the geometry of the Bloch bundle
Maximally localized Wannier functions for periodic Schrödinger operators: the MV functional and the geometry of the Bloch bundle Gianluca Panati Università di Roma La Sapienza 4ème Rencontre du GDR Dynamique
More informationKitaev honeycomb lattice model: from A to B and beyond
Kitaev honeycomb lattice model: from A to B and beyond Jiri Vala Department of Mathematical Physics National University of Ireland at Maynooth Postdoc: PhD students: Collaborators: Graham Kells Ahmet Bolukbasi
More informationLecture 4: Equations of motion and canonical quantization Read Sakurai Chapter 1.6 and 1.7
Lecture 4: Equations of motion and canonical quantization Read Sakurai Chapter 1.6 and 1.7 In Lecture 1 and 2, we have discussed how to represent the state of a quantum mechanical system based the superposition
More informationLecture 8 Nature of ensemble: Role of symmetry, interactions and other system conditions: Part II
Lecture 8 Nature of ensemble: Role of symmetry, interactions and other system conditions: Part II We continue our discussion of symmetries and their role in matrix representation in this lecture. An example
More informationCompressed representation of Kohn-Sham orbitals via selected columns of the density matrix
Lin Lin Compressed Kohn-Sham Orbitals 1 Compressed representation of Kohn-Sham orbitals via selected columns of the density matrix Lin Lin Department of Mathematics, UC Berkeley; Computational Research
More informationarxiv: v1 [physics.comp-ph] 25 Jan 2018
VARIATIONAL FORMULATION FOR WANNIER FUNCTIONS WITH ENTANGLED BAND STRUCTURE ANIL DAMLE, ANTOINE LEVITT, AND LIN LIN arxiv:1801.08572v1 [physics.comp-ph] 25 Jan 2018 Abstract. Wannier functions provide
More informationarxiv: v1 [cond-mat.mes-hall] 21 Oct 2010
Inversion Symmetric Topological Insulators arxiv:.458v [cond-mat.mes-hall] Oct Taylor L. Hughes, Emil Prodan, and B. Andrei Bernevig 3 Department of Physics, University of Illinois, West Green St, Urbana
More informationThree Most Important Topics (MIT) Today
Three Most Important Topics (MIT) Today Electrons in periodic potential Energy gap nearly free electron Bloch Theorem Energy gap tight binding Chapter 1 1 Electrons in Periodic Potential We now know the
More informationBloch Wilson Hamiltonian and a Generalization of the Gell-Mann Low Theorem 1
Bloch Wilson Hamiltonian and a Generalization of the Gell-Mann Low Theorem 1 Axel Weber 2 arxiv:hep-th/9911198v1 25 Nov 1999 Instituto de Física y Matemáticas, Universidad Michoacana de San Nicolás de
More information2 The Density Operator
In this chapter we introduce the density operator, which provides an alternative way to describe the state of a quantum mechanical system. So far we have only dealt with situations where the state of a
More informationMeasuring many-body topological invariants using polarons
1 Anyon workshop, Kaiserslautern, 12/15/2014 Measuring many-body topological invariants using polarons Fabian Grusdt Physics Department and Research Center OPTIMAS, University of Kaiserslautern, Germany
More informationLecture #1. Review. Postulates of quantum mechanics (1-3) Postulate 1
L1.P1 Lecture #1 Review Postulates of quantum mechanics (1-3) Postulate 1 The state of a system at any instant of time may be represented by a wave function which is continuous and differentiable. Specifically,
More information1 Planck-Einstein Relation E = hν
C/CS/Phys C191 Representations and Wavefunctions 09/30/08 Fall 2008 Lecture 8 1 Planck-Einstein Relation E = hν This is the equation relating energy to frequency. It was the earliest equation of quantum
More informationWe start with some important background material in classical and quantum mechanics.
Chapter Basics We start with some important background material in classical and quantum mechanics.. Classical mechanics Lagrangian mechanics Compared to Newtonian mechanics, Lagrangian mechanics has the
More informationKeywords: Bloch-Floquet Theorem, Periodic Media, Photonic Crystals, Plasmonic Crystals, Electronic Crystals, Phononic Crystals
On the Bloch Theorem and Orthogonality Relations Sina Khorasani School of Electrical Engineering Sharif University of Technology P. O. Box 11365-9363 Tehran, Iran Email: khorasani@sina.sharif.edu Abstract
More informationBerry Phases and Curvatures in Electronic-Structure Theory. David Vanderbilt Rutgers University
Berry Phases and Curvatures in Electronic-Structure Theory David Vanderbilt Rutgers University Rahman prize for: Theory of polarization (King-Smith & Vanderbilt) Ultrasoft pseudopotentials Three quick
More information3.14. The model of Haldane on a honeycomb lattice
4 Phys60.n..7. Marginal case: 4 t Dirac points at k=(,). Not an insulator. No topological index...8. case IV: 4 t All the four special points has z 0. We just use u I for the whole BZ. No singularity.
More information