Methanol diffusion and dynamics in zeolite H-ZSM-5 probed by quasi-elastic neutron scattering and classical molecular dynamics simulations

Zeolite ZSM-5 is a key catalyst in commercially relevant processes including the widely studied methanol to hydrocarbon reaction, and molecular diffusion in zeolite pores is known to be a crucial factor in controlling catalytic reactions. Here, we present critical analyses of recent quasi-elastic neutron scattering (QENS) data and complementary molecular dynamics (MD) simulations. The QENS experiments show that the nature of methanol diffusion dynamics in ZSM-5 pores is dependent both on the Si/Al ratio (11, 25, 36, 40 and 140), which determines the Brønsted acid site density of the zeolite, and that the nature of the type of motion observed may vary qualitatively over a relatively small temperature range. At 373 K, on increasing the ratio from 11 to 140, the observed mobile methanol fraction increases and the nature of methanol dynamics changes from rotational (in ZSM-5 with Si/Al of 11) to translational diffusion. The latter is either confined localized diffusion within a pore in zeolites with ratios up to 40 or non-localized, longer-range diffusion in zeolite samples with the ratio of 140. The complementary MD simulations conducted over long time scales (1 ns), which are longer than those measured in the present study by QENS (≈1-440 ps), at 373 K predict the occurrence of long-range translational diffusion of methanol in ZSM-5, independent of the Si/Al ratios (15, 47, 95, 191 and siliceous MFI). The rate of diffusion increases slightly by increasing the ratio from 15 to 95 and thereafter does not depend on zeolite composition. Discrepancies in the observed mobile methanol fraction between the MD simulations (100% methanol mobility in ZSM-5 pores across all Si/Al ratios) and QENS experiments (for example, ≈80% immobile methanol in ZSM-5 with Si/Al of 11) are attributed to the differences in time resolutions of the techniques. This perspective provides comprehensive information on the effect of acid site density on methanol dynamics in ZSM-5 pores and highlights the complementarity of QENS and MD, and their advantages and limitations. This article is part of the theme issue 'Exploring the length scales, timescales and chemistry of challenging materials (Part 2)'.

Highly Efficient Proton Conduction in the Metal-Organic Framework Material MFM-300(Cr)·SO4(H3O)2

The development of materials showing rapid proton conduction with a low activation energy and stable performance over a wide temperature range is an important and challenging line of research. Here, we report confinement of sulfuric acid within porous MFM-300(Cr) to give MFM-300(Cr)·SO₄(H₃O)₂, which exhibits a record-low activation energy of 0.04 eV, resulting in stable proton conductivity between 25 and 80 °C of >10^–2 S cm^–1. In situ synchrotron X-ray powder diffraction (SXPD), neutron powder diffraction (NPD), quasielastic neutron scattering (QENS), and molecular dynamics (MD) simulation reveal the pathways of proton transport and the molecular mechanism of proton diffusion within the pores. Confined sulfuric acid species together with adsorbed water molecules play a critical role in promoting the proton transfer through this robust network to afford a material in which proton conductivity is almost temperature-independent.

Methanol dynamics in H-ZSM-5 with Si/Al ratio of 25: a quasi-elastic neutron scattering (QENS) study

Methanol dynamics in zeolite H-ZSM-5 (Si/Al of 25) with a methanol loading of ~ 30 molecules per unit cell has been studied at 298, 323, 348 and 373 K by incoherent quasi-elastic neutron scattering (QENS). The elastic incoherent structure factor (EISF) reveals that the majority of methanol is immobile, in the range between 70 and 80%, depending on the measurement temperature. At 298 K, ≈ 20% methanol is mobile on the instrumental timescale, exhibiting isotropic rotational dynamics with a rotational diffusion coefficient (DR) of 4.75 × 1010 s−1. Upon increasing the measurement temperature from 298 to 323 K, the nature of the methanol dynamics changes from rotational to translational diffusion dynamics. Similar translational diffusion rates are measured at 348 and 373 K, though with a larger mobile fraction as temperature increases. The translational diffusion is characterised as jump diffusion confined to a sphere with a radius close to that of a ZSM-5 channel. The diffusion coefficients may be calculated using either the Volino–Dianoux (VD) model of diffusion confined to a sphere, or the Chudley–Elliot (CE) jump diffusion model. The VD model gives rise to a self-diffusion co-efficient (Ds) of methanol in the range of 7.8–8.4 × 10–10 m2 s−1. The CE model gives a Ds of around 1.2 (± 0.1) × 10–9 m2 s−1 with a jump distance of 2.8 (either + 0.15 or − 0.1) Å and a residence time (τ) of ~ 10.8 (either + 0.1 or − 0.2) ps. A correlation between the present and earlier studies that report methanol dynamics in H-ZSM-5 with Si/Al of 36 is made, suggesting that with increasing Si/Al ratio, the mobile fraction of methanol increases while DR decreases.

