A novel decomposition of torsional potentials into pairwise interactions: A study of energy second derivatives

A general method of analyzing intramolecular torsional potentials in terms of energy second derivatives that couple the rotating atoms is presented. The method offers a rigorous decomposition of the total torsional potential into pairwise (dihedral) interactions and enables one to derive nonbonded torsional interactions between 1–4 atoms as well as between more distant atoms and sites. The method is demonstrated on ethane, propane and acetaldehyde. It is shown that the 1–4 H…H dihedral potentials in ethane and propane are very similar, thereby supporting the notion of transferable force field potential functions. However, the dihedral potentials that are obtained differ from 1–4 potentials that are used in current force fields. Intramolecular three body effects are clearly seen in this method and are found to be relatively large for the dihedral interactions, although in the one case studied (propane) the overall effect on the methyl‐methyl interaction is negligible due to cancellation of terms. The analysis explicitly shows that the barrier in acetaldehyde is due mainly to the dihedral H…H interaction.

[1]  R. Wyatt,et al.  Origin of the Barrier Hindering Internal Rotation in Ethane , 1964 .

[2]  H. Scheraga,et al.  METHOD FOR CALCULATION INTERNAL ROTATION BARRIERS. , 1965, The Journal of chemical physics.

[3]  Leo Radom,et al.  Molecular orbital theory of the electronic structure of organic compounds. XIII. Fourier component analysis of internal rotation potential functions in saturated molecules , 1972 .

[4]  D. E. Williams,et al.  Coulombic interactions in crystalline hydrocarbons , 1974 .

[5]  Norman L. Allinger,et al.  Conformational analysis. 125. The importance of twofold barriers in saturated molecules , 1977 .

[6]  L. Bartell.  REPRESENTATIONS OF MOLECULAR FORCE FIELDS. 3. GAUCHE CONFORMATIONAL ENERGY , 1977 .

[7]  Peter Pulay,et al.  Systematic AB Initio Gradient Calculation of Molecular Geometries, Force Constants, and Dipole Moment Derivatives , 1979 .

[8]  Yasuki Endo,et al.  Microwave spectra of deuterated ethanes: Internal rotation potential function and rz structure , 1981 .

[9]  Force field, dipole moment derivatives, and vibronic constants of benzene from a combination of experimental and ab initio quantum chemical information , 1981 .

[10]  S. Lifson,et al.  Born–Oppenheimer energy surfaces of similar molecules: Interrelations between bond lengths, bond angles, and frequencies of normal vibrations in alkanes , 1982 .

[11]  U. Singh,et al.  A NEW FORCE FIELD FOR MOLECULAR MECHANICAL SIMULATION OF NUCLEIC ACIDS AND PROTEINS , 1984 .

[12]  R. Amos.  Dipole moment derivatives of H2O and H2S , 1984 .

[13]  Lennart Nilsson,et al.  Empirical energy functions for energy minimization and dynamics of nucleic acids , 1986 .

[14]  P. Kollman,et al.  An all atom force field for simulations of proteins and nucleic acids , 1986, Journal of computational chemistry.

[15]  Y. Miwa,et al.  Molecular mechanics simulations of thermodynamic functions and infrared spectra of alkanes , 1988 .

[16]  A. T. Hagler,et al.  Determination of atomic point charges and point dipoles from the Cartesian derivatives of the molecular dipole moment and second moments, and from energy second derivatives of planar dimers. I. Theory , 1989 .

[17]  A. Hagler,et al.  Determination of atomic point charges and point dipoles from the Cartesian derivatives of the molecular dipole moment and second moments, and from energy second derivatives of planar dimers. II. Applications to model systems , 1989 .

[18]  Arnold T. Hagler,et al.  Direct evaluation of nonbonding interactions from ab initio calculations , 1989 .