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.

Explicit calculation of nuclear-magnetic-resonance relaxation rates in small pores to elucidate molecular-scale fluid dynamics.

Nuclear-magnetic-resonance (NMR) spin-lattice (T_{1}^{-1}) and spin-spin (T_{2}^{-1}) relaxation rate measurements can act as effective nondestructive probes of the nanoscale dynamics of ^{1}H spins in porous media. In particular, fast-field-cycling T_{1}^{-1} dispersion measurements contain information on the dynamics of diffusing spins over time scales spanning many orders of magnitude. Previously published experimental T_{1}^{-1} dispersions from a plaster paste, synthetic saponite, mortar, and oil-bearing shale are reanalyzed using a model and associated theory which describe the relaxation rate contributions due to the interaction between spin ensembles in quasi-two-dimensional pores. Application of the model yields physically meaningful diffusion correlation times for all systems. In particular, the surface diffusion correlation time and the surface desorption time take similar values for each system, suggesting that surface mobility and desorption are linked processes. The bulk fluid diffusion correlation time is found to be two to five times the value for the pure liquid at room temperature for each system. Reanalysis of the oil-bearing shale yields diffusion time constants for both the oil and water constituents. The shale is found to be oil wetting and the water T_{1}^{-1} dispersion is found to be associated with aqueous Mn^{2+} paramagnetic impurities in the bulk water. These results escalate the NMR T_{1}^{-1} dispersion measurement technique as the primary probe of molecular-scale dynamics in porous media yielding diffusion parameters and a wealth of information on pore morphology.

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.

Solvation pressure in chloroform.

Molecular dynamics (MD) simulations of chloroform vapor and liquid at normal temperature and pressure and liquid under hydrostatic pressure are presented, giving bond lengths and vibrational frequencies as functions of pressure. The change in bond lengths between vapor and liquid at normal temperature and pressure is consistent with a pressure equivalent to the cohesive energy density (CED) of the liquid, supporting the solvation pressure model which predicts that solvated molecules or nanoparticles experience a pressure equal to the CED of the liquid. Experimental data for certain Raman frequencies of chloroform in the vapor phase, in the liquid, and in the liquid under pressure are presented and compared to MD. Results for C-Cl vibrational modes are in general agreement with the solvation pressure model whereas frequencies associated with the C-H bond are not. The results demonstrate that masking interactions exist in the real liquid that can be reduced or eliminated in simplified simulations.