Mechanisms of Hydrogen Exchange of Methane with H-Zeolite Y: An ab Initio Embedded Cluster Study

Interfacing the Core-Shell or the Drude Polarizable Force Field With Car-Parrinello Molecular Dynamics for QM/MM Simulations

We report a quantum mechanics/polarizable–molecular mechanics (QM/p–MM) potential based molecular dynamics (MD) technique where the core–shell (or the Drude) type polarizable MM force field is interfaced with the plane-wave density functional theory based QM force field which allows Car–Parrinello MD for the QM subsystem. In the QM/p-MM Lagrangian proposed here, the shell (or the Drude) MM variables are treated as extended degrees of freedom along with the Kohn–Sham (KS) orbitals describing the QM wavefunction. The shell and the KS orbital degrees of freedom are then adiabatically decoupled from the nuclear degrees of freedom. In this respect, we also present here the Nosé–Hoover Chain thermostat implementation for the dynamical subsystems. Our approach is then used to investigate the effect of MM polarization on the QM/MM results. Especially, the consequence of MM polarization on reaction free energy barriers, defect formation energy, and structural and dynamical properties are investigated. A low point charge polarizable potential (p–MZHB) for pure siliceous systems is also reported here.

Quantum chemical ab initio prediction of proton exchange barriers between CH4 and different H-zeolites

A hybrid MP2:DFT (second-order Møller-Plesset perturbation theory-density functional theory) method that combines MP2 calculations for cluster models with DFT calculations for the full periodic structure is used to localize minima and transition structures for proton jumps at different Brønsted sites in different frameworks (chabazite, faujasite, ferrierite, and ZSM-5) and at different crystallographic positions of a given framework. The MP2 limit for the periodic structures is obtained by extrapolating the results of a series of cluster models of increasing size. A coupled-cluster (CCSD(T)) correction to MP2 energies is calculated for cluster models consisting of three tetrahedra. For the adsorption energies, this difference is small, between 0.1 and 0.9 kJ/mol, but for the intrinsic proton exchange barriers, this difference makes a significant (10.85 ± 0.25 kJ/mol) and almost constant contribution across different systems. The total values of the adsorption energies vary between 22 and 34 kJ/mol, whereas the total proton exchange energy barriers fall in the narrow range of 152-156 kJ/mol. After adding nuclear motion contributions (harmonic approximation, 298 K), intrinsic enthalpy barriers between 134 and 141 kJ/mol and apparent energy barriers between 105 and 118 kJ/mol are predicted for the different sites examined for the different frameworks. These predictions are consistent with experimental results available for faujasite, ferrierite, and ZSM-5.

Significant influence of Zn on activation of the C-H bonds of small alkanes by Brønsted acid sites of zeolite

Herein, we analyze earlier obtained and new data about peculiarities of the H/D hydrogen exchange of small C(1)-n-C(4) alkanes on Zn-modified high-silica zeolites ZSM-5 and BEA in comparison with the exchange for corresponding purely acidic forms of these zeolites. This allows us to identify an evident promoting effect of Zn on the activation of C-H bonds of alkanes by zeolite Brønsted sites. The effect of Zn is demonstrated by observing the regioselectivity of the H/D exchange for propane and n-butane as well as by the increase in the rate and a decrease in the apparent activation energy of the exchange for all C(1)-n-C(4) alkanes upon modification of zeolites with Zn. The influence of Zn on alkane activation has been rationalized by dissociative adsorption of alkanes on Zn oxide species inside zeolite pores, which precedes the interaction of alkane with Brønsted acid sites.

EPR and DFT study on the stabilization of radiation-generated methyl radicals in dehydrated Na-A zeolite

Electron paramagnetic resonance (EPR) spectroscopy was applied to study paramagnetic species stabilized in Na-A zeolite exposed to gaseous methane and gamma-irradiated at 77 K. Two types of EPR spectra were recorded during thermal annealing of zeolite up to room temperature. Owing to the results for the zeolite exposed to (13)CH(4) the multiplet observed at 110 K was assigned to a (.-)CH(3)...Na(+) complex. After decay of the multiplet, the isotropic quartet of methyl radical was recorded in the temperature range of 170-280 K. On the basis of the EPR parameters it is postulated that (.-)CH(3) radicals in this temperature region are able to freely rotate inside the zeolite cage. The structures of the (.-)CH(3)...Na(+) adsorption complex and respective hyperfine coupling constants were calculated by applying DFT quantum chemical methods. Two different models were applied to represent the zeolite framework: the 6T structure of one six-membered ring and the 3T cluster. The hyperfine coupling constants calculated for the (.-)CH(3)...Na(+) adsorption complex for both applied models show very good agreement with those obtained experimentally.

Neopentane and Solid Acids: Direct Hydron Exchange before Cracking

