Advanced Statistical Mechanics

Size: px
Start display at page:

Download "Advanced Statistical Mechanics"

Transcription

1 Advanced Statistical Mechanics Victor Gurarie Fall 0 Week 3: D and D transverse field Ising models, exact solution

2 D Ising Model: Peierls Argument D Ising model is defined on a D lattice. Each spin of the D Ising model lives on a -dimensional space, thus it has two coordinates. Let s call them horizontal and vertical coordinates, and let s denote them x and τ. So each spin can be written as σ τ x. Suppose the Ising model has length N h in the horizontal direction and N v in the vertical direction, so all together it has N h N v spins. The energy of the Ising model can be written as E = and the partition function is x=,...,n h ; τ=,...,n v Z = σ i =± [ ] Jh σxσ τ x+ τ + J v σxσ τ x τ+. e E T.. For generality we introduce two distinct couplings, J h and J v in the horizontal and vertical directions, although in the classic Ising model J h = J v = J. We also need to specify what happens at the boundary. Let s take the periodic boundary conditions. To reduce clutter, let us adopt this common notations: σ Nv+ x σ x, σ τ N h + σ τ..3 In this way, the periodic boundary conditions are automatically encoded in.. Peierls came up with a hand-waving argument which supports that at temperature below some critical temperature T c, the Ising spins order. Let us go over this argument quickly in case where J h = J v = J, and J/T = K. A droplet of negative spins in the sea of positive spins has energy E 0 + JL,.4 where E 0 is the energy of all spins pointing up, or E 0 = JN v N h N v N h is the total number of bonds on the square lattice of our size. Here L is the length of the boundary between a droplet of spins pointing down, and the rest of the spins pointing up each bond where spin at one end points up and a spin at the other point points down has energy +J, which differs from its energy when both spins point up which is J by J. So J times the number of bonds with opposite spins at ends is the energy of the droplet, relative to E 0. There are L such bonds, so the total is E 0 + JL. Its contribution to Z is then e E 0 T KL. On the other hand, there are c L ways to position a droplet. This produces a contribution to Z of the order of E E0/T L ln c KL e

3 which, at large K small temperature, as K = J/T, will always be small. In other words, the free energy of the droplet is F = E T S = JL LT ln c E 0.5 and its minimum at T < J/ ln c is achieved at zero L. This means at T < J/ ln c, there is going to be no droplets around and the spins will all be positive. At T > J/ ln c the minimum of free energy is at very large L, in other words, these droplets will be around, they will be large, and eventually they will blur the difference between spin-ups and spin-downs, so the spins will not be ordered. The transition happens at critical temperature T c which this argument says is J/ ln c. This is just a rough estimate, in practice this temperature is proportional to J with some numerical coefficient in front. If J h J v, there is still a transition temperature T c which can be expressed in terms of these. This expression is known exactly, but its form is not illuminating and so won t be given here. Transfer Matrix and the Hamiltonian We rewrite the partition function as a transfer matrix. This is a little more involved than in D. In one dimension, the transfer matrix has two indices, σ and σ. In D, we designate all spins along one row of spins as one index of the transfer matrix, and all spins along an adjacent row as the other index. Thus a transfer matrix is a matrix of the size N h by N h. The transfer matrix takes the form T σ,σ,...,σ N h σ,σ,...,σ N h Nh = e x=k vσx σ x+ +Kvσ x σ x.. Here we introduced the natural notations K h = J h T, K v = J v T.. Again, it has two indices, labeling its rows and columns. Each index goes from to N h. However, we denote the index not by one integer going between these two values, but rather by a string of integers σ, σ, σ 3,..., σ N lower index and σ, σ, σ 3,..., σ N upper index, each σ taking values ±, all for Nv distinct strings of integers.

4 Then we can write the partition function in the following way σ=± e x=,...,nh ; τ=,...,nv [K hσxσ τ x+ τ +Kvστ xσx τ+ ] = T σ...σ N h T σ3...σ3 Nh σ σ=±...σ N σ h...σ N h This in turn can be interpreted as... T σ...σ Nh σ Nv...σ Nv N h.3 Z = Tr T Nv = n λ Nv n λ Nv 0,.4 where λ n are the eigenvalues of T there are going to be N h of those because this matrix is N h by N h and λ0 is the largest of these eigenvalues. So the task is to diagonalize the matrix T, written explicitly in.. This is a complicated matrix and it s not easy to diagonalize. However, there is a particular limit in which the matrix considerably simplifies. This is the limit K v, K h. Let us derive the transfer matrix in this limit. To do that, we first of all split the transfer matrix into a product of two matrices, the one responsible for the horizontal part of the transfer matrix, and the one responsible for the vertical part, in the following way T σ,σ,...,σ N h σ,σ,...,σ N h = σ=± [T h ] σ, σ,..., σ N h σ,σ,...,σ N h [T v ] σ,σ,...,σ Nh σ, σ,..., σ N h..5 Here and [T h ] σ, σ,..., σ N h σ,σ,...,σ N h = e K Nh h [T v ] σ,σ,...,σ N h σ, σ,..., σ N h x= σ xσx+ δ σ δ σ σ σ... δ σ Nv σ Nv.6 = e Kv Nh x= σ xσ x..7 It should be clear that multiplying T h by T v, as is done in.5, gives T, as in.. Next we employ the following strategy. We will show that, under conditions K v and K h, each of these two matrices can be written as T h = e H h, T v = e Hv.8 where H h and H v are very small matrices. Now in general for two matrices A and B, e A e B e A+B. But for two very small matrices e H h e Hv H h H v H h H v e H h H v..9 So we will be able to derive T in this way. 3

5 Let s follow this route. Take T h first. Suppose K v. Then we can Taylor-expand the exponential: [T ] σ, σ,..., σ N h h = δ σ σ,σ,...,σ N σ h δ σ... δ σ N h σ σ N h +K h x= [ δ σ δ σ σ σ... δ σ x σ x σ x δ σ x σ x ] σ x+δ σ x+ δ σ x+... δ σ N h σx+ σx+ σn h.0 We rewrite this with the help of the Pauli matrices [ τ 3 ] σ x 0 x = = σ xδ σ x 0 σ.. x This gives σ x [T ] σ, σ,..., σ N h h = δ σ σ,σ,...,σ N σ h δ σ... δ σ N h σ σn h + K h x= [ δ σ δ σ σ σ ]... δ σ x τ σ xτ 3 x+δ σ 3 x+... δ σ N h x σx+ σn h. Now the first term in this sum is nothing but an identity matrix a very large, N h N h identity matrix. The second term represents a sum of terms, each is a product of many little identities acting on all spins except σ x and σ x+, which are acted upon by Pauli matrices instead. A shorthand notation is often introduced, where all these little Kronecker deltas are omitted just as multiplication by is usually omitted, which leads to the following simplified notation yet completely equivalent. [T h ] σ, σ,..., σ N h σ,σ,...,σ N h Here I stands for this very large identity matrix. The final step is to recall that K h, so we can write N = I + K h τxτ 3 x x= T h = e K h N x= τ 3 xτ 3 x+..4 So we accomplished the task of writing T h in the form e H h, where H h = K h N x= by τ 3 xτ 3 x+..5 Now on to T v. T v will be transformed in the following way. We first write [T v ] σ,σ,...,σ N h σ, σ,..., σ N h = Let us study each matrix in this product. N h x= e Kv σ x σ x..6 e Kv σ xσ x.7 4

