Molecular simulation of the vapor-liquid equilibria of xylene mixtures: Force field performance, and Wolf vs. Ewald for electrostatic interactions

This article explores how well vapor-liquid equilibria of pure components and binary mixtures of xylenes can be predicted using different force fields in molecular simulations. The accuracy of the Wolf method and the Ewald summation is evaluated. Monte Carlo simulations in the Gibbs ensemble are performed at conditions comparable to experimental data, using four different force fields. Similar results using the Wolf and the Ewald methods can be obtained for the prediction of densities and the phase compositions of binary mixtures. With the Wolf method, up to 50% less CPU time is used compared to the Ewald method, at the cost of accuracy and additional parameter calibration. The densities of p-xylene and m-xylene can be well estimated using the TraPPE-UA and AUA force fields. The largest differences of VLE with experiments are observed for o-xylene. The p-xylene/o-xylene binary mixtures at 6.66 and 81.3 kPa are simulated, leading to an excellent agreement in the predictions of the composition of the liquid phase compared to experiments. The composition of the vapor phase is dominated by the properties of the component with the largest mole fraction in the liquid phase. The accuracy of the predictions of the phase composition are related to the quality of the density predictions of the pure component systems. The phase composition of the binary system of xylenes is very sensitive to slight differences in vapor phase density of each xylene isomer, and how well the differences are captured by the force fields.