Inﬂuence of Topology and Brønsted Acid Site Presence on Methanol Di ﬀ usion in Zeolites Beta and MFI

: Detailed insight into molecular di ﬀ usion in zeolite frameworks is crucial for the analysis of the factors governing their catalytic performance in methanol-to-hydrocarbons (MTH) reactions. In this work


Introduction
The zeolite-catalyzed conversion of methanol to hydrocarbons (MTH) has received considerable attention over the last four decades as a reliable source of light olefins with diverse applications in the chemical industry [1].These aluminosilicate crystals combine high surface area, great thermal stability, and shape/size selectivity, performing as excellent catalysts in the conversion process [2][3][4][5].Despite several proposed reaction mechanisms [6][7][8][9][10][11], a clear understanding of how structural features affect the catalytic performance is still missing.Diffusion studies are critical to improve our insight, since the movement of methanol through the pore system is affected by the interaction with the active sites, variations in temperature and modifications in the zeolite topology and micropore dimension, which are all factors that also control the efficiency of zeolites as catalysts in the MTH process.
Earlier experiments have shown that the MTH process over acidic zeolites MFI and Beta have similar conversion efficiencies but differ in molecular shape and size selectivity [12].The role of Brønsted acid sites in the analysis of multicomponent diffusion studies, such as CH 3 OH/(CH 3 ) 2 O [13] or CH 3 OH/H 2 O [14], is ignored, with scarce knowledge on the most favorable adsorption conformation for the methylation of zeolites.The majority of the reported mechanisms propose the C-O cleavage in the methanol molecule at the Brønsted acid sites of the zeolite framework [15].The complexity of the suggested mechanisms further expands with recent studies demonstrating the possible formation of a carbene-like intermediate when the framework methoxy species reacts with hydrocarbon-pool compounds [16,17].The gap in our knowledge of the effect of the Brønsted acid sites on the mobility of methanol within the zeolite frameworks can be resolved by comparing the diffusivities of methanol in all-silica and acidic frameworks to quantify the influence of the Brønsted acid sites on a crucial part of the MTH process.Diffusion studies serve as a useful tool in determining and understanding the methanol-zeolite interactions, which ultimately enhances our existing knowledge of the catalytic activity, selectivity, and lifespan of the zeolites [18,19].
It is well-established that the MTH process in zeolites is a temperature-sensitive reaction.However, considerable uncertainty remains on the optimum range of temperature within which the reaction occurs.Neutron spectroscopy experiments of methanol in H-ZSM-5 [20] have shown that methylation can occur under ambient conditions, while other experimental studies have indicated that methylation occurs at high temperatures between 473 and 523 K [21].It is worth noting that methoxy species formed at temperatures below 373 K have low stability and are susceptible to post-dehydration hydrolysis, which results in the reformation of methanol [22].Some proposed kinetic models have shown the deactivation of ZSM-5 in the MTH process by coke formation observed at a relatively higher temperature range of 573-723 K.At temperatures above 723 K, the irreversible dealumination of the zeolite catalyst is observed in the production of light olefins from chloromethane [23].Dimethyl ether formation, which is presumed to be an intermediate in the MTH process, shows high yield with a reduction in selectivity as the temperature increases [24].From this information, a tentative temperature range for the MTH process can be set with lower and upper limits of 373-573 K.In this work, we present the diffusion analysis of methanol with the highest temperature set at 473 K, which should provide sufficient information on the initial stages of the dehydration process as a function of temperature.
Reports centered around the topological features and their relevance in catalysis have shown that zeolites with relatively small pores, specifically 10-membered channels, have better selectivity toward light olefins [25].Relatively larger pore zeolites, such as MOR and Beta, are popular in reactions involving the aromatization of olefins, although they tend to deactivate easily and they also require comparatively higher concentrations of substrates [26][27][28].Hence, the pore dimension and framework features are crucial for the conversion process and product selectivity.In order to analyze the correlation between framework topology and catalytic activity in the MTH reaction, other factors that govern the performance of the catalysts have to remain unchanged, e.g., the acid site distribution and methanol concentration [29,30].Only then can one distinctly recognize the effect of different pore dimensions on the catalytic activity [31].In this regard, computer simulations are a natural choice to perform this analysis, allowing us to systematically tune the features of the model to subsequently observe the effects of the parameters of interest.
In the present work, we report the investigation of methanol diffusion in all-silica and acidic zeolites MFI and Beta, using molecular dynamics (MD) simulations to analyze the effect of topology, temperature, and molecular concentration on the mobility of methanol and its interaction with the Brønsted acid sites.We inspect experimentally relevant parameters, such as diffusion barriers and mean square displacement of the adsorbates, over the range of temperature 373-473 K and at methanol concentrations between 2 and 5 molecules per unit cell.The calculated results demonstrate the relevance of methanol diffusion studies as they inform and validate some proposed reaction mechanisms involved in zeolite methylation, which is a key aspect of the MTH reactions.

