|q⟩ Bad Qubits

intermediate · Physics · Time-Independent Perturbation Theory

Checkpoint: Perturb a System

Bring it together

This checkpoint combines everything from the module on a single concrete system. You will compute the ground-state energy of a two-level Hamiltonian to second order in perturbation theory and compare against the exact answer.

The system

The unperturbed Hamiltonian is diagonal,

H^(0)=(2005),\hat{H}^{(0)} = \begin{pmatrix} 2 & 0 \\ 0 & 5 \end{pmatrix},

with ground state 0|0\rangle at E0(0)=2E_0^{(0)} = 2 and excited state 1|1\rangle at E1(0)=5E_1^{(0)} = 5. The perturbation, in the same basis, is

H^=(0.40.30.30.2).\hat{H}' = \begin{pmatrix} 0.4 & 0.3 \\ 0.3 & -0.2 \end{pmatrix}.

The recipe, step by step

First order — the diagonal expectation value:

E0(1)=0H^0=0.4.E_0^{(1)} = \langle 0|\hat{H}'|0\rangle = 0.4.

Second order — sum over the other state (here just 1|1\rangle):

E0(2)=1H^02E0(0)E1(0)=(0.3)225=0.093=0.03.E_0^{(2)} = \frac{|\langle 1|\hat{H}'|0\rangle|^2}{E_0^{(0)} - E_1^{(0)}} = \frac{(0.3)^2}{2 - 5} = \frac{0.09}{-3} = -0.03.

Adding everything, the perturbative ground-state energy is

E0E0(0)+E0(1)+E0(2)=2+0.40.03=2.37.E_0 \approx E_0^{(0)} + E_0^{(1)} + E_0^{(2)} = 2 + 0.4 - 0.03 = 2.37.

Checking against the exact answer

This 2×22\times 2 problem can be diagonalized exactly. The full ground Hamiltonian is H^=(2.40.30.34.8)\hat{H} = \begin{pmatrix} 2.4 & 0.3 \\ 0.3 & 4.8 \end{pmatrix}, whose lower eigenvalue is

E0exact=2.4+4.82(4.82.42)2+0.32=3.61.22+0.092.3631.E_0^{\text{exact}} = \frac{2.4 + 4.8}{2} - \sqrt{\left(\frac{4.8 - 2.4}{2}\right)^2 + 0.3^2} = 3.6 - \sqrt{1.2^2 + 0.09} \approx 2.3631.

Our second-order estimate 2.372.37 sits within about 0.0070.007 of the exact 2.36312.3631 — the residual is the third-order-and-beyond contribution. The agreement confirms the perturbation is genuinely small here (coupling 0.30.3 versus gap 2.6\approx 2.6), exactly the validity condition from earlier in the module.

Try it

For the system above, compute the ground-state energy to second order, E0=E0(0)+E0(1)+E0(2)E_0 = E_0^{(0)} + E_0^{(1)} + E_0^{(2)}, 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.