A DFT Study on Isomorphously Substituted MCM-22 Zeolite

Acidity of Isomorphic Substituted Zeolites with B, Al and Ga Revisited

Isomorphic substitution of zeolites with B, Al and Ga is a widely used approach in catalysis. The experimentally reported trend of their acidities decreases in the order: Al>Ga>B. However, a consistent explanation is still lacking in the literature. To bring more understanding of this trend, density functional theory computations were conducted on several model systems. First, the acidity of small clusters with two (2T) and five (5T) tetrahedral sites was analyzed. These systems were then projected onto three large void structures: H-[A]-BEA (52T), H-[A]-FAU (84T) and H-[A]-MOR (112T) with A=B, Al, Ga. Our electron density and Interacting Quantum Atom analyses show that the acidity of Al zeolites originates from the much stronger O-Al bond, which is dominated by the electrostatic attraction. The bridging hydroxyl therefore donates more charge density to the metal, the proton becomes more positive and consequently more acidic. Ga zeolites are more acidic than B zeolites due to the greater covalent nature on the O-Ga bond. The resulting acidity, as seen by ammonia, depends on both the acidic oxygen and the charge distribution of the surrounding oxygens exerted by the substituents.

DFT stimulation and experimental insights of chiral Cu(II)-salen scaffold within the pocket of MWW-zeolite and its catalytic study

A Cu(II)-salen complex encapsulated in MWW-framework as an efficient chiral organocatalyst was developed for the synthesis of 3,4-dihydropyrimidin-2-(1H)-one (DHPMs) derivatives via an asymmetric pathway. In order to confirm its structural properties, single-crystal X-ray diffraction, powder XRD, BET, XPS, FE-SEM, EDX, UV-Vis, and FTIR spectra were used. Using computer-assisted DFT calculations, the Cu(II)-salen complex has been fine-tuned to fit into the pocket of the porous MWW support while keeping its chirality. This organocatalyst was shown to be a potent catalyst for the formation of the desired DHPMs product under short reaction times. Furthermore, this green protocol allows rapid and simple isolation of active MWW-trapped Cu(II)-salen scaffolds and its reusability in at least five consecutive runs without losing much of its activity.

Interlayer spacing in pillared and grafted MCM-22 type silicas: density functional theory analysis versus experiment

Pillaring of synthetic layered crystalline silicates and aluminosilicates provides a strategy to enhance their adsorption and separation performance, and can facilitate the understanding of such behavior in more complex natural clays. We perform the first-principles density functional theory calculations for the pillaring of the pure silica polymorph of an MCM-22 type molecular sieve. Starting with a precursor material MCM-22P with fully hydroxylated layers, a pillaring agent, (EtO)₃SiR, can react with hydroxyl groups (-OH) on adjacent internal surfaces, 2(-OH) + (EtO)₃SiR + H₂O → (-O)₂SiOHR + 3EtOH, to form a pillar bridging these surfaces, or with a single hydroxyl, -OH + (EtO)₃SiR + 2H₂O → (-O)Si(OH)₂R + 3EtOH, grafting to one surface. For computational efficiency, we replace the experimental organic ligand, R, by a methyl group. We find that the interlayer spacing in MCM-22 is reduced by 2.66 Å relative to weakly bound layers in the precursor MCM-22P. Including (-O)₂SiR bridges for 50% (100%) of the hydroxyl sites in MCM-22P increases the interlayer spacing relative to MCM-22 by 2.52 Å (2.46 Å). For comparison, we also analyze the system where all -OH groups in MCM-22P are replaced by non-bridging grafted (-O)Si(OH)₂R which results in a smaller interlayer spacing expansion of 2.17 Å relative to MCM-22. Our results for the interlayer spacing in the pillared materials are compatible with experimental observations for a similar MCM-22 type material with low Al content (Si : Al = 51 : 1) of an expansion relative to MCM-22 of roughly 2.8 Å and 2.5 Å from our x-ray diffraction and scanning transmission electron microscopy analyses, respectively. The latter analysis reveals significant variation in individual layer spacings.

Adsorption and dehydration of ethanol on isomorphously B, Al, and Ga substituted H-ZSM-5 zeolite: an embedded ONIOM study

