Mutual solubilities study for binary mixtures of dipropylene glycol dimethyl ether and water via mol

Abstract It is a challenging task to accurately compute the structures and physicochemical properties of partially soluble liquid mixtures using classical force field simulation methods due to the limitations of both sampling and intermolecular potential. In the present study, we employed replica exchange molecular dynamics in combination with the AMOEBA polarizable force field to simulate partially miscible liquid mixtures of dipropylene glycol dimethyl ether and water. We demonstrated that our strategy is promising in studying such complicated binary mixture systems without involving empirical fitting of any available experimental data of mixtures. We correctly predicted the trend of liquid concentration distribution as the function of temperature, and the concentration drifts found in the simulations are consistent with the concentration difference between experimental values and initial configurations. Based on self-diffusion coefficient analysis, the lower critical solution temperature of the system was identified at about 290 K, which is in satisfactory agreement with experimental evidence. Although accurate concentration prediction of such a partially miscible system using direct molecular dynamics method remains an intractable task, we expect that such atomistic simulations will play a more important role in future studies of liquid mixtures.

[1]  H. G. Petersen,et al.  Error estimates on averages of correlated data , 1989 .

[2]  A. Mark,et al.  Convergence and sampling efficiency in replica exchange simulations of peptide folding in explicit solvent. , 2007, The Journal of chemical physics.

[3]  H. C. Andersen Rattle: A “velocity” version of the shake algorithm for molecular dynamics calculations , 1983 .

[4]  Alexander D. MacKerell,et al.  All-atom empirical potential for molecular modeling and dynamics studies of proteins. , 1998, The journal of physical chemistry. B.

[5]  K. Gubbins,et al.  Vapor-liquid and vapor-liquid-liquid equilibria of carbon dioxide/n-perfluoroalkane/n-alkane ternary mixtures. , 2005, The journal of physical chemistry. B.

[6]  J. I. Siepmann,et al.  Pressure dependence of the vapor-liquid-liquid phase behavior in ternary mixtures consisting of n-alkanes, n-perfluoroalkanes, and carbon dioxide. , 2005, The journal of physical chemistry. B.

[7]  Monica H Lamm,et al.  Molecular simulation of complete phase diagrams for binary mixtures , 2001 .

[8]  P. Jungwirth,et al.  Propensity of heavier halides for the water/vapor interface revisited using the Amoeba force field , 2005 .

[9]  M. Mazutti,et al.  Phase equilibrium data and thermodynamic modelling of the system (propane + DMF + methanol) at high pressures , 2011 .

[10]  T. Darden,et al.  Particle mesh Ewald: An N⋅log(N) method for Ewald sums in large systems , 1993 .

[11]  G. Pellicane,et al.  A comprehensive study of the phase diagram of symmetrical hard-core Yukawa mixtures , 1998 .

[12]  Alexander D. MacKerell,et al.  Atomistic Models and Force Fields , 2001 .

[13]  Pengyu Y. Ren,et al.  Calculation of protein–ligand binding free energy by using a polarizable potential , 2008, Proceedings of the National Academy of Sciences.

[14]  W. L. Jorgensen,et al.  Development and Testing of the OPLS All-Atom Force Field on Conformational Energetics and Properties of Organic Liquids , 1996 .

[15]  Carsten Kutzner,et al.  GROMACS 4:  Algorithms for Highly Efficient, Load-Balanced, and Scalable Molecular Simulation. , 2008, Journal of chemical theory and computation.

[16]  A. Roitberg,et al.  Coupling of replica exchange simulations to a non-Boltzmann structure reservoir. , 2007, The journal of physical chemistry. B.

[17]  C. Geyer Markov Chain Monte Carlo Maximum Likelihood , 1991 .

[18]  Yuko Okamoto,et al.  Replica-exchange Monte Carlo method for the isobaric isothermal ensemble , 2001 .

[19]  Claire S. Adjiman,et al.  Understanding the fluid phase behaviour of crude oil: Asphaltene precipitation , 2011 .

[20]  Pengyu Y. Ren,et al.  Ion solvation thermodynamics from simulation with a polarizable force field. , 2003, Journal of the American Chemical Society.

[21]  W. L. Jorgensen,et al.  Comparison of simple potential functions for simulating liquid water , 1983 .

[22]  A. Leach Molecular Modelling: Principles and Applications , 1996 .

[23]  Jean-Philip Piquemal,et al.  Polarizable molecular dynamics simulation of Zn(II) in water using the AMOEBA force field. , 2010, Journal of Chemical Theory and Computation.

[24]  E. Maginn,et al.  Isomolar semigrand ensemble molecular dynamics: development and application to liquid-liquid equilibria. , 2005, The Journal of chemical physics.

[25]  Peter A. Kollman,et al.  AMBER, a package of computer programs for applying molecular mechanics, normal mode analysis, molecular dynamics and free energy calculations to simulate the structural and energetic properties of molecules , 1995 .

[26]  Margaret E. Johnson,et al.  Current status of the AMOEBA polarizable force field. , 2010, The journal of physical chemistry. B.

[27]  Pengyu Y. Ren,et al.  Consistent treatment of inter‐ and intramolecular polarization in molecular mechanics calculations , 2002, J. Comput. Chem..

[28]  David M. Eike,et al.  Atomistic simulation of solid-liquid coexistence for molecular systems: application to triazole and benzene. , 2006, The Journal of chemical physics.

[29]  W. L. Jorgensen,et al.  Prediction of Properties from Simulations: Free Energies of Solvation in Hexadecane, Octanol, and Water , 2000 .

[30]  D. van der Spoel,et al.  A temperature predictor for parallel tempering simulations. , 2008, Physical chemistry chemical physics : PCCP.

[31]  I. Marrucho,et al.  Liquid-liquid equilibrium of (perfluoroalkane + alkane) binary mixtures , 2006 .

[32]  K. Shukla Total and excess properties of nonideal ternary fluid mixtures via isothermal-isobaric molecular dynamics simulation: size and energy parameter effects , 1997 .

[33]  T. Kristóf,et al.  The extrapolation of phase equilibrium curves of mixtures in the isobaric—isothermal Gibbs ensemble , 2002 .

[34]  José Elguero,et al.  Study of the methanol trimer potential energy surface , 1997 .

[35]  M. Seibert,et al.  Reproducible polypeptide folding and structure prediction using molecular dynamics simulations. , 2005, Journal of molecular biology.