Skip to content

Nondegenerate Perturbation Theory

Nondegenerate time-independent perturbation theory computes how an isolated eigenvalue and its eigenstate change when a Hamiltonian is changed slightly.

Write

H(λ)=H0+λV,H(\lambda) = H_0+\lambda V,

where H0H_0 is exactly solved and λ\lambda is a bookkeeping parameter. The unperturbed eigenvalue equation is

H0∣n(0)⟩=En(0)∣n(0)⟩.H_0\lvert n^{(0)}\rangle = E_n^{(0)} \lvert n^{(0)}\rangle.

The method assumes the level En(0)E_n^{(0)} is isolated: no other unperturbed state has the same energy, and the perturbation does not mix it too strongly with nearby levels.

Seek an eigenvalue and eigenstate of

H(λ)∣n(λ)⟩=En(λ)∣n(λ)⟩H(\lambda)\lvert n(\lambda)\rangle = E_n(\lambda)\lvert n(\lambda)\rangle

as formal expansions:

En(λ)=En(0)+λEn(1)+λ2En(2)+⋯ ,E_n(\lambda) = E_n^{(0)} + \lambda E_n^{(1)} + \lambda^2E_n^{(2)} + \cdots,

and

∣n(λ)⟩=∣n(0)⟩+λ∣n(1)⟩+λ2∣n(2)⟩+⋯ .\lvert n(\lambda)\rangle = \lvert n^{(0)}\rangle + \lambda\lvert n^{(1)}\rangle + \lambda^2\lvert n^{(2)}\rangle + \cdots.

The physical perturbation may be obtained by setting λ=1\lambda=1 after the order counting is clear, but this is safe only if the corresponding dimensionless ratios are small.

A common convention is intermediate normalization:

⟨n(0)∣n(λ)⟩=1.\langle n^{(0)}|n(\lambda)\rangle=1.

This implies

⟨n(0)∣n(1)⟩=0.\langle n^{(0)}|n^{(1)}\rangle=0.

This condition fixes the component of ∣n(1)⟩\lvert n^{(1)}\rangle parallel to ∣n(0)⟩\lvert n^{(0)}\rangle. Other normalization conventions are possible. Energy corrections do not depend on this convention, but state-correction formulas do.

Substituting the expansions into the eigenvalue equation and matching powers of λ\lambda gives, at first order,

(H0−En(0))∣n(1)⟩=(En(1)−V)∣n(0)⟩.(H_0-E_n^{(0)})\lvert n^{(1)}\rangle = (E_n^{(1)}-V) \lvert n^{(0)}\rangle.

Projecting with ⟨n(0)∣\langle n^{(0)}| gives the first-order energy correction:

En(1)=⟨n(0)∣V∣n(0)⟩.E_n^{(1)} = \langle n^{(0)}|V|n^{(0)}\rangle.

Using the shorthand

Vmn=⟨m(0)∣V∣n(0)⟩,V_{mn} = \langle m^{(0)}|V|n^{(0)}\rangle,

this is

En(1)=Vnn.E_n^{(1)}=V_{nn}.

For symmetry tests, sign checks, parameter-response interpretation, and focused examples, see First-Order Energy Corrections.

Project the first-order equation onto ⟨m(0)∣\langle m^{(0)}| with m≠nm\ne n:

(Em(0)−En(0))⟨m(0)∣n(1)⟩=−Vmn.(E_m^{(0)}-E_n^{(0)}) \langle m^{(0)}|n^{(1)}\rangle = -V_{mn}.

Thus

∣n(1)⟩=∑m≠nVmnEn(0)−Em(0)∣m(0)⟩.\lvert n^{(1)}\rangle = \sum_{m\ne n} \frac{V_{mn}} {E_n^{(0)}-E_m^{(0)}} \lvert m^{(0)}\rangle.

This formula is the first place energy-gap denominators appear. Small denominators are a warning sign, not a detail to ignore.

For normalization and phase conventions, observable corrections, selection-rule diagnostics, exact checks, and the geometric interpretation of this vector, see First-Order State Corrections.

The second-order energy correction is

En(2)=∑m≠n∣Vmn∣2En(0)−Em(0).E_n^{(2)} = \sum_{m\ne n} \frac{\lvert V_{mn}\rvert^2} {E_n^{(0)}-E_m^{(0)}}.

This formula is one of the main workhorses of bound-state perturbation theory. Its sign depends on whether the coupled states lie above or below the state being corrected.

For a ground state, all other unperturbed energies satisfy Em(0)>E0(0)E_m^{(0)}\gt E_0^{(0)}, so every term in E0(2)E_0^{(2)} is nonpositive.

