neg_hessian = np.zeros((num_par, num_par))   # Hessian matrix