Interatomic Lennard-Jones potentials of linear and branched alkanes calibrated by Gibbs ensemble simulations for vapor-liquid equilibria

We propose Lennard-Jones potential parameters for interatomic interactions of linear and branched alkanes based on matching the results of Gibbs ensemble simulations of vapor-liquid equilibria to experimental data. The alkane model is similar as in the OPLS-AA [W. L. Jorgensen, D. S. Maxwell, and J. Tirado-Rives, J. Am. Chem. Soc. 118, 11225 (1996)], but multiple atom types for carbon based on the number of covalently bonded hydrogen atoms are necessary to accurately reproduce liquid densities and enthalpies of vaporization with the errors of 2.1% and 3.3%, respectively, for hydrocarbons of various chain lengths and structures. We find that the attraction energies of the carbon atoms are almost proportional to the number of covalent hydrogen atoms with each increasing the carbon energy parameter by approximate to0.033 kcal/mol. Though the present force field outperforms the OPLS-AA force field for alkanes we studied, systematic deviations for vapor pressures are still observed with errors of 15%-30%, and critical temperatures are slightly underestimated. We think that these shortcomings are probably due to the inadequacy of the two-parameter Lennard-Jones potential, and especially its behavior at short distances. (C) 2004 American Institute of Physics.