COSMIC(90): an improved molecular mechanics treatment of hydrocarbons and conjugated systems

Four modifications to the COSMIC molecular mechanics force field are described, which greatly increase both its versatility and the accuracy of calculated conformational energies. The Hill non-bonded van der Waals potential function has been replaced by a two-parameter Morse curve and a new H-H potential, similar to that in MM3, incorporated. Hydrocarbon energies in particular are much improved. A simple iterative Hückel pi-electron molecular orbital calculation allows modelling of conjugated systems. Calculated bond lengths and rotational barriers for a series of conjugated hydrocarbons and nitrogen heterocycles are shown to be as accurate as those determined by the MM2 SCF method. Explicit hydrogen-bonding potentials for H-bond acceptor-donor atom pairs have been included to give better hydrogen bond energies and lengths. The van der Waals radii of protonic hydrogens are reduced to 0.5 A and the energy well depth is increased to 1.0 kcal mol-1. Two new general atom types, N+sp2 and O-sp3, have been introduced which allow a wide variety of charged conjugated systems to be studied. A minimum of parameterisation is required, as the new types are easily included in the Hückel scheme which automatically adjusts bond and torsional parameters according to the defined bond-order relationships.