We investigate the dissociation of a Ca -Cl pair in water using classical molecular dynamics simulations with a polarizable interaction potential, parameterized from ab initio calculations. By computing the potential of mean force as a function not only of the interionic distance but also of the coordination numbers by water molecules, we show that it is necessary to use a collective variable describing the cation hydration in order to capture the dissociation mechanism. In the contact ion pair, the Ca cation has a first coordination sphere containing 5 or 6 water molecules.
View Article and Find Full Text PDFWe report the diffusion coefficient and viscosity of popular rigid water models: two non-polarizable ones (SPC/E with three sites, and TIP4P/2005 with four sites) and a polarizable one (Dang-Chang, four sites). We exploit the dependence of the diffusion coefficient on the system size (Yeh and Hummer 2004 J. Phys.
View Article and Find Full Text PDFWe present a new polarizable force field for aqueous ions (Li(+), Na(+), K(+), Rb(+), Cs(+), Mg(2 +), Ca(2 +), Sr(2 +), and Cl(-)) derived from condensed phase ab initio calculations. We use maximally localized Wannier functions together with a generalized force and dipole-matching procedure to determine the whole set of parameters. Experimental data are then used only for validation purposes and a good agreement is obtained for structural, dynamic, and thermodynamic properties.
View Article and Find Full Text PDFDipole polarizabilities of a series of ions in aqueous solutions are computed from first-principles. The procedure is based on the study of the linear response of the maximally localized Wannier functions to an applied external field, within density functional theory. For most monoatomic cations (Li(+), Na(+), K(+), Rb(+), Mg(2+), Ca(2+) and Sr(2+)) the computed polarizabilities are the same as in the gas phase.
View Article and Find Full Text PDFThe potential-driven ordering transition of a LiCl layer adsorbed on the (100) surface of a metallic aluminum electrode is studied by molecular dynamics simulations. The transition causes a sharp peak in the potential dependence of the differential capacitance of the interface. This result is in qualitative agreement with recently reported experimental work on the interface between a room temperature ionic liquid and a well-defined Au(100) surface.
View Article and Find Full Text PDF