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 145 records · Page 8
All sourcesNASA NTRSDOE OSTI
DOE OSTI2018

Selective Benzyl Alcohol Oxidation over Pd Catalysts

In the last decades, the selective liquid phase oxidation of alcohols to the corresponding carbonyl compounds has been a subject of growing interest. Research has focused on green methods that use “clean” oxidants such as O 2 in combination with supported metal nanoparticles as the catalyst. Among the alcohols, benzyl alcohol is one of the most studied substrates. Indeed, benzyl alcohol can be converted to benzaldehyde, largely for use in the pharmaceutical and agricultural industries. This conversion serves as model reaction in testing new potential catalysts, that can then be applied to other systems. Pd based catalysts have been extensively studied as active catalytic metals for alcohol oxidation for their high activity and selectivity to the corresponding aldehyde. Several catalytic materials obtained by careful control of the morphology of Pd nanoparticles, (including bimetallic systems) and by tuning the support properties have been developed. Moreover, reaction conditions, including solvent, temperature, pressure and alcohol concentration have been investigated to tune the selectivity to the desired products. Different reaction mechanisms and microkinetic models have been proposed. The aim of this review is to provide a critical description of the recent advances on Pd catalyzed benzyl alcohol oxidation.

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

Theoretical Modeling of Electrochemical Proton-Coupled Electron Transfer

Proton-coupled electron transfer (PCET) plays an essential role in a wide range of electrocatalytic processes. Here, a vast array of theoretical and computational methods have been developed to study electrochemical PCET. These methods can be used to calculate redox potentials and pKa values for molecular electrocatalysts, proton-coupled redox potentials and bond dissociation free energies for PCET at metal and semiconductor interfaces, and reorganization energies associated with electrochemical PCET. Periodic density functional theory can also be used to compute PCET activation energies and perform molecular dynamics simulations of electrochemical interfaces. Various approaches for maintaining a constant electrode potential in electronic structure calculations and modeling complex interactions in the electric double layer (EDL) have been developed. Theoretical formulations for both homogeneous and heterogeneous electrochemical PCET spanning the adiabatic, nonadiabatic, and solvent-controlled regimes have been developed and provide analytical expressions for the rate constants and current densities as functions of applied potential. The quantum mechanical treatment of the proton and inclusion of excited vibronic states have been shown to be critical for describing experimental data, such as Tafel slopes and potential dependent kinetic isotope effects. The calculated rate constants can be used as input to microkinetic models and voltammogram simulations to elucidate complex electrocatalytic processes.

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

Toward Carbon Monoxide Methanation at Mild Conditions on Dual-Site Catalysts

The catalytic carbon monoxide (CO) methanation is an ideal model reaction for the fundamental understanding of catalysis on the gas–solid interface and is crucial for various industrial processes. However, the harsh operating conditions make the reaction unsustainable, and the limitations set by the scaling relations between the dissociation energy barrier and dissociative binding energy of CO further increase the difficulty in designing high-performance methanation catalysts operating under milder conditions. In this work, we proposed a theoretical strategy to circumvent the limitations elegantly and achieve both facile CO dissociation and C/O hydrogenation on the catalyst containing a confined dual site. The DFT-based microkinetic modeling (MKM) reveals that the designed Co-Cr 2 /G dual-site catalyst could provide 4–6 orders of magnitude higher turnover frequency for CH 4 production than the cobalt step sites. We believe that the proposed strategy in the current work will provide essential guidance for designing state-of-the-art methanation catalysts under mild conditions.

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

How the Bioinspired Fe 2 Mo 6 S 8 Chevrel Breaks Electrocatalytic Nitrogen Reduction Scaling Relations

The nitrogen reduction reaction (NRR) is a renewable alternative to the energy- and CO 2 -intensive Haber–Bosch NH 3 synthesis process but is severely limited by the low activity and selectivity of studied electrocatalysts. The Chevrel phase Fe 2 Mo 6 S 8 has a surface Fe–S–Mo coordination environment that mimics the nitrogenase FeMo-cofactor and was recently shown to provide state-of-the-art activity and selectivity for NRR. Here, we elucidate the previously unknown NRR mechanism on Fe 2 Mo 6 S 8 via grand-canonical density functional theory (GC-DFT) that realistically models solvated and biased surfaces. Fe sites of Fe 2 Mo 6 S 8 selectively stabilize the key *NNH intermediate via a narrow band of free-atom-like surface d-states that selectively hybridize with p-states of *NNH, which results in Fe sites breaking NRR scaling relationships. These sharp d-states arise from an Fe–S bond dissociation during N 2 adsorption that mimics the mechanism of the nitrogenase FeMo-cofactor. Furthermore, we developed a new GC-DFT-based approach for calculating transition states as a function of bias (GC-NEB) and applied it to produce a microkinetic model for NRR at Fe 2 Mo 6 S 8 that predicts high activity and selectivity, in close agreement with experiments. Furthermore, our results suggest new design principles that may identify effective NRR electrocatalysts that minimize the barriers for *N 2 protonation and *NH 3 desorption and that may be broadly applied to the rational discovery of stable, multinary electrocatalysts for other reactions where narrow bands of surface d-states can be tuned to selectively stabilize key reaction intermediates and guide selectivity toward a target product. Furthermore, our results highlight the importance of using GC-DFT and GC-NEB to accurately model electrocatalytic reactions.

10 SYNTHETIC FUELS↗
DOE OSTI2022

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↗
DOE OSTI2023

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↗
DOE OSTI2021

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↗
DOE OSTI2025

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↗
DOE OSTI2018

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↗
DOE OSTI2020

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↗
DOE OSTI2023

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↗
DOE OSTI2019

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↗
DOE OSTI2022

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↗
DOE OSTI2017

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↗
DOE OSTI2024

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↗
DOE OSTI2018

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↗
DOE OSTI2021

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↗