Using operators to expand the block matrices forming the Hessian of a molecular potential

J Comput Chem. 2014 Jun 5;35(15):1149-58. doi: 10.1002/jcc.23609. Epub 2014 Apr 16.

Abstract

We derive compact expressions of the second-order derivatives of bond length, bond angle, and proper and improper torsion angle potentials, in terms of operators represented in two orthonormal bases. Hereby, simple rules to generate the Hessian of an internal coordinate or a molecular potential can be formulated. The algorithms we provide can be implemented efficiently in high-level programming languages using vectorization. Finally, the method leads to compact expressions for a second-order expansion of an internal coordinate or a molecular potential.

Keywords: Hessian; internal coordinates; molecular potentials; operators; second-order derivatives.