Estimation of Transition-Metal Empirical Parameters for Molecular Mechanical Force Fields

Force-field parameters of the first row transition metals together with a few additional common elements such as those from the second (Rh, Ru) and third (Hg, Pt) rows of elements in ligated forms were determined based on the density functional theory calculations. Bonding characteristics were determined by averaging metal–ligand force constants in optimal geometries from several chosen complexes of each metal in the most common oxidation numbers and structural arrangements. Parameters of Lennard-Jones potential were determined based on a supermolecular model. Our determined molecular mechanical parameters are compared with presently available parameters published by other groups. We performed two different kinds of testing in order to demonstrate the reliability of these parameters in the case of ligated metallo complexes. First, the nonbonding potential was constructed for an additional set of 19 larger systems containing common complexes with organic molecules. The second test compares the Pt–O and Pt–H radial distribution functions for cisplatin in a box of TIP3P water with lately published studies