Molecular modeling of aqueous electrolytes at interfaces: effects of long-range dispersion forces and of ionic charge rescaling

Abstract

Molecular dynamics simulations of aqueous electrolytes generally rely on empirical force fields, combining dispersion interactions - described by a truncated Lennard-Jones (LJ) potential - and electrostatic interactions - described by a Coulomb potential computed with a long-range solver. Recently, force fields using rescaled ionic charges (electronic continuum correction, ECC), possibly complemented with rescaling of LJ parameters (electronic continuum correction rescaled, ECCR), have shown promising results in bulk, but their performance at interfaces has been less explored. Here we started by exploring the impact of the LJ potential truncation on the surface tension of a sodium chloride aqueous solution. We show a discrepancy between the numerical predictions for truncated LJ interactions with a large cutoff and for untruncated LJ interactions computed with a long-range solver, which can bias comparison of force field predictions with experiments. Using a long-range solver for LJ interactions, we then show that an ionic charge rescaling factor chosen to correct long-range electrostatic interactions in bulk also describes accurately image charge repulsion at the liquid-vapor interface, and that the rescaling of LJ parameters in ECCR models - aimed at capturing local ion-ion and ion-water interactions in bulk - also describes well the formation of an ionic double layer at the liquid-vapor interface. Overall, these results suggest that the molecular modeling of aqueous electrolytes at interfaces would benefit from using long-range solvers for dispersion forces, and from using ECCR models, where the charge rescaling factor should be chosen to correct long-range electrostatic interactions.

arXiv commentMain text: 6 pages, 3 figures; Supplemental material: 9 pages, 3 figures; to be published in J. Chem. Phys

References in corpus (12)

• How to mesh up Ewald sums (I): A theoretical and numerical comparison of various particle mesh routines
• Accounting for electronic polarization in nonpolarizable force fields
• Understanding water's anomalies with locally favored structures
• A force field of Li${}^{+}$ , Na${}^{+}$ , K${}^{+}$, Mg${}^{2+}$, Ca${}^{2+}$, Cl${}^{-}$, and SO${}_4^{2-}$ in aqueous solution based on the TIP4P/2005 water model and scaled charges for the ions
• Anomalies in water as obtained from computer simulations of the TIP4P/2005 model: density maxima, and density, isothermal compressibility and heat capacity minima
• Development and application of a particle-particle particle-mesh Ewald method for dispersion interactions
• Two-structure thermodynamics for the TIP4P/2005 model of water covering supercooled and deeply stretched regions
• Recent Progress in Molecular Simulation of Aqueous Electrolytes: Force Fields, Chemical Potentials and Solubility
• Viscosity and self-diffusion of supercooled and stretched water from molecular dynamics simulations
• Water: one molecule, two surfaces, one mistake
• Critical Drying of Liquids
• Communication: Truncated non-bonded potentials can yield unphysical behavior in molecular dynamics simulations of interfaces

Cited by in corpus (9)

• Molecular Simulation of Electrode-Solution Interfaces
• Scaled charges for ions: an improvement but not the final word for modeling electrolytes in water
• Computation of Electrical Conductivities of Aqueous Electrolyte Solutions: Two Surfaces , One Property
• In Silico Seawater
• How important is the dielectric constant in water modeling? Evaluation of the performance of the TIP4P/$\epsilon$ force field and its compatibility with the Joung-Cheatham NaCl model
• Impact of the sodium and calcium chlorides uptake on the interfacial behavior of ice: premelting, structure, and dynamics
• Effect of Succinonitrile on Ion Transport in PEO-based Lithium Ion Battery Electrolytes
• Molecular insight on ultra-confined ionic transport in wetting films: the key role of friction
• Molecular dynamics study of perchloric acid using the extended Madrid-2019 force field