Effect of steam de-alumination on the interactions of propene with H-ZSM-5 zeolites

Steam de-alumination is used to prepare a H-ZSM-5 material representative of industrial acid zeolite catalysts. Characterisation shows extensive loss of zeolite acidity but minimal loss of framework crystallinity in the treated material. The material's interaction with propene is probed by means of inelastic and quasielastic neutron scattering, providing information on the reactivity and mobility of the propene respectively. These results are compared to those previously obtained for propene in the untreated zeolite. The steaming treatment resulted in decreased reactivity of the zeolite toward olefin oligomerization, higher temperatures for reaction initiation, and increased mobility of the propene in the zeolite at all temperatures. Analysis of the motions of the propene revealed by QENS shows the mobility to be comparable to those previously reported for propane in similar materials but occurring at slower velocities due to the greater rigidity and polarisation of the propene molecule.

Onset of Propene Oligomerization Reactivity in ZSM-5 Studied by Inelastic Neutron Scattering Spectroscopy

The techniques of quasi-elastic and inelastic neutron scattering (QENS and INS) are applied to investigate the oligomerization of propene over a ZSM-5 zeolite. Investigations are performed at low temperatures, allowing identification of the onset of the oligomerization reaction and observation of the low-energy spectral changes due to intermediate formation that are difficult to observe by optical methods. Oligomerization proceeds via formation of a hydrogen-bonded precursor by an interaction of the propene with an internal acid site followed by protonation and chain growth with protonation being the rate-limiting step. The use of quasi-elastic neutron scattering to observe changes in system mobility with temperature via the elastic window scan technique allows identification of the active temperature range where catalyst activity commences and permits targeting of the more time-consuming INS investigations to conditions of interest. From examination of the product's spectrum, the structure of the resulting oligomer is deduced to be primarily linear.

The Application of Quasi-Elastic Neutron Scattering to Investigate Hydrogen Diffusion in an Iron-Based Fischer–Tropsch Synthesis Catalyst

Iron-based Fischer–Tropsch synthesis (FTS) catalysts evolve in situ on exposure to synthesis gas (CO & H2) forming a mixture of iron oxides, iron carbides and carbonaceous deposits. Recently, the application of inelastic neutron scattering has shown the progressive formation of a hydrocarbonaceous overlayer during this catalyst conditioning period. The evolving nature of the catalyst alters the proportion of phases present within the catalyst, which may influence the transport of hydrogen within the reaction system. Preliminary quasi-elastic neutron scattering (QENS) measurements are used to investigate hydrogen diffusion within an un-promoted iron FTS catalyst that has experienced varying levels of time-on-stream (0, 12 and 24 h) of ambient pressure CO hydrogenation at 623 K. Measurements on the catalyst samples in the absence of hydrogen show the unreacted sample (t = 0 h) to exhibit little increase in motion over the temperature range studied, whereas the t = 12 and 24 h samples exhibit a pronounced change in motion with temperature. The contrast is attributed to the presence of the afore-mentioned hydrocarbonaceous overlayer. Measurements on the samples in the presence of liquid hydrogen show hydrogen diffusional characteristics to be modified as a function of the catalyst conditioning process but, due to the complexity of the evolving catalyst matrix, the hydrogen motion cannot be attributed to a particular phase or component of the catalyst. Problems in the use of hydrogen as a probe molecule in this instance are briefly considered. Coincident neutron diffraction studies undertaken alongside the QENS measurements confirm the transition from hematite pre-catalyst to that of Hägg carbide during the course of extended times-on-stream.

Molecular behaviour of methanol and dimethyl ether in H-ZSM-5 catalysts as a function of Si/Al ratio: a quasielastic neutron scattering study

Qualitative and quantitative differences are found in methanol and dimethyl ether mobility in H-ZSM-5 catalysts of varying Si/Al ratios (Brønsted acid site concentrations) using quasielastic neutron scattering.

Probing the dynamics and structure of confined benzene in MCM-41 based catalysts

Combined MD simulations and QENS experiments on benzene in MCM-41 provide insight into the dynamics and structure of benzene

SAPO-34 Framework Contraction on Adsorption of Ammonia: A Neutron Scattering Study

