Next Article in Journal
Role of ROS/RNS in Preeclampsia: Are Connexins the Missing Piece?
Next Article in Special Issue
Polymorphic Forms of Valinomycin Investigated by NMR Crystallography
Previous Article in Journal
Conjugates of Ciprofloxacin and Levofloxacin with Cell-Penetrating Peptide Exhibit Antifungal Activity and Mammalian Cytotoxicity
Previous Article in Special Issue
Monitoring the Site-Specific Solid-State NMR Data in Oligopeptides
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Review

A Multi-Scale Study of Water Dynamics under Confinement, Exploiting Numerical Simulations in Relation to NMR Relaxometry, PGSE and NMR Micro-Imaging Experiments: An Application to the Clay/Water Interface

Interfaces, Confinement, Matériaux et Nanostructures (ICMN), UMR 7374, CNRS and Université d’Orléans, 1b rue de la Férollerie, CS 40059, CEDEX 2, F-45071 Orléans, France
*
Authors to whom correspondence should be addressed.
Int. J. Mol. Sci. 2020, 21(13), 4697; https://doi.org/10.3390/ijms21134697
Submission received: 27 May 2020 / Revised: 19 June 2020 / Accepted: 22 June 2020 / Published: 30 June 2020

Abstract

:
Water mobility within the porous network of dense clay sediments was investigated over a broad dynamical range by using 2H nuclear magnetic resonance spectroscopy. Multi-quanta 2H NMR spectroscopy and relaxation measurements were first performed to identify the contributions of the various relaxation mechanisms monitoring the time evolution of the nuclear magnetisation of the confined heavy water. Secondly, multi-quanta spin-locking NMR relaxation measurements were then performed over a broad frequency domain, probing the mobility of the confined water molecules on a time-scale varying between microseconds and milliseconds. Thirdly, 1H NMR pulsed-gradient spin-echo attenuation experiments were performed to quantify water mobility on a time-scale limited by the NMR transverse relaxation time of the confined NMR probe, typically a few milliseconds. Fourthly, the long living quantum state of the magnetisation of quadrupolar nuclei was exploited to probe a two-time correlation function at a time-scale reaching one second. Finally, magnetic resonance imaging measurements allow probing the same dynamical process on time-scales varying between seconds and several hours. In that context, multi-scale modelling is required to interpret these NMR measurements and extract information on the influences of the structural properties of the porous network on the apparent mobility of the diffusing water molecules. That dual experimental and numerical approach appears generalizable to a large variety of porous networks, including zeolites, micelles and synthetic or biological membranes.

1. Introduction

In the last few decades, confined fluids were the subject of numerous experimental [1], theoretical [2] and numerical [3] studies because confinement strongly modifies the physico-chemical properties of bulk liquids. In that framework, the clay/water interface was frequently investigated [4,5,6,7,8,9,10,11,12] because clays exhibit a large variety of physico-chemical properties, such as water absorption, swelling, gelling, thixotropy, ionic exchange capacity and surface acidity exploited in numerous industrial applications, including waste storing [11], cosmetic and paint industries, drilling and heterogeneous catalysis [12]. Furthermore, clays are synthetic or natural materials, easily purified, with a well-controlled structure and atomic composition, leading to ideal models of solid/liquid interfaces. Optimising these numerous applications requires a detailed analysis of the influence of the clay composition on the structural and dynamical properties of the confined fluids.
In that context, a large number of experimental studies, including inelastic neutron scattering (INS) [13,14,15], quasi-elastic neutron scattering (QENS) [16,17], neutron spin echo [18], multi-quanta nuclear magnetic resonance (NMR) relaxometry (i.e., the frequency variation of the NMR relaxation rates) [19,20,21,22,23,24], NMR pulsed-gradient spin-echo (PGSE) attenuation [25,26,27,28,29,30] and NMR micro-imaging [29,30], were performed to quantify the mobility of confined water molecules over a broad range of diffusing time. The short times of mobility for the confined water molecules (from pico-seconds up to 100 nanoseconds) were easily probed by classical neutron scattering experiments (INS, QENS and neutron spin echo). Field-cycling NMR relaxometry is also a powerful tool to investigate the mobility of a large class of confined fluids [31,32,33,34,35,36,37,38,39,40] by probing the frequency variation of the relaxation rate. However, that method becomes useless when the NMR relaxation times of the confined fluids become smaller than the time required to switch the strength of the magnetic field. In that context, multi-quanta NMR relaxometry was successfully applied to probe water mobility in the intermediate time-scale (between microseconds and milliseconds). That time-scale may be extended up to seconds by measuring two-time correlation functions [22,23,24,41] of the long-living quantum state implied in the nuclear magnetisation of confined quadrupolar nuclei. Finally, NMR PGSE attenuation and NMR micro-imaging were also used to investigate the long-time water mobility (between seconds and hours). Because of the large time-scale covered by these numerous experiments (see Figure 1), multi-scale modelling is required to interpret and analyse these various experimental data by using quantum molecular dynamics (for INS) [15], classical molecular dynamics (for QENS and neutron spin echo) [17,18] or Brownian dynamics and differential equations [19,20,21,22,23,24,25,26,27,28,29,30,41] (for NMR relaxometry, two-time correlation function, PGSE and micro-imaging). The same approaches were applied to investigate the mobility of various cations (7Li [42], 23Na [43] and 133Cs [44]) neutralising the negative charge of the clay platelets. While this study focuses on the clay water interfaces, it may be successfully applied to other interfacial systems, including cementous material, zeolite, micro-porous silica [45,46,47,48], synthetic and natural macromolecules.

2. Materials and Methods

Three different classes of clay materials were used in that study: aqueous diluted dispersions [49], dense sediments [23] with a well-controlled hydration state and compacted aggregates [29]. While swelling clays (montmorillonie [50], beidellite [23]) are implied in the two first studies, the samples of the third study are prepared with kaolinite, i.e., unswelling clay [29]. The swelling clays result from the sandwiching of one octahedral layer of aluminum or magnesium oxides between two layers of tetrahedral silicon oxides. Because of some substitutions of these metallic atoms by less charged atoms, the individual clay platelet bears negative electric charge neutralised by solvated and exchangeable counterions. The charge locations and the nature of these counterions play central roles in the swelling behaviour of these clay platelets [50,51,52].
Montmorillonite was used to study dilute aqueous clay dispersions. It is a natural clay purchased from the Source Clays Mineral repository at Purdue University. It was prepared and purified according to classical procedure [50]. A sequence of centrifugation was also performed to reduce the size dispersion of the clay which is quantified by analysing TEM micrographs [49]. Two clay samples were selected for that study, with average lateral extents of ( 70 ± 10 ) nm and ( 300 ± 50 ) nm for the samples labelled Size3 and Size2 respectively. The general formula per unit cell of the purified montmorillonite is: (Al 3.06 Mg 0.48 Fe 0.46 )(Si 7.76 Al 0.24 )O 20 (OH) 4 Na 0.77 [49].
Beidellite is natural clay used to prepare self-supporting clay films with well characterised degrees of hydration. It is purchased from the Source Clays Mineral repository at Purdue University and purified according the same classical procedures [50]. After purification, its general formula per unit cell is: (Al 3.77 Mg 0.21 Fe 0.11 )(Si 7.27 Al 0.73 )O 20 (OH) 4 Na 0.67 [23]. Clay films are obtained by oedometric compression of diluted aqueous dispersions and dried under nitrogen flux before equilibration with a reservoir of heavy water (from Aldrich, purity 99.9%) at fixed water chemical potential by using saturated salt solutions [23]. We fix the water partial pressure p/ p 0 = 0.92 by using saturated KNO 3 aqueous solution of pure heavy water. This partial pressure was selected because it corresponds to a well characterised hydration state with an interlayer separation of 15.6 Å [53,54,55] corresponding to two hydration layers, as illustrated in Figure 2b. In order to avoid clay drying during the set of NMR experiments, a small reservoir of heavy water at fixed chemical potential is inserted in the cap of the NMR tube, out of the detection coil.
Kaolinite is the non-swelling clay selected for the last study. It results from the stacking of platelets composed from the superposition of one layer of tetrahedral silicon with one layer of octahedral aluminium oxides. The cohesion between individual platelets is insured by a dense network of hydrogen bonds between the hydroxides of octahedral aluminium oxides of one platelet and the oxygen atoms from the tetrahedral layer of silicon oxides pertaining to the neighbouring platelet. The kaolinite sample was purchased from the Source Clays Mineral repository at Purdue University. Its general formula per unit cell is: (Al 1.9 Ti 0.065 Fe 0.035 )(Si 1.92 Al 0.08 )O 5 (OH) 4 [29,30]. Compacted clay samples are prepared by uniaxial compaction of the dry clay powder [29] and its porosity is determined by helium pycnometry [29].
2H NMR spectra and relaxation measurements were performed on a DSX360 Bruker spectrometer operating at a field of 8.465 Tesla. On this spectrometer, typical pulse duration for the total inversion of the magnetisation is equal to 23 microseconds. Spectra were recorded by using a fast acquisition mode with a typical sampling time of 0.25 microseconds, corresponding to a spectral width of 4 MHz. The NMR spectra are classically recorded by using a single detection pulse. Rectangular clay lamellae (typically 30 × 5.5 mm 2 ) are cut in the self-supporting clay films and are inserted in a sealed cylindrical tube (see Figure 3a) fitting the gap of the home-made solenoidal coil used to perform the various NMR measurements [20,21,22,23,24,27,28,29,30,41,42,43,44]. As a consequence, these NMR measurements can be performed at different orientations of the clay film by reference with the static magnetic field B 0 (see Figure 3a). In the case of aqueous clay dispersions, this setup was slightly modified: the sample holder was a tube fitting the diameter of the coil, and its axis was oriented perpendicularly to the coil axis (see Figure 3b). It can be rotated by 90 ° while maintaining the sample in the magnet, allowing one to perform dynamical studies of the clay mobility into the static magnetic field B 0 .
1H NMR PGSE self-diffusion and magnetic resonance imaging (MRI) measurements are performed on a Bruker DSX100 spectrometer operating at a static magnetic field B 0 of 2.351 Tesla. This spectrometer is equipped with a saddle coil, with the axis of the sample holder parallel to the static magnetic field B 0 . Typical pulse duration for a complete inversion of the proton magnetisation is 10 microseconds.

3. Results and Discussion

3.1. NMR Spectra

In contrast with bulk water, splitting of the 2H NMR resonance line of heavy water (see Figure 2a) is a general trend of the confined water molecules in the presence of clay platelets [56,57,58]. It results from a partial average of the instantaneous quadrupolar coupling (see Equation (A5) in the Appendix A) felt by the deuterium atom. Since its principal axis is oriented along the O D director [59], any asymmetry of the water reorientation motions induces such residual quadrupolar coupling. The NMR signal of bulk water in the outer reservoir (see Section 2) induces the distortion of the NMR spectra (see Figure 2a). Therefore, the data displayed in Figure 2d were compared [23] with the residual quadrupolar coupling, noted ν Q , extracted from the oscillations of the transverse magnetisation detected by using a Hahn echo pulse sequence (see Figure 4a). As displayed in Figure 2d, the agreement between both set of measurements is quantitative. As predicted by molecular modelling of the clay/water interface [54,55,60], water molecules exhibit specific orientation in the vicinity of the clay surface. As illustrated in Figure 2b, the configurations of the water molecules obtained by grand canonical Monte Carlo simulations confirm the specific orientation of confined water molecules in contact with the clay surface (see Figure 6c in the reference [23]) and are responsible for the splitting of the 2H NMR spectra. As a consequence, the observed residual splitting results from an average over the distribution of the clay directors (see Equation (A6) in the Appendix A) inside the clay film in addition to the specific organisation of the water molecules confined in the interlamellar space of the clay platelets (see Figure 2c), leading to an apparent residual splitting varying according to the film orientation inside the static magnetic field B 0 (see Figure 2d).
Such residual quadrupolar coupling was also detected (Figure 5a,b) in the presence of dilute and isotropic aqueous clay dispersions [49,56,57,58] resulting from the individual platelet orientation induced by the static magnetic field B 0 . Such a platelet alignment in the static magnetic field results from the large anisotropy of their magnetic susceptibility [49]. As previously illustrated [25], the water self-diffusion tensor may be used to identify the preferential orientation of the clay directors. As shown in Figure 5c, the clay directors of dilute aqueous dispersions of montmorillonite are preferentially oriented perpendicular to the static magnetic field. Figure 5a also exhibits reduced asymmetry of the water doublet. It results from the paramagnetic impurities (iron) of the clay platelets. As detailed in the Appendix, the resulting paramagnetic coupling interferes with the quadrupolar coupling responsible for the deuterium relaxation and induces such asymmetry of the deuterium doublet [57,61]. As displayed in Figure 5a, the platelet orientation is not instantaneous and depends upon the clay concentration and the strength of the static magnetic field (see Figure 5b). The strong enhancement of the time-scale characterising the clay reorientation in the static magnetic field reported in Figure 5d is the fingerprint of their mutual repulsion resulting from the overlap of their ionic diffuse layers of neutralising counterions [1,49,61]. It corresponds to a microscopic detection of the sol–gel transition occurring at slightly lower clay concentrations than those detected by macroscopic rheological measurements [49].

3.2. NMR Relaxometry

