Building better water models using the shape of the charge distribution of a water molecule

Exploring optimization strategies for improving explicit water models: Rigid n-point model and polarizable model based on Drude oscillator

Rigid n-point water models are widely used in atomistic simulations, but have known accuracy drawbacks. Increasing the number of point charges, as well as adding electronic polarizability, are two common strategies for accuracy improvements. Both strategies come at considerable computational cost, which weighs heavily against modest possible accuracy improvements in practical simulations. In an effort to provide guidance for model development, here we have explored the limiting accuracy of "electrostatically globally optimal" n-point water models in terms of their ability to reproduce properties of water dimer-a mimic of the condensed state of water. For a given n, each model is built upon a set of reference multipole moments (e.g. ab initio) and then optimized to reproduce water dimer total dipole moment. The models are then evaluated with respect to the accuracy of reproducing the geometry of the water dimer. We find that global optimization of the charge distribution alone can deliver high accuracy of the water model: for n = 4 or n = 5, the geometry of the resulting water dimer can be almost within 50 of the ab initio reference, which is half that of the experimental error margin. Thus, global optimization of the charge distribution of classical n-point water models can lead to high accuracy models. We also find that while the accuracy improvement in going from n = 3 to n = 4 is substantial, the additional accuracy increase in going from n = 4 to n = 5 is marginal. Next, we have explored accuracy limitations of the standard practice of adding electronic polarizability (via a Drude particle) to a "rigid base"-pre-optimization rigid n-point water model. The resulting model (n = 3) shows a relatively small improvement in accuracy, suggesting that the strategy of merely adding the polarizability to an inferior accuracy water model used as the base cannot fix the defects of the latter. An alternative strategy in which the parameters of the rigid base model are globally optimized along with the polarizability parameter is much more promising: the resulting 3-point polarizable model out-performs even the 5-point optimal rigid model by a large margin. We suggest that future development efforts consider 3- and 4-point polarizable models where global optimization of the "rigid base" is coupled to optimization of the polarizability to deliver globally optimal solutions.

Ab Initio Molecular Dynamics Investigation of the Electronic and Structural Stability of Anionic O2-(H2O)n, n = 1-16 Clusters

We report an ab initio molecular dynamics investigation of the electronic and structural stability of water molecules binding to a nucleation O₂^- particle, O₂^-(H₂O)_n with n = 1-16, to ascertain the factors that create particularly stable species. Our results compare well with previous experimental and theoretical reports for clusters with less water content, find three new geometries for species with 7, 9, and 10 water molecules, and determine that 8, 11, 13, and 15 water molecules form remarkably stable structures around O₂^-. These special clusters correspond to well-defined compact structures formed by cubes and four-member rings made of water's hydrogen bonds interacting with a negative kernel formed by O₂^- with five water molecules, O₂^-(H₂O)₅, in which the negative charge is localized in the first four water molecules, while the fifth molecule provides geometrical stability. We assess the clusters' energetic stability based on dissociation energies, analyze electron detachment energies to understand its geometrical evolution, and investigate its charge distribution based upon isosurfaces of the highest occupied molecular orbital (HOMO). This research can help provide theoretical insight into the starting steps of nucleation of water clusters around ionic particles.