\begin{equation}
E_{v}=(v+\frac{1}{2})\hbar\omega-(v+\frac{1}{2})^{2}\hbar\omega x_{e}
\end{equation}