surface interactions of liquids in porous media are of great importance , particularly in the field of heterogeneous catalysis,[15 ] and the ability to understand surface interactions is essential for efficient and rational catalyst design . however ,
probing liquid  surface interactions in liquid - saturated porous media is particularly challenging . established techniques for probing the interaction of molecules at solid surfaces include isosteric heat of adsorption , temperature - programmed desorption ( tpd ) , infra - red ( ir ) spectroscopy , and nuclear magnetic resonance ( nmr ) chemical shift (  ) . however , all these measurements have limitations , and none are able to probe , non - destructively , the behaviour of molecules on catalyst surfaces at realistic reaction conditions .
for example , whilst isosteric heat of adsorption measurements are non - invasive , they require the determination of adsorption isotherms at different temperatures , and are therefore extremely time - consuming ; characterisation of co - adsorbed systems may take in excess of well over ten hours .
tpd is known to cause in situ reactions of some organic molecules leading to dehydration , isomerisation and decomposition , and is also limited in its ability to probe co - adsorption.[9 , 10 ] ir is well known to be limited to liquids and solids that are transparent in the frequency range of interest .
finally , nmr chemical shift interpretation is robust only for monolayer surface coverage.[12 , 13 ] it has recently been reported that under certain conditions chemical shift measurements can be used to probe surface interactions over pvp - stabilised metal nanoparticles . in this example , the chemical shift of formic acid ( adsorbate ) on metal colloid catalysts , measured by c nmr spectroscopy in aqueous suspension , was used to infer the strength of surface interaction . however , to achieve this measurement it was necessary to avoid direct contact between the c atom of the adsorbate and the metal surface by inserting oxygen - atom spacers , in order to eliminate spectral line broadening . whilst this is clearly an elegant measurement ,
it can not readily be transferred to measurements of liquid  surface interactions in saturated porous materials of industrial relevance . in general , the use of nmr methods to probe adsorption is made challenging because of the broadening of the nmr lineshape associated with the adsorbed species , which arises because of the magnetic susceptibility differences between adsorbate and adsorbent . there exist many reports of nmr experiments probing different aspects of adsorption .
c nmr is often the method of choice because of the wider chemical shift range , and hence spectral resolution attainable , relative to h observation . examples of c nmr to study adsorption include chemical shift studies to investigate the structure of the adsorbed species and the determination of acid site strength ( mostly in zeolites ) using suitable probe molecules such as acetonitrile or acetone.[16 , 17 ] c nmr approaches usually require the use of enriched c species due to the low signal - to - noise ratio of c signal at natural abundance .
however , it is noted that c studies of adsorbed species at natural isotopic abundance using the distortionless enhancement by polarisation transfer ( dept ) technique have been reported , although the strength of interaction was not probed using this technique .
h nmr has been used by various researchers to study acidic strength of hydroxyl groups in zeolites , using deuterated probe molecules such as acetonitrile and pyridine.[1719 ] for example , zheng et al .
have conducted experimental and theoretical studies of deuterated pyridine adsorption to investigate acidic strength of solid acids and found a linear correlation between the h chemical shift of adsorbed pyridine and the proton affinity .
thus , whilst various implementations of nmr spectroscopy have been used to study adsorption phenomena in porous media , a robust , generic technique readily applicable to co - adsorption and with the ability to study liquid - saturated pore spaces , as opposed to monolayer coverages , has not been presented .
the method proposed in the present work exploits the modification of the nuclear spin relaxation time characteristics of the adsorbate that result from its interaction with the surface of the pore .
the advantage of using so - called nmr relaxation measurements is that the characterisation of the adsorption interaction does not rely on the nmr lineshape and the 
peak positionthe actual peak position associated with liquid confined within a porous medium , or chemical shift , may be influenced by factors other than the adsorbate 
moreover , measurements based on chemical shift will be limited in terms of the systems that can be studied when multi - component liquid systems inside the pore space are of interest , because of the broad , overlapping , line shapes associated with each component .
in contrast , the relaxation times of different chemical species in a multi - component system are often quite distinct even when their respective chemical shifts overlap .
thus far[2124 ] it has been used as a qualitative probe yielding the relative strength of interactions between different adsorbate  adsorbent systems .
it is the purpose of this paper to demonstrate that the ratio of the spin - lattice to spin - spin ( or transverse ) relaxation time ( t1/t2 ) can be related directly to the activation energy of desorption characterising the strongest adsorption sites on the surface of the adsorbent , as determined by temperature programmed desorption ( tpd ) . in recent years
, nmr relaxation has emerged as a non - invasive , chemically sensitive technique for studying surface interactions of liquids in saturated porous media . following radio - frequency ( rf ) excitation ,
longitudinal t1 relaxation processes drive the longitudinal magnetisation to the equilibrium position , that is , aligned with the external magnetic field , whilst transverse t2 relaxation processes determine the rate of loss of phase coherence of the magnetisation in the transverse plane . reduced t1 and t2 relaxation times
are observed when liquid molecules adsorb on a solid surface due to a change in the molecular mobility ; in bulk liquids , t1t2 . both t1 and
t2 are affected by changes in the rotational correlation time of the adsorbate molecules . however , t2 is further influenced by a translational correlation time associated with surface diffusion.[27 , 28 ] consequently , when molecules adsorb on surfaces , changes in their translational and rotational dynamics influence t2 more than t1 , resulting in t1>t2 . in porous materials with a high surface - to - volume ratio ,
s / v , the observed relaxation rates are proportional to s / v , so absolute t1 and t2 measurements can not be readily used to compare interactions between materials with differing pore geometry , pore size and density of adsorption sites . however , the ratio of relaxation times t1/t2 is ( to leading order ) independent of these characteristics .
it has been observed empirically that the ratio t1/t2 provides an indication of the relative strength of surface interaction for different liquids in the same catalyst .
for example , in the aforementioned work , it was observed that water imbibed in pd / al2o3 catalyst trilobes exhibits a much larger t1/t2 ratio ( stronger surface interaction ) than 2-butanone on the same surface ; water is known to poison this palladium catalyst in the hydrogenation of 2-butanone by preferential adsorption on the active sites .
other studies , which exploit t1/t2 as a probe of surface interaction strength and the ability of these measurements to help understand catalytic performance , have recently been reported.[2123 ] it is important to note that besides the advantage of being independent of pore geometry , the use of t1 and t2 relaxation times is particularly advantageous if co - adsorption in porous materials is to be probed .
this is because the method has the advantage that different species can often be separated based on their relaxation time values even when the two species can not be separated in the chemical shift domain , which is often the case in porous materials .
this was shown when studying co - adsorption of various binary mixtures in porous al2o3 and sio2 supported catalysts .
the aim of the present work is to confirm theoretically that the ratio t1/t2 can be used as a qualitative descriptor of surface affinity , and extend the measurement to provide a quantitative metric of the same .
moreover , if this is the case , can the t1/t2 measurement be shown to correlate with an accepted laboratory characterisation of the adsorption interaction , namely tpd ?
indeed , if the theoretical interpretation summarised below is appropriate , we might expect the value of t1/t2 to correlate with the strongest adsorption energies characterising a particular system , since nuclear spin relaxation is dominated by the strongest relaxation sinks present . in this section
we establish proportionality between nmr relaxation and adsorption energy for a molecule interacting with the surface of a pore ; the relaxation of the adsorbed molecule will be modulated by the reduction in molecular mobility and the dipole 
the present analysis develops from the established theory of surface relaxation , which assumes the pore surface contains paramagnetic impurities that act as relaxation sinks and it is dipole - electron coupling that determines t1/t2 . in recent work we demonstrated that the analysis of godefroy et al . for proton
electron dipole coupling is extensible to proton  proton dipole coupling , albeit with a much weaker coupling constant . in the absence of paramagnetic species , it is the active binding sites on the pore surface that are considered to be the relaxation sinks .
molecules adsorbed onto a surface undergo two - dimensional ( 2d ) translational motion governed by the activation energy em .
the diffusion coefficient associated with this surface motion is : in which dm0 is a temperature independent contribution to surface diffusion , r is the ideal gas constant , and t is the temperature . in liquid - saturated porous materials , adsorbed molecules
the effective surface diffusion coefficient deff ( observed ) is modified from dm due to the finite residence time on the surface .
therefore , the effective diffusion coefficient is : in which e = emes is an activation energy for surface diffusion .
the diffusion coefficient has an associated correlation time m that describes the time for diffusion between active binding sites , such that deff(t)/(4m ) , where  is the thickness of the adsorbed surface layer .
thus , the surface correlation time is defined as : there will also be a surface residence time for adsorbed molecules , which is similarly defined as : the observed relaxation times for liquid saturated porous media are : in which 1,2 are the surface relaxivity values governing longitudinal and transverse relaxation processes at the pore surface , respectively , t1,2,bulk are the respective relaxation times of the bulk liquid , and p is the fraction of spins on the pore surface . for the catalyst supports of interest , the pore diameters are on the order of nanometres , such that the observed relaxation times are determined only by surface relaxation .
the ratio of observed ( surface ) relaxation times is defined by the spectral density function j(0 ) as : at high field , this ratio is known to be insensitive to larmor frequency .
based on typical measured values of m , s for liquids on an alumina surface , and given that 0 is large , we assume that ( m/s)(0m)1 . through a series of trivial algebraic steps
suggest that t1/t2 is sensitive to small changes in ln(m ) but insensitive to large changes in ln(s/m ) as m , s are not independent variables .
therefore , for the purpose of establishing a relationship between surface relaxation and adsorption energy it follows from equation ( 3 ) and the reduced form of equation ( 6 ) , given 0 is constant , that : we define a dimensionless surface interaction parameter esurf=t2/t1 for the purpose of interpreting surface relaxation as a surface energy . for the materials of interest in this study , we are concerned only with the dynamics of the adsorbed liquid that give rise to altered nmr relaxation times .
however , it is important to consider the extent to which other properties of the sample will influence the observed relaxation times , and the t1/t2 ratio .
the most important of these are 1 ) the presence of paramagnetic species on the pore surface , and 2 ) the presence of magnetic susceptibility contrast between the solid and liquid .
first , considering the case of presence of paramagnetic species on the pore surface ; such species will alter the sensitivity of the relaxation mechanism to m and s .
empirical evidence suggests the ratio t1/t2 remains independent of the density of surface paramagnetic species ( such as fe ) when the concentration of the paramagnetic impurity is less than 5000 ppm,[36 , 37 ] even though the individual values of t1 and t2 are affected by impurities at concentrations greater than 200 ppm . at paramagnetic impurity contents > 5000 ppm , re - calibration of the correlation between esurf ( nmr ) and emax ( tpd )
, the presence of magnetic susceptibility contrast between the solid and liquid will distort the magnetic field such that the diffusion - sensitive t2 measurement is no longer determined by surface adsorption .
we have shown recently that the correct t2 relaxation time can be extracted in the presence of pore - scale magnetic field inhomogeneities .
neither of these conditions apply to the materials studied here , so the observed t1/t2 ratios can be interpreted directly in terms of a surface adsorption energy .
in this section we establish proportionality between nmr relaxation and adsorption energy for a molecule interacting with the surface of a pore ; the relaxation of the adsorbed molecule will be modulated by the reduction in molecular mobility and the dipole  dipole coupling between protons in the adsorbed molecule and those bound to the surface .
the present analysis develops from the established theory of surface relaxation , which assumes the pore surface contains paramagnetic impurities that act as relaxation sinks and it is dipole - electron coupling that determines t1/t2 . in recent work we demonstrated that the analysis of godefroy et al . for proton
 electron dipole coupling is extensible to proton  proton dipole coupling , albeit with a much weaker coupling constant . in the absence of paramagnetic species , it is the active binding sites on the pore surface that are considered to be the relaxation sinks .
molecules adsorbed onto a surface undergo two - dimensional ( 2d ) translational motion governed by the activation energy em .
the diffusion coefficient associated with this surface motion is : in which dm0 is a temperature independent contribution to surface diffusion , r is the ideal gas constant , and t is the temperature . in liquid - saturated porous materials , adsorbed molecules exchange with molecules not directly interacting with the pore surface .
the effective surface diffusion coefficient deff ( observed ) is modified from dm due to the finite residence time on the surface .
therefore , the effective diffusion coefficient is : in which e = emes is an activation energy for surface diffusion .
the diffusion coefficient has an associated correlation time m that describes the time for diffusion between active binding sites , such that deff(t)/(4m ) , where  is the thickness of the adsorbed surface layer .
thus , the surface correlation time is defined as : there will also be a surface residence time for adsorbed molecules , which is similarly defined as : the observed relaxation times for liquid saturated porous media are : in which 1,2 are the surface relaxivity values governing longitudinal and transverse relaxation processes at the pore surface , respectively , t1,2,bulk are the respective relaxation times of the bulk liquid , and p is the fraction of spins on the pore surface . for the catalyst supports of interest , the pore diameters are on the order of nanometres , such that the observed relaxation times are determined only by surface relaxation .
the ratio of observed ( surface ) relaxation times is defined by the spectral density function j(0 ) as : at high field , this ratio is known to be insensitive to larmor frequency .
based on typical measured values of m , s for liquids on an alumina surface , and given that 0 is large , we assume that ( m/s)(0m)1 . through a series of trivial algebraic steps it can be shown that t1/t2 is proportional to ln(s/m)/ln(0m ) .
theoretical calculations previously reported by mcdonald et al . suggest that t1/t2 is sensitive to small changes in ln(m ) but insensitive to large changes in ln(s/m ) as m , s are not independent variables .
therefore , for the purpose of establishing a relationship between surface relaxation and adsorption energy it follows from equation ( 3 ) and the reduced form of equation ( 6 ) , given 0 is constant , that : we define a dimensionless surface interaction parameter esurf=t2/t1 for the purpose of interpreting surface relaxation as a surface energy . for the materials of interest in this study ,
we are concerned only with the dynamics of the adsorbed liquid that give rise to altered nmr relaxation times .
however , it is important to consider the extent to which other properties of the sample will influence the observed relaxation times , and the t1/t2 ratio .
the most important of these are 1 ) the presence of paramagnetic species on the pore surface , and 2 ) the presence of magnetic susceptibility contrast between the solid and liquid .
first , considering the case of presence of paramagnetic species on the pore surface ; such species will alter the sensitivity of the relaxation mechanism to m and s .
empirical evidence suggests the ratio t1/t2 remains independent of the density of surface paramagnetic species ( such as fe ) when the concentration of the paramagnetic impurity is less than 5000 ppm,[36 , 37 ] even though the individual values of t1 and t2 are affected by impurities at concentrations greater than 200 ppm . at paramagnetic impurity contents > 5000 ppm , re - calibration of the correlation between esurf ( nmr ) and emax ( tpd )
, the presence of magnetic susceptibility contrast between the solid and liquid will distort the magnetic field such that the diffusion - sensitive t2 measurement is no longer determined by surface adsorption .
we have shown recently that the correct t2 relaxation time can be extracted in the presence of pore - scale magnetic field inhomogeneities .
neither of these conditions apply to the materials studied here , so the observed t1/t2 ratios can be interpreted directly in terms of a surface adsorption energy .
zirconia ( zro2 ) , -alumina ( -al2o3 ) , -alumina ( -al2o3 ) and silica ( sio2 ) were supplied by johnson matthey plc ; anatase titania ( tio2-a ) and rutile titania ( tio2-r ) were supplied by evonik - degussa .
all materials were used as - received in the form of 34 mm diameter cylindrical and 45 mm length extrudate pellets .
samples for nmr relaxation and tpd measurements were prepared by impregnating the different porous materials in deionised water for at least 24 h. excess surface water was removed from the pellets by gently contacting them on a pre - soaked filter paper .
the suppliers of the alumina , silica and zirconia samples do not report any measureable paramagnetic impurities .
the average pore diameter determined was as follows : tio2-a ( 23 nm ) ; tio2-r ( 38 nm ) ; -al2o3 ( 15 nm ) ; -al2o3 ( 7 nm ) ; sio2 ( 5 nm ) and zro2 ( 9 nm ) .
nmr experiments were carried out on a bruker biospec ( horizontal bore ) av 85 mhz spectrometer .
samples were prepared by soaking the different solid materials in the liquid for at least 24 h ; the pellets of solid material were dried on a pre - soaked filter paper , in order to remove any excess liquid on the external surface , and finally transferred into 20 mm diameter glass vials , which were then sealed .
a standard t1t2 pulse sequence was used to acquire 2d correlation data . the specific sequence of rf pulses is illustrated in figure 1 .
a 180 rf inversion pulse rotates the spin magnetisation , initially at equilibrium along the z axis , onto the z axis .
longitudinal t1 relaxation then occurs for a time 1 ; during this time the spin system recovers along the z axis and is aligned with the static magnetic field of magnitude b0 . after this recovery time , a 90 rf excitation pulse is applied and rotates the recovered spin magnetisation into the xy plane .
a series of 180 rf refocusing pulses are then applied to generate a train of n spin echoes , each separated in time by te=22 .
the amplitude of each echo is recorded as a single point , and all echo amplitudes are recorded in a single scan ( no chemical resolution ) .
the amplitude of the initial echo is determined by the degree of recovery during 1 ; the envelope of the echo train is described by the transverse t2 relaxation of the spin ensemble . by repeating the experiment for different 1 recovery times , a 2d data matrix is constructed . here ,
16 t1 recovery delays were used , ranging from 1=1 ms to 10 s ; n=1024 spin echoes were acquired in a single shot with an echo time spacing te=0.8 ms .
a recycle delay of 10 s was included between each scan to ensure maximum signal was obtained at all times .
the total experiment duration was 1.5 h , and included 32 repeat scans to accommodate the rf phase cycle and provide signal averaging to improve the signal - to - noise ratio ( snr ) of the data .
the t1t2 pulse sequence , showing the rf pulses as thin ( 90 ) and thick ( 180 ) vertical bars .
the amplitude of each of the n spin echoes is recorded as a single datum .
the 2d data are inverted numerically to form a 2d distribution of t2 correlated against t1 .
the nmr data are described by the first kind fredholm integral equation : in which the nmr signal is b ,  is the experimental error ( noise ) , f(t1,t2 ) is the required 2d correlation , and the kernel function k(1,t1,nte , t2 ) represents the expected form of the data so that : in which the first exponent describes the t1 relaxation and the second exponent describes the t2 relaxation . as the two exponents do not share a common time base , the kernel function is separable , that is , the expected behaviour of the signal amplitude for a given set of ( 1,t1 ) can be determined separately from the expected behaviour of the signal amplitude for a given set of ( nte , t2 ) .
this separation allows us to solve the fredholm integral equation in ( 9 ) efficiently in vector 
the result is biased to be positive , smooth , and bounded by pre - determined limits in t1 and t2 .
a stable distribution , obtained in the presence of noise , is found using tikhonov regularisation with the smoothing parameter chosen using the generalised cross validation method ; a review of data inversion techniques is presented elsewhere .
the inversion is highly susceptible to noise fluctuations and the degree of smoothing increases as the snr decreases .
variations in distribution shape of less than an order of magnitude on the relaxation time axes are therefore considered to be determined by the signal quality and not the sample .
tpd experiments were performed on a catlab - pcs ( hiden analytical ) , comprising a microreactor module with integrated mass spectrometer .
pellets saturated with water were placed into the glass microreactor in a high purity helium flow at a constant rate of 40 ml min and left for one hour at 45 c until all the physisorbed water was removed .
tpd curves of water ( m=18 amu ) were recorded in the range 451000 c with a heating rate of 20 c min .
each tpd measurement lasted 3 h. the amount of water desorbed from each porous material during the tpd experiments was calculated by using a calibration factor , which was obtained by measuring the area under the tpd profile corresponding to a known amount of desorbed water from a calibration sample .
the analysis of the tpd curves was carried out according to the condensation approximation method reported by barrie in order to obtain the distribution function of the activation energy of desorption .
having derived equation ( 7 ) , we now demonstrate that t1/t2 for the systems studied can indeed be related to the strongest relaxation sinks present ( i.e. , the strongest adsorption sites present)the strength of the strongest adsorption sites being determined by tpd analysis .
we determine the strength of water adsorption in six porous oxides ( extruded pellets ) used as supports and catalysts : tio2-a ( anatase ) , tio2-r ( rutile ) ( supplied by evonik - degussa ) , zro2 , -al2o3 , -al2o3 and sio2 . we chose water as the adsorbate because it is central to aqueous - phase heterogeneous catalytic processes and , in the context of this work , has the additional advantages of not decomposing during thermally driven desorption .
the tpd energy distributions for water desorbing from the porous oxides are shown in figure 2 . due to the uncertainty associated with the tail of the tpd curves ,
we determine the maximum activation energy of desorption ( emax ) where the integral area of the desorption curves is 95 % of the total , see figure 2 .
energy distribution functions obtained from tpd analysis of water in a ) tio2-a , b ) tio2-r , c ) -al2o3 , d ) sio2 , e ) -al2o3 and f ) zro2 .
the vertical dashed lines indicate the activation energy of desorption ( emax ) at which 95 % of species have desorbed .
it is clearly seen that the shape of the tpd curves differs widely between the oxides studied . ignoring any detailed structure in these plots
, we see that the tio2 materials ( figure 2a , b ) are associated with significantly lower maximum desorption energies than the other materials studied .
based on the values of emax determined from the tpd the strength of the interaction of water with the porous solid decreases in the order : the t1t2 correlations for water in the porous oxides are shown in figure 3 . in each case , a single relaxation time component is observed ; the relaxation times are dominated by surface adsorbed species due to the small pore size and high surface area in these oxides .
the ratio t1/t2 is obtained from the logarithmic mean of the individual t1,2 dimensions , which corresponds almost exactly with the maximum intensity of the 2d peak in these mono - modal distributions . in these porous materials with narrow , mono - modal pore size distributions , details of the peak shape are determined predominantly by the raw data quality ( degree of smoothing on inversion ) and are not considered representative of physical sample properties .
the 2d correlations provide a straightforward visual comparison and are required to interpret relaxation results obtained from complicated systems with multiple liquids or diffusive exchange.[24 , 29 ] t1t2 correlation plots for water in : a ) tio2-a , b ) tio2-r , c ) -al2o3 , d ) sio2 , e ) -al2o3 and f ) zro2 .
the solid diagonal line indicates t1=t2 ; the dashed line indicates t1/t2 at the maximum of the peak .
a comparison between the nmr relaxation analysis and tpd is achieved by converting the t1/t2 ratio into an effective surface interaction parameter , esurf , using equation ( 7 ) .
the comparison between emax ( tpd ) and esurf ( nmr ) is given in figure 4 , in which an excellent correlation is obtained between the two measurements .
notably , the tio2 samples have a markedly lower esurf value compared to the other oxides , consistent with their tpd spectra being significantly different to those of the other oxides . elsewhere ,
a comparison between zro2 and tio2 surfaces has been carried out by ignatchenko and co - workers , who concluded , through isotopic exchange experiments and density functional theory ( dft ) calculations , that zro2 has a greater affinity for water than tio2 .
our results agree with these previous findings , and the high value of esurf obtained for water on zro2 correlates with the known strong affinity of water for this oxide surface .
values of the isosteric heat of adsorption for water on sio2[47 , 48 ] have been reported to be comparable to those measured over -al2o3 .
again , this is in line with the finding of this work , which shows very similar esurf values for sio2 and -al2o3 , with sio2 showing a slightly higher value .
overall , our results show that esurf for a given liquid adsorbed within a range of chemically similar materials does correlate directly with the maximum activation energy of desorption obtained by tpd analysis , as shown in figure 4 .
solid systems under investigation can be obtained , the esurf values can be converted to an estimate of the absolute maximum desorption energies characteristic of that system .
comparison of esurf ( nmr ) against emax ( tpd ) for water on oxides .
symbols correspond to ( left to right ) tio2-a , tio2-r , -al2o3 , sio2 , -al2o3 and zro2 .
error bars represent the uncertainty in determining the peak position in the t1t2 distributions ( esurf ) and 5 % integral area under the tpd curves ( emax ) .
the adsorption site densities , calculated from the area under the tpd curve according to the procedure described in the experimental section , were : 0.11 ( tio2-a and tio2-r ) , 0.84 ( -al2o3 ) , 1.38 ( sio2 ) , 2.95 ( -al2o3 ) , and 0.60 mmol m ( zro2 ) .
no strong correlation is found between these values and the t1/t2 ratios , as expected from the earlier discussion[29 , 30 ] regarding the relatively low sensitivity of t1/t2 to site density compared to the changing nature of interactions between liquids and surfaces of different chemistry .
zirconia ( zro2 ) , -alumina ( -al2o3 ) , -alumina ( -al2o3 ) and silica ( sio2 ) were supplied by johnson matthey plc ; anatase titania ( tio2-a ) and rutile titania ( tio2-r ) were supplied by evonik - degussa .
all materials were used as - received in the form of 34 mm diameter cylindrical and 45 mm length extrudate pellets .
samples for nmr relaxation and tpd measurements were prepared by impregnating the different porous materials in deionised water for at least 24 h. excess surface water was removed from the pellets by gently contacting them on a pre - soaked filter paper .
the suppliers of the alumina , silica and zirconia samples do not report any measureable paramagnetic impurities .
the average pore diameter determined was as follows : tio2-a ( 23 nm ) ; tio2-r ( 38 nm ) ; -al2o3 ( 15 nm ) ; -al2o3 ( 7 nm ) ; sio2 ( 5 nm ) and zro2 ( 9 nm ) .
nmr experiments were carried out on a bruker biospec ( horizontal bore ) av 85 mhz spectrometer .
samples were prepared by soaking the different solid materials in the liquid for at least 24 h ; the pellets of solid material were dried on a pre - soaked filter paper , in order to remove any excess liquid on the external surface , and finally transferred into 20 mm diameter glass vials , which were then sealed . a standard t1t2 pulse sequence was used to acquire 2d correlation data .
a 180 rf inversion pulse rotates the spin magnetisation , initially at equilibrium along the z axis , onto the z axis .
longitudinal t1 relaxation then occurs for a time 1 ; during this time the spin system recovers along the z axis and is aligned with the static magnetic field of magnitude b0 . after this recovery time , a 90 rf excitation pulse is applied and rotates the recovered spin magnetisation into the xy plane . a series of 180 rf refocusing pulses
are then applied to generate a train of n spin echoes , each separated in time by te=22 .
the amplitude of each echo is recorded as a single point , and all echo amplitudes are recorded in a single scan ( no chemical resolution ) .
the amplitude of the initial echo is determined by the degree of recovery during 1 ; the envelope of the echo train is described by the transverse t2 relaxation of the spin ensemble . by repeating the experiment for different 1 recovery times , a 2d data matrix is constructed . here , 16 t1 recovery delays were used , ranging from 1=1 ms to 10 s ; n=1024 spin echoes were acquired in a single shot with an echo time spacing te=0.8 ms .
a recycle delay of 10 s was included between each scan to ensure maximum signal was obtained at all times .
the total experiment duration was 1.5 h , and included 32 repeat scans to accommodate the rf phase cycle and provide signal averaging to improve the signal - to - noise ratio ( snr ) of the data .
the t1t2 pulse sequence , showing the rf pulses as thin ( 90 ) and thick ( 180 ) vertical bars .
the amplitude of each of the n spin echoes is recorded as a single datum .
the 2d data are inverted numerically to form a 2d distribution of t2 correlated against t1 .
the nmr data are described by the first kind fredholm integral equation : in which the nmr signal is b ,  is the experimental error ( noise ) , f(t1,t2 ) is the required 2d correlation , and the kernel function k(1,t1,nte , t2 ) represents the expected form of the data so that : in which the first exponent describes the t1 relaxation and the second exponent describes the t2 relaxation . as the two exponents do not share a common time base , the kernel function is separable , that is , the expected behaviour of the signal amplitude for a given set of ( 1,t1 ) can be determined separately from the expected behaviour of the signal amplitude for a given set of ( nte , t2 ) .
this separation allows us to solve the fredholm integral equation in ( 9 ) efficiently in vector 
the result is biased to be positive , smooth , and bounded by pre - determined limits in t1 and t2 .
a stable distribution , obtained in the presence of noise , is found using tikhonov regularisation with the smoothing parameter chosen using the generalised cross validation method ; a review of data inversion techniques is presented elsewhere .
the inversion is highly susceptible to noise fluctuations and the degree of smoothing increases as the snr decreases .
variations in distribution shape of less than an order of magnitude on the relaxation time axes are therefore considered to be determined by the signal quality and not the sample .
tpd experiments were performed on a catlab - pcs ( hiden analytical ) , comprising a microreactor module with integrated mass spectrometer .
pellets saturated with water were placed into the glass microreactor in a high purity helium flow at a constant rate of 40 ml min and left for one hour at 45 c until all the physisorbed water was removed .
tpd curves of water ( m=18 amu ) were recorded in the range 451000 c with a heating rate of 20 c min .
each tpd measurement lasted 3 h. the amount of water desorbed from each porous material during the tpd experiments was calculated by using a calibration factor , which was obtained by measuring the area under the tpd profile corresponding to a known amount of desorbed water from a calibration sample .
the analysis of the tpd curves was carried out according to the condensation approximation method reported by barrie in order to obtain the distribution function of the activation energy of desorption .
having derived equation ( 7 ) , we now demonstrate that t1/t2 for the systems studied can indeed be related to the strongest relaxation sinks present ( i.e. , the strongest adsorption sites present)the strength of the strongest adsorption sites being determined by tpd analysis .
we determine the strength of water adsorption in six porous oxides ( extruded pellets ) used as supports and catalysts : tio2-a ( anatase ) , tio2-r ( rutile ) ( supplied by evonik - degussa ) , zro2 , -al2o3 , -al2o3 and sio2 .
we chose water as the adsorbate because it is central to aqueous - phase heterogeneous catalytic processes and , in the context of this work , has the additional advantages of not decomposing during thermally driven desorption .
the tpd energy distributions for water desorbing from the porous oxides are shown in figure 2 . due to the uncertainty associated with the tail of the tpd curves ,
we determine the maximum activation energy of desorption ( emax ) where the integral area of the desorption curves is 95 % of the total , see figure 2 .
energy distribution functions obtained from tpd analysis of water in a ) tio2-a , b ) tio2-r , c ) -al2o3 , d ) sio2 , e ) -al2o3 and f ) zro2 .
the vertical dashed lines indicate the activation energy of desorption ( emax ) at which 95 % of species have desorbed .
it is clearly seen that the shape of the tpd curves differs widely between the oxides studied . ignoring any detailed structure in these plots
, we see that the tio2 materials ( figure 2a , b ) are associated with significantly lower maximum desorption energies than the other materials studied . based on the values of emax determined from the tpd the strength of the interaction of water with the porous solid decreases in the order : the t1t2 correlations for water in the porous oxides are shown in figure 3 . in each case , a single relaxation time component is observed ; the relaxation times are dominated by surface adsorbed species due to the small pore size and high surface area in these oxides .
the ratio t1/t2 is obtained from the logarithmic mean of the individual t1,2 dimensions , which corresponds almost exactly with the maximum intensity of the 2d peak in these mono - modal distributions . in these porous materials with narrow , mono - modal pore size distributions , details of the peak shape are determined predominantly by the raw data quality ( degree of smoothing on inversion ) and are not considered representative of physical sample properties .
the 2d correlations provide a straightforward visual comparison and are required to interpret relaxation results obtained from complicated systems with multiple liquids or diffusive exchange.[24 , 29 ] t1t2 correlation plots for water in : a ) tio2-a , b ) tio2-r , c ) -al2o3 , d ) sio2 , e ) -al2o3 and f ) zro2 .
the solid diagonal line indicates t1=t2 ; the dashed line indicates t1/t2 at the maximum of the peak .
a comparison between the nmr relaxation analysis and tpd is achieved by converting the t1/t2 ratio into an effective surface interaction parameter , esurf , using equation ( 7 ) .
the comparison between emax ( tpd ) and esurf ( nmr ) is given in figure 4 , in which an excellent correlation is obtained between the two measurements .
notably , the tio2 samples have a markedly lower esurf value compared to the other oxides , consistent with their tpd spectra being significantly different to those of the other oxides . elsewhere ,
a comparison between zro2 and tio2 surfaces has been carried out by ignatchenko and co - workers , who concluded , through isotopic exchange experiments and density functional theory ( dft ) calculations , that zro2 has a greater affinity for water than tio2 .
our results agree with these previous findings , and the high value of esurf obtained for water on zro2 correlates with the known strong affinity of water for this oxide surface . measurements of isosteric heat of adsorption have also been used to study porous oxides .
values of the isosteric heat of adsorption for water on sio2[47 , 48 ] have been reported to be comparable to those measured over -al2o3 .
again , this is in line with the finding of this work , which shows very similar esurf values for sio2 and -al2o3 , with sio2 showing a slightly higher value .
overall , our results show that esurf for a given liquid adsorbed within a range of chemically similar materials does correlate directly with the maximum activation energy of desorption obtained by tpd analysis , as shown in figure 4 .
solid systems under investigation can be obtained , the esurf values can be converted to an estimate of the absolute maximum desorption energies characteristic of that system .
comparison of esurf ( nmr ) against emax ( tpd ) for water on oxides .
symbols correspond to ( left to right ) tio2-a , tio2-r , -al2o3 , sio2 , -al2o3 and zro2 .
error bars represent the uncertainty in determining the peak position in the t1t2 distributions ( esurf ) and 5 % integral area under the tpd curves ( emax ) .
the adsorption site densities , calculated from the area under the tpd curve according to the procedure described in the experimental section , were : 0.11 ( tio2-a and tio2-r ) , 0.84 ( -al2o3 ) , 1.38 ( sio2 ) , 2.95 ( -al2o3 ) , and 0.60 mmol m ( zro2 ) .
no strong correlation is found between these values and the t1/t2 ratios , as expected from the earlier discussion[29 , 30 ] regarding the relatively low sensitivity of t1/t2 to site density compared to the changing nature of interactions between liquids and surfaces of different chemistry .
in this work we have presented a theoretical analysis to show that the ratio of nmr relaxation times , t1/t2 , can be related directly to the adsorbate 
adsorbent interaction energy characterising adsorption in a porous material , and we introduce a quantitative metric esurf ( based on the relaxation time ratio ) characterising the strength of this surface interaction .
we confirm that the t1/t2 ratio is insensitive to pore geometry , allowing adsorption of liquids to be compared between different materials .
given that the relaxation time characteristics of the nuclear spin system will be dominated by the strongest adsorption sites , the hypothesis that t2/t1 should correlate directly with the maximum activation energy of desorption observed in a temperature - programmed desorption experiment was tested ; the data were found to support this hypothesis . given this result
, we now have a physical interpretation of the nmr relaxation data in terms of a surface activation energy of desorption .
it follows from the analysis that for a series of materials for which at least two calibration values are obtainable by tpd , the esurf parameter yields a direct estimate of the maximum activation energy of desorption from the surface .
although the method has been demonstrated for water , the nmr relaxation analysis can be extended to probe non - invasively any adsorbed species , notably organic molecules .
further , given that nmr relaxation is readily extended to measurements at elevated conditions of temperature and pressure , these results suggest that the technique might be considered a useful addition to the toolkit of techniques employed to characterise adsorption in porous media used for applications in catalysis and separations processes . beyond these applications , demonstrating that t1/t2 , when expressed as esurf , is a quantitative measure of interaction strength has implications for determining surface wettability in oil - field reservoir rocks , an important parameter for predicting oil recovery .
an improved correlation between relaxation time and surface adsorption energy could be obtained by measuring t1/t2 as a function of temperature or field strength , although such procedures for providing measurements of surface interaction strength would be considerably more time consuming .
nmr relaxation has been validated as a robust and versatile tool for quantitatively comparing liquid 
solid interactions in porous materials . in principle , the method is applicable to the characterisation of both liquid and gas phase adsorption , although signal - to - noise considerations may limit its applicability to gas  solid interactions .
this limitation and the consideration of chemically different surfaces ( e.g. , carbon - based ) , or materials with much smaller pore size ( e.g. , zeolites ) are the subject of further consideration .
 nuclear magnetic resonance ( nmr ) relaxation times are shown to provide a unique probe of adsorbate  adsorbent interactions in liquid - saturated porous materials . 
 a short theoretical analysis is presented , which shows that the ratio of the longitudinal to transverse relaxation times ( t1/t2 ) is related to an adsorbate  adsorbent interaction energy , and we introduce a quantitative metric esurf ( based on the relaxation time ratio ) characterising the strength of this surface interaction . 
 we then consider the interaction of water with a range of oxide surfaces ( tio2 anatase , tio2 rutile , -al2o3 , sio2 , -al2o3 and zro2 ) and show that esurf correlates with the strongest adsorption sites present , as determined by temperature programmed desorption ( tpd ) . 
 thus we demonstrate that nmr relaxation measurements have a direct physical interpretation in terms of the characterisation of activation energy of desorption from the surface . 
 further , for a series of chemically similar solid materials , in this case a range of oxide materials , for which at least two calibration values are obtainable by tpd , the esurf parameter yields a direct estimate of the maximum activation energy of desorption from the surface . 
 the results suggest that t1/t2 measurements may become a useful addition to the methods available to characterise liquid - phase adsorption in porous materials . 
 the particular motivation for this work is to characterise adsorbate  
 surface interactions in liquid - phase catalysis . 