Fick and Maxwell-Stefan diffusion of the liquid mixture cyclohexane + toluene + acetone + methanol and its subsystems

Transport diffusion of the quaternary mixture cyclohexane + toluene + acetone + methanol and most of its binary and ternary subsystems at ambient conditions of temperature and pressure is studied by molecular dynamics simulation. Liquid-liquid phase separation of cyclohexane + methanol extends into the according ternary and quaternary mixtures, rendering them highly non-ideal. Maxwell-Stefan diffusion and intradiffusion coefficients are predicted with the Green-Kubo formalism, while the thermodynamic factor is sampled with Kirkwood-Buff integration. From these data, the Fick diffusion coefficient is determined and thoroughly compared with experimental literature values. Given that the simulation results for all binaries, one of the ternaries and the quaternary mixture are successfully validated, it can be assumed that the results for the three remaining ternary mixtures, for which no experimental data exist, are sound, too. These three ternaries exhibit non-idealities that entail Maxwell-Stefan diffusion coefficient maxima which are accompanied by Fick diffusion coefficient minima due to the thermodynamic factor. Particularly the very different nature of cyclohexane and methanol molecules, where the latter strongly self-associate, leads to anomalies. The predictive models by Li et al. [2001, Fluid Phase Equilibr., 187, 195] and Allie-Ebrahim et al. [2018, Phys. Chem. Chem. Phys., 20, 18436] for binary and ternary mixtures, respectively, are found to yield good results for predicting the Fick diffusion coefficient of these challenging systems.