Comparison between ab initio and polarizable molecular dynamics simulations of 1-butyl-3-methylimidazolium tetrafluoroborate and chloride in water

Lessons Learned on Obtaining Reliable Dynamic Properties for Ionic Liquids

Ionic liquids are nowadays investigated with respect to their use as electrolytes for high-performance energy storage materials. In this study, we provide a tutorial on how to calculate dynamic properties such as self-diffusion coefficients, ionic conductivities, transference numbers, as well as ion pair and ion cage dynamics, that all play a role in judging the applicability of ionic liquids as electrolytes. For the case of the ionic liquid[ C 2 C 1 Im ] [ NTf 2 ] ${[{\rm{C}}_2 {\rm{C}}_1 {\rm{Im}}][{\rm{NTf}}_2 ]}$ , we investigate the performance of different force fields. Amongst them are non-polarizable models employing unity charges, a charge-scaled version of a non-polarizable model, a polarizable model and another non-polarizable model with refined Lennard-Jones parameters. We also study the influence of the system size on the dynamic properties. While all studied force field models capture qualitatively correct trends, only the polarizable force field and the non-polarizable force field with refined Lennard-Jones parameters provide quantitative agreement to reference data, making the latter model very attractive for the reason of lower computational costs.

Ion-Specific Effects on Ion and Polyelectrolyte Solvation

Ion-specific effects on aqueous solvation of monovalent counter ions, Na+ ${^+ }$ , K+ ${^+ }$ , Cl- ${^- }$ , and Br- ${^- }$ , and two model polyelectrolytes (PEs), poly(styrene sulfonate) (PSS) and poly(diallyldimethylammonium) (PDADMA) were here studied with ab initio molecular dynamics (AIMD) and classical molecular dynamics (MD) simulations based on the OPLS-aa force-field which is an empirical fixed point-charge force-field. Ion-specific binding to the PE charge groups was also characterized. Both computational methods predict similar response for the solvation of the PEs but differ notably in description of ion solvation. Notably, AIMD captures the experimentally observed differences in Cl- ${^- }$ and Br- ${^- }$ anion solvation and binding with the PEs, while the classical MD simulations fail to differentiate the ion species response. Furthermore, the findings show that combining AIMD with the computationally less costly classical MD simulations allows benefiting from both the increased accuracy and statistics reach.

Efficient Molecular Dynamics Simulations of Deep Eutectic Solvents with First-Principles Accuracy Using Machine Learning Interatomic Potentials

In recent years, deep eutectic solvents emerged as highly tunable and ecofriendly alternatives to common organic solvents and liquid electrolytes. In the present work, the ability of machine learning (ML) interatomic potentials for molecular dynamics (MD) simulations of these liquids is explored, showcasing a trained neural network potential for a 1:2 ratio mixture of choline chloride and urea (reline). Using the ML potentials trained on density functional theory data, MD simulations for large systems of thousands of atoms and nanosecond-long time scales are feasible at a fraction of the computational cost of the target first-principles simulations. The obtained structural and dynamical properties of reline from MD simulations using our machine learning models are in good agreement with the first-principles MD simulations and experimental results. Running a single MD simulation is highlighted as a general shortcoming of typical first-principles studies if the dynamic properties are investigated. Furthermore, velocity cross-correlation functions are employed to study the collective dynamics of the molecular components in reline.

Recent Developments in Polarizable Molecular Dynamics Simulations of Electrolyte Solutions

Polarizable molecular dynamics simulations are a fast progressing field in the scientific research of ionic liquids. The fundamentals of polarizable simulations, as well as their application to ionic liquids, were summarized in a review [Bedrov, D.; Piquemal, J.-P.; Borodin, O.; MacKerell, Jr., A. D.; Roux, B.; Schröder, C. Molecular Dynamics Simulations of Ionic Liquids and Electrolytes Using Polarizable Force Fields. Chem. Rev. 2019, 119, 7940–7995] in 2019. Since then, new methods to treat intermolecular interaction of induced dipoles in these highly charged systems were developed. This concerns the damping of these interactions and additional charge transfer as well as the prediction of ionic materials with ultrahigh refractive indices. In addition to the progress of the polarizable force fields, also thermostats and barostats for polarizable simulations evolved recently.

Comparison between Ab Initio Molecular Dynamics and OPLS-Based Force Fields for Ionic Liquid Solvent Organization

OPLS-based force fields (FFs) have been shown to provide accurate bulk-phase properties for a wide variety of imidazolium-based ionic liquids (ILs). However, the ability of OPLS to reproduce an IL solvent structure is not as well-validated given the relative lack of high-level theoretical or experimental data available for comparison. In this study, ab initio molecular dynamics (AIMD) simulations were performed for three widely used ILs: the 1-butyl-3-methylimidazolium cation with chloride, tetrafluoroborate, or hexafluorophosphate anions, that is, [BMIM][Cl], [BMIM][BF₄], and [BMIM][PF₆], respectively, as a basis for further assessment of two unique IL FFs: the ±0.8 charge-scaled OPLS-2009IL FF and the OPLS-VSIL FF. The OPLS-2009IL FF employs a traditional all-atom functional form, whereas the OPLS-VSIL FF was developed using a virtual site that offloads negative charge to inside the plane of the ring with careful attention given to reproducing hydrogen bonding. Detailed comparisons between AIMD and the OPLS FFs were made based on radial distribution functions (RDFs), combined distribution functions (CDFs), and spatial distribution functions (SDFs) to examine cation-anion interactions and π⁺-π⁺ stacking between the imidazolium rings. While both FFs were able to correctly capture the general solvent structure of these popular ILs, the OPLS-VSIL FF quantitatively reproduced interaction distances more accurately. In addition, this work provides further insights into the different short- and long-range structure patterns of these popular ILs.

