move open_potential() from pair.cpp to force.cpp and extend to bonds/angles