The unprecedented ability of computations to probe atomic-level details of catalytic systems holds immense promise for the fundamentals-based bottom-up design of novel heterogeneous catalysts, which are at the heart of the chemical and energy sectors of industry. Here, we critically analyze recent advances in computational heterogeneous catalysis. First, we will survey the progress in electronic structure methods and atomistic catalyst models employed, which have enabled the catalysis community to build increasingly intricate, realistic, and accurate models of the active sites of supported transition-metal catalysts. We then review developments in microkinetic modeling, specifically mean-field microkinetic models and kinetic Monte Carlo simulations, which bridge the gap between nanoscale computational insights and macroscale experimental kinetics data with increasing fidelity. We finally review the advancements in theoretical methods for accelerating catalyst design and discovery. Throughout the review, we provide ample examples of applications, discuss remaining challenges, and provide our outlook for the near future.

Computational Methods in Heterogeneous Catalysis

Reaction kinetics were studied to quantify the effects of polar aprotic organic solvents on the acid-catalyzed conversion of xylose into furfural. A solvent of particular importance is γ-valerolactone (GVL), which leads to significant increases in reaction rates compared to water in addition to increased product selectivity. GVL has similar effects on the kinetics for the dehydration of 1,2-propanediol to propanal and for the hydrolysis of cellobiose to glucose. Based on results obtained for homogeneous Brønsted acid catalysts that span a range of pKa values, we suggest that an aprotic organic solvent affects the reaction kinetics by changing the stabilization of the acidic proton relative to the protonated transition state. This same behavior is displayed by strong solid Brønsted acid catalysts, such as H-mordenite and H-beta.

Solvent Effects in Acid-Catalyzed Biomass Conversion Reactions

Aqueous-phase reactions within microporous Bronsted acids occur at active centers comprised of water-reactant-clustered hydronium ions, solvated within extended hydrogen-bonded water networks that tend to stabilize reactive intermediates and transition states differently The effects of these diverse clustered and networked structures were disentangled here by measuring turnover rates of gas-phase ethanol dehydration to diethyl ether (DEE) on H-form zeolites as water pressure was increased to the point of intrapore condensation, causing protons to become solvated in larger clusters that subsequently become solvated by extended hydrogen-bonded water networks, according to in situ IR spectra Measured first-order rate constants in ethanol quantify the stability of SN2 transition states that eliminate DEE relative to (C2H5OH)(H+)(H2O)n clusters of increasing molecularity, whose structures were respectively determined using metadynamics and ab initio molecular dynamics simulations At low water pressures (2–10 kPa H2O), rate inhibition by water (−1 reaction order) reflects the need to displace one water by ethanol in the cluster en route to the DEE-formation transition state, which resides at the periphery of water–ethanol clusters At higher water pressures (10–75 kPa H2O), water–ethanol clusters reach their maximum stable size ((C2H5OH)(H+)(H2O)4–5), and water begins to form extended hydrogen-bonded networks; concomitantly, rate inhibition by water (up to −3 reaction order) becomes stronger than expected from the molecularity of the reaction, reflecting the more extensive disruption of hydrogen bonds at DEE-formation transition states that contain an additional solvated non-polar ethyl group compared to the relevant reactant cluster, as described by non-ideal thermodynamic formalisms of reaction rates Microporous voids of different hydrophilic binding site density (Beta; varying H+ and Si–OH density) and different size and shape (Beta, MFI, TON, CHA, AEI, FAU), influence the relative extents to which intermediates and transition states disrupt their confined water networks, which manifest as different kinetic orders of inhibition at high water pressures The confinement of water within sub-nanometer spaces influences the structures and dynamics of the complexes and extended networks formed, and in turn their ability to accommodate the evolution in polarity and hydrogen-bonding capacity as reactive intermediates become transition states in Bronsted acid-catalyzed reactions

/pdf/structure-and-solvation-of-confined-water-and-water-ethanol-58b9o9lpn0.pdf

Structure and solvation of confined water and water-ethanol clusters within microporous Brønsted acids and their effects on ethanol dehydration catalysis.

Advances in Understanding the Selective Hydrogenolysis of Biomass Derivatives

Solvent molecules interact with reactive species and alter the rates and selectivities of catalytic reactions by orders of magnitude. Specifically, solvent molecules can modify the free energies of liquid phase and surface species via solvation, participating directly as a reactant or co-catalyst, or competitively binding to active sites. These effects carry consequences for reactions relevant for the conversion of renewable or recyclable feedstocks, the development of distributed chemical manufacturing, and the utilization of renewable energy to drive chemical reactions. First, we describe the quantitative impact of these effects on steady-state catalytic turnover rates through a rate expression derived for a generic catalytic reaction (A → B), which illustrates the functional dependence of rates on each category of solvent interaction. Second, we connect these concepts to recent investigations of the effects of solvents on catalysis to show how interactions between solvent and reactant molecules at solid–liquid interfaces influence catalytic reactions. This discussion demonstrates that the design of effective liquid phase catalytic processes benefits from a clear understanding of these intermolecular interactions and their implications for rates and selectivities.

Influence of solvent structure and hydrogen bonding on catalysis at solid-liquid interfaces.

The composition of the liquid phase can alter the rates of individual reaction steps and thus alter the selectivity of acid-catalyzed reactions, but these solvent effects are difficult to anticipate for design purposes. Herein, we report the kinetics and selectivity of Bronsted acid-catalyzed 1,2-propanediol dehydration in pure water and in aqueous mixtures of the polar aprotic cosolvents γ-valerolactone, 1,4-dioxane, tetrahydrofuran, N-methyl-2-pyrrolidone, tetramethylene sulfoxide, and dimethyl sulfoxide at 433 K. We find that the major product of 1,2-propanediol dehydration is propanal in most mixed-solvent environments with selectivities between 1 and 68 mol %. In contrast, 1,2-propanediol dehydration in aqueous mixtures of dimethyl sulfoxide affords acetone as the major product with up to 48% selectivity with minimal propanal formation. We use classical molecular dynamics simulations to probe these solvent effects by computing the difference between the solvation free energies of 1,2-propanediol and ...

Liam Witteman

Papers

Effect of Mixed-Solvent Environments on the Selectivity of Acid-Catalyzed Dehydration Reactions