This is probably a very trivial question… why, in pair_tersoff.cpp, are the analytical derivatives of the different parts of the potential divided by r?
It’s so that the components of force can be just something
like fx = delx * fforce, instead of fx = delx/r * force. I.e.
the 1/r is pulled into fforce. E.g. in pair_lj_cut.cpp: