Search NASA⌕ Search

SEARCH · Search NASA

Results for “Microkinetic Modeling”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 91 records · Page 5
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Effects of surface diffusion in electrocatalytic CO 2 reduction on Cu revealed by kinetic Monte Carlo simulations

We report Kinetic Monte Carlo (KMC) methods are frequently used for mechanistic studies of thermally driven heterogeneous catalysis systems but are underused for electrocatalysis. Here, we develop a lattice KMC approach for electrocatalytic CO 2 reduction. The work is motivated by a prior experimental report that performed electroreduction of a mixed feed of 12 CO 2 and 13 CO on Cu; differences in the 13 C content of C2 products ethylene and ethanol (Δ 13 C) were interpreted as evidence of site selectivity. The lattice KMC model considers the effect of surface diffusion on this system. In the limit of infinitely fast diffusion (mean-field approximation), the key intermediates 12 CO* and 13 CO* would be well mixed on the surface and no evidence of site selectivity could have been observed. Using a simple two-site model and adapting a previously reported microkinetic model, we assess the effects of diffusion on the relative isotope fractions in the products using the estimated surface diffusion rate of CO* from literature reports. We find that the size of the active sites and the total surface adsorbate coverage can have a large influence on the values of Δ 13 C that can be observed. Δ 13 C is less sensitive to the CO* diffusion rate as long as it is within the estimated range. We further offer possible methods to estimate surface distribution of intermediates and to predict intrinsic selectivity of active sites based on experimental observations. This work illustrates the importance of considering surface diffusion in the study of electrochemical CO 2 reduction to multi-carbon products. Our approach is entirely based on a freely available open-source code, so will be readily adaptable to other electrocatalytic systems.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2024

Investigation of Ethane Dehydrogenation and Hydrogenolysis on Pt(111), Pt(211), and Pt(100): Bayesian Quantification and Correction of DFT-Based Enthalpic and Entropic Uncertainties

Computational investigations of heterogeneously catalyzed reactions using density functional theory (DFT) are often inaccurate, largely due to uncertainties in the choice of DFT functional (enthalpic uncertainty) and approximations for modeling adsorbate movement along the catalyst surface (entropic uncertainty). This work illustrates that both uncertainties are significant in the investigation of ethane dehydrogenation (EDH) and hydrogenolysis on Pt catalysts by considering the complete deconstruction of ethane on Pt(111), Pt(211), and Pt(100) using microkinetic modeling (MKM). Hence, this work uses both noncalibrated and Bayesian-calibrated MKMs to quantify and correct inaccuracies in macroscopic properties due to both uncertainties. A Bayesian approach to the correction of entropic errors was introduced using a “Modified Fermi Function (MFF)” to calibrate between the two bounds of entropy represented by the harmonic oscillator (HO) and free translator (FT) approximations. Regardless of enthalpic and entropic uncertainties, all three surfaces are capable of ethane activation; however, Pt(211) was found to be the most active and is largely responsible for methane production. Next, Pt(111) is largely responsible for acetylene production, and Pt(100) has the highest ethylene selectivity but is most susceptible to coking. By comparison of different calibrated models, the FT entropy approximation was found to better describe EDH under typical experimental conditions. Statistical evidence was found to support Pt(111) as the active site for EDH, assuming that one single site is responsible for the chemistry. On the three surfaces, competing second dehydrogenations to CH 2 CH 2 and CH 3 CH were observed as well as isomerization of CH 3 CH back to CH 2 CH 2 and deeper dehydrogenation of CH 3 CH. In conclusion, C–C cleavage was found to largely proceed via the CH 3 C intermediate on Pt(100) and Pt(111), while on Pt(211), it was via both CHC and CH 3 C.

Bayesian model selection↗
DOE OSTI2020

Kinetic Isotope Effects: Interpretation and Prediction Using Degrees of Rate Control

