The manuscript analyzes the accuracy of our recently developed reciprocal polarizable embedding scheme, where a density functional theory model of the QM region is coupled to a dipole- and quadrupole polarizable water potential of the MM region. We present calculations of water clusters and liquid water where we analyze the energy, atomic forces and total polarization to demonstrate that artifacts in energy and polarization introduced by the QM/MM coupling are small and well-behaved. Furthermore, our methodology improves the consistency of the structure of optimized water hexamer geometries when compared to results obtained with models that neglect polarization. Additionally, the manuscript provides evidence that our coupling scheme eliminates artifacts in the structure of liquid water obtained with simpler electrostatic embedding models.
Large vortex-like structure of dipole field in computer models of liquid water and dipole-bridge between biomolecules
