2nd order pertubation theory for harmonic oscillator I'm having some trouble calculating the 2nd order energy shift in a problem.
I am given the pertubation:
$\hat{H}'=\alpha \hat{p}$,
where $\alpha$ is a constant, and $\hat{p}$ is given by:
$p=i\sqrt{\frac{\hbar m\omega }{2}}\left( {{a}_{+}}-{{a}_{-}} \right)$,
where ${a}_{+}$ and ${a}_{-}$ are the usual ladder operators.
Now, according to my book, the 2nd order energy shift is given by:
$E_{n}^{2}=\sum\limits_{m\ne n}{\frac{{{\left| \left\langle  \psi _{m}^{0} \right|H'\left| \psi _{n}^{0} \right\rangle  \right|}^{2}}}{E_{n}^{0}-E_{m}^{0}}}$
Now, what I have tried to do is to calculate the term inside the power of 2. And so far I have done this:
$\begin{align}
  & E_{n}^{1}=\alpha i\sqrt{\frac{\hbar m\omega }{2}}\int{\psi _{m}^{*}\left( {{{\hat{a}}}_{+}}-{{{\hat{a}}}_{-}} \right)}\,{{\psi }_{n}}\,dx=\alpha i\sqrt{\frac{\hbar m\omega }{2}}\left( \int{\psi _{m}^{*}\,{{{\hat{a}}}_{+}}{{\psi }_{n}}\,dx-\int{\psi _{m}^{*}\,{{{\hat{a}}}_{-}}{{\psi }_{n}}\,dx}} \right) \\ 
 & =\alpha i\sqrt{\frac{\hbar m\omega }{2}}\left( \sqrt{n+1}\int{\psi _{m}^{*}\,{{\psi }_{n+1}}\,dx-\sqrt{n}\int{\psi _{m}^{*}\,{{\psi }_{n-1}}\,dx}} \right)  
\end{align}
$
As you can see, I end up with the two integrals. But I don't know what to do next ? 'Cause if $m > n$, and only by 1, then the first integral will be 1, and the other will be zero. And if $n > m$, only by 1, then the second integral will be 1, and the first will be zero. Otherwise both will be zero.
And it seems wrong to have to make two expressions for the energy shift for $n > m$ and $m > n$.
So am I on the right track, or doing it totally wrong ?
Thanks in advance.
Regards
 A: Well, you end up with integrals but those are very, very easy to solve for the harmonic oscillator!
Since your problem is already formulated in terms of the raising and lowering operators $a_+$ and $a_-$. Recall that
$$a_+ | n \rangle = \sqrt{n+1} | n+1 \rangle$$
$$a_- | n \rangle = \sqrt{n} | n - 1 \rangle$$
$$ a_+ a_- | n \rangle = n | n \rangle$$
where $|n\rangle$ is short-hand for $|\psi_n^0 \rangle$.
These relations make it almost trivial to compute matrix elements involving eigenstates of the harmonic oscillator and those operators. 
Just as an example, let's prove that all eigenstates have zero expectation value for $x$: We know $x$ is proportional to $a_+ + a_-$. Inserting that into the matrix element gives us 
$$\langle n | x | n \rangle \propto \langle n | a_+ | n \rangle
+ \langle n | a_- | n \rangle = \sqrt{n+1} \langle n | n+1\rangle + \sqrt{n} \langle n | n- 1 \rangle = 0$$
because eigenstates are orthogonal.
EDIT: Continuing with the derivation where you left off, you see that you get a non-zero contribution only if $m = n+1$ or $m = n-1$. So in the infinite sum over all states $m \not= n$, only two terms will contribute, making it possible to easily carry out that sum: Just add those two non-zero terms.