The second source of difference between isotropic and confined liquids is given by the orders of magnitude of the longitudinal ( R 1 = 1 / T 1 ) and transverse ( R 2 = 1 / T 2 ) NMR relaxation rates: while both relaxation rates are general equal in bulk liquids at room temperature, confinement induces a large enhancement of the NMR transverse relaxation rate by reference to its longitudinal counterpart (see Appendix). In the framework of the classical BPP model [61,62], this behaviour was previously assigned [57,63] to a large increase of the time-scale characterising the decorrelation of the relaxation mechanism monitoring the NMR relaxation process of the confined molecules by reference to their dynamical behaviour in bulk liquid. While the classical BPP model adequately interprets the NMR relaxation process in bulk liquids, it becomes useless in analysing NMR relaxation of the same confined fluid. As an example, the 2H NMR relaxation rate of bulk heavy water ( R 2 = R 1 = 2.4 s 1 ) is compatible with an orientational correlation time of a few pico-seconds, in agreement with QENS measurements. By contrast, the large enhancement of 2H transverse relaxation time of heavy water in the presence of swelling clays ( R 2 = 13,000 s 1 and R 1 = 110 s 1 ) [64] requires, in the framework of the BPP model, an orientational correlation time of 20 nanoseconds while its order of magnitude reported by QENS measurements still fits a few pico-seconds! As a consequence, one must definitively reject any interpretation of the NMR relaxation rates of confined fluids on the basis of this classical BPP model. In that context, numerous analytical [36] and numerical [65,66,67] models were developed to interpret such enhancement of the NMR transverse relaxation detected for confined fluids. In contrast with the isotropic BPP model, that interpretation of the NMR relaxation of confined fluids focuses on the diffusion time [36,65,66,67] required by the confined molecular probe to explore different local environments and lose the memory of their initial configuration.
In the case of dilute aqueous dispersions [19,49,57], that diffusion time corresponds to the average residence time τ c of the water molecules at the surface of the clay particle. For dense clay sediments, the same diffusion time corresponds to the average residence time τ c of the water molecules confined between two clay platelets [20,21,22,23,24]. NMR relaxometry measurements are then performed to investigate the frequency variation of the NMR relaxation rates and extract the average residence time τ c of the diffusing probes. The order of magnitude of the residence time τ c is simply evaluated by inverting the characteristic angular velocity ω c . In that context, determining the average size of these clay platelets is a simple way to estimate the order of magnitude of the mobility of these confined water molecules [20,21,22,23,24].
One may be tempted to interpret the large enhancement of the transverse relaxation rate on the basis of iron clusters dispersed in the clay sediments and inducing heterogeneities of the magnetic susceptibility within the mineral. This analysis is not valid in our case since the longitudinal relaxation rate, a priori not affected by heterogeneities of the chemical shift [68]), is also enhanced compared to its value in the bulk. To avoid such heterogeneities of the magnetic susceptibility, we used either synthetic or purified natural clays (see Section 2) in which iron atoms are only present in limited quantity and dispersed within the clay matrix by substituting metallic atoms of the clay network (see Section 2). Furthermore, the paramagnetic contribution of these atoms to the relaxation mechanisms is fully described in our analysis (see Appendix A and Appendix D) and the multi-quanta relaxation measurements were performed to separately quantify the contributions of the quadrupolar and paramagnetic relaxation mechanisms. Nevertheless, we also observed large enhancement of transverse relaxation rate of the confined water molecules in comparison with their longitudinal relaxation rate, resulting from the slow modulation of these relaxation mechanisms. This was previously demonstrated by measuring, at various strengths of the static magnetic field, the relaxation rates of confined quadrupolar nuclei in the presence of synthetic and natural but purified clay sediments [69]. While the transverse relaxation rate appears independent of the static magnetic field B 0 , the longitudinal relaxation rate increases with a decrease of the magnetic field. A completely opposite trend should be detected for relaxation enhancements induced by the heterogeneities of the magnetic field; i.e., no influence of the strength of the static magnetic field on the longitudinal relaxation rates and an increase of the transverse relaxation rate as a function of the increase of the static magnetic field!
The relaxation of the deuterium nuclei is mainly monitored by the quadrupolar coupling. As explained in the Appendix (see Appendix A), this quadrupolar coupling is modulated only by the reorientation of the water molecule (see Equations (A4) and (A5)), more precisely by the fluctuations of the O D director into the static magnetic field B 0 , since the principal axis of the electrostatic coupling felt by the deuterium atom of heavy water is oriented along the O D director [59]. This is the reason why we should expect a short correlation time (i.e., a few ps), in agreement with QENS measurements [17]. This situation differs strongly from the modulation of the dipolar coupling, which includes also a contribution from the separation between the pair in interacting dipolar nuclei; see Equation (A9). As a consequence, molecular diffusion of the confined dipolar probes becomes an important source of NMR relaxation, inducing a slow modulation of the NMR dipolar coupling in the case of confined fluids [33,35,36,39,67,70]. However, because of the specific orientation of the water molecules confined in the interlamellar space between two clay platelets, the orientation fluctuations of the O D director of these confined water molecules is strongly restricted until their desorption. This point is well illustrated by numerical simulations [21]. As a consequence, the residence time of the water molecules adsorbed within in clay porous network monitors the decorrelation of their quadrupolar coupling, inducing the abovementioned tremendous increase of the correlation time.
For that purpose, a set of multi-quanta NMR relaxation measurements are performed with the pulse sequences illustrated in Figure 4a–c. Figure 6a displays the corresponding relaxation rates from which we can extract the apparent spectral densities measured in the laboratory frame (Figure 6b); see Equations (A21) and (A22a)–(A22e) in the Appendix D. In that approach, we focus on the dominant contributions to the NMR relaxation mechanisms. As detailed in the Appendix B in Equations (A10)–(A15), the so-called spectral densities J ( ω ) are used to evaluate the contributions of the various relaxation mechanisms to the NMR relaxation rates. Since under confinement each spectral density (see Equation (A23)) evaluated at zero angular velocity is always larger than its counterpart evaluated at high angular velocity,
J ( 0 ) J ( ω 0 ) J ( 2 ω 0 )
we obtain a limited number of parameters characterising all the multi-quanta NMR relaxation rates (see Equation (A24) in the Appendix D). By using Wigner rotation matrices (cf. Equations (A16a)–(A16c) in the Appendix B), we can easily extract from these apparent spectral densities their intrinsic counterparts evaluated in the frame of the oriented clay film in which the water diffusional process occurs. As shown in the Appendix D in the Equation (A24), a limited number of parameters are required to fully quantify the angular variation of the multi-quanta relaxation rates. That analysis allows not only distinguishing the relative contributions from the quadrupolar and paramagnetic relaxation mechanisms but also determines the intrinsic spectral densities evaluated at zero angular velocity in the frame of the clay film (see Figure 6c,d and Equations (A16a)–(A16c) in the Appendix B). As illustrated in Figure 6c, the apparent spectral density of the quadrupolar coupling J 0 Q ( 0 ) exhibits a large angular variation because its intrinsic contributions significantly differ
J 1 Q , i n t r i n s i c ( 0 ) J 0 Q , i n t r i n s i c ( 0 ) > J 2 Q , i n t r i n s i c ( 0 )
while the angular variation of the apparent spectral density of the paramagnetic coupling ( J 0 D ( 0 ) ) is less noticeable (Figure 6d) because of the limited difference between its intrinsic components
J 0 D , i n t r i n s i c ( 0 ) J 1 D , i n t r i n s i c ( 0 ) > J 2 D , i n t r i n s i c ( 0 )
These relationships are not specific to water molecules confined within beidellite sediments but were also reported for various clay sediments [20,21,22,23,24].
Figure 7 illustrates the angular variation of the relative weight of the various intrinsic contributions to the m = 0 component of the apparent spectral density monitoring the quadrupolar and dipolar relaxation mechanisms (see Equation (A16a) in the Appendix B). By coupling these results with the data displayed in Figure 6c,d, it becomes possible to extract all the required dynamical information by performing a limited number of multi-quanta relaxometry measurements. For that purpose, we selected only three different orientations of the clay film into the static magnetic field; i.e., 0 ° , 90 ° and 30 ° . At 0 ° , only the m = 0 intrinsic component contributes to the apparent relaxation process. Then, measurements performed at 90 ° , still remove the contribution from the m = 1 intrinsic component, but simultaneously optimise the contribution of the m = 2 intrinsic component. Finally, an intermediate angle (30 ° ) is selected to add a reduced contribution of the m = 1 intrinsic component. We avoid using the angle value 45 ° , i.e., the angle optimising this contribution, because of the large enhancement expected to occur from the quadrupolar relaxation mechanism (see Figure 6c).
Figure 8 illustrates the pulse sequence used for these multi-quanta NMR relaxometry measurements. Relaxometry measurements were already performed for various confined fluids [19] by using only NMR locking of the longitudinal and transverse magnetisations of dipolar nuclei (spin I = 1/2). The use of quadrupolar nuclei (spin I > 1/2) allows performing here multi-quanta relaxometry, increasing significantly the range of probed angular velocities. This point is well illustrated In Table 1 by listing the broad range of probed angular velocities with a limited number of irradiation powers corresponding to the values of ω 1 . Figure 9 displays the typical agreement between experimental data and fitted curves. Fourier transforms of these results better illustrate the specific resonance frequencies (see Equation (A20) in the Appendix) probed by these multi-quanta spin-locking relaxation measurements. Finally, Figure 10 displays the variation of the reduced spectral densities probed by these multi-quanta NMR experiments. Two frequency domains are clearly displayed, with a plateau at low frequency and a unique monotonous decrease at high frequency. The transition between these two dynamical regimes occurs at ( 6 ± 1 ) × 10 4 rad/s, leading to a characteristic residence time of ( 17 ± 3 ) microseconds for the water molecules confined within the interlamellar space of beidellite clay, in agreement with the data obtained by QENS measurements [17] and MD numerical simulations [17] performed at a much smaller time-scale. Such agreement was also obtained for various heavy water molecules confined between clay lamellae under controlled hydration conditions [20,21,22,24].

3.3. Two-Time Correlation Function

Dynamical information on the long-time mobility of the confined water molecules may be extracted by measuring two-time correlation function. The pulse sequence used to perform said measurement is displayed in Figure 11: it is based on echo attenuation measurement by saving the magnetisation during the long evolution period ( τ M ) on the long-living T 20 coherence or quantum spin state (cf. Equations (A21)–(A24) in the Appendix). That sequence exploits the heterogeneity of the platelet orientations within the clay sediment (see Figure 2c). During the first evolution period, the first-order coherence T 10 oscillates according to the initial local value of the residual quadrupolar coupling (noted ω Q ( 0 ) ) felt by the confined water molecules. After a mixing period τ M , the same water molecules diffuse through the clay film and may reach a new clay aggregate with a different orientation within the film, corresponding to another residual quadrupolar coupling (noted ω Q ( τ M ) )). As a consequence, the measured magnetisation evolves according to [41,71,72]:
I ( t e , τ M ) = cos ω Q i ( 0 ) t e ) × cos ω Q j ( τ M ) t e ) × e ( R 20 τ M + 2 R 11 t e )
If the mixing time τ M is shorter than the time-scale characterising water exchange between platelets pertaining to different orientations, both residual quadrupolar couplings are the same and Equation (4) predicts damped oscillations of a squared cosinus function. By contrast, a strong attenuation of these oscillations is expected to occur if the mixing time exceeds the residence time of the water molecules confined within their initial aggregate.
Figure 12a results from measurements performed on the same beidellite clay sample. After normalisation of the first point, Figure 12b displays a noticeable decrease of the intensity of the abovementioned oscillations of the initial squared cosinus function resulting from water exchange between aggregates with different orientations within the clay film. The exponent α = 1.5 is heuristic and has no physical meaning. It results from a fit of both experimental and numerical data and is used to quantify the order of magnitude of the time τ e x c h . The time-scale corresponding to said exchange τ e x c h is quite long ( 33 ± 5 ) milliseconds, much longer than the average residence time τ c of the same water molecules in the interlamellar space inside the same aggregate ( 17 ± 3 ) microseconds. By using the self-diffusion coefficient of bulk water D b u l k , one can estimate the upper limit of the correlation length characterising the propagation of the orientation ordering of the clay aggregates within the film
L m a x = 2 D b u l k τ e x c h = 11 ± 2 μ m
that is one order of magnitude larger than the average size of the beidellite platelets used in that study.
The decrease of the two-time correlation function displayed in Figure 12b is modelled by a set of differential equations describing the exchange between 3375 elementary cubic sub-cells within the clay sediment, each characterised by its residual quadrupolar coupling [41]. A Gaussian random field, with a priori field-field correlation [73,74], was used to generate the distribution of the 3375 local values of the residual quadrupolar coupling. The standard deviation of this Gaussian field–field correlation function [41] was selected to match the length of two elementary sub-cells. Each sub-cell is indexed by three indices ( i , j , k ), allowing us to describe water exchange between the magnetisation pertaining to neighbouring sub-cells:
d σ i , j , k d t = R i , j , k 6 k e x c h t σ i , j , k + k e x c h t σ i + 1 , j , k + σ i 1 , j , k + σ i , j + 1 , k + σ i , j 1 , k + σ i , j , k + 1 + σ i , j , k 1
where R i , j , k contains contributions from pulse, local residual quadrupolar coupling and relaxation mechanisms. Equation (6) is a generalised Bloch equation, modified to describe spin-exchange between neighbouring elementary sub-cells [68]. The set of differential equations described by Equation (6) was solved numerically by using an iterative method [41]. As displayed in Figure 12b, the agreement with experimental data is quantitative, validating our interpretation of this two-time correlation function.

