Probing the Acid Strength of Brønsted Acidic Zeolites with Acetonitrile: An Atomistic and Quantum Chemical Study

Simulating Crystal Structure, Acidity, Proton Distribution, and IR Spectra of Acid Zeolite HSAPO-34: A High Accuracy Study

It is a challenge to characterize the acid properties of microporous materials in either experiments or theory. This study presents the crystal structure, acid site, acid strength, proton siting, and IR spectra of HSAPO-34 from the SCAN + rVV10 method. The results indicate: the crystal structures of various acid sites of HSAPO-34 deviate from the space group of R3¯; the acid strength inferred from the DPE value likely decreases with the proton binding sites at O(2), O(4), O(1),and O(3), contrary to the stability order in view of the internal energy; the calculated ensemble-averaged DPE is about 1525 kJ/mol at 673.15 K; and the proton siting and the proton distribution are distinctly influenced by the temperature: at low temperatures, the proton is predominantly located at O(3), while it prefers O(2) at high temperatures, and the proton at O(4) assumedly has the least distribution at 273.15-773.15 K. In line with the neutron diffraction experiment, a correction factor of 0.979 is needed to correct for the calculated hydroxyl stretching vibration (ν(O-H)) of HSAPO-34. It seems that the SCAN meta-GGA method, compensating for some drawbacks of the GGA method, could provide satisfying results regarding the acid properties of HSAPO-34.

Identifying and controlling the acid site distributions in mordenite zeolite for dimethyl ether carbonylation reaction by means of selective ion-exchange

The identification and regulation of acid site distributions in mordenite zeolites are accomplished by using tetramethylammonium ions to selectively exchange with the counter ions in 12-membered channels.

Can Hammett indicators accurately measure the acidity of zeolite catalysts with confined space? Insights into the mechanism of coloration

The acidic properties of zeolite catalysts play a crucial role in governing catalytic performances, which makes the acidity characterization an important subject in the field of zeolite catalysis.

Identifying the effective phosphorous species over modified P-ZSM-5 zeolite: a theoretical study

In this work, a density functional theory (DFT) study was carried out to address the fundamental description of the effective phosphorous species that could improve the framework stability and reduce the coke deposition formation on the P-ZSM-5 zeolite.

Scaling Relations for Acidity and Reactivity of Zeolites

Zeolites are widely applied as solid acid catalysts in various technological processes. In this work we have computationally investigated how catalytic reactivity scales with acidity for a range of zeolites with different topologies and chemical compositions. We found that straightforward correlations are limited to zeolites with the same topology. The adsorption energies of bases such as carbon monoxide (CO), acetonitrile (CH3CN), ammonia (NH3), trimethylamine (N(CH3)3), and pyridine (C5H5N) give the same trend of acid strength for FAU zeolites with varying composition. Crystal orbital Hamilton populations (COHP) analysis provides a detailed molecular orbital picture of adsorbed base molecules on the Brønsted acid sites (BAS). Bonding is dominated by strong σ donation from guest molecules to the BAS for the adsorbed CO and CH3CN complexes. An electronic descriptor of acid strength is constructed based on the bond order calculations, which is an intrinsic parameter rather than adsorption energy that contains additional contributions due to secondary effects such as van der Waals interactions with the zeolite walls. The bond order parameter derived for the CH3CN adsorption complex represents a useful descriptor for the intrinsic acid strength of FAU zeolites. For FAU zeolites the activation energy for the conversion of π-adsorbed isobutene into alkoxy species correlates well with the acid strength determined by the NH3 adsorption energies. Other zeolites such as MFI and CHA do not follow the scaling relations obtained for FAU; we ascribe this to the different van der Waals interactions and steric effects induced by zeolite framework topology.

External or internal surface of H-ZSM-5 zeolite, which is more effective for the Beckmann rearrangement reaction?

The influence of Brønsted acid site location on the Beckmann rearrangement reaction over H-ZSM-5 zeolite has been explored.

