Veronique Van Speybroeck

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.

Reproducibility in density functional theory calculations of solids

A comparison of DFT methods Density functional theory (DFT) is now routinely used for simulating material properties. Many software packages are available, which makes it challenging to know which are the best to use for a specific calculation. Lejaeghere et al. compared the calculated values for the equation of states for 71 elemental crystals from 15 different widely used DFT codes employing 40 different potentials (see the Perspective by Skylaris). Although there were variations in the calculated values, most recent codes and methods converged toward a single value, with errors comparable to those of experiment. Science, this issue p. 10.1126/science.aad3000; see also p. 1394 A survey of recent density functional theory methods shows a convergence to more accurate property calculations. [Also see Perspective by Skylaris] INTRODUCTION The reproducibility of results is one of the underlying principles of science. An observation can only be accepted by the scientific community when it can be confirmed by independent studies. However, reproducibility does not come easily. Recent works have painfully exposed cases where previous conclusions were not upheld. The scrutiny of the scientific community has also turned to research involving computer programs, finding that reproducibility depends more strongly on implementation than commonly thought. These problems are especially relevant for property predictions of crystals and molecules, which hinge on precise computer implementations of the governing equation of quantum physics. RATIONALE This work focuses on density functional theory (DFT), a particularly popular quantum method for both academic and industrial applications. More than 15,000 DFT papers are published each year, and DFT is now increasingly used in an automated fashion to build large databases or apply multiscale techniques with limited human supervision. Therefore, the reproducibility of DFT results underlies the scientific credibility of a substantial fraction of current work in the natural and engineering sciences. A plethora of DFT computer codes are available, many of them differing considerably in their details of implementation, and each yielding a certain “precision” relative to other codes. How is one to decide for more than a few simple cases which code predicts the correct result, and which does not? We devised a procedure to assess the precision of DFT methods and used this to demonstrate reproducibility among many of the most widely used DFT codes. The essential part of this assessment is a pairwise comparison of a wide range of methods with respect to their predictions of the equations of state of the elemental crystals. This effort required the combined expertise of a large group of code developers and expert users. RESULTS We calculated equation-of-state data for four classes of DFT implementations, totaling 40 methods. Most codes agree very well, with pairwise differences that are comparable to those between different high-precision experiments. Even in the case of pseudization approaches, which largely depend on the atomic potentials used, a similar precision can be obtained as when using the full potential. The remaining deviations are due to subtle effects, such as specific numerical implementations or the treatment of relativistic terms. CONCLUSION Our work demonstrates that the precision of DFT implementations can be determined, even in the absence of one absolute reference code. Although this was not the case 5 to 10 years ago, most of the commonly used codes and methods are now found to predict essentially identical results. The established precision of DFT codes not only ensures the reproducibility of DFT predictions but also puts several past and future developments on a firmer footing. Any newly developed methodology can now be tested against the benchmark to verify whether it reaches the same level of precision. New DFT applications can be shown to have used a sufficiently precise method. Moreover, high-precision DFT calculations are essential for developing improvements to DFT methodology, such as new density functionals, which may further increase the predictive power of the simulations. Recent DFT methods yield reproducible results. Whereas older DFT implementations predict different values (red darts), codes have now evolved to mutual agreement (green darts). The scoreboard illustrates the good pairwise agreement of four classes of DFT implementations (horizontal direction) with all-electron results (vertical direction). Each number reflects the average difference between the equations of state for a given pair of methods, with the green-to-red color scheme showing the range from the best to the poorest agreement. The widespread popularity of density functional theory has given rise to an extensive range of dedicated codes for predicting molecular and crystalline properties. However, each code implements the formalism in a different way, raising questions about the reproducibility of such predictions. We report the results of a community-wide effort that compared 15 solid-state codes, using 40 different potentials or basis set types, to assess the quality of the Perdew-Burke-Ernzerhof equations of state for 71 elemental crystals. We conclude that predictions from recent codes and pseudopotentials agree very well, with pairwise differences that are comparable to those between different high-precision experiments. Older methods, however, have less precise agreement. Our benchmark provides a framework for users and developers to document the precision of new applications and methodological improvements.

Electronic effects of linker substitution on Lewis acid catalysis with metal-organic frameworks.

