Accurate diffusion coefficients of the excess proton and hydroxide in water via extensive ab initio simulations with different schemes

Abstract

[ This article may be downloaded for personal use only. Any other use requires prior permission of the author and AIP Publishing. This article appeared in The Journal of Chemical Physics, 157, 024504 (2022) and may be found at: https://doi.org/10.1063/5.0093958 ] Despite its simple molecular formula, obtaining an accurate in silico description of water is far from straightforward. Many of its very peculiar properties are quite elusive, and in particular obtaining good estimations of the diffusion coefficients of the solvated proton and hydroxide at a reasonable computational cost has been an unsolved challenge until now. Here I present extensive results of several unusually long ab initio MD simulations employing different combinations of the Born--Oppenheimer and 2${}^{\mathrm{nd}}$-generation Car--Parrinello MD propagation methods with different ensembles (NVE and NVT) and thermostats, which show that these methods together with the RPBE-D3 functional provide a very accurate estimation of the diffusion coefficients of the solvated H${}_3$O${}^{+}$ and OH${}^{-}$ ions, together with an extremely accurate description of several properties of neutral water (such as the structure of the liquid and its diffusion and shear viscosity coefficients). In addition, I show that the estimations of $D_{{\mathrm{H}}_3{\mathrm{O}}^{+}}$ and $D_{\mathrm{O}{\mathrm{H}}^{-}}$ depend dramatically on the simulation length, being necessary to reach timescales in the order of hundreds of picoseconds to obtain reliable results.

arXiv comment18 pages including Supporting Information

References in corpus (19)

• Canonical sampling through velocity-rescaling
• Separable Dual Space Gaussian Pseudo-potentials
• Van der Waals Density Functional for General Geometries
• Relativistic separable dual-space Gaussian Pseudopotentials from H to Rn
• CP2K: An Electronic Structure and Molecular Dynamics Software Package -- Quickstep: Efficient and Accurate Electronic Structure Calculations
• Perspective: How good is DFT for water?
• 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
• The Individual and Collective Effects of Exact Exchange and Dispersion Interactions on the Ab Initio Structure of Liquid Water
• Why does hydronium diffuse faster than hydroxide in liquid water?
• Density, structure and dynamics of water: the effect of Van der Waals interactions
• Ab-Initio Molecular Dynamics
• Network equilibration and first-principles liquid water
• Stochastic thermostats: comparison of local and global schemes
• Room temperature compressibility and diffusivity of liquid water from first principles
• Ab-initio Structure and Thermodynamics of the RPBE-D3 Water/Vapor Interface by Neural-Network Molecular Dynamics
• Decoding the spectroscopic features and timescales of aqueous proton defects
• The Hopping Mechanism of the Hydrated Excess Proton and Its Contribution to Proton Diffusion in Water

Cited by in corpus (2)

• Mechanistic insights into water autoionization
• When Theory Meets Experiment: What Does it Take to Accurately Predict ${}^1$H NMR Dipolar Relaxation Rates in Neat Liquid Water from Theory?