Search NASA⌕ Search

SEARCH · Search NASA

Results for “microkinetics”

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 OSTI2020

Oxidative dehydrogenation of propane on the oxygen adsorbed edges of boron nitride nanoribbons

Metal-free boron nitride has recently been reported to exhibit high catalytic activity and selectivity towards oxidative dehydrogenation of propane (ODHP). Despite several experimental studies, the exact nature and function of the active sites in this emerging catalyst are unknown. In the present work, density functional theory calculations are combined with microkinetic modeling to systematically explore the ODHP on oxygen passivated boron nitride nanoribbons (BNNRs). Here, the relative stabilities of different edge structures of BNNR were examined using ab initio thermodynamic analysis, and the most stable oxygen passivated edge structure was used for a mechanistic study. Microkinetic analysis revealed that a N 2 O or NO x -type active site is active and selective for the conversion of propane to propene. In addition to a heterogeneous catalytic cycle, the proposed N 2 O/NO x -type active sites are expected to be able to generate gas phase C 3 H 7 ˙, which in turn can trigger gas phase reactions of ODHP as experimentally speculated.

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

Passive NO x adsorption on Pd/H-ZSM-5: Experiments and modeling

A combined experimental and modeling study of the passive NO x adsorber (PNA) is presented that advances the understanding and prediction of the effects of various operating parameters and material properties of H-ZSM-5 and Pd/H-ZSM-5. Experiments reveal the NO x uptake process to be kinetically limited. Two microkinetic schemes are developed to explain the NO x uptake and release and NO 2 generation data. Both schemes are incorporated into a 1+1 D monolith model to predict and validate NO x uptake and temperature programmed desorption (TPD) data. Scheme I involves Z – [PdOH] + , Z – Pd 2+ Z – and Z – Pd + with Z – Pd + as the strongest NO binding site. Scheme II involves reduction of PdO 2 to PdO and nitrate formation. Model tuning utilizes a combination of infrared spectroscopy measurements and density functional theory (DFT) estimates of energy barriers. Lastly, the model is validated at different conditions including feed parameters and Pd loading.

42 ENGINEERING↗
DOE OSTI2021

NO x adsorption with CO and C 2 H 4 on Pd/SSZ-13: Experiments and modeling

A transient monolith model containing microkinetic schemes for NO uptake and release over Pd/SSZ-13 without and with reductants CO or C2H4 is presented. The scheme involves three cationic Pd species (Z – [PdOH] + , Z – Pd 2+ Z – , Z – Pd + ) as the active sites during uptake of NO, CO and C 2 H 4 , and their desorption and conversion at higher temperature. Kinetic parameters are estimated through a combination of density functional theory (DFT) estimates and a fit of uptake, desorption and conversion data. The tuned model is validated at different uptake temperatures, ramp rates, and flowrates. A “degree of uptake control” parameter is defined that helps to identify the step(s) that are uptake controlling. The model helps to interpret the data features and is used to identify operating conditions to meet application-relevant performance metrics, including NO trapping efficiency and NO release temperature. Furthermore, the model demonstrates higher NO uptake on Pd/SSZ-13 compared to Pd/ZSM-5.

42 ENGINEERING↗
DOE OSTI2021

Efficient Screening of Bi–Metallic Electrocatalysts for Glycerol Valorization

