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 145 records · Page 8

Catlas: an automated framework for catalyst discovery demonstrated for direct syngas conversion

Catalyst discovery is paramount to support access to energy and key chemical feedstocks in a post fossil fuel era. Exhaustive computational searches of large material design spaces using ab initio methods like density functional theory (DFT) are infeasible. We seek to explore large design spaces at relatively low computational cost by leveraging large, generalized, graph-based machine learning (ML) models, which are pretrained and therefore require no upfront data collection or training. We present Catlas, a framework that distributes and automates the generation of adsorbate-surface configurations and ML inference of DFT energies to achieve this goal. Catlas is open source, making ML assisted catalyst screenings easy and available to all. To demonstrate its efficacy, we use Catlas to explore catalyst candidates for the direct conversion of syngas to multi-carbon oxygenates. For this case study, we explore 947 stable/metastable binary, transition metal intermetallics as possible catalyst candidates. On this subset of materials, we are able to predict the adsorption energy of key descriptors, *CO and *OH, with good accuracy (0.16, 0.14 eV MAE, respectively). Using the projected selectivity towards C 2+ oxygenates from an existing microkinetic model, we identified 144 candidate materials. For 10 promising candidates, DFT calculations reveal a good correlation with our assessment using ML. Finally, among the top elemental combinations were Pt–Ti, Pd–V, Ni–Nb, and Ti–Zn, all of which appear unexplored experimentally.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Solvation of furfural at metal–water interfaces: Implications for aqueous phase hydrogenation reactions

Metal–water interfaces are central to understanding aqueous-phase heterogeneous catalytic processes. However, the explicit modeling of the interface is still challenging as it necessitates extensive sampling of the interfaces’ degrees of freedom. Herein, we use ab initio molecular dynamics (AIMD) simulations to study the adsorption of furfural, a platform biomass chemical on several catalytically relevant metal–water interfaces (Pt, Rh, Pd, Cu, and Au) at low coverages. We find that furfural adsorption is destabilized on all the metal–water interfaces compared to the metal–gas interfaces considered in this work. This destabilization is a result of the energetic penalty associated with the displacement of water molecules near the surface upon adsorption of furfural, further evidenced by a linear correlation between solvation energy and the change in surface water coverage. To predict solvation energies without the need for computationally expensive AIMD simulations, we demonstrate OH binding energy as a good descriptor to estimate the solvation energies of furfural. Using microkinetic modeling, we further explain the origin of the activity for furfural hydrogenation on intrinsically strong-binding metals under aqueous conditions, i.e., the endothermic solvation energies for furfural adsorption prevent surface poisoning. In conclusion, our work sheds light on the development of active aqueous-phase catalytic systems via rationally tuning the solvation energies of reaction intermediates.

aqueous-phase hydrogenation↗

Catalytic CO Oxidation on MgAl 2 O 4 -Supported Iridium Single Atoms: Ligand Configuration and Site Geometry

Probing and understanding the local chemical environment of an active site is essential for designing high-performance single-atom catalysts (SACs). Density functional theory (DFT) calculations were performed to investigate the ligand configuration and site geometry of MgAl 2 O 4 -supported iridium single atoms (Ir-1) toward catalytic carbon monoxide (CO) oxidation. We employed MgAl 2 O 4 (111) and MgAl 2 O 4 (211) as the model substrates with adsorbed Ir single atoms of different site geometries. DFT calculations revealed that the Mg-site on MgAl 2 O 4 (111) and the step site on MgAl 2 O 4 (211) are the most stable adsorption sites for Ir single atoms. Irrespective of site choices, CO oxidation on supported Ir single atoms follows the Eley-Rideal (E-R) mechanism, in which the surface oxygen vacancies close to the Ir single atoms activate molecular O 2 with a negligible barrier and the rate-limiting step is the gas-phase CO directly attacking the O-Ir species that is modulated by a CO ligand. Furthermore, first-principles X-ray absorption near-edge spectra of reaction intermediates along with in situ/operando X-ray absorption spectroscopy (XAS) suggest that Ir single atoms adsorb primarily on the step sites of MgAl 2 O 4 . However, microkinetic modeling predicts that a higher activity can be attained on the equally stable Mg-site, maximizing the population of which in catalyst synthesis might prove fruitful in future studies. Electronic structure analysis indicates that the CO ligand increases the reactivity of adsorbed oxygen atoms bound to Ir single atoms by increasing the antibonding nature of the O-Ir bond.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Surface intermediates identified using IR spectroscopy during the catalytic oxidation of C 2 H 6 on IrO 2 (1 1 0) at near-ambient pressures

