Developing machine-learned potentials to simultaneously capture the dynamics of excess protons and hydroxide ions in classical and path integral simulations

Abstract

The transport of excess protons and hydroxide ions in water underlies numerous important chemical and biological processes. Accurately simulating the associated transport mechanisms ideally requires utilizing ab initio molecular dynamics simulations to model the bond breaking and formation involved in proton transfer and path-integral simulations to model the nuclear quantum effects relevant to light hydrogen atoms. These requirements result in a prohibitive computational cost, especially at the time and length scales needed to converge proton transport properties. Here, we present machine-learned potentials (MLPs) that can model both excess protons and hydroxide ions at the generalized gradient approximation and hybrid density functional theory levels of accuracy and use them to perform multiple nanoseconds of both classical and path-integral proton defect simulations at a fraction of the cost of the corresponding ab initio simulations. We show that the MLPs are able to reproduce ab initio trends and converge properties such as the diffusion coefficients of both excess protons and hydroxide ions. We use our multi-nanosecond simulations, which allow us to monitor large numbers of proton transfer events, to analyze the role of hypercoordination in the transport mechanism of the hydroxide ion and provide further evidence for the asymmetry in diffusion between excess protons and hydroxide ions.

References in corpus (27)

• Canonical sampling through velocity-rescaling
• Separable Dual Space Gaussian Pseudo-potentials
• CP2K: An Electronic Structure and Molecular Dynamics Software Package -- Quickstep: Efficient and Accurate Electronic Structure Calculations
• Perspective: How good is DFT for water?
• Competing quantum effects in the dynamics of a flexible water model
• How van der Waals interactions determine the unique properties of water
• Nuclear quantum effects enter the mainstream
• An Efficient and Accurate Car-Parrinello-like Approach to Born-Oppenheimer Molecular Dynamics
• Efficient stochastic thermostatting of path integral molecular dynamics
• i-PI 2.0: A Universal Force Engine for Advanced Molecular Simulations
• Ab initio thermodynamics of liquid and solid water
• The Individual and Collective Effects of Exact Exchange and Dispersion Interactions on the Ab Initio Structure of Liquid Water
• How to remove the spurious resonances from ring polymer molecular dynamics
• Why does hydronium diffuse faster than hydroxide in liquid water?
• 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
• i-PI: A Python interface for ab initio path integral molecular dynamics simulations
• Automated Fitting of Neural Network Potentials at Coupled Cluster Accuracy: Protonated Water Clusters as Testing Ground
• The Interplay of Structure and Dynamics in the Raman Spectrum of Liquid Water over the Full Frequency and Temperature Range
• 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
• Unraveling quantum mechanical effects in water using isotopic fractionation
• Ab initio molecular dynamics with nuclear quantum effects at classical cost: ring polymer contraction for density functional theory
• Effect of quantum nuclear motion on hydrogen bonding
• Isobaric first-principles molecular dynamics of liquid water with nonlocal van der Waals interactions
• Accurate molecular dynamics and nuclear quantum effects at low cost by multiple steps in real and imaginary time: using density functional theory to accelerate wavefunction methods
• Probing defects and correlations in the hydrogen-bond network of ab initio water
• Decoding the spectroscopic features and timescales of aqueous proton defects
• Unified description of hydrogen bonding and of proton transfer by a two-state effective Hamiltonian

Cited by in corpus (4)

• Entropy governs the structure and reactivity of water dissociation under electric fields
• Aqueous Solution Chemistry In Silico and the Role of Data Driven Approaches
• Self-learning path integral hybrid Monte Carlo with mixed ab initio and machine learning potentials for modeling nuclear quantum effects in water
• Active Learning of a Neural Network Potential for Large-Scale Atomistic Simulations of Polymer Electrolyte Membranes