Engineering Papers⌕ Search

SEARCH · Engineering Papers

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 73 records · Page 4
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Interfacial Water Manipulation with Ionic Liquids for the Oxygen Reduction Reaction

The role that interfacial water plays in both promoting and inhibiting the electrochemical oxygen reduction reaction (ORR) remains somewhat unclear and controversial. Here, we use electroanalytical chemistry, spectroscopy, and microkinetic modeling to probe the impact of hydrophobic ionic liquid (IL) thin films on interfacial water structure and its role in promoting or inhibiting the ORR. Through the use of in situ ATR-SEIRAS, we find that the IL thin films limit the content of water at the interface and prevent the formation of hydrogen bond stabilized water organization. The impact of this exclusion of water on the promotion of the ORR is through the reduced coverage of OH ad spectator species. The decreased solvation of “active” OH ad species weakens its interaction with the catalyst surface, lowering the barrier to the last step in the ORR mechanism. This “destabilized” OH ad impact on ORR kinetics is confirmed through a microkinetic model. Furthermore, the results presented here highlight the mechanistic pathway through which hydrophobic ILs enhance ORR kinetics and point to pathways to both further improve this performance enhancement and the potential for integration of hydrophobic ILs into other technologically relevant elementary electrochemical reactions.

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

Non-Idealities in Lab-Scale Kinetic Testing: A Theoretical Study of a Modular Temkin Reactor

The Temkin reactor can be applied for industrial relevant catalyst testing with unmodified catalyst particles. It was assumed in the literature that this reactor behaves as a cascade of continuously stirred tank reactors (CSTR). However, this assumption was based only on outlet gas composition or inert residence time distribution measurements. The present work theoretically investigates the catalytic CO2 methanation as a test case on different catalyst geometries, a sphere, and a ring, inside a single Temkin reaction chamber under isothermal conditions. Axial gas-phase species profiles from detailed computational fluid dynamics (CFD) are compared with a CSTR and 1D plug-flow reactor (PFR) model using a sophisticated microkinetic model. In addition, a 1D chemical reactor network (CRN) model was developed, and model parameters were adjusted based on the CFD simulations. Whereas the ideal reactor models overpredict the axial product concentrations, the CRN model results agree well with the CFD simulations, especially under low to medium flow rates. This study shows that complex flow patterns greatly influence species fields inside the Temkin reactor. Although residence time measurements suggest CSTR-like behavior, the reactive flow cannot be described by either a CSTR or PFR model but with the developed CRN model.

Wehinger, Gregor D. (ORCID:0000000217743391)↗
DOE OSTI2023

A Comparative Study of Electrical Double Layer Effects for CO Reduction Reaction Kinetics

Solvation models describe how the interactions between the solutes and solvents affect the reactivity and selectivity in electrochemical processes. In this study, we developed a framework for evaluating the effects of applied potential and electrical double layer on CO reduction (COR), comparing fully explicit, implicit, and hybrid solvation models at the standard hydrogen electrode (SHE) scale. We analyzed all crucial intermediates leading to the production of C 1 and C 2 products and found good agreement across these models. Some notable differences were observed in the implicit description of *C and *CO adsorption at higher overpotentials and overall trends in the adsorption energies within the hybrid model. Using this unified SHE framework, we built comparable microkinetic models for COR kinetics and rates. Despite small differences in thermodynamics, solvation-model-based microkinetic simulations showed good agreement for onset potentials against the benchmark experiment. Only qualitative difference was observed for C 1+ versus hydrogen evolution at high overpotentials for the implicit model. Finally, we constructed a generalized C 2 selectivity map in descriptor space (U SHE , ΔG electrolyte CO ), which highlights the limitations of a copper-based COR catalyst and guides the search for optimal descriptor parameters to maximize C 2 selectivity. In conclusion, these findings demonstrate the importance of considering electrical double layer effects in reduction reactions and offer a useful framework for comparing solvation models and predicting optimal electrochemical conditions for specific applications.

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

Theoretical Investigation of Solvent Effects on the Hydrodeoxygenation of Propionic Acid over a Ni(111) Catalyst Model

The effect of two solvents, liquid water and 1,4-dioxane, has been studied from first principles on the hydrodeoxygenation of propionic acid over a Ni (111) catalyst surface model. A mean-field microkinetic model was developed to investigate these effects at a temperature of 473 K. Under all reaction conditions, a decarbonylation mechanism is favored significantly over a decarboxylation pathway. Although no significant solvent effects were observed on the decarbonylation rate, a substantial solvent stabilization of two key surface intermediates in the decarboxylation mechanism, CH 3 CCOO and CH 3 CHCOO, lead to a notable increase of the decarboxylation rate by two orders of magnitude in liquid water and by one order of magnitude in liquid 1,4-dioxane. Furthermore, a significant solvent stabilization of the transition state of C-H bond cleavage of the α-carbon of CH 3 CHCO, relative to the stabilization of the C-C bond cleavage of the α-carbon of CH 3 CHCO, leads to a change in dominant pathway in the liquid phase environments. Finally, a sensitivity analysis shows that the C-OH bond cleavage of propionic acid and C-C bond cleavage of the α-carbon of CH 3 CHCO are the most rate controlling states in the gas phase. In contrast, in solvents the dehydrogenation of CH 3 CHCO becomes the most influential step. This shift in rate controlling state is attributed to the solvent effect on the dehydrogenation of CH 3 CHCO, which is facilitated in aqueous phase. Altogether, it is likely that the investigated (111) facet of Ni is not active for the hydrodeoxygenation of propionic acid in neither the gas nor liquid phase and other Ni facets or phases must be responsible for the experimentally observed kinetics.

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

