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

Thiago J. Pinheiro dos Santos,Betul Orcan-Ekmekci,Walter G. Chapman,Philip M. Singer,Dilipkumar N. Asthagiri

2024-01-06

Abstract:Traditional theories of the NMR autocorrelation function for intramolecular dipole pairs assume single-exponential decay, yet the calculated autocorrelation of realistic systems display 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 autocorrelation using simple, physical models within a rigorous statistical mechanical development that encompasses both rotational 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 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.

Chemical Physics,Soft Condensed Matter

What problem does this paper attempt to address?

This paper mainly discusses the theoretical and modeling problems of molecular patterns in nuclear magnetic resonance (NMR) relaxation. The traditional NMR autocorrelation function theory assumes single exponential decay, but actual system calculations show multi-exponential behavior, leading to anomalous NMR relaxation dispersion. The researchers proposed a method to simulate and explain the multi-exponential autocorrelation using a simple physical model within the framework of statistical mechanics, while considering rotational and translational diffusion. They transformed the problem into the average of the diffusion propagator, which is described by the Fokker-Planck equation. By expanding the characteristic function of the propagator, each mode has spatial and temporal components related to specific correlation times, similar to the stationary states of a quantum harmonic oscillator.

The researchers validated the statistical mechanics description using molecular dynamics simulations, without using any relaxation models or adjustable parameters to explain the simulation results. They found that the main poles of the simulated autocorrelation function were consistent with the eigenvalues predicted by the theory. The paper also discusses the limitations of early theories, such as the Bloembergen-Purcell-Pound model, and later attempts to explain multi-exponential behavior.

In this way, the paper establishes a theoretical framework for extending the autocorrelation function to multiple "modes," where each mode corresponds to specific molecular motions, aiding in the understanding and prediction of NMR dipole-dipole relaxation dispersion in fluids. For non-interacting dipoles at fixed distances, the traditional BPP theory predicts single exponential decay. However, for dipoles with variable distances, they observed multi-exponential decay, which is related to the different conformations and dynamics of the molecules.

In summary, the paper addresses how to better understand and simulate the multi-exponential behavior in NMR relaxation processes in complex liquids using statistical mechanics methods, providing insights into molecular-level phenomena.