\begin{equation}
\label{eq:ham}
H= \sum_{i=1}^N \left( \frac{{\bf p_i^2}}{2m} - \frac{e^{2}}{4\pi\epsilon_{0}}
\frac{Z}{r} \right)
+ \sum_{i