Kinetic isotope effects (KIEs) have been used for decades in catalysis research as a tool for clarifying reaction mechanisms. Significant primary kinetic isotope effects have usually been interpreted as being a result of isotope substitution at a site of bond breaking (or forming) in the rate-determining step in the reaction mechanism. However, quantitative analysis of the magnitude of the KIE in complex multistep reaction mechanisms is seldom reported. We prove here that the logarithm of the rate ratio for two isotopes is the weighted average over all species in the mechanism of their standard-state free energy difference between the two isotopes, divided by RT. The weighting factor is the degree of rate control (DRC) for that species (e.g., transition state, intermediate, reactant) when the rate is measured separately for each isotope. It is instead the degree of selectivity ratio control (DSRC) when the KIE is measured as the product selectivity in a parallel competition between two isotopes within the same reactant molecule. Since only a few species have non-zero DRCs (or DSRCs) for most reactions, only these few contribute to this weighted average and the KIE. We show that this provides a simple way to interpret and quantitatively predict kinetic isotope effects that is powerful in the insights it provides, allowing one to evaluate directly which species contribute most to the KIE. By applying it to H/D KIEs in several example mechanisms, we further show that the traditional way of interpreting KIEs that focuses only on the rate-determining step can easily lead to misunderstanding of KIE and the reaction mechanism. Finally, this highlights the importance to consider the effect of isotope substitution on the energies of all species with large DRCs (i.e, those whose energies are kinetically relevant). This method also offers opportunities for quantitative validation of mechanism-based microkinetic models.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2022

Modified Energy Span Analysis Reveals Heterogeneous Catalytic Kinetics

Mechanistic modeling is a cornerstone of catalyst development generally conducted with microkinetic models or density functional theory-based energy profiles. We extend the energy span model of homogeneous catalysis to heterogeneous systems by introducing the modified energy span analysis (MESA) model by implementing collision theory and gas-phase concentration effects. We determine analytically turnover frequencies, coverages, rate-determining steps, apparent activation energies, and reaction orders in agreement with microkinetic and kinetic Monte Carlo simulations. The model applies to discrete catalyst state systems, including single-atom catalysts, homogeneous systems, and microporous materials. A generalizable rate expression is derived for reactor modeling or mechanistic insights from solely experimentally measured reaction orders. Here, we illustrate MESA on three published mechanisms and reveal unexpected phenomena absent in traditional energy span profiles.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2023

Effect of Co/SiO 2 Single-Site Heterogeneity on Small Alkane Dehydrogenation Kinetics

Atomically dispersed, high-spin Co(II) atoms in distorted tetrahedral coordination to an amorphous silica (am-SiO 2 ) support, and more recently in zeolite frameworks, are active and selective for light alkane dehydrogenation. This paper investigates how variations in the geometry of the active sites affect the ethane dehydrogenation activity of atomically dispersed Co(II) on an am-SiO 2 support. We generate a distribution of sites and determine the geometric parameters that exhibit the strongest correlation with the coordination geometry and activity of the metal atom by means of linear dimensionality reduction techniques. We perform electronic structure calculations and microkinetic modeling and deduce the mechanism and kinetics for a representative sample of sites. Irrespective of the active site geometry, the rate of ethane dehydrogenation is governed by the β-hydride elimination, which involves quartet-doublet spin-crossing and proceeds adiabatically due to strong spin-orbit coupling. Informed by the complete microkinetic analysis of the sites, we derive the reduced rate expression as a function of three site-dependent quantities. We show that these site-dependent quantities correlate with the energy of formation of the ethyl intermediate that forms via C-H bond activation. This correlation allows us to derive the site-averaged rate for the entire distribution of sites. Among various sites, the tri-coordinate and planar tetra-coordinate Co sites exhibit higher Lewis acidity than the tetrahedral sites, and consequently, higher initial rates. Finally, we discuss the implications for more active catalysts.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2021

Activation of N 2 on Manganese Nitride-Supported Ni 3 and Fe 3 Clusters and Relevance to Ammonia Formation