6 is a two by two matrix, which explicitly takes the familiar from D Ising model form [ ] e K v σ xσ x σ x = σ x e K v e Kv e Kv This matrix is close to identity if K v. In this case we can write it as e Kv.8 e K v e Kv e K v = e Kv = e Kv δ + e K v τ e Kv e e Kv τ..9 e Kv e Kv e Kv Here the natural notations are Therefore, we can write τ = 0, δ = [ ] e K v σ x σ σ x x = e Kv e e Kv τx.. σ x Here the index x on τ x serves as a reminder that τ 3 is applied to the x-th spin. Now we can write [T v ] σ,σ,...,σ N h σ, σ,..., σ N h = N h x= Finally, we recall that K v, so e Kv into a big sum. So we succeeded in writing [T v ] σ,σ,...,σ N h σ, σ,..., σ N h e Kv σ x σ x = e N h K v N h x= e e Kv τ x.., so all these exponentials can be combined = e N hk v e e Lv N h x= τ x..3 T v = e KvN h e Hv, H v = e Lv x= τ x..4 This completes the derivation of T v. Finally, we put everything together. We introduce convenient notations γ = e Kv, β = K h..5 This gives T = e KvN h e H,.6 where H = γ x= τ x β 5 x= τ 3 xτ 3 x+..7

7 This Hamiltonian is called the D transverse field Ising model. Indeed, this is a quantum mechanical Hamiltonian, which acts on quantum spins aligned in a D chain. τ x are spin operators Pauli matrices acting on x-th specific spin. We have a term τxτ 3 x+ 3 which describes two spins trying to align in the same direction zdirection and a term γτx which corresponds to a magnetic field acting in the x direction. So it s a D quantum spin chain with a magnetic field in the transverse direction. Our goal is now to study this H and find its ground and excited states. This will solve the D Ising model, at least in the regime where K v, K h. 3 Ordered and Disordered Regimes 3. Ordered regime If γ β, then the Hamiltonian is approximately neglect the γ term Then it has two ground states, 0 = N h x= H = β 0 x x=, or 0 = τ 3 xτ 3 x+. 3. N h x= 0 i. 3. The notations here correspond to a wave function a spinor for spins on sites x. These are spins all pointing up or down. In the spirit of spontaneous symmetry breaking the spins choose one of these two states and ignore the other. suppose the ground state corresponds to all spins up. The energy is E 0 = βn h. 3.3 An excited states would correspond to taking one spin and making it point down. This would correspond to the energy E = βn h + β. 3.4 Of course, there are many such states, since any spin out of N h spins, can be made to point down. But this does not matter. We see that E E 0 = β. Now since T = e H, this leads to the eigenvalues of the transfer matrix λ 0 = e E 0, λ = e E, 3.5 6

8 and the correlation length is It is finite. ξ = ln λ 0 λ = = E E 0 β Disordered regime Suppose now γ β. Neglect β to obtain The ground state takes the form H = γ 0 = N h x= x= τ x. 3.7 x. 3.8 This is because, is the eigenvector for τ with the eigenvalue. This corresponds to the energy E = γn h. 3.9 An excited state would correspond to taking one spin and replacing its spinor by, again, there are many such excited states. eigenvalue. So this gives This is an eigenvector of τ but with an E = γn h + γ. 3.0 So we see that ξ = = E E 0 γ. 3. In the next section, we are going to solve our problem exactly. We will see that the exact expression is ξ = γ β. 3. This matches 3.6 and 3. in the limits γ β and β γ respectively. 4 Jordan-Wigner transformation Now we will solve the Hamiltonian.7 exactly. First of all we do a unitary transformation so that H = γ x= τ 3 x β 7 x= τ xτ x+. 4.

9 Corresponding transformation is given by U = + iτ, and Then, it is convenient to introduce U τ 3 U = τ, U τ U = τ b x = τ x + iτx, b x = τ x iτx 4.3 They satisfy b x = b x = 0, b xb x + b x b x = 4.4 In this way they look like creation and annihilation operators. The transformations expressing τ in terms of b read τx = b x + b x, τx = ib x b x, τx 3 = i [ ] τ x, τx = bx b x b xb x. 4.5 Now we introduce the Jordan-Wigner string a x = e iπ y<x b yb y b x, a x = e iπ y<x b yb y b x. 4.6 Then not only a i obey creation and annihilation relations on a given site, but also a xa y + a y a x = 0, x y, 4.7 and so do a x, a y and a x, a y. Therefore, these are now fermions, called Jordan-Wigner fermions. Therefore, the Hamiltonian becomes taking into account the Jordan-Wigner string which leads to the minus sign in the term a x a x H = γ x= a x a x a x a x N h β x= a x a x a x+ + a x+ 4.8 This is the so-called Bogoliubov Hamiltonian, studied in superconductivity. This Hamiltonian can now be diagonalized using standard methods. 5 Diagonalizing the Hamiltonian - the straightforward method To diagonalize this Hamiltonian we go to the plane wave basis a x = Nh k e ikx a k, a x = Nh 8 k e ikx a k 5.

10 Substituting this into 4.8 gives H = γ k a k a k a k a k β k a k a k a k + a k e ik. 5. Here the sum over k goes from π to π, as usual this is called the Brillouin zone, in steps of π/n h. This Hamiltonian clearly splits into terms coupling the modes k and k for each positive k. Therefore, we can write it as H = H k, 5.3 k>0 where H k = γ β cos k a k a k a k a k H k can be written in a convenient matrix form iβ sin k a k a k + iβ sin k a ka k. 5.4 H k = a γ β cos k iβ sin k ak k a k iβ sin k γ + β cos k a. 5.5 k This matrix has the Bogoliubov-de-Gennes form from the theory of superconductivity. It is a little inconvenient that it has imaginary entries, but these imaginary entries can be completely eliminated if one does a transformation a k = iã k, a k = iã k. 5.6 The new variable ã k are as good a fermionic variable as the old one, and the Hamiltonian becomes H k = a γ β cos k β sin k ak k ã k β sin k γ + β cos k ã. 5.7 k Now the Hamiltonian can be diagonalized by a Bogoliubov transformation. This transformation reads a k = cosφ k c k + sinφ k c k, ã k = sinφ k c k + cosφ k c k. 5.8 Here k > 0. It is clear from this that φ k = φ k. Importantly, these transformations preserve the fermionic anticommutation relations. For example, { ak, a } { k = cosφk c k + sinφ k c k, sinφ k c k + cosφ k c } k = cos φ k + sin φ k =, 5.9 9