Functionalized linkers can greatly increase the activity of metal-organic framework (MOF) catalysts with coordinatively unsaturated sites. A clear linear free-energy relationship (LFER) was found between Hammett σ(m) values of the linker substituents X and the rate k(X) of a carbonyl-ene reaction. This is the first LFER ever observed for MOF catalysts. A 56-fold increase in rate was found when the substituent is a nitro group.

Unraveling the reaction mechanisms governing methanol-to-olefins catalysis by theory and experiment

The conversion of methanol to olefins (MTO) over a heterogeneous nanoporous catalyst material is a highly complex process involving a cascade of elementary reactions. The elucidation of the reaction mechanisms leading to either the desired production of ethene and/or propene or undesired deactivation has challenged researchers for many decades. Clearly, catalyst choice, in particular topology and acidity, as well as the specific process conditions determine the overall MTO activity and selectivity; however, the subtle balances between these factors remain not fully understood. In this review, an overview of proposed reaction mechanisms for the MTO process is given, focusing on the archetypal MTO catalysts, H-ZSM-5 and H-SAPO-34. The presence of organic species, that is, the so-called hydrocarbon pool, in the inorganic framework forms the starting point for the majority of the mechanistic routes. The combination of theory and experiment enables a detailed description of reaction mechanisms and corresponding reaction intermediates. The identification of such intermediates occurs by different spectroscopic techniques, for which theory and experiment also complement each other. Depending on the catalyst topology, reaction mechanisms proposed thus far involve aromatic or aliphatic intermediates. Ab initio simulations taking into account the zeolitic environment can nowadays be used to obtain reliable reaction barriers and chemical kinetics of individual reactions. As a result, computational chemistry and by extension computational spectroscopy have matured to the level at which reliable theoretical data can be obtained, supplying information that is very hard to acquire experimentally. Special emphasis is given to theoretical developments that open new perspectives and possibilities that aid to unravel a process as complex as methanol conversion over an acidic porous material.

Error Estimates for Solid-State Density-Functional Theory Predictions: An Overview by Means of the Ground-State Elemental Crystals

Predictions of observable properties by density-functional theory calculations (DFT) are used increasingly often by experimental condensed-matter physicists and materials engineers as data. These predictions are used to analyze recent measurements, or to plan future experiments in a rational way. Increasingly more experimental scientists in these fields therefore face the natural question: what is the expected error for such a first-principles prediction? Information and experience about this question is implicitly available in the computational community, scattered over two decades of literature. The present review aims to summarize and quantify this implicit knowledge. This eventually leads to a practical protocol that allows any scientist—experimental or theoretical—to determine justifiable error estimates for many basic property predictions, without having to perform additional DFT calculations. A central role is played by a large and diverse test set of crystalline solids, containing all ground-state elemental crystals (except most lanthanides). For several properties of each crystal, the difference between DFT results and experimental values is assessed. We discuss trends in these deviations and review explanations suggested in the literature. A prerequisite for such an error analysis is that different implementations of the same first-principles formalism provide the same predictions. Therefore, the reproducibility of predictions across several mainstream methods and codes is discussed too. A quality factor Δ expresses the spread in predictions from two distinct DFT implementations by a single number. To compare the PAW method to the highly accurate APW+lo approach, a code assessment of VASP and GPAW (PAW) with respect to WIEN2k (APW+lo) yields Δ-values of 1.9 and 3.3 meV/atom, respectively. In both cases the PAW potentials recommended by the respective codes have been used. These differences are an order of magnitude smaller than the typical difference with experiment, and therefore predictions by APW+lo and PAW are for practical purposes identical.

First Principle Kinetic Studies of Zeolite-Catalyzed Methylation Reactions

