Experimental and Theoretical Evidence for the Promotional Effect of Acid Sites on the Diffusion of Alkenes through Small-Pore Zeolites.

The diffusion of saturated and unsaturated hydrocarbons is of fundamental importance for many zeolite-catalyzed processes. Transport of small alkenes in the confined zeolite pores can become hindered, resulting in a significant impact on the ultimate product selectivity and separation. Herein, intracrystalline light olefin/paraffin diffusion through the 8-ring windows of zeolite SAPO-34 is characterized by a complementary set of first-principle molecular dynamics simulations, PFG-NMR experiments, and pulse-response temporal analysis of products measurements, yielding information at different length and time scales. Our results clearly show a promotional effect of the presence of Brønsted acid sites on the diffusion rate of ethene and propene, whereas transport of alkanes is found to be insensitive to the presence of acid sites. The enhanced diffusivity of unsaturated hydrocarbons is ascribed to the formation of favorable π-H interactions with acid protons, as confirmed by IR spectroscopy measurements. The acid site distribution is proven to be an important design parameter for optimizing product distributions and separations.

Light Olefin Diffusion during the MTO Process on H-SAPO-34: A Complex Interplay of Molecular Factors.

The methanol-to-olefins process over H-SAPO-34 is characterized by its high shape selectivity toward light olefins. The catalyst is a supramolecular system consisting of nanometer-sized inorganic cages, decorated by Brønsted acid sites, in which organic compounds, mostly methylated benzene species, are trapped. These hydrocarbon pool species are essential to catalyze the methanol conversion but may also clog the pores. As such, diffusion of ethene and propene plays an essential role in determining the ultimate product selectivity. Enhanced sampling molecular dynamics simulations based on either force fields or density functional theory are used to determine how molecular factors influence the diffusion of light olefins through the 8-ring windows of H-SAPO-34. Our simulations show that diffusion through the 8-ring in general is a hindered process, corresponding to a hopping event of the diffusing molecule between neighboring cages. The loading of different methanol, alkene, and aromatic species in the cages may substantially slow down or facilitate the diffusion process. The presence of Brønsted acid sites in the 8-ring enhances the diffusion process due to the formation of a favorable π-complex host-guest interaction. Aromatic hydrocarbon pool species severely hinder the diffusion and their spatial distribution in the zeolite crystal may have a significant effect on the product selectivity. Herein, we unveil how molecular factors influence the diffusion of light olefins in a complex environment with confined hydrocarbon pool species, high olefin loadings, and the presence of acid sites by means of enhanced molecular dynamics simulations under operating conditions.

Structural and Photophysical Properties of Various Polypyridyl Ligands: A Combined Experimental and Computational Study.

Covalent triazine frameworks (CTFs) with polypyridyl ligands are very promising supports to anchor photocatalytic complexes. Herein, we investigate the photophysical properties of a series of ligands which vary by the extent of the aromatic system, the nitrogen content and their topologies to aid in selecting interesting building blocks for CTFs. Interestingly, some linkers have a rotational degree of freedom, allowing both a trans and cis structure, where only the latter allows anchoring. Therefore, the influence of the dihedral angle on the UV-Vis spectrum is studied. The photophysical properties are investigated by a combined computational and experimental study. Theoretically, both static and molecular dynamics simulations are performed to deduce ground- and excited state properties based on density functional theory (DFT) and time-dependent DFT. The position of the main absorption peak shifts towards higher wavelengths for an increased size of the π-system and a higher π-electron deficiency. We found that the position of the main absorption peak among the different ligands studied in this work can amount to 271 nm; which has a significant impact on the photophysical properties of the ligands. This broad range of shifts allows modulation of the electronic structure by varying the ligands and may help in a rational design of efficient photocatalysts.

A Supramolecular View on the Cooperative Role of Brønsted and Lewis Acid Sites in Zeolites for Methanol Conversion.