Glycerol is a byproduct of biodiesel production and, as such, it is of limited economic value. By means of electrooxidation, glycerol can be used as a feedstock for scalable hydrogen production, in addition to conversion to value-added products. The development of novel and efficient catalytic electrode materials for the anodic side of the reaction is a key towards a hydrogen-based energy economy. In the present study, a computational screening protocol combining DFT, scaling relations, and microkinetic modeling allows for a rational selection of novel catalysts that can deliver efficient glycerol electrooxidation, low cost of production, and environmental sustainability. Activity and chemical selectivity towards hydrogen production on pure metal catalysts is discussed in terms of volcano-shaped plots. We find that the selectivity in the glycerol oxidation reaction is influenced by a different energy landscape when in the presence of water and best classified by a comparison of O—H and C—H bond-breaking barriers. In addition, we screened 3570 bi-metallic catalysts in the AB (L1 0 ) and A 3 B (L1 2 ) ordered structures for activity, stability, price, and toxicity. By filtering based on the criteria for toxicity, resistance to oxidation, miscibility, and price, we have identified 5 L1 0 structured catalysts (AgPd, AuPd, PtSb, CuPt, and AgPt) and 20 L1 2 catalysts (Ga 3 Ta, In 3 Ta, Ir 3 W, Ir 3 Mo, Cu 3 Pt, Ir 3 Ta, Ir 3 Re, Pd 3 Bi, Pd 3 Cu, Pd 3 W, Pd 3 Co, Pd 3 Sn, Pd 3 Mo, Pd 3 Ag, Pd 3 Ga, Pd 3 Ta, Au 3 Ru, Pd 3 In, Au 3 Ir, and Pd 3 Au) that are all predicted to show high activity. We also identify an additional 37 L1 0 and 92 L1 2 structured electrocatalysts with an anticipated medium-high activity.

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

Mechanistic understanding of pH effects on the oxygen evolution reaction

The oxygen-evolution reaction (OER) is pivotal in many energy-conversion technologies as it is an important counter reaction to others that convert stable chemicals to higher-value products using electrochemistry. The local microenvironment and pH for the anode OER can vary from acidic to neutral to alkaline depending on the system being explored, making definitive mechanistic insights difficult. In this paper, we couple experiments, first-principles calculations based on density functional theory, microkinetics, and transport modeling to explore the entire pH range of the OER. At low current densities, neutral pH values unexpectedly perform better than the acidic and alkaline conditions, and this trend is reversed at higher current densities (> 20 mA cm -2 ). Using multiscale modeling, this switch is rationalized by a change from a dual-reaction mechanism to a single rate-determining step. The model also shows how the alkaline reaction rates dominate in the middle to high pH range. Furthermore, we explore that the local pH for near-neutral conditions is much different (e.g., 2.4 at the reaction surface vs. 9 in the bulk) than the pH extremes, demonstrating the criticality that transport phenomena plays in kinetic activity.

36 MATERIALS SCIENCE↗
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↗
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 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

Importance of Site Diversity and Connectivity in Electrochemical CO Reduction on Cu

Electrochemical CO 2 reduction on Cu is a promising approach to produce value-added chemicals using renewable feedstocks, yet various Cu preparations have led to differences in activity and selectivity toward single and multicarbon products. Here, we find, surprisingly, that the effective catalytic activity toward ethylene improves when there is a larger fraction of less active sites acting as reservoirs of *CO on the surface of Cu nanoparticle electrocatalysts. In an adaptation of chemical transient kinetics to electrocatalysis, we measure the dynamic response of a gas diffusion electrode (GDE) cell when the feed gas is abruptly switched between Ar (inert) and CO. When switching from Ar to CO, CO reduction (COR) begins promptly, but when switching from CO to Ar, COR can be maintained for several seconds (delay time) despite the absence of the CO reactant in the gas phase. A three-site microkinetic model captures the observed dynamic behavior and shows that Cu catalysts exhibiting delay times have a less active *CO reservoir that exhibits fast diffusion to active sites. The observed delay times and the estimated *CO reservoir sizes are affected by catalyst preparation, applied potential, and microenvironment (electrolyte cation identity, electrolyte pH, and CO partial pressure). Notably, we estimate that the *CO reservoir surface coverage can be as high as 88 ± 7% on oxide-derived Cu (OD-Cu) at high overpotentials (–1.52 V vs SHE) and this increases in reservoir coverage coincide with increased turnover frequencies to ethylene. We also estimate that *CO can travel substantial distances (up to 10s of nm) prior to desorption or reaction. It appears that active C–C coupling sites by themselves do not control selectivity to C 2+ products in electrochemical COR; the supply of CO to those sites is also a crucial factor. More generally, the overall activity of Cu electrocatalysts cannot be approximated from linear combinations of individual site activities. Future designs must consider the diversity of the catalyst network and account for intersite transportation pathways.

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

