Ca2+ -Cl- Association in Water Revisited: the Role of Cation Hydration.

We investigate the dissociation of a Ca2+ -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 Ca2+ cation has a first coordination sphere containing 5 or 6 water molecules. The minimum free-energy path for dissociation involves a two-step process: First one or two additional water molecules enter the cation coordination shell, increasing the coordination number up to 7 with an almost fixed interionic distance. Then the dissociation of the ionic pair occurs at this fixed coordination number.

Diffusion coefficient and shear viscosity of rigid water models.

We 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. Chem. B 108 15873) to obtain the size-independent value. This also provides an estimate of the viscosity of all water models, which we compare to the Green-Kubo result. In all cases, a good agreement is found. The TIP4P/2005 model is in better agreement with the experimental data for both diffusion and viscosity. The SPC/E and Dang-Chang models overestimate the diffusion coefficient and underestimate the viscosity.

A transferable ab initio based force field for aqueous ions.

We 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. The same procedure applied to crystalline phases allows to parametrize the interaction between cations and the chloride anion. Finally, we illustrate the good transferability of the force field to other thermodynamic conditions by investigating concentrated solutions.

Ions in solutions: Determining their polarizabilities from first-principles.

Dipole 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. For Cs(+) and a series of anions (F(-), Cl(-), Br(-) and I(-)), environmental effects are observed, which reduce the polarizabilities in aqueous solutions with respect to their gas phase values. The polarizabilities of H((aq)) (+), OH((aq)) (-) have also been determined along an ab initio molecular dynamics simulation. We observe that the polarizability of a molecule instantaneously switches upon proton transfer events. Finally, we also computed the polarizability tensor in the case of a strongly anisotropic molecular ion, UO(2) (2+). The results of these calculations will be useful in building interaction potentials that include polarization effects.

Potential-induced ordering transition of the adsorbed layer at the ionic liquid/electrified metal interface.

The 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. In the LiCl/Al simulations, the transition occurs when the interaction model includes the induction of dipoles on the ions of the liquid by their mutual interaction and their interaction with the electrode surface as well as the polarization of the metal by the charges and dipoles on the ions ("image" interactions). When the electrode or ion polarization effects are not included, the transition is no longer observed. The interaction between the induced charges on the metal atoms and the induced dipoles on the ions creates an additional screening, which stabilizes the formation of a crystalline layer at the anode. When the crystallographic plane of the metal is changed to (110) instead of (100), the two first adsorbed layer are crystalline on both the anode and the cathode, but the structure is different: the crystal is formed through an epitaxial mechanism to adapt to the electrode surface structure. In the case of the (110) crystallographic plane, the charging of the adsorbed layer occurs through the formation of nonstoichiometric crystalline layers.