This is a suite of MATLAB codes to calculate effective Eckart inertias for large amplitude torsions published by Bryan M. Wong, Ryan L. Thom, and Robert W. Field in The Journal of Physical Chemistry A 110, 7406 - 7413 (2006). The Eckart inertias are obtained by solving a system of transcendental equations using the Powell dogleg method. Since this system is highly nonlinear, we have implemented analytical Jacobians in the dogleg method to maximize computational efficiency. The zip file includes some sample input files.
Create scripts with code, output, and formatted text in a single executable document.