A systematic molecular level and spectroscopic investigation is presented to show the cooperative role of Brønsted acid and Lewis acid sites in zeolites for the conversion of methanol. Extra-framework alkaline-earth metal containing species and aluminum species decrease the number of Brønsted acid sites, as protonated metal clusters are formed. A combined experimental and theoretical effort shows that postsynthetically modified ZSM-5 zeolites, by incorporation of extra-framework alkaline-earth metals or by demetalation with dealuminating agents, contain both mononuclear [MOH]+ and double protonated binuclear metal clusters [M(µ-OH)2M]2+ (M = Mg, Ca, Sr, Ba, and HOAl). The metal in the extra-framework clusters has a Lewis acid character, which is confirmed experimentally and theoretically by IR spectra of adsorbed pyridine. The strength of the Lewis acid sites (Mg > Ca > Sr > Ba) was characterized by a blue shift of characteristic IR peaks, thus offering a tool to sample Lewis acidity experimentally. The incorporation of extra-framework Lewis acid sites has a substantial influence on the reactivity of propene and benzene methylations. Alkaline-earth Lewis acid sites yield increased benzene methylation barriers and destabilization of typical aromatic intermediates, whereas propene methylation routes are less affected. The effect on the catalytic function is especially induced by the double protonated binuclear species. Overall, the extra-framework metal clusters have a dual effect on the catalytic function. By reducing the number of Brønsted acid sites and suppressing typical catalytic reactions in which aromatics are involved, an optimal propene selectivity and increased lifetime for methanol conversion over zeolites is obtained. The combined experimental and theoretical approach gives a unique insight into the nature of the supramolecular zeolite catalyst for methanol conversion which can be meticulously tuned by subtle interplay of Brønsted and Lewis acid sites.

Reliably Modeling the Mechanical Stability of Rigid and Flexible Metal-Organic Frameworks.

Over the past two decades, metal-organic frameworks (MOFs) have matured from interesting academic peculiarities toward a continuously expanding class of hybrid, nanoporous materials tuned for targeted technological applications such as gas storage and heterogeneous catalysis. These oft-times crystalline materials, composed of inorganic moieties interconnected by organic ligands, can be endowed with desired structural and chemical features by judiciously functionalizing or substituting these building blocks. As a result of this reticular synthesis, MOF research is situated at the intriguing intersection between chemistry and physics, and the building block approach could pave the way toward the construction of an almost infinite number of possible crystalline structures, provided that they exhibit stability under the desired operational conditions. However, this enormous potential is largely untapped to date, as MOFs have not yet found a major breakthrough in technological applications. One of the remaining challenges for this scale-up is the densification of MOF powders, which is generally achieved by subjecting the material to a pressurization step. However, application of an external pressure may substantially alter the chemical and physical properties of the material. A reliable theoretical guidance that can presynthetically identify the most stable materials could help overcome this technological challenge. In this Account, we describe the recent research the progress on computational characterization of the mechanical stability of MOFs. So far, three complementary approaches have been proposed, focusing on different aspects of mechanical stability: (i) the Born stability criteria, (ii) the anisotropy in mechanical moduli such as the Young and shear moduli, and (iii) the pressure-versus-volume equations of state. As these three methods are grounded in distinct computational approaches, it is expected that their accuracy and efficiency will vary. To date, however, it is unclear which set of properties are suited and reliable for a given application, as a comprehensive comparison for a broad variety of MOFs is absent, impeding the widespread use of these theoretical frameworks. Herein, we fill this gap by critically assessing the performance of the three computational models on a broad set of MOFs that are representative for current applications. These materials encompass the mechanically rigid UiO-66(Zr) and MOF-5(Zn) as well as the flexible MIL-47(V) and MIL-53(Al), which undergo pressure-induced phase transitions. It is observed that the Born stability criteria and pressure-versus-volume equations of state give complementary insight into the macroscopic and microscopic origins of instability, respectively. However, interpretation of the Born stability criteria becomes increasingly difficult when less symmetric materials are considered. Moreover, pressure fluctuations during the simulations hamper their accuracy for flexible materials. In contrast, the pressure-versus-volume equations of state are determined in a thermodynamic ensemble specifically targeted to mitigate the effects of these instantaneous fluctuations, yielding more accurate results. The critical Account presented here paves the way toward a solid computational framework for an extensive presynthetic screening of MOFs to select those that are mechanically stable and can be postsynthetically densified before their use in targeted applications.

Dynamic Interplay between Defective UiO-66 and Protic Solvents in Activated Processes.

UiO-66, composed by Zr-oxide inorganic bricks [Zr6 (µ3 -O)4 (µ3 -OH)4 ] and organic terephthalate linkers, is one of the most studied metal-organic frameworks (MOFs) due to its exceptional thermal, chemical, and mechanical stability. Thanks to its high connectivity, the material can withstand structural deformations during activation processes such as linker exchange, dehydration, and defect formation. These processes do alter the zirconium coordination number in a dynamic way, creating open metal sites for catalysis and thus are able to tune the catalytic properties. In this work, it is shown, by means of first-principle molecular-dynamics simulations at operating conditions, how protic solvents may facilitate such changes in the metal coordination. Solvent can induce structural rearrangements in the material that can lead to undercoordinated but also overcoordinated metal sites. This is demonstrated by simulating activation processes along well-chosen collective variables. Such enhanced MD simulations are able to track the intrinsic dynamics of the framework at realistic conditions.

