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

Abstract

In this contribution, we compute the ${}^1$H nuclear magnetic resonance (NMR) relaxation rate of liquid water at ambient conditions. We are using structural and dynamical information from Coupled Cluster Molecular Dynamics (CCMD) trajectories generated at CCSD(T) electronic structure accuracy while considering also nuclear quantum effects in addition to consulting information from X-ray and neutron scattering experiments. Our analysis is based on a recently presented computational framework for determining the frequency-dependent NMR dipole-dipole relaxation rate of spin $1/2$ nuclei from Molecular Dynamics (MD) simulations, which allows for an effective disentanglement of its structural and dynamical contributions, and is including a correction for finite-size effects inherent to MD simulations with periodic boundary conditions. A close to perfect agreement with experimental relaxation data is achieved if structural and dynamical informations from CCMD trajectories are considered including a re-balancing of the rotational and translational dynamics, according to the product of the self-diffusion coefficient and the reorientational correlation time of the H-H vector $D_0\times {\tau }_{\mathrm{HH}}$. The simulations show that this balance is significantly altered when nuclear quantum effects are taken into account. Our analysis suggests that the intermolecular and intramolecular contribution to the ${}^1$H NMR relaxation rate of liquid water are almost similar in magnitude, unlike to what was predicted earlier from classical MD simulations.

13 pages, 8 figures, added references, corrected authors list in metadata

References in corpus (10)

• Efficient stochastic thermostatting of path integral molecular dynamics
• 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
• High-dimensional neural network potentials for solvation: The case of protonated water clusters in helium
• Quadrupolar ${}^{23}$Na${}^{+}$ NMR Relaxation as a Probe of Subpicosecond Collective Dynamics in Aqueous Electrolyte Solutions
• Nuclear Quantum Effects in Liquid Water Are Negligible for Structure but Significant for Dynamics
• OrthoBoXY: A Simple Way to Compute True Self-Diffusion Coefficients from MD Simulations with Periodic Boundary Conditions Without Prior Knowledge of the Viscosity
• Accurate diffusion coefficients of the excess proton and hydroxide in water via extensive ab initio simulations with different schemes
• NMR relaxation rates of quadrupolar aqueous ions from classical molecular dynamics using force-field specific Sternheimer factors
• Computing the Frequency-Dependent NMR Relaxation of ${}^1$H Nuclei in Liquid Water
• An OrthoBoXY-Method for Various Alternative Box Geometries