Dual-site models were constructed to represent manganese nitride (Mn 4 N)-supported Ni 3 and Fe 3 clusters for NH 3 synthesis. Density functional theory calculations produced an energy barrier of approximately 0.55 eV for N–N bond activation at the interfacial nitrogen vacancy sites (N v ); also, the hydrogenation and removal of interfacial N is promoted by earth-abundant Ni and Fe metals. Steady-state microkinetic modeling revealed that the turnover frequencies of NH 3 production follow an order of Fe 3 @Mn 4 N ≈ Ni 3 @Mn 4 N > Mn 4 N > Fe >> Ni. Moreover, we present clear evidence that, before NH 3 formation, NH migrates from N v onto the metallic sites. Using N binding energy (BE N ) and the transition-state energy of N 2 activation (E TS ) as descriptors, we concluded that the beneficial effects owing to interfacial N v sites are the most pronounced when BE N is either too strong or too weak while E TS is high; otherwise, excessive N v sites may hinder catalyst performance.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2018

Computational Design of Active Site Structures with Improved Transition-State Scaling for Ammonia Synthesis

Here, the Haber–Bosch process for the reduction of atmospheric nitrogen to ammonia is one of the most optimized heterogeneous catalytic reactions, but there are aspects of the industrial process that remain less than ideal. It has been shown that the activity of metal catalysts is limited by a Brønsted–Evans–Polanyi (BEP) scaling relationship between the reaction and transition-state energies for N 2 dissociation, leading to a negligible production rate at ambient conditions and a modest rate under harsh conditions. In this study, we use density functional theory (DFT) calculations in conjunction with mean-field microkinetic modeling to study the rate of NH 3 synthesis on model active sites that require the singly coordinated dissociative adsorption of N atoms onto transition metal atoms. Our results demonstrate that this ”on-top” binding of nitrogen exhibits significantly improved scaling behavior, which can be rationalized in terms of transition-state geometries and leads to considerably higher predicted activity. While synthesis of these model systems is likely challenging, the stabilization of such an active site could enable thermochemical ammonia synthesis under more benign conditions.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2019

pH Effects on Hydrogen Evolution and Oxidation over Pt(111): Insights from First-Principles

A long-standing question in electrocatalysis is how the electrolyte pH affects the catalytic activity of proton–electron transfer reactions. A prime example for this is the hydrogen evolution/oxidation reaction (HER/HOR) over metal catalysts. While it has long been established that alkaline conditions result in a more sluggish reaction kinetics than in acidic conditions, the underlying reason for this trend remains contentious. Here, we apply density functional calculations to evaluate prevailing hypotheses for the origin of this effect: shifts in hydrogen binding, proton donor, and water reorganization energy. We present a microkinetic model, based on ab initio reaction energetics of all possible elementary steps. Our model shows a good agreement with experimental trends. We find that with increasing pH, the proton donor changes from hydronium to water. In conclusion, our model suggests that the intrinsically larger barriers for the splitting of water with respect to hydronium are the cause of HER kinetics being slower in alkaline than in acidic media.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2024

Machine Learning Accelerated First-Principles Study of the Hydrodeoxygenation of Propanoic Acid

The complex reaction network of catalytic biomass conversions often involves hundreds of surface intermediates and thousands of reaction steps, greatly hindering the rational design of metal catalysts for these conversions. Here, we present a framework of machine learning (ML)-accelerated first-principles studies for the hydrodeoxygenation (HDO) of propanoic acid over transition metal surfaces. The microkinetic model (MKM) is initially parametrized by ML-predicted energies and iteratively improved by identifying the rate-determining species and steps (RDS), computing their energies by density functional theory (DFT), and reparameterizing the MKM until all the RDS are computed by DFT. The Gaussian process (GP) model performs significantly better than the linear ridge regression model for predicting both the adsorption free energies and transition state free energies. Parameterized with energies from the GP model, only 5–20% of the full reaction network has to be computed by DFT for the MKM to possess DFT-level accuracy for the TOF and dominant reaction pathway. While the linear ridge regression model performs worse than the GP model, its performance is greatly improved when only transition states are predicted by the regression model and adsorption energies are computed by DFT. Overall, we find that a high accuracy in adsorption free energies is more important for a reliable MKM than a high accuracy in TS free energies. Lastly, based on the GP model with GOH and GCHCHCO as catalyst descriptors, we build two-dimensional volcano plots in activity and selectivity that can help design promising alloy catalysts for HDO reactions of organic acids.

