Monte Carlo predictions of phase equilibria and structure for dimethyl ether + sulfur dioxide and dimethyl ether + carbon dioxide

A new force field for dimethyl ether (DME) based on the Lennard-Jones (LJ) 12-6 plus point charge functional form is presented in this work. This force field reproduces experimental saturated liquid and vapor densities, vapor pressures, heats of vaporization, and critical properties to within the statistical uncertainty of the combined experimental and simulation measurements for temperatures between the normal boiling and critical point. Critical parameters and normal boiling point are predicted to within 0.1% of experiment. This force field is used in grand canonical histogram reweighting Monte Carlo simulations to predict the pressure composition diagrams for the binary mixtures DME + SO(2) at 363.15 K and DME + CO(2) at 335.15 and 308.15 K. For the DME + SO(2) mixture, simulation is able to qualitatively reproduce the minimum pressure azeotropy observed experimentally for this mixture, but quantitative errors exist, suggesting that multibody effects may be important in this system. For the DME + CO(2) mixture, simulation is able to predict the pressure-composition behavior within 1% of experimental data. Simulations in the isobaric-isothermal ensemble are used to determine the microstructure of DME + SO(2) and DME + CO(2) mixtures. The DME + SO(2) shows weak pairing between DME and SO(2) molecules, while no specific pairing or aggregation is observed for mixtures of DME + CO(2).