Dielectric constant and low-frequency infrared spectra for liquid water and ice Ih within the E3B model

Two intrinsic difficulties in modeling condensed-phase water with conventional rigid non-polarizable water models are: reproducing the static dielectric constants for liquid water and ice Ih, and generating the peak at about 200 cm(-1) in the low-frequency infrared spectrum for liquid water. The primary physical reason for these failures is believed to be the missing polarization effect in these models, and consequently various sophisticated polarizable water models have been developed. However, in this work we pursue a different strategy and propose a simple empirical scheme to include the polarization effect only on the dipole surface (without modifying a model's intermolecular interaction potential). We implement this strategy for our explicit three-body (E3B) model. Our calculated static dielectric constants and low-frequency infrared spectra are in good agreement with experiment for both liquid water and ice Ih over wide temperature ranges, albeit with one fitting parameter for each phase. The success of our modeling also suggests that thermal fluctuations about local minima and the energy differences between different proton-disordered configurations play minor roles in the static dielectric constant of ice Ih. Our analysis shows that the polarization effect is important in resolving the two difficulties mentioned above and sheds some light on the origin of several features in the low-frequency infrared spectra for liquid water and ice Ih.