adsorption↗
DOE OSTI2018

Selectivity of Synthesis Gas Conversion to C 2+ Oxygenates on fcc(111) Transition-Metal Surfaces

Using a combined density functional theory and descriptor based microkinetic model approach, we predict production rate volcanos for higher oxygenate formation on (111) transition-metal surfaces. Despite their lower activity for CO conversion compared to stepped surfaces, (111) transition metal surfaces bring the potential for selectivity toward C 2+ oxygenates. As a result, the volcano plots can be used to rationalize and predict activity and selectivity trends for transition-metal-based catalysts.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2021

Achieving industrial ammonia synthesis rates at near-ambient conditions through modified scaling relations on a confined dual site

The production of ammonia through the Haber–Bosch process is regarded as one of the most important inventions of the 20th century. Despite significant efforts in optimizing the process, it still consumes 1 to 2% of the worldwide annual energy for the high working temperatures and pressures. The design of a catalyst with a high activity at milder conditions represents another challenge for this reaction. Herein, we combine density functional theory and microkinetic modeling to illustrate a strategy to facilitate low-temperature and -pressure ammonia synthesis through modified energy-scaling relationships using a confined dual site. Our results suggest that an ammonia synthesis rate two to three orders of magnitude higher than the commercial Ru catalyst can be achieved under the same reaction conditions with the introduction of confinement. Such strategies will open pathways for the development of catalysts for the Haber–Bosch process that can operate at milder conditions and present more economically viable alternatives to current industrial solutions.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2017

Titania-Supported Single-Atom Platinum Catalyst for Water-Gas Shift Reaction

Density functional theory (DFT) and microkinetic modeling are used to examine the water-gas shift (WGS) reaction on TiO 2 (110)-supported single Pt atoms. A CO-adsorbed Pt atom supported on a TiO 2 surface exhibits a lower activity for the WGS than the Pt 2+ ion supported on a doped TiO 2 surface, although DFT predicts the CO-adsorbed Pt on an undoped TiO 2 surface to be thermodynamically more stable. Finally, trends observed in WGS activity for different active sites of TiO 2 -supported Pt catalysts can be explained in terms of CO adsorption strength and the ability to form an oxygen vacancy on the oxide support.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2020

DFT and microkinetic comparison of Pt, Pd and Rh-catalyzed ammonia oxidation

Ammonia oxidation is the heart of the Ostwald process and is important in emissions control. Catalytic behaviors are a function of conditions and are observed to vary across the platinum group metals (PGMs) Pt, Pd, Rh. Here, we combine density functional theory computations and microkinetic modeling to rationalize these dependencies. We compute reactions over model (2 1 1) and (1 1 1) surfaces of PGMs. Binding energies are similar on Pd and Pt and generally greater on Rh, while activation energies vary across all metals. Rates on (2 1 1) surfaces are greater than (1 1 1) surfaces. Additionally, the stepped Pt is most active and stepped Rh most selective to N 2 at ammonia slip conditions, while at Ostwald process conditions, stepped Pd is most active and stepped Pt most selective to NO. Degree of rate and selectivity control analysis provides insights into the reactions limiting performance of PGMs. Both activation barriers and surface coverages influence rates and selectivities.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2018

First principles, microkinetic, and experimental analysis of Lewis acid site speciation during ethanol dehydration on Sn-Beta zeolites