For the derivation, reduced-resolvent form, sign theorem, static-polarizability connection, continuum terms, and convergence diagnostics, see Second-Order Energy Corrections.

For the all-order recursion, third-order subtraction structure, convergence radius, asymptotic truncation, and the switch to subspace diagonalization, see Higher-Order Structure.

For the historical identity, normalization conventions, repartitioning, comparison with Brillouin–Wigner theory, and atomic and molecular uses of the order-by-order method, see Rayleigh–Schrödinger Perturbation Theory.

A useful diagnostic is

ϵmn=∣λVmnEn(0)−Em(0)∣.\epsilon_{mn} = \left| \frac{\lambda V_{mn}} {E_n^{(0)}-E_m^{(0)}} \right|.

For nondegenerate perturbation theory to be reliable, the relevant ϵmn\epsilon_{mn} should be much smaller than 11. This is not a rigorous universal bound, but it is the basic physical test.

If the perturbation has nonzero matrix elements connecting nearly degenerate levels, diagonalize the coupled subspace instead.

Consider

H=(E1λvλv∗E2),H = \begin{pmatrix} E_1&\lambda v\\ \lambda v^*&E_2 \end{pmatrix},

with E1<E2E_1\lt E_2 and Δ=E2−E1>0\Delta=E_2-E_1\gt 0. For the lower level, V11=0V_{11}=0 and

E1(2)=∣v∣2E1−E2=−∣v∣2Δ.E_1^{(2)} = \frac{\lvert v\rvert^2}{E_1-E_2} = -\frac{\lvert v\rvert^2}{\Delta}.

The exact lower eigenvalue is

E1+E22−(Δ2)2+λ2∣v∣2.\frac{E_1+E_2}{2} - \sqrt{ \left(\frac{\Delta}{2}\right)^2 + \lambda^2\lvert v\rvert^2 }.

Expanding for λ∣v∣≪Δ\lambda\lvert v\rvert\ll\Delta gives

E1−λ2∣v∣2Δ+O(λ4),E_1 - \lambda^2 \frac{\lvert v\rvert^2}{\Delta} + O(\lambda^4),

matching the perturbative result.

Nondegenerate perturbation theory is excellent for small shifts of isolated levels, weak mixing between well-separated states, and expectation-value corrections when state corrections are included consistently.

It is not automatically reliable for wavefunction details, near-degenerate levels, tunneling splittings, or observables sensitive to high-energy tails.

  • Using the formula when Em(0)=En(0)E_m^{(0)}=E_n^{(0)} or nearly so.
  • Forgetting that the component of ∣n(1)⟩\lvert n^{(1)}\rangle along ∣n(0)⟩\lvert n^{(0)}\rangle depends on the normalization convention.
  • Treating state corrections themselves as observables.
  • Setting λ=1\lambda=1 before identifying the actual small dimensionless ratio.
  • Ignoring symmetry, which can force VnnV_{nn} or entire blocks of VV to vanish.
  • R. Shankar, Principles of Quantum Mechanics, 2nd ed., Springer, 1994.
  • J. J. Sakurai and J. Napolitano, Modern Quantum Mechanics, 3rd ed., Cambridge University Press, 2020.
  • C. Cohen-Tannoudji, B. Diu, and F. Laloë, Quantum Mechanics, Wiley, 1977.
  • A. Messiah, Quantum Mechanics, Dover, 1999.
  1. Show that if Vnn=0V_{nn}=0, the first nonzero energy correction may be second order.
Solution

The first-order correction is

En(1)=Vnn.E_n^{(1)}=V_{nn}.

If Vnn=0V_{nn}=0, this term vanishes. The second-order correction is

En(2)=∑m≠n∣Vmn∣2En(0)−Em(0).E_n^{(2)} = \sum_{m\ne n} \frac{\lvert V_{mn}\rvert^2} {E_n^{(0)}-E_m^{(0)}}.

This can be nonzero if the perturbation couples ∣n(0)⟩\lvert n^{(0)}\rangle to other unperturbed states.

  1. For the two-level Hamiltonian above, explain why the perturbation expansion fails when Δ→0\Delta\to0.
Solution

The second-order correction contains

−∣v∣2Δ.-\frac{\lvert v\rvert^2}{\Delta}.

As Δ→0\Delta\to0, this correction diverges, indicating that the assumed unperturbed basis is no longer a good starting point. One must diagonalize the coupled two-state subspace, which is degenerate perturbation theory in the Δ=0\Delta=0 limit.