The hydrogen/deuterium exchange reaction of 2,2-dimethylpropane (neopentane) over D(2)O-exchanged zeolites (MOR, FAU, BEA, MFI) using a batch recirculation reactor was studied by means of gas chromatography coupled with mass spectrometer. In the temperature range 473-573 K, H/D exchange proceeds without side reaction such as cracking at short contact times. Indeed the C-H bond has appeared favorably involved in the activation of neopentane compared to the less accessible C-C bond. The transition state allowing hydron exchange is most likely a carbonium species (pentacoordinated carbon) as in the case of the H/D exchange between methane and solid acid. The activation energies of the H/D exchange between neopentane and zeolites are the same for all zeolites indicating a common carbonium ion type transition state. On the basis of previous results in the case of the exchange between methane and liquid superacids, the deuterium exchange rates in neopentane were tentatively related to the acidity of the solids. However the order of activity MOR > MFI > BEA > FAU seems to be related to the size of the pores, which may suggest the involvement of a confinement effect in the zeolites cavities. Moreover we found that H/D exchange takes also place between neopentane and deuterated sulfated zirconia (SZ) emphasizing its strong acidity.

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.

DFT Investigation of Alkoxide vs Alkylammonium Formation in Amine-Substituted Zeolites

Density functional theory (DFT) cluster calculations were used to describe bifunctional acid-base properties of amine-substituted zeolites containing a Brønsted acid site. Preliminary results (J. Am. Chem. Soc. 2004, 126, 9162) indicated that efficient use of both functional groups might lead to a substantial lowering of activation barriers. In this paper, comparison is made between the alkoxide formation in zeolites containing only oxygen bridges and alkylammonium formation on the bridging NH groups in amine-functionalized zeolites for various guest species, such as methanol, ethene, and chloromethane. The amine functionalization only lowers barriers for SN2 type reactions with otherwise highly strained transition states, as is the case for chloromethane. In these new materials more basic sites are introduced into the zeolite framework, enabling optimal linear SN2 type transition states incorporating various T sites.

The initial stages of solid acid-catalyzed reactions of adsorbed propane. A mechanistic study by in situ MAS NMR

In situ solid-state NMR spectroscopy was employed to study the kinetics of hydrogen/deuterium exchange and scrambling as well as (13)C scrambling reactions of labeled propane over Al(2)O(3)-promoted sulfated zirconia (SZA) catalyst under mild conditions (30-102 degrees C). Three competitive pathways of isotope redistribution were observed during the course of the reaction: (1) a regioselective H/D exchange between acidic protons of the solid surface and the deuterons of the methyl group of propane-1,1,1,3,3,3-d(6), monitored by in situ (1)H MAS NMR; (2) an intramolecular H/D scrambling between methyl deuterons and protons of the methylene group, without exchange with the catalyst surface, monitored by in situ (2)H MAS NMR; (3) a intramolecular (13)C scrambling, by skeletal rearrangement process, favored at higher temperatures, monitored by in situ (13)C MAS NMR. The activation energy of (13)C scrambling was estimated to be very close to that of (2)H scrambling, suggesting that these two processes imply a common transition state, responsible for both vicinal hydride migration and protonated cyclopropane formation. All pathways are consistent with a classical carbenium ion-type mechanism.

A periodic DFT study of intramolecular isomerization reactions of toluene and xylenes catalyzed by acidic mordenite

A periodic density functional theory (DFT) study of the isomerization reactions of toluene and xylene catalyzed by acidic mordenite is reported. Monomolecular isomerization reactions have been considered and analyzed. The different reaction pathways have been discussed in detail. The use of periodic structure calculations allows consideration and analysis of zeolite electrostatic contributions and steric constraints that occur within zeolite micropores. Major differences in the details of protonation reaction pathways are found when periodic structures are used rather than small cluster models of the Brønsted acidic site. Complex relationships are found between zeolite topology and reaction pathways.

A theoretical study of the alkylation reaction of toluene with methanol catalyzed by acidic mordenite

A theoretical study of the alkylation reaction of toluene with methanol catalyzed by the acidic Mordenite (Si/Al = 23) is reported. Cluster DFT as well as periodical structure DFT calculations have been performed. Full reaction energy diagrams of the elementary reaction steps that lead to the formation of the three xylene isomers are given. The use of periodical structure calculations allows one to account for zeolite framework electrostatic contributions and steric constraints that take place in zeolitic catalysts. Especially the steric constraint energy contribution has a significant effect on the energies and bond formation paths. The activation energy barrier of p-xylene formation is found to be approximately 20 kJ/mol lower than the corresponding values for the formation of its isomers. Computed host-guest binding energies according to the DFT method need a correction due to the absence of the dispersive interaction with the zeolite wall. Apparent activation energies obtained with this correction are in good agreement with experimental data.