Density functional theory calculations are combined with kinetic measurements of ethanol dehydration to diethyl ether to identify the relative catalytic contributions of structurally distinct speciations of Sn sites in zeolite Beta frameworks. The structural complexities of the Beta framework require nonstandard techniques for entropy and energy calculations, including consideration of anharmonic effects in vibrational modes, employment of quasi-harmonic densities of states methods to evaluate entropies, and use of hybrid density functionals to evaluate binding energies. Calculated energies and entropies are used to construct a microkinetic model that is iteratively refined to identify all kinetically and thermodynamically sensitive reaction steps and intermediates which are subsequently treated with the higher-level methods. The rate and equilibrium constants obtained from this tiered approach agree well with measured reaction orders in ethanol and water. Site balances provide evidence for the interconversion of Sn sites between “closed” configurations that are tetra-coordinated to the framework, open configurations formed by hydrolysis that are tri-coordinated to the framework and contain a hydroxyl ligand and proximal silanol group (“hydroxy-open”), and open tri-coordinate configurations formed by reaction with ethanol, yielding an ethoxy ligand and proximal silanol group (“ethoxy-open”). Closed, hydroxy-open, and ethoxy-open Sn sites adsorb ethanol via distinct modes and react via distinct pathways, with the consequence that prevalent dehydration pathways depend on the speciation of Sn sites under reaction conditions. The kinetic modeling indicates that, under the conditions studied (404 K, 0.5–35 kPa ethanol, 0.1–50 kPa water), bimolecular dehydration on the closed Sn site is the sole kinetically-relevant step, and ethanol, ethanol-ethanol dimers, and ethanol-water dimers are the most abundant surface intermediates. In conclusion, these results highlight the importance of considering the distribution of heteroatom coordination modes to zeolite frameworks under reaction conditions, as well as pathways for their interconversion in the presence of reacting molecules, to obtain more robust mechanistic and kinetic interpretations of catalytic pathways in Lewis acid zeolites.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2018

Mechanistic Insights into the Synthesis of Higher Alcohols from Syngas on CuCo Alloys

Synthesis gas (CO + H 2 ) conversion is an important process in the transformation of coal, natural gas, or biomass into higher-value products. The explicit conversion into C 2+ oxygenates on transition-metal-based catalysts suffers from a low selectivity, being a consequence of an imperative integration of C–O bond splitting and C–C coupling reactions. Recently, it has been demonstrated that a bimetallic CuCo catalyst has high higher alcohol selectivity under mild reaction conditions, but the details of the reaction mechanism on the surface are still elusive. In this work, we studied the formation of methane, methanol, and ethanol from syngas on a close-packed (111) and a stepped (211) CuCo surface combining density functional theory (DFT) and microkinetic modeling. We found the CuCo alloy to be a promising candidate catalyst, displaying the required coverage of CO and CH x on the surface to facilitate C–C coupling. In addition, we found the selectivity to be very structure sensitive: the CuCo (211) surface is selective toward ethanol under certain reaction conditions, while the (111) surface is selective toward methanol. As a result, we identified the much lower C–O dissociation barrier and the higher rate of CH x –CO coupling as the reason for the high activity and selectivity toward ethanol on the (211) surface.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Rhodium Single–Atom Catalysts on Titania for Reverse Water Gas Shift Reaction Explored by First Principles Mechanistic Analysis and Compared to Nanoclusters