On the importance of anharmonicities and nuclear quantum effects in modelling the structural properties and thermal expansion of MOF-5.

In this article, we investigate the influence of anharmonicities and nuclear quantum effects (NQEs) in modelling the structural properties and thermal expansion of the empty MOF-5 metal-organic framework. To introduce NQEs in classical molecular dynamics simulations, two different methodologies are considered, comparing the approximate, but computationally cheap, method of generalised Langevin equation thermostatting to the more advanced, computationally demanding path integral molecular dynamics technique. For both methodologies, similar results were obtained for all the properties under investigation. The structural properties of MOF-5, probed by means of radial distribution functions (RDFs), show some distinct differences with respect to a classical description. Besides a broadening of the RDF peaks under the influence of quantum fluctuations, a different temperature dependence is also observed due to a dominant zero-point energy (ZPE) contribution. For the thermal expansion of MOF-5, by contrast, NQEs appear to be only of secondary importance with respect to an adequate modelling of the anharmonicities of the potential energy surface (PES), as demonstrated by the use of two differently parametrised force fields. Despite the small effect in the temperature dependence of the volume of MOF-5, NQEs do however significantly affect the absolute volume of MOF-5, in which the ZPE resulting from the intertwining of NQEs and anharmonicities plays a crucial role. A sufficiently accurate description of the PES is therefore prerequisite when modelling NQEs.

Extension of the QuickFF force field protocol for an improved accuracy of structural, vibrational, mechanical and thermal properties of metal-organic frameworks.

QuickFF was originally launched in 2015 to derive accurate force fields for isolated and complex molecular systems in a quick and easy way. Apart from the general applicability, the functionality was especially tested for metal-organic frameworks (MOFs), a class of hybrid materials consisting of organic and inorganic building blocks. Herein, we launch a new release of the QuickFF protocol which includes new major features to predict structural, vibrational, mechanical and thermal properties with greater accuracy, without compromising its robustness and transparent workflow. First, the ab initio data necessary for the fitting procedure may now also be derived from periodic models for the molecular system, as opposed to the earlier cluster-based models. This is essential for an accurate description of MOFs with one-dimensional metal-oxide chains. Second, cross terms that couple internal coordinates (ICs) and anharmonic contributions for bond and bend terms are implemented. These features are essential for a proper description of vibrational and thermal properties. Third, the fitting scheme was modified to improve robustness and accuracy. The new features are tested on MIL-53(Al), MOF-5, CAU-13 and NOTT-300. As expected, periodic input data are proven to be essential for a correct description of structural, vibrational and thermodynamic properties of MIL-53(Al). Bulk moduli and thermal expansion coefficients of MOF-5 are very accurately reproduced by static and dynamic simulations using the newly derived force fields which include cross terms and anharmonic corrections. For the flexible materials CAU-13 and NOTT-300, the transition pressure is accurately predicted provided cross terms are taken into account. © 2018 Wiley Periodicals, Inc.

Influence of a Confined Methanol Solvent on the Reactivity of Active Sites in UiO-66.

UiO-66, composed of Zr-oxide bricks and terephthalate linkers, is currently one of the most studied metal-organic frameworks due to its exceptional stability. Defects can be introduced in the structure, creating undercoordinated Zr atoms which are Lewis acid sites. Here, additional Brønsted sites can be generated by coordinated protic species from the solvent. In this Article, a multilevel modeling approach was applied to unravel the effect of a confined methanol solvent on the active sites in UiO-66. First, active sites were explored with static periodic density functional theory calculations to investigate adsorption of water and methanol. Solvent was then introduced in the pores with grand canonical Monte Carlo simulations, followed by a series of molecular dynamics simulations at operating conditions. A hydrogen-bonded network of methanol molecules is formed, allowing the protons to shuttle between solvent methanol, adsorbed water, and the inorganic brick. Upon deprotonation of an active site, the methanol solvent aids the transfer of protons and stabilizes charged configurations via hydrogen bonding, which could be crucial in stabilizing reactive intermediates. The multilevel modeling approach adopted here sheds light on the important role of a confined solvent on the active sites in the UiO-66 material, introducing dynamic acidity in the system at finite temperatures by which protons may be easily shuttled from various positions at the active sites.

