NMR spin-rotation relaxation and diffusion of methane

NMR Relaxation of Gas Adsorbed in Microporous Material

NMR relaxometry has been widely applied to characterize fluid confined in porous media because of its versatility, chemical selectivity, and noninvasive nature. Here we extend its usage to gas adsorbed in microporous materials by establishing a new quantitative model based on the molecular level NMR relaxation mechanism revealed by the molecular simulation of a prototypical adsorption system, CH4 adsorbed in ZIF-8. The model enables new NMR relaxometry-based characterization methods for thermodynamic, dynamic, and structural properties of adsorption systems, as demonstrated and validated by the experiments where the adsorption capacity and self-diffusivity of H2, CH4, and small alcohols adsorbed in ZIF-8 are deduced from the NMR relaxation data. The findings can serve for a better understanding of the composition-structure-properties relationships of a wide range of adsorption systems which is essential for the development and application of new functional microporous materials.

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.

MD Simulation of Water Using a Rigid Body Description Requires a Small Time Step to Ensure Equipartition

In simulations of aqueous systems, it is common to freeze the bond vibration and angle bending modes in water to allow for a longer time step δt for integrating the equations of motion. Thus, δt = 2 fs is often used in simulating rigid models of water. We simulate the SPC/E model of water using δt from 0.5 to 3.0 fs and up to 4 fs using hydrogen mass repartitioning. In these simulations, we find that for all but δt = 0.5 fs, equipartition is not obtained between translational and rotational modes, with the rotational modes exhibiting a lower temperature than the translation modes. To probe the reasons for the lack of equipartition, we study the autocorrelation of the translational velocity of the center of mass and the angular velocity of the rigid water molecule, respectively. We find that the rotational relaxation occurs on a timescale comparable to vibrational periods, calling into question the original motivations for freezing the vibrations. Furthermore, a time step with δt ≥ 1 fs is not able to capture accurately the fast rotational relaxation, which reveals its impact as an effective slowing-down of rotational relaxation. The fluctuation-dissipation relation then leads to the conclusion that the rotational temperature should be cooler for δt greater than the reference value of 0.5 fs. Consideration of fluctuation-dissipation in equilibrium molecular dynamics simulations also emphasizes the need to capture the temporal evolution of fluctuations with fidelity and the role of δt in this regard. The time step also influences the solution thermodynamic properties: both the mean system potential energies and the excess entropy of hydration of a soft repulsive cavity are sensitive to δt.

Kinetics of Guest-Induced Structural Transitions in Metal-Organic-Framework MIL-53(Al)-NH2 Probed by High-Pressure Nuclear Magnetic Resonance

A nuclear magnetic resonance (NMR) study of a pore opening in amino-functionalized metal-organic framework (MOF) MIL-53(Al) in response to methane pressure variation is presented. Variations of both NMR signal intensities and transversal relaxation rates for methane are found to reveal hysteretic structural transitions in the MOF material, which are smeared out over broad pressure ranges. Experiments with pressure reversals upon an incomplete adsorption/desorption gave deeper insight into the microscopic transition mechanisms. These experiments have unequivocally proven that the non-stepwise pore opening/closing transitions observed in the experiments are governed by a distribution of the opening/closing pressures over different MOF crystallites, for example, due to a distribution of the crystal sizes or shapes. The slow kinetics of the structural transitions measured in the hysteresis regime revealed a complex free energy landscape for the phase transition process.

Unified Description of Proton NMR Relaxation in Water, Acetonitrile, and Methane from Molecular Dynamics Simulations in the Liquid, Supercritical, and Gas Phases

A comprehensive calculation of proton NMR relaxation in water, acetonitrile, and methane across a wide range of the phase diagram is provided via ab initio and force-field-based molecular dynamics simulations. The formalism used for the spin-rotation (SR) contribution to relaxation is developed for use with any molecular symmetry and utilizes the full molecular SR tensors, which are calculated from first-principles via Kohn-Sham (KS) DFT. In combination with calculations of the dipolar contribution, near-quantitative agreement with total measured relaxation rates is achieved.

Effect of Nanoconfinement on NMR Relaxation of Heptane in Kerogen from Molecular Simulations and Measurements

Kerogen-rich shale reservoirs will play a key role during the energy transition, yet the effects of nanoconfinement on the NMR relaxation of hydrocarbons in kerogen are poorly understood. We use atomistic MD simulations to investigate the effects of nanoconfinement on the ¹H NMR relaxation times T₁ and T₂ of heptane in kerogen. In the case of T₁, we discover the important role of confinement in reducing T₁ by ∼3 orders of magnitude from that of bulk heptane, in agreement with measurements of heptane dissolved in kerogen from the Kimmeridge Shale, without any models or free parameters. In the case of T₂, we discover that confinement breaks spatial isotropy and gives rise to residual dipolar coupling which reduces T₂ by ∼5 orders of magnitude from the value for bulk heptane. We use the simulated T₂ to calibrate the surface relaxivity and thence predict the pore-size distribution of the organic nanopores in kerogen, without additional experimental data.

Surface Interactions and Nanoconfinement of Methane and Methane plus CO2 Revealed by High-Pressure Magic Angle Spinning NMR Spectroscopy and Molecular Dynamics