Neutron scattering data was recorded from SAPO-34 using the OSIRIS instrument before and after repeated ammonia adsorption at pressures up to 8 bar. Coherent scattering from the zeolite framework provides the neutron powder diffraction pattern and gave evidence for anisotropic contraction on ammonia dosing. Incoherent quasielastic scattering from the hydrogen of the ammonia showed that mobile ammonia was present in the framework. The quasielastic data was fitted to a model where the ammonia was confined within the chabazite cage in the c direction of the crystal lattice, with diffusion solely occurring through the perpendicular 8-membered rings. The calculated diffusion constant reached a maximum of 6.3×10^-8  m²  s^-1 at 5 bar.

Modulating proton diffusion and conductivity in metal-organic frameworks by incorporation of accessible free carboxylic acid groups

Three multi-carboxylic acid functionalised ligands have been designed, synthesised and utilised to prepare the new barium-based MOFs, MFM-510, -511, and -512, which show excellent stability to water-vapour. MFM-510 and MFM-511 show moderate proton conductivities (2.1 × 10-5 and 5.1 × 10-5 S cm-1, respectively) at 99% RH and 298 K, attributed to the lack of free protons or hindered proton diffusion within the framework structures. In contrast, MFM-512, which incorporates a pendant carboxylic acid group directed into the pore of the framework, shows a two orders of magnitude enhancement in proton conductivity (2.9 × 10-3 S cm-1). Quasi-elastic neutron scattering (QENS) suggests that the proton dynamics of MFM-512 are mediated by "free diffusion inside a sphere" confirming that incorporation of free carboxylic acid groups within the pores of MOFs is an efficient albeit synthetically challenging strategy to improve proton conductivity.

Molecular behaviour of phenol in zeolite Beta catalysts as a function of acid site presence: a quasielastic neutron scattering and molecular dynamics simulation study

The dynamic behaviour of phenol in zeolite Beta is strongly influenced by the presence of Brønsted acid sites.

Comparing ammonia diffusion in NH3-SCR zeolite catalysts: a quasielastic neutron scattering and molecular dynamics simulation study

The diffusion of ammonia in the small pore zeolite and potential commercial NH3-SCR catalyst levynite (LEV) was measured and compared with its mobility in the chabazite (CHA) topology (more established in NOx abatement catalysis), using quasielastic neutron scattering (QENS) and molecular dynamics (MD) simulations at 273, 323 and 373 K. The QENS experiments suggest that mobility in LEV is dominated by jump diffusion through the 8-ring windows between cages (as previously observed in CHA) which takes place at very similar rates in the two zeolites, yielding similar experimental self-diffusion coefficients (Ds). After confirming that the same characteristic motions are observed between the MD simulations and the QENS experiments on the picosecond scale, the simulations suggest that on the nanoscale, the diffusivity is higher by a factor of ∼2 in the CHA framework than in LEV. This difference between zeolites is primarily explained by the CHA cages having six 8-ring windows in the building unit, compared to only three such windows in the LEV cage building unit, thereby doubling the geometric opportunities to perform jump diffusion between cages (as characterised by the QENS experiments) leading to the corresponding increase in the MD calculated Ds. The techniques illustrate the importance of probing both pico- and nanoscale dynamics when studying intracrystalline diffusion in both NH3-SCR catalyst design, and in porous materials generally, where notable consistencies and differences may be found on either scale.

Ammonia mobility in chabazite: insight into the diffusion component of the NH3-SCR process

The diffusion of ammonia in commercial NH3-SCR catalyst Cu-CHA was measured and compared with H-CHA using quasielastic neutron scattering (QENS) and molecular dynamics (MD) simulations to assess the effect of counterion presence on NH3 mobility in automotive emission control relevant zeolite catalysts. QENS experiments observed jump diffusion with a jump distance of 3 Å, giving similar self-diffusion coefficient measurements for both Cu- and H-CHA samples, in the range of ca. 5-10 × 10(-10) m(2) s(-1) over the measured temperature range. Self-diffusivities calculated by MD were within a factor of 6 of those measured experimentally at each temperature. The activation energies of diffusion were also similar for both studied systems: 3.7 and 4.4 kJ mol(-1) for the H- and Cu-chabazite respectively, suggesting that counterion presence has little impact on ammonia diffusivity on the timescale of the QENS experiment. An explanation is given by the MD simulations, which showed the strong coordination of NH3 with Cu(2+) counterions in the centre of the chabazite cage, shielding other molecules from interaction with the ion, and allowing for intercage diffusion through the 8-ring windows (consistent with the experimentally observed jump length) to carry on unhindered.