Advances in theory and their application within the field of zeolite chemistry.

Zeolites are versatile and fascinating materials which are vital for a wide range of industries, due to their unique structural and chemical properties, which are the basis of applications in gas separation, ion exchange and catalysis. Given their economic impact, there is a powerful incentive for smart design of new materials with enhanced functionalities to obtain the best material for a given application. Over the last decades, theoretical modeling has matured to a level that model guided design has become within reach. Major hurdles have been overcome to reach this point and almost all contemporary methods in computational materials chemistry are actively used in the field of modeling zeolite chemistry and applications. Integration of complementary modeling approaches is necessary to obtain reliable predictions and rationalizations from theory. A close synergy between experimentalists and theoreticians has led to a deep understanding of the complexity of the system at hand, but also allowed the identification of shortcomings in current theoretical approaches. Inspired by the importance of zeolite characterization which can now be performed at the single atom and single molecule level from experiment, computational spectroscopy has grown in importance in the last decade. In this review most of the currently available modeling tools are introduced and illustrated on the most challenging problems in zeolite science. Directions for future model developments will be given.

QuickFF: A program for a quick and easy derivation of force fields for metal-organic frameworks from ab initio input.

QuickFF is a software package to derive accurate force fields for isolated and complex molecular systems in a quick and easy manner. Apart from its general applicability, the program has been designed to generate force fields for metal-organic frameworks in an automated fashion. The force field parameters for the covalent interaction are derived from ab initio data. The mathematical expression of the covalent energy is kept simple to ensure robustness and to avoid fitting deficiencies as much as possible. The user needs to produce an equilibrium structure and a Hessian matrix for one or more building units. Afterward, a force field is generated for the system using a three-step method implemented in QuickFF. The first two steps of the methodology are designed to minimize correlations among the force field parameters. In the last step, the parameters are refined by imposing the force field parameters to reproduce the ab initio Hessian matrix in Cartesian coordinate space as accurate as possible. The method is applied on a set of 1000 organic molecules to show the easiness of the software protocol. To illustrate its application to metal-organic frameworks (MOFs), QuickFF is used to determine force fields for MIL-53(Al) and MOF-5. For both materials, accurate force fields were already generated in literature but they requested a lot of manual interventions. QuickFF is a tool that can easily be used by anyone with a basic knowledge of performing ab initio calculations. As a result, accurate force fields are generated with minimal effort. © 2015 Wiley Periodicals, Inc.

PPV Polymerization through the Gilch Route: Diradical Character of Monomers.

Despite various studies on the polymerization of poly(p-phenylene vinylene) (PPV) through different precursor routes, detailed mechanistic knowledge on the individual reaction steps and intermediates is still incomplete. The present study aims to gain more insight into the radical polymerization of PPV through the Gilch route. The initial steps of the polymerization involve the formation of a p-quinodimethane intermediate, which spontaneously self-initiates through a dimerization process leading to the formation of diradical species; chain propagation ensues on both sides of the diradical or chain termination occurs by the formation of side products, such as [2.2]paracyclophanes. Furthermore, different p-quinodimethane systems were assessed with respect to the size of their aromatic core as well as the presence of heteroatoms in/on the conjugated system. The nature of the aromatic core and the specific substituents alter the electronic structure of the p-quinodimethane monomers, affecting the mechanism of polymerization. The diradical character of the monomers has been investigated with several advanced methodologies, such as spin-projected UHF, CASSCF, CASPT2, and DMRG calculations. It was shown that larger aromatic cores led to a higher diradical character in the monomers, which in turn is proposed to cause rapid initiation.

Possibility of [1,5] sigmatropic shifts in bicyclo[4.2.0]octa-2,4-dienes.

