Given the two orthogonal states for $H_0$ , $|n(t)>_I, |m(t)>_I$, in the interaction picture, we want to find the probability of transforming from one to the other after time t, aka:
$ \ (1) \ |<n(t)|m(t)>|^2$.
Naively, I would think to do this:
$ \ (2) \ |<n(t)|m(t)>|^2 = |<n(0)|U_I^\dagger U_I|m(0)>|^2 = |<n(0)|m(0)>|^2 = \delta_{nm}$
Where
$U_I = exp[\frac{-i}{\hbar} \int_{t_0}^{t} d\tau e^{\frac{i H_0 \tau}{\hbar}}V(\tau)e^{\frac{-i H_0 \tau}{\hbar}}]$
and V is the perturbation (and is hermitian). (U is unitary)
Using the perturbation theory I would get $\delta_{nm}$ for the 0th and 1st order, but for the second order, I get something completely different. Which way is the correct way to go, and why?
Thanks in advance!