Part VI, Chapter 1 | Page 2 of 3

Worked Examples

Anharmonic oscillator, Stark effect, and Zeeman effect in detail

Now we apply non-degenerate perturbation theory to several classic problems that illustrate the power and limitations of the method. Each example demonstrates different aspects of the formalism developed on the previous page.

Example 1: Anharmonic Oscillator

Consider a quantum harmonic oscillator perturbed by a quartic anharmonic term. The full Hamiltonian is:

$$\hat{H} = \underbrace{\frac{\hat{p}^2}{2m} + \frac{1}{2}m\omega^2\hat{x}^2}_{\hat{H}_0} + \underbrace{\lambda\hat{x}^4}_{\hat{H}'}$$

The unperturbed system is the harmonic oscillator with well-known energy levels $E_n^{(0)} = \hbar\omega(n + \tfrac{1}{2})$ and ladder operator algebra.

Step 1: Express in Ladder Operators

The position operator in terms of ladder operators is:

$$\hat{x} = \sqrt{\frac{\hbar}{2m\omega}}(\hat{a} + \hat{a}^\dagger)$$

Therefore the perturbation becomes:

$$\hat{H}' = \lambda\hat{x}^4 = \lambda\left(\frac{\hbar}{2m\omega}\right)^2 (\hat{a} + \hat{a}^\dagger)^4$$

Step 2: First-Order Energy Correction to Ground State

We need $E_0^{(1)} = \langle 0|\hat{H}'|0\rangle = \lambda(\hbar/2m\omega)^2\langle 0|(\hat{a}+\hat{a}^\dagger)^4|0\rangle$.

Expanding $(\hat{a}+\hat{a}^\dagger)^4$, only terms with equal numbers of $\hat{a}$ and $\hat{a}^\dagger$ survive (to return to $|0\rangle$). Using normal ordering and the commutation relation $[\hat{a},\hat{a}^\dagger]=1$:

$$\langle 0|(\hat{a}+\hat{a}^\dagger)^4|0\rangle = \langle 0|\hat{a}^2(\hat{a}^\dagger)^2 + \hat{a}\hat{a}^\dagger\hat{a}\hat{a}^\dagger + \hat{a}(\hat{a}^\dagger)^2\hat{a} + \hat{a}^\dagger\hat{a}^2\hat{a}^\dagger + \hat{a}^\dagger\hat{a}\hat{a}^\dagger\hat{a} + (\hat{a}^\dagger)^2\hat{a}^2|0\rangle$$

Evaluating each term carefully using $\hat{a}|0\rangle = 0$ and $\hat{a}^\dagger|n\rangle = \sqrt{n+1}|n+1\rangle$:

$$\langle 0|(\hat{a}+\hat{a}^\dagger)^4|0\rangle = 3$$

Therefore the first-order energy correction is:

$$\boxed{E_0^{(1)} = \frac{3\lambda\hbar^2}{4m^2\omega^2}}$$

The anharmonic term raises the ground state energy, consistent with the positive quartic potential adding confinement.

Step 3: Second-Order Energy Correction

For the second-order correction, we need matrix elements $\langle m|\hat{x}^4|0\rangle$ for $m \neq 0$. The operator $\hat{x}^4$ can connect $|0\rangle$ to $|2\rangle$ and $|4\rangle$ only (by selection rules from ladder operators):

$$\langle 2|(\hat{a}+\hat{a}^\dagger)^4|0\rangle = 6\sqrt{2}, \quad \langle 4|(\hat{a}+\hat{a}^\dagger)^4|0\rangle = 2\sqrt{6}$$

Applying the second-order formula:

$$E_0^{(2)} = \lambda^2\left(\frac{\hbar}{2m\omega}\right)^4\left[\frac{|\langle 2|(\hat{a}+\hat{a}^\dagger)^4|0\rangle|^2}{E_0^{(0)}-E_2^{(0)}} + \frac{|\langle 4|(\hat{a}+\hat{a}^\dagger)^4|0\rangle|^2}{E_0^{(0)}-E_4^{(0)}}\right]$$
$$E_0^{(2)} = \lambda^2\left(\frac{\hbar}{2m\omega}\right)^4\left[\frac{72}{-2\hbar\omega} + \frac{24}{-4\hbar\omega}\right] = -\frac{21\lambda^2\hbar^3}{8m^4\omega^5}$$

As expected, the second-order correction is negative (lowering the energy) and proportional to $\lambda^2$.

Example 2: Stark Effect (Non-Degenerate Levels)

The Stark effect describes the splitting of atomic energy levels in an external uniform electric field $\vec{\mathcal{E}} = \mathcal{E}_0\hat{z}$. The perturbation Hamiltonian is:

$$\hat{H}' = e\mathcal{E}_0 z = e\mathcal{E}_0 r\cos\theta$$

For hydrogen, the first-order correction vanishes for all states with definite parity (since $z$ is odd under parity, $\langle nlm|z|nlm\rangle = 0$ for all states with well-defined $l$). This is the quadratic Stark effect:

$$E_n^{(2)} = e^2\mathcal{E}_0^2 \sum_{n'\ell' m' \neq n\ell m} \frac{|\langle n'\ell' m'|r\cos\theta|n\ell m\rangle|^2}{E_n - E_{n'}}$$

For the hydrogen ground state (1s), the exact second-order result is:

$$E_1^{(2)} = -\frac{9}{4}a_0^3 \mathcal{E}_0^2 = -\frac{1}{2}\alpha_d \mathcal{E}_0^2$$

where $\alpha_d = \frac{9}{2}a_0^3$ is the static electric polarizability of hydrogen. This energy shift is always negative and quadratic in the field strength.

Important Distinction

For non-degenerate levels, the Stark effect is quadratic in the electric field ($\Delta E \propto \mathcal{E}_0^2$). For degenerate levels like $n=2$ in hydrogen, there is a linear Stark effect ($\Delta E \propto \mathcal{E}_0$) because states with opposite parity ($|2s\rangle$ and $|2p\rangle$) can mix. The linear effect is much larger and requires degenerate perturbation theory, covered in Chapter 2.

Example 3: Weak-Field Zeeman Effect

An atom placed in a weak magnetic field $\vec{B} = B\hat{z}$ experiences the perturbation:

$$\hat{H}' = -\hat{\vec{\mu}}\cdot\vec{B} = \frac{eB}{2m_e}(\hat{L}_z + 2\hat{S}_z)$$

In the coupled basis $|j, m_j\rangle$, the first-order energy correction involves the Lande g-factor:

$$E^{(1)} = g_J \mu_B B m_j$$

where the Lande g-factor is:

$$g_J = 1 + \frac{j(j+1) + s(s+1) - l(l+1)}{2j(j+1)}$$

and $\mu_B = e\hbar/(2m_e) \approx 5.79 \times 10^{-5}$ eV/T is the Bohr magneton. Each $j$-level splits into $2j+1$ equally spaced sublevels.

Worked Example: Perturbed Infinite Square Well

Problem: A particle in an infinite square well (width $a$) is perturbed by $\hat{H}' = V_0\sin(\pi x/a)$. Find the second-order energy correction to the ground state.

Step 1: Unperturbed ground state

$$\psi_1^{(0)}(x) = \sqrt{\frac{2}{a}}\sin\left(\frac{\pi x}{a}\right), \quad E_1^{(0)} = \frac{\pi^2\hbar^2}{2ma^2}$$

Step 2: First-order correction

$$E_1^{(1)} = \frac{2V_0}{a}\int_0^a \sin^2\!\left(\frac{\pi x}{a}\right)\sin\!\left(\frac{\pi x}{a}\right)dx = V_0 \cdot \frac{8}{3\pi} \approx 0.849\, V_0$$

Step 3: Matrix elements for second order

Using $\sin^2\theta = (1-\cos 2\theta)/2$ and product-to-sum formulas, the only non-vanishing matrix element is with $n=3$:

$$\langle 3|\hat{H}'|1\rangle = \frac{V_0}{2}$$

Step 4: Second-order formula

$$E_1^{(2)} = \frac{|\langle 3|\hat{H}'|1\rangle|^2}{E_1^{(0)} - E_3^{(0)}} = \frac{(V_0/2)^2}{\frac{\pi^2\hbar^2}{2ma^2}(1 - 9)} = -\frac{V_0^2 ma^2}{32\pi^2\hbar^2}$$

The second-order correction is negative and scales as $V_0^2$, as expected from perturbation theory.

Convergence and Breakdown

The convergence of the perturbation series depends critically on the ratio of perturbation matrix elements to energy level spacings. Several important scenarios arise:

  • Good convergence: When $|\langle m|\hat{H}'|n\rangle| \ll |E_n^{(0)} - E_m^{(0)}|$ for all states, low-order perturbation theory gives excellent results.
  • Asymptotic series: The perturbation series for the anharmonic oscillator ($\lambda x^4$) actually diverges for any $\lambda > 0$! The series is asymptotic: the first few terms approximate well, but adding more terms eventually makes the result worse. The reason is that for $\lambda < 0$, the potential is unbounded below, so the series has zero radius of convergence.
  • Near-degeneracy: When two levels are close in energy, the corresponding terms in the sum become very large. The cure is to treat the near-degenerate subspace exactly (quasi-degenerate perturbation theory).

Practical Rule of Thumb

A useful criterion for when to trust perturbation theory to order $k$: the $(k+1)$-th order correction should be significantly smaller than the $k$-th order correction. If you compute:

$$\frac{|E_n^{(k+1)}|}{|E_n^{(k)}|} \lesssim 0.1$$

then the perturbation series is converging well at order $k$, and the truncation error is roughly the size of the last term kept.

Applications in Physics

The examples on this page illustrate the broad applicability of non-degenerate perturbation theory:

  • Anharmonic oscillator: Models vibrations of molecules beyond harmonic approximation
  • Quadratic Stark effect: Gives atomic polarizabilities, critical for understanding van der Waals forces
  • Zeeman effect: Foundation of magnetic resonance spectroscopy (NMR, ESR) and astrophysical spectral analysis
  • Crystal field theory: Explains optical properties of transition metal complexes
  • Lamb shift: First evidence for quantum electrodynamic corrections in hydrogen