Results and Discussion
The comparison between the acidic and all-silica structures allows us to assess the effect of the Brønsted acid sites on the diffusion of methanol through the micropore system of the zeolites.As a general observation, diffusion is slower in the presence of acid sites, owing to H-bonding interactions between methanol OH groups and the acidic protons; this translational motion slows down for both acidic and all-silica zeolites when the molecular loading is increased, as shown in the mean square displacement (MSD) plots of Figure 1.The slowest diffusion of methanol in zeolite Beta is observed for a molecular loading of 4 molecules per unit cell (mpuc) in the acidic structure at a temperature of 373 K, with a diffusion coefficient of 89.9 × 10 −10 m 2 s −1 .This value steadily increases to 143.0 × 10 −10 m 2 s −1 at a temperature of 473 K, yielding an activation energy for the translational diffusion of 6.8 kJ•mol −1 (see Table 1).When the concentration of methanol is halved to 2 mpuc, the coefficients increase by factors of between 1.2 and 1.4, with the energy barrier increasing by 2 kJ•mol −1 .The diffusion in the all-silica structure increases by a factor of 1.2 to 1.3 compared to the acidic framework, with a marginal decrease in the activation energy, down to 8.5 and 5.3 kJ•mol −1 for concentrations of 2 and 4 mpuc, respectively.Zeolite MFI shows a considerable reduction in translational diffusion compared to zeolite Beta, with diffusion coefficients lower by factors between 20 and 40 for the acidic structure and between 10 and 30 for the all-silica framework.This observation is a direct consequence of the reduction of the pore size in zeolite MFI, with a dimension of 10 MR compared to zeolite Beta, which features 12 MR pore windows.Thus, the methanol agglomeration is more pronounced in MFI, which affects more markedly the H-bonding interactions between methanol and the acidic protons.The effect of the micropore dimension can be observed clearly in the trajectories of randomly selected molecules at 473 K, as shown in Figure 2. In zeolite Beta, the molecule covers practically all the micropore volumes, irrespective of whether the zeolite is in its acidic or all-silica form.In the 10 MR pores of zeolite MFI, the methanol movement is highly reduced where Brønsted acid sites are present, showing less hindered diffusion in the all-silica structure.It is worth noting that the localized motion of the methanol molecules in the spherical intersections of the straight and sinusoidal channels, even in all-silica zeolite MFI, is consistent with recent quasielastic neutron scattering (QENS) observations [18].We observe that for a concentration of 3 mpuc and at a temperature of 373 K, the diffusion coefficient in acidic MFI can be as low as 3.4 × 10 −10 m 2 s −1 , increasing to 9.7 × 10 −10 m 2 s −1 for the highest probed temperature of 473 K similar to other diffusion studies for methanol in zeolite HY [32].This variation yields an activation energy of 15.4 kJ•mol −1 .When the molecular loading is increased to 5 mpuc in acidic MFI, the coefficients decrease by factors of between 1.5 and 1.8 as a direct consequence of higher molecular agglomeration, although the activation energy of diffusion is reduced from 15.4 kJ•mol −1 for 3 mpuc to 13.0 kJ•mol −1 for 5 mpuc.This behavior is in contrast to the all-silica MFI framework, where the activation energy of diffusion increases with the concentration of methanol, from 4.3 to 9.5 kJ•mol −1 , as shown in Table 1.It has been reported previously that when the molecular concentration increases, re-adsorption events at the acid sites play a more prominent role than molecular agglomeration in the variation of the diffusion barrier in acidic zeolite; it is in the all-silica structure where molecular agglomeration has a more determining impact [33].These trends are explained by first considering that methanol adsorbs preferentially at the acid sites in the acidic zeolites Beta and MFI, where the molecule resides for most of the simulation time, as shown in the MD trajectories of Figure 2. The micropore volume covered by the diffusion of each molecule is hence much lower compared to the all-silica structure, and therefore an increment in molecular loading will not notably affect a translation that is already strongly reduced by the H-bonding interaction with the acidic protons.This effect does not occur in the all-silica MFI framework, where the methanol molecules roam more freely in the supercell, but they are primarily located within the intersections of the channels without a preferential adsorption site; consequently, additional molecules will more clearly disrupt their diffusion, which explains why the diffusion barrier increases with concentration in the all-silica zeolite MFI.We also have to consider the re-adsorption of methanol at the acid sites [34].At the lowest concentration of 3 mpuc, more acidic protons are freely available for a methanol molecule to re-adsorb when the H-bonding interaction with a current acid site is broken; therefore, re-adsorption events will reduce the diffusion and yield a higher energy barrier.When the number of freely available sites is decreased as a consequence of the increase in molecular concentration, re-adsorption is less probable, hence having a less significant impact on the energy barrier.Accordingly, the activation energy of diffusion decreases in acidic zeolites Beta and MFI as the concentration of methanol increases.This trend should hold as long as the acid sites are not saturated by methanol molecules, in which case agglomeration should start to play a more important role [33].It is apparent that the larger pores of zeolite Beta, which exert less steric hindrance on the relatively small methanol molecules, diminish the effects of molecular agglomeration in its all-silica framework, since its diffusion barrier decreases from 8.5 to 5.3 kJ•mol −1 with the increase of the molecular concentration.
The importance of the pore dimension in the catalytic activity of zeolites cannot be over-stated, as product selectivity is highly dependent on the size and arrangement of the microporous system.Zeolite MFI is distinguished by 10 MR channels, and it is known to be more selective toward C 5 -C 6 carbon chains, with reduced preference for aromatics in the MTH process.Observed coefficients for the translational diffusion of methanol in ZSM-5 from both nuclear magnetic resonance [35] and quasi-elastic neutron scattering methods [36] have yielded values within the range of 2.0-3.5 × 10 −10 m 2 s −1 at 300 and 335 K respectively, which are in agreement with the values calculated in this work for similar temperatures (see Table 1).In contrast, high methanol mobility is observed for zeolite Beta, which presents a larger pore diameter of approximately 6.7 Å [37].Zeolite Beta generally shows higher diffusion coefficients for methanol with selectivity toward aromatics as the intra-crystalline porosity significantly increases the access to acid sites, hence improving the conversion rate and selectivity for aromatic compounds [38,39].
From the analysis of the MSD plots, we can conclude that methanol has a higher protonation probability in zeolite MFI compared to zeolite Beta.This analysis can be expanded further by examining the radial distribution function (RDF) and their integrals for the atomic pair formed by the acidic proton and the oxygen atom of the methanol molecule (*O m , H a ), as shown in Figure 3.The integration of the RDFs adds information with regard to the number of atoms within the coordination sphere.In Beta zeolite, a narrow and intense peak is observed at an *O m -H a distance of 2.6 Å, whose intensity is reduced as the temperature increases.Minor peaks can also be detected within the region from 4.0 to 6.0 Å.In comparison, zeolite MFI shows a similar profile, but the intensity of its peak around 2.5 Å is higher.In fact, although this peak also diminishes as the temperature increases, it is as intense at 473 K in zeolite MFI as it is at 373 K in zeolite Beta, which highlights the higher protonation probability of methanol in the 10 MR structure.The molecular confinement in the micropore system plays a critical role in strengthening the H-bonding interaction between the catalyst and the adsorbed molecule [26].This can be corroborated further by considering the average residence time τ of methanol molecules within a distance of 3.5 Å from the acidic proton (plots are available in the Supplementary Data Figure S1) (see Table 2).In zeolite Beta, τ shows a maximum value of about 3 ps for a concentration of 2 mpuc and a temperature of 373 K; the residence time decreases to values of approximately 2 ps at the temperatures 423 and 473 K.As expected, in zeolite MFI, methanol spends almost twice as much time engaged in H-bonding interactions with the acidic protons when compared to zeolite Beta, with τ values between 5 and 6 ps for a loading of 3 mpuc at 373 K, becoming steadily smaller for higher temperatures, down to a value of 2-3 ps for 473 K.However, although a higher retention time is a factor in defining the catalytic performance in the MTH reaction, it also increases the probability of secondary reactions, for instance the formation of coke, which is detrimental to the lifetime of the catalyst [40,41].Some experimental studies have suggested that framework methoxy species are the starting point for the C-C bond formation in MTH reactions, giving much relevance to the methanol-acidic proton interaction [22,42].Many variations in the proposed mechanism for MTH stem from the different orientations that the methanol molecule can assume during adsorption and protonation [43].Once the methanol is H-bonded to the acidic proton, its methyl group can either point to the centre of the pore (end-on configuration) or lean toward the framework O atoms in the vicinity of the acid site (side-on configuration).Density functional theory (DFT) studies on the different methanol adsorption complexes formed at Brønsted acid sites have shown that methanol is adsorbed at the acid site preferably through the end-on adsorption mode, which shows higher thermochemical stability.However, the side-on adsorption configuration directly precedes the dehydration step leading to the cleavage of the C-O bond in methanol, followed by the methoxylation of the zeolite [21].The RDFs of four atomic pairs are examined in order to elucidate the predominant methanol adsorption complex in both zeolites Beta and MFI, as seen in Figure 4.In the case of zeolite Beta, the atomic pair formed by *O m from the hydroxyl OH end of methanol and the acidic proton H a of the zeolite framework (*O m , H a ) shows an intense peak at an *O m -H a distance of 2.6 Å, which we have already discussed earlier in this section.Additionally, three distinct peaks are observed for the atomic pairs composed of the carbon atom *C of the methanol and the non-acidic framework O at (*C, O) distances of 4.6, 6.4, and 8.5 Å.Each peak represents the methanol *C interaction with the three nearest neighboring framework O atoms to the Brønsted acid site.Furthermore, the pair formed by methyl hydrogen *H and the non-acidic oxygen O of the zeolite (*H, O) shows two peaks at the distances 5.4 and 9.1 Å.This is expected, as the spatial arrangement of the three methyl *H atoms will always cause one of them to be closer to the non-acidic framework O relative to the other two, hence producing two peaks.The RDF profiles also show the pair *H m from the hydroxyl of methanol and the non-acidic O of the framework (*H m , O), which yields what appears to be three overlapping peaks of low intensities at distances of 3.6, 5.1, and 7.2 Å.The typical end-on adsorption mode is characterized by a primary H-bond *O m -H a and a secondary H-bond interaction *H m -O while the side-on adsorption complex can be described by a single H-bond *O m -H a [44].The observed peak intensities for the atomic pair (*H m , O) suggest very weak interactions between the *H m atom and the non-acidic O atoms of the zeolite, matching the description of a side-on configuration.In Figure 5, we have created an illustrative image to depict the side-on adsorbed methanol, with average distances for the selected atomic pairs (*O m , H a ), (*C, O), (*H, O), and (*H m , O) based on our RDF analysis.Similar trends are observed in zeolite MFI, although significantly higher signals are recorded for each of the selected atomic pairs at approximately the same distances, indicating higher occurrence of side-on configuration.Thus, our analysis demonstrates higher susceptibility for methoxylation in zeolite MFI compared to zeolite Beta.As a general observation, the peak intensities for atomic pairs (*O m -H a ) and (*C-O) decrease with increasing concentration and temperature in zeolite Beta.Zeolite MFI displays similar trends, except that the peak observed for the (*C-O) atomic pair becomes rather more intense as the methanol concentration increases.The observed peaks for the pairs (*H, O) and (*H m , O) exhibit very little change when either the temperature or loading is varied.These results are confirmed by data from hybrid DFT calculations of active species in ZSM-5 for the MTH process [45], which report on the transient nature of the side-on conformations, leading to methanol dehydration.This provides insight into geometry transitions that take place during proton abstraction and the sequential methoxylation of the zeolite frameworks during the MTH process.
The RDF for the atomic pair formed by the carbon atoms of two different methanol molecules (*C, *C) provides valuable information concerning the molecular clustering inside the micropore system, as shown in Figure 6.We have discussed the impact of the micropore size on the protonation probability of methanol, which also extends to the molecular clustering.The RDF profile of zeolite Beta shows an intense peak at approximately 4.1 Å, with a lesser signal above 8.0 Å.In comparison, zeolite MFI shows what seems to be two overlapping signals of lower intensity within a *C-*C distance region between 4.0 and 7.0 Å.Therefore, we can conclude that the larger pore size in zeolite Beta allows a more intense clustering of the methanol molecules.We have observed that the retention time of methanol by the acidic protons is reduced in both zeolites when the methanol concentration is increased, which means that a higher number of molecules adsorbed in the systems disrupts more easily the interaction with the acid sites due to more frequent methanol-methanol contacts.The methanol-methanol interactions have a significant impact on the adsorption, desorption, and diffusion characteristics of methanol.XRD reports show that at high loadings, methanol clusters of side-on adsorption conformation up to 6 molecules are formed at the intersections of channels in ZSM-5, where the free space diameter (≈10 Å) is approximately twice the channel diameter [46].Larger clusters of methanol molecules in zeolite MFI should be more weakly bound, and their occurrence may be limited to channel intersections, which explains its overlapping RDF peaks and lower integration numbers observed for the atomic pair (*C, *C) as seen in Figure 6.In contrast, zeolite Beta, whose 12 MR channels allow for stronger methanol-methanol interactions, exhibits more intense methanol clustering with relatively higher integration numbers compared to zeolite MFI.These findings may be of further relevance when considering the mixed diffusion of methanol with other polar molecules that are also involved in the MTH, such as dimethyl ether, ethanol, and water.Under industrially significant conditions, water is used as a regulatory tool for methanol diffusivity, based on the relatively higher polarity of water causing an increase in molecular clustering [47].Ultimately, methanol clustering facilitates proton transfer from Brønsted acid sites to create methoxonium ions, due to the increased H-bonding interactions from multiple methanol molecules.This may result in an almost barrier-less reaction step for methanol dehydration in the MTH process.

