Gibbs ensemble Monte Carlo simulations of binary mixtures of methane, difluoromethane, and carbon dioxide

Gibbs ensemble Monte Carlo simulations were used to study the vapor-liquid equilibrium of binary mixtures of carbon dioxide + methane and carbon dioxide + difluoromethane. The potential forcefields we employ are all atomistic models, and have not previously been mixed together to study the vapor-liquid equilibrium of the binary mixtures. In addition, we characterize the microscopic structure of these liquid mixtures. In carbon dioxide + methane at 230 K and 56 bar, the microscopic structure of carbon dioxide in the mixture is the same as that in the pure liquid. In carbon dioxide + difluoromethane at 283 K and 56 bar, the presence of carbon dioxide does not noticeably perturb the liquid structure of difluoromethane, but the structure of carbon dioxide is subtly changed, due to a strong interaction between it and difluoromethane. The simulations in the isobaric-isothermal (NPT) ensemble agree well with the experimental data, except at the two extreme regions of the pressure range. The good agreement of most simulated state points with experimental data encourages one to develop more accurate potentials for predicting the thermodynamic properties of these systems as well as other complicated systems, which are less amenable to measurement by experiment.