|q⟩ Bad Qubits

intermediate · Physics · Time-Dependent Perturbation & Fermi's Golden Rule

Transition Probabilities

The amplitude cf(1)(t)c_f^{(1)}(t) is complex and not directly observable. By the Born rule, the physically measurable quantity is its squared magnitude — the probability of finding the system in state f|f\rangle.

From amplitude to probability

The first-order transition probability is

Pif(t)=cf(1)(t)2.P_{i\to f}(t) = \big|c_f^{(1)}(t)\big|^2 .

For the constant perturbation H^=V^\hat H' = \hat V switched on at t=0t=0, we found

cf(1)(t)=Vfi2sin(ωfit/2)ωfi,\big|c_f^{(1)}(t)\big| = \frac{|V_{fi}|}{\hbar}\,\frac{2\,|\sin(\omega_{fi}t/2)|}{|\omega_{fi}|} ,

so squaring gives

  Pif(t)=Vfi224sin2(ωfit/2)ωfi2  \boxed{\;P_{i\to f}(t) = \frac{|V_{fi}|^2}{\hbar^2}\,\frac{4\sin^2(\omega_{fi}t/2)}{\omega_{fi}^2}\;}

where Vfi=fV^iV_{fi}=\langle f|\hat V|i\rangle and ωfi=(EfEi)/\omega_{fi}=(E_f-E_i)/\hbar.

Reading the formula

Two factors control the probability:

The small-energy-gap limit

Near ωfi=0\omega_{fi}=0 the factor has a finite limit. Using sinxx\sin x \approx x,

4sin2(ωfit/2)ωfi2    ωfi0    t2,\frac{4\sin^2(\omega_{fi}t/2)}{\omega_{fi}^2} \;\xrightarrow{\;\omega_{fi}\to 0\;}\; t^2 ,

so for transitions to (near-)degenerate states the probability grows quadratically in time:

Pif(t)    Vfi22t2.P_{i\to f}(t) \;\approx\; \frac{|V_{fi}|^2}{\hbar^2}\,t^2 .

A useful shape: the sinc function

Writing Ω=ωfi\Omega = \omega_{fi}, the energy-mismatch factor is t2sinc2(Ωt/2)t^2\,\mathrm{sinc}^2(\Omega t/2) with sinc(x)=sinx/x\mathrm{sinc}(x)=\sin x / x. As tt grows the central peak at Ω=0\Omega=0 becomes taller (height t2t^2) and narrower (width 2π/t\sim 2\pi/t). This narrowing is the seed of energy conservation: only final states with EfEiE_f \approx E_i accumulate appreciable probability, and the longer the perturbation acts, the more strictly that condition is enforced — a time–energy uncertainty relation ΔEt\Delta E\,t \sim \hbar.

Try it

For a constant perturbation with Vfi=0.3V_{fi}=0.3, ωfi=1.5\omega_{fi}=1.5, t=2.0t=2.0, and =1\hbar=1, compute Pif(t)=(Vfi2/2)4sin2(ωfit/2)/ωfi2P_{i\to f}(t) = (|V_{fi}|^2/\hbar^2)\,4\sin^2(\omega_{fi}t/2)/\omega_{fi}^2 and return it. (Check that your answer lies between 0 and 1.)

Run your code to see the quantum state.

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