Nuclear-magnetic-resonance relaxation due to the translational diffusion of fluid confined to quasi-two-dimensional pores.

Nuclear-magnetic-resonance (NMR) relaxation experimentation is an effective technique for nondestructively probing the dynamics of proton-bearing fluids in porous media. The frequency-dependent relaxation rate T_{1}^{-1} can yield a wealth of information on the fluid dynamics within the pore provided data can be fit to a suitable spin diffusion model. A spin diffusion model yields the dipolar correlation function G(t) describing the relative translational motion of pairs of ^{1}H spins which then can be Fourier transformed to yield T_{1}^{-1}. G(t) for spins confined to a quasi-two-dimensional (Q2D) pore of thickness h is determined using theoretical and Monte Carlo techniques. G(t) shows a transition from three- to two-dimensional motion with the transition time proportional to h^{2}. T_{1}^{-1} is found to be independent of frequency over the range 0.01-100 MHz provided h≳5 nm and increases with decreasing frequency and decreasing h for pores of thickness h<3 nm. T_{1}^{-1} increases linearly with the bulk water diffusion correlation time τ_{b} allowing a simple and direct estimate of the bulk water diffusion coefficient from the high-frequency limit of T_{1}^{-1} dispersion measurements in systems where the influence of paramagnetic impurities is negligible. Monte Carlo simulations of hydrated Q2D pores are executed for a range of surface-to-bulk desorption rates for a thin pore. G(t) is found to decorrelate when spins move from the surface to the bulk, display three-dimensional properties at intermediate times, and finally show a bulk-mediated surface diffusion (Lévy) mechanism at longer times. The results may be used to interpret NMR relaxation rates in hydrated porous systems in which the paramagnetic impurity density is negligible.

Model for the interpretation of nuclear magnetic resonance relaxometry of hydrated porous silicate materials.

Nuclear magnetic resonance (NMR) relaxation experimentation is an effective technique for probing the dynamics of proton spins in porous media, but interpretation requires the application of appropriate spin-diffusion models. Molecular dynamics (MD) simulations of porous silicate-based systems containing a quasi-two-dimensional water-filled pore are presented. The MD simulations suggest that the residency time of the water on the pore surface is in the range 0.03-12 ns, typically 2-5 orders of magnitude less than values determined from fits to experimental NMR measurements using the established surface-layer (SL) diffusion models of Korb and co-workers [Phys. Rev. E 56, 1934 (1997)]. Instead, MD identifies four distinct water layers in a tobermorite-based pore containing surface Ca2+ ions. Three highly structured water layers exist within 1 nm of the surface and the central region of the pore contains a homogeneous region of bulklike water. These regions are referred to as layer 1 and 2 (L1, L2), transition layer (TL), and bulk (B), respectively. Guided by the MD simulations, a two-layer (2L) spin-diffusion NMR relaxation model is proposed comprising two two-dimensional layers of slow- and fast-moving water associated with L2 and layers TL+B, respectively. The 2L model provides an improved fit to NMR relaxation times obtained from cementitious material compared to the SL model, yields diffusion correlation times in the range 18-75 ns and 28-40 ps in good agreement with MD, and resolves the surface residency time discrepancy. The 2L model, coupled with NMR relaxation experimentation, provides a simple yet powerful method of characterizing the dynamical properties of proton-bearing porous silicate-based systems such as porous glasses, cementitious materials, and oil-bearing rocks.

Nuclear magnetic resonance relaxometry of water in two and quasi-two dimensions.

Molecular dynamics (MD) and Monte Carlo (MC) methods are used to determine the spin-pair correlation function G(*)(t) for the diffusion of bulk water in three dimensions (3D) and pore water in two dimensions (2D) and quasi-two dimensions (Q2D). The correlation function is required for the determination of the nuclear magnetic resonance spin-lattice and spin-spin relaxation times T(1) and T(2). It is shown that the analytic form of the powder-average correlation function, introduced by Sholl [Sholl, J. Phys. C: Solid State Phys. 7, 3378 (1974)] for the diffusion of spins on a 3D lattice, is of general validity. An analytic expression for G(*)(t) for a uniform spin fluid is derived in 2D. An analytic expression for the long-time behavior of G(*)(t) is derived for spins diffusing on 3D, 2D, and Q2D lattices. An analytic correction term, which accounts for spin pairs outside the scope of the numerical simulations, is derived for 3D and 2D and shown to improve the accuracy of the simulations. The contributions to T(1) due to translational and rotational motion obtained from the MD simulation of bulk water at 300 K are 7.4 s and 10±1 s, respectively, at 150 MHz, leading to an overall time of 4.3±0.4 s compared to the experimental value of 3.8 s. In Q2D systems, in which water is confined by alpha-quartz surfaces to thicknesses of 1-5 nm, T(1) for both translational and rotational relaxation is reduced due to the orientation and adsorption of spins at the surfaces. A method of parametrizing the MC lattice-diffusion simulations in 3D, 2D, and Q2D systems is presented. MC results for G(*)(t) for 3D and 2D systems are found to be consistent with an analytic uniform fluid model for t~/>40 ps. The value of TT(1) for translational diffusion obtained from the MC simulation of bulk water is found to be 4.8 s at 15 MHz. G(*)(t) obtained from MC simulations of Q2D systems, where water is confined by hard walls, is found to execute a distinct transition from 3D to 2D behavior. The T(1) is found to be similar to the 3D bulk water result at all pore thicknesses.