This study explores the fundamental, molecular- to microscopic-level behavior of methane gas confined into nanoporous silica proxies with different pore diameters and surface-to-volume (S/V) ratios. Surfaces and pore walls of nanoporous silica matrices are decorated with hydroxyl (-OH) groups, resembling natural heterogeneity. High-pressure MAS NMR was utilized to characterize the interactions between methane and the engineered nanoporous silica proxies under various temperature and pressure regimes. There was a change in the chemical shift position of confined methane slightly in the mixtures with nanoporous silica up to 393 K, as shown by high-pressure ¹³C-NMR. The ¹³C-NMR chemical shift of methane was changed by pressure, explained by the densification of methane inside the nanoporous silica materials. The influence of pore diameter and S/V of the nanoporous silica materials on the behaviors and dynamics of methane were studied. The presence of CO₂ in mixtures of silica and methane needs analysis with caution because CO₂ in a supercritical state and gaseous CO₂ change the original structure of nanoporous silica and change surface area and pore volume. According to simulation, the picosecond scale dynamics of methane confined in larger pores of amorphous silica is faster. In the 4 nm pore, the diffusivity obtained from MD simulations in the pore with a higher S/V ratio is slower due to the trapping of methane molecules in adsorbed layers close to the corrugated pore surface. In contrast, relaxation measured with NMR for smaller pores (higher S/V) exhibits larger T₁, indicating slower relaxation.

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

Atomistic molecular dynamics simulations are used to predict ¹H NMR T₁ relaxation of water from paramagnetic Gd³⁺ ions in solution at 25 °C. Simulations of the T₁ relaxivity dispersion function r₁ computed from the Gd³⁺-¹H dipole-dipole autocorrelation function agree within ≃8% of measurements in the range f₀ ≃ 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 f₀ ≲ 5 MHz, the simulation overestimates r₁ compared to measurements, which is used to estimate the zero-field electron-spin relaxation time. The simulations show potential for predicting r₁ at high frequencies in chelated Gd³⁺ contrast-agents used for clinical MRI.

NMR 1H-1H Dipole Relaxation in Fluids: Relaxation of Individual 1H-1H Pairs versus Relaxation of Molecular Modes

The intramolecular ¹H NMR dipole-dipole relaxation of molecular fluids has traditionally been interpreted within the Bloembergen-Purcell-Pound (BPP) theory of NMR intramolecular relaxation. The BPP theory draws upon Debye's theory for describing the rotational diffusion of the ¹H-¹H pair and predicts a monoexponential decay of the ¹H-¹H dipole-dipole autocorrelation function between distinct spin pairs. Using molecular dynamics (MD) simulations, we show that for both n-heptane and water this is not the case. In particular, the autocorrelation function of individual ¹H-¹H intramolecular pairs itself evinces a rich stretched-exponential behavior, implying a distribution in rotational correlation times. However, for the high-symmetry molecule neopentane, the individual ¹H-¹H intramolecular pairs do conform to the BPP description, suggesting an important role of molecular symmetry in aiding agreement with the BPP model. The intermolecular autocorrelation functions for n-heptane, water, and neopentane also do not admit a monoexponential behavior of individual ¹H-¹H intermolecular pairs at distinct initial separations. We suggest expanding the autocorrelation function in terms of modes, provisionally termed molecular modes, that do have an exponential relaxation behavior. With care, the resulting Fredholm integral equation of the first kind can be inverted to recover the probability distribution of the molecular modes. The advantages and limitations of this approach are noted.

Critical Role of Confinement in the NMR Surface Relaxation and Diffusion of n-Heptane in a Polymer Matrix Revealed by MD Simulations

The mechanism behind the NMR surface-relaxation times (T_1S,2S) and the large T_1S/T_2S ratio of light hydrocarbons confined in the nanopores of kerogen remains poorly understood and consequently has engendered much debate. Toward bringing a molecular-scale resolution to this problem, we present molecular dynamics (MD) simulations of ¹H NMR relaxation and diffusion of n-heptane in a polymer matrix. The high-viscosity polymer is a model for kerogen and bitumen that provides an organic "surface" for heptane. Diffusion of n-heptane shows a power-law dependence on the concentration of n-heptane (ϕ_C7) in the polymer matrix, consistent with Archie's model of tortuosity. We calculate the autocorrelation function G(t) for ¹H-¹H dipole-dipole interactions of n-heptane in the polymer matrix and use this to generate the NMR frequency (f₀) dependence of T_1S,2S as a function of ϕ_C7. We find that increasing molecular confinement increases the correlation time, which decreases the surface-relaxation times for n-heptane in the polymer matrix. For weak confinement (ϕ_C7 > 50 vol %), we find that T_1S/T_2S ≃ 1. Under strong confinement (ϕ_C7 ≲ 50 vol %), we find that T_1S/T_2S ≳ 4 increases with decreasing ϕ_C7 and that the dispersion relation T_1S ∝ f₀ is consistent with previously reported measurements of polydisperse polymers and bitumen. Such frequency dependence in bitumen has been previously attributed to paramagnetism; instead, our studies suggests that ¹H-¹H dipole-dipole interactions enhanced by organic nanopore confinement dominate the NMR response in saturated organic-rich shales.