Characterizing surface intermediates during catalysis is essential for developing mechanistic models of catalytic reactions and for identifying rate-controlling steps. Here, in this work, we used polarization-dependent reflection absorption infrared spectroscopy (PD-RAIRS) to investigate surface intermediates that form on IrO 2 (1 1 0) during the catalytic oxidation of C 2 H 6 at pressures near 0.5 Torr. Our results show that adsorbed CO and HCO 2 are the only intermediates detectable with PD-RAIRS across a wide range of temperatures and reactant compositions during complete C 2 H 6 oxidation, giving rise to dominant peaks at about 2093 and 1310 cm −1 , respectively. These assignments were corroborated by RAIRS measurements following CO and HCOOH adsorption in ultrahigh vacuum, RAIRS with C 2 H 6 – 18 O 2 mixtures during catalysis and vibrational frequencies and IR intensities computed using density functional theory. We find that increasing the C 2 H 6 mole fraction in the reaction mixture raises the CO surface coverage relative to HCO 2 , correlating with increasing catalytic oxidation rates. The formation of CO and HCO 2 surface species at temperatures well below catalytic ignition indicates that the overall rate of C 2 H 6 oxidation on IrO 2 (1 1 0) is governed by steps late in the reaction sequence, such as H 2 O formation, CO oxidation and HCO 2 dehydrogenation. These findings provide key insight into the mechanism of alkane oxidation on IrO 2 (1 1 0) and valuable input for first-principles microkinetic modeling.