Correction: In situ spectroscopic investigations of MoOx/Fe2O3 catalysts for the selective oxidation of methanol

Correction for ‘In situ spectroscopic investigations of MoO_x/Fe₂O₃ catalysts for the selective oxidation of methanol’ by Catherine Brookes et al., Catal. Sci. Technol., 2016, 6, 722–730.

Spectroscopic investigation into the design of solid-acid catalysts for the low temperature dehydration of ethanol

The increased demand for bulk hydrocarbons necessitates research into increasingly sustainable, energy-efficient catalytic processes. Owing to intricately designed structure-property correlations, SAPO-34 has become established as a promising material for the low temperature ethanol dehydration to produce ethylene. However, further optimization of this process requires a precise knowledge of the reaction mechanism at a molecular level. In order to achieve this a range of spectroscopic characterization techniques are required to probe both the interaction with the active site, and also the wider role of the framework. To this end we employ a combination of in situ infra-red and neutron scattering techniques to elucidate the influence of the surface ethoxy species in the activation of both diethyl ether and ethanol, towards the improved formation of ethylene at low temperatures. The combined conclusions of these studies is that the formation of ethylene is the rate determining step, which is of fundamental importance towards the development of this process and the introduction of bio-ethanol as a viable feedstock for ethylene production.

Methanol oxidation on Fe2O3 catalysts and the effects of surface Mo

The adsorption of methanol on haematite has been investigated using temperature programmed methods, combined with in situ DRIFTS. Model catalysts based on this material have then been made with a shell-core configuration of molybdenum oxide monolayers on top of the haematite core. These are used as models of industrial iron molybdate catalysts, used to selectively oxidise methanol to formaldehyde, one of the major chemical outlets for methanol. Haematite itself is completely ineffective in this respect since it oxidises it to CO2 and the DRIFTS shows that this occurs by oxidation of methoxy to formate at around 200 °C. The decomposition behaviour is affected by the absence or presence of oxygen in the gas phase; oxygen destabilises the methoxy and enhances formate production. In contrast, when a monolayer of molybdena is placed onto the surface by incipient wetness, and it remains there after calcination, the pathway to formate production is blocked and formaldehyde is the main gas phase product in TPD after methanol dosing.

The application of inelastic neutron scattering to investigate the interaction of methyl propanoate with silica

A modern industrial route for the manufacture of methyl methacrylate involves the reaction of methyl propanoate and formaldehyde over a silica-supported Cs catalyst. Although the process has been successfully commercialised, little is known about the surface interactions responsible for the forward chemistry. This work concentrates upon the interaction of methyl propanoate over a representative silica. A combination of infrared spectroscopy, inelastic neutron scattering, DFT calculations, X-ray diffraction and temperature-programmed desorption is used to deduce how the ester interacts with the silica surface.

The reaction of formic acid with RaneyTM copper

The interaction of formic acid with Raney^TM Cu proves to be complex. Rather than the expected generation of a monolayer of bidentate formate, we find the formation of a Cu(II) compound. This process occurs by direct reaction of copper and formic acid; in contrast, previous methods are by solution reaction. This is a rare example of formic acid acting as an oxidant rather than, as more commonly found, a reductant. The combination of diffraction, spectroscopic and computational methods has allowed this unexpected process to be characterized.

Proton Conduction in a Phosphonate-Based Metal-Organic Framework Mediated by Intrinsic "Free Diffusion inside a Sphere"

Understanding the molecular mechanism of proton conduction is crucial for the design of new materials with improved conductivity. Quasi-elastic neutron scattering (QENS) has been used to probe the mechanism of proton diffusion within a new phosphonate-based metal–organic framework (MOF) material, MFM-500(Ni). QENS suggests that the proton conductivity (4.5 × 10^–4 S/cm at 98% relative humidity and 25 °C) of MFM-500(Ni) is mediated by intrinsic “free diffusion inside a sphere”, representing the first example of such a mechanism observed in MOFs.

Correction: Evidence for a surface gold hydride on a nanostructured gold catalyst

Correction for 'Evidence for a surface gold hydride on a nanostructured gold catalyst' by I. P. Silverwood et al., Chem. Commun., 2016, 52, 533-536.

Towards microfluidic reactors for in situ synchrotron infrared studies