Models
Zeolite framework types MFI and Beta are selected for this study.Zeolite Beta is characterized by interconnected, 3-dimensional 12-membered ring (12 MR) channels with tetragonal P4 1 22 symmetry, while MFI has 10 MR channels, which are also 3-dimensionally interconnected, and an orthorhombic Pnma symmetry [48,49].Both all-silica and acidic frameworks are employed in this study, which allows us to determine the effect of the acid sites on the diffusion of methanol.The unit cell of zeolite Beta with lattice parameters a = 12.46 Å, b = 12.46 Å, and c = 26.22Å is subjected to an expansion to a 4 × 4 × 2 system, leading to a supercell size of approximately 50 Å.The unit cell of MFI, with lattice parameters a = 20.07Å, b = 19.69Å and c = 13.34Å, is replicated to create a 3 × 3 × 4 supercell of approximately 60 Å.
The micropore system of the MFI framework is characterized by sinusoidal channels, running along the a direction, in addition to straight channels, which are aligned parallel to the b direction, plus the intersections between sinusoidal and straight channels.In the present work, we have replaced a number of Si by Al and placed the acidic protons in such a way that each of these structural features is represented in the model.The Al substitution in the sinusoidal channel takes place at the T7 site, binding the acidic proton to the oxygen O18.The site T7 is also used for the straight channel, but this time binding the proton to the oxygen O17.The intersection region is covered by replacing Si by Al at site T6, with the proton bound to O13.Unlike zeolite MFI, the framework of zeolite Beta has the same straight 12 MR channels along the a and b directions, requiring substitution only at the tetrahedral site T6, as shown in Figure 7.A Si/Al ratio of 31 is employed for both zeolites, with the acid sites having an adequate spatial separation consistent with Dempsey's rule [50].

