|q⟩ Bad Qubits

intermediate · Physics · Time-Independent Perturbation Theory

First-Order State Corrections

How the state shifts

In the previous lesson we projected the first-order equation onto n(0)\langle n^{(0)}| to extract the energy. To find the state correction n(1)|n^{(1)}\rangle we project onto a different unperturbed state m(0)\langle m^{(0)}| with mnm \neq n. Starting from

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

and using m(0)H^(0)=Em(0)m(0)\langle m^{(0)}|\hat{H}^{(0)} = E_m^{(0)}\langle m^{(0)}| together with m(0)n(0)=0\langle m^{(0)}|n^{(0)}\rangle = 0, the En(1)E_n^{(1)} term drops out and we get

(Em(0)En(0))m(0)n(1)+m(0)H^n(0)=0.\big(E_m^{(0)} - E_n^{(0)}\big)\langle m^{(0)}|n^{(1)}\rangle + \langle m^{(0)}|\hat{H}'|n^{(0)}\rangle = 0.

Solving for the overlap of n(1)|n^{(1)}\rangle with each m(0)|m^{(0)}\rangle gives the coefficients. We choose the normalization n(0)n(1)=0\langle n^{(0)}|n^{(1)}\rangle = 0 (the correction has no component along the original state — that would only rescale it). Expanding n(1)|n^{(1)}\rangle in the unperturbed basis:

n(1)=mnm(0)H^n(0)En(0)Em(0)m(0).\boxed{\,|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\,.}

Reading the formula

Each unperturbed state m(0)|m^{(0)}\rangle mixes into n|n\rangle by an amount set by two factors:

This is why a vanishing gap (degeneracy) is dangerous: the coefficient diverges, and ordinary perturbation theory fails. It also explains the convergence condition from the first lesson — the ratio of coupling to gap must be small for the correction to be a genuine small admixture.

A two-level example

Take H^(0)=diag(E0(0),E1(0))\hat{H}^{(0)} = \mathrm{diag}(E_0^{(0)}, E_1^{(0)}) with E0(0)=1E_0^{(0)} = 1, E1(0)=4E_1^{(0)} = 4, and a perturbation with off-diagonal coupling 1(0)H^0(0)=0.6\langle 1^{(0)}|\hat{H}'|0^{(0)}\rangle = 0.6. The only state that can mix into 0|0\rangle is 1|1\rangle, with coefficient

c1=1(0)H^0(0)E0(0)E1(0)=0.614=0.63=0.2,c_1 = \frac{\langle 1^{(0)}|\hat{H}'|0^{(0)}\rangle}{E_0^{(0)} - E_1^{(0)}} = \frac{0.6}{1 - 4} = \frac{0.6}{-3} = -0.2,

so 0(1)=0.21|0^{(1)}\rangle = -0.2\,|1\rangle and the corrected state is approximately 00.21|0\rangle - 0.2\,|1\rangle (before renormalization).

Try it

Using the example above (E0(0)=1E_0^{(0)} = 1, E1(0)=4E_1^{(0)} = 4, coupling 1(0)H^0(0)=0.6\langle 1^{(0)}|\hat{H}'|0^{(0)}\rangle = 0.6), compute the mixing coefficient c1c_1 of 1|1\rangle in the first-order correction to the ground state, 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.