Tuning Active Site Flexibility by Defect Engineering of Graphene Ribbon Edge‐hosted Fe−N 3 Sites

Abstract Nitrogen‐doped, carbon‐supported transition metal catalysts are excellent for several reactions. Structural engineering of M−N x sites to boost catalytic activity is rarely studied. Here, we demonstrate that the structural flexibility of Fe−N 3 site is vital for tuning the electronic structure of Fe atoms and regulating the catalytic transfer hydrogenation (CTH) activity. By introducing carbon defects, we construct Fe−N 3 sites with varying Fe−N bond lengths distinguishable by X‐ray absorption spectroscopy. We investigate the CTH activity by density‐functional theory and microkinetic calculations and reveal that the vertical displacement of the Fe atom out of the plane of the support, induced by the Fe−N 3 distortion, raises the Fe orbital and strengthens binding. We propose that the activity is controlled by the relaxation of the reconstructed site, which is further affected by Fe−N bond length, an excellent activity descriptor. We elucidate the origin of the CTH activity and principles for high‐performing Fe−N−C catalysts by defect engineering.

Yang, Piaoping↗
DOE OSTI2023

Tuning Active Site Flexibility by Defect Engineering of Graphene Ribbon Edge–hosted Fe–N3 Sites

Nitrogen-doped, carbon-supported transition metal catalysts are excellent for several reactions. Structural engineering of metal-Nx sites to boost catalytic activity is rarely studied. Here, we demonstrate that the structural flexibility of Fe-N3 site is vital for tuning the electronic structure of Fe atoms and regulating the catalytic transfer hydrogenation (CTH) activity. By introducing carbon defects, we construct Fe-N3 sites with varying Fe-N bond lengths distinguishable by X-ray absorption spectroscopy. We investigate the CTH activity by density functional theory and microkinetic calculations and reveal that the vertical displacement of the Fe atom out of the plane of the support, induced by the Fe-N3 distortion, raises the Fe [[EQUATION]]orbital and strengthens binding. We propose that the activity is controlled by the relaxation of the reconstructed site, which is further affected by Fe-N bond length, an excellent activity descriptor. Furthermore, we elucidate the origin of the CTH activity and principles for high-performing Fe-N-C catalysts by defect engineering.

09 BIOMASS FUELS↗
DOE OSTI2024

Theoretical Investigation of the Adsorbate and Potential–Induced Stability of Cu Facets During Electrochemical CO 2 and CO Reduction

The activity and product selectivity of electrocatalysts for reactions like the carbon dioxide reduction reaction (CO 2 RR) are intimately dependent on the catalyst's structure and composition. While engineering catalytic surfaces can improve performance, discovering the key sets of rational design principles remains challenging due to limitations in modeling catalyst stability under operating conditions. Herein, we perform first-principles density functional calculations adopting implicit solvation methods with potential control to study the influence of adsorbates and applied potential on the stability of different facets of model Cu electrocatalysts. Using coverage dependencies extracted from microkinetic models, we describe an approach for calculating potential and adsorbate-dependent contributions to surface energies under reaction conditions, where Wulff constructions are used to understand the morphological evolution of Cu electrocatalysts under CO 2 RR conditions. Here we identify that CO*, a key reaction intermediate, exhibits higher kinetically and thermodynamically accessible coverages on (100) relative to (111) facets, which can translate into an increased relative stabilization of the (100) facet during CO 2 RR. Our results support the known tendency for increased (111) faceting of Cu nanoparticles under more reducing conditions and that the relative increase in (100) faceting observed under CO 2 RR conditions is likely attributed to differences in CO* coverage between these facets.

30 DIRECT ENERGY CONVERSION↗