Methylations of ethene, propene, and butene by methanol over the acidic microporous H-ZSM-5 catalyst are studied by means of state of the art computational techniques, to derive Arrhenius plots and rate constants from first principles that can directly be compared with the experimental data. For these key elementary reactions in the methanol to hydrocarbons (MTH) process, direct kinetic data became available only recently [J. Catal.2005, 224, 115-123; J. Catal.2005, 234, 385-400]. At 350 °C, apparent activation energies of 103, 69, and 45 kJ/mol and rate constants of 2.6 × 10(-4), 4.5 × 10(-3), and 1.3 × 10(-2) mol/(g h mbar) for ethene, propene, and butene were derived, giving following relative ratios for methylation k(ethene)/k(propene)/k(butene) = 1:17:50. In this work, rate constants including pre-exponential factors are calculated which give very good agreement with the experimental data: apparent activation energies of 94, 62, and 37 kJ/mol for ethene, propene, and butene are found, and relative ratios of methylation k(ethene)/k(propene)/k(butene) = 1:23:763. The entropies of gas phase alkenes are underestimated in the harmonic oscillator approximation due to the occurrence of internal rotations. These low vibrational modes were substituted by manually constructed partition functions. Overall, the absolute reaction rates can be calculated with near chemical accuracy, and qualitative trends are very well reproduced. In addition, the proposed scheme is computationally very efficient and constitutes significant progress in kinetic modeling of reactions in heterogeneous catalysis.

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.

The remarkable catalytic activity of the saturated metal organic framework V-MIL-47 in the cyclohexene oxidation

The remarkable catalytic activity of the saturated metal organic framework MIL-47 in the epoxidation of cyclohexene is elucidated by means of both experimental results and theoretical calculations.

Theoretical Insights on Methylbenzene Side-Chain Growth in ZSM-5 Zeolites for Methanol-to-Olefin Conversion

The key step in the conversion of methane to polyolefins is the catalytic conversion of methanol to light olefins. The most recent formulations of a reaction mechanism for this process are based on the idea of a complex hydrocarbon-pool network, in which certain organic species in the zeolite pores are methylated and from which light olefins are eliminated. Two major mechanisms have been proposed to date-the paring mechanism and the side-chain mechanism-recently joined by a third, the alkene mechanism. Recently we succeeded in simulating a full catalytic cycle for the first of these in ZSM-5, with inclusion of the zeolite framework and contents. In this paper, we will investigate crucial reaction steps of the second proposal (the side-chain route) using both small and large zeolite cluster models of ZSM-5. The deprotonation step, which forms an exocyclic double bond, depends crucially on the number and positioning of the other methyl groups but also on steric effects that are typical for the zeolite lattice. Because of steric considerations, we find exocyclic bond formation in the ortho position to the geminal methyl group to be more favourable than exocyclic bond formation in the para position. The side-chain growth proceeds relatively easily but the major bottleneck is identified as subsequent de-alkylation to produce ethene. These results suggest that the current formulation of the side-chain route in ZSM-5 may actually be a deactivating route to coke precursors rather than an active ethene-producing hydrocarbon-pool route. Other routes may be operating in alternative zeotype materials like the silico-aluminophosphate SAPO-34.

Active site engineering in UiO-66 type metal–organic frameworks by intentional creation of defects: a theoretical rationalization

The catalytic activity of the Zr-benzenedicarboxylate (Zr-BDC) UiO-66 can be drastically increased if some BDC linkers are missing, as this removes the full coordination of the framework metal ions. As a result, metal centers become more accessible and thus more active for Lewis acid catalysed reactions. Addition of modulators (MDL) to the synthesis mixture can create more linker deficiencies (Vermoortele et al., J. Am. Chem. Soc., 2013, 135, 11465) and leads to a significant increase in the catalytic activity due to the creation of a larger number of open sites. In this paper, we rationalize the function of the modulators under real synthesis conditions by the construction of free energy diagrams. The UiO-66 type materials form a very appropriate test case as the effect of addition of modulators hydrochloric acid (HCl) and trifluoroacetate (TFA) has been intensively investigated experimentally for the synthesis process and post-synthetic thermal activation. Under synthesis conditions, direct removal of BDC linkers requires a high free energy, but replacement of such linker by one or more TFA species might occur especially at high TFA : BDC ratios in the reaction mixture. Post-synthesis activation procedures at higher temperatures lead to substantial removal of the species coordinated to the Zr bricks, creating open metal sites. A mechanistic pathway is presented for the dehydroxylation process of the hexanuclear Zr cluster. For the citronellal cyclization, we show that the presence of some residual TFA in the structure may lead to faster reactions in complete agreement with the experiment. Hirshfeld-e partial charges for the Zr ions have been computed to investigate their sensitivity to substituent effects; a strong correlation with the experimental Hammett parameters and with the rates of the citronellal cyclization is found. The theoretical rationalization may serve as a basis for detailed active site engineering studies.