Fick Diffusion Coefficients in Ternary Liquid Systems from Equilibrium Molecular Dynamics Simulations

An approach for computing Fick diffusivities directly from equilibrium molecular dynamics (MD) simulations is presented and demonstrated for a ternary chloroform–acetone–methanol liquid mixture. In our approach, Fick diffusivities are calculated from the Maxwell–Stefan (MS) diffusivities and the so-called matrix of thermodynamic factors. MS diffusivities describe the friction between different molecular species and can be directly computed from MD simulations. The thermodynamic factor describes the deviation from ideal mixing behavior and is difficult to extract from both experiments and simulations. Here, we show that the thermodynamic factor in ternary systems can be obtained from density fluctuations in small subsystems embedded in a larger simulation box. Since the computation uses the Kirkwood–Buff coefficients, the present approach provides a general route toward the thermodynamics of the mixture. In experiments, Fick diffusion coefficients are measured, while previously equilibrium molecular dynami...

[1]  A. Karr,et al.  Analysis of Binary and Ternary Mixtures. System Acetone-Chloroform-Methyl Isobutyl Kentone , 1951 .

[2]  J. Kirkwood,et al.  The Statistical Mechanical Theory of Solutions. I , 1951 .

[3]  J. Timmermans.  The physico-chemical constants of binary systems in concentrated solutions , 1959 .

[4]  E. M. Kartzmark,et al.  Thermodynamic and other properties of methanol + acetone, carbon disulphide + acetone, carbon disulphide + methanol, and carbon disulphide + methanol + acetone , 1973 .

[5]  R. Reid,et al.  The Properties of Gases and Liquids , 1977 .

[6]  J. Gmehling.  Vapor-Liquid Equilibrium Data Collection , 1977 .

[7]  J. Gaube.  J. Gmehling, U. Onken: Vapor‐Liquid Equilibrium Data Collection, Aqueous‐Organic Systems, in der Reihe: Chemistry Data Series, Vol. I, Part. 1. DECHEMA, Frankfurt 1977. 698 Seiten, Preis: DM 120,‐ , 1978 .

[8]  Katsuji Noda,et al.  Viscosities and densities at 298.15 K for mixtures of methanol, acetone, and water , 1982 .

[9]  Richard L. Rowley,et al.  Ternary liquid mixture viscosities and densities , 1984 .

[10]  M. Schoen,et al.  The shear viscosity of a Lennard-Jones fluid calculated by equilibrium molecular dynamics , 1985 .

[11]  William L. Jorgensen,et al.  Optimized intermolecular potential functions for liquid alcohols , 1986 .

[12]  H. Kooijman,et al.  COMPOSITION DERIVATIVES OF ACTIVITY COEFFICIENT MODELS (FOR THE ESTIMATION OF THERMODYNAMIC FACTORS IN DIFFUSION) , 1991 .

[13]  Alexis T. Bell,et al.  Transport diffusivity of methane in silicalite from equilibrium and nonequilibrium simulations , 1993 .

[14]  Ross Taylor,et al.  Multicomponent mass transfer , 1993 .

[15]  G. Kuiken.  Thermodynamics of Irreversible Processes: Applications to Diffusion and Rheology , 1994 .

[16]  R. Krishna,et al.  The Maxwell-Stefan approach to mass transfer , 1997 .

[17]  D. C. Rapaport,et al.  The Art of Molecular Dynamics Simulation , 1997 .

[18]  P. Oracz,et al.  Vapour-liquid equilibria. XII. The ternary system methanolchloroformacetone at 303.15 K , 1997 .

[19]  V. Zanden,et al.  Using molecular dynamics to obtain Maxwell-Stefan diffusion coefficients in liquid systems , 1998 .

[20]  E. Ruckenstein,et al.  Entrainer effect in supercritical mixtures , 2001 .

[21]  D. Frenkel,et al.  Understanding molecular simulation : from algorithms to applications. 2nd ed. , 2002 .

[22]  Stanley I. Sandler,et al.  A Priori Phase Equilibrium Prediction from a Segment Contribution Solvation Model , 2002 .

[23]  A. Shapiro,et al.  Verifying reciprocal relations for experimental diffusion coefficients in multicomponent mixtures , 2003 .

[24]  H. Hasse,et al.  Self Diffusion and Binary Maxwell–Stefan Diffusion in Simple Fluids with the Green–Kubo Method , 2009, 0904.4629.

[25]  G. Grest,et al.  Interdiffusion of solvent into glassy polymer films: a molecular dynamics study. , 2004, The Journal of chemical physics.

[26]  D. Keffer,et al.  The Composition Dependence of Self and Transport Diffusivities from Molecular Dynamics Simulations , 2004 .

[27]  A. Shapiro,et al.  Modeling diffusion coefficients in binary mixtures , 2004 .

[28]  B. Edwards,et al.  Determination of Statistically Reliable Transport Diffusivities from Molecular Dynamics Simulation , 2004 .

[29]  D. Wheeler,et al.  Molecular Dynamics Simulations of Multicomponent Diffusion. 2. Nonequilibrium Method , 2004 .