Multi‐Scale Model‐Informed Deep Learning for Plasma‐Nanoparticle Interaction

The Overarching Goal of this proposed research is to understand and quantitively determine the interactions between non-thermal plasma (hot electrons, reactive radicals, vibrationally excited species) and surface reactions on influencing the activity and selectivity of the desired reactions via developing multi-scale model informed deep learning algorithm. Investigating non-thermal plasma-surface interaction is feasible due to the low bulk temperature in the discharge region. To investigate the role of plasma-nanoparticle interaction on enhancing the reaction kinetics, we will focus on ammonia cracking to generate clean hydrogen over earth-abundant, non-critical metallic nanoparticles, which is of great significance for decarbonization. We hypothesize that (1) reactive radicals interacting with surface reaction species via Eley–Rideal mechanism will significantly lower the energetics of the potential rate-limiting step of nitrogen formation; (2) the surface will be charged heterogeneously under non-thermal plasma conditions and the charged site will lower the energetics of ammonia cracking through Langmuir– Hinshelwood mechanism; (3) vibrationally excited ammonia will further promote the initial N-H bond cleavage. To access the hypothesis, we will (1) reveal the surface charge effects on tunning the reaction energetics via interpretable, physics-informed deep learning accelerated density functional theory (DFT) calculations; (2) determine the reactive radicals interacting with surface reaction species on tuning the reaction energetics via DFT; (3) reveal the surface charge effects on tunning the reaction energetics via DFT and deep learning models, (4) quantify how vibrationally excited species, reactive radicals, and surface charging effects on enhancing the catalysis via developing DFT-based microkinetic modeling (MKM) and active learning. Deep and active learning of plasma-nanoparticle interactions effects on enhancing ammonia cracking to generate hydrogen represents a new paradigm for designing high performance plasma materials. The fundamental science of how plasma-nanoparticle interactions will change the plasma kinetics and will improve the energy efficiency for decarbonization and sustainability. The interpretable and physics-informed machine learning model will accelerate low temperature plasma chemistry and material discovery with physics rules and model interpretation.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2026

Group-Additivity–Embedded Multiscale Modeling for Electric Field-Enhanced Nanocatalysis

Elucidating structure-performance relationships remains a central challenge in field-enhanced catalysis, where nanoparticles exhibit nonuniform surface sites with site-dependent responses to electric fields. Low-coordination sites (edges, corners, and tips) are particularly electric field-sensitive (EF), leading to nonuniform charge distribution, adsorption energies, and catalytic activity. Here, using ammonia decomposition on a ruthenium cluster as a model system, we develop a transferable multiscale framework integrating density functional theory, group additivity (GA), Brønsted-Evans-Polanyi scaling, and microkinetic modeling to predict EF-dependent activity across nonuniform cluster sites. Across sites and fields, the nitrogen adsorption energy (E N ) emerges as the governing descriptor, yielding robust volcano relationships whose optimum shifts systematically with field: negative fields strengthen N binding via electron accumulation, while positive fields weaken N binding via charge depletion, moving the optimal E N toward weaker binding. Microkinetic analysis shows that N≡N bond formation remains the key kinetic bottleneck over most conditions; positive fields lower the effective barrier and, critically, increase the fraction of near-optimal active sites, leading to a net enhancement in overall activity relative to zero-field and negative-field cases. By capturing EF- and site-dependent energetics with high accuracy and low computational cost, this GA-embedded multi-scale simulation workflow provides a physically interpretable route to predict and design field-enhanced nanocatalysis.

ammonia decomposition↗
DOE OSTI2026

Interpretable Deep Learning for Advancing Field-Enhanced Catalysis

This DOE Early Career project developed a physics-informed, interpretable AI-and-modeling framework to understand and exploit electric-field effects in heterogeneous catalysis, with ammonia cracking and synthesis as a representative pathway. The team built and validated methods to map local electric fields on metal surfaces and nanoparticles, showing that low-coordination features (tips/edges/corners) can concentrate fields by several-fold relative to flat facets. Using DFT-generated datasets, the project created physics-guided machine learning models that rapidly predict local electric fields and field-dependent adsorption energetics with near-DFT accuracy while reducing computational cost by orders of magnitude. These predictions were integrated with microkinetic modeling to quantify how field-dipole interactions reshape reaction energetics and mechanisms, enabling large increases in predicted catalytic rates and substantial reductions in operating temperature under favorable field conditions. To accelerate discovery of earth-abundant catalysts, the project combined interpretable ML screening (with electronic-structure descriptors identified as key drivers) with a generative inverse-design workflow based on diffusion models and physics constraints. The resulting closed-loop approach, linking simulation, mechanistic modeling, and AI, provides reusable tools and datasets for designing catalysts and operating conditions in field-enhanced catalysis, with broad relevance to electrostatic catalysis, plasma catalysis, electrocatalysis, and other energy-related chemical transformations.

30 DIRECT ENERGY CONVERSION↗
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 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 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 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 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↗