Pope, Connor [Univ. of Florida, Gainesville, FL (U↗

Theoretical and Experimental Studies of CoGa Catalysts for the Hydrogenation of CO 2 to Methanol

We report methanol is an important chemical compound which is used both as a fuel and as a platform molecule in chemical production. Synthesizing methanol, as well as dimethyl ether, directly from carbon dioxide and hydrogen produced using renewable electricity would be a major step forward in enabling an environmentally sustainable economy. We utilize density functional theory combined with microkinetic modeling to understand the methanol synthesis reaction mechanism on a model CoGa catalyst. A series of catalysts with varying Ga content are synthesized and experimentally tested for catalytic performance. The performance of these catalysts is sensitive to the Co:Ga ratio, whereby increased Ga content results in increased methanol and dimethyl ether selectivity and increased Co content results in increased selectivity towards methane. We find that the most active catalysts have up to 95% CO-free selectivity towards methanol and dimethyl ether during CO 2 hydrogenation and are comparable in performance to a commercial CuZn catalyst. Finally, using in situ DRIFTS we experimentally verify the presence of a surface formate intermediate during CO 2 hydrogenation in support of our theoretical calculations.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Microkinetic simulation and fitting of the temperature programmed reaction of methanol on CeO 2 (111): H 2 and H 2 O + V production

The kinetics and mechanism for temperature programmed reaction following adsorption of an adsorbate can be better understood by simulation and fitting with comparison to experiment. A case study is presented here for the chemistry of adsorption of methanol on a CeO 2 (111) surface followed by heating. The gas products observed are CH 3 OH, CH 2 O, H 2 , H 2 O, CO, CO 2 . At low temperatures (< 500 K), there is formation of H 2 and H 2 O, where the H 2 O formation is accompanied by lattice oxygen vacancy (V) formation and is thus important in determining the selectivity towards different products. Microkinetic modeling was performed using a recently published method for fitting to gain mechanistic knowledge of the H 2 and H 2 O + V formation at < 500 K. In the kinetic models used here, most of the H 2 and H 2 O + V V formation can be explained by a mechanism in which a metastable state of hydrogen on the surface (H*) acts as an intermediate. Two possibilities were investigated for the source of the metastable H* intermediate: H* from CH bond breaking of methoxies, or promotion of H+ to H* via electron transfer from ionic methoxies absorbed in oxygen vacancies (CH 3 O - /V). From this study, we consider the latter to be more likely at < 500 K. For the H 2 O formation, it was found to be critical that H 2 O cannot dissociate directly on oxygen vacancies. Finally, catalytic chemistry was observed in simulations, including catalytic formation of oxygen vacancies. Various features of the experimental results were reproduced, including methoxies being the major carbon containing species on the surface at < 500 K.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Mechanistic Insights into the Catalytic Condensation of Methyl Ketones on MgO Surfaces

Ketone coupling via aldol condensation is one of the promising routes to produce cyclic and value-added precursors for renewable hydrocarbon biofuels. A first-principles-based microkinetic modeling is performed to evaluate the surface-mediated reaction mechanisms and the role of water molecules in the observed activities for 2-pentanone and 3-pentanone aldol condensation on dehydroxylated MgO(111) surface and hydroxylated terminated surface[OH-MgO(111)]. Here, we have identified the enhancement of the surface OH group to MgO(111) surface catalytic activity by destabilizing the binding strength of reaction intermediates and reducing the energy barriers of rate-determining steps(proton transfer and dehydration steps). The 2-pentanone has one elementary step less in the complete reaction mechanism of aldol condensation and preferable energy barrier for proton transfer and dehydration steps, revealing 2-pentanone as terminal ketone is more reactive than 3-pentanone as central ketone. The water molecules dominated the OH-MgO( 111) surface after further addition of water, leading to the reduction of turnover frequency of the aldol condensation dimer product as the loss of aldol condensation reaction intermediates in competitive adsorption with water molecules.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Bromine and iodine for selective partial oxidation of propane and methane

The oxidative dehydrogenation of propane to propylene can be enhanced through the addition of a halogen to a feed of oxygen and propane. Here, we examine halogen-assisted oxidative dehydrogenation by a mixture of I 2 and Br 2 . We find that adding a small amount of Br 2 to a feed of I 2 , C 3 H 8 , and O 2 produces propylene with 71% propane conversion and 83% propylene selectivity. Good performance is also obtained if a small amount of I 2 is added to a feed of Br 2 , C 3 H 8 , and O 2 . In both cases, the propylene yield is greater than the sum of the yields obtained with iodine and bromine separately. A variety of studies were performed to help understand the reaction mechanism. Methane halogenation was also investigated by using combinations of iodine and bromine to produce methyl iodide, and a microkinetic model was also constructed for methane halogenation in the gas-phase to help understand the mechanism.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Kinetics-Informed Neural Networks

Chemical kinetics and reaction engineering consists of the phenomenological framework for the disentanglement of reaction mechanisms, optimization of reaction performance and the rational design of chemical processes. Here, we utilize feed-forward artificial neural networks as basis functions to solve ordinary differential equations (ODEs) constrained by differential algebraic equations (DAEs) that describe microkinetic models (MKMs). We present an algebraic framework for the mathematical description and classification of reaction networks, types of elementary reaction, and chemical species. Under this framework, we demonstrate that the simultaneous training of neural nets and kinetic model parameters in a regularized multi-objective optimization setting leads to the solution of the inverse problem through the estimation of kinetic parameters from synthetic experimental data. We analyze a set of scenarios to establish the extent to which kinetic parameters can be retrieved from transient kinetic data, and assess the robustness of the methodology with respect to statistical noise. Furthermore, this approach to inverse kinetic ODEs can assist in the elucidation of reaction mechanisms based on transient data.

36 MATERIALS SCIENCE↗

Mechanism and kinetics of 1-dodecanol etherification over tungstated zirconia

Growing interest in finding renewable alternatives to conventional fossil fuels and petroleum-derived specialty chemicals has motivated the investigation of biomass-derived alcohols to make ethers as diesel additives or lubricants. To optimize the direct etherification of long chain alcohols in the liquid phase, it is necessary to develop an understanding of the kinetics and mechanism of etherification and dehydration reactions. In this study, tungstated zirconia was identified as a selective solid-acid catalyst for the liquid-phase etherification of 1-dodecanol. Investigations of the mechanism and kinetics of this reaction suggest that cooperation between Brønsted- and Lewis-acid sites on tungstated zirconia enhances the selectivity to ether by increasing the surface concentration of adsorbed alcohol, thereby promoting bi-molecular ether formation relative to unimolecular alcohol dehydration. The suggested rate limiting step for etherification is the formation of a C–O bond between two adsorbed alcohol molecules, and the suggested rate-limiting step for dehydration is the cleavage of the C–H bond of the β-carbon atom in an adsorbed alcohol. Measurements of the kinetic isotope effects for etherification and dehydration support the proposed mechanism. Finally, a microkinetic model based on the proposed mechanism for dodecanol etherification and dehydration over tungstated zirconia accurately describes the observed effects of alcohol concentration and product inhibition.

di-dodecylether↗

Catalytic resonance theory: Circumfluence of programmable catalytic loops

Chemical reactions on heterogeneous catalyst surfaces exhibit complex networks of elementary reactions with multiple pathways to fluid phase products, sometimes leading to surface reaction loops consisting of a closed cycle reaction pathway. While conventional catalysts at steady state exhibit zero net flux in either direction around a catalytic loop, the loop turnover frequency of three-species surface loops was evaluated in this work via microkinetic modeling to assess the reaction loop behavior resulting from a catalytic surface oscillating between two or more surface energy states. For dynamic heterogeneous catalysts undergoing applied oscillations of surface energy (i.e., programs), surface reaction loops of three species were shown to exhibit non-zero net flow of molecules around the loop, with the extent of loop turnover frequency varying with the applied frequency and amplitude of a sinusoidal or square wave oscillation. Alternatively, some dynamic surface reactions exhibited oscillations only between two surface species or resulted in surfaces covered by a single species at all times. Here, more complex behavior was observed for dynamic surface programs of three distinct electronic states, with the temporal order of applied states controlling the direction of flow of molecules within a three-molecule surface loop. Catalytic loops have the potential to limit overall catalytic reaction rates and use energy in programmable catalysts, while some applications could purposely impose non-zero loop turnover frequency for improved surface reaction control.

Catalysis↗

Activity Trends for Catalytic CO and NO Co-Oxidation at Low Temperature Diesel Emission Conditions

The diesel oxidation catalyst (DOC) is an essential component of modern vehicle emissions control systems. The pervasive challenge for low temperature oxidation of engine exhaust gas is the mutual inhibition between the various pollutants, causing a marked increase in light-off temperature. Using a combination of density functional theory and descriptor-based microkinetic modeling, we have screened catalysts for low temperature co-oxidation of CO and NO with specific emphasis on minimizing inhibition effects. Compared to standard Pt–Pd alloys, we find that coinage metal alloys, i.e., Cu, Ag, and Au with at least one oxophilic constituent, should possess more robust low temperature activity with minimal inhibition. We attribute this remarkable performance to high surface concentrations of oxygen due to the oxophilic component and less competitive adsorption between CO and NO to the exposed coinage metal sites. We believe that these fundamental insights provide valuable design principles for improved low temperature oxidation catalysts.

02 PETROLEUM↗

Computational Investigation of the Catalytic Hydrodeoxygenation of Propanoic Acid over a Cu(111) Surface

Cu-based alloy catalysts have recently been investigated experimentally for the hydrodeoxygenation (HDO) of biomass-derived organic acids. Here, the HDO of propanoic acid (PAc) has been studied over Cu(111) by mean-field microkinetic modeling based on parameters obtained from first-principles calculations. Models were developed for the gas- and liquid-phase HDO in condensed water and 1,4-dioxane. In agreement with experimental observations, the gas-phase PAc conversion rate is low at 573 K and increases in liquid water by 1 order of magnitude. In all reaction environments, the decarboxylation mechanism is dominant at low hydrogen partial pressures less than 0.1 bar, and the C–COO bond dissociation is the rate-controlling elementary step. This observation contrasts with the rate-controlling step identified over most group VIII metal surfaces, which is the C–OH bond dissociation in the decarbonylation mechanism. At high hydrogen (H 2 ) partial pressures greater than 10 bar, the HDO of PAc produces propionaldehyde that can readsorb and further react through decarbonylation to produce C 2 alkane products, which is conceptually different from the low H 2 partial pressure scenario. At high H 2 partial pressures, the initial hydrogenation at the carbonyl carbon of PAc becomes the rate-controlling elementary step.

Catalysts↗

Insights into the Hydrogen Evolution Reaction on 2D Transition-Metal Dichalcogenides

Understanding hydrogen evolution reaction (HER) behaviors over two-dimensional transition-metal dichalcogenides (2D-TMDs) is critical for the development of nonprecious HER electrocatalysts with better activity. Here, in this work, by combining density functional theory calculations with microkinetic modeling, we thoroughly investigated the HER mechanism on 2D-TMDs. We find an important dependence of simulated cell size on the calculated hydrogen adsorption energy and the activation barrier for MoS 2 . Distinct from previous “H migration” mechanisms proposed for the Heyrovsky reaction, the rate-determining step for MoS 2 , we propose that the Mo site only serves as the stabilized transition state rather than H adsorption. In comparison to transition-metal electrocatalysts, we find that the activation barrier of the Heyrovsky reaction on 2D-TMDs scales with the hydrogen adsorption energy exactly as for transition metals except that all activation energies are displaced upward by ca. 0.4 eV. This higher Heyrovsky activation barrier is responsible for the substantially lower activity of 2D-TMDs. We further show that this higher activation barrier stems from the more positively charged adsorbed hydrogen on the chalcogenides interacting repulsively with the incoming proton. Based on these insights, we discuss potential strategies for the design of nonprecious HER catalysts with activity comparable to Pt.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Point Defects Control Guest Molecule Diffusion in the 1D Pores of Zn(tbip)

Molecular diffusion plays a critical role in metal-organic frameworks (MOFs) within the application of kinetic chemical separations. We carefully study in this work the unexpected role of point defects for short-chain alkanes diffusing in Zn(tbip), an MOF with rigid one-dimensional (1D) channels. Inside a defect-free Zn(tbip), guest molecules are expected to follow single-file diffusion along 1D channels. It has been found previously by Heinke et al. that these parallel 1D channels are connected for molecular diffusion. Our density functional theory (DFT) calculations suggest that linker vacancy defects could arise under experimentally relevant conditions by removing a pair of linkers. Further climbing-image nudged elastic band (cNEB) DFT calculations indicate that hopping of short-chain alkanes between adjacent 1D channels over defect windows can occur at moderate temperatures. In addition to providing insights into connected adjacent 1D channels in Zn(tbip), Heinke et al. also inferred that most 1D pores are blocked from a microkinetic model. Additionally, we explored the influence of hydrolyzed linker created by the formation of linker removal inside 1D pores. Our DFT calculations show that the linkers can effectively block the pores and the linker diffusion in 1D channels is slow. Our results, for the first time, offer a mechanistic explanation of the unexpected molecular diffusion behavior in this MOF with 1D channels.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Unraveling 1-Hexene Hydrogenation over Dilute Pd-in-Au Alloys

Dilute Pd-in-Au alloys are valuable materials for selectively hydrogenating alkynes and isomerizing alkenes. By diluting Pd in a Au host, the selectivity toward semihydrogenated alkene isomers can be significantly enhanced and the unfavorable overhydrogenation to alkanes is suppressed. However, a detailed mechanistic study on the origin of the enhanced alkene selectivity over dilute alloy catalysts is still missing. Here, we combine experiment and theory to unravel the reaction mechanism, identifying rate-limiting and selectivity-controlling steps in 1-hexene hydrogenation over dilute Pd-in-Au catalysts. Using isotope-exchange hydrogenation experiments, we show that 1-hexene and hydrogen over a bimetallic Pd4Au96 in silica catalyst preferentially form 1-hexene isomers, (trans and cis) 2- and 3-hexene and only a small amounts of hexane. The reaction is consistent with a Horiuti–Polanyi mechanism, similar to a monometallic Pd nanoparticle catalyst. Computation of the free-energy profiles for 1-hexene hydrogenation and isomerization over a single Pd atom in a Au surface using first principles calculations indicated that the isomerization of 1-hexene to 2-hexene is energetically favorable due to the relatively large barrier for H2 dissociation preventing hydrogenation to n-hexane. Microkinetic modeling established that H2 dissociation on the single-atom Pd sites and H spillover from these sites onto the Au host are rate-limiting and key in steering the selectivity of dilute Pd-in-Au alloys toward the hexene isomers. Furthermore, the mechanistic insights from this study contribute to the rational design of optimized dilute alloy catalysts for selective alkene isomerization.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Modifying Interface Solvation and Oxygen Reduction Electrocatalysis with Hydrophobic Species

Ionic liquid interlayers improve the oxygen reduction reaction (ORR) kinetics on bulk and nanostructured catalysts for both Pt and alloyed-Pt materials. Despite the demonstrated performance enhancement at the half-cell and membrane electrode assembly level, the mechanism of the improvement is not fully understood. In this work, we combine single-crystal experiments with microkinetic modeling to uncover the origin of the ORR kinetic improvement on Pt(111) in the presence of ionic liquids and hydrophobic cations. With the incorporation of a modified Frumkin isotherm, our model accurately simulates the disorder–order transition observed in hydroxyl and bisulfate adsorption on Pt(111) under acidic conditions. Voltametric analysis shows that ionic liquids impact solvation to break so-called scaling relations between the adsorption strength of OH ad and O ad , but these effects have little impact on ORR activity. Instead, destabilized OH ad reduces the overall hydroxyl (spectator) coverage, resulting in higher availability of active sites.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