Determination of an ethane intermolecular potential model for use in molecular simulations from ab initio calculations

Counterpoise-corrected, supermolecule, ab initio energies obtained at the M P2/6-311+G(2df,2pd) level were computed for 22 different relative orientati ons of two ethane molecules as a function of the separation distance betwee n the molecular centers. These energies were used to regress the parameters in several simple, analytical, interatomic or site-site models that can be used for implementation in molecular simulations. Sensitivity analysis ind icates that the intermolecular potential surface is insensitive to C-C inte ractions and that the parameters in the C-C model are coupled and unobtaina ble from the dimer energies. Representation of the potential surface can be made in terms of C-H and H-H interatomic potentials if the C-C interaction s are treated as shielded. Simple Lennard-Jones and exp-6 models do not ade quately represent the potential surface using these shielded models, nor do they produce the anticipated physics for the interatomic potentials. The e xp-6 model with a damping function and the modified-Morse interatomic poten tials both reproduce the intermolecular potential surface well with physica lly realistic intersite potentials suitable for use in molecular dynamics s imulations. (C) 2001 American Institute of Physics.