A comparison of classical interatomic potentials applied to highly concentrated aqueous lithium chloride solutions

Abstract

Aqueous lithium chloride solutions up to very high concentrations were investigated in classical molecular dynamics simulations. Various force fields based on the 12-6 Lennard-Jones model, parametrized for non-polarizable water solvent molecules (SPC/E, TIP4P, TIP4PEw), were inspected. Twenty-nine combinations of ion-water interaction models were examined at four different salt concentrations. Densities, static dielectric constants and self-diffusion coefficients were calculated. Results derived from the different force fields scatter over a wide range of values. Neutron and X-ray weighted structure factors were also calculated from the radial distribution functions and compared with experimental data. It was found that the agreement between calculated and experimental curves is rather poor for several investigated potential models, even though some of them have previously been applied in computer simulations. None of the investigated models yield satisfactory results for all the tested quantities. Only two parameter sets provide acceptable predictions for the structure of highly concentrated aqueous LiCl solutions. Some approaches for adjusting potential parameters, such as those of Aragones [Aragones et al., J. Phys. Chem. B 118 (2014) 7680] and Pluharova [Pluharova et al, J. Phys. Chem. A 117 (2013) 11766], were tested as well; the simulations presented here underline their usefulness. These refining methods are suited to obtain more appropriate ion/water potentials.

65 pages together with Supplementary informations, submitted to J.Mol.Liq

References in corpus (7)

• Dependence of the dielectric constant of electrolyte solutions on ionic concentration - a microfield approach
• On the accurate calculation of the dielectric constant and the diffusion coefficient from molecular dynamics simulations: the case of SPC/E water
• Recent Progress in Molecular Simulation of Aqueous Electrolytes: Force Fields, Chemical Potentials and Solubility
• 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
• Reverse Monte Carlo investigations concerning recent isotopic substitution neutron diffraction data on liquid water
• Ion clusters and ion-water potentials in MD-simulations

Cited by in corpus (8)

• 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
• The Madrid-2019 force field for electrolytes in water using TIP4P/2005 and scaled charges: extension to the ions F${}^{-}$, Br${}^{-}$, I${}^{-}$, Rb${}^{+}$, Cs${}^{+}$
• 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
• Chloride ions as integral parts of hydrogen bonded networks in aqueous salt solutions: the appearance of solvent separated anion pairs
• The structure of aqueous lithium chloride solutions at high concentrations as revealed by a comparison of classical interatomic potential models
• Experimental and molecular dynamics study of the ionic conductivity in aqueous LiCl electrolytes
• Towards the correct microscopic structure of aqueous CsCl solutions with a comparison of classical interatomic potential models