Dehydration reactions are important in the petroleum and petrochemical industries, especially for the feedstock production. In this work, the catalytic activity of zeolites with different acidities for the dehydration of ethanol to ethylene and diethylether is investigated by density functional calculations on cluster models of three isomorphous B, Al, and Ga substituted H-ZSM-5 zeolites. Both unimolecular and bimolecular mechanisms are investigated. Detailed reaction profiles for the dehydration reaction, assuming either a stepwise or a concerted mechanism, were calculated by using the ONIOM(MP2:M06-2X) + SCREEP method. The adsorption energies of ethanol are -21.6, -28.1, and -27.7 kcal mol^-1 on H-[B]-ZSM-5, H-[Al]-ZSM-5, and H-[Ga]-ZSM-5 zeolites, respectively. The activation energies for the rate-determining step of the unimolecular concerted mechanism for the ethylene formation are 48.5, 42.6, and 43.6 kcal mol^-1 on H-[B]-ZSM-5, H-[Al]-ZSM-5, and H-[Ga]-ZSM-5 zeolites, respectively. The activation energies for the ethoxy formation as the rate-determining step for the bimolecular formation of diethylether are 42.3, 40.0, and 41.1 kcal mol^-1 on H-[B]-ZSM-5, H-[Al]-ZSM-5, and H-[Ga]-ZSM-5 zeolites, respectively. The results indicate that the catalytic activities for the dehydration of ethanol decrease in the order H-[Al]-ZSM-5 ~ H-[Ga]-ZSM-5 > H-[B]-ZSM-5. Besides the acid strength, the zeolite framework affects the reaction by stabilizing the reaction intermediates, leading to more stable adsorption complexes and lower activation barriers.

Deboronation-assisted construction of defective Ti(OSi)3OH species in MWW-type titanosilicate and their enhanced catalytic performance

Regulating the state of titanium species via the deboronation-assisted route is a facile strategy to construct highly efficient titanosilicate catalysts.

A DFT study on catalytic oxidative desulfurization with H2O2 over Ti-MWW zeolite

The catalytic mechanism of Ti-MWW in oxidative desulfurization with H₂O₂ was investigated by quantum chemical calculations. A defect model (Ti-d) and a perfect model (Ti-p) were proposed for Ti-MWW, and two possible reaction pathways starting from Ti-d and Ti-p were considered. On Ti-d, the hydroperoxy bidentate intermediate TiOOH (η²) was formed by activating H₂O₂ at the Ti center. Afterwards, aromatic sulfides were oxidized to sulfoxides and to ultimate sulfones by TiOOH (η²). The order of oxidation reactivity was benzothiophene > dibenzothiophene > thiophene, conforming to experimental observations. The Ti-p pathway proposed for oxidation of sulfides with H₂O₂ resulted in higher energy barriers compared to the Ti-d pathway. Natural bond orbital charge analysis was carried out to understand the charge distribution. This work showed that the defective Ti-MWW model for oxidative desulfurization was more active than the perfect model. Graphical abstract Catalytic oxidative desulfurizationwith H₂O₂ over Ti-MWW 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.

Location, distribution and acidity of Al substitution in ZSM-5 with different Si/Al ratios – a periodic DFT computation

Systematic periodic density functional theory computations including dispersion correction (GGA-PBE-D3) have been carried out to characterize the location, distribution and acidity of Al substitution in ZSM-5 with different Si/Al ratios (up to 8 Al atoms).

A comparative first principles study on trivalent ion incorporated SSZ-13 zeolites

The dispersion-corrected density functional theory has been used to study the trivalent ions B, Al, Ga, and Fe incorporated SSZ-13-type zeolites.

The structure and activity of potassium supported on SBA-15 molecular sieve: Density functional theory study

The geometries, vibrational frequencies, electronic properties and reactivity of potassium supported on SBA-15 have been theoretically investigated by the density functional theory (DFT) method. The structural model of the potassium supported on SBA-15 was constructed based on our previous work [Wang ZX, Wang DX, Zhao Z, Chen Y, Lan J, A DFT study of the structural units in SBA-15 mesoporous molecular sieve, Comput. Theor. Chem.963, 403, 2011]. This paper is the extension of our previous work. The most favored location of potassium atom was obtained by the calculation of substitution energy. The calculated vibrational frequencies of K /SBA-15 are in good agreement with the experimental results. By analyzing the properties of electronic structure, we found that the O atom of Si - O (2)- K group acts as the Lewis base center and the K atom acts as the Lewis acid center. The reactivity of K /SBA-15 was investigated by calculating the activation of oxygen molecule. The oxygen molecule can be activated by K /SBA-15 with an energy barrier of 103.2 kJ/mol. In the final state, the activated oxygen atoms become new Lewis acid centers, which are predicted to act as the active sites in the catalytic reactions. This study provides a deep insight into the properties of supported potassium catalysts and offers fundamental information for further research.