11 {a k, ã k } = { cosφ k c k + sinφ k c k, cosφ k c k sinφ k c k, } = cos φ k sin φ k +cos φ k sin φ k = 0, 5.0 and so on to compute the left hand side, we assume that c k satisfy the fermionic anticommutation relations, in other words, if c k are fermions, then so are a k. This Bogoliubov transformation is nothing but the orthogonal rotation of the matrix, in other words, ak cos φk sin φ ã = k ck k sin φ k cos φ k c k The Hamiltonian, when expressed in terms of c k, becomes H k = c γ β cos k β sin k k c k U ck U β sin k γ + β cos k c, U = k. 5. cos φk sin φ k. sin φ k cos φ k 5. We can now choose the orthogonal matrix U in such a way that the Hamiltonian becomes diagonal. Its eigenvalues are h k and h k where h k = γ β cos k + β sin k = γ + β γβ cos k = γ β + 4γβ sin k. This brings the Hamiltonian to the form Finally, this means H = k>0 5.3 H k = c hk 0 ck k c k 0 h k c. 5.4 k γ β + 4γβ sin k c k c k c k c k. 5.5 This is a very interesting Hamiltonian. In its ground state all fermionic states are empty. The ground state energy is then E 0 = k>0 γ β + 4γβ sin k = N h π π dk π γ β + 4γβ sin k. 5.6 To create an excitation, we populate one fermionic state, either with positive momentum q or with negative momentum q. That increases the energy by E q E 0 = γ β + 4γβ sin q. 5.7 The minimum energy increase occurs if we populate the state with q = 0. This gives E q=0 E 0 = γ β, 5.8 0

12 and leads to the previously advertised result for the correlation length. ξ = E q=0 E 0 = γ β. 5.9 The point where γ = β is the transition point. At this point the correlation length is infinity. Interestingly, at small q the single fermion excitation energy has a relativistic form E q E 0 /4 = γ β + γβq. 5.0 Here γ β plays the role of a relativistic mass. That s why the critical point γ = β is often called massless, while the phases with gaps in the spectrum are often called massive. Finally, we observe that this Hamiltonian 5.5 represents the Majorana fermions, not the usual Dirac fermions. That s because it represents the positive momentum particles whose energy increases with k, so they move to the right. But also it has negative momentum particles whose energy increases as k becomes more negative, so they move to the left. Real-world Dirac fermions would ve been moving either left and/or right at all k, not just at positive k. These kinds of special fermions which have their momentum restricted to either positive values c k moving to the right or negative values c k moving to the left are called Majorana fermions. 6 Diagonalizing the Hamiltonian - alternative scheme Here is another way to diagonalize the Hamiltonian, faster but involving a more convoluted mathematical construction. To diagonalize it, we employ the Bogoliubov transformations, c a A B a c = U a = a, c c = a a U. 6. B A the anticommutation relations. This is shorthand notations, in components this looks like this: U ij is a matrix, and N h [ ] c x = Axy a y + B xy a N h [ ] y, c x = B xy a y + A xya y. 6. y= Let us prove that U must be a unitary matrix. This is a consequence of the anticommutation relations. Let us introduce ψ j = y= a a. 6.3

13 This is again a shorthand notation. Here j goes from to N h. If j N h, then ψ j = a j. If N h < j N h, then ψ j = a j N h. We also introduce φ j = c c. 6.4 Now it should be clear that ψ jψ m + ψ m ψ j = δ ij. 6.5 Indeed, if m = j N h, this amounts to anticommutations between a m and a j. If m > N h, n > N h, this is again an anticommutator between a m N h, a n Nh. If m > N h and j < N h, this is an anticommutator between a m N h and a j, which is zero. Finally, if m N h, j > N h then this is an anticommutator between a m and a j Nh, which is also zero. In the same way, φ jφ m + φ m φ j = δ mj. 6.6 Now let us use that φ j = m U jm ψ m, φ j = m ψ mu mj. 6.7 Substitute these into 6.6, we find [ ψ n U nju ml ψ l + U ml ψ l ψ nu ] nj = δmj. 6.8 nl Now we note that ψ nψ l + ψ l ψ n = δ nl. 6.9 This gives This means UU = I, or U is a unitary matrix. Let us write the Hamiltonian 4.8 in the following way H = a a U mn U nj = δ mj. 6.0 n γ + β β β Γ Γ γ β a a = ψ Hψ. 6. Here H is a N h by N h matrix. We can now employ the Bogoliubov transformation, to write ψ = U φ. 6. Then H = φ UHU φ. 6.3

14 We choose U so that UHU is diagonal. Then in terms of c the Hamiltonian becomes H = H = k h k [ c k c k c k c k], 6.4 Here h k are eigenvalues of H. So the program is then to find the eigenvalues of H. We write γ + β H = where and Γ are matrices such that and Σ 3, Σ are Pauli matrices. Σ 3 + i β ΓΣ, 6.5 xy = δ x,y+ + δ x,y, Γ xy = δ x,y δ x,y+, 6.6 To find h k, we will use a series of unitary transformations. We first rotate Σ 3 into Σ by employing U = + iσ and H = U HU. We find H = γ + β Σ + i β ΓΣ 0 γ + β + Γ = γ + β. 6.7 Γ 0 The eigenvalue equations become These can be solved in terms of plane waves, hs x = γt x + βt x+, 6.8 ht x = γs x + βs x. 6.9 s x = s k e ikx, t x = t k e ikx. 6.0 Here x varies from to N h. k must vary from 0 to π because x is integer if k > π, then e ixk = e ixk π, so k larger than π is equivalent to k π. We now find This is equivalent to finding the eigenvalues of k = 0, π N h, 4π N h,..., π π N h. 6. hs k = γ + βe ik t k, ht k = γs k + βe ik s k γ + βe ik γ + βe ik

15 The eigenvalues are then h = γ + βe ik γ + βe ik = γ + β + γβ cosk = γ β + γβ cosk k hk = ± γ β + 4γβ cos. 6.5 It is convenient to introduce q = π k. As k varies from 0 to π, q varies from π to π. Then hq ± Thus the Hamiltonian becomes H = q γ β + 4γβ sin q γ β + 4γβ sin q Now we can construct the ground and excited states. fermions, then. 6.6 c qc q c q c q. 6.7 If we fill all the states with c qc q c q c q =. 6.8 This gives the minimal possible energy, that is the ground state. The ground state energy is thus this q, E 0 = q γ β + 4γβ sin q = N h π π dq q γ β π + 4γβ sin. 6.9 To create an excitation, we take a fermion with momentum q and remove it. Then for This increases energy by E q E 0 = c qc q c q c q = γ β + 4γβ sin q. 6.3 This gives excitation energy. It s also possible to remove more than one fermion, to get even higher energy. Interesting to note that at small q it has a relativistic form, This is the spectrum of particles of mass m = γ β. h q = γ β + γβq

