ELSEVIER Journal of Molecular Structure (Theochem) 311 (1994) 9-17 THEO CHEM Investigation of electron correlation effects on molecular geometries Brian J. Teppen", David M. Miller", Ming Cao", Regina F. Frey", Susan Q. Newton", Frank A. Momany", Michael Ramek", Lothar Schafer'v" "Department of Agronomy, University of Arkansas, Fayetteville, AR 72701, USA "Depanment of Chemistry and Biochemistry. University of Arkansas, Fayetteville, AR 72701, USA clBM Corporation, 500 Maryville College Drive. St Louis. MO. 63141. USA dMolecular Simulations Inc., Burlington, MA 01803-5297. USA "Institutefur Physikalische und Theoretische Chemie, Technische Universitdt Graz . A-80lO Graz , Austria (Received 19 July 1993; accepted 9 August 1993) Abstract A survey is presented as part of an ongoing program in which the effects of electron correlation on molecular structures are studied for the purpose of obtaining data that can be used in developing force field parameters for empirical molecular modeling. In our approach, molecular structures devoid of dispersion forces, i.e. optimized at the Hartree- Fock (HF) level, are compared with the same structures optimized at the MP2 level, in which the dispersion forces are "switched on". So far, the investigations have included hydrocarbons, to study aliphatic non-bonded interactions. systems with internal hydrogen bonding (ethylene glycol, glycerol, glycine and the model dipeptide N-formyl alanine amide), and silicates, which are studied in order to derive force field parameters applicable to modeling clays. In hydrocarbons, electron correlation effects lead to contractions in I,S-non-bonded distances in all rotamers with gauche (G) torsions, particularly in GG sequences, in which small changes in torsional angles « 5°) can cause large changes (a few tenths of an angstrom) in non-bonded distances in the attractive region of the van der Waals potential. In ethylene glycol and glycerol, hydrogen bonds in MP2 structures are O. I to 0.2 A shorter than the same bonds in the correspond- ing HF structures. In silicates, significant correlation effects are found both in primary structural parameters and torsional angles, and again long-range non-bonded distances are contracted in structures in which correlation effects are active. In the context of the widely used single-point energy calculations at the correlation level (for example single- point MP2 energies calculated for HF geometries), we find that the errors connected with the lack of geometry optimization at the MP2 level increase with the size of a system and are larger for folded conformations than stretched ones. 1. Introduction Developing force field parameters for empirical molecular modeling (MM) is a difficult task for many reasons, including the sparsity of the experi- "Corresponding author. mental database and the frequently contradictory nature of the experimental data. Therefore, non- experimental sources of information which can be used in this context are increasingly being explored. Among these, ab initio calculations have been applied in various ways [1]. One of the most important contributions to molecular stability derives from intramolecular 0166-1280/94/$07.00 © 1994 Elsevier Science B.V. All rights reserved SSDl 0166-1280(93)03505-2