Interatomic Potentials
The Born model of ionic solids, which uses a series of energy equations to describe the interactions between ions, is employed in the description of the zeolite structures and adsorbates [51].The charges for methanol were obtained by a Mulliken population analysis of ab initio DFT calculations by Blanco and Auerbach [52].These parameters were scaled by Plant and associates [53] to attain a more polarized methyl group corresponding to methanol at acidic sites in zeolite frameworks.Methanol intramolecular harmonic potentials for covalent bonds, angles, and dihedrals were determined using ab initio DFT calculation at the quadratic configuration interaction calculation with single and double excitations (QCISD) level of theory to attain the ground state energy and electron density distribution [52].The data obtained were in good agreement with results from electron diffraction [54] and millimeter-wave spectroscopic analysis of methanol [55].The selected Lennard-Jones potentials listed in Table 3, which describe the methanol-methanol interactions, are empirically fitted values from the consistent valence force field (CVFF) [56] to reproduce diffusion coefficients similar to that of liquid-phase methanol [57].4.
Table 5 shows the Lennard-Jones potentials employed to represent the dispersion forces involved in the methanol-zeolite interactions.These parameters originate from a combination of studies from Shubbin and associates [60] on the diffusion of iso-butanol in H-ZSM-5 and methanol diffusion in all-silica faujasite from Blanco and Auerbach's work [52].The potentials employed in this work have been applied extensively in previous studies of micro-and meso-porous materials and have been shown to accurately describe the adsorption behavior in zeolites [61][62][63][64].