16 7 Critical exponent The lowest excited state is at q = 0, giving E q=0 E 0 = γ β. 7. This is the gap in the spectrum. The correlation length is the inverse gap and is given by ξ = γ β. 7. It diverges at γ = β. It is customary, if the correlation length is divergent at some critical T c, to denote the power of its divergence as ν: ξ We see that for the D Ising model ν =. Indeed, recall that β = K h = J h T, γ = e Kv T T c ν. 7.3 = e Jv T. 7.4 So β γ = J h T Jv e T. 7.5 The critical temperature T c is when J h e Jv Tc = T c This is when the gap vanishes and the correlation length is infinity. If T is close to T c, then with some slope α. So β γ J h T c + J v e Jv Tc T T c = αt T c 7.7 T c ξ ν = is one of Onsager s famous exact results. T T c. 7.8 Another result for the free energy easily follows as well. The free energy is given by F = T log Z. 7.9 Since Z is proportional to exp N v E 0 where E 0 is the ground state energy and N v is the length of the Ising model in the time direction, that gives F = T N v E

17 In other words, the ground state energy in quantum mechanics is the same as the free energy in statistical mechanics. In this problem π dq q E 0 N h γ β π π + 4γβ sin 7. So π dq q F = T N v N h γ β π π + 4γβ sin. 7. Notice that F is proportional to N v N h, which is the total number of sites in our D Ising model, which is like a volume of the system. So F is extensive proportional to volume, as we should have expected. The integral in the expression for F is a known integral, whose value is related to elliptic functions. We are not interested in its exact value. Rather we are only interested in the specific heat c = T F dt 7.3 and only at T close to T c. To find it we need to differentiate F with respect to T, remembering that γ and β are both temperature dependent. We recall that at T close to T c, β γ αt T c. So we can write π dq q F = T N v N h α π π T T c + 4γT βt sin. 7.4 Differentiating with respect to T once gives F π T = T N dq α T T c + vn γβ T sin q h π π α T T c + 4γT βt sin q π dq q N v N h α π π T T c + 4γT βt sin. 7.5 Computing the second derivative gives F = T N vn π h dq α T T c + γβ T sin q T π π α T T c + 4γT βt sin 3 q π dq α + γβ sin q T T N v N h π π α T T c + 4γT βt sin q N v N h π π dq π α T T c + T γβ sin q α T T c + 4γT βt sin q

18 As T T c, the first and the third term remain finite. The second term can be estimated to behave as log T T c. So we find c log. 7.7 T T c This is another of Onsager s famous results. A standard notation in statistical mechanics is From here it follows that α = 0. c T T c α 7.8 7

Magnets, 1D quantum system, and quantum Phase transitions

Magnets, 1D quantum system, and quantum Phase transitions 134 Phys620.nb 10 Magnets, 1D quantum system, and quantum Phase transitions In 1D, fermions can be mapped into bosons, and vice versa. 10.1. magnetization and frustrated magnets (in any dimensions) Consider

More information

Introduction to the Mathematics of the XY -Spin Chain

Introduction to the Mathematics of the XY -Spin Chain Introduction to the Mathematics of the XY -Spin Chain Günter Stolz June 9, 2014 Abstract In the following we present an introduction to the mathematical theory of the XY spin chain. The importance of this

More information

The 1+1-dimensional Ising model

The 1+1-dimensional Ising model Chapter 4 The 1+1-dimensional Ising model The 1+1-dimensional Ising model is one of the most important models in statistical mechanics. It is an interacting system, and behaves accordingly. Yet for a variety

More information

Lecture notes on topological insulators

Lecture notes on topological insulators Lecture notes on topological insulators Ming-Che Chang Department of Physics, National Taiwan Normal University, Taipei, Taiwan Dated: May 8, 07 I. D p-wave SUPERCONDUCTOR Here we study p-wave SC in D

More information

Physics 505 Homework No. 8 Solutions S Spinor rotations. Somewhat based on a problem in Schwabl.

Physics 505 Homework No. 8 Solutions S Spinor rotations. Somewhat based on a problem in Schwabl. Physics 505 Homework No 8 s S8- Spinor rotations Somewhat based on a problem in Schwabl a) Suppose n is a unit vector We are interested in n σ Show that n σ) = I where I is the identity matrix We will

More information

Physics 239/139 Spring 2018 Assignment 2 Solutions

Physics 239/139 Spring 2018 Assignment 2 Solutions University of California at San Diego Department of Physics Prof. John McGreevy Physics 39/139 Spring 018 Assignment Solutions Due 1:30pm Monday, April 16, 018 1. Classical circuits brain-warmer. (a) Show

More information

Chapter 17 The bilinear covariant fields of the Dirac electron. from my book: Understanding Relativistic Quantum Field Theory.

Chapter 17 The bilinear covariant fields of the Dirac electron. from my book: Understanding Relativistic Quantum Field Theory. Chapter 17 The bilinear covariant fields of the Dirac electron from my book: Understanding Relativistic Quantum Field Theory Hans de Vries November 10, 008 Chapter Contents 17 The bilinear covariant fields

More information

ANTICOMMUTING INTEGRALS AND FERMIONIC FIELD THEORIES FOR TWO-DIMENSIONAL ISING MODELS. V.N. Plechko

ANTICOMMUTING INTEGRALS AND FERMIONIC FIELD THEORIES FOR TWO-DIMENSIONAL ISING MODELS. V.N. Plechko ANTICOMMUTING INTEGRALS AND FERMIONIC FIELD THEORIES FOR TWO-DIMENSIONAL ISING MODELS V.N. Plechko Bogoliubov Theoretical Laboratory, Joint Institute for Nuclear Research, JINR-Dubna, 141980 Dubna, Moscow

More information

Tutorial 5 Clifford Algebra and so(n)

Tutorial 5 Clifford Algebra and so(n) Tutorial 5 Clifford Algebra and so(n) 1 Definition of Clifford Algebra A set of N Hermitian matrices γ 1, γ,..., γ N obeying the anti-commutator γ i, γ j } = δ ij I (1) is the basis for an algebra called

More information

QM and Angular Momentum

QM and Angular Momentum Chapter 5 QM and Angular Momentum 5. Angular Momentum Operators In your Introductory Quantum Mechanics (QM) course you learned about the basic properties of low spin systems. Here we want to review that

More information

Physics 127b: Statistical Mechanics. Renormalization Group: 1d Ising Model. Perturbation expansion

Physics 127b: Statistical Mechanics. Renormalization Group: 1d Ising Model. Perturbation expansion Physics 17b: Statistical Mechanics Renormalization Group: 1d Ising Model The ReNormalization Group (RNG) gives an understanding of scaling and universality, and provides various approximation schemes to

