Polarizable Water Model with Ab Initio Neural Network Dynamic Charges and Spontaneous Charge Transfer

Abstract

Simulating water accurately has been a challenge due to the complexity of describing polarization and intermolecular charge transfer. Quantum mechanical (QM) electronic structures provide an accurate description of polarization in response to local environments, which is nevertheless too expensive for large water systems. In this study, we have developed a polarizable water model integrating Charge Model 5 atomic charges at the level of the second-order Mo̷ller-Plesset perturbation theory, predicted by an accurate and transferable charge neural network (ChargeNN) model. The spontaneous intermolecular charge transfer has been explicitly accounted for, enabling a precise treatment of hydrogen bonds and out-of-plane polarization. Our ChargeNN water model successfully reproduces various properties of water in gas, liquid, and solid phases. For example, ChargeNN correctly captures the hydrogen-bond stretching peak and bending-libration combination band, which are absent in the spectra using fixed charges, highlighting the significance of accurate polarization and charge transfer. Finally, the molecular dynamical simulations using ChargeNN for liquid water and a large water droplet with a ∼4.5 nm radius reveal that the strong interfacial electric fields are concurrently induced by the partial collapse of the hydrogen-bond network and surface-to-interior charge transfer. Our study paves the way for QM-polarizable force fields, aiming for large-scale molecular simulations with high accuracy.