The thermal equilibration of the methyl esters of endiandric acids D and E was subject to a computational study. An electrocyclic pathway via an electrocyclic ring opening followed by a ring flip and a subsequent electrocyclization proposed by Nicolaou [ Nicolaou , K. C. ; Chen , J. S. Chem. Soc. Rev. 2009 , 38 , 2993 ], was computationally explored. The free-energy barrier for this electrocyclic route was shown to be very close to the bicyclo[4.2.0]octa-2,4-diene reported by Huisgen [ Huisgen , R. ; Boche , G. ; Dahmen , A. ; Hechtl , W. Tetrahedron Lett. 1968 , 5215 ]. Furthermore, the possibility of a [1,5] sigmatropic alkyl group shift of bicyclo[4.2.0]octa-2,4-diene systems at high temperatures was explored in a combined computational and experimental study. Calculated reaction barriers for an open-shell singlet biradical-mediated stepwise [1,5] sigmatropic alkyl group shift were shown to be comparable with the reaction barriers for the bicyclo[4.1.0]hepta-2,4-diene (norcaradiene) walk rearrangement. However, the stepwise sigmatropic pathway is suggested to only be feasible for appropriately substituted compounds. Experiments conducted on a deuterated analogous diol derivative confirmed the calculated (large) differences in barriers between electrocyclic and sigmatropic pathways.

First principle chemical kinetics in zeolites: the methanol-to-olefin process as a case study.

To optimally design next generation catalysts a thorough understanding of the chemical phenomena at the molecular scale is a prerequisite. Apart from qualitative knowledge on the reaction mechanism, it is also essential to be able to predict accurate rate constants. Molecular modeling has become a ubiquitous tool within the field of heterogeneous catalysis. Herein, we review current computational procedures to determine chemical kinetics from first principles, thus by using no experimental input and by modeling the catalyst and reacting species at the molecular level. Therefore, we use the methanol-to-olefin (MTO) process as a case study to illustrate the various theoretical concepts. This process is a showcase example where rational design of the catalyst was for a long time performed on the basis of trial and error, due to insufficient knowledge of the mechanism. For theoreticians the MTO process is particularly challenging as the catalyst has an inherent supramolecular nature, for which not only the Brønsted acidic site is important but also organic species, trapped in the zeolite pores, must be essentially present during active catalyst operation. All these aspects give rise to specific challenges for theoretical modeling. It is shown that present computational techniques have matured to a level where accurate enthalpy barriers and rate constants can be predicted for reactions occurring at a single active site. The comparison with experimental data such as apparent kinetic data for well-defined elementary reactions has become feasible as current computational techniques also allow predicting adsorption enthalpies with reasonable accuracy. Real catalysts are truly heterogeneous in a space- and time-like manner. Future theory developments should focus on extending our view towards phenomena occurring at longer length and time scales and integrating information from various scales towards a unified understanding of the catalyst. Within this respect molecular dynamics methods complemented with additional techniques to simulate rare events are now gradually making their entrance within zeolite catalysis. Recent applications have already given a flavor of the benefit of such techniques to simulate chemical reactions in complex molecular environments.

Solved? The reductive radiation chemistry of alanine.

The structural changes throughout the entire reductive radiation-induced pathway of l-α-alanine are solved on an atomistic level with the aid of periodic DFT and nudged elastic band (NEB) simulations. This yields unprecedented information on the conformational changes taking place, including the protonation state of the carboxyl group in the "unstable" and "stable" alanine radicals and the internal transformation converting these two radical variants at temperatures above 220 K. The structures of all stable radicals were verified by calculating EPR properties and comparing those with experimental data. The variation of the energy throughout the full radiochemical process provides crucial insight into the reason why these structural changes and rearrangements occur. Starting from electron capture, the excess electron quickly localizes on the carbon of a carboxyl group, which pyramidalizes and receives a proton from the amino group of a neighboring alanine molecule, forming a first stable radical species (up to 150 K). In the temperature interval 150-220 K, this radical deaminates and deprotonates at the carboxyl group, the detached amino group undergoes inversion and its methyl group sustains an internal rotation. This yields the so-called "unstable alanine radical". Above 220 K, triggered by the attachment of an additional proton on the detached amino group, the radical then undergoes an internal rotation in the reverse direction, giving rise to the "stable alanine radical", which is the final stage in the reductive radiation-induced decay of alanine.

Exploring the vibrational fingerprint of the electronic excitation energy via molecular dynamics.

A Fourier-based method is presented to relate changes of the molecular structure during a molecular dynamics simulation with fluctuations in the electronic excitation energy. The method implies sampling of the ground state potential energy surface. Subsequently, the power spectrum of the velocities is compared with the power spectrum of the excitation energy computed using time-dependent density functional theory. Peaks in both spectra are compared, and motions exhibiting a linear or quadratic behavior can be distinguished. The quadratically active motions are mainly responsible for the changes in the excitation energy and hence cause shifts between the dynamic and static values of the spectral property. Moreover, information about the potential energy surface of various excited states can be obtained. The procedure is illustrated with three case studies. The first electronic excitation is explored in detail and dominant vibrational motions responsible for changes in the excitation energy are identified for ethylene, biphenyl, and hexamethylbenzene. The proposed method is also extended to other low-energy excitations. Finally, the vibrational fingerprint of the excitation energy of a more complex molecule, in particular the azo dye ethyl orange in a water environment, is analyzed.