More information

8.334: Statistical Mechanics II Problem Set # 4 Due: 4/9/14 Transfer Matrices & Position space renormalization

8.334: Statistical Mechanics II Problem Set # 4 Due: 4/9/14 Transfer Matrices & Position space renormalization 8.334: Statistical Mechanics II Problem Set # 4 Due: 4/9/14 Transfer Matrices & Position space renormalization This problem set is partly intended to introduce the transfer matrix method, which is used

More information

4. Two-level systems. 4.1 Generalities

4. Two-level systems. 4.1 Generalities 4. Two-level systems 4.1 Generalities 4. Rotations and angular momentum 4..1 Classical rotations 4.. QM angular momentum as generator of rotations 4..3 Example of Two-Level System: Neutron Interferometry

More information

G : Quantum Mechanics II

G : Quantum Mechanics II G5.666: Quantum Mechanics II Notes for Lecture 5 I. REPRESENTING STATES IN THE FULL HILBERT SPACE Given a representation of the states that span the spin Hilbert space, we now need to consider the problem

More information

The quantum state as a vector

The quantum state as a vector The quantum state as a vector February 6, 27 Wave mechanics In our review of the development of wave mechanics, we have established several basic properties of the quantum description of nature:. A particle

More information

129 Lecture Notes More on Dirac Equation

129 Lecture Notes More on Dirac Equation 19 Lecture Notes More on Dirac Equation 1 Ultra-relativistic Limit We have solved the Diraction in the Lecture Notes on Relativistic Quantum Mechanics, and saw that the upper lower two components are large

More information

1 Infinite-Dimensional Vector Spaces

1 Infinite-Dimensional Vector Spaces Theoretical Physics Notes 4: Linear Operators In this installment of the notes, we move from linear operators in a finitedimensional vector space (which can be represented as matrices) to linear operators

More information

Quantum Mechanics II

Quantum Mechanics II Quantum Mechanics II Prof. Boris Altshuler March 8, 011 1 Lecture 19 1.1 Second Quantization Recall our results from simple harmonic oscillator. We know the Hamiltonian very well so no need to repeat here.

More information

Problem set for the course Skálázás és renormálás a statisztikus fizikában, 2014

Problem set for the course Skálázás és renormálás a statisztikus fizikában, 2014 1 Problem set for the course Skálázás és renormálás a statisztikus fizikában, 014 Rules: You can choose at wish from problems having the same main number (i.e. from a given section), but you can collect

More information

Maxwell s equations. electric field charge density. current density

Maxwell s equations. electric field charge density. current density Maxwell s equations based on S-54 Our next task is to find a quantum field theory description of spin-1 particles, e.g. photons. Classical electrodynamics is governed by Maxwell s equations: electric field

More information

VI.D Self Duality in the Two Dimensional Ising Model

VI.D Self Duality in the Two Dimensional Ising Model VI.D Self Duality in the Two Dimensional Ising Model Kramers and Wannier discovered a hidden symmetry that relates the properties of the Ising model on the square lattice at low and high temperatures.

More information

Linear algebra. 1.1 Numbers. d n x = 10 n (1.1) n=m x

