I know that if a quantum mechanical particle of mass $m$ is subject to the potential $V(x)=\frac{1}{2}m\omega^2x^2$ then the wave functions $\phi_n(x)=\frac{1}{\sqrt{n!}}(a^t)^n\phi_0(x)$ are eigenfunctions of the Hamiltonian, with eigenvalues $E_n=\hbar \omega(\frac{1}{2}+n)$. Here $a=\frac{1}{\sqrt{2\hbar m\omega}}(m\omega x+ip)$ and $a^t=\frac{1}{\sqrt{2\hbar m\omega}}(m\omega x-ip)$
But what if we modify the potential $V(x)=\frac{1}{2}m\omega^2x^2+kx$ for some constant $k$? How do the eigenvalues change?
I also have the relation $[a,(a^t)^n]=n(a^t)^{n-1}$ and know how to express the position and momentum operators $x$ and $p$ in terms of $a,a^t$. All I can get to is this:
$H\phi_n=E_n\phi_n+kx\phi_n=E_n\phi_n + k\sqrt{\frac{\hbar}{2m\omega}}(a+a^t)\phi_n$
but don't know where to go from here since $a\phi_n=\sqrt{n}\phi_{n-1}$ and $a^t\phi_n=\sqrt{n+1}\phi_{n+1}$ don't seem to be of any use.