\mathbf{F}=-\partial U/\partial \mathbf{r}_a
