The solvent water molecules are described by the polarizable model described by Chang (DC) , which is compatible with PIM. DC water is a rigid 4-site rigid model, with an additional virtual site M along the symmetry axis of the molecule, which carries a negative partial charge and the induced dipole. The short range repulsion and dispersion interactions are described by Lennard-Jones (LJ) interaction on the oxygen atom only. The parameters of the DC model are summarized in Table . The functional form of the water potential is included in the SI. This model has been shown to correctly reproduce the trends in several thermophysical properties of bulk water from ambient to high temperature conditions . The PIM potential parameters are determined sequentially. First, we derive the ionic polarizability (α I in Equation ) of La 3+ and Cl -ions. For La 3+ , DFT calculations are done with an isolated gas phase cation. The dipole induced on the cation by an electric field perturbation is analyzed by constructing maximally localized Wannier functions (MLWF) using the method by Molina et al. . Since the effect of the coordination environment on the polarizability of anions is significant , α I of Cl -ion is determined from analogous simulations of the solvated anion with box #3 in Table .