Ion Correlation in Choline Chloride-Urea Deep Eutectic Solvent (Reline) from Polarizable Molecular Dynamics Simulations

In recent years, deep eutectic solvents (DESs) emerged as highly tunable and environmentally friendly alternatives to common ionic liquids and organic solvents. In this work, a polarizable model based on the CHARMM Drude polarizable force field is developed for a 1:2 ratio mixture of choline chloride/urea (reline) DES. To successfully reproduce the structure of the liquid as compared to first-principles molecular dynamics simulations, a damping factor was introduced to correct the observed over-binding between the chloride and the hydrogen bonding site of choline. Investigated radial distributions reveal the formation of hydrogen bonds between all the constituents of reline and similar interactions for chloride and urea's oxygen atoms, which could contribute to the melting point depression of the mixture. Predicted dynamic properties from our polarizable force field are in good agreement with experiments, showing significant improvements over nonpolarizable models. Similar to some ionic liquids, an oscillatory behavior in the velocity autocorrelation function of the anion is visible, which can be interpreted as a rattling motion of the lighter anion surrounded by the heavier cations. The obtained results for ionic conductivity of reline show some degree of correlated ion motion in this DES. However, a joint diffusion of ion pairs cannot be observed during the simulations.

Charge transfer and polarisability in ionic liquids: a case study

The practical use of ionic liquids (ILs) is benefiting from a growing understanding of the underpinning structural and dynamic properties, facilitated through classical molecular dynamics (MD) simulations. The predictive and explanatory power of a classical MD simulation is inextricably linked to the underlying force field. A key aspect of the forcefield for ILs is the ability to recover charge based interactions. Our focus in this paper is on the description and recovery of charge transfer and polarisability effects, demonstrated through MD simulations of the widely used 1-butyl-3-methylimidazolium bis(trifluoromethanesulfonyl)imide [C₄C₁im][NTf₂] IL. We study the charge distributions generated by a range of ab initio methods, and present an interpolation method for determining atom-wise scaled partial charges. Two novel methods for determining the mean field (total) charge transfer from anion to cation are presented. The impact of using different charge models and different partial charge scaling (unscaled, uniformly scaled, atom-wise scaled) are compared to fully polarisable simulations. We study a range of Drude particle explicitly polarisable potentials and shed light on the performance of current approaches to counter known problems. It is demonstrated that small changes in the charge description and MD methodology can have a significant impact; biasing some properties, while leaving others unaffected within the structural and dynamic domains.

Chemistry Dissolved in Ionic Liquids. A Theoretical Perspective

The theoretical treatment of ionic liquids must focus now on more realistic models while at the same time keeping an accurate methodology when following recent ionic liquids research trends or allowing predictability to come to the foreground. In this Perspective, we summarize in three cases of advanced ionic liquid research what methodological progress has been made and point out difficulties that need to be overcome. As particular examples to discuss we choose reactions, chirality, and radicals in ionic liquids. All these topics have in common that an explicit or accurate treatment of the electronic structure and/or intermolecular interactions is required (accurate methodology), while at the same time system size and complexity as well as simulation time (realistic model) play an important role and must be covered as well.

Computer Simulations of Deep Eutectic Solvents: Challenges, Solutions, and Perspectives

Deep eutectic solvents (DESs) are one of the most rapidly evolving types of solvents, appearing in a broad range of applications, such as nanotechnology, electrochemistry, biomass transformation, pharmaceuticals, membrane technology, biocomposite development, modern 3D-printing, and many others. The range of their applicability continues to expand, which demands the development of new DESs with improved properties. To do so requires an understanding of the fundamental relationship between the structure and properties of DESs. Computer simulation and machine learning techniques provide a fruitful approach as they can predict and reveal physical mechanisms and readily be linked to experiments. This review is devoted to the computational research of DESs and describes technical features of DES simulations and the corresponding perspectives on various DES applications. The aim is to demonstrate the current frontiers of computational research of DESs and discuss future perspectives.

Liquid Ethylene Glycol: Prediction of Physical Properties, Conformer Population and Interfacial Enrichment with a Refined Non-Polarizable Force Field

Periodic density functional theory based molecular dynamics simulations confirm the fraction of molecules in neat liquid ethylene glycol with their central OCCO dihedral in the trans conformation to be 21%,...

Fine-Tuning the Polarizable CL&Pol Force Field for the Deep Eutectic Solvent Ethaline

Polarizable force fields are gradually becoming a common choice for ionic soft matter, in particular, for molecular dynamics (MD) simulations of ionic liquids (ILs) and deep eutectic solvents (DESs). The CL&Pol force field introduced in 2019 is the first general, transferable, and polarizable force field for MD simulations of different types of DESs. The original formulation contains, however, some problems that appear in simulations of ethaline and may also have a broader impact. First, the originally proposed atomic diameter parameters are unbalanced, resulting in too weak interactions between the chlorides and the hydroxyl groups of the ethylene glycol molecules. This, in turn, causes an artificial phase separation in long simulations. Second, there is an overpolarization of chlorides due to strong induced dipoles that give rise to the presence of peaks and antipeaks at very low q-vector values (2.4 nm^-1) in the partial components of the structure factors. In physical terms, this is equivalent to overestimated spatial nanoscale heterogeneity. To correct these problems, we adjusted the chloride-hydroxyl radial distribution functions against ab initio data and then extended the use of the Tang-Toennis damping function for the chlorides' induced dipoles. These adjustments correct the problems without losing the robustness of the CL&Pol force field. The results were also compared with the nonpolarizable version, the CL&P force field. We expect that the corrections will facilitate reliable use of the CL&Pol force field for other types of DESs.