We ran molecular dynamics (MD) and Monte Carlo (MC) simulations using the RASPA software. We used a combination of Lennard-Jones and electrostatic potentials. Water molecules are represented with the rigid TIP4P model and the electrolytes are modeled as charged particles that differs by their Lennard-Jones parameters and charges. We considered the following monovalent electrolytes which are: LiCl, NaCl, KCl, RbCl, CsCl, KF, KBr and MgCl 2 , CaCl 2 , SrCl 2 and BaCl 2 . The details of the force field are given in the SI. We used Lorentz-Berthelot mixing rules for cross-terms in Lennard-Jones potential, and Ewald summation was used to account efficiently for electrostatic interactions. We used a cutoff of 8.0 Å for both the Lennard-Jones potential and the separation between real space and Fourier space in the Ewald summation. MD simulations were performed with a timestep of 1 fs and all simulations used a Nosé-Hoover thermostat, setting the temperature at 300 K.