3.4. Macroscopic Water Diffusion

Kaolinite compacted sample (dry sample porosity ϕ = 0.52) was selected to investigate water mobility within unswelling aggregates. After saturation of the clay aggregate by pure water, 1 H NMR PGSE attenuation measurements were performed to measure the water diffusion tensor inside the sediment (see Section 2). Figure 13 illustrates the pulse sequence used to perform the PGSE NMR attenuation measurements [29]. The encoding of transverse magnetisation is performed by a pulse sequence composed of an initial π / 2 pulse used to create transverse magnetisation, a set of bipolar field gradients (duration δ , strength G) applied along any desired direction and a final π / 2 pulse to transfer back magnetisation in the direction parallel to the static magnetic field B 0 (called longitudinal direction). After the evolution period ( Δ ), magnetisation is refocused by using a pulse sequence similar to the encoding one, but with opposite field gradients. That sequence allows probing water mobility along any direction.
Because of the presence of iron (see Section 2), the water relaxation times in saturated sediment ( T 1 = 5.3 × 10 2 s and T 2 = 1.2 × 10 3 s) are much shorter than their value in bulk water ( T 1 = T 2 = 3 s). As a consequence, optimisation of the pulse sequence requires limiting the mixing time Δ to 150 ms, i.e., no more than three T 1 values (see Figure 13). Another limitation of this PGSE pulse sequence originates from the time required to build the magnetic field gradients (typically 60 microseconds). As a consequence, the total duration ( 4 τ = 3.04 ms) of the coding and refocusing pulse sequence (see Figure 13) reaches milliseconds, i.e., the order of magnitude of the transverse relaxation time T 2 (see above), leading again to important attenuation of the detectable NMR signal. Finally, the intensity of the NMR signal varies according to [75,76]
E ( G , Δ ) = C exp 4 π 2 q 2 e G T D e G Δ + 3 2 τ δ 6 4 τ T 2 Δ T 1
where q = γ δ G / π ; γ is the 1H gyromagnetic ratio ( γ = 2.6752 × 10 8 rad/s); G is the strength of the field gradient and e G its director; D is the water self-diffusion tensor; and δ , Δ and τ are time delays illustrated in Figure 13. Since all measurements are with a unique set of time delays ( δ = 500 μ s, Δ = 20 ms and τ = 760 μ s), the attenuation of the NMR signal satisfies
E ( G , Δ ) = E ( 0 , Δ ) exp 4 π 2 q 2 e G T D e G Δ + 3 2 τ δ 6
reducing the signal attenuation to the contribution from the self-diffusion propagator, i.e., the intermediate scattering function probed by QENS [77]. The maximum resolution of these PGSE NMR measurements is limited by the maximum strength of the magnetic field gradient ( G m a x = 1.6 Tesla/m) [76]:
resolution m a x = 1 q m a x = π γ δ G m a x 15 μ m
Figure 14 illustrates the results obtained by 1H PGSE NMR attenuation measurements within water saturated kaolinite aggregates. The size of the individual kaolinite platelets (∼100 nm) was determined by TEM [29]. It was two orders of magnitude smaller than the maximum resolution of these PGSE measurements. We did not detect significant variation of the self-diffusion coefficient as a function of diffusion time Δ [29].
As displayed in Figure 14, the attenuation perfectly fit a Gaussien relationship, whatever the direction probed by the applied field gradient. Six non-collinear directions of the field gradient were selected and are noted e 1 = ( 1 , 0 , 0 ) , e 2 = ( 0 , 1 , 0 ) , e 3 = ( 0 , 0 , 1 ) , e 4 = ( 1 / 2 , 1 / 2 , 0 ) , e 5 = ( 1 / 2 , 0 , 1 / 2 ) and e 6 = ( 0 , 1 / 2 , 1 / 2 ) , with e 3 oriented parallel to the cylindrical axis of the clay sample.
Said procedure allows determining the principal axes of the self-diffusion tensor and the corresponding eigenvalues, D l o n g = ( 0.92 ± 0.02 ) × 10 9 m 2 /s and D t r a n s = ( 1.05 ± 0.02 ) × 10 9 m 2 /s, exhibiting reduced anisotropy of water mobility within the clay sample. Since the longitudinal direction coincides with the compression axis of the clay sample, we can deduce that the small decrease of the water mobility in that direction is induced by the sample preparation. From the water mobility of bulk water D b u l k = ( 2.10 ± 0.02 ) × 10 9 m 2 /s, we can evaluate the tortuosity factor θ i of the sample resulting from the geometry and tortuosity of its porous network by evaluating the ratio θ i = D b u l k / D i , leading to θ l o n g = 2.28 ± 0.05 and θ t r a n s = 2.00 ± 0.05 .
Nuclear magnetic imaging may be also used to investigate water mobility within the kaolinite porous network in a very long time-scale (up to a few days), better corresponding to macroscopic scales implied in industrial or geological processes. For that purpose, a limited amount of heavy water (0.23 mL) is added on the top of the water saturated kaolinite compacted sample and MRI is used to probe the macroscopic exchange between heavy and natural water by analysing the time evolution of H 2 O concentration profiles. The MRI measurements were performed on the same Bruker DSX100 spectrometer (see Section 2). The pulse sequence used to perform magnetic resonance imaging (MRI) experiments is illustrated in Figure 15. After an initial π / 2 excitation pulse, a magnetic field gradient ( G z ) is applied during a period δ / 2 in the vertical direction, parallel to the static magnetic field, thereby corresponding to the cylindrical axis of the sample holder. Inversion pulse (noted π ) is then applied and the same field gradient ( G z ) is applied during the period δ with a simultaneous acquisition of the NMR signal. This pulse sequence generates an echo after the evolution period δ / 2 leading to a noticeable attenuation ( exp δ / ( 2 T 2 ) 0.34 ) of the NMR signal pertaining to the confined water molecules without altering the NMR signal from bulk water. These MRI water concentration profiles are recorded with 256 elementary time steps of 10 μ s during the evolution period δ , leading to water concentration profiles composed of 256 water layers along the O z direction. In order to simultaneously detect the water molecules confined inside the clay aggregate (length L 1 = 12.3 mm) and free above it (length L 2 = 7 mm), we probe a total spatial window ( L T ) of 25 mm, see Figure 16a,b, leading to a spatial resolution ( L T / 256) of 97.7 μ m. For that purpose, a magnetic field gradient ( G z ) of 0.094 Tesla/m is applied, corresponding to a wave vector q = 1.024 × 10 4 m 1 (see above) leading to the required spatial resolution. In order to improve the signal/noise ratio, the pulse sequence displayed in Figure 15 is repeated 16 times with a time delay of 0.5 s. That repetition delay is much longer than the longitudinal relaxation time of confined water molecules ( T 1 = 5.3 × 10 2 s), thus optimising acquisition of their NMR signal. By contrast, the NMR signal of bulk water ( T 1 = 3 s) is then saturated, leading to non-quantitative detection of the bulk water concentration profiles displayed on the right side of Figure 16a.
Figure 16a displays the water concentration profiles detected after addition of 0.23 mL of bulk water on the top of the kaolinite sample, where the compression was initially applied (see Section 2). As displayed in Figure 16a, significant signal attenuations occur at the lower and upper limits of the spectral window because it exceeds the spatial domain corresponding to complete excitation and detection of the NMR coil. For that purpose, we restrict our analysis by integrating the water concentration profiles over a limited area of 5 mm centred on the clay sample (see Figure 16a). Nevertheless, Figure 16a still exhibits a linear decrease of the initial water concentration profile probably due to compaction implied in the sample preparation (see Section 2), leading to lower porosity at the top of the clay sample, where the compression was applied (see above). Figure 17 exhibits the long-time variation of the partially-integrated water concentration profiles, after normalisation by the initial value. The limiting asymptotic value may be used to evaluate the porosity (noted φ ) of the hydrated sample since
asymptotic value = total H 2 O total H 2 O + D 2 O = L 1 φ L 2 + L 1 φ 0.49 ± 0.01
where L 1 (∼12.3 mm) and L 2 (∼7 mm) are the lengths occupied by the confined hydrated sample and the added bulk water respectively (see above). The corresponding porosity φ = 0.55 ± 0.01 is fully compatible with the value determined by picnometry for the compacted sample of dry kaolinite ( φ = 0.52 ), since no mechanical constraint is applied during hydration of the kaolinite sample. These measurements are fully validated by numerical simulations of Brownian dynamics—modelling the macroscopic exchange between bulk water and added heavy water (Figure 16b). For that purpose, the water self-diffusion coefficient of the confined water molecules was set equal to the value measured by PGSE in the longitudinal direction; i.e., 0.92 × 10 9 m 2 /s (see above). During these BD simulations, the water exchange at the boundary between free and confined water was constricted to insure conservation of the water local density, thereby describing, at a microscopic scale, the incompressibility of liquid water. As displayed in Figure 16a,b and Figure 17, semi-quantitative agreements are obtained between experimental and numerical data, validating the use of MRI concentration profiles to probe such long-time interfacial water exchanges between macroscopic samples of confined and bulk water.

4. Conclusions

Multi-scale analysis of the structural and dynamical properties of water molecules confined by clay porous network was determined by a set of complementary nuclear magnetic resonance experiments exploiting spectroscopy, multi-quanta relaxometry, two-time correlation functions, pulsed gradient spin-echo attenuation and magnetic resonance imaging. This whole set of experimental investigations allows one to cover a broad range of characteristic times, varying between microseconds and days. In that context, a complementary set of numerical models, including molecular dynamics, Brownian dynamics and resolution of macroscopic differential equations, was used to validate the interpretation of these multi-scale NMR measurements performed at such various time-scales. For that purpose, bottom-up approaches are applied to fix the values of parameters required by coarse grain modelling from numerical results obtained by simulating more detailed systems. While this study was performed for various clay/water interfaces, including dilute aqueous dispersions and dense sediments or aggregates, it can be extended to other deuterated molecules or quadrupolar ions (7Li, 23Na, 133Cs) confined by various interfacial systems, including zeolites and porous silicates, micelles and membranes, polyelectrolytes and macromolecules.

Author Contributions

Conceptualization, writing—original draft preparation, writing—review and editing, P.P. and A.D.; funding acquisition, A.D. All authors have read and agreed to the published version of the manuscript.

Funding

The NMR spectrometers used in that study were purchased thanks to grants from Région Centre (France). This study was partially supported by grants from the CNRS Interdisciplinary Project “Nucléaire, Energie, Environnement, Déchets, Société (NEEDS)” through its “Milieux Poreux (MiPor)” Program, which provided financial support (projects MultiDyn and TransReac).

Acknowledgments

The authors cordially thank A.M. Faugère for technical assistance with the NMR spectrometers, Badets for devising the NMR specific probe, F. Condé for pycnometry and F. Warmont for TEM micrographs (ICMN, Orléans, France). We are also grateful to L.J. Michot (PHENIX, Paris, France), and E. Ferrage, F. Hubert and Th. Dabat (IC2MP-Hydrasa, Poitiers, France) for providing purified clay samples. The authors gratefully acknowledge the anonymous reviewers for their constructive comments to improve the manuscript.

Conflicts of Interest

The authors declare no conflict of interest.

Abbreviations

The following abbreviations are used in this manuscript:
EFGElectric Field Gradient
GCMCGrand Canonical Monte Carlo
INSInelastic Neutron scattering
MRIMagnetic Resonance Imaging
PGSEPulsed-Gradient Spin-Echo
QENSQuasi Elastic Neutron Scattering
TEMTransmission Electron Microscopy

Appendix A. Quadrupolar and Heteronuclear Dipolar Hamiltonian

