Molecular Modelling of Ionic Liquids: General Guidelines on Fixed-Charge Force Fields for Balanced Descriptions

30 June 2022, Version 2
This content is a preprint and has not undergone peer review at the time of posting.

Abstract

It has been increasingly common to investigate dynamic and thermodynamic properties of green solvents at atomistic scales with molecular simulation. These designed solvents such as ionic liquids are often highly charged species, which pose a problem for molecular modelling with classical fixed-charge force fields. Simulation outcomes with atomic charges derived from ab initio calculations often show significant deviations from the experimental reference, and charge scaling is widely applied as a simple yet efficient solution to achieve satisfactory experiment-simulation accordance. Although the massive emphasis on the reproduction of bulk properties such as density of the liquids as charge scaling criteria, a more crucial thermodynamic observable is the solvation behavior of external agents in these green solvents. Astonishingly, our recent large-scale benchmark calculation on solvation thermodynamics suggests that the solvation-free-energy-derived scaling factor is generally slightly larger (~0.1) than the bulk-property-derived estimate. This phenomenon is rather not unexpected, as the density-matching estimate only considers the solvent-solvent interaction (overfitting), while accurate calculations of solvation free energies require balanced descriptions of solute-solvent and solvent-solvent interactions. A more interesting observation is that different solute-solvent pairs exhibit different responses to the variation of the scaling parameter, which arises from the competing electrostatic and vdW contributions. Another perspective provided in our previous extensive benchmark is about the suitability of general-purpose force fields for bonded and vdW interactions in the modelling of ionic liquids. The bond stretching and angle bending terms in pre-fitted GAFF derivatives are often problematic, while the torsional potential shows satisfactory reproduction of ab initio results. In the current work, we expect to accumulate more experiences from large-scale fast-growth solvation free energy simulations in ionic liquids with compositions different from our previous benchmark. The obtained new results are combined with our previous dataset to form a large solvation set (and also partition or water-ionic-liquids transfer set), from which a universally applicable charge scaling factor with at least half-optimal performance is derived. Aside from the force-field issue, another extremely important modelling detail considered in the current work is the finite-size effect. It is observed that the finite-size artifacts in solvation thermodynamics are much more severe than mass density, which emphasizes the use of a sufficiently large ionic solvent box in molecular simulations of ionic liquids derivatives. Finally, combining the extensive computational perspectives accumulated in our series works, general guidelines for molecular modelling of ionic liquids with fixed-charge force fields are summarized.

Keywords

Ionic Liquids
Free Energy Calculation
GAFF
Fast Growth
Force Field
Solvation Free Energy
Charge Scaling
Electrostatic Potential
Finite-Size Effects
ab initio Calculation
Property-matching
Force Field Refitting
Bond Stretching
Angle Bending
Torsional
Large-scale Calculation
Partition Coefficient
Mass Density
Water-ionic-liquids Transfer Free Energy
Restrained Electrostatic Potential
AMBER
GROMACS
Slow Growth
Gaussian Approximation
1-ethyl-3-methylimidazolium
dicyanamide
bis(trifluoromethylsulfonyl)imide
trifluoromethylsulfonate
Imidazolium
B97-3c
PACKMOL
ORCA
Nonequilibrium Alchemical Transformation
Dipole Moments
Error Metrics
Ranking Coefficients
Force-matching
Density-matching
Water-Ionic-Liquids Partition
BP86
CAM-B3LYP
B3LYP
M05-2X
M06-2X
MN15
PW6B95
TPSSh
wB97X-D
MP2
B2PLYP
mPW2PLYP
HF
def2-TZVPP
ESP RRMSE

Comments

Comments are not moderated before they are posted, but they can be removed by the site moderators if they are found to be in contravention of our Commenting Policy [opens in a new tab] - please read this policy before you post. Comments should be used for scholarly discussion of the content in question. You can find more information about how to use the commenting feature here [opens in a new tab] .
This site is protected by reCAPTCHA and the Google Privacy Policy [opens in a new tab] and Terms of Service [opens in a new tab] apply.