The direct transformation of CO2 into high-value-added hydrocarbons (i.e., olefins and aromatics) has the potential to make a decisive impact in our society. However, despite the efforts of the ...scientific community, no direct synthetic route exists today to synthesize olefins and aromatics from CO2 with high productivities and low undesired CO selectivity. Herein, we report the combination of a series of catalysts comprising potassium superoxide doped iron oxide and a highly acidic zeolite (ZSM-5 and MOR) that directly convert CO2 to either light olefins (in MOR) or aromatics (in ZSM-5) with high space–time yields (STYC2‑C4= = 11.4 mmol·g–1·h–1; STYAROM = 9.2 mmol·g–1·h–1) at CO selectivities as low as 12.8% and a CO2 conversion of 49.8% (reaction conditions: T = 375 °C, P = 30 bar, H2/CO2 = 3, and 5000 mL·g–1·h–1). Comprehensive solid-state nuclear magnetic resonance characterization of the zeolite component reveals that the key for the low CO selectivity is the formation of surface formate species on the zeolite framework. The remarkable difference in selectivity between the two zeolites is further rationalized by first-principles simulations, which show a difference in reactivity for crucial carbenium ion intermediates in MOR and ZSM-5.
The mobility of the copper cations acting as active sites for the selective catalytic reduction of nitrogen oxides with ammonia in Cu-CHA catalysts varies with temperature and feed composition. ...Herein, the migration of Cu(NH3)2+ complexes between two adjacent cavities of the chabazite structure, including other reactant molecules (NO, O2, H2O, and NH3), in the initial and final cavities is investigated using ab initio molecular dynamics (AIMD) simulations combined with enhanced sampling techniques to describe hopping events from one cage to the other. We find that such diffusion is only significantly hindered by the presence of excess NH3 or NO in the initial cavity, since both reactants form with Cu(NH3)2+ stable intermediates which are too bulky to cross the 8-ring windows connecting the cavities. The presence of O2 modifies strongly the interaction of NO with Cu+. At low temperatures, we observe NO detachment from Cu+ and increased mobility of the Cu(NH3)2+ complex, while at high temperatures, NO reacts spontaneously with O2 to form NO2. The present simulations give evidence for recent experimental observations, namely, an NH3 inhibition effect on the SCR reaction at low temperatures, and transport limitations of NO and NH3 at high temperatures. Our first principle simulations mimicking operating conditions support the existence of two different reaction mechanisms operating at low and high temperatures, the former involving dimeric Cu(NH3)2-O2-Cu(NH3)2 species and the latter occurring by direct NO oxidation to NO2 in one single cavity.
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.
Within this Perspective, we critically reflect on the role of first-principles molecular dynamics (MD) simulations in unraveling the catalytic function within zeolites under operating conditions. ...First-principles MD simulations refer to methods where the dynamics of the nuclei is followed in time by integrating the Newtonian equations of motion on a potential energy surface that is determined by solving the quantum-mechanical many-body problem for the electrons. Catalytic solids used in industrial applications show an intriguing high degree of complexity, with phenomena taking place at a broad range of length and time scales. Additionally, the state and function of a catalyst critically depend on the operating conditions, such as temperature, moisture, presence of water, etc. Herein we show by means of a series of exemplary cases how first-principles MD simulations are instrumental to unravel the catalyst complexity at the molecular scale. Examples show how the nature of reactive species at higher catalytic temperatures may drastically change compared to species at lower temperatures and how the nature of active sites may dynamically change upon exposure to water. To simulate rare events, first-principles MD simulations need to be used in combination with enhanced sampling techniques to efficiently sample low-probability regions of phase space. Using these techniques, it is shown how competitive pathways at operating conditions can be discovered and how broad transition state regions can be explored. Interestingly, such simulations can also be used to study hindered diffusion under operating conditions. The cases shown clearly illustrate how first-principles MD simulations reveal insights into the catalytic function at operating conditions, which could not be discovered using static or local approaches where only a few points are considered on the potential energy surface (PES). Despite these advantages, some major hurdles still exist to fully integrate first-principles MD methods in a standard computational catalytic workflow or to use the output of MD simulations as input for multiple length/time scale methods that aim to bridge to the reactor scale. First of all, methods are needed that allow us to evaluate the interatomic forces with quantum-mechanical accuracy, albeit at a much lower computational cost compared to currently used density functional theory (DFT) methods. The use of DFT limits the currently attainable length/time scales to hundreds of picoseconds and a few nanometers, which are much smaller than realistic catalyst particle dimensions and time scales encountered in the catalysis process. One solution could be to construct machine learning potentials (MLPs), where a numerical potential is derived from underlying quantum-mechanical data, which could be used in subsequent MD simulations. As such, much longer length and time scales could be reached; however, quite some research is still necessary to construct MLPs for the complex systems encountered in industrially used catalysts. Second, most currently used enhanced sampling techniques in catalysis make use of collective variables (CVs), which are mostly determined based on chemical intuition. To explore complex reactive networks with MD simulations, methods are needed that allow the automatic discovery of CVs or methods that do not rely on a priori definition of CVs. Recently, various data-driven methods have been proposed, which could be explored for complex catalytic systems. Lastly, first-principles MD methods are currently mostly used to investigate local reactive events. We hope that with the rise of data-driven methods and more efficient methods to describe the PES, first-principles MD methods will in the future also be able to describe longer length/time scale processes in catalysis. This might lead to a consistent dynamic description of all stepsdiffusion, adsorption, and reactionas they take place at the catalyst particle level.
Display omitted
•Enhanced sampling methods are used to study the methylation of ethene in H-ZSM-5.•Transformations of free energy profiles between collective variables were validated.•Proper reaction ...rates independent of the collective variable were obtained.•Near chemical accuracy was obtained with the proposed protocol.
The methylation reaction of ethene with methanol over the Brønsted acidic ZSM-5 catalyst is one of the prototype reactions within zeolite catalysis for which experimental kinetic data is available. It is one of the premier reactions within the methanol-to-olefins process and has been the subject of extensive theoretical testing to predict the reaction rates. Herein, we apply, for the first time, first principle molecular dynamics methods to determine the intrinsic reaction kinetics taking into account the full configurational entropy. As chemical reactions are rare events, enhanced sampling methods are necessary to obtain sufficient sampling of the configurational space at the activated region. A plethora of methods is available which depend on specific choices like the selection of collective variables along which the dynamics is enhanced. Herein, a thorough first principle molecular dynamics study is presented to determine the reaction kinetics via various enhanced MD techniques on an exemplary reaction within zeolite catalysis for which reference theoretical and experimental data are available.
Catalytic alkene cracking on H-ZSM-5 involves a complex reaction network with many possible reaction routes and often elusive intermediates. Herein, advanced molecular dynamics simulations at 773 K, ...a typical cracking temperature, are performed to clarify the nature of the intermediates and to elucidate dominant cracking pathways at operating conditions. A series of C4–C8 alkene intermediates are investigated to evaluate the influence of chain length and degree of branching on their stability. Our simulations reveal that linear, secondary carbenium ions are relatively unstable, although their lifetime increases with carbon number. Tertiary carbenium ions, on the other hand, are shown to be very stable, irrespective of the chain length. Highly branched carbenium ions, though, tend to rapidly rearrange into more stable cationic species, either via cracking or isomerization reactions. Dominant cracking pathways were determined by combining these insights on carbenium ion stability with intrinsic free energy barriers for various octene β-scission reactions, determined via umbrella sampling simulations at operating temperature (773 K). Cracking modes A (3° → 3°) and B2 (3° → 2°) are expected to be dominant at operating conditions, whereas modes B1 (2° → 3°), C (2° → 2°), D2 (2° → 1°), and E2 (3° → 1°) are expected to be less important. All β-scission modes in which a transition state with primary carbocation character is involved have high intrinsic free energy barriers. Reactions starting from secondary carbenium ions will contribute less as these intermediates are short living at the high cracking temperature. Our results show the importance of simulations at operating conditions to properly evaluate the carbenium ion stability for β-scission reactions and to assess the mobility of all species in the pores of the zeolite.
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.
The influence of the acid site density on the diffusivities of alkanes and alkenes through the 8‐ring windows of SAPO‐34 was probed using a complementary set of ab initio MD simulations and experimental techniques (PFG‐NMR and TAP pulse‐response measurements). The presence of Brønsted acid sites was found to have a clear promotional effect on alkene diffusion but it does not influence alkane diffusion.
Display omitted
•Acidity can be tuned by isomorphic substitution of H-MeAlPO-5 materials.•Static DFT calculations indicate a linear correlation between benzene methylation barriers and acid ...strength.•Experimental results corroborate computational trends with Mg-, Si- and ZrAlPO-5 catalysts.•Benzene methylation is hindered at high temperature and DME partial pressure.•MgAlPO-5 promotes a selective methylation of alkenes compared to SAPO-5 counterpart.
Methylation of aromatic compounds is a key reaction step in various industrial processes such as the aromatic cycle of methanol-to-hydrocarbons chemistry. The study of isolated methylation reactions and of the influence of catalyst acidity on their kinetics is a challenging task. Herein, we have studied unidirectional metal-substituted H-MeAlPO-5 materials to evaluate the effect of acid strength on the kinetics of benzene methylation with DME. First-principle simulations showed a direct correlation between the methylation barrier and acid site strength, which depends on the metal substituent. Three H-MeAlPO-5 catalysts with high (Me = Mg), moderate (Me = Si) and low acidity (Me = Zr) were experimentally tested, confirming a linear relationship between the methylation activation energy and acid strength. The effects of temperature and reactant partial pressure were evaluated, showing significant differences in the byproduct distribution between H-MgAlPO-5 and H-SAPO-5. Comparison with propene methylation suggested that the Mg substituted catalyst is also the most active for the selective methylation of alkenes.