Liquid–vapour equilibrium of n-alkanes using interface simulations

Abstract

Molecular Dynamics simulations were performed to calculate liquid–vapour coexisting properties of n-alkane chains up to 16 carbon atoms using interface simulations. The lattice sum or Ewald method on the dispersion forces of the Lennard–Jones potential was applied to calculate the full interaction. The liquid and vapour coexisting densities were obtained for two flexible force field models, NERD and TraPPE-UA, where the intermolecular interactions are of the Lennard–Jones type. We have recently shown [P. Orea, J. López-Lemus, and J. Alejandre, J. Chem. Phys. 123, 114702 (2005)] that the liquid–vapour densities for simple fluids do not depend on interfacial area and therefore it is possible to use a small number of molecules in a simulation. We show that the same trend is found on the simulation of these hydrocarbon molecules. The phase diagram of ethane/n-decane binary mixtures is also obtained at 410.95 K for the NERD model. The simulation results from this work were compared with those obtained using methods with interfaces using large cut-off distances and with methods without interfaces for the same potential model. In both comparisons, excellent agreement was found. The results of liquid density from the TraPPE-UA model are in good agreement with experimental data while those from the NERD model are underestimated at low temperatures. Our findings are consistent with results published by other authors for small hydrocarbons.

Acknowledgements

We wish to acknowledge the Consejo Nacional de Ciencia y Tecnología (CONACyT), México for financial support. JLL also thanks to PROMEP-México, for Grant FE03/2005. JA also thanks Mark E. Tuckerman for helpful discussions about the MTS and Nosé–Hoover chain methods. We thank the referee for helpful suggestions to improve this work.