[30]  Molecular dynamics simulation of solvent-polymer interdiffusion: Fickian diffusion. , 2003, The Journal of chemical physics.

[31]  B. Edwards,et al.  On the relationship between Fickian diffusivities at the continuum and molecular levels. , 2005, The journal of physical chemistry. B.

[32]  H. Hasse,et al.  Self-Diffusion and Binary Maxwell–Stefan Diffusion Coefficients of Quadrupolar Real Fluids from Molecular Simulation , 2005, 0904.4633.

[33]  Rajamani Krishna,et al.  The Darken Relation for Multicomponent Diffusion in Liquid Mixtures of Linear Alkanes: An Investigation Using Molecular Dynamics (MD) Simulations , 2005 .

[34]  André Bardow,et al.  Ternary diffusivities by model-based analysis of Raman spectroscopy measurements , 2006 .

[35]  Arieh Ben-Naim,et al.  Molecular Theory of Solutions , 2006 .

[36]  Wei Shi,et al.  Continuous Fractional Component Monte Carlo:  An Adaptive Biasing Method for Open System Atomistic Simulations. , 2007, Journal of chemical theory and computation.

[37]  Alberto Striolo,et al.  Hydrogen-bond dynamics for water confined in carbon tetrachloride-acetone mixtures. , 2008, The journal of physical chemistry. B.

[38]  D. Bedeaux,et al.  Non-equilibrium Thermodynamics of Heterogeneous Systems , 2008, Series on Advances in Statistical Mechanics.

[39]  C. Nieto‐Draghi,et al.  Diffusion coefficients in CO2/n-alkane binary liquid mixtures by molecular simulation. , 2008, The journal of physical chemistry. B.

[40]  Wei Shi,et al.  Improvement in molecule exchange efficiency in Gibbs ensemble Monte Carlo: Development and implementation of the continuous fractional component move , 2008, J. Comput. Chem..

[41]  W. Marquardt,et al.  Prediction of multicomponent mutual diffusion in liquids: Model discrimination using NMR data , 2009 .

[42]  Randall Q. Snurr,et al.  A new perspective on the order-n algorithm for computing correlation functions , 2009 .

[43]  Stan G Moore,et al.  Improved implementation of Kirkwood-Buff solution theory in periodic molecular simulations. , 2009, Physical review. E, Statistical, nonlinear, and soft matter physics.

[44]  Stanley I. Sandler,et al.  Improvements of COSMO-SAC for vapor–liquid and liquid–liquid equilibrium predictions , 2010 .

[45]  T. Vlugt,et al.  Multicomponent Maxwell−Stefan Diffusivities at Infinite Dilution , 2011 .

[46]  André Bardow,et al.  Calculating thermodynamic properties from fluctuations at small scales. , 2011, The journal of physical chemistry. B.

[47]  T. Vlugt,et al.  Maxwell-Stefan diffusivities in liquid mixtures: using molecular dynamics for testing model predictions , 2011 .

[48]  André Bardow,et al.  Predictive Darken Equation for Maxwell-Stefan Diffusivities in Multicomponent Mixtures , 2011 .

[49]  T. Vlugt,et al.  Maxwell-Stefan diffusivities in binary mixtures of ionic liquids with dimethyl sulfoxide (DMSO) and H2O. , 2011, The journal of physical chemistry. B.

[50]  A. Chandra,et al.  Structural, single-particle and pair dynamical properties of acetone–chloroform mixtures with dissolved solutes , 2011 .

[51]  G. Peters,et al.  Pair correlation function integrals: computation and use. , 2011, The Journal of chemical physics.

[52]  Larisa Zoranic,et al.  A comparative Molecular Dynamics study of water–methanol and acetone–methanol mixtures , 2011 .

[53]  Jean-Marc Simon,et al.  Thermodynamics of a small system in a μT reservoir , 2011 .

[54]  Hans Hasse,et al.  Prediction of self-diffusion coefficient and shear viscosity of water and its binary mixtures with methanol and ethanol by molecular simulation. , 2011, The Journal of chemical physics.

[55]  D. Bedeaux,et al.  Fick diffusion coefficients of liquid mixtures directly obtained from equilibrium molecular dynamics. , 2011, The journal of physical chemistry. B.

[56]  Jean-Marc Simon,et al.  Thermodynamics of small systems embedded in a reservoir: a detailed analysis of finite size effects , 2012 .

[57]  D. Bedeaux,et al.  Correction to “Fick Diffusion Coefficients of Liquid Mixtures Directly Obtained from Equilibrium Molecular Dynamics” , 2012 .

[58]  H. Hasse,et al.  On the prediction of transport properties of monomethylamine, dimethylamine, dimethylether and hydrogen chloride by molecular simulation , 2012 .

[59]  H. Hasse,et al.  Prediction of Transport Properties of Liquid Ammonia and Its Binary Mixture with Methanol by Molecular Simulation , 2012 .

[60]  Kurt Kremer,et al.  Kirkwood-Buff Analysis of Liquid Mixtures in an Open Boundary Simulation. , 2012, Journal of chemical theory and computation.