Adsorption in zeolites using mechanically embedded ONIOM clusters

We have explored mechanically embedded three-layer QM/QM/MM ONIOM models for computational studies of binding in Al-substituted zeolites.

Probing zeolites by vibrational spectroscopies

This review addresses the most relevant aspects of vibrational spectroscopies (IR, Raman and INS) applied to zeolites and zeotype materials. Surface Brønsted and Lewis acidity and surface basicity are treated in detail. The role of probe molecules and the relevance of tuning both the proton affinity and the steric hindrance of the probe to fully understand and map the complex site population present inside microporous materials are critically discussed. A detailed description of the methods needed to precisely determine the IR absorption coefficients is given, making IR a quantitative technique. The thermodynamic parameters of the adsorption process that can be extracted from a variable-temperature IR study are described. Finally, cutting-edge space- and time-resolved experiments are reviewed. All aspects are discussed by reporting relevant examples. When available, the theoretical literature related to the reviewed experimental results is reported to support the interpretation of the vibrational spectra on an atomic level.

Fourier self-deconvolution of the IR spectra as a tool for investigation of distinct functional groups in porous materials: Brønsted acid sites in zeolites

For many decades, IR and FT-IR spectroscopy has generated valuable information about different functional groups in zeolites, metal-organic frameworks (MOFs), and other porous materials. However, this technique cannot distinguish between functional groups in different local environments. Our study demonstrates that this limitation could be overcome by using Fourier self-deconvolution of infrared spectra (FSD-IR). We apply this method to study three acidic mordenite zeolites and show (i) that these zeolites contain six distinct Brønsted acid sites (BAS) as opposed to 2-4 different BAS previously considered in literature and (ii) that the relative amounts of these BAS are different in the three zeolites examined. We then analyze possible locations of six BAS in the mordenite structure and explain a number of conflicting results in literature. On this basis, we conclude that the FSD-IR method allows direct visualization and examination of distributions of distinct BAS in zeolites, thus providing a unique research opportunity, which no other method can provide. Given the similarities in the IR analysis of different functional groups in solids, we expect that the FSD-IR method will be also instrumental in the research into other porous materials, such as solid oxides and MOFs. The latter point is illustrated by FSD of the IR spectrum of hydroxyl groups in a sample of α-alumina.

Effects of the zeolite framework on the adsorptions and hydrogen-exchange reactions of unsaturated aliphatic, aromatic, and heterocyclic compounds in ZSM-5 zeolite: a combination of perturbation theory (MP2) and a newly developed density functional theory (M06-2X) in ONIOM scheme

The confinement effect on the adsorption and reaction mechanism of unsaturated aliphatic, aromatic and heterocyclic compounds on H-ZSM-5 zeolite has been investigated by the four ONIOM methods (MP2:M06-2X), (MP2:B3LYP), (MP2:HF), and (MP2:UFF). The H-ZSM-5 'nanoreactor' porous intersection, where chemical reactions take place, is represented by a quantum cluster of 34 tetrahedral units. Ethene, benzene, ethylbenzene, and pyridine are chosen to represent reactions of various adsorbates of aliphatic, aromatic and heterocyclic compounds. Among the four combined methods, (MP2:M06-2X) outperforms the others. The results confirm that the method that takes weak interactions, especially the van der Waals interaction, into account is essential for describing the confinement effect from the zeolite framework. The effects of the infinite zeolitic framework on the cluster model are also included by a set of point charges generated by the embedded ONIOM model. The energies for the adsorption of ethene, benzene, ethylbenzene, and pyridine on H-ZSM-5 from an embedded ONIOM(MP2:M06-2X) calculation are predicted to be -14.0, -19.8, -24.7, and -48.4 kcal/mol, respectively, which are very close to available experimental observations. The adsorption energy of pyridine agrees well with the experiment data of -47.6 kcal/mol. We also applied the same computational methodology on the systematic investigation of the H/H exchange reaction of benzene and ethylbenzene with the acidic H-ZSM-5 zeolite. The H/H exchange reaction was found to take place in a single concerted step. The calculated apparent activation energies for benzene and ethylbenzene are 12.6 and 4.9 kcal/mol, which can be compared to the experimental estimates of 11.0 and 6.9 kcal/mol, respectively. The confinement effect of the extended zeolite framework has been clearly demonstrated not only to stabilize the adsorption complexes but also to improve their corresponding activation energies to approach the experimental benchmark.