Synthesis modulation as a tool to increase the catalytic activity of metal-organic frameworks: the unique case of UiO-66(Zr).

The catalytic activity of the zirconium terephthalate UiO-66(Zr) can be drastically increased by using a modulation approach. The combined use of trifluoroacetic acid and HCl during the synthesis results in a highly crystalline material, with partial substitution of terephthalates by trifluoroacetate. Thermal activation of the material leads not only to dehydroxylation of the hexanuclear Zr cluster but also to post-synthetic removal of the trifluoroacetate groups, resulting in a more open framework with a large number of open sites. Consequently, the material is a highly active catalyst for several Lewis acid catalyzed reactions.

Design of zeolite by inverse sigma transformation.

Although the search for new zeolites has traditionally been based on trial and error, more rational methods are now available. The theoretical concept of inverse σ transformation of a zeolite framework to generate a new structure by removal of a layer of framework atoms and contraction has for the first time been achieved experimentally. The reactivity of framework germanium atoms in strong mineral acid was exploited to selectively remove germanium-containing four-ring units from an UTL type germanosilicate zeolite. Annealing of the leached framework through calcination led to the new all-silica COK-14 zeolite with intersecting 12- and 10-membered ring channel systems. An intermediate stage of this inverse σ transformation with dislodged germanate four-rings still residing in the pores could be demonstrated. Inverse σ transformation involving elimination of germanium-containing structural units opens perspectives for the synthesis of many more zeolites.

Enthalpy and entropy barriers explain the effects of topology on the kinetics of zeolite-catalyzed reactions.

The methylation of ethene, propene, and trans-2-butene on zeolites H-ZSM-58 (DDR), H-ZSM-22 (TON), and H-ZSM-5 (MFI) is studied to elucidate the particular influence of topology on the kinetics of zeolite-catalyzed reactions. H-ZSM-58 and H-ZSM-22 are found to display overall lower methylation rates compared to H-ZSM-5 and also different trends in methylation rates with increasing alkene size. These variations may be rationalized based on a decomposition of the free-energy barriers into enthalpic and entropic contributions, which reveals that the lower methylation rates on H-ZSM-58 and H-ZSM-22 have virtually opposite reasons. On H-ZSM-58, the lower methylation rates are caused by higher enthalpy barriers, owing to inefficient stabilization of the reaction intermediates in the large cage-like pores. On the other hand, on H-ZSM-22, the methylation rates mostly suffer from higher entropy barriers, because excessive entropy losses are incurred inside the narrow-channel structure. These results show that the kinetics of crucial elementary steps hinge on the balance between proper stabilization of the reaction intermediates inside the zeolite pores and the resulting entropy losses. These fundamental insights into their inner workings are indispensable for ultimately selecting or designing better zeolite catalysts.

Identification of intermediates in zeolite-catalyzed reactions by in situ UV/Vis microspectroscopy and a complementary set of molecular simulations.

The optical absorption properties of (poly)aromatic hydrocarbons occluded in a nanoporous environment were investigated by theoretical and experimental methods. The carbonaceous species are an essential part of a working catalyst for the methanol-to-olefins (MTO) process. In situ UV/Vis microscopy measurements on methanol conversion over the acidic solid catalysts H-SAPO-34 and H-SSZ-13 revealed the growth of various broad absorption bands around 400, 480, and 580 nm. The cationic nature of the involved species was determined by interaction of ammonia with the methanol-treated samples. To determine which organic species contribute to the various bands, a systematic series of aromatics was analyzed by means of time-dependent density functional theory (TDDFT) calculations. Static gas-phase simulations revealed the influence of structurally different hydrocarbons on the absorption spectra, whereas the influence of the zeolitic framework was examined by using supramolecular models within a quantum mechanics/molecular mechanics framework. To fully understand the origin of the main absorption peaks, a molecular dynamics (MD) study on the organic species trapped in the inorganic host was essential. During such simulation the flexibility is fully taken into account and the effect on the UV/Vis spectra is determined by performing TDDFT calculations on various snapshots of the MD run. This procedure allows an energy absorption scale to be provided and the various absorption bands determined from in situ UV/Vis spectra to be assigned to structurally different species.