The quadrupolar Hamiltonian [33,78,79] is defined by:
H Q = C Q p = 2 2 ( 1 ) p F 2 , p Q , L T 2 , p I R with C Q = 3 2 e Q ( 1 + γ ) I ( 2 I 1 )
where e is the electron charge, Q is the quadrupolar moment of the nuclei [80] and ( 1 + γ ) is the Steinhermer antishielding factor [80]. In the above equation,
F 2 , 0 Q , L = 1 2 V z z L , F 2 , ± 1 Q , L = 1 6 V x z L ± i V y z L , F 2 , ± 2 Q , L = 1 2 6 V x x L V y y L ± 2 i V x y L
and
T 2 , 0 Q , I R = 1 6 3 I z 2 I ( I + 1 ) , T 2 , ± 1 Q , I R = 1 2 I z I ± + I ± I z , T 2 , ± 2 Q , I R = 1 2 I ± 2
where V α β L are the components of the electric field gradient (EFG) evaluated in the laboratory frame (noted L); T 2 , ± p Q , I R (for p = −2 to 2) are the second order irreducible tensors operators; I x , I y and I z are the spin operators; and I ± = I x ± i I y .
In the presence of a static quadrupolar coupling, the equidistant Zeeman energy levels are modified by the residual quadrupolar coupling, leading to a quadrupolar splitting of the resonance lines according to:
ω m 1 , m = 3 8 C Q V z z L ( 1 2 m )
for m varying between I and I + 1 .
During a change of frame, the components of the EFG transform like the second order spherical harmonics [33]:
F 2 , q Q , L = p = 2 2 F 2 , p Q , P D p , q L P ( θ , φ , ψ )
with D p , q L P ( θ , φ , ψ ) , the components of the Wigner rotation matrices [33] where the set of ( θ , φ , ψ ) Euler angles defines the orientation, into the static magnetic field, of the principal axis of the tensor describing the EFG felt by the quadrupolar nucleus. Three sets of frames are useful to describe the orientation of the principal component of the EFG: the laboratory frame (noted L), a frame attached to the dense clay sediment (noted F) and a frame attached to the individual quadrupolar nucleus (noted P; see Figure 2c)). The e z directors of these different frames are, respectively, the direction of the static magnetic field B 0 (laboratory frame L), the normal to the clay sediment n (sediment frame F) and the director of the principal component of the EFG, noted V z z P (particle frame P; see Figure 2c)).
The measured quadrupolar splitting is derived from Equations (A3) and (A4):
ω m 1 , m a p p = A m V z z P p = 2 2 D p , 0 L F ( θ L F , φ L F , ψ L F ) D 0 , p F P ( θ F P , φ F P , ψ F P ) with A m = 3 e Q ( 1 + γ ) ( 1 2 m ) 4 I ( 2 I 1 )
The angular average is evaluated in Equation (A5) over all the orientations of EFG principal component within the sediment. The first set of Wigner rotation matrix describes the orientation of the macroscopic clay film with respect to the magnetic field B 0 and the Wigner rotation matrix in the bracket characterizes the average orientation of EFG principal component within the sediment (see Figure 2c). For clay sediments with cylindrical symmetry, only the component p = 0 contributes to Equation (A5) which reduces to:
ω m 1 , m a p p = A m V z z P 3 cos 2 θ L F 1 2 3 cos 2 θ F P 1 2
For deuterium atoms pertaining to heavy water, the electrostatic field gradient V z z P is oriented along the O D director of the water molecule and corresponds to quadrupolar coupling constant ω Q = 990 × 10 3 rad/s [59].
In addition to the quadrupolar coupling, the heteronuclear dipolar coupling may also be responsible for the NMR relaxation of the confined probes because of the presence of paramagnetic impurities. The corresponding heteronuclear dipolar Hamiltonian [33,78,79] becomes:
H D ( t ) = C D m = 2 2 ( 1 ) m T 2 , m D , I R F 2 , m D , L ( t ) with C D = μ 0 4 π γ I γ S
where the C D is the dipolar coupling constant and the spin operators become [33,78,79]:
T 2 , 0 D , I R = 1 6 2 I z S z 1 2 ( I + S + I S + ) , T 2 , ± 1 D , I R = 1 2 ( I z S ± + I ± S z ) and T 2 , ± 2 D , I R = 1 2 I ± S
The functions F 2 , m D , L ( t ) in Equation (A7) are related to the second order spherical harmonics describing the reorientation of the vector joining the two coupled spins (noted r I S ( t ) ) by reference with the static magnetic field [33,78,79]:
F 2 , m D , L ( t ) = 24 π 5 Y 2 , m ( θ , φ ) r I S 3

Appendix B. NMR Relaxation Theory

In the framework of the Redfield theory [81], the time evolution of the spin quantum states, also called coherences, is described by the master equation [33,78,79]:
d σ * d t = i H S * , σ * + f ( σ * )
As noted by the asterisk ( * ), all terms are evaluated in the Larmor frequency rotating frame. The commutator describes the contribution from the static Hamiltonians H S * , including the excitation pulses and the residual quadrupolar Hamiltonian. The second term describes the contribution from the fluctuating parts of the quadrupolar and dipolar Hamiltonians:
H X F * ( t ) = C X m = 2 2 ( 1 ) m T 2 , m X , I R e i m ω 0 t F 2 , m X , L ( t ) F 2 , m X , L ( t )
where the index X stands for the various relaxation mechanisms (i.e., Q or D). This last contribution to the master equation is given by [33,78,79,81]:
f ( σ * ) = 0 t s u p H X F * ( t ) , e i H S * τ H X F * + ( t τ ) e i H S * τ , σ * ( t ) d τ
If the time-scales characterising the decorrelation of the various Hamiltonians are much smaller than the time-scale sampled by the evolution of the coherences, the upper limit of the integral t s u p in Equation (A12) may be set equal to infinity. This hypothesis restricts the validity of the Redfield theory applied to NMR relaxation [33,78,79].
Let us introduce the autocorrelation functions of the fluctuating components of the Hamiltonian
G m X , L ( τ ) = F 2 , m X , L ( 0 ) < F 2 , m X , L > × F 2 , m X , L ( τ ) < F 2 , m X , L > + F 2 , m X , L ( 0 ) < F 2 , m X , L > × F 2 , m X , L ( τ ) < F 2 , m X , L > with m { 0 , 1 , 2 }
By neglecting the time evolution of the coherences during the irradiation pulses, Equation (A12) becomes then [33,78,79]:
f ( σ * ) = m = 0 2 T 2 , m X , I R , T 2 , m X , I R , σ * J m X , L ( m ω 0 )
where the so-called spectral densities J m X , L ( m ω 0 ) satisfy the relationship [33,78,79]:
J m X , L ( m ω 0 ) = C X 2 0 G m X , L ( t ) e i m ω 0 t d t
The above-mentioned approximation is generally valid for classical relaxation measurements because the duration of the detection pulses (typically a few μ s) is much shorter than the time evolution of the coherences. Finally, a complete basis set of coherences is required to translate Equation (A14) into a matrix form [64,82,83,84,85].
In the case of spin-locking relaxation measurements, one cannot neglect the time evolution of the coherences during the irradiation power, since it is applied during all the evolution period of the coherences. As a consequence, the time evolution of the dipolar and quadrupolar couplings under the influence of the static Hamiltonian H S * must be taken into account in Equation (A12) as described implicitly by the term e i H S * τ H X F * + ( t τ ) e i H S * τ in the double commutator (see Equation (A12)). For that purpose, the static Hamiltonian is also formulated in a matrix form by using the complete basis set of coherences [20,42,44]. After evaluating its eigenvalues (noted ± i λ p ) and corresponding eigenvectors (noted v p ), one obtains a new complete basis set. The problem is then easily solved by projecting, into this eigenvectors basis set, the initial basis set of the coherences used to describe the T 2 , m X , I R spin operators.
By using the Wigner rotation matrices [35] (cf. Equation (A4)), it is possible to relate the derivation of the apparent correlation functions G m X , L ( τ ) , evaluated in the laboratory frame (noted L), with their intrinsic value evaluated in the frame attached to the clay sediment (noted F) [86]:
G 0 X , L ( τ ) = ( 1 3 cos 2 θ L F ) 2 4 G 0 X , F ( τ ) + 3 cos 2 θ L F sin 2 θ L F G 1 X , F ( τ ) + 3 ( 1 3 cos 2 θ L F ) 2 4 G 2 X , F ( τ )
G 1 X , L ( τ ) = 3 cos 2 θ L F sin 2 θ L F 2 G 0 X , F ( τ ) + 1 3 cos 2 θ L F + 4 cos 4 θ L F 2 G 1 X , F ( τ ) + 1 cos 4 θ L F 2 G 2 X , F ( τ )
G 2 X , L ( τ ) = 3 ( 1 cos 2 θ L F ) 2 8 G 0 X , F ( τ ) + 1 cos 4 θ L F 2 G 1 X , F ( τ ) + 1 + 6 cos 2 θ L F + cos 4 θ L F 8 G 2 X , F ( τ )
Obviously, the same relationship may be deduced for the corresponding spectral densities thanks to the linearity of the Fourier transform (see Equation (A15)).

Appendix C. Matrix Representation of the Irreducible Tensor Operators

Depending on the spin state, complete orthonormal basis sets may be constructed by using the irreducible tensor operators, also called coherences [87,88,89]. Symmetric and antisymmetric combinations of the coherences [83,84,85] are also introduced:
T l p ( s ) = 1 2 ( T l p + T l p ) and T l p ( a ) = 1 2 ( T l p T l p )
By using these new coherences, the three spin operators, I x , I y and I z , become proportional to T 11 ( a ) , T 11 ( s ) and T 10 respectively, simplifying the formulation of the Hamiltonians describing the irradiation pulse and the heteronuclear dipolar coupling (see Equation (A8)).
As displayed in Figure A1 and Figure A2, the size of the basis set increases significantly as a function of the spin of the nucleus.
Figure A1. Complete orthogonal basis set (four elements) describing the evolution of the spin I = 1/2 (Pauli matrices): (a) symbolic representation T i j ; (b) explicit matrix representation of irreducible tensor operators.
Figure A1. Complete orthogonal basis set (four elements) describing the evolution of the spin I = 1/2 (Pauli matrices): (a) symbolic representation T i j ; (b) explicit matrix representation of irreducible tensor operators.
Ijms 21 04697 g0a1
Figure A2. Complete orthogonal basis set (nine elements) describing the evolution of the spin I = 1: (a) symbolic representation T i j ; (b) explicit matrix representation of irreducible tensor operators.
Figure A2. Complete orthogonal basis set (nine elements) describing the evolution of the spin I = 1: (a) symbolic representation T i j ; (b) explicit matrix representation of irreducible tensor operators.
Ijms 21 04697 g0a2

Appendix D. Application to the Relaxation of Quadrupolar Nuclei