Computational and FTIR spectroscopic studies on carbon monoxide and dinitrogen adsorption on a high-silica H-FER zeolite

Adsorption (at a low temperature) of carbon monoxide and dinitrogen on a high-silica ferrierite-type zeolite (H-FER, Si : Al = 27.5 : 1) was investigated by means of variable temperature infrared spectroscopy and theoretical calculations at the periodic DFT level. This combined experimental and computational approach led to detailed characterization of several types of hydrogen-bonded OHCO and OHN(2) complexes, formed by interaction between the adsorbed molecules and the Brønsted acid OH groups of the zeolite. CO or N(2), forming linear complexes with OH groups pointing towards a sufficiently ample void space, show the largest adsorption enthalpy which was found to be in the (approximate) range of -25 to -29 kJ mol(-1) for CO and -15 to -19 kJ mol(-1) for N(2). Less stable OHCO or OHN(2) complexes can be formed when either the Brønsted acid OH group is involved in intra-zeolite hydrogen bonding or when the free space available is too small to allow formation of linear complexes without previous re-location of the proton of the OH group involved. The details of experimental IR spectra in the O-H, C-O, and N-N stretching regions could be interpreted on the basis of good agreement between experimental and calculated results.

Combined DFT Theoretical Calculation and Solid-State NMR Studies of Al Substitution and Acid Sites in Zeolite MCM-22

The distributions of Brönsted acidic protons and their acid strengths in zeolite H-MCM-22 have been characterized by density functional theory (DFT) calculations as well as magic angle spinning (MAS) NMR experiments. The embedded scheme (ONIOM) that combines the quantum mechanical (QM) description of active sites and semiempirical AM1 treatment of the neighboring environment was applied to predict the aluminum substitution mechanism and proton affinity (PA), as well as adsorption behaviors of acetone and trimethylphosphine oxide (TMPO) onto the zeolite. Our theoretical results indicate that the Al substitution takes place in the order of Al1-OH-Si2 > Al8-OH-Si8 > Al5-OH-Si7. The DFT theoretical and NMR experimental results suggest that the acid strength of the three Brönsted acid sites in H-MCM-22 zeolite is slightly lower than that of H-ZSM-5 zeolite and the accessible Brönsted acidic protons most likely reside in both the supercages (at the Al8-OH-Si8 and Al1-OH-Si2 sites) and external surface pocket (at the Al8-OH-Si8 site) rather than in the sinusoidal channels (Al5-OH-Si7), with the Al1-OH-Si2 site having the strongest acid strength (as probed by TMPO). This may partially explain the special selectivity of acid-catalyzed reactions occurring inside the channels of H-MCM-22 zeolite.

Design of a "green" one-step catalytic production of epsilon-caprolactam (precursor of nylon-6)

The ever-increasing industrial demand for nylon-6 (polycaprolactam) necessitates the development of environmentally benign methods of producing its precursor, epsilon-caprolactam, from cyclohexanone. It is currently manufactured in two popular double-step processes, each of which uses highly aggressive reagents, and each generates substantial quantities of largely unwanted ammonium sulfate as by-product. Here we describe a viable laboratory-scale, single-step, solvent-free process of producing epsilon-caprolactam using a family of designed bifunctional, heterogeneous, nanoporous catalysts containing isolated acidic and redox sites, which smoothly convert cyclohexanone to epsilon-caprolactam with selectivities in the range 65-78% in air and ammonia at 80 degrees C. The catalysts are microporous (pore diameter 7.3 A) aluminophosphates in which small fractions of the Al(III)O4(5-) and P(V)O4(3-) tetrahedra constituting the 4-connected open framework are replaced by Co(III)PO4(5-) and Si(IV)O4(4-) tetrahedra, which become the loci of the redox and acidic centers, respectively. The catalysts may be further optimized, and already may be so designed as to generate selectivities of approximately 80% for the intermediate oxime, formed from NH2OH, which is produced in situ within the pore system. The advantages of such designed heterogeneous catalysts, and their application to a range of other chemical conversions, are also adumbrated.

