where ∆q i = qi -q are the errors in internal coordinates, and f q i,i are the diagonal elements of the Hq Hessian in internal coordinates (force constants). In this way, Eq. 9 enables us to partition GEO into the contributions from the errors in geometric parameters in internal coordinates (bond lengths, bond angles, torsion angles, etc.). We also define: