A modification to the COSMIC parameterisation using ab initio constrained potential functions

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.

Conformational analysis. Part 16. Conformational free energies in substituted piperidines and piperidinium salts

The conformational free energies (-delta G degree) of a number of 4-substituted piperdines and piperidinium salts have been determined by the J-value method. For the 4-substituted piperdines (R = Me, Phenyl, CO2Et, Br, OH, F) the relative conformer energies are almost identical to those of the analogous cyclohexanes. The methyl and phenyl compounds showed no change in the couplings on protonation, implying no change in the conformer energies. In contrast, in the remaining compounds with polar 4-substituents an almost constant stabilisation of the axial conformer of ca. 0.7 -0.8 kcal mol-1 was observed on protonation. In three cases (R = F, OH and Br) the conformational preference is reversed on protonation and the axial form is favored. The conformer energies of both the free bases and the piperidinium salts can be quantitatively predicted by molecular mechanics calculations using the COSMIC force-field, in which the electrostatic interactions are calculated by a simple Coulombic model with the partial atomic charges in the molecules given by the CHARGE2 routine, and an effective dielectric constant of five. The precise agreement obtained demonstrates conclusively that the electrostatic interactions between the substituents and the protonated nitrogen are the cause of the conformational changes on protonation, and that these can be modelled successfully using existing force-fields.

Molecular modelling of poly(aryl ether ketones). I. Aryl..aryl interactions in crystal structures

Non-bonded potentials for the aryl..aryl interaction have been derived using crystal structure data of a number of small aromatic molecules. The potentials, based on atom-centered interactions, give an accurate reproduction of the benzene crystal geometry and sublimation energy when used in conjunction with coulombic energies evaluated using point atomic charges. An examination of the charge distribution on benzene suggested values of 0.13e (H) and -0.13e (C) to be suitable. The transferability of the potentials has been shown by prediction of crystal geometries and sublimation energies of other hydrocarbon molecules and, with additional interactions for the oxygen atom included, preliminary polymer crystal structure calculations have been carried out. These demonstrate the validity of the derived parameters by successfully predicting crystallographic unit cell dimensions and ring conformations in the poly(phenylene oxide) and poly(aryl ether ketone) crystals.