Ammonia IRMS-TPD Study on the Distribution of Acid Sites in Mordenite

Using an IRMS-TPD (temperature programmed desorption) of ammonia, we studied the nature, strength, crystallographic location, and distribution of acid sites of mordenite. In this method, infrared spectroscopy (IR) and mass spectroscopy (MS) work together to follow the thermal behavior of adsorbed and desorbed ammonia, respectively; therefore, adsorbed species were identified, and their thermal behavior was directly connected with the desorption of ammonia during an elevation of temperature. IR-measured TPD of the NH4(+) cation was similar to MS-measured TPD, thus showing the nature of Brønsted acidity. From the behavior of OH bands, it was found that the Brønsted acid sites consisted of two kinds of OH bands at high and low wavenumbers, ascribable to OH bands situated on 12- and 8-member rings (MR) of mordenite structure, respectively. The amount and strength of these Brønsted hydroxyls were measured quantitatively based on a theoretical equation using a curve fitting method. Up to ca. 30% of the exchange degree, NH4(+) was exchanged with Na+ on the 12-MR to arrive at saturation; therefore, in this region, the Brønsted acid site was situated on the large pore of 12-MR. The NH4(+) cation was then exchanged with Na+ on 8-MR, and finally exceeded the amount on 12-MR. In the 99% NH4-mordenite, Brønsted acid sites were located predominantly on the 8-MR more than on the 12-MR. Irrespective of the NH4(+) exchange degree, the strengths deltaH of Brønsted OH were 145 and 153 kJ mol(-1) on the 12- and 8-MR, respectively; that is, the strength of Brønsted acid site on the 8-MR was larger than that on the 12-MR. A density functional theory (DFT) calculation supported the difference in the strengths of the acid sites. Catalytic cracking activity of the Brønsted acid sites on the 8-MR declined rapidly, while that on the 12-MR was remarkably kept. The difference in strength and/or steric capacity may cause such a difference in the life of a catalyst.

Probing the Structural and Electronic Factors Affecting the Adsorption and Reactivity of Alkenes in Acidic Zeolites Using DFT Calculations and Multivariate Statistical Methods

Quantum mechanical (QM) cluster calculations have been performed on a model of ZSM-5 at DFT and MP2 levels. We investigated how the adsorption energies and the energetics of alkoxide intermediate formation of six different alkene substrates, ethene, propene, 1-butene, cis/trans butene, and isobutene, vary in this zeolite model. An analysis of the DFT geometric, electronic, and energetic parameters of the zeolite-substrate complexes, transition states, and alkoxide intermediates is performed using principal components analysis (PCA) and partial least squares (PLS). These deliver an insight into the correlated changes that occur between molecular structure and energy along the reaction coordinate between the physisorbed and chemisorbed species within the zeolite. To the best of our knowledge, this is the first occasion multivariate techniques such as PCA or PLS have been employed to profile the changes in electronics, distances, and angles in QM calculations of catalytic systems such as zeolites. We find the calculated adsorption and the alkoxide intermediate energies correlate strongly with the absolute charge on the substrate and the length of the substrate double bond. The transition states' energies are not affected by the zeolite framework as modeled, which explains why they correlate strongly with the gas-phase substrate protonation energy. Our cluster results show that for ethene, propene, 1-butene, and isobutene, the relative energetics associated with the formation of the alkoxide intermediate in ZSM-5 follow the same trends as calculations where the effects of the framework are included.