E=E'+\mathbf{V}\cdot\mathbf{P}'+\frac{1}{2}\mu V^2
