Exploring Ion Polarizabilities and Their Correlation with Van der Waals Radii: A Theoretical Investigation

Polarizability (⍺) is a fundamental property which measures the tendency of the electron cloud of an atom, ion, or molecule to be distorted by electric field. Polarizability contributes to important physical properties such as molecular interactions or dielectric constants; thus, it is essential to have accurate polarizabilities in molecular simulations. However, it remains a challenge to develop polarizable force fields for ions in computational chemistry. In particular, a comprehensive set of polarizabilities for ions has not been derived. Herein, we derived a systematic set of polarizabilities for atoms and ions across the periodic table based on high-level quantum mechanics calculations. These values have excellent agreement with experimental data. Furthermore, we examined the relationship between the obtained polarizabilities and the van der Waals radii (R_VDW) that we previously determined (J. Chem. Theory Comput., 2023, 19, 2064). Two relationships, R_VDW∝α^(1/7) and R_VDW∝α^(1/3), proposed in previous studies were examined in the present work. Our results indicated the former relationship, which was derived based on the quantum harmonic oscillator model, prevails for atoms and cations, but neither relationship provides a satisfactory fit for anions. This is consistent with the tight-binding nature of the electrons in atoms and cations, while it is more challenging to quantify the polarizabilities of anions because of their more dispersed electron clouds. Moreover, we compared different approaches to determine the dispersion coefficients, including the London equation, Slater-Kirkwood equation, SAPT calculations, and TD-DFT method, along with the approach based on van der Waals constants. Our results indicated that although different approaches predict deviated magnitudes for the dispersion coefficients, their predictions are highly correlated, implying that each of these approaches can be used to evaluate dispersion interactions after proper scaling. Finally, we have developed a parameterization strategy for the 12-6-4 model based on our study. We specifically compared the performance of the 12-6-4 model with SAPT and SobEDA analyses to model interactions involving Na+/Mg2+ and various ligands containing He, Ne, Ar, H2O, NH3, [H2PO4]−, and [HPO4]2−. The results demonstrate that the 12-6-4 parameters effectively replicate both the total interaction energy and the individual energy components (electrostatics, exchange-repulsion, dispersion, and induction), highlighting the physical relevance of the 12-6-4 model and the effectiveness of our parameterization approach. This study has significant implications for advancing the development of next-generation ion models and polarizable force fields.