Simulations
All simulations are carried out with the DL_POLY 4.0 code [65].Visualizations and plots for this study are constructed using visualization for electronic and structural analysis (VESTA) and visual molecular dynamics (VMD) [66,67].The MD simulations are initiated by an equilibration, which is conducted first under an NVE ensemble over 2 ns, followed by 2 ns employing an NVT ensemble.Afterward, the production run consists of 12 ns of an NVE ensemble.A time step of 0.5 fs is used in the simulation, saving the atomic coordinates every 1 ps.The simulations are carried out at the temperatures 373, 423, and 473 K.The Berendsen thermostat is applied to regulate the temperature during the equilibration with the NVT ensemble, with a time constant of 1 ps [68].Two different loadings of methanol molecules are employed to analyze the effect of the molecular concentration on the diffusion.At the lower loading, 64 and 93 methanol molecules are adsorbed in the supercells of zeolites Beta and MFI, i.e., 2 and ≈3 mpuc, respectively, while a twofold increase is used for the higher concentrations, i.e., 128 and 186 molecules, resulting in 4 and ≈5 molecules per unit cell, respectively.For each loading, the two zeolites both have the same ratio between the number of molecules and the volume of the supercell.The multiple initial times t 0 method is employed to average the trajectory from a total of 12 ns of production run into 1 ns, shifting the initial time t 0 every 30 ps.This allows a reduction in statistical noise that may be generated in the mean square displacement (MSD) plots.The MSD plots allow the determination of self-diffusion coefficients (D s ) from the Einstein relationship: The calculated activation energies for diffusion are subsequently determined using the Arrhenius equation: where k b is the Boltzmann constant and A is a pre-exponential factor.The value of activation energy E a is derived from the slope of the linear regression obtained from the plot of ln(D s ) vs. 1 T .The good linearity observed in the MSD plots demonstrates that 1 ns is enough to achieve true self-diffusion.Radial distribution function (RDF) plots are used to elucidate the average kinetic relationship between the adsorbate and the zeolite framework, which allows an easy identification of H-bonding.
The contact correlation function C(t) between methanol and the acidic proton is calculated with the following equations: where N is the total number of methanol molecules.The function n i (t) acquires the value of 1 if methanol is within a chosen distance from the acidic proton (in the present work, we use a cut-off of 3.5 Å) or it equals 0 otherwise.Function C(t) provides the probability of finding methanol near the acidic proton at time t if the molecule was in this region at times t 0 and t − ∆t.The residence time t r is calculated from the integration of the function C(t): When C(t) displays an exponential decay, it can be expressed as C(t) = e −t/τ , which allows to equal t r to τ [69].

