Abstract:Tröster et al. (J. Phys. Chem B 2013, 117, 9486-9500) recently suggested a mixed computational and empirical approach to the optimization of polarizable molecular mechanics (PMM) water models. In the empirical part the parameters of Buckingham potentials are optimized by PMM molecular dynamics (MD) simulations. The computational part applies hybrid calculations, which combine the quantum mechanical description of a H2O molecule by density functional theory (DFT) with a PMM model of its liquid phase environment generated by MD. While the static dipole moments and polarizabilities of the PMM water models are fixed at the experimental gas phase values, the DFT/PMM calculations are employed to optimize the remaining electrostatic properties. These properties cover the width of a Gaussian inducible dipole positioned at the oxygen and the locations of massless negative charge points within the molecule (the positive charges are attached to the hydrogens). The authors considered the cases of one and two negative charges rendering the PMM four- and five-point models TL4P and TL5P. Here we extend their approach to three negative charges, thus suggesting the PMM six-point model TL6P. As compared to the predecessors and to other PMM models, which also exhibit partial charges at fixed positions, TL6P turned out to predict all studied properties of liquid water at p0 = 1 bar and T0 = 300 K with a remarkable accuracy. These properties cover, for instance, the diffusion constant, viscosity, isobaric heat capacity, isothermal compressibility, dielectric constant, density, and the isobaric thermal expansion coefficient. This success concurrently provides a microscopic physical explanation of corresponding shortcomings of previous models. It uniquely assigns the failures of previous models to substantial inaccuracies in the description of the higher electrostatic multipole moments of liquid phase water molecules. Resulting favorable properties concerning the transferability to other temperatures and conditions like the melting of ice are also discussed.

Inclusion of Water Multipoles into the Implicit Solvation Framework Leads to Accuracy Gains

Including Explicit Water Molecules as Part of the Protein Structure in MM/PBSA Calculations.

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

Molecular Theories and Simulation of Ions and Polar Molecules in Water

The AGBNP2 Implicit Solvation Model

Completely Multipolar Model for Many-Body Water-Ion and Ion-Ion Interactions

Implicit Solvation Using the Superposition Approximation (IS-SPA): Extension to Polar Solutes in Chloroform

Completely Multipolar Model as a General Framework for Many-Body Interactions as Illustrated for Water

Polarizable Continuum Models and Green's Function GW Formalism: On the Dynamics of the Solvent Electrons

Polarizable six-point water models from computational and empirical optimization

Accurate Calculation of Many-Body Energies in Water Clusters Using a Classical Geometry-Dependent Induction Model

Transferable Ion Force Fields in Water from a Simultaneous Optimization of Ion Solvation and Ion–Ion Interaction

Building Water Models, A Different Approach

Polarizable Continuum Models and Green's Function $\bf{GW}$ Formalism: On the Dynamics of the Solvent Electrons

Comparison of permutationally invariant polynomials, neural networks, and Gaussian approximation potentials in representing water interactions through many-body expansions

Improving the Efficiency of Non-equilibrium Sampling in the Aqueous Environment via Implicit-Solvent Simulations.

Interplay of local hydrogen-bonding and long-ranged dipolar forces in simulations of confined water

Condensed Phase Water Molecular Multipole Moments from Deep Neural Network Models Trained on Ab Initio Simulation Data

PW-SMD: A Plane-Wave Implicit Solvation Model Based on Electron Density for Surface Chemistry and Crystalline Systems in Aqueous Solution

Development of an Advanced Force Field for Water Using Variational Energy Decomposition Analysis

Ion clusters and ion-water potentials in MD-simulations