Methodologies for creating reactive potential energy surfaces from molecular mechanics force-fields are becoming increasingly popular. To date, molecular mechanics force-fields in biochemistry and small molecule organic chemistry tend to use harmonic expressions to treat bonding stretches, which is a poor approximation in reactive and non-equilibirum molecular dynamics simulations since bonds are often displaced significantly from their equilibrium positions. For such applications there is need for a better treatment of anharmonicity. In this contribution, Morse bonding potentials have been extensively parameterized for the atom types in the MM3 force field of Allinger and co-workers using high level CCSD(T)(F12*) energies. To our knowledge this is amongst the first instances of a comprehensive parametrization of Morse potentials in a popular organic chemistry force field. In the context of molecular dynamics simulations, this data will: (1) facilitate the fitting of reactive potential energy surfaces using empirical valence bond approaches, and (2) enable more accurate treatments of energy transfer.
23.Hattig, C.; Tew, D. P.; Kohn, A., Communications: Accurate and efficient approximations to explicitly correlated coupled-cluster singles and doubles, CCSD-F12.