implemented hessian in pair-wise potential, hessian function needs to be written for all potentials