in c++ potential interactions change hessian function

hessian(...) -> energy_gradient_hessian(...)
add hessian to LJCut