Anodically bonded etched silicon microfluidic devices that allow infrared spectroscopic measurement of solutions are reported. These extend spatially well-resolved in situ infrared measurement to higher temperatures and pressures than previously reported, making them useful for effectively time-resolved measurement of realistic catalytic processes. A data processing technique necessary for the mitigation of interference fringes caused by multiple reflections of the probe beam is also described.

Structural and spectroscopic characterisation of C4 oxygenates relevant to structure/activity relationships of the hydrogenation of α,β-unsaturated carbonyls

In the present work, we have investigated the conformational isomerism and calculated the vibrational spectra of the C4 oxygenates: 3-butyne-2-one, 3-butene-2-one, 2-butanone and 2-butanol using density functional theory. The calculations are validated by comparison to structural data where available and new, experimental inelastic neutron scattering and infrared spectra of the compounds. We find that for 3-butene-2-one and 2-butanol the spectra show clear evidence for the presence of conformational isomerism and this is supported by the calculations. Complete vibrational assignments for all four molecules are provided and this provides the essential information needed to generate structure/activity relationships for the sequential catalytic hydrogenation of 3-butyne-2-one to 2-butanol.

Room temperature methoxylation in zeolites: insight into a key step of the methanol-to-hydrocarbons process

Neutron scattering methods studying mobility and vibrational spectra observed complete room temperature methoxylation in a commercial sample of methanol-to-hydrocarbons (MTH) catalyst H-ZSM-5.

Neutron scattering in catalysis and energy materials

This themed collection focuses on the use of neutron scattering to interrogate catalysts and materials used in energy storage. Neutrons are a powerful and unique probe that allow characterisation of these experimentally challenging systems in multiple ways.

In situ spectroscopic investigations of MoOx/Fe2O3 catalysts for the selective oxidation of methanol

Methoxy adsorbed on MoO_x/Fe₂O₃ surfaces.

Methanol diffusion in zeolite HY: a combined quasielastic neutron scattering and molecular dynamics simulation study

The diffusion of methanol in zeolite HY is studied using tandem quasielastic neutron scattering (QENS) experiments and molecular dynamics (MD) simulations at 300–400 K.

An attenuated total reflection Fourier transform infrared (ATR FT-IR) spectroscopic study of gas adsorption on colloidal stearate-capped ZnO catalyst substrate

Attenuated total reflection Fourier transform infrared (ATR FT-IR) spectroscopy has been applied in situ to study gas adsorption on a colloidal stearate-capped zinc oxide (ZnO) surface. Infrared spectra of a colloidal stearate-capped ZnO catalyst substrate were assigned at room temperature using zinc stearate as a reference compound. Heating was shown to create a monodentate species that allowed conformational change to occur, leading to altered binding geometry of the stearate ligands upon cooling. CO2 and H2 adsorption measurements demonstrated that the ligand shell was permeable and did not cover the entire surface, allowing adsorption and reaction with at least some portion of the ZnO surface. It has been demonstrated that stearate ligands did not prevent the usual chemisorption processes involved in catalytic reactions on a model ZnO catalyst system, yet the ligand-support system is dynamic under representative reaction conditions.

Experimental arrangements suitable for the acquisition of inelastic neutron scattering spectra of heterogeneous catalysts

Inelastic neutron scattering (INS) is increasingly being used for the characterization of heterogeneous catalysts. As the technique is uniquely sensitive to hydrogen atoms, vibrational spectra can be obtained that emphasize a hydrogenous component or hydrogen-containing moieties adsorbed on to an inorganic support. However, due to sensitivity constraints, the technique typically requires large sample masses (∼10 g catalyst). A reaction system is hereby described that enables suitable quantities of heterogeneous catalysts to be appropriately activated and operated under steady-state conditions for extended periods of time prior to acquisition of the INS spectrum. In addition to ex situ studies, a cell is described which negates the need for a sample transfer stage between reaction testing and INS measurement. This cell can operate up to temperatures of 823 K and pressures up to 20 bar. The apparatus is also amenable to adsorption experiments at the gas-solid interface.

A microwave-heated infrared reaction cell for the in situ study of heterogeneous catalysts

A transmission infrared microreactor cell which holds a pressed disc in a controlled atmosphere and allows microwave and conventional heating up to 423 K is demonstrated using the oxidation of carbon monoxide over the standard catalyst EUROPT-1. Optical characteristics are determined by the choice of CaF2 as the window material, allowing transmission from 77,000-1000 cm(-1). An oscillating microwave power regime with a peak height of 200 W is used and time-resolved infrared spectra and mass spectrometry show oscillations in the reaction which correspond to the microwave heating.