Improving Condensed Phase Water Dynamics with Explicit Nuclear Quantum Effects: the Polarizable Q-AMOEBA Force Field

Abstract

We introduce a new parametrization of the AMOEBA polarizable force field for water denoted Q-AMOEBA, for use in simulations that explicitly account for nuclear quantum effects (NQEs). This study is made possible thanks to the recently introduced adaptive Quantum Thermal Bath (adQTB) simulation technique which computational cost is comparable to classical molecular dynamics. The flexible Q-AMOEBA model conserves the initial AMOEBA functional form, with an intermolecular potential including an atomic multipole description of electrostatic interactions (up to quadrupole), a polarization contribution based on the Thole interaction model and a buffered 14-7 potential to model van der Waals interactions. It has been obtained by using a Force Balance fitting strategy including high-level quantum chemistry reference energies and selected condensed phase properties targets. The final Q-AMOEBA model is shown to accurately reproduce both gas phase and condensed phase properties, notably improving the original AMOEBA water model. This development allows the fine study of NQEs on water liquid phase properties such as the average H-O-H angle compared to its gas phase equilibrium value, isotope effects etc... Q-AMOEBA also provides improved infrared spectroscopy prediction capabilities compared to AMOEBA03. Overall, we show that the impact of NQEs depends on the underlying model functional form and on the associated strength of hydrogen bonds. Since adQTB simulations can be performed at near classical computational cost using the Tinker-HP package, Q-AMOEBA can be extended to organic molecules, proteins and nucleic acids opening the possibility for the large scale study of the importance of NQEs in biophysics.

References in corpus (9)

• Competing quantum effects in the dynamics of a flexible water model
• What ice can teach us about water interactions: a critical comparison of the performance of different water models
• How to remove the spurious resonances from ring polymer molecular dynamics
• On the accuracy of the MB-pol many-body potential for water: Interaction energies, vibrational frequencies, and classical thermodynamic and dynamical properties from clusters to liquid water and ice
• Accelerating the convergence of path integral dynamics with a generalized Langevin equation
• Tinker-HP : Accelerating Molecular Dynamics Simulations of Large Complex Systems with Advanced Point Dipole Polarizable Force Fields using GPUs and Multi-GPUs systems
• Quantum contributions in the ice phases: the path to a new empirical model for water -- TIP4PQ/2005
• Development of the Quantum Inspired SIBFA Many-Body Polarizable Force Field: Enabling Condensed Phase Molecular Dynamics Simulations
• Nuclear Quantum Effects in liquid water at near classical computational cost using the adaptive Quantum Thermal Bath

Cited by in corpus (3)

• Force-Field-Enhanced Neural Network Interactions: from Local Equivariant Embedding to Atom-in-Molecule properties and long-range effects
• Routine Molecular Dynamics Simulations Including Nuclear Quantum Effects: from Force Fields to Machine Learning Potentials
• The Q-AMOEBA (CF) Polarizable Potential