The differential equation describing the time evolution of the coherences (Equations (A10)–(A15)) may be written in a matrix form. To simplify the derivations of these matrices, we selected to use symmetric and antisymmetric combinations of the coherences [83,84,85]; see Equation (A17). By using these coherences, the three spin operators, I x , I y and I z , become proportional to T 11 ( a ) , T 11 ( s ) and T 10 respectively, simplifying the formulation of the Hamiltonians describing the irradiation pulse and the heteronuclear dipolar coupling (see Equations (A7)–(A9)).
For spin I = 1 nuclei, the time evolution of the coherences under the influence of the static Hamiltonian (noted H S * in Equation (A10)) including the residual quadrupolar coupling ( ω Q ) and the irradiation pulse ( ω 1 ) becomes [64,83]:
d d t T 20 T 11 ( a ) T 21 ( s ) T 22 ( s ) T 10 T 11 ( s ) T 21 ( a ) T 22 ( a ) = i 0 0 3 ω 1 0 0 0 0 0 0 0 ω Q 0 0 0 0 0 3 ω 1 ω Q 0 ω 1 0 0 0 0 0 0 ω 1 0 0 0 0 0 0 0 0 0 0 ω 1 0 0 0 0 0 0 ω 1 0 ω Q 0 0 0 0 0 0 ω Q 0 ω 1 0 0 0 0 0 0 ω 1 0 T 20 T 11 ( a ) T 21 ( s ) T 22 ( s ) T 10 T 11 ( s ) T 21 ( a ) T 22 ( a )
leading to two independent sub-sets of coherences. Analytical derivation of the eigenvalues ± i λ p with p { 0 , , 3 } and their corresponding eigenvectors v p are used to derive the general solutions of Equation (A18) [20,83]. Among others, two coherences are of practical interest for the interpretation of spin-locking relaxation measurements:
e i H S * τ T 20 e i H S * τ = 3 ω 1 ω Q 1 cos ( λ 1 τ ) λ 1 2 T 11 ( a ) + ω Q 2 + ω 1 2 ( 1 3 cos ( λ 1 τ ) ) λ 1 2 T 20 i 3 ω 1 sin ( λ 1 τ ) λ 1 T 21 ( s ) + 3 ω 1 2 cos ( λ 1 τ ) 1 λ 1 2 T 22 ( s ) and
e i H S * τ T 10 e i H S * τ = λ 2 cos ( λ 3 τ ) + λ 3 cos ( λ 2 τ ) λ 1 T 10 i ω 1 sin ( λ 3 τ ) + sin ( λ 2 τ ) λ 1 T 11 ( s ) + ω 1 cos ( λ 3 τ ) cos ( λ 2 τ ) λ 1 T 21 ( a ) i λ 2 sin ( λ 3 τ ) λ 3 sin ( λ 2 τ ) λ 1 T 22 ( a )
where the characteristic angular velocities λ 0 , λ 1 , λ 2 and λ 3 are defined by:
λ 0 = 0 , λ 1 = ω Q 2 + 4 ω 1 2 , λ 2 = λ 1 + ω Q 2 , and λ 3 = λ 1 ω Q 2 respectively .
As a consequence, under the simultaneous influences of the residual quadrupolar coupling ( ω Q ) and the irradiation pulse ( ω 1 ), the T 20 coherence oscillates according to the angular velocity λ 1 while two angular velocities ( λ 2 and λ 3 ) drive the oscillations of the T 10 coherence.
Straightforward calculations of the set of Equations (A14) and (A15) lead to the contributions of the quadrupolar [83] and heteronuclear dipolar [64] relaxation mechanisms to the time evolution of the coherences.
d d t T 20 T 11 ( a ) T 21 ( s ) T 22 ( s ) T 10 T 11 ( s ) T 21 ( a ) T 22 ( a ) = diag ( A , B , C , D , E , B , C , D ) × T 20 T 11 ( a ) T 21 ( s ) T 22 ( s ) T 10 T 11 ( s ) T 21 ( a ) T 22 ( a )
where A = A Q + A D , B = B Q + B D , C = C Q + C D , D = D Q + D D and E = E Q + E D , with
A Q = 3 J 1 Q ( ω 1 ) and A D = J 0 D ( ω 0 ω S ) 3 + J 1 D ( ω 0 ) + 2 J 2 D ( ω 0 + ω S )
B Q = 3 J 0 Q ( 0 ) 2 + 5 J 1 Q ( ω 0 ) 2 + J 2 Q ( 2 ω 0 ) and B D = 2 J 0 D ( 0 ) 9 + J 0 D ( ω 0 ω S ) 18 + J 1 D ( ω 0 ) 6 + J 1 D ( ω S ) 3 + J 2 D ( ω 0 + ω S ) 3
C Q = 3 J 0 Q ( 0 ) 2 + J 1 Q ( ω 0 ) 2 + J 2 Q ( 2 ω 0 ) and C D = 2 J 0 D ( 0 ) 9 + 5 J 0 D ( ω 0 ω S ) 18 + 5 J 1 D ( ω 0 ) 6 + J 1 D ( ω S ) 3 + 5 J 2 D ( ω 0 + ω S ) 3
D Q = J 1 Q ( ω 0 ) + 2 J 2 Q ( 2 ω 0 ) and D D = 8 J 0 D ( 0 ) 9 + J 0 D ( ω 0 ω S ) 9 + J 1 D ( ω 0 ) 3 + 4 J 1 D ( ω S ) 3 + 2 J 2 D ( ω 0 + ω S ) 3
E Q = J 1 Q ( ω 0 ) + 4 J 2 Q ( 2 ω 0 ) and E D = J 0 D ( ω 0 ω S ) 9 + J 1 D ( ω 0 ) 3 + 2 J 2 D ( ω 0 + ω S ) 3
Under slow modulation of the quadrupolar and heteronuclear dipolar couplings, i.e., when:
J 0 X ( 0 ) > > J m X ( ω 0 ) > > J m X ( ω S ) with m { 1 , 2 } and X { Q , D }
the set of Equations (A22a)–(A22e) reduces to [21,22]:
A = 3 U Q + U D B = C = 3 J 0 Q ( 0 ) 2 + 5 U Q 2 + 2 J 0 D ( 0 ) 9 + U D 2 D = 3 U Q + 8 J 0 D ( 0 ) 9 + U D E = 5 U Q + U D 3 where U Q = J 1 Q ( ω 0 ) J 2 Q ( 2 ω 0 ) and U D = J 0 D ( ω S ω 0 ) 3 + J 1 D ( ω 0 ) + 2 J 1 D ( ω S + ω 0 )
Four independent measurements of the relaxation of the T 20 , T 11 , T 22 and T 10 coherences lead then to the four dominant contributions ( U Q , J 0 Q ( 0 ) , U D and J 0 D ( 0 ) ) quantifying the quadrupolar and heteronuclear dipolar relaxation mechanisms [21,22].
The derivation of the time evolution of the coherences under spin-locking condition requires taking into account the evolution of the fluctuating part of the quadrupolar [83] and dipolar [20] Hamiltonians under the influence of the static Hamiltonians, as described by the term e i H S * τ H X F * + ( t τ ) e i H S * τ in Equation (A12). In the next approximation, we focus only on the m = 0 component of the fluctuating Hamiltonians because their m = 1 and m = 2 components oscillate at angular velocities ( ω 0 and 2 ω 0 ) (cf. Equation (A14)) much larger than the characteristic angular velocities ( λ i ) (cf. Equations (A19a)–(A20)). By using Equations (A19a) and (A19b), that approximation leads to:
d d t T 20 T 11 ( a ) T 21 ( s ) T 22 ( s ) T 10 T 11 ( s ) T 21 ( a ) T 22 ( a ) = A 3 K 0 0 0 0 0 0 0 B 0 2 K D 0 0 0 0 0 0 C 0 0 0 0 0 0 K Q + K D 0 D 0 0 0 0 0 0 0 0 E 0 K 0 0 0 0 0 0 L 0 2 K D 0 0 0 0 0 0 M 0 0 0 0 0 0 K Q + K D 0 D T 20 T 11 ( a ) T 21 ( s ) T 22 ( s ) T 10 T 11 ( s ) T 21 ( a ) T 22 ( a )
with
A Q = 3 J 1 Q ( ω 0 ) and A D = J 0 D ( ω 0 ω S ) 3 + J 1 D ( ω 0 ) + 2 J 2 D ( ω 0 + ω S )
B Q = 3 ω Q 2 J 0 Q ( 0 ) + 4 ω 1 2 J 0 Q ( λ 1 ) 2 λ 1 2 + 5 J 1 Q ( ω 0 ) 2 + J 2 Q ( 2 ω 0 ) and B D = 2 λ 2 J 0 D ( λ 3 ) + λ 3 J 0 D ( λ 2 ) 9 λ 1 + J 0 D ( ω 0 ω S ) 18 + J 1 D ( ω 0 ) 6 + J 1 D ( ω S ) 3 + J 2 D ( ω 0 + ω S ) 3
C Q = 3 ω Q 2 J 0 Q ( 0 ) + 4 ω 1 2 J 0 Q ( λ 1 ) 2 λ 1 2 + J 1 Q ( ω 0 ) 2 + J 2 Q ( 2 ω 0 ) and C D = 2 λ 2 J 0 D ( λ 3 ) + λ 3 J 0 D ( λ 2 ) 9 λ 1 + 5 J 0 D ( ω 0 ω S ) 18 + 5 J 1 D ( ω 0 ) 6 + J 1 D ( ω S ) 3 + 5 J 2 D ( ω 0 + ω S ) 3
D Q = J 1 Q ( ω 0 ) + 2 J 2 Q ( 2 ω 0 ) and D D = 8 J 0 D ( 0 ) 9 + J 0 D ( ω 0 ω S ) 9 + J 1 D ( ω 0 ) 3 + 4 J 1 D ( ω S ) 3 + 2 J 2 D ( ω 0 + ω S ) 3
E Q = J 1 Q ( ω 0 ) + 4 J 2 Q ( 2 ω 0 ) and E D = J 0 D ( ω 0 ω S ) 9 + J 1 D ( ω 0 ) 3 + 2 J 2 D ( ω 0 + ω S ) 3
L Q = 3 ω Q 2 J 0 Q ( 0 ) + 6 ω 1 2 J 0 Q ( 0 ) + J 0 Q ( λ 1 ) 2 λ 1 2 + 5 J 1 Q ( ω 0 ) 2 + 2 J 2 Q ( 2 ω 0 ) and L D = 2 λ 2 J 0 D ( λ 3 ) + λ 3 J 0 D ( λ 2 ) 9 λ 1 + J 0 D ( ω 0 ω S ) 18 + J 1 D ( ω 0 ) 6 + J 1 D ( ω S ) 3 + J 2 D ( ω 0 + ω S ) 3
M Q = 3 ω Q 2 J 0 Q ( 0 ) + 2 ω 1 2 J 0 Q ( 0 ) + J 0 Q ( λ 1 ) 2 λ 1 2 + J 1 Q ( ω 0 ) 2 + J 2 Q ( 2 ω 0 ) and M D = 2 λ 2 J 0 D ( λ 3 ) + λ 3 J 0 D ( λ 2 ) 9 λ 1 + 5 J 0 D ( ω 0 ω S ) 18 + 5 J 1 D ( ω 0 ) 6 + J 1 D ( ω S ) 3 + 5 J 2 D ( ω 0 + ω S ) 3
K Q = 3 ω 1 ω Q ( J 0 Q ( 0 ) J 0 Q ( λ 1 ) ) 2 λ 1 2 and K D = 2 ω 1 J 0 D ( λ 3 ) J 0 D ( λ 2 ) 9 λ 1
As mentioned above, the quadrupolar relaxation mechanism (implying the T 20 coherence) samples the spectral densities at the angular velocity λ 1 , while the heteronuclear dipolar relaxation mechanism (implying the T 10 coherence) samples the two other angular velocities ( λ 2 and λ 3 ), extending notably, the dynamical range probed by spin-locking relaxation measurements.

