Theory and modeling of molecular modes in the NMR relaxation of fluids.

Traditional theories of the nuclear magnetic resonance (NMR) autocorrelation function for intra-molecular dipole pairs assume a single-exponential decay, yet the calculated autocorrelation of realistic systems displays a rich, multi-exponential behavior, resulting in anomalous NMR relaxation dispersion (i.e., frequency dependence). We develop an approach to model and interpret the multi-exponential intra-molecular autocorrelation using simple, physical models within a rigorous statistical mechanical development that encompasses both rotational diffusion and translational diffusion in the same framework. We recast the problem of evaluating the autocorrelation in terms of averaging over a diffusion propagator whose evolution is described by a Fokker-Planck equation. The time-independent part admits an eigenfunction expansion, allowing us to write the propagator as a sum over modes. Each mode has a spatial part that depends on the specified eigenfunction and a temporal part that depends on the corresponding eigenvalue (i.e., correlation time) with a simple, exponential decay. The spatial part is a probability distribution of the dipole pair, analogous to the stationary states of a quantum harmonic oscillator. Drawing inspiration from the idea of inherent structures in liquids, we interpret each of the spatial contributions as a specific molecular mode. These modes can be used to model and predict the NMR dipole-dipole relaxation dispersion of fluids by incorporating phenomena on the molecular level. We validate our statistical mechanical description of the distribution in molecular modes with molecular dynamics simulations interpreted without any relaxation models or adjustable parameters: the most important poles in the Padé-Laplace transform of the simulated autocorrelation agree with the eigenvalues predicted by the theory.

Thermal and concentration effects on 1H NMR relaxation of Gd3+-aqua using MD simulations and measurements.

Gadolinium-based contrast agents are key in clinical MRI for enhancing the longitudinal NMR relativity (r1) of hydrogen nuclei (1H) in water and improving the contrast among different tissues. The importance of MRI in clinical practice cannot be gainsaid, yet the interpretation of MRI relies on models with severe assumptions, reflecting a poor understanding of the molecular-scale relaxation processes. In a step towards building a clearer understanding of the relaxation processes, here we investigate thermal and concentration effects on r1 of the Gd3+-aqua complex using both semi-classical molecular dynamics (MD) simulations and measurements. We follow the MD simulation approach recently introduced by [Singer et al., Phys. Chem. Chem. Phys., 2021, 23, 20974], in which no NMR relaxation model or free-parameter is assumed to predict r1, thereby bringing new insights into the physics of r1 on a molecular scale. We expand the autocorrelation function G(t) in terms of molecular modes and determine the thermal activation energies of the two largest modes, both of which are consistent with the range of literature values for rotational diffusion. We also determine the activation energies for translational diffusion and low-field electron-spin relaxation, both of which are consistent with the literature. Furthermore, we validate the MD simulations at human body temperature and concentrations of the paramagnetic ion used in clinical MRI, and we quantify the uncertainties in both simulations and measurements.

Predicting 1H NMR relaxation in Gd3+-aqua using molecular dynamics simulations.

Atomistic molecular dynamics simulations are used to predict 1H NMR T1 relaxation of water from paramagnetic Gd3+ ions in solution at 25 °C. Simulations of the T1 relaxivity dispersion function r1 computed from the Gd3+-1H dipole-dipole autocorrelation function agree within ≃8% of measurements in the range f0 ≃ 5 ↔ 500 MHz, without any adjustable parameters in the interpretation of the simulations, and without any relaxation models. The simulation results are discussed in the context of the Solomon-Bloembergen-Morgan inner-sphere relaxation model, and the Hwang-Freed outer-sphere relaxation model. Below f0 ≲ 5 MHz, the simulation overestimates r1 compared to measurements, which is used to estimate the zero-field electron-spin relaxation time. The simulations show potential for predicting r1 at high frequencies in chelated Gd3+ contrast-agents used for clinical MRI.