Second Quantization
Total Page:16
File Type:pdf, Size:1020Kb
Physics 127c: Statistical Mechanics Second Quantization Ladder Operators in the SHO It is useful to first review the use of ladder operators in the simple harmonic oscillator. Here I present the bare bones—review your favorite Quantum textbook for more details. The Hamiltonian for the SHO is 1 1 H = p2 + Kx2. (1) 2m 2 For this section only, I will scale variables so that m = 1,K = 1,and also h¯ = 1, to save writing. So the Hamiltonian is 1 1 H = p + x2 (2) 2 2 and using p =−i∂/∂x (remember h¯ = 1) Schrodinger’s equation Hψ = Eψ gives eigenstates ψn and energy eigenvalues εn d2 1 1 2 − ψn(x) + x ψn(x) = εnψ(x). (3) 2 dx2 2 This differential equation can be solved with all the technology of orthogonal polynomials to find the eigen- 1 εn = n + values 2 , etc. We can instead solve the problem algebraically by introducing up and down ladder operators b+ and b 1 b = √ (p − ix) (4) 2 1 b+ = (b)+ = √ (p + ix) (5) 2 where the + superscript denotes the Hermitian conjugate. We know that the commutator [x,p] ≡ [xp − px] = i (6) is equivalent to the coordinate representation p →−i∂/∂x, and so it is natural to look at the commutation properties of b, b+: 1 [b, b+] = [p − ix,p + ix] (7) 2 =−i[x,p] = 1. (8) The Hamiltonian is 1 H = b+b + , (9) 2 and we are led to look at [b, b+b] = bb+b − b+bb (10) = [b, b+]b = b, (11) and [b+,b+b] = b+b+b − b+bb+ (12) =−b+[b, b+] =−b+. (13) 1 Collecting the results: [b, b+] = 1, (14a) [b, b+b] = b, (14b) [b+,b+b] =−b+. (14c) Since H is the sum of squares, the eigenvalues must be positive, and so there is lowest energy ε0 and corresponding normalized eigenfunction ψ0. Let’s consider the state bψ0, and ask how the Hamiltonian acts on this state using the commutation rules to transfer the action of H onto ψ0: 1 H(bψ ) = (b+b + )bψ (15) 0 2 0 1 = [b(b+b + ) − b]ψ (16) 2 0 = (ε0 − 1)(bψ0). (17) This suggests that bψ0 is an eigenfunction with energy ε0 − 1, but since ε0 was defined as the lowest energy this is not possible, and the only way out is bψ0 = 0. (18) Thus the ground state eigenfunction is defined by the action of the down ladder operator b giving zero. Similarly + + H(b ψ0) = (ε0 + 1)(b ψ0) (19) + so that b ψ0 is an eigenfunction with energy ε0 + 1. Continuing in this way we find a “ladder” of eigenstates + n 1 ψn = (b ) ψ with energy εn = n + . (20) 0 2 We can argue that this generates the complete set of states, etc., etc. (see your Quantum text book). Note that + b bψn = nψn (21) so the operator b+b acts as a level counter. We have not addressed the normalization of the ψn created in this way. The normalization of the ψn can be related: + + hψn|ψni = ψn−1|bb |ψn−1 = ψn−1|b b + 1|ψn−1 = n hψn−1|ψn−1i . (22) Since ψ0 is normalized, this gives hψn|ψni = n!. Then the functions −1/2 un = (n!) ψn (23) are normalized eigenfunctions, hun|uni = 1. Rewriting the effect of the ladder operators on ψn in terms of the normalized states un gives √ + b un = n + 1un+ (24) √ 1 bun = nun−1 (25) + b bun = nun. (26) Remembering that the level n can be thought of as n bosons in the single particle state of energy 1 we therefore have the following dictionary. We have introduced a raising operator or boson creation operator b+, and a lowering operator or boson annihilation operator b. The operator b+b = n counts the level or the 2 number of bosons. The Hamiltonian is just b+b plus a constant. (Going back to unscaled units, it would be b+b times the harmonic oscillator energy hω¯ , plus a constant.) The lowering or annihilation operator gives zero when acting on the lowest or zero particle state. The analysis in terms of ode theory has been replaced by an analysis in terms of the algebra of the b, b+ operators, a considerable simplification. Other expectation values can also be calculated e.g. 1 u |x2|u =− u |(b+ − b)(b+ − b)|u (27) 0 0 2 0 0 1 =− u |(b+b+ − b+b − bb+ + bb|u (28) 2 0 0 1 1 = u |(bb+|u = (29) 2 0 0 2 + where b |u0i = 0 = hu0| b has been used. Again this is considerably easier than integrating over the Gaussian form of u0(x). Many Particle Systems Occupation number representation For a many particle system we choose as the basis any symmetrized (bosons) or antisymmetrized (fermions) product of complete orthonormal set of single particle states φi. Often these will be chosen as momentum states, position states, or eigenstates of the noninteracting Hamiltonian. Bosons: The N particle basis state is X ψ (rE , rE ...rE ) = C φ (P rE )φ (P rE )...φ (P rE )...φ (P rE ) B 1 2 N 1 1 1 2 1 n1 Q N (30) P where φq ,q = 1 ...Q occurs nq times (“nq particles in the state φq ”). The sum runs over N! permutations P and the normalization constant is −1/2 C = (N!n1!n2! ...nQ!) . (31) Last term we introduced the occupation number representation where this state is represented ψB = n1n2 ...nQ ... (32) + where all the nq are zero for the empty states q>Q. We define boson creation operator bq and annihilation operator bq that adds or subtracts one particle from the single particle state φq + 1/2 b n n ...nq ...nQ = (nq + 1) n n ...(nq + 1)...nQ , (33a) q 1 2 1 2 1/2 bq n1n2 ...nq ...nQ = (nq ) n1n2 ...(nq − 1)...nQ . (33b) The prefactors are motivated by the ladder-operator results, and can be checked to preserve the nor- malization. Of course if nq = 0, the second result is just zero, consistent with the action of b on a zero particle state. It is easy to check that b, b+ satisfy + [br ,bs ] = δrs, (34) + + [br ,bs] = [br ,bs ] = 0. (35) 3 Fermions The antisymmetrized Fermi basis state is the Slater determinant X P ψF (rE1, rE2 ...rEN ) = C (−1) φ1(P rE1)φ2(P rE2)...φN (P rEN ) (36) P where now all the φi have to be different for a nonzero wave function (equivalent to the Pauli exclusion principle). The normalization constant is C = (N!)−1/2. (37) The occupation number representation is ψF = |n1n2 ...i (38) where each ni is one if the state appears in the Slater determinant and zero otherwise. Note that the order of the states affects the overall sign through the (−1)P , so we must have some convention for ordering the list of one particle states. What are the appropriate creation and annihilation operators? First consider the single particle state. Possible occupation numbers are 0, 1, i.e. the states |0i and |1i. Define creation operator a+ and annihilation operator a such that a+ |0i = |1i a+ |1i = 0 (39) a |0i = 0 a |1i = |0i . (Take care of the difference between 0 and |0i.) It is easy to check that a,a+ satisfy anticommutation rules a,a+ = 1 (40a) {a,a} = 0 = a+,a+ (40b) where {A, B} = AB + BA. (41) For a multiparticle system S (−1) |n1n2 ...0s ...i if ns = 1 as |n1n2 ...ns ...i = (42a) 0ifns = 0 + 0ifns = 1 as |n1n2 ...ns ...i = S (42b) (−1) |n1n2 ...1s ...i if ns = 0 where S = n1 + n2 +···ns−1 (43) is to keep track of the ordering convention. The full algebra is then a ,a+ = δ , r s rs (44a) + + {ar ,as} = 0 = ar ,as . (44b) Note, for example that a+a+ | , ,...i =−a+a+ | , ,...i 1 2 0 0 2 1 0 0 (45) consistent with the definition in terms of Slater determinants. 4 In either case the number operator that counts the number of particles is the same + bs bs BOSONS ns = + (46) as as FERMIONS and has eigenvalues corresponding to the possible occupation numbers of the state s (0, 1 for Fermions, 0, 1,...for Bosons). It follows from the commutation or anticommutation rules that in either case [ns,nr ] = 0 (47) so that the set {nr } are simultaneously diagonalizable, as we expect from the definition of occupation number states. A general many particle wave function ψ can be expressed as a linear combination of the occupation number basis states. The expectation value hψ |ns| ψi may then be nonintegral. We have replaced complicated symmetric product or determinant wave functions by states in the simpler occupation number representation, and have developed the algebra of operators connecting these states. Any occupation number state can be generated by acting on the no-particle state |0i (the “vacuum state”) with appropriate combinations of creation operators c+. This notation is called “second quantization” as opposed to the original “first quantization” notation. It is the notation of quantum field theory. The next task is to write physical operators in terms of these elementary operators. Physical Operators We are usually interested in few-body operators, i.e. those involving a fewP particle coordinates in the T =−h2 ∇2 “usual”P first quantization notation, such as the kinetic energy 2m i i or the potential energy V = u(rE −Er ) i<j i j . P For a one body operator O1 = fi the second quantized expression is X + O1 = hφr |f | φsi cr cs (48) r,s where hφr |f | φsi is the matrix element Z 3 ∗ hφr |f | φsi = d rφr (r)f(E r)φE s(r)E (49) and c represents a for Fermions andPb for Bosons.