NMR relaxation rates of quadrupolar aqueous ions from classical molecular dynamics using force-field specific Sternheimer factors

Abstract

The nuclear magnetic resonance (NMR) relaxation of quadrupolar nuclei is governed by the electric field gradient (EFG) fluctuations at their position. In classical molecular dynamics (MD), the electron cloud contribution to the EFG can be included via the Sternheimer approximation, in which the full EFG at the nucleus that can be computed using quantum DFT is considered to be proportional to that arising from the external, classical charge distribution. In this work, we systematically assess the quality of the Sternheimer approximation as well as the impact of the classical force field (FF) on the NMR relaxation rates of aqueous quadrupolar ions at infinite dilution. In particular, we compare the rates obtained using an ab initio parametrized polarizable FF, a recently developed empirical FF with scaled ionic charges and a simple empirical non-polarizable FF with formal ionic charges. Surprisingly, all three FFs considered yield good values for the rates of smaller and less polarizable solutes (Li + , Na + , K + , Cl − ), provided that a model-specific Sternheimer parametrization is employed. Yet,