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

Abstract

In this work, a force field for several ions in water is proposed. In particular, we consider the cations Li${}^{+}$ , Na${}^{+}$ , K${}^{+}$, Mg${}^{2+}$, Ca${}^{2+}$ and the anions Cl${}^{-}$, and SO${}_4^{2-}$ . These ions were selected as they appear in the composition of seawater, and they are also found in biological systems. The force field proposed (denoted as Madrid-2019) is nonpolarizable, and both water molecules and sulfate anions are rigid. For water, we use the TIP4P/2005 model. The main idea behind this work is to further explore the possibility of using scaled charges for describing ionic solutions. Monovalent and divalent ions are modeled using charges of 0.85 and 1.7, respectively (in electron units). The model allows a very accurate description of the densities of the solutions up to high concentrations. It also gives good predictions of viscosities up to 3 m concentrations. Calculated structural properties are also in reasonable agreement with the experiment. We have checked that no crystallization occurred in the simulations at concentrations similar to the solubility limit. A test for ternary mixtures shows that the force field provides excellent performance at an affordable computer cost. In summary, the use of scaled charges, which could be regarded as an effective and simple way of accounting for polarization (at least to a certain extend), improves the overall description of ionic systems in water. However, for purely ionic systems, scaled charges will not adequately describe neither the solid nor the melt.

References in corpus (9)

• Canonical sampling through velocity-rescaling
• What ice can teach us about water interactions: a critical comparison of the performance of different water models
• A route to explain water anomalies from results on an aqueous solution of salt
• Polarizable Mean-Field Model of Water for Biological Simulations with Amber and Charmm force fields
• Crystal lattice properties fully determine short-range interaction parameters for alkali and halide ions
• Polarizable molecular interactions in condensed phase and their equivalent nonpolarizable models
• Water: one molecule, two surfaces, one mistake
• A comparison of classical interatomic potentials applied to highly concentrated aqueous lithium chloride solutions
• Stabilization of AgI's polar surfaces by the aqueous environment, and its implications for ice formation

Cited by in corpus (8)

• The Madrid-2019 force field for electrolytes in water using TIP4P/2005 and scaled charges: extension to the ions F${}^{-}$, Br${}^{-}$, I${}^{-}$, Rb${}^{+}$, Cs${}^{+}$
• A review on laser-induced crystallization from solution
• Computation of Electrical Conductivities of Aqueous Electrolyte Solutions: Two Surfaces , One Property
• How Grain Boundaries and Interfacial Electrostatic Interactions Modulate Water Desalination Via Nanoporous Hexagonal Boron Nitride
• Understanding the Anomalous Diffusion of Water in Aqueous Electrolytes Using Machine Learned Potentials
• Three phase equilibria of the methane hydrate in NaCl solutions: A simulation study
• In Silico Seawater
• NMR relaxation rates of quadrupolar aqueous ions from classical molecular dynamics using force-field specific Sternheimer factors