Improved Force Field Model for the Deep Eutectic Solvent Ethaline: Reliable Physicochemical Properties

In this work, we combined various parameters found in the literature for the choline cation, chloride anion, and ethylene glycol to set up force field models (FFMs) for a eutectic mixture, namely, ethaline (1:2 choline chloride/ethylene glycol (ChCl:2EG)). The validation of these models was carried out on the basis of physical and chemical properties, such as the density, expansion coefficient, enthalpy of vaporization, self-diffusion coefficients, isothermal compressibility, surface tension, and shear viscosity. After the initial evaluation of the FFMs, a refinement was found necessary and accomplished by taking into account polarization effects in a mean-field manner. This was achieved by rescaling the electrostatic charges of the ions based on partial charges derived from ab initio molecular dynamics (MD) simulations of the bulk system. Classical all-atom MD simulations performed over a large range of temperatures (298.15-373.15 K) using the refined FFMs clearly showed improved results, allowing a better prediction of experimental properties. Specific structural properties (radial distribution functions and hydrogen bonding) were then analyzed in order to support the adequacy of the proposed refinement. The final selected FFM leads to excellent agreement between simulated and experimental data on dynamic and structural properties. Moreover, compared to the previously reported force field model (Perkins, S. L.; Painter, P.; Colina, C. M. Experimental and Computational Studies of Choline Chloride-Based Deep Eutectic Solvents. J. Chem. Eng. Data 2014, 59, 3652-3662), a 10% improvement in simulated transport properties, i.e., self-diffusion coefficients, was achieved. The isothermal compressibility, surface tension, and shear viscosity for ethaline are accessed in MD simulations for the first time.