The thermocatalytic reduction of CO2 by H2 often proceeds via two competing reaction mechanisms—the reverse water gas shift reaction (rWGSR, CO2 + H2 ⇌ CO + H2O) and methanation (CO2 + 4H2 ⇌ CH4 + 2H2O). Atomically dispersed Rh1 catalysts on TiO2 show high selectivity toward the rWGSR compared with larger Rh nanoclusters, but the origin of this size-dependent selectivity has not been fully explained. Here we report density functional theory (DFT) calculations and microkinetic simulations that clarify the Rh1 active sites and rWGSR pathway on anatase TiO2(101), as well as the high rWGSR selectivity of Rh1 compared with supported Rhx (x = 2–8 atoms) nanoclusters. DFT-computed formation energies, vibrational frequency analysis, and microkinetic modeling suggest three plausible active sites: Rh1 on titania (Rh1/TiO2(101)), Rh1 with a nearby hydroxyl group (Rh1OH/TiO2(101)), and Rh1 near an oxygen vacancy at a three-fold coordinated site (Rh1 near O3cvac). Predicted turnover frequencies and apparent activation barriers for Rh1 indicate a faster reaction involving CO2 dissociation assisted by a support oxygen vacancy via Rh1 near O3cvac, as well as slower reactions involving Rh1OH/TiO2(101) or Rh1/TiO2(101) through a COOH intermediate. Furthermore, these Rh1 sites are selective toward CO rather than CH4 because of the weak adsorption of CO, large barrier for C-O bond dissociation, and the lack of nearby metal sites for H2 dissociation, in contrast to Rhx nanoclusters, including Rh2 dimers.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2021

Identification of Active Metal Carbide and Nitride Catalytic Facets for Hydrodeoxygenation Reactions

The catalytic hydrodeoxygenation (HDO) reaction is of considerable interest for biomass conversion to valuable chemicals and fuels, where one of the critical bottlenecks is the lack of cost-effective and efficient catalysts. To discover cost-efficient catalysts for the HDO reaction, we employed a density functional theory-based hierarchical catalyst design strategy based on catalytic descriptors, reaction energy profiles, and microkinetic modeling (MKM). We focused on the carbide and nitride catalyst space, for which we calculated 121 catalyst surfaces of Mo 2 C, MoC, Mo 2 N, W 2 C, NbC, VC, VN, and NbN catalysts. Based on the computed surface energies, reaction energies of oxygen removal, carbon binding strength, and the surface area of nanoparticles, the likely active facets are the Mo 2 C(111), MoC(011), VN(100), Mo 2 N(001), Mo 2 N(011), and Mo 2 N(100) surfaces. Further, detailed energy profiles were obtained, and MKM was performed for a model reaction (glycolaldehyde + 2H 2 . ethylene + 2H 2 O) on the Mo 2 C(111), VN(100), and MoC(100) surfaces. Based on the computed volcano map obtained from MKM, the predicted active facets for this HDO reaction are the Mo 2 C(111), MoC(011), VN(011), Mo 2 N(001), Mo 2 N(011), and Mo 2 N(100) surfaces. Additionally, none of the carbide and nitride catalyst surfaces are located in the optimal catalytic activity part. Therefore, it is essential to modify the catalyst via adding dopants or alloying to improve the catalytic activity. Catalytic modifications that can destabilize the surface adsorption of O*/H 2 O* and decrease the energy barriers of O-H bond formation are recommended to facilitate the HDO on the carbide and nitride catalysts. These a priori investigations provide guidelines for future low-cost HDO catalyst development.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2020

Investigation of the reaction mechanism of the hydrodeoxygenation of propionic acid over a Rh(1 1 1) surface: A first principles study

Microkinetic models based on first principles calculations have been used to study the vapor and liquid phase hydrodeoxygenation of propionic acid on a Rh(111) surface. Calculations suggest that both decarboxylation and decarbonylation do not occur at an appreciable rate under all reaction environments. Propanol and propionaldehyde are the main products on this surface and they are produced at similar rates in both vapor and liquid phase environments. While a condensed phase can shift the reaction rate, the dominant pathways and selectivity towards the various products are hardly affected. At 473 K, the turnover frequency is increased by about a factor of 1.5 in liquid water relative to the gas phase. In liquid 1,4-dioxane, the turnover frequency is also slightly increased relative to the gas phase. Here, given the uncertainty in Rh cavity radius in the liquid phase calculations, computations with different cavity radii have been performed. With larger Rh cavity radius, the promotional effect of the solvents on the turnover frequency becomes more significant, while practically no changes are observed for a smaller cavity radius.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