Abstract: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.

Building Water Models, A Different Approach

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

A New Coarse-Grained Model for Water: the Importance of Electrostatic Interactions

Building Water Models Compatible with Charge Scaling Molecular Dynamics

Parameterization of Monovalent Ions for the OPC3, OPC, TIP3P-FB, and TIP4P-FB Water Models

Bottom-up Derived Flexible Water Model with Dipole and Quadrupole Moments for Coarse-Grained Molecular Simulations.

Molecular Simulation of Water and Hydration Effects in Different Environments: Challenges and Developments for Dftb Based Models

A QM/MM Derived Polarizable Water Model for Molecular Simulation

Deconstructing Classical Water Models at Interfaces and in Bulk

Competing quantum effects in the dynamics of a flexible water model

Simulating a flexible water model as rigid: Best practices and lessons learned

On the transferability of three water models developed by adaptive force matching

Simple and Accurate One-Body Energy and Dipole Moment Surfaces for Water and Beyond

Assessment of an Electrostatic Energy-Based Charge Model for Modeling the Electrostatic Interactions in Water Solvent.

An Alternative Model for Simulating Water Between Two Monolayer Surfaces

Modeling Molecular Interactions in Water: From Pairwise to Many-Body Potential Energy Functions

A Three-Point Coarse-Grained Model of Five-Water Cluster with Permanent Dipoles and Quadrupoles

A Transferrable Range-Separated Force Field for Water: Combining the Power of Both Physically-Motivated Models and Machine Learning Techniques

q-AQUA: a many-body CCSD(T) water potential, including 4-body interactions, demonstrates the quantum nature of water from clusters to the liquid phase

Seven-Site Effective Pair Potential for Simulating Liquid Water

Three‐site and five‐site fixed‐charge water models compatible with AMOEBA force field