I consider the following Hamiltonian $$H=\frac{p^2}{2m}+\frac{m\omega^2}{2}x^2+\Theta(t)Fx,$$ where $F$ is an external constant force. So the Hamiltonian describes an unperturbed harmonic oscillator if $t<0$ and a constant perturbed harmonic oscillator if $t>0$.
My aim is to calculate the exact transition probability from the initial unperturbed groundstate to a final new excited state.
What I already know: The ecxited eigenvalues ar $\bar{E}_n=E_n-\frac{F^2}{2m\omega^2}$
I need the transition amplitude $\langle \bar{n}|0\rangle$ to calculate the probability.
The time evolution for the new excited states are $|\bar{n}\rangle=e^{-\frac{i}{\hbar}\bar{E}_nt}|n\rangle$. The state $|n\rangle$ I can express with the ladder operators.
Now I have some problems. I calculate the amplitude, $e^{\frac{i}{\hbar}\bar{E}_nt}\langle n|0\rangle=0$, because $|n\rangle$ and $|0\rangle$ are orthogonal. Of course this schould be nonzero.
So where is the mistake?