[molpro-user] precision problems in calculationg Hessian matrix

[张世扬]

When we calculate harmonic vibration frequencies and normal modes, the default step size of the numerical differentiation of the energy is rstep=0.001 [bohr]. And the default convergence threshold for energy is 1.d-6. In this case, we cannot get precision of Hessian matrix as the output file shows to us. Why the elements of Hessian matrix in the molpro output file can achieve this precision?