add numerical gradient and hessian to frozen atom potential
add thermodynamics to base system