|q⟩ Bad Qubits

intermediate · Physics · Time-Independent Perturbation Theory

Second-Order Energy Corrections

Going to second order

The first-order energy used only the unperturbed states. The second-order correction is the first place the state correction feeds back into the energy. Return to the order-by-order expansion and collect the λ2\lambda^2 terms of the eigenvalue equation:

H^(0)n(2)+H^n(1)=En(0)n(2)+En(1)n(1)+En(2)n(0).\hat{H}^{(0)}|n^{(2)}\rangle + \hat{H}'|n^{(1)}\rangle = E_n^{(0)}|n^{(2)}\rangle + E_n^{(1)}|n^{(1)}\rangle + E_n^{(2)}|n^{(0)}\rangle.

Project onto n(0)\langle n^{(0)}|. As before the H^(0)\hat{H}^{(0)} term cancels, and with the normalization choice n(0)n(1)=0\langle n^{(0)}|n^{(1)}\rangle = 0 the En(1)E_n^{(1)} term drops too. What survives is

En(2)=n(0)H^n(1).E_n^{(2)} = \langle n^{(0)}|\hat{H}'|n^{(1)}\rangle.

Now substitute the explicit first-order state from the previous lesson, n(1)=mnm(0)H^n(0)En(0)Em(0)m(0)|n^{(1)}\rangle = \sum_{m\neq n} \frac{\langle m^{(0)}|\hat{H}'|n^{(0)}\rangle}{E_n^{(0)}-E_m^{(0)}}|m^{(0)}\rangle. Using n(0)H^m(0)=m(0)H^n(0)\langle n^{(0)}|\hat{H}'|m^{(0)}\rangle = \langle m^{(0)}|\hat{H}'|n^{(0)}\rangle^* (Hermiticity), the numerator becomes a squared modulus:

En(2)=mnm(0)H^n(0)2En(0)Em(0).\boxed{\,E_n^{(2)} = \sum_{m \neq n} \frac{\big|\langle m^{(0)}|\hat{H}'|n^{(0)}\rangle\big|^2}{E_n^{(0)} - E_m^{(0)}}\,.}

What the formula tells you

Two features are worth memorizing:

  1. The numerator is always non-negative because it is 2|\cdot|^2. The sign of each term is therefore set entirely by the energy denominator.
  2. The ground state always shifts down. For the lowest level, every other state has higher energy, so E0(0)Em(0)<0E_0^{(0)} - E_m^{(0)} < 0 for all mm, making every term negative. Second-order perturbation theory lowers the ground-state energy — a manifestation of level repulsion: nearby states "push apart," and the lower one is pushed further down.

The single-term case

For a two-level system there is only one other state, so the sum collapses to one term. With H^(0)=diag(E0(0),E1(0))\hat{H}^{(0)} = \mathrm{diag}(E_0^{(0)}, E_1^{(0)}) and real coupling V=1(0)H^0(0)V = \langle 1^{(0)}|\hat{H}'|0^{(0)}\rangle, the ground state's second-order shift is

E0(2)=V2E0(0)E1(0).E_0^{(2)} = \frac{|V|^2}{E_0^{(0)} - E_1^{(0)}}.

If E0(0)=1E_0^{(0)} = 1, E1(0)=4E_1^{(0)} = 4, and V=0.6V = 0.6, then E0(2)=0.36/(14)=0.12E_0^{(2)} = 0.36 / (1 - 4) = -0.12 — negative, as expected for a ground state.

Try it

Compute the second-order ground-state shift E0(2)E_0^{(2)} for E0(0)=1E_0^{(0)} = 1, E1(0)=4E_1^{(0)} = 4, and real coupling V=0.6V = 0.6, and return it as a number.

Run your code to see the quantum state.

Sign in on the full site to ask questions and join the discussion.