Density functional calculations on the distribution, acidity, and catalysis of Ti(IV) and Ti(III) ions in MCM-22 zeolite

Isolated Ti species in zeolites show unique catalytic activities for a variety of chemical reactions. In this work, density functional calculations were used to explore three current concerns: 1) the distributions of Ti(IV) and Ti(III) ions in the MCM-22 zeolite; 2) the Lewis acidity of the Ti(IV) and Ti(III) sites; and 3) activation of alkane C-H bonds by photocatalysis with Ti-doped zeolites. Neither the Ti(IV) nor Ti(III) ions are randomly distributed in the MCM-22 zeolite. The orders of relative stability are very close for the eight Ti(IV) and Ti(III) sites, and the T3 site is the most probable in both cases. The wavelengths for Ti(IV)-Ti(III) excitations were calculated to lie in the range λ=246.9-290.2 nm. The Ti3(IV) site shows Lewis acidity toward NH(3) in two different modes, and these two modes can coexist with each other. The calculated Ti(IV) coordination numbers, Ti(IV)-O bond elongations, and charge transfers caused by NH(3) adsorption are in good agreement with previous results. Similarly, two different NH(3) adsorption modes exist for the Ti3(III) site; the site that exhibits radical transfer from the lattice O to N atoms is preferred due to the higher adsorption energy. This indicates that the Ti3(III) site does not show Lewis acidity, in contrast to the Ti3(IV) site. At the Ti3(III) site, the energy barrier for activating the methane C-H bond was calculated to be 33.3 kJ mol(-1) and is greatly reduced by replacing the hydrogen atoms with methyl groups. In addition, the reactivity is improved when switching from MCM-22 to TS-1 zeolite. The studies on the various Ti species reveal that lattice O atoms rather than Ti(III) radicals are crucial to the activation of alkane C-H bonds. This work provides new insights into and aids understanding of the catalysis by isolated Ti species in zeolites.

Brønsted sites on acid-treated montmorillonite: a theoretical study with probe molecules

The effect of isomorphic substitution on Brønsted acid sites, generated during acid treatment of montmorillonite, was studied using the density functional theory (DFT) and the cluster model method. Cluster models describing acid montmorillonite were constructed and used as a matrix for the protonation of several test molecules under different types of isomorphic substitution. Our results concerning geometries, energies, and infrared vibrational frequencies indicate that isomorphic substitution of aluminum with magnesium in the octahedral sheet leads to stronger acid sites and more effective protonation of test molecules than those produced by substitution of silicon with aluminum in the tetrahedral sheet.

Behavior of Extraframework Fe Sites in MFI and MCM-22 Zeolites upon Interaction with N2O and NO

We report on the characterization of an isomorphously substituted Fe-MCM-22 sample containing both Fe and Al in framework positions (Si/Fe = 44, Si/Al = 25). XANES spectroscopy was used to study the evolution of Fe sites as a consequence of thermal activation at high temperature (1073 K) and subsequent oxidation with N2O. The results were compared to those obtained in the same conditions on a well-known Fe-silicalite sample (Si/Fe = 68, Si/Al = infinity). In both samples, thermal activation causes migration of a fraction of Fe ions from framework to extraframework positions, this migration being accompanied by a reduction of Fe3+ to Fe2+. Upon oxidation with N2O at 523 K, the two samples show a different behavior. While in Fe-silicalite practically all of the Fe2+ sites formed by thermal activation are reoxidized to Fe3+, in Fe-MCM-22 only a fraction of the extraframework iron sites is involved in the reoxidation process. The accessibility of the extraframework Fe sites was also investigated by using the NO molecule as a surface probe. Upon NO dosage on the sample, the modification of the pre-edge peak and of the edge position suggests an important charge release from the extraframework Fe2+ ions to the adsorbed molecules. This could be formalized with the formation of Fe3+(NO-) complexes, compatible (on the basis of the simple molecular orbital theory) with a bent NO geometry. The formation of a complex family of Fe2+ mono-, di-, and trinitrosyl complexes was also confirmed by FTIR spectroscopy. Similarly to what was observed in the oxidation experiments, the fraction of extraframework Fe sites able to interact with NO in Fe-MCM-22 sample is smaller than that in Fe-silicalite treated in the same conditions. This trend is explained with a major clustering of extraframework Fe sites in Fe-MCM-22 sample, as was also suggested by FTIR experiments. These results suggest that the dispersion of iron in zeolitic matrixes prepared by isomorphous substitution could also depend on the zeolitic structure.