Polarization response of water and methanol investigated by a polarizable force field and density functional theory...

Electronic polarization response in hydrogen-bond clusters and liquid configurations of water and methanol has been studied by density functional theory (DFT) and by a polarizable force field based on the chemical potential equalization (CPE) principle. It has been shown that an accurate CPE parametrization based on isolated molecular properties is not completely transferable to strongly interacting hydrogen-bond clusters with discrepancies between CPE and DFT overall dipole moments as large as 15%. This is due to the lack of intermolecular charge transfer in the standard CPE implementation. A CPE scheme for evaluating the amount of transferred charge has been developed. The charge transfer parameters are determined with the aid of accurate DFT calculations using only hydrogen-bond dimer configurations. The amount of transferred charge is found to be of the order of few hundredths of electrons, as already found in recent studies on hydrogen-bond systems. The parameters of the model are then used, without further adjustment, to different hydrogen-bond clustered forms of water and methanol (oligomer and liquid configurations). In agreement with different approaches proposed in literature for studying charge transfer effects, the transferred charge in hydrogen-bond dimers is found to decrease exponentially with the hydrogen-bond distance. When allowance is made for charge transfer according to the proposed scheme, the CPE dipole moments are found to reproduce satisfactorily the DFT data.