E=\sum_a\frac{1}{2}m_av_a^2+U(\mathbf{r}_1,\mathbf{r}_2,\ldots)
