Formula_2887

\[H_{\omega_{s}}=-\sum_{k \in \omega_{s}} \frac{1}{2 m_{k}} \Delta r_{k}+\sum_{i, j \in \omega_{s}} v_{i j}\left(r_{i}-r_{j}\right),\]