Linear algebra. 1.1 Numbers. d n x = 10 n (1.1) n=m x 1 Linear algebra 1.1 Numbers The natural numbers are the positive integers and zero. Rational numbers are ratios of integers. Irrational numbers have decimal digits d n d n x = 10 n (1.1 n=m x that do

More information

conventions and notation

conventions and notation Ph95a lecture notes, //0 The Bloch Equations A quick review of spin- conventions and notation The quantum state of a spin- particle is represented by a vector in a two-dimensional complex Hilbert space

More information

in-medium pair wave functions the Cooper pair wave function the superconducting order parameter anomalous averages of the field operators

in-medium pair wave functions the Cooper pair wave function the superconducting order parameter anomalous averages of the field operators (by A. A. Shanenko) in-medium wave functions in-medium pair-wave functions and spatial pair particle correlations momentum condensation and ODLRO (off-diagonal long range order) U(1) symmetry breaking

More information

Plan for the rest of the semester. ψ a

Plan for the rest of the semester. ψ a Plan for the rest of the semester ϕ ψ a ϕ(x) e iα(x) ϕ(x) 167 Representations of Lorentz Group based on S-33 We defined a unitary operator that implemented a Lorentz transformation on a scalar field: and

More information

Attempts at relativistic QM

Attempts at relativistic QM Attempts at relativistic QM based on S-1 A proper description of particle physics should incorporate both quantum mechanics and special relativity. However historically combining quantum mechanics and

More information

Second quantization. Emmanuel Fromager

Second quantization. Emmanuel Fromager Institut de Chimie, Strasbourg, France Page 1 Emmanuel Fromager Institut de Chimie de Strasbourg - Laboratoire de Chimie Quantique - Université de Strasbourg /CNRS M2 lecture, Strasbourg, France. Institut

More information

where a is the lattice constant of the triangular Bravais lattice. reciprocal space is spanned by

where a is the lattice constant of the triangular Bravais lattice. reciprocal space is spanned by Contents 5 Topological States of Matter 1 5.1 Intro.......................................... 1 5.2 Integer Quantum Hall Effect..................... 1 5.3 Graphene......................................

More information

Lecture 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 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 information

Representations of Lorentz Group

Representations of Lorentz Group Representations of Lorentz Group based on S-33 We defined a unitary operator that implemented a Lorentz transformation on a scalar field: How do we find the smallest (irreducible) representations of the

More information

Potentially useful reading Sakurai and Napolitano, sections 1.5 (Rotation), Schumacher & Westmoreland chapter 2

Potentially useful reading Sakurai and Napolitano, sections 1.5 (Rotation), Schumacher & Westmoreland chapter 2 Problem Set 2: Interferometers & Random Matrices Graduate Quantum I Physics 6572 James Sethna Due Friday Sept. 5 Last correction at August 28, 2014, 8:32 pm Potentially useful reading Sakurai and Napolitano,

More information

H&M Chapter 5 Review of Dirac Equation

H&M Chapter 5 Review of Dirac Equation HM Chapter 5 Review of Dirac Equation Dirac s Quandary Notation Reminder Dirac Equation for free particle Mostly an exercise in notation Define currents Make a complete list of all possible currents Aside

More information

Electrons in a periodic potential

Electrons 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 information

Problem Set No. 3: Canonical Quantization Due Date: Wednesday October 19, 2018, 5:00 pm. 1 Spin waves in a quantum Heisenberg antiferromagnet

Problem Set No. 3: Canonical Quantization Due Date: Wednesday October 19, 2018, 5:00 pm. 1 Spin waves in a quantum Heisenberg antiferromagnet Physics 58, Fall Semester 018 Professor Eduardo Fradkin Problem Set No. 3: Canonical Quantization Due Date: Wednesday October 19, 018, 5:00 pm 1 Spin waves in a quantum Heisenberg antiferromagnet In this

More information

Integrable structure of various melting crystal models

Integrable structure of various melting crystal models Integrable structure of various melting crystal models Kanehisa Takasaki, Kinki University Taipei, April 10 12, 2015 Contents 1. Ordinary melting crystal model 2. Modified melting crystal model 3. Orbifold

More information

Lecture 4 - Relativistic wave equations. Relativistic wave equations must satisfy several general postulates. These are;

Lecture 4 - Relativistic wave equations. Relativistic wave equations must satisfy several general postulates. These are; Lecture 4 - Relativistic wave equations Postulates Relativistic wave equations must satisfy several general postulates. These are;. The equation is developed for a field amplitude function, ψ 2. The normal

More information

C/CS/Phy191 Problem Set 6 Solutions 3/23/05

C/CS/Phy191 Problem Set 6 Solutions 3/23/05 C/CS/Phy191 Problem Set 6 Solutions 3/3/05 1. Using the standard basis (i.e. 0 and 1, eigenstates of Ŝ z, calculate the eigenvalues and eigenvectors associated with measuring the component of spin along

More information

The Klein-Gordon Equation Meets the Cauchy Horizon

The Klein-Gordon Equation Meets the Cauchy Horizon Enrico Fermi Institute and Department of Physics University of Chicago University of Mississippi May 10, 2005 Relativistic Wave Equations At the present time, our best theory for describing nature is Quantum

More information

SSH Model. Alessandro David. November 3, 2016

SSH 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 information

Quantum Physics 2006/07

Quantum Physics 2006/07 Quantum Physics 6/7 Lecture 7: More on the Dirac Equation In the last lecture we showed that the Dirac equation for a free particle i h t ψr, t = i hc α + β mc ψr, t has plane wave solutions ψr, t = exp

More information

PY 351 Modern Physics - Lecture notes, 3

PY 351 Modern Physics - Lecture notes, 3 PY 351 Modern Physics - Lecture notes, 3 Copyright by Claudio Rebbi, Boston University, October 2016. These notes cannot be duplicated and distributed without explicit permission of the author. Time dependence

More information

Physics 110. Electricity and Magnetism. Professor Dine. Spring, Handout: Vectors and Tensors: Everything You Need to Know

Physics 110. Electricity and Magnetism. Professor Dine. Spring, Handout: Vectors and Tensors: Everything You Need to Know Physics 110. Electricity and Magnetism. Professor Dine Spring, 2008. Handout: Vectors and Tensors: Everything You Need to Know What makes E&M hard, more than anything else, is the problem that the electric

More information

3 Symmetry Protected Topological Phase

3 Symmetry Protected Topological Phase Physics 3b Lecture 16 Caltech, 05/30/18 3 Symmetry Protected Topological Phase 3.1 Breakdown of noninteracting SPT phases with interaction Building on our previous discussion of the Majorana chain and

More information

Enhancing Superconductivity by Disorder

Enhancing Superconductivity by Disorder UNIVERSITY OF COPENHAGEN FACULTY OF SCIENCE Enhancing Superconductivity by Disorder Written by Marie Ernø-Møller 16.01.19 Supervised by Brian Møller Andersen Abstract In this thesis an s-wave superconductor

More information

Relativistic Dirac fermions on one-dimensional lattice

Relativistic Dirac fermions on one-dimensional lattice Niodem Szpa DUE, 211-1-2 Relativistic Dirac fermions on one-dimensional lattice Niodem Szpa Universität Duisburg-Essen & Ralf Schützhold Plan: 2 Jan 211 Discretized relativistic Dirac fermions (in an external

More information

The 3 dimensional Schrödinger Equation

The 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 information

Systems of Linear ODEs

Systems of Linear ODEs P a g e 1 Systems of Linear ODEs Systems of ordinary differential equations can be solved in much the same way as discrete dynamical systems if the differential equations are linear. We will focus here

More information

Kitaev honeycomb lattice model: from A to B and beyond

Kitaev 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 information

Quantum Mechanics C (130C) Winter 2014 Assignment 7

Quantum Mechanics C (130C) Winter 2014 Assignment 7 University of California at San Diego Department of Physics Prof. John McGreevy Quantum Mechanics C (130C) Winter 014 Assignment 7 Posted March 3, 014 Due 11am Thursday March 13, 014 This is the last problem

More information

Section 10: Many Particle Quantum Mechanics Solutions

Section 10: Many Particle Quantum Mechanics Solutions Physics 143a: Quantum Mechanics I Section 10: Many Particle Quantum Mechanics Solutions Spring 015, Harvard Here is a summary of the most important points from this week (with a few of my own tidbits),

More information

7.1 Creation and annihilation operators

7.1 Creation and annihilation operators Chapter 7 Second Quantization Creation and annihilation operators. Occupation number. Anticommutation relations. Normal product. Wick s theorem. One-body operator in second quantization. Hartree- Fock

More information

E 2 = p 2 + m 2. [J i, J j ] = iɛ ijk J k

E 2 = p 2 + m 2. [J i, J j ] = iɛ ijk J k 3.1. KLEIN GORDON April 17, 2015 Lecture XXXI Relativsitic Quantum Mechanics 3.1 Klein Gordon Before we get to the Dirac equation, let s consider the most straightforward derivation of a relativistically

More information

Topological insulator with time-reversal symmetry

Topological 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 information

Classical field theory 2012 (NS-364B) Feynman propagator

Classical field theory 2012 (NS-364B) Feynman propagator Classical field theory 212 (NS-364B Feynman propagator 1. Introduction States in quantum mechanics in Schrödinger picture evolve as ( Ψt = Û(t,t Ψt, Û(t,t = T exp ı t dt Ĥ(t, (1 t where Û(t,t denotes the

More information

A. Time dependence from separation of timescales

A. Time dependence from separation of timescales Lecture 4 Adiabatic Theorem So far we have considered time independent semiclassical problems. What makes these semiclassical is that the gradient term (which is multiplied by 2 ) was small. In other words,

More information

d 3 r d 3 vf( r, v) = N (2) = CV C = n where n N/V is the total number of molecules per unit volume. Hence e βmv2 /2 d 3 rd 3 v (5)

d 3 r d 3 vf( r, v) = N (2) = CV C = n where n N/V is the total number of molecules per unit volume. Hence e βmv2 /2 d 3 rd 3 v (5) LECTURE 12 Maxwell Velocity Distribution Suppose we have a dilute gas of molecules, each with mass m. If the gas is dilute enough, we can ignore the interactions between the molecules and the energy will

More information

PHYS 508 (2015-1) Final Exam January 27, Wednesday.

PHYS 508 (2015-1) Final Exam January 27, Wednesday. PHYS 508 (2015-1) Final Exam January 27, Wednesday. Q1. Scattering with identical particles The quantum statistics have some interesting consequences for the scattering of identical particles. This is

More information

Topological invariants for 1-dimensional superconductors

Topological invariants for 1-dimensional superconductors Topological invariants for 1-dimensional superconductors Eddy Ardonne Jan Budich 1306.4459 1308.soon SPORE 13 2013-07-31 Intro: Transverse field Ising model H TFI = L 1 i=0 hσ z i + σ x i σ x i+1 σ s:

More information

Time dependent perturbation theory 1 D. E. Soper 2 University of Oregon 11 May 2012

Time dependent perturbation theory 1 D. E. Soper 2 University of Oregon 11 May 2012 Time dependent perturbation theory D. E. Soper University of Oregon May 0 offer here some background for Chapter 5 of J. J. Sakurai, Modern Quantum Mechanics. The problem Let the hamiltonian for a system

More information

Second Quantization: Quantum Fields

Second Quantization: Quantum Fields Second Quantization: Quantum Fields Bosons and Fermions Let X j stand for the coordinate and spin subscript (if any) of the j-th particle, so that the vector of state Ψ of N particles has the form Ψ Ψ(X

More information

Topological Insulator Surface States and Electrical Transport. Alexander Pearce Intro to Topological Insulators: Week 11 February 2, / 21

Topological Insulator Surface States and Electrical Transport. Alexander Pearce Intro to Topological Insulators: Week 11 February 2, / 21 Topological Insulator Surface States and Electrical Transport Alexander Pearce Intro to Topological Insulators: Week 11 February 2, 2017 1 / 21 This notes are predominately based on: J.K. Asbóth, L. Oroszlány

More information

(1.1) In particular, ψ( q 1, m 1 ; ; q N, m N ) 2 is the probability to find the first particle

(1.1) In particular, ψ( q 1, m 1 ; ; q N, m N ) 2 is the probability to find the first particle Chapter 1 Identical particles 1.1 Distinguishable particles The Hilbert space of N has to be a subspace H = N n=1h n. Observables Ân of the n-th particle are self-adjoint operators of the form 1 1 1 1

More information

MODEL WITH SPIN; CHARGE AND SPIN EXCITATIONS 57

MODEL WITH SPIN; CHARGE AND SPIN EXCITATIONS 57 56 BOSONIZATION Note that there seems to be some arbitrariness in the above expressions in terms of the bosonic fields since by anticommuting two fermionic fields one can introduce a minus sine and thus

More information

3 Quantization of the Dirac equation

3 Quantization of the Dirac equation 3 Quantization of the Dirac equation 3.1 Identical particles As is well known, quantum mechanics implies that no measurement can be performed to distinguish particles in the same quantum state. Elementary

More information

Linearization of Differential Equation Models

Linearization of Differential Equation Models Linearization of Differential Equation Models 1 Motivation We cannot solve most nonlinear models, so we often instead try to get an overall feel for the way the model behaves: we sometimes talk about looking

More information

16.1. PROBLEM SET I 197

16.1. PROBLEM SET I 197 6.. PROBLEM SET I 97 Answers: Problem set I. a In one dimension, the current operator is specified by ĵ = m ψ ˆpψ + ψˆpψ. Applied to the left hand side of the system outside the region of the potential,

More information

Exact diagonalization methods

Exact diagonalization methods Summer School on Computational Statistical Physics August 4-11, 2010, NCCU, Taipei, Taiwan Exact diagonalization methods Anders W. Sandvik, Boston University Representation of states in the computer bit

More information

Donoghue, Golowich, Holstein Chapter 4, 6

Donoghue, Golowich, Holstein Chapter 4, 6 1 Week 7: Non linear sigma models and pion lagrangians Reading material from the books Burgess-Moore, Chapter 9.3 Donoghue, Golowich, Holstein Chapter 4, 6 Weinberg, Chap. 19 1 Goldstone boson lagrangians

More information

Introduction to Group Theory

Introduction to Group Theory Chapter 10 Introduction to Group Theory Since symmetries described by groups play such an important role in modern physics, we will take a little time to introduce the basic structure (as seen by a physicist)

More information

VI.D Self Duality in the Two Dimensional Ising Model

VI.D Self Duality in the Two Dimensional Ising Model VI.D Self Duality in the Two Dimensional Ising Model Kramers and Wannier discovered a hidden symmetry that relates the properties of the Ising model on the square lattice at low and high temperatures.

More information

1 Mathematical preliminaries

1 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 information

Dirac Equation. Chapter 1

Dirac Equation. Chapter 1 Chapter Dirac Equation This course will be devoted principally to an exposition of the dynamics of Abelian and non-abelian gauge theories. These form the basis of the Standard Model, that is, the theory

More information

III. Particle Physics and Isospin

III. Particle Physics and Isospin . Particle Physics and sospin Up to now we have concentrated exclusively on actual, physical rotations, either of coordinates or of spin states, or both. n this chapter we will be concentrating on internal

More information

Physics 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 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 information

Quantum 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. 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 information

Construction of spinors in various dimensions

Construction of spinors in various dimensions Construction of spinors in various dimensions Rhys Davies November 23 2011 These notes grew out of a desire to have a nice Majorana representation of the gamma matrices in eight Euclidean dimensions I

More information

221B Lecture Notes Spontaneous Symmetry Breaking

221B Lecture Notes Spontaneous Symmetry Breaking B Lecture Notes Spontaneous Symmetry Breaking Spontaneous Symmetry Breaking Spontaneous Symmetry Breaking is an ubiquitous concept in modern physics, especially in condensed matter and particle physics.

More information

Symmetries, Groups, and Conservation Laws

Symmetries, Groups, and Conservation Laws Chapter Symmetries, Groups, and Conservation Laws The dynamical properties and interactions of a system of particles and fields are derived from the principle of least action, where the action is a 4-dimensional

More information

Summary of free theory: one particle state: vacuum state is annihilated by all a s: then, one particle state has normalization:

Summary of free theory: one particle state: vacuum state is annihilated by all a s: then, one particle state has normalization: The LSZ reduction formula based on S-5 In order to describe scattering experiments we need to construct appropriate initial and final states and calculate scattering amplitude. Summary of free theory:

More information

x(τ) y(τ) = x(τ) f(τ), where f(τ) is the solution of the variational equation δs(f) = 0 f(τ) = 0

x(τ) y(τ) = x(τ) f(τ), where f(τ) is the solution of the variational equation δs(f) = 0 f(τ) = 0 Calculation of the path integral We have introduced the formal concept of path integral. The importance of the concept is that path integrals can often be calculated without referring to a limiting process,

More information

Chapter 2. Linear Algebra. rather simple and learning them will eventually allow us to explain the strange results of

Chapter 2. Linear Algebra. rather simple and learning them will eventually allow us to explain the strange results of Chapter 2 Linear Algebra In this chapter, we study the formal structure that provides the background for quantum mechanics. The basic ideas of the mathematical machinery, linear algebra, are rather simple

More information

On the Random XY Spin Chain

On the Random XY Spin Chain 1 CBMS: B ham, AL June 17, 2014 On the Random XY Spin Chain Robert Sims University of Arizona 2 The Isotropic XY-Spin Chain Fix a real-valued sequence {ν j } j 1 and for each integer n 1, set acting on

More information

C/CS/Phys C191 Particle-in-a-box, Spin 10/02/08 Fall 2008 Lecture 11

C/CS/Phys C191 Particle-in-a-box, Spin 10/02/08 Fall 2008 Lecture 11 C/CS/Phys C191 Particle-in-a-box, Spin 10/0/08 Fall 008 Lecture 11 Last time we saw that the time dependent Schr. eqn. can be decomposed into two equations, one in time (t) and one in space (x): space

More information

LS coupling. 2 2 n + H s o + H h f + H B. (1) 2m

LS coupling. 2 2 n + H s o + H h f + H B. (1) 2m LS coupling 1 The big picture We start from the Hamiltonian of an atomic system: H = [ ] 2 2 n Ze2 1 + 1 e 2 1 + H s o + H h f + H B. (1) 2m n e 4πɛ 0 r n 2 4πɛ 0 r nm n,m Here n runs pver the electrons,

More information

lattice QCD and the hadron spectrum Jozef Dudek ODU/JLab

lattice QCD and the hadron spectrum Jozef Dudek ODU/JLab lattice QCD and the hadron spectrum Jozef Dudek ODU/JLab the light meson spectrum relatively simple models of hadrons: bound states of constituent quarks and antiquarks the quark model empirical meson

More information

We already came across a form of indistinguishably in the canonical partition function: V N Q =

We already came across a form of indistinguishably in the canonical partition function: V N Q = Bosons en fermions Indistinguishability We already came across a form of indistinguishably in the canonical partition function: for distinguishable particles Q = Λ 3N βe p r, r 2,..., r N ))dτ dτ 2...

More information

S j H o = gµ o H o. j=1

S j H o = gµ o H o. j=1 LECTURE 17 Ferromagnetism (Refs.: Sections 10.6-10.7 of Reif; Book by J. S. Smart, Effective Field Theories of Magnetism) Consider a solid consisting of N identical atoms arranged in a regular lattice.

More information

Mean-Field Theory: HF and BCS

Mean-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 information

Chapter 16 The relativistic Dirac equation. from my book: Understanding Relativistic Quantum Field Theory. Hans de Vries

Chapter 16 The relativistic Dirac equation. from my book: Understanding Relativistic Quantum Field Theory. Hans de Vries Chapter 6 The relativistic Dirac equation from my book: Understanding Relativistic Quantum Field Theory Hans de Vries November, 28 2 Chapter Contents 6 The relativistic Dirac equation 6. The linearized

More information

GROUP THEORY PRIMER. New terms: so(2n), so(2n+1), symplectic algebra sp(2n)

GROUP THEORY PRIMER. New terms: so(2n), so(2n+1), symplectic algebra sp(2n) GROUP THEORY PRIMER New terms: so(2n), so(2n+1), symplectic algebra sp(2n) 1. Some examples of semi-simple Lie algebras In the previous chapter, we developed the idea of understanding semi-simple Lie algebras

More information

Potts And XY, Together At Last

Potts And XY, Together At Last Potts And XY, Together At Last Daniel Kolodrubetz Massachusetts Institute of Technology, Center for Theoretical Physics (Dated: May 16, 212) We investigate the behavior of an XY model coupled multiplicatively

More information

Chm 331 Fall 2015, Exercise Set 4 NMR Review Problems

Chm 331 Fall 2015, Exercise Set 4 NMR Review Problems Chm 331 Fall 015, Exercise Set 4 NMR Review Problems Mr. Linck Version.0. Compiled December 1, 015 at 11:04:44 4.1 Diagonal Matrix Elements for the nmr H 0 Find the diagonal matrix elements for H 0 (the

More information

SUSY Breaking in Gauge Theories

SUSY Breaking in Gauge Theories SUSY Breaking in Gauge Theories Joshua Berger With the Witten index constraint on SUSY breaking having been introduced in last week s Journal club, we proceed to explicitly determine the constraints on

More information

Physics 218. Quantum Field Theory. Professor Dine. Green s Functions and S Matrices from the Operator (Hamiltonian) Viewpoint

Physics 218. Quantum Field Theory. Professor Dine. Green s Functions and S Matrices from the Operator (Hamiltonian) Viewpoint Physics 28. Quantum Field Theory. Professor Dine Green s Functions and S Matrices from the Operator (Hamiltonian) Viewpoint Field Theory in a Box Consider a real scalar field, with lagrangian L = 2 ( µφ)

More information

Elementary Linear Algebra

Elementary Linear Algebra Matrices J MUSCAT Elementary Linear Algebra Matrices Definition Dr J Muscat 2002 A matrix is a rectangular array of numbers, arranged in rows and columns a a 2 a 3 a n a 2 a 22 a 23 a 2n A = a m a mn We

More information

Graduate Quantum Mechanics I: Prelims and Solutions (Fall 2015)

Graduate Quantum Mechanics I: Prelims and Solutions (Fall 2015) Graduate Quantum Mechanics I: Prelims and Solutions (Fall 015 Problem 1 (0 points Suppose A and B are two two-level systems represented by the Pauli-matrices σx A,B σ x = ( 0 1 ;σ 1 0 y = ( ( 0 i 1 0 ;σ

More information

Adiabatic particle pumping and anomalous velocity

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 information

Spins and spin-orbit coupling in semiconductors, metals, and nanostructures

Spins and spin-orbit coupling in semiconductors, metals, and nanostructures B. Halperin Spin lecture 1 Spins and spin-orbit coupling in semiconductors, metals, and nanostructures Behavior of non-equilibrium spin populations. Spin relaxation and spin transport. How does one produce

More information

Lecture 10: Solving the Time-Independent Schrödinger Equation. 1 Stationary States 1. 2 Solving for Energy Eigenstates 3

Lecture 10: Solving the Time-Independent Schrödinger Equation. 1 Stationary States 1. 2 Solving for Energy Eigenstates 3 Contents Lecture 1: Solving the Time-Independent Schrödinger Equation B. Zwiebach March 14, 16 1 Stationary States 1 Solving for Energy Eigenstates 3 3 Free particle on a circle. 6 1 Stationary States

More information

Lecture 4 - Dirac Spinors

Lecture 4 - Dirac Spinors Lecture 4 - Dirac Spinors Schrödinger & Klein-Gordon Equations Dirac Equation Gamma & Pauli spin matrices Solutions of Dirac Equation Fermion & Antifermion states Left and Right-handedness Non-Relativistic

More information