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

The possibility of obtaining an accurate site-site potential model suitable for use in molecular dynamics (MD) simulations of methane from ab initio calculations has been explored. Counterpoise-corrected (CPC), supermolecule, ab initio energies at the MP2/6-311+G(2df,2pd) level were computed for eleven relative orientations of two methane molecules as a function of C-C separation distance. C-C, C-H, and H-H interaction parameters in a pairwise-additive, site-site potential model for rigid methane molecules were regressed from the ab initio energies, and the resultant model accurately reproduced the ab initio energies. The model suggests that C-H attractions are dominant in weakly binding the methane dimer. CPC energies for methane trimers, tetramers, and a pentamer were also calculated at the same level. The results indicate that the n-mer energy per pair of interactions monotonically converges with increasing n, but that the assumption of pairwise additivity commonly used in MD simulations is reasonably ...