References

  1. Israelachvili, J.N. Intermolecular and Surface Forces; Academic Press: New York, NY, USA, 1985. [Google Scholar]
  2. Henderson, D. Fundamentals of Inhomogeneous Fluids; M. Dekker: New York, NY, USA, 1992. [Google Scholar]
  3. Mansoori, G.A.; Rice, S.A. Confined Fluids: Structure, Properties and Phase Behavior. In Advances in Chemical Physics; John Wiley & Sons, Ltd.: Hoboken, NJ, USA, 2014; Volume 156, Chapter 5; pp. 197–294. [Google Scholar] [CrossRef]
  4. Zax, D.B.; Yang, D.K.; Santos, R.A.; Hegemann, H.; Giannelis, E.P.; Manias, E. Dynamical heterogeneity in nanoconfined poly(styrene) chains. J. Chem. Phys. 2000, 112, 2945–2951. [Google Scholar] [CrossRef] [Green Version]
  5. Jobbagy, M.; Iyi, N. Interplay of charge density and relative humidity on the structure of nitrate layered double hydroxides. J. Phys. Chem. C 2010, 114, 18153–18158. [Google Scholar] [CrossRef]
  6. Lee, S.S.; Fenter, P.; Park, C.; Sturchio, N.C.; Nagy, K.L. Hydrated cation speciation at the muscovite (001)-water interface. Langmuir 2010, 26, 16647–16651. [Google Scholar] [CrossRef]
  7. Bowers, G.M.; Singer, J.W.; Bish, D.L.; Kirkpatrick, R.J. Alkali metal and H2O dynamics at the smectite/water interface. J. Phys. Chem. C 2011, 115, 23395–23407. [Google Scholar] [CrossRef]
  8. Boily, J.F. Water structure and hydrogen bonding at goethite/water interfaces: Implications for proton affinities. J. Phys. Chem. C 2012, 116, 4714–4724. [Google Scholar] [CrossRef]
  9. Ho, T.A.; Argyris, D.; Cole, D.R.; Striolo, A. Aqueous NaCl and CsCl solutions confined in crystalline slit-shaped silica nanopores of varying degree of protonation. Langmuir 2012, 28, 1256–1266. [Google Scholar] [CrossRef]
  10. Malani, A.; Ayappa, K.G. Relaxation and jump dynamics of water at the mica interface. J. Chem. Phys. 2012, 136, 194701. [Google Scholar] [CrossRef]
  11. Sato, K.; Fujimoto, K.; Dai, W.; Hunger, M. Molecular mechanism of heavily adhesive Cs: Why radioactive Cs is not decontaminated from soil. J. Phys. Chem. C 2013, 117, 14075–14080. [Google Scholar] [CrossRef]
  12. Buyukdagli, S.; Blossey, R. Dipolar correlations in structured solvents under nanoconfinement. J. Chem. Phys. 2014, 140, 234903. [Google Scholar] [CrossRef] [Green Version]
  13. Jiménez-Ruiz, M.; Ferrage, E.; Delville, A.; Michot, L.J. Anisotropy on the collective dynamics of water confined in swelling clay minerals. J. Phys. Chem. A 2012, 116, 2379–2387. [Google Scholar] [CrossRef]
  14. Michot, L.J.; Ferrage, E.; Delville, A.; Jiménez-Ruiz, M. Influence of layer charge, hydration state and cation nature on the collective dynamics of interlayer water in synthetic swelling clay minerals. Appl. Clay Sci. 2016, 119, 375–384. [Google Scholar] [CrossRef]
  15. Jiménez-Ruiz, M.; Ferrage, E.; Blanchard, M.; Fernandez-Castanon, J.; Delville, A.; Johnson, M.R.; Michot, L.J. Combination of inelastic neutron scattering experiments and ab initio quantum calculations for the study of the hydration properties of oriented saponites. J. Phys. Chem. C 2017, 121, 5029–5040. [Google Scholar] [CrossRef]
  16. Michot, L.J.; Delville, A.; Humbert, B.; Plazanet, M.; Levitz, P. Diffusion of water in a synthetic clay with tetrahedral charges by combined neutron time-of-flight measurements and molecular dynamics simulations. J. Phys. Chem. C 2007, 111, 9818–9831. [Google Scholar] [CrossRef]
  17. Michot, L.J.; Ferrage, E.; Jiménez-Ruiz, M.; Boehm, M.; Delville, A. Anisotropic features of water and ion dynamics in synthetic Na- and Ca-smectites with tetrahedral layer charge. A combined Quasi-Elastic Neutron-Scattering and Molecular Dynamics simulations study. J. Phys. Chem. C 2012, 116, 16619–16633. [Google Scholar] [CrossRef]
  18. Marry, V.; Dubois, E.; Malikova, N.; Durand-Vidal, S.; Longeville, S.; Breu, J. Water dynamics in hectorite clays: Influence of temperature studied by coupling neutron spin echo and molecular dynamics. Environ. Sci. Technol. 2011, 45, 2850–2855. [Google Scholar] [CrossRef]
  19. Delville, A.; Letellier, M. Structure and dynamics of simple liquids in heterogeneous condition: An NMR study of the clay-water interface. Langmuir 1995, 11, 1361–1367. [Google Scholar] [CrossRef]
  20. Porion, P.; Michot, L.J.; Faugère, A.M.; Delville, A. Influence of confinement on the long-range mobility of water molecules within clay aggregates: A 2H NMR analysis using spin-locking relaxation rates. J. Phys. Chem. C 2007, 111, 13117–13128. [Google Scholar] [CrossRef]
  21. Porion, P.; Michot, L.J.; Warmont, F.; Faugère, A.M.; Delville, A. Long-time dynamics of confined water molecules probed by 2H NMR multiquanta relaxometry: An application to dense clay sediments. J. Phys. Chem. C 2012, 116, 17682–17697. [Google Scholar] [CrossRef]
  22. Porion, P.; Faugère, A.M.; Delville, A. Multiscale water dynamics within dense clay sediments probed by 2H multiquanta NMR relaxometry and two-time stimulated echo NMR spectroscopy. J. Phys. Chem. C 2013, 117, 26119–26134. [Google Scholar] [CrossRef]
  23. Porion, P.; Faugère, A.M.; Delville, A. Structural and dynamical properties of water molecules confined within clay sediments probed by deuterium NMR spectroscopy, multiquanta relaxometry, and two-time stimulated echo attenuation. J. Phys. Chem. C 2014, 118, 20429–20444. [Google Scholar] [CrossRef]
  24. Porion, P.; Faugère, A.M.; Rollet, A.L.; Dubois, E.; Marry, V.; Michot, L.J.; Delville, A. Influence of strong confinement on the structure and dynamics of liquids: A study of the clay/water interface exploiting 2H NMR spectroscopy and spin-locking relaxometry. J. Phys. Chem. C 2018, 122, 16830–16841. [Google Scholar] [CrossRef]
  25. Porion, P.; Rodts, S.; Al-Mukhtar, M.; Faugère, A.M.; Delville, A. Anisotropy of the solvent self-diffusion tensor as a probe of nematic ordering within dispersion of nano-composite. Phys. Rev. Lett. 2001, 87, 208302. [Google Scholar] [CrossRef] [PubMed]
  26. Porion, P.; Al-Mukhtar, M.; Faugère, A.M.; Pellenq, R.J.M.; Meyer, S.; Delville, A. Water self-diffusion within nematic dispersion of nanocomposites: A multiscale analysis of 1H pulsed gradient spin-echo NMR measurements. J. Phys. Chem. B 2003, 107, 4012–4023. [Google Scholar] [CrossRef]
  27. Porion, P.; Al-Mukhtar, M.; Faugère, A.M.; Delville, A. 23Na Nuclear Magnetic Resonance and 1H Pulsed Gradient Spin-Echo detection of the critical concentration corresponding to the Isotrope/Nematic transition within aqueous dispersions of charged anisotropic nanoparticles. J. Phys. Chem. B 2004, 108, 10825–10831. [Google Scholar] [CrossRef]
  28. Porion, P.; Faugère, A.M.; Delville, A. 1H and 7Li NMR pulsed gradient spin echo measurements and multiscale modeling of the water and ionic mobility within aqueous dispersions of charged anisotropic nanoparticles. J. Phys. Chem. C 2008, 112, 11893–11900. [Google Scholar] [CrossRef]
  29. Porion, P.; Ferrage, E.; Hubert, F.; Tertre, E.; Dabat, T.; Faugère, A.M.; Condé, F.; Warmont, F.; Delville, A. Water mobility within compacted clay samples: Multi-scale analysis exploiting 1H NMR Pulsed Gradient Spin Echo and Magnetic Resonance Imaging of water density profiles. ACS Omega 2018, 3, 7399–7406. [Google Scholar] [CrossRef] [Green Version]
  30. Dabat, T.; Porion, P.; Hubert, F.; Paineau, E.; Dazas, B.; Grégoire, B.; Tertre, E.; Delville, A.; Ferrage, E. Influence of preferred orientation of clay particles on the diffusion of water in kaolinite porous media at constant porosity. Appl. Clay Sci. 2020, 184, 105354. [Google Scholar] [CrossRef]
  31. Kimmich, R. Field cycling in NMR relaxation spectroscopy: Applications in biological, chemical and polymer physics. Bull. Magn. Reson. 1979, 1, 195–218. [Google Scholar]
  32. Stapf, S.; Kimmich, R.; Seitter, R.O. Proton and deuteron field-cycling NMR relaxometry of liquids in porous glasses: Evidence for Lévy-walk statistics. Phys. Rev. Lett. 1995, 75, 2855–2858. [Google Scholar] [CrossRef]
  33. Kimmich, R. NMR: Tomography, Diffusometry, Relaxometry; Springer: Berlin/Heidelberg, Germany, 1997. [Google Scholar]
  34. Korb, J.P.; Whaley-Hodges, M.; Bryant, R.G. Translational diffusion of liquids at surfaces of microporous materials: Theoretical analysis of field-cycling magnetic relaxation measurements. Phys. Rev. E 1997, 56, 1934–1945. [Google Scholar] [CrossRef]
  35. Kimmich, R.; Anoardo, E. Field-cycling NMR relaxometry. Prog. Nucl. Magn. Reson. Spectrosc. 2004, 44, 257–320. [Google Scholar] [CrossRef]
  36. Korb, J.P. Nuclear magnetic relaxation of liquids in porous media. New J. Phys. 2011, 13, 035016. [Google Scholar] [CrossRef] [Green Version]
  37. Kruk, D.; Herrmann, A.; Rössler, E.A. Field-cycling NMR relaxometry of viscous liquids and polymers. Prog. Nucl. Magn. Reson. Spectrosc. 2012, 63, 33–64. [Google Scholar] [CrossRef]
  38. Steele, R.M.; Korb, J.P.; Ferrante, G.; Bubici, S. New applications and perspectives of fast field cycling NMR relaxometry. Magn. Reson. Chem. 2016, 54, 502–509. [Google Scholar] [CrossRef]
  39. Korb, J.P. Multiscale nuclear magnetic relaxation dispersion of complex liquids in bulk and confinement. Prog. Nucl. Magn. Reson. Spectrosc. 2018, 104, 12–55. [Google Scholar] [CrossRef]
  40. Faux, D.; Kogon, R.; Bortolotti, V.; McDonald, P. Advances in the interpretation of frequency-dependent Nuclear Magnetic Resonance measurements from porous material. Molecules 2019, 24, 3688. [Google Scholar] [CrossRef] [Green Version]
  41. Porion, P.; Faugère, A.M.; Delville, A. Long-distance water exchange within dense clay sediments probed by two-time 2H stimulated echo NMR spectroscopy. J. Phys. Chem. C 2013, 117, 9920–9931. [Google Scholar] [CrossRef]
  42. Porion, P.; Faugère, A.M.; Delville, A. Long-time scale ionic dynamics in dense clay sediments measured by the frequency variation of the 7Li multiple-quantum NMR relaxation rates in relation with a multiscale modeling. J. Phys. Chem. C 2009, 113, 10580–10597. [Google Scholar] [CrossRef]
  43. Porion, P.; Al-Mukhtar, M.; Meyer, S.; Faugère, A.M.; van der Maarel, J.R.C.; Delville, A. Nematic ordering of suspensions of charged anisotropic colloids detected by 23Na Nuclear Quadrupolar Spectroscopy. J. Phys. Chem. B 2001, 105, 10505–10514. [Google Scholar] [CrossRef]
  44. Porion, P.; Warmont, F.; Faugère, A.M.; Rollet, A.L.; Dubois, E.; Marry, V.; Michot, L.J.; Delville, A. 133Cs Nuclear Magnetic Resonance relaxometry as a probe of the mobility of cesium cations confined within dense clay sediments. J. Phys. Chem. C 2015, 119, 15360–15372. [Google Scholar] [CrossRef]
  45. Hwang, D.W.; Chu, C.C.; Sinha, A.K.; Hwang, L.P. Dynamics of supercooled water in various mesopore sizes. J. Chem. Phys. 2007, 126, 044702. [Google Scholar] [CrossRef]
  46. Sattig, M.; Vogel, M. Dynamic crossovers and stepwise solidification of confined water: A 2H NMR study. J. Phys. Chem. Lett. 2014, 5, 174–178. [Google Scholar] [CrossRef]
  47. Demuth, D.; Sattig, M.; Steinrücken, E.; Weigler, M.; Vogel, M. 2H NMR studies on the dynamics of pure and mixed hydrogen-bonded liquids in confinement. Z. Phys. Chem. 2018, 232, 1059–1087. [Google Scholar] [CrossRef]
  48. Weigler, M.; Brodrecht, M.; Buntkowsky, G.; Vogel, M. Reorientation of deeply cooled water in mesoporous silica: NMR studies of the pore-size dependence. J. Phys. Chem. B 2019, 123, 2123–2134. [Google Scholar] [CrossRef]
  49. Porion, P.; Faugère, A.M.; Michot, L.J.; Paineau, E.; Delville, A. Orientational microdynamics and magnetic-field-induced ordering of clay platelets detected by 2H NMR spectroscopy. Langmuir 2010, 26, 7035–7044. [Google Scholar] [CrossRef]
  50. Michot, L.J.; Bihannic, I.; Porsch, K.; Maddi, S.; Baravian, C.; Mougel, J.; Levitz, P. Phase diagrams of Wyoming Na-montmorillonite clay. Influence of particle anisotropy. Langmuir 2004, 20, 10829–10837. [Google Scholar] [CrossRef]
  51. Sato, T.; Watanabe, T.; Otsuka, R. Effects of layer charge, charge location, and energy change on expansion properties of dioctahedral smectites. Clays Clay Miner. 1992, 40, 103–113. [Google Scholar] [CrossRef]
  52. Delville, A.; Gasmi, N.; Pellenq, R.J.M.; Caillol, J.M.; Van Damme, H. Correlations between the stability of charged interfaces and ionic exchange capacity: A Monte Carlo study. Langmuir 1998, 14, 5077–5082. [Google Scholar] [CrossRef]
  53. Ferrage, E.; Lanson, B.; Malikova, N.; Plançon, A.; Sakharov, B.A.; Drits, V.A. New insights on the distribution of interlayer water in bi-hydrated smectite from X-ray diffraction profile modeling of 00l reflections. Chem. Mater. 2005, 17, 3499–3512. [Google Scholar] [CrossRef] [Green Version]
  54. Rinnert, E.; Carteret, C.; Humbert, B.; Fragneto-Cusani, G.; Ramsay, J.D.F.; Delville, A.; Robert, J.L.; Bihannic, I.; Pelletier, M.; Michot, L.J. Hydration of a synthetic clay with tetrahedral charges: A multidisciplinary experimental and numerical study. J. Phys. Chem. B 2005, 109, 23745–23759. [Google Scholar] [CrossRef]
  55. Ferrage, E.; Sakharov, B.A.; Michot, L.J.; Delville, A.; Bauer, A.; Lanson, B.; Grangeon, S.; Frapper, G.; Jiménez-Ruiz, M.; Cuello, G.J. Hydration properties and interlayer organization of water and ions in synthetic Na-smectite with tetrahedral layer charge. Part 2. Toward a precise coupling between molecular simulations and diffraction data. J. Phys. Chem. C 2011, 115, 1867–1881. [Google Scholar] [CrossRef]
  56. Woessner, D.E.; Snowden, B.S. NMR doublet splitting in aqueous montmorillonite gels. J. Chem. Phys. 1969, 50, 1516–1523. [Google Scholar] [CrossRef]
  57. Delville, A.; Grandjean, J.; Laszlo, P. Order acquisition by clay platelets in a magnetic field. NMR study of the structure and microdynamics of the adsorbed water layer. J. Phys. Chem. 1991, 95, 1383–1392. [Google Scholar] [CrossRef]
  58. Porion, P.; Faugère, A.M.; Michot, L.J.; Paineau, E.; Delville, A. 2H NMR spectroscopy and multiquantum relaxometry as a probe of the magnetic-field-induced ordering of clay nanoplatelets within aqueous dispersions. J. Phys. Chem. C 2011, 115, 14253–14263. [Google Scholar] [CrossRef]
  59. Edmonds, D.T.; Mackay, A.L. Pure quadrupole-resonance of deuteron in ice. J. Magn. Reson. 1975, 20, 515–519. [Google Scholar] [CrossRef]
  60. Delville, A. Structure of liquids at a solid interface: An application to the swelling of clay by water. Langmuir 1992, 8, 1796–1805. [Google Scholar] [CrossRef]
  61. Petit, D.; Korb, J.P.; Delville, A.; Grandjean, J.; Laszlo, P. Theory of nuclear spin relaxation in heterogeneous media and application to the cross correlation between quadrupolar and dipolar fluctuations of deuterons in clay gels. J. Magn. Reson. 1992, 96, 252–279. [Google Scholar] [CrossRef]
  62. Bloembergen, N.; Purcell, E.M.; Pound, R.V. Relaxation effects in Nuclear Magnetic Resonance adsorption. Phys. Rev. 1948, 73, 679–712. [Google Scholar] [CrossRef]
  63. Piculell, L. Water spin relaxation in colloidal systems. Part 1, 17O and 2H relaxation in dispersions of colloidal silica. J. Chem. Soc. Faraday Trans. 1 1986, 82, 387–399. [Google Scholar] [CrossRef]
  64. Porion, P.; Michot, L.J.; Faugère, A.M.; Delville, A. Structural and dynamical properties of the water molecules confined in dense clay sediments: A study combining 2H NMR spectroscopy and multiscale numerical modeling. J. Phys. Chem. C 2007, 111, 5441–5453. [Google Scholar] [CrossRef]
  65. Korb, J.P.; Delville, A.; Xu, S.; Demeulenaere, G.; Costa, P.; Jonas, J. Relative role of surface interactions and topological effects in nuclear magnetic resonance of confined liquids. J. Chem. Phys. 1994, 101, 7074–7081. [Google Scholar] [CrossRef]
  66. Pasquier, V.; Levitz, P.; Delville, A. 129Xe NMR as a probe of gas diffusion and relaxation in disordered porous media: An application to Vycor. J. Phys. Chem. 1996, 100, 10249–10256. [Google Scholar] [CrossRef]
  67. Levitz, P.E. Confined dynamics, forms and transitions in colloidal systems: From clay to DNA. Magn. Reson. Imaging 2005, 23, 147–152. [Google Scholar] [CrossRef]
  68. Woessner, D.E. Nuclear transfer effects in Nuclear Magnetic Resonance pulse experiments. J. Chem. Phys. 1961, 35, 41–48. [Google Scholar] [CrossRef]
  69. Delville, A.; Porion, P.; Faugère, A.M. Ion diffusion within charged porous network as probed by nuclear quadrupolar relaxation. J. Phys. Chem. B 2000, 104, 1546–1551. [Google Scholar] [CrossRef]
  70. Faux, D.A.; McDonald, P.J. Nuclear-magnetic-resonance relaxation rates for fluid confined to closed, channel, or planar pores. Phys. Rev. E 2018, 98, 063110. [Google Scholar] [CrossRef]
  71. Böhmer, R.; Jeffrey, K.R.; Vogel, M. Solid-state Li NMR with applications to the translational dynamics in ion conductors. Prog. Nucl. Magn. Reson. Spectrosc. 2007, 50, 87–174. [Google Scholar] [CrossRef]
  72. Böhmer, R.; Qi, F. Spin relaxation and ultra-slow Li motion in an aluminosilicate glass ceramic. Solid State Nucl. Magn. Reson. 2007, 31, 28–34. [Google Scholar] [CrossRef] [PubMed]
  73. Berk, N.F. Scattering properties of the leveled-wave model of random morphologies. Phys. Rev. A 1991, 44, 5069–5079. [Google Scholar] [CrossRef] [PubMed]
  74. Levitz, P. Off-lattice reconstruction of porous media: Critical evaluation, geometrical confinement and molecular transport. Adv. Colloid Interface Sci. 1998, 77, 71–106. [Google Scholar] [CrossRef]
  75. Cotts, R.M.; Hoch, M.J.R.; Sun, T.; Markert, J.T. Pulsed field gradient stimulated echo methods for improved NMR diffusion measurements in heterogeneous systems. J. Magn. Reson. 1989, 83, 252–266. [Google Scholar] [CrossRef]
  76. Callaghan, P.T. Principles of Nuclear Magnetic Resonance Microscopy; Clarendon Press: Oxford, UK, 1991. [Google Scholar]
  77. Hansen, J.P.; McDonald, I.R. Theory of Simple Liquids; Academic Press: London, UK, 1986. [Google Scholar]
  78. Abragam, A. The Principles of Nuclear Magnetism; Clarendon Press: Oxford, UK, 1961. [Google Scholar]
  79. Mehring, M. Principles of High Resolution NMR in Solids, 2nd ed.; Springer: Berlin/Heidelberg, Germany, 1983. [Google Scholar]
  80. Hertz, H.G. Magnetic relaxation by quadrupole interaction of ionic nuclei in electrolyte solutions. Part I: Limiting values for infinite dilution. Ber. Bunsen-Ges. Phys. Chem. 1973, 77, 531–540. [Google Scholar] [CrossRef]
  81. Redfield, A.G. On the theory of relaxation processes. IBM J. Res. Dev. 1957, 1, 19–31. [Google Scholar] [CrossRef]
  82. Porion, P.; Faugère, A.M.; Delville, A. 7Li NMR spectroscopy and multiquantum relaxation as a probe of the microstructure and dynamics of confined Li+ cations: An application to dense clay sediments. J. Phys. Chem. C 2008, 112, 9808–9821. [Google Scholar] [CrossRef]
  83. Van der Maarel, J.R.C. The relaxation dynamics of spin I = 1 nuclei with a static quadrupolar coupling and a radio-frequency field. J. Chem. Phys. 1993, 99, 5646–5653. [Google Scholar] [CrossRef]
  84. Van der Maarel, J.R.C. Thermal relaxation and coherence dynamics of spin 3/2. I. Static and fluctuating quadrupolar interactions in the multipole basis. Concepts Magn. Reson. Part A 2003, 19A, 97–116. [Google Scholar] [CrossRef]
  85. Van der Maarel, J.R.C. Thermal relaxation and coherence dynamics of spin 3/2. II. Strong radio-frequency field. Concepts Magn. Reson. Part A 2003, 19A, 117–133. [Google Scholar] [CrossRef]
  86. Barbara, T.M.; Vold, R.R.; Vold, R.L. A determination of individual spectral densities in a smectic liquid-crystal from angle dependent nuclear spin relaxation measurements. J. Chem. Phys. 1983, 79, 6338–6340. [Google Scholar] [CrossRef]
  87. Bowden, G.J.; Hutchison, W.D. Tensor operator formalism for multiple-quantum NMR. 1. Spin-1 nuclei. J. Magn. Reson. 1986, 67, 403–414. [Google Scholar] [CrossRef]
  88. Bowden, G.J.; Hutchison, W.D.; Khachan, J. Tensor operator formalism for multiple-quantum NMR. 2. Spins-3/2, spin-2, and spin-5/2 and general I. J. Magn. Reson. 1986, 67, 415–437. [Google Scholar] [CrossRef]
  89. Müller, N.; Bodenhausen, G.; Ernst, R.R. Relaxation-induced violations of coherence transfer selection rules in Nuclear Magnetic Resonance. J. Magn. Reson. 1987, 75, 297–334. [Google Scholar] [CrossRef]
Figure 1. Schematic view of the different time-scales explored by various experimental and numerical dynamical studies.
Figure 1. Schematic view of the different time-scales explored by various experimental and numerical dynamical studies.
Ijms 21 04697 g001
Figure 2. (a) 2H NMR spectra as a function of the film orientation θ L F in the static magnetic field B 0 ; (b) Snapshot illustrating one grand canonical Monte Carlo (GCMC) equilibrium configuration of confined water molecules and neutralising sodium counterions from which one can deduce concentration profiles of sodium counterions with oxygen and hydrogen atoms pertaining to water molecules confined between two beidellite clay lamellae; (c) Schematic view of multi-scale organisation of the clay sediment resulting from the coexistence of clay aggregates with various orientations of the platelet directors; (d) Variation of the residual quadrupolar coupling ν Q extracted from the 2H NMR spectra and the Hahn echo sequence as a function of the film orientation θ L F in the static magnetic field B 0 . Reprinted with permission from [23]. Copyright (2014) American Chemical Society.
Figure 2. (a) 2H NMR spectra as a function of the film orientation θ L F in the static magnetic field B 0 ; (b) Snapshot illustrating one grand canonical Monte Carlo (GCMC) equilibrium configuration of confined water molecules and neutralising sodium counterions from which one can deduce concentration profiles of sodium counterions with oxygen and hydrogen atoms pertaining to water molecules confined between two beidellite clay lamellae; (c) Schematic view of multi-scale organisation of the clay sediment resulting from the coexistence of clay aggregates with various orientations of the platelet directors; (d) Variation of the residual quadrupolar coupling ν Q extracted from the 2H NMR spectra and the Hahn echo sequence as a function of the film orientation θ L F in the static magnetic field B 0 . Reprinted with permission from [23]. Copyright (2014) American Chemical Society.
Ijms 21 04697 g002
Figure 3. Schematic view of: (a) the clay film within the NMR tube inserted into the detection coil (see text) and (b) the clay suspensions into the NMR tube. The X L and Z L directors are parallel to the axes of the NMR detection coil and the static magnetic field B 0 respectively. The θ L F Euler angle and the director n L F characterise the orientation of the clay sample by reference to the static magnetic field B 0 . Reprinted with permission from [23]. Copyright (2014) American Chemical Society.
Figure 3. Schematic view of: (a) the clay film within the NMR tube inserted into the detection coil (see text) and (b) the clay suspensions into the NMR tube. The X L and Z L directors are parallel to the axes of the NMR detection coil and the static magnetic field B 0 respectively. The θ L F Euler angle and the director n L F characterise the orientation of the clay sample by reference to the static magnetic field B 0 . Reprinted with permission from [23]. Copyright (2014) American Chemical Society.
Ijms 21 04697 g003
Figure 4. Pulse sequences and coherence transfer pathways used to measure the NMR relaxation rates of: (a) T 11 ( a , s ) ; (b) T 20 and (c) T 22 ( a , s ) coherences (see the Appendix C). The delay δ o p t is selected to optimise the coherence tranfer, and the delay ϵ is set equal to 10 μ s. Reprinted with permission from [23]. Copyright (2014) American Chemical Society.
Figure 4. Pulse sequences and coherence transfer pathways used to measure the NMR relaxation rates of: (a) T 11 ( a , s ) ; (b) T 20 and (c) T 22 ( a , s ) coherences (see the Appendix C). The delay δ o p t is selected to optimise the coherence tranfer, and the delay ϵ is set equal to 10 μ s. Reprinted with permission from [23]. Copyright (2014) American Chemical Society.
Ijms 21 04697 g004
Figure 5. (a) 2H NMR spectra recorded during the recovery of the residual quadrupolar splitting ν Q after a 90 ° rotation of a previously equilibrated clay sample; (b) Concentration variations of the limiting equilibrium residual quadrupolar splitting ν Q measured at B 0 = 2.351 and 8.465 Tesla for the Size2 clay samples and at B 0 = 8.465 Tesla for the Size3 clay samples. Additional measurements are performed at B 0 = 8.465 Tesla for the Size2 samples in the presence of salt (I = 10 4 M); (c) Determination of the various components of self-diffusion tensor D of the water molecules in the Size2 clay sample (c = 21.5 g/L) measured by the 2H NMR PGSE method at B 0 = 2.351 Tesla. The PGSE measurements are performed after equilibration of the clay platelet orientation in the static magnetic field B 0 parallel to the e z director; (d) Concentration variation of the characteristic time t c measured at B 0 = 2.351 and 8.465 Tesla for the Size2 clay samples. Reprinted with permission from [49]. Copyright (2010) American Chemical Society.
Figure 5. (a) 2H NMR spectra recorded during the recovery of the residual quadrupolar splitting ν Q after a 90 ° rotation of a previously equilibrated clay sample; (b) Concentration variations of the limiting equilibrium residual quadrupolar splitting ν Q measured at B 0 = 2.351 and 8.465 Tesla for the Size2 clay samples and at B 0 = 8.465 Tesla for the Size3 clay samples. Additional measurements are performed at B 0 = 8.465 Tesla for the Size2 samples in the presence of salt (I = 10 4 M); (c) Determination of the various components of self-diffusion tensor D of the water molecules in the Size2 clay sample (c = 21.5 g/L) measured by the 2H NMR PGSE method at B 0 = 2.351 Tesla. The PGSE measurements are performed after equilibration of the clay platelet orientation in the static magnetic field B 0 parallel to the e z director; (d) Concentration variation of the characteristic time t c measured at B 0 = 2.351 and 8.465 Tesla for the Size2 clay samples. Reprinted with permission from [49]. Copyright (2010) American Chemical Society.
Ijms 21 04697 g005
Figure 6. Influence of film orientation θ L F in the static magnetic field B 0 : (a) on the apparent multi-quanta relaxation rates of T 10 , T 11 ( a , s ) , T 20 and T 22 ( a , s ) coherences, noted R 10 , R 11 , R 20 and R 22 , respectively; and (b) on the on the apparent spectral densities J 0 Q ( 0 ) , J 0 D ( 0 ) , U Q and U D extracted from these R i j values; see Equation (A24) in Appendix D. The intrinsic spectral densities J 0 X , i n t r ( 0 ) , J 1 X , i n t r ( 0 ) and J 2 X , i n t r ( 0 ) with X { D , Q } , extracted from angular variation of the apparent spectral densities J 0 X ( 0 ) quantifying (c) heterogeneous quadripolar ( X = Q ) and (d) dipolar ( X = D ) relaxation mechanisms (see the Appendix). Reprinted with permission from [23]. Copyright (2014) American Chemical Society.
Figure 6. Influence of film orientation θ L F in the static magnetic field B 0 : (a) on the apparent multi-quanta relaxation rates of T 10 , T 11 ( a , s ) , T 20 and T 22 ( a , s ) coherences, noted R 10 , R 11 , R 20 and R 22 , respectively; and (b) on the on the apparent spectral densities J 0 Q ( 0 ) , J 0 D ( 0 ) , U Q and U D extracted from these R i j values; see Equation (A24) in Appendix D. The intrinsic spectral densities J 0 X , i n t r ( 0 ) , J 1 X , i n t r ( 0 ) and J 2 X , i n t r ( 0 ) with X { D , Q } , extracted from angular variation of the apparent spectral densities J 0 X ( 0 ) quantifying (c) heterogeneous quadripolar ( X = Q ) and (d) dipolar ( X = D ) relaxation mechanisms (see the Appendix). Reprinted with permission from [23]. Copyright (2014) American Chemical Society.
Ijms 21 04697 g006
Figure 7. Angular variations of the relative weights of the various intrinsic contributions J m X , i n t r ( ω ) to the m = 0 component of the apparent spectral density J 0 X ( ω ) monitoring the quadrupolar and dipolar relaxation mechanisms (see Equation (A16a) in the Appendix B).
Figure 7. Angular variations of the relative weights of the various intrinsic contributions J m X , i n t r ( ω ) to the m = 0 component of the apparent spectral density J 0 X ( ω ) monitoring the quadrupolar and dipolar relaxation mechanisms (see Equation (A16a) in the Appendix B).
Ijms 21 04697 g007
Figure 8. Pulse sequences and coherences transfer pathways used to measure the multi-quanta spin-locking relaxation rates of the (a) T 11 ( s ) ; (b) T 21 ( a ) and T 21 ( s ) ; and (c) T 22 ( a ) coherences, noted T 11 ρ ( s ) , T 21 ρ ( a ) , T 21 ρ ( s ) and T 22 ρ ( a ) respectively. The delay ϵ is set equal to 10 μ s. The pulse durations ψ 10 21 and ψ 10 22 are selected to optimise magnetisation transfer between the corresponding coherences; i.e., T 10 T 21 and T 10 T 22 respectively. Reprinted with permission from [23]. Copyright (2014) American Chemical Society.
Figure 8. Pulse sequences and coherences transfer pathways used to measure the multi-quanta spin-locking relaxation rates of the (a) T 11 ( s ) ; (b) T 21 ( a ) and T 21 ( s ) ; and (c) T 22 ( a ) coherences, noted T 11 ρ ( s ) , T 21 ρ ( a ) , T 21 ρ ( s ) and T 22 ρ ( a ) respectively. The delay ϵ is set equal to 10 μ s. The pulse durations ψ 10 21 and ψ 10 22 are selected to optimise magnetisation transfer between the corresponding coherences; i.e., T 10 T 21 and T 10 T 22 respectively. Reprinted with permission from [23]. Copyright (2014) American Chemical Society.
Ijms 21 04697 g008
Figure 9. Time evolution of (a) T 11 ( s ) ; (b) T 21 ( a ) ; (c) T 21 ( s ) and (d) T 22 ( a ) coherences measured under spin-locking conditions, noted T 11 ρ ( s ) , T 21 ρ ( a ) , T 21 ρ ( s ) and T 22 ρ ( a ) , respectively; and their Fourier transforms for (e) T 11 ( s ) ; (f) T 21 ( a ) ; (g) T 21 ( s ) and (h) T 22 ( a ) coherences (2H NMR). Reprinted with permission from [23]. Copyright (2014) American Chemical Society.
Figure 9. Time evolution of (a) T 11 ( s ) ; (b) T 21 ( a ) ; (c) T 21 ( s ) and (d) T 22 ( a ) coherences measured under spin-locking conditions, noted T 11 ρ ( s ) , T 21 ρ ( a ) , T 21 ρ ( s ) and T 22 ρ ( a ) , respectively; and their Fourier transforms for (e) T 11 ( s ) ; (f) T 21 ( a ) ; (g) T 21 ( s ) and (h) T 22 ( a ) coherences (2H NMR). Reprinted with permission from [23]. Copyright (2014) American Chemical Society.
Ijms 21 04697 g009
Figure 10. Variation of the intrinsic spectral densities J m Q , i n t r ( λ 1 ) with m { 0 , 1 } , describing quadrupolar coupling, and J m D , i n t r ( λ i ) with m { 0 , 1 } and i { 2 , 3 } , describing the heteronuclear dipolar coupling, as a function of the corresponding averaged angular velocities < λ i > (2H NMR). Reprinted with permission from [23]. Copyright (2014) American Chemical Society.
Figure 10. Variation of the intrinsic spectral densities J m Q , i n t r ( λ 1 ) with m { 0 , 1 } , describing quadrupolar coupling, and J m D , i n t r ( λ i ) with m { 0 , 1 } and i { 2 , 3 } , describing the heteronuclear dipolar coupling, as a function of the corresponding averaged angular velocities < λ i > (2H NMR). Reprinted with permission from [23]. Copyright (2014) American Chemical Society.
Ijms 21 04697 g010
Figure 11. Pulse sequence and coherence pathway used to measure the attenuation of the two-time 2H NMR stimulated echo I ( t e , τ M ) as a function of the evolution period t e and the mixing time τ M . The fourth pulse duration, noted ψ , is optimised in order to maximise the transfer from T 20 to T 22 ( s ) coherence and simultaneously minimise the transfer from T 10 to T 22 ( a ) coherence. Reprinted with permission from [23]. Copyright (2014) American Chemical Society.
Figure 11. Pulse sequence and coherence pathway used to measure the attenuation of the two-time 2H NMR stimulated echo I ( t e , τ M ) as a function of the evolution period t e and the mixing time τ M . The fourth pulse duration, noted ψ , is optimised in order to maximise the transfer from T 20 to T 22 ( s ) coherence and simultaneously minimise the transfer from T 10 to T 22 ( a ) coherence. Reprinted with permission from [23]. Copyright (2014) American Chemical Society.
Ijms 21 04697 g011
Figure 12. (a) Variation of the two-time stimulated echo attenuation I ( t e , τ M ) as a function of mixing time τ M (2H NMR). The data are normalised to take into account the relaxation of the T 20 coherence during the mixing time τ M (see Equation (4)); (b) Two-time correlation function extracted from the normalised stimulated echo attenuation as a function of the mixing time τ M . The red line corresponds to the best fit of a stretched exponential function, f ( t ) = A exp ( t / τ e x c h ) α , to determine the exchange time τ e x c h ( τ e x c h = 33 ± 5 ms with an exponent α set equal to 1.5). The green line dots are obtained by numerical modelling (see text). Reprinted with permission from [23]. Copyright (2014) American Chemical Society.
Figure 12. (a) Variation of the two-time stimulated echo attenuation I ( t e , τ M ) as a function of mixing time τ M (2H NMR). The data are normalised to take into account the relaxation of the T 20 coherence during the mixing time τ M (see Equation (4)); (b) Two-time correlation function extracted from the normalised stimulated echo attenuation as a function of the mixing time τ M . The red line corresponds to the best fit of a stretched exponential function, f ( t ) = A exp ( t / τ e x c h ) α , to determine the exchange time τ e x c h ( τ e x c h = 33 ± 5 ms with an exponent α set equal to 1.5). The green line dots are obtained by numerical modelling (see text). Reprinted with permission from [23]. Copyright (2014) American Chemical Society.
Ijms 21 04697 g012
Figure 13. Schematic view of the pulse sequence used to perform pulsed gradient spin-echo (PGSE-NMR) attenuation measurements. Details on the time delays, pulse durations and strengths of the magnetic field gradients G are given in the text. Reprinted with permission from [29]. Copyright (2018) American Chemical Society.
Figure 13. Schematic view of the pulse sequence used to perform pulsed gradient spin-echo (PGSE-NMR) attenuation measurements. Details on the time delays, pulse durations and strengths of the magnetic field gradients G are given in the text. Reprinted with permission from [29]. Copyright (2018) American Chemical Society.
Ijms 21 04697 g013
Figure 14. Water self-diffusion propagators measured within the saturated clay sample. More details on the selected diffusion directors e i are given (see text above), leading to the components of the self- diffusion tensor D 1 = 1.041 × 10 9 m 2 /s, D 2 = 1.062 × 10 9 m 2 /s, D 3 = 0.918 × 10 9 m 2 /s, D 4 = 1.067 × 10 9 m 2 /s, D 5 = 0.991 × 10 9 m 2 /s, and D 6 = 0.992 × 10 9 m 2 /s. Reprinted with permission from [29]. Copyright (2018) American Chemical Society.
Figure 14. Water self-diffusion propagators measured within the saturated clay sample. More details on the selected diffusion directors e i are given (see text above), leading to the components of the self- diffusion tensor D 1 = 1.041 × 10 9 m 2 /s, D 2 = 1.062 × 10 9 m 2 /s, D 3 = 0.918 × 10 9 m 2 /s, D 4 = 1.067 × 10 9 m 2 /s, D 5 = 0.991 × 10 9 m 2 /s, and D 6 = 0.992 × 10 9 m 2 /s. Reprinted with permission from [29]. Copyright (2018) American Chemical Society.
Ijms 21 04697 g014
Figure 15. Schematic view of the pulse sequence used to create the water concentration profiles inside the compacted clay sample by magnetic resonance imaging (MRI) technique. Reprinted with permission from [29]. Copyright (2018) American Chemical Society.
Figure 15. Schematic view of the pulse sequence used to create the water concentration profiles inside the compacted clay sample by magnetic resonance imaging (MRI) technique. Reprinted with permission from [29]. Copyright (2018) American Chemical Society.
Ijms 21 04697 g015
Figure 16. Time evolution of the water concentration profiles: (a) measured by MRI after the addition of bulk heavy water on the top of the water saturated sample and (b) obtained by numerical simulations of Brownian dynamics (see the text). Reprinted with permission from [29]. Copyright (2018) American Chemical Society.
Figure 16. Time evolution of the water concentration profiles: (a) measured by MRI after the addition of bulk heavy water on the top of the water saturated sample and (b) obtained by numerical simulations of Brownian dynamics (see the text). Reprinted with permission from [29]. Copyright (2018) American Chemical Society.
Ijms 21 04697 g016
Figure 17. Direct comparison of the time evolutions of the measured and simulated numbers of confined water molecules. Experimental and numerical data are fitted by a biexponential law: f ( t ) = m 1 exp ( t / m 2 ) + ( 1 m 1 m 4 ) exp ( t / m 3 ) + m 4 , where m i are constants. Reprinted with permission from [29]. Copyright (2018) American Chemical Society.
Figure 17. Direct comparison of the time evolutions of the measured and simulated numbers of confined water molecules. Experimental and numerical data are fitted by a biexponential law: f ( t ) = m 1 exp ( t / m 2 ) + ( 1 m 1 m 4 ) exp ( t / m 3 ) + m 4 , where m i are constants. Reprinted with permission from [29]. Copyright (2018) American Chemical Society.
Ijms 21 04697 g017
Table 1. Set of characteristic angular velocities ( λ 1 , λ 2 , λ 3 ) detected by multi-quanta spin-locking relaxometry for 2H NMR experiments obtained by varying the irradiation power ω 1 and the angle  θ L F  (see text).
Table 1. Set of characteristic angular velocities ( λ 1 , λ 2 , λ 3 ) detected by multi-quanta spin-locking relaxometry for 2H NMR experiments obtained by varying the irradiation power ω 1 and the angle  θ L F  (see text).
θ LF 0 ° 30 ° 90 ° 0 ° 30 ° 90 ° 0 ° 30 ° 90 °
ω 1 (10 5  rad/s) λ 1 (10 5  rad/s) λ 2 (10 5  rad/s) λ 3 (10 5  rad/s)
1.122 2.462.282.22 1.601.361.23 0.860.860.86
0.561 1.481.291.36 1.110.860.86 0.370.370.43
0.280 1.050.800.80 0.920.680.62 0.180.180.09
0.140 0.860.550.55 0.800.620.49 0.060.060.09
0.070 0.800.550.43 0.800.550.43 0.060.060.03

Share and Cite

MDPI and ACS Style

Porion, P.; Delville, A. A Multi-Scale Study of Water Dynamics under Confinement, Exploiting Numerical Simulations in Relation to NMR Relaxometry, PGSE and NMR Micro-Imaging Experiments: An Application to the Clay/Water Interface. Int. J. Mol. Sci. 2020, 21, 4697. https://doi.org/10.3390/ijms21134697

AMA Style

Porion P, Delville A. A Multi-Scale Study of Water Dynamics under Confinement, Exploiting Numerical Simulations in Relation to NMR Relaxometry, PGSE and NMR Micro-Imaging Experiments: An Application to the Clay/Water Interface. International Journal of Molecular Sciences. 2020; 21(13):4697. https://doi.org/10.3390/ijms21134697

Chicago/Turabian Style

Porion, Patrice, and Alfred Delville. 2020. "A Multi-Scale Study of Water Dynamics under Confinement, Exploiting Numerical Simulations in Relation to NMR Relaxometry, PGSE and NMR Micro-Imaging Experiments: An Application to the Clay/Water Interface" International Journal of Molecular Sciences 21, no. 13: 4697. https://doi.org/10.3390/ijms21134697

APA Style

Porion, P., & Delville, A. (2020). A Multi-Scale Study of Water Dynamics under Confinement, Exploiting Numerical Simulations in Relation to NMR Relaxometry, PGSE and NMR Micro-Imaging Experiments: An Application to the Clay/Water Interface. International Journal of Molecular Sciences, 21(13), 4697. https://doi.org/10.3390/ijms21134697

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop