ABSTRACT
Deep eutectic solvents (DES) have gained a reputation as inexpensive and easy to handle ionic liquid analogues. This work employs molecular dynamics (MD) to simulate a variety of DES. The hydrogen bond acceptor (HBA) choline chloride was paired with the hydrogen bond donors (HBD) glycerol, 1,4-butanediol, and levulinic acid. Levulinic acid was also paired with the zwitterionic HBA betaine. In order to evaluate the reliability of data MD simulations can provide for DES, two force fields were compared: the Merck Molecular Force Field and the General Amber Force Field with two different sets of partial charges for the latter. The force fields were evaluated by comparing available experimental thermodynamic and transport properties against simulated values. Structural analysis was performed on the eutectic systems and compared to non-eutectic compositions. All force fields could be validated against certain experimental properties, but performance varied depending on the system and property in question. While extensive hydrogen bonding was found for all systems, details about the contribution of individual groups strongly varied among force fields. Interaction potentials revealed that HBA–HBA interactions weaken linearly with increasing HBD ratio, while HBD–HBD interactions grew disproportionally in magnitude, which might hint at the eutectic composition of a system.
Acknowledgments
The authors thank Prof. Wolfgang Arlt for supporting the experimental part of this study and Dr.-Ing. Maximilian Greiner for valuable discussions on the molecular dynamics simulations. This work has made use of the computer resources provided by the Leibniz Supercomputing Centre.
Disclosure statement
No potential conflict of interest was reported by the authors.