Conclusion
We have used molecular dynamics simulations for an in-depth investigation of methanol diffusion in zeolites MFI and Beta.Measurements at temperatures of between 373 and 473 K yielded diffusion coefficients within the range of 2.0-10.0× 10 −10 m 2 s −1 for acidic zeolite MFI, with coefficients of two orders of magnitude larger in zeolite Beta, suggesting greater methanol mobility in the larger microporous system of zeolite Beta.The activation energy of translational diffusion was calculated for both zeolites, giving values in acidic zeolite MFI of between 13.0 and 15.4 kJ•mol −1 , compared to zeolite Beta with values of between 6.0 and 8.0 kJ•mol −1 .The analysis of methanol diffusion in the all-silica zeolites indicates that the presence of acid sites results in a significant reduction in methanol diffusion, especially in zeolite MFI.Examination of the MSD and RDF plots showed a higher probability of methanol protonation in zeolite MFI, whereas zeolite Beta exhibited lower peak intensities between the hydroxyl O atom of the methanol molecule and the protons of the Brønsted acid sites.

Figure 2 .
Figure 2. Trajectory of a single methanol molecule at 473 K in (a) H-Beta (b) all-silica Beta c direction at a loading of 4 mpuc, and in (c) H-MFI and (d) all-silica MFI c direction at a loading of 5 mpuc.Element color code: O (red), Si (orange), Al (blue), center of mass for CH 3 OH (green).

