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.

Mean-Field Description Of Ionic Size Effects With Nonuniform Ionic Sizes: A Numerical Approach

Effect of orientational ordering of water dipoles on stratification of counterions of different size in multicomponent electrolyte solution near charged surface -- a mean field approach

Molecular Mean-Field Theory of Ionic Solutions: a Poisson-Nernst-Planck-Bikerman Model

Ion Size Effects on Individual Fluxes Via Poisson-Nernst-Planck Systems with Bikerman's Local Hard-Sphere Potential: Analysis Without Electroneutrality Boundary Conditions

Finite ion size effects on I-V relations via Poisson-Nernst-Planck systems with two cations: A case study

Multi-scale theory in the molecular simulation of electrolyte solutions

A Molecular Simulation Study on the Role of Ion Sizes and Dielectric Images in Near-Surface Ion Distribution Far from the Strong Coupling Limit

Ion size effects on the electric double layer of a spherical particle in a realistic salt-free concentrated suspension

Classical density functional analysis of the ionic size effects on the properties of charge regulating electric double layers

Electrokinetic properties of NaCl solution via molecular dynamics simulations with scaled-charge electrolytes

Particles Inside Electrolytes with Ion-Specific Interactions, Their Effective Charge Distributions, and Effective Interactions

Ion Distributions, Exclusion Coefficients, and Separation Factors of Electrolytes in a Charged Cylindrical Nanopore: A Partially Perturbative Density Functional Theory Study

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

Density-functional Theory of Spherical Electric Double Layers and Zeta Potentials of Colloidal Particles in Restricted-Primitive-model Electrolyte Solutions.

Size Dependent Lithium-Ion Conductivity of Solid Electrolytes in Machine Learning Molcular Dynamics Simulations

Influence of solvent polarization and non-uniform ion size on electrostatic properties between charged surfaces in an electrolyte solution

Effective charges and virial pressure of concentrated macroion solutions

Scalable molecular simulation of electrolyte solutions with quantum chemical accuracy

Ionic Interactions in Biological and Physical Systems: a Variational Treatment

A local approximation of fundamental measure theory incorporated into three dimensional Poisson-Nernst-Planck equations to account for hard sphere repulsion among ions

Effective Electrostatic Interactions in Colloid-Nanoparticle Mixtures