NMR Investigation into the Influence of Surface Interactions on Liquid Diffusion in a Mesoporous Catalyst Support

Pulsed field gradient NMR diffusion measurements provide a non-invasive measure of the mass transport (self-diffusion) characteristics of liquids confined to porous catalyst materials. Here we explore the ability of this technique to probe the diffusive behaviour of a series of short-chain primary alcohols within a mesoporous catalyst support material; through the comparison of our results with highly surface-sensitive NMR relaxation data, we show that the evaluation of bulk-pore diffusion dynamics may provide a simple and indirect method to access and explore surface interaction phenomena occurring at the catalyst-liquid interface.


Introduction
It is well-established that the implementation of heterogeneous catalytic process is fundamental to the future of sustainable chemistry [1,2]. The development of liquid-phase reaction processes is of particular importance regarding the enhanced utilisation of biomass-waste organic matter derived from food refuse, animal waste and used vegetable oils-which may be processed to produce renewable platform chemicals and fuels [3][4][5]. While it is clear that relevant catalytic processes are necessary to facilitate such transformations [6][7][8][9], such reaction systems can differ significantly from the gas-phase processes encountered in more traditional catalyst applications [10]. Most notably, biomassderived compounds typically possess multiple polar functional groups; the resulting intermolecular interactions lead these materials to exhibit notably low volatilities and vapour pressures, such that effective transportation relies heavily on the addition of solvents [11]. It follows that the relevant chemistry and mass transport phenomena are dominated by the liquid-phase; while this is favourable for efficient macrokinetic transport processes, the high molecular density within such systems can lead to complex and competitive dynamics at the catalyst surface [11]. As a result, the development of robust approaches for the investigation of liquid-phase dynamics within catalytically relevant porous media has become a significant goal within modern catalytic research.
Nuclear magnetic resonance (NMR) relaxation and diffusion measurements have shown significant promise in this area [12]. Obtaining useful insight from traditional NMR chemical shift phenomena is challenging when considering liquid-saturated porous material as a result of significant line-broadening effects due to adsorption interactions and magnetic susceptibility differences at the solid-liquid interface. While such effects may be mitigated through the observation of nuclei known to exhibit a wide range of chemical shift values (e.g. 13 C), the resolution of such measurements is often limited by their inherently low natural abundance. The measurement of nuclear spin relaxation rates and molecular self-diffusion coefficients are, however, broadly independent of observable chemical shift phenomena. Rather, such measurements depend on the rate of decay of relevant NMR signals as a result of molecular dynamics. NMR relaxation studies have been particularly successful in determining adsorption phenomena within liquid-saturated 1 3 catalyst materials; such measurements probe the longitudinal ( T 1 ) and/or transverse ( T 2 ) relaxation time constants of the confined liquid. Notably, the ratio of these time constants T 1 ∕T 2 is now well-established as a measure of relative surface affinity exhibited by liquids and liquid mixtures at the catalyst pore surface [13][14][15], and we have recently extended such analysis to demonstrate that this ratio may be interpreted as a quantitative indicator of adsorption energetics [16,17].
Alternatively, pulsed field gradient (PFG) NMR diffusion measurements rely on the application of a series of short magnetic field gradient pulses to encode and then decode the position of nuclear spins either side of an observation period (see Sect. 2). Such measurements have been widely applied to the elucidation of mass transport phenomena in materials of relevance to heterogeneous catalysis [18], including mass-transfer limitations [19] and the influence of pore structure design [20][21][22]. Intriguingly, we note that liquids confined to porous materials will be subject to significant and repeated interactions with the pore walls, such that the measurement of their diffusive characteristics may present a potential method for the elucidation of surface interaction phenomena. While PFG NMR diffusion measurements have previously been applied to obtain direct insight into surface diffusion phenomena within liquid-saturated catalyst materials [23], the direct observation of dynamics associated with the adsorbed surface layer is challenging due to its small population and self-diffusion coefficient, as well as the inherently rapid rates of nuclear spin relaxation of species near the pore surface; such characteristics significantly reduce the observable signal from this layer, relative to that obtained from the bulk-pore population. It follows that this approach is only possible with access to NMR hardware capable of producing both large static magnetic fields and field gradient pulses ranging in the several T m −1 . Herein, with the aid of preliminary diffusion data obtained from a series primary alcohols within a silica catalyst support, we discuss how the observation of self-diffusion phenomena associated with the bulk-pore population of confined liquids may be utilised as an indirect method for the elucidation of surface phenomena within liquid-saturated mesoporous materials. This approach is predicted to be applicable to experiments using a range of NMR hardware, including low magnetic field benchtop systems [24].

Pulsed Field Gradient NMR
Given the general relationship between static magnetic field strength B 0 , gyromagnetic ratio and Larmor frequency 0 , the application of a spatially varying magnetic field (a gradient) = g x , g y , g z causes the precessional frequency of the spin system to vary with position [25], For a gradient applied in only the z-direction the frequency is and so when viewed from the rotating frame of referencethat being a frame of reference which precesses at 0results in a spatially-dependent offset [26] The time-evolution of the Larmor precession therefore exhibits a clear dependence on the spatial coordinate z . Indeed, provided that the sample magnetisation has been excited into the transverse plane through the application of an appropriate radio frequency (RF) pulse sequence, the spatial dependence of Ω(z) results in rotation of the transverse magnetisation by a phase factor exp (i ) . The phase is given by where t g defines the length of the applied gradient pulse. Neglecting any details of relaxation effects or chemical shift phenomena, and assuming the absence of any phase offset, the complex transverse magnetisation following this gradient pulse may then be expressed as where M 0 is the magnitude of the transverse magnetisation immediately prior to the pulse. Importantly, this expression describes a spatially varying magnetisation helix of pitch 2 g z t g −1 which encodes the position of the nuclear spins under study [27]. PFG NMR diffusion experiments extend the ability to encode the spatial positioning of nuclear spins to the measurement of molecular displacement rates. Such measurements are performed by employing multiple field gradient pulses separated by a diffusion observation time t Δ . Consider, for example, the displacement of a single spin by a factor Δz = z � − z in which the z-coordinate changes from z at time t = 0 to z ′ at t = t Δ . Displacement of the spin during the observation time results in a phase shift of the form The application of a second gradient pulse of equal magnitude but opposite polarity aims to remove the spatial encoding defined by Eq. (6). In the absence of any diffusive displacement Δz = 0 and the transverse magnetisation returns to its original state. If molecular motion has occurred, however, Δz ≠ 0 and the original magnetisation is not recovered. For a spin system of arbitrary size the corresponding NMR signal is given by where S(0) is the signal obtained in the absence of any applied field gradient. The propagator P Δ, t Δ describes the probability that a randomly selected spin will experience a displacement Δz during t Δ , and for normal diffusion is described by the Gaussian function where D is the self-diffusion coefficient which quantifies the rate of displacement. It follows that for Δz ≠ 0 diffusive attenuation of the signal S(g z ) will be observed. Combining Eqs. (8) and (9) we obtain the general expression [28] where t eff is an effective observation time which accounts for diffusion during the gradient pulses, and depends on the NMR pulse sequence employed. It follows that D may be obtained by incrementing g z through a range of values and fitting the acquired attenuation data to Eq. (10).

Diffusion in Restricted Systems
The diffusive behaviour of liquids confined to porous media is complicated by restrictions arising from the geometrical features of the accessible pore structure [29]. The application of PFG NMR measurements to such systems provides an effective self-diffusion coefficient D eff , which is dependent on the characteristics of the confining pore system and on the diffusion observation time employed [27,30]. Differing degrees of diffusive restriction may be characterised by the dimensionless parameter [25] where D 0 is the self-diffusion coefficient of the unrestricted liquid and is a characteristic length scale defining the pore structure. This expression provides a simple comparison between pore size and the mean squared displacement (MSD) of the diffusing species (Δz) 2 , which, for the evolution of a system over a time period t may be related to the self-diffusion coefficient of the system according to the Einstein expression, For restricted diffusion observed with short t Δ a shorttime diffusion limit emerges for which Eq. (11) gives Ξ ≪ 1 . Within this diffusive regime the majority of molecular displacements are not large enough for significant restrictions to occur from interactions with the pore walls, such that the observed effective self-diffusion coefficient tends towards D 0 . An intermediate regime emerges for longer t Δ values, and is characterised by Ξ ≈ 1 . This regime corresponds to a substantial portion of the diffusing molecules experiencing restrictions due to encounters with the pore walls, such that the MSD is reduced relative to that associated with unrestricted diffusion occurring over the same t Δ . An effective self-diffusion coefficient is observed which depends on both t Δ and the surface-to-volume ( S∕V ) of the confining pore structure according to [27] In the present work we concentrate on diffusive characteristics observed at long t Δ values, wherein all of the diffusing molecules are subject to significant interactions with the pore walls. This long-time diffusion limit is characterised by Ξ ≫ 1 . For isolated pores the MSD is limited to the square of the pore diameter and D eff = 0 ; for a well-connected pore network, however, the effective self-diffusion coefficient D eff t Δ → ∞ ≡ D ∞ is reduced from that of the unrestricted liquid according to [28] where is the tortuosity of the accessible pore network.
We note here that while Eq. (14) suggests that knowledge of the pore structure tortuosity may simply be obtained through the ratio D 0 ∕D ∞ , this is typically far from the case. Indeed, deviations in this ratio across a range of molecular liquids within the same pore structure are common, implying sensitivity to interactions between the probe liquid and porous matrix. The results of such measurements are instead described in terms of the so-called PFG interaction parameter [31], For weakly interacting non-viscous liquids (such as short-chain alkanes) restricted within mesoporous media this ratio is often considered equal to the tortuosity of the accessible pore network, = [32]. The observation that > may be facilitated by non-negligible adsorption interactions at the solid liquid interface, or by the occurrence of hindered diffusion, whereby region of the porous materials exhibit pore diameters of a similar size to the diffusion probe molecule employed. Prevalent examples of such are configurational diffusion through the pore network of zeolites, and the effects of coke deposition [33,34]. Conversely, the observation that < has been proposed to originate from the disruption of dynamic hydrogen bonding networks by the presence of the pore walls [20,32]. Herein, we present a brief investigation into the PFG interaction parameter presented by a homologous series of short-chain primary alcohols within an industrial silica catalyst support; the diffusion of cyclohexane is also investigated as a weakly interacting reference, and the assumption that such liquids may be used to provide an estimate of the pore structure tortuosity addressed.

Sample Preparation
A commercial G57 silica catalyst support was obtained from Johnson Matthey; material textural properties are summarised in Table 1. Methanol, ethanol, 1-propanol, 1-butanol and cyclohexane ( ≥ 99% purity, Sigma Aldrich) were used as received. The silica support material was first dried at 105 °C for 12 h. Imbibed samples were then prepared by soaking in excess liquid for at least 24 h under ambient conditions. The silica was then separated from each liquid and rolled over a pre-soaked filter paper to remove any liquid on the outer surface of the material. Finally, the imbibed material was transferred to sealed 5 mm NMR tubes to a height of approximately 15 mm, which were sealed with a plug of soaked filter paper so as to saturate the atmosphere within each NMR tube.

Pulsed Field Gradient NMR Measurements
A summary of the PFG NMR acquisition parameters employed in this work is provided in Table 2; all measurements were performed at 20.0 ± 0.1 °C and at ambient pressure. 1 H PFG NMR diffusion experiments were performed using a Bruker DMX spectrometer equipped with a 7.1 T superconducting magnet (corresponding to a 1 H frequency of 300.13 MHz) and a Bruker Biospin Diff-30 diffusion probe capable of producing magnetic field gradient pulses of up to 11.76 T m −1 . The diffusion of unrestricted liquids was analysed using the pulsed gradient stimulated echo (PGSTE) sequence illustrated in Fig. 1a [35]. Unrestricted self-diffusion coefficients D 0 were obtained from this analysis by fitting the acquired signal attenuation data to the Stejskal-Tanner equation [35], where Our PGSTE measurements were performed by holding t g = 1 ms constant and varying the magnetic field gradient strength; 16 linearly spaced g z values were employed while the observation time was set to t Δ = 50 ms. Trapezoidal gradient pulses of area t g g z were employed to ensure consistent pulse shaped across all values of g z . A homospoil   16 32 gradient of magnitude g max ∕3 (where g max is the magnitude of the maximum applied gradient) and length t H = 5 ms was also applied during the storage interval T to remove any remaining coherent transverse magnetisation. The echo time was e = 3.2 ms. The diffusion of liquids restricted to mesoporous silica was analysed using the alternating pulsed gradient stimulated echo (APGSTE) sequence illustrated in Fig. 1b [36]; this sequence was employed to minimise the influence of any background magnetic field gradients resulting from susceptibility differences at the solid-liquid interface [37]. The effective self-diffusion coefficients of restricted liquids D eff were therefore obtained by fitting the acquired signal attenuation data to where [36] In analogy to our PGSTE experiments, APGSTE analysis was carried out by holding t g = 1 ms constant and varying g z across 16 linearly spaced values; in this case, however, the observation time was set to t Δ = 100 ms to allow for sufficient self-diffusion throughout the silica pore network. Trapezoidal gradient pulse lobes (in this case of area t g g z ∕2 ) and a homospoil gradient of magnitude g max ∕3 and length t H = 10 ms were employed, while the echo time was e = 2.7 ms.

Results and Discussion
Our acquired PGSTE and APGSTE diffusion attenuation data are show in Fig. 2a and b, respectively. The resulting self-diffusion coefficients are detailed in Table 3. In all cases the diffusion observation time for our APGSTE analysis was set to t Δ = 100 ms; given the resultant effective selfdiffusion coefficients D eff obtained from these experiments we may readily calculate that the corresponding root mean squared displacements ( RMSD = √ 2t Δ D eff ) of the confined liquids fall within the range 0.64-1.59 µm. Clearly these displacements are far larger than the average pore size of the silica material employed here ( d pore = 15 nm), resulting in molecular self-diffusion across a significant number of pore diameters during this observation time. In turn, this suggests that we may define D eff t Δ = 100 ms ≈ D ∞ . Calculation of the restricted diffusion parameter Ξ according to Eq. (11) supports this assumption, with each of the five liquids investigated here exhibiting Ξ ≫ 1 , indicative of the long-time diffusion limit. It follows that our acquired D eff values may be used to calculate the PFG interaction parameter according to Eq. (15); these value are also provided in Table 3.
It is clear from our calculated values that the primary alcohols investigated here exhibit larger D 0 ∕D eff values than expected from the tortuosity of our silica material, which may be approximated from our cyclohexane diffusion data, such that ≈ ≈ 1.5 . This value is in good agreement with previous investigations into the tortuosity of similar mesoporous catalyst support materials [32]. The observed increase in away from this value suggest that adsorption interactions at the solid-liquid interface have a non-negligible effect on the self-diffusion of our primary alcohols within this material. Furthermore, there is a clear increase in with increasing alcohol carbon chain length; notably, this trend is similar to that previously observed when considering the nuclear spin relaxation characteristics ( T 1 ∕T 2 ratio) of these same adsorbate/adsorbent systems at intermediate magnetic field [17]. This observation is of particular interest as we have shown that the T 1 ∕T 2 ratio obtained under such conditions may be directly related to adsorption energetics [16]. As such, we provide in Fig. 3 a direct comparison of our calculated data with these T 1 ∕T 2 values. We find that there exists a remarkably strong, positive correlation between these two measurements. Given the well-established sensitivity of the ratio T 1 ∕T 2 to surface phenomena, this observation implies that the restricted diffusion observed through our APGSTE analysis exhibits significant sensitivity to adsorption interactions at the pore surface. The solid diagonal line within Fig. 3 illustrates a linear fit to the data and suggests that D 0 ∕D ∞ ≈ T 1 ∕T 2 ∕10 + 1.3 . Of particular interest is the solution to this equation for T 1 ∕T 2 = 1 , indicative of unrestricted-and therefore totally Fig. 1 PFG NMR diffusion pulse sequence diagrams for a PGSTE and b APGSTE analysis. Each pulse sequence consists of a radio frequency (RF) and gradient pulse ( g ) axis. All remaining symbols are defined within the main text and Table 2 1 3 non-interacting-liquids [38]. In this case our data fit predicts that in the absence of any adsorption interactions D 0 ∕D ∞ ≈ 1.4 . This value is clearly in excellent agreement with the tortuosity implied through our cyclohexane diffusion data, providing validation for the use of short chain alkanes as a measure of tortuosity in such systems.
It is appropriate here to briefly interrogate how our selfdiffusion measurements might exhibit such sensitivity to adsorption phenomena. To a simple approximation our liquid-saturated silica may be treated as a biphasic system, consisting of an adsorbed surface layer and a bulk-like liquid population towards the centre of the pores. The observed diffusion characteristics of such systems is highly dependent upon the rate of exchange between the two populations, relative to the diffusion observation time employed. If the observation time is short with respect to this exchange process (biphasic slow-exchange) then the resulting PFG data will attenuate as (20) S g z S(0) = P 1 exp −bD 1 + P 2 exp −bD 2 .
Here S g z ∕S(0) is the normalised signal attenuation, D 1 and D 2 are the self-diffusion coefficients associated with Fig. 2 PFG NMR signal attenuation data for short-chain primary alcohols and cyclohexane in a the unrestricted bulk and b restricted within mesoporous silica. Solid lines in a represent a fit to Eqs. (16) and (17) while solid lines in b represent a fit to Eqs. (18) and (19). The resultant self-diffusion coefficients are detailed in Table 3   Table 3 Summary of the PFG NMR diffusion characteristics obtained from Fig. 2 D eff values correspond to the effective self-diffusion of these liquids in mesoporous silica ( d pore = 15 nm) with the observation time set to t Δ = 100 ms. An estimate of the restricted diffusion parameter Ξ is given according to Eq. (11) assuming = d pore ∕2 . The PFG interaction parameter is calculated according to Eq. (15) assuming D eff t Δ = 100 ms ≡ D ∞ . T 1 ∕T 2 values are taken from Ref [17] Liquid D 0 × 10 10 (m 2 s −1 )  . 3 Comparison of the PFG interaction parameter with ratio of NMR relaxation time constants T 1 ∕T 2 (a measure of surface affinity) obtained for these same adsorbate/adsorbent systems. NMR relaxation data is taken from Ref [17] populations P 1 and P 2 respectively (with P 1 = 1 − P 2 ), and b depends on the NMR pulse sequence employed. Data of this form will be highly apparent in the resulting logattenuation plot as a curve with limiting gradients −D 1 and −D 2 . [23] Conversely, if the diffusion observation time is long with respect to the exchange process (biphasic fastexchange) then the resulting PFG attenuation data will follow with In this case a single self-diffusion coefficient D is observed, defined by the weighted average of the rates of diffusion within the two populations. Given the long diffusion observation time employed in our APGSTE measurements it is reasonable to assume that this biphasic fast-exchange approach is applicable to the self-diffusion processes explored here; this interpretation is supported by the log-attenuation plots in Fig. 2, which present a single diffusion coefficient as characterised by a straight line on the log-scale.
In the present case we must also take account of the influence of the pore structure tortuosity on our observed self-diffusion characteristics. In particular, the effective self-diffusion coefficient of our restricted liquids may be defined according to Here D 0 is again the self-diffusion coefficient of the unrestricted bulk, while D surf is the self-diffusion coefficient for species within an adsorbed surface layer of population P ; is again the tortuosity of the overall pore structure while surf is an apparent surface tortuosity which influences diffusion through the adsorbed surface layer [39]. This tortuosity parameter relates to the topology of the pore surface on the length scale s ∼ √ D surf s [30], where s is the mean surface residence time of the adsorbate at the pore surface. It follows from this expression that for (1 − P)D 0 ∕ ≫ PD surf ∕ surf we obtain D 0 ∕D ∞ ≈ . However, if the surface term PD surf ∕ surf is non-negligible then D 0 ∕D ∞ > . It follows that our observed trend in PFG interaction parameter might be explained through the increased prominence of this surface term with increasing alcohol carbon chain length. Furthermore, we note that for a series of similar adsorbates (e.g. short-chain primary alcohols) within the same mesoporous material, the prominence of this surface term will be dominated by the surface diffusion coefficient D surf . This diffusion coefficient provides a potential link to adsorption energetics through the usual Arrhenius form where E a is the activation energy for the surface diffusion process, typically considered some fraction of the adsorption energy [40], R is the gas constant, T the absolute temperature, and D surf ,0 the Arrhenius pre-exponential factor for activated surface diffusion. While this expression does not provide a direct explanation for the strong linear relationship observed between the two data sets in Fig. 3, perturbation of the observed self-diffusivity behaviour by activated surface dynamics clearly presents a promising link between these different measurements. Indeed, recalling that the surface diffusion coefficient may also be expressed in terms of the average adsorption site residence time m and the mean distance between adsorption sites as D surf = 2 ∕4 m [41], and that the expression T 1 ∕T 2 ∝ −1∕ ln m has been derived elsewhere [16], some simple algebraic manipulation suggests that the relationship between surface relaxation and self-diffusion should be exponential in nature, rather than the linear expression obtained from the data in Fig. 3. We note, however, that a far more detailed evaluation of this relationship over a wider range of data points is necessary to determine whether such an expression is robust; such experiments are the subject of future investigation.

Conclusion
This work has detailed a brief evaluation of the diffusive characteristics of a series of short-chain primary alcohols and cyclohexane within a commercial mesoporous silica catalyst support material. Through the application of PFG NMR diffusion measurements we have revealed a distinct inequality in the ratio of unrestricted-to-restricted selfdiffusion coefficients across the series of liquids investigated. Indeed, while this ratio is often defined as a measure of the tortuosity of the accessible pore space, our results illustrate a notable increase in this parameter with increasing carbon chain length. A direct comparison of this data with the results of nuclear spin relaxation measurements reported previously reveal a strong, positive correlation, and suggests that our NMR diffusion data is sensitive to adsorption interactions occurring at the solid-liquid interface. A biphasic fast-exchange theory has been presented which suggests that the sensitivity of bulk-pore diffusivity to surface phenomena occurs as a result of the influence of a surface diffusion term. Overall, the results described here demonstrate that such measurements may be a potential (24) D surf = D surf ,0 exp −E a RT , source of additional insight for providing rational connections between adsorption phenomena and mass transfer characteristics within systems of relevance to heterogeneous catalysis.