Figure 3 .
Figure 3.The radial distribution function and its integration for *O m -H a showing H-bonding between the acid site and methanol as a function of temperature in H-Beta for loadings of (a) 2 and (b) 4 mpuc, and in H-MFI for loadings of (c) 3 and (d) 5 mpuc.

Figure 6 .
Figure 6.The radial distribution function for *C-*C showing the probability for methanol agglomeration as a function of temperature in H-Beta with loadings of (a) 2 and (b) 4 mpuc; in all-silica Beta with loadings of (c) 2 and (d) 4 mpuc; in H-MFI with loadings of (e) 3 and (f) 5 mpuc; and all-silica MFI with loadings of (g) 3 and (h) 5 mpuc.

Figure 7 .
Figure 7. (a) Zeolite Beta with substituted Al at tetrahedral site T6 and proton on O12.Zeolite MFI showing the substituted Al at (b) T7 and T6 with protons on O17 and O13 for straight and the intersection between channels, respectively and (c) site T7 with proton on O18 for the sinusoidal channel.Element color code: O (red), Si (orange), Al (blue), H (gray).
The zeolite atoms Si, Al, and non-acidic O are assigned formal charges +4.00, +3.00, and −2.00 respectively, with fractional charges of −1.426 for the acidic oxygen bridging Si and Al atoms (labelled O a ) and +0.426 for the proton of the Brønsted acid site (labeled as H a ).Buckingham potentials are used to describe the short-range interactions between the atomic pairs Si-O, O-O and Al-O based on parameters reported by Sanders et al. [58] and Catlow et al. [59], while the bridging hydroxyl group is represented by a Morse potential derived by Schroder and collaborators [55].Three-body harmonic potentials are utilized to describe the bends of the angles O-Si-O, O-Al-O, O-Si-O a , and O-Al-O a (labeled as O-T-O).The full set of potentials is listed in Table

2
potentials U θ ijk = k 2 θ ijk − θ 0 Angles k (eV•rad −2 ) θ 0 ( • ) O-T-O b 2.09724 109.47aO a and H a denote the bridging hydroxyl O and the proton of the Brønsted acid sites respectively.b The same three-body potential is used for the atomic triads O-Si-O, O-Al-O, O-Si-O a , and O-Al-O a.

Table 2 .
Residence times τ obtained after the exponential fit of the contact correlation function C(t) at 373, 423 and 473 K.

Table 3 .
Potential parameters describing the intra-and inter-molecular interactions for methanol a .Harmonic dihedral U ϕ ijkl =k [1 + cos βϕ ijkl − α O m and H m denote the O and hydroxyl H atom of methanol, respectively. a

Table 4 .
Potentials describing the interatomic interactions within the zeolite framework and the Brønsted acid sites a .Buckingham Potentials U r ij = A• exp −

Table 5 .
Lennard-Jones potential describing the zeolite-methanol a interaction.
a The asterisk denotes the atoms of the methanol molecule.