G Quantum Mechanical Perturbation Theory
Quantum mechanical perturbation theory is a widely used method in solid- state physics. Without the details of derivation, we shall list a number of basic formulas of time-independent (stationary) and time-dependent perturbation theory below. For simplicity, we shall use the Dirac notation for wavefunctions and matrix elements.
G.1 Time-Independent Perturbation Theory
Assume that the complete solution (eigenfunctions and eigenvalues) of the Schrödinger equation " " H (0) (0) (0) 0 ψi = Ei ψi (G.1.1) is known for a system described by a simple Hamiltonian H0. If the system is subject to a time-independent (stationary) perturbation described by the Hamiltonian H1 – which can be an external perturbation or the interaction between the components of the system –, the eigenvalues and eigenfunctions change. The method for determining the new ones depends on whether the unperturbed energy level in question is degenerate or not.
G.1.1 Nondegenerate Perturbation Theory
We now introduce a fictitious coupling constant λ, whose value will be treated as a parameter in the calculations and set equal to unity in the final result, and write the full Hamiltonian H = H0 + H1 as
H = H0 + λH1 . (G.1.2)
The parameter λ is purely a bookkeeping device to keep track of the relative order of magnitude of the various terms, since the energy eigenvalues and eigenfunctions will be sought in the form of an expansion in powers of λ: 580 G Quantum Mechanical Perturbation Theory " " ∞ " (0) n (n) ψi = ψi + λ ψi , n=1 ∞ (G.1.3) (0) n (n) Ei = Ei + λ Ei . n=1 The series is convergent if the perturbation is weak, that is, in addition to the formally introduced parameter λ, the interaction Hamiltonian itself contains a small parameter, the physical coupling constant. By substituting this expansion into the Schrödinger equation and collect- ing the same powers of λ from both sides, we obtain " " (0) (0) (0) H0 ψ = E ψ , " i " i i " " (1) (0) (0) (1) (1) (0) H0 ψ + H1 ψ = E ψ + E ψ , (G.1.4) i " i " i i " i i " " H (2) H (1) (0) (2) (1) (1) (2) (0) 0 ψi + 1 ψi = Ei ψi + Ei ψi + Ei ψi and similar equations for higher-order corrections. The corrections to the en- ergy and wavefunction of any order are related to the lower-order ones by the recursion formula " " H − (0) (n) H − (1) (n−1) ( 0 Ei ) ψi +( 1 Ei ) ψi " " (G.1.5) − (2) (n−2) − − (n) (0) Ei ψi ... Ei ψi =0.
Multiplying the second! equation in (G.1.4) (which comes from the terms (0) that are linear in λ)by ψi from the left, the first-order correction to the energy is ! " (1) (0)H (0) Ei = ψi 1 ψi . (G.1.6) To determine! the correction to the wavefunction, the same equation is multi- (0) plied by ψj (j = i): ! " ! ! " (0) (0) (1) (0)H (0) (0) (0) (1) Ej ψj ψi + ψj 1 ψi = Ei ψj ψi . (G.1.7) " H (n) Since the eigenfunctions of 0 make up a complete set, the functions ψi can be expanded in terms of them: " " (n) (n) (0) ψi = Cij ψj . (G.1.8) j
(n) The coefficients Cii are not determined by the previous equations: their val- ues depend on the normalization of the perturbed wavefunction. Substituting the previous formula into (G.1.7), we have ! (0) (n) (0)H (0) (0) (n) Ej Cij + ψj 1 ψi = Ei Cij , (G.1.9) G.1 Time-Independent Perturbation Theory 581 and hence ! " (0) (0) (1)" ψ H1 ψ (0)" ψ = j i ψ . (G.1.10) i (0) − (0) j j =i Ei Ej The second-order correction to the energy is then
! "2 (0) (0) H1 (2) ψj ψi E = , (G.1.11) i (0) − (0) j =i Ei Ej and the second-order correction to the wavefunction is ! "! " (0) (0) (0) (0) (2)" ψ H1 ψ ψ H1 ψ (0) ψ = j k k i ψ i (0) − (0) (0) − (0) j j =i k =i Ei Ej Ei Ek ! "! " (G.1.12) (0) (0) (0) (0) ψ H1 ψ ψ H1 ψ (0) − j i i i ψ . (0) − (0) 2 j j =i Ei Ej
However, this wavefunction is not normalized to unity. Proper normalization is ensured by the choice ! "! " (0) (0) (0) (0) (2)" ψ H1 ψ ψ H1 ψ (0)" ψ = j k k i ψ i (0) − (0) (0) − (0) j j =i k =i Ei Ej Ei Ek ! "! " (0)H (0) (0)H (0) " ψj 1 ψi ψi 1 ψi (0) − ψ (G.1.13) (0) − (0) 2 j j =i Ei Ej ! "! " (0) (0) (0) (0) ψ H1 ψ ψ H1 ψ (0)" − 1 i j j i ψ . 2 (0) − (0) 2 i j =i Ei Ej
Finally, the third-order correction to the energy is ! "! "! " (0) (0) (0) (0) (0) (0) (3) ψ H1 ψ ψ H1 ψ ψ H1 ψ E = i j j k k i i (0) − (0) (0) − (0) j =i k =i Ei Ej Ei Ek ! "! "! " (G.1.14) (0) (0) (0) (0) (0) (0) ψ H1 ψ ψ H1 ψ ψ H1 ψ − i j j i i i . (0) − (0) 2 j =i Ei Ej
In this Rayleigh–Schrödinger perturbation theory the explicit form of higher-order corrections becomes increasingly complicated. A relatively simple recursion formula can be obtained by introducing the projection operator "! "! − (0) (0) (0) (0) Pi =1 ψi ψi = ψj ψj (G.1.15) j =i 582 G Quantum Mechanical Perturbation Theory " (0) which projects onto the subspace that is orthogonal to the state ψi .The nth-order energy correction can then be written as ! " (n) (0)H (n−1) Ei = ψi 1 ψi , (G.1.16) where the matrix element is to be taken with the wavefunction " # " (n) 1 H − (1) (n−1) ψi = (0) Pi 1 Ei ψi E −H0 i " "$ (G.1.17) − (2) (n−2) − − (n−1) (1) Ei ψi ... Ei ψi , which is in the subspace mentioned above. Formally simpler expressions can be obtained when the Brillouin–Wigner perturbation theory is used. The perturbed wavefunction in the Schrödinger equation " " H0 + H1 ψi = Ei ψi (G.1.18) is then chosen in the form " " " (0) ψi = C0 ψi + Δψi , (G.1.19) " " (0) where Δψi is orthogonal to ψi , and C0 takes care of the appropriate normalization. After some algebra, the eigenvalue equation reads " " " H − H − (0) (0) 0 Ei Δψi + 1 ψi = C0 Ei Ei ψi . (G.1.20)
By applying the projection operator Pi, and exploiting the relations " " " (0) Pi ψi =0,Pi ψi = Δψi (G.1.21) as well as the commutation of Pi and H0, " " Ei −H0 Δψi = PiH1 ψi (G.1.22) is obtained. Its formal solution is " " " (0) Pi H ψi = C0 ψi + 1 ψi . (G.1.23) Ei −H0 Iteration then yields " ∞ n " Pi H (0) ψi = C0 1 ψi , (G.1.24) E −H0 n=0 i and ∞ ! n " (0)H Pi H (0) ΔEi = ψi 1 1 ψi (G.1.25) E −H0 n=0 i G.1 Time-Independent Perturbation Theory 583 for the energy correction. In this method the energy denominator contains the (0) perturbed energy Ei rather than the unperturbed one Ei . To first order in the interaction, ! " (0) (0)H (0) Ei = Ei + ψi 1 ψi , (G.1.26) while to second order, ! " (0) (0) (0) Ei = E + ψ H1 ψ i ! i "!i " (0) (0) (0) (0) ψ H1 ψ ψ H1 ψ (G.1.27) + i j j i + ... . − (0) j =i Ei Ej
It is easy to show that by rearranging the energy denominator and expanding it as ∞ n 1 1 1 −ΔEi = (0) = (0) (0) , (G.1.28) Ei −H0 −H −H −H Ei 0 +ΔEi Ei 0 n=0 Ei 0 the results of the Rayleigh–Schrödinger perturbation theory are recovered. The formulas of time-dependent perturbation theory can also be used to determine the ground-state energy and wavefunction of the perturbed sys- tem, provided the interaction is assumed to be turned on adiabatically. The appropriate formulas are given in Section G.2.
G.1.2 Degenerate Perturbation Theory
In the previous subsection we studied the shift of nondegenerate energy levels due to the perturbation. For degenerate levels a slightly different method has to be used because the formal application of the previous formulas would yield vanishing energy denominators. Assuming that the ith energy level of the unperturbed system is p-fold (0) (0)" degenerate – that is, the same energy E belongs to each of the states ψ , " " i i1 (0) (0) ψi2 , ..., ψip –, any linear combination of these degenerate eigenstates is also an eigenstate of H0 with the same energy. We shall use such linear combinations to determine the perturbed states. We write the wavefunctions of the states of the perturbed system that arise from the degenerate states as " " " (0) (0) ψ = cik ψik + cn ψn , (G.1.29) k n =i
where the cik are of order unity, whereas the other coefficients cn that specify the mixing with the unperturbed eigenstates whose energy is different from (0) Ei are small, proportional to the perturbation. By substituting! this form (0) into the Schrödinger equation, and multiplying both sides by ψij from the left, 584 G Quantum Mechanical Perturbation Theory ! " ! " (0)H (0) (0)H (0) ΔEcij = ψij 1 ψik cik + ψij 1 ψn cn (G.1.30) k n =i is obtained. Since the coefficients cn are small, the second term on the right- hand side can be neglected in calculating the leading-order energy correction, which is given by #! " $ (0)H (0) − ψij 1 ψik δjkΔE cik =0. (G.1.31) k This homogeneous system of equations has nontrivial solutions if the deter- minant of the coefficient matrix vanishes: ! " (0)H (0) − det ψij 1 ψik δjkΔE =0. (G.1.32)
The solutions of this pth-order equation! – that" is, the eigenvalues of the matrix (0)H (0) made up of the matrix elements ψij 1 ψik – specify the eventual splitting of the initially p-fold degenerate level, i.e., the shift of the perturbed levels with respect to the unperturbed one. Thus the interaction Hamiltonian needs to be diagonalized on the subspace of the degenerate states of H0. In general, the degeneracy is lifted at least partially by the perturbation. As discussed in Appendix D on group theory, the symmetry properties of the full Hamiltonian determine which irreducible representations appear, and what the degree of degeneracy is for each new level.
G.2 Time-Dependent Perturbation Theory
If the perturbation depends explicitly on time, no stationary states can arise. We may then be interested in the evolution of" the system: What states can (0) be reached at time t from an initial state ψi if the perturbation is turned on suddenly at time t0? The answer lies in the solution of the time-dependent Schrödinger equation " " ∂ H0 + λH1(t) ψ (t) = − ψ (t) . (G.2.1) i i ∂t i " The wavefunction ψi(t) is sought in the form
" (0)" (0) −iEj t/ ψi(t) = cij (t) ψj e , (G.2.2) j subject to the initial condition
cij (t0)=δij . (G.2.3) Since the time dependence of the unperturbed state has been written out explicitly, the functions cij (t) are expected to vary slowly in time. Expanding the coefficients once again into powers of λ, G.2 Time-Dependent Perturbation Theory 585 ∞ (0) r (r) cij (t)=cij (t)+ λ cij (t) , (G.2.4) r=1 where, naturally, the zeroth-order term is a constant:
(0) cij (t)=δij . (G.2.5) Substituting this series expansion into the Schrödinger equation, we find
(0) (0) ! " ∂ (r) i(E −E )t/ (0) (0) (r−1) − c (t)= e j k ψ H1(t) ψ c (t) . (G.2.6) i ∂t ij j k ik k The explicit formulas for the first two terms obtained through iteration are