Search NASA⌕ Search

SEARCH · Search NASA

Results for “computed reaction energetics”

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 19 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Connecting cation site location to alkane dehydrogenation activity in Ni/BEA catalysts

Ni-modified beta zeolite (Ni/BEA) catalysts activate carbon-hydrogen bonds in light alkanes, as demonstrated through isobutane reaction testing. Controlled synthesis of Ni/BEA allows for efficient introduction of ion-exchanged Ni sites at varying Ni loadings (0.43% - 1.8%). These catalysts exhibit site time yields (STY) for H2 production that increase with increasing Ni loading. A detailed analysis of secondary reactions and carbon deposition based on the relative molar flowrates of product C and H indicates that the observed increase in H2 STY with increasing Ni loading is likely attributed to both increasing alkane activation activity and increasing formation of hydrogen-deficient aromatic products retained within the catalyst pores. In situ diffuse-reflectance UV-visible-NIR absorbance and X-ray absorption spectroscopies indicate isolated, 4-coordinate Ni(2+) species for all loadings. Quantum mechanics/molecular mechanics modeling identifies two distinct Ni(2+) sites consistent with the structural characterization, but with differing relative stabilities due to their coordination environment. Computed reaction energetics for isobutane dehydrogenation demonstrate that the more stable Ni(2+) species at a six-membered 4Si-2Al ring, Ni-6MR, is less active for isobutane dehydrogenation than the less stable Ni(2+) at the five-membered 3Si-2Al ring, Ni-5MR. The differing local structure of the isolated cationic Ni sites in Ni/BEA offers a possible rationalization for the increased H2 STY observed at greater Ni loadings.

catalysts↗
DOE OSTI2021

BETO 2021 Peer Review - Electrocatalytic CO2 Utilization

The goal of the ChemCatBio DataHub project is to accelerate the catalyst and process development cycle by developing transformational tools for prediction and collaboration in catalyst R&D. The project is currently focused on the development of the Catalyst Property Database (CPD), a free and public resource released in September 2020. The CPD was designed to advance the state of the art for application of computational data. When computational data, such as computed reaction energetics, is used in catalyst design, it is almost always generated by the researchers seeking to use it, even if similar data has been published previously. One barrier to data reuse that results in this duplication of effort is the difficult process of finding and applying published data, which can be slow, error-prone, and manual. The CPD seeks to overcome these challenges by creating a centralized, searchable database of quality catalyst property data. At present, the CPD contains computed adsorption energies for intermediates along catalytic pathways. During FY21 and FY22, development of the CPD continues with a focus on external users and meeting their requirements. A batch upload capability, training and curation procedures, user interviews, and a demonstration of the CPD's utility in accelerating catalyst research are planned. Overall, the Data Hub project and CPD aim to reduce the time and cost of catalyst research by harnessing the power of data in catalyst discovery.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗
DOE OSTI2021

BETO 2021 Peer Review - ChemCatBio Data Hub 2.6.2.500

The goal of the ChemCatBio DataHub project is to accelerate the catalyst and process development cycle by developing transformational tools for prediction and collaboration in catalyst R&D. The project is currently focused on the development of the Catalyst Property Database (CPD), a free and public resource released in September 2020. The CPD was designed in response to the observation that when data, such as computed reaction energetics, is used in catalyst design, it is almost always generated by the researchers seeking to use it, even if similar data has been published previously. One barrier to data reuse that results in this duplication of effort is the difficult process of finding and applying published data, which can be slow, error-prone, and manual. The CPD seeks to overcome these challenges by creating a centralized, searchable database of quality catalyst property data. At present, the CPD contains computed adsorption energies for intermediates along catalytic pathways. During FY21 and FY22, development of the CPD continues with a focus on external users and meeting their requirements. A batch upload capability, training and curation procedures, user interviews, and a demonstration of the CPD's utility in accelerating catalyst research are planned. Overall, the Data Hub project and CPD aim to reduce the time and cost of catalyst research by harnessing the power of data in catalyst discovery.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗
DOE OSTI2022

Pathways for the Formation of C 2+ Products under Alkaline Conditions during the Electrochemical Reduction of CO 2

Establishing how Cu facilitates the electrochemical CO 2 reduction reaction (CO2RR) to C 2+ products remains a critical challenge. Under typical reaction conditions, the pH near the electrode is considerably more alkaline than that in the bulk due to mass transport limitations. Challenges with probing alkaline pathways using computational methods have limited understanding of the CO2RR under experimentally relevant conditions. In this work, using the Volmer reaction on Cu (100), we demonstrate that predicted activation barriers can substantially differ between acidic and alkaline pathways. We compute reaction energetics for alkaline *CO protonation and find that, while the formation of *CHO is preferred thermodynamically, the formation of *COH is favored kinetically at high overpotential. However, we find formation of *CHO via reaction of *H and *CO feasible at room temperature. Further, we report potential-dependent energetics for forming the first C-C bond in CO2RR and find that CO dimerization likely dominates. Finally, we investigate how long-range van der Waals interactions impact our results by comparing to the meta-GGA B97M-rV.

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

Constrained nuclear–electronic orbital method for periodic density functional theory: Application to H 2 chemisorption on Si(001) surfaces

The nuclear–electronic orbital (NEO) method provides a powerful computational framework for incorporating nuclear quantum effects (NQE) in electronic structure calculations beyond the Born–Oppenheimer approximation. By incorporating additional constraints to the position operator on quantum particles like protons, the NEO method enables calculation of effective potential that accounts for NQE. Here, in this work, we present a new constrained NEO (cNEO) formulation for density functional theory (cNEO-DFT) calculations in the context of extended periodic systems. Using the nudged elastic band method, we discuss an application of the cNEO-DFT approach to studying the adsorption of a hydrogen molecule on the Si(001) surfaces. The calculation shows how NQE impacts the reaction energetics. The proton density changes are computed along the reaction pathways. This work demonstrates the capability of the new cNEO-DFT method to study a wide range of chemical processes, such as surface reactions where the quantum nature of light atoms like protons is non-negligible.

Chemical processes↗
DOE OSTI2021

Hydrolysis of Small Oxo/Hydroxo Molecules Containing High Oxidation State Actinides (Th, Pa, U, Np, Pu): A Computational Study

The energetics of hydrolysis reactions for high oxidation states of oxo/hydroxo monomeric actinide species (Th IV O 2 , Pa IV O 2 , U IV O 2 , Pa V O 2 (OH), U V O 2 (OH), U VI O 3 , Np VI O 3 , Np VII O 3 (OH), and Pu VII O 3 (OH)) were calculated at the CCSD(T) level. The first step is the formation of a Lewis acid/base adduct with H 2 O (hydration), followed by a proton transfer to form a dihydroxide molecule (hydrolysis); this process is repeated until all oxo groups are hydrolyzed. The physisorption (hydration) for each H 2 O addition was predicted to be exothermic, ca. -20 kcal/mol. The hydrolysis products are preferred energetically over the hydration products for the +IV and +V oxidation states. The compounds with An VI are a turning point in terms of favoring hydration over hydrolysis. For An VII O 3 (OH), hydration products are preferred, and only two waters can bind; the complete hydrolysis process is now endothermic, and the oxidation state for the An in An(OH) 7 is +VI with two OH groups each having one-half an electron. The natural bond order charges and the reaction energies provide insights into the nature of the hydrolysis/hydration processes. The actinide charges and bond ionicity generally decrease across the period. Here, the ionic character decreases as the oxidation state and coordination number increase so that covalency increases moving to the right in the actinide period.

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

Description of reaction and vibrational energetics of CO 2 –NH 3 interaction using quantum computing algorithms

CO 2 capture is critical to solving global warming. Amine-based solvents are extensively used to chemically absorb CO 2 . Thus, it is crucial to study the chemical absorption of CO 2 by amine-based solvents to better understand and optimize CO 2 capture processes. Here, we use quantum computing algorithms to quantify molecular vibrational energies and reaction pathways between CO 2 and a simplified amine-based solvent model—NH 3 . Molecular vibrational properties are important to understanding kinetics of reactions. However, the molecule size correlates with the strength of anharmonicity effect on vibrational properties, which can be challenging to address using classical computing. Quantum computing can help enhance molecular vibrational calculations by including anharmonicity. We implement a variational quantum eigensolver (VQE) algorithm in a quantum simulator to calculate ground state vibrational energies of reactants and products of the CO 2 and NH 3 reaction. The VQE calculations yield ground vibrational energies of CO 2 and NH 3 with similar accuracy to classical computing. In the presence of hardware noise, Compact Heuristic for Chemistry (CHC) ansatz with shallower circuit depth performs better than Unitary Vibrational Coupled Cluster. The “Zero Noise Extrapolation” error-mitigation approach in combination with CHC ansatz improves the vibrational calculation accuracy. Excited vibrational states are accessed with quantum equation of motion method for CO 2 and NH 3 . Using quantum Hartree–Fock (HF) embedding algorithm to calculate electronic energies, the corresponding reaction profile compares favorably with Coupled Cluster Singles and Doubles while being more accurate than HF. Our research showcases quantum computing applications in the study of CO 2 capture reactions.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Prediction of Probabilistic Shock Initiation Thresholds of Energetic Materials Through Evolution of Thermal-Mechanical Dissipation and Reactive Heating

The ignition threshold of an energetic material (EM) quantifies the macroscopic conditions for the onset of self-sustaining chemical reactions. The threshold is an important theoretical and practical measure of material attributes that relate to safety and reliability. Historically, the thresholds are measured experimentally. In this work, we present a new Lagrangian computational framework for establishing the probabilistic ignition thresholds of heterogeneous EM out of the evolutions of coupled mechanical-thermal-chemical processes using mesoscale simulations. Furthermore, the simulations explicitly account for microstructural heterogeneities, constituent properties, and interfacial processes and capture processes responsible for the development of material damage and the formation of hotspots in which chemical reactions initiate. The specific mechanisms tracked include viscoelasticity, viscoplasticity, fracture, post-fracture contact, frictional heating, heat conduction, reactive chemical heating, gaseous product generation, and convective heat transfer. To determine the ignition threshold, the minimum macroscopic loading required to achieve self-sustaining chemical reactions with a rate of reactive heat generation exceeding the rate of heat loss due to conduction and other dissipative mechanisms is determined. Probabilistic quantification of the processes and the thresholds are obtained via the use of statistically equivalent microstructure sample sets (SEMSS). The predictions are in agreement with available experimental data.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Modernist materials synthesis: Finding thermodynamic shortcuts with hyperdimensional chemistry

Synthesis remains a challenge for advancing materials science. A key focus of this challenge is how to enable selective synthesis, particularly as it pertains to metastable materials. This perspective addresses the question: how can “spectator” elements, such as those found in double ion exchange (metathesis) reactions, enable selective materials synthesis? By observing reaction pathways as they happen (in situ) and calculating their energetics using modern computational thermodynamics, we observe transient, crystalline intermediates that suggest that many reactions attain a local thermodynamic equilibrium dictated by local chemical potentials far before achieving a global equilibrium set by the average composition. Furthermore, using this knowledge, one can thermodynamically “shortcut” unfavorable intermediates by including additional elements beyond those of the desired target, providing access to a greater number of intermediates with advantageous energetics and selective phase nucleation. Ultimately, data-driven modeling that unites first-principles approaches with experimental insights will refine the accuracy of emerging predictive retrosynthetic models for complex materials synthesis.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Solvation at metal/water interfaces: An ab initio molecular dynamics benchmark of common computational approaches

Determining the influence of the solvent on electrochemical reaction energetics is a central challenge in our understanding of electrochemical interfaces. To date, it is unclear how well existing methods predict solvation energies at solid/liquid interfaces, since they cannot be assessed experimentally. Ab initio molecular dynamics (AIMD) simulations present a physically highly accurate, but also a very costly approach. In this work, we employ extensive AIMD simulations to benchmark solvation at charge-neutral metal/water interfaces against commonly applied continuum solvent models. We consider a variety of adsorbates including *CO, *CHO, *COH, *OCCHO, *OH, and *OOH on Cu, Au, and Pt facets solvated by water. The surfaces and adsorbates considered are relevant, among other reactions, to electrochemical CO 2 reduction and the oxygen redox reactions. We determine directional hydrogen bonds and steric water competition to be critical for a correct description of solvation at the metal/water interfaces. As a consequence, we find that the most frequently applied continuum solvation methods, which do not yet capture these properties, do not presently provide more accurate energetics over simulations in vacuum. We find most of the computed benchmark solvation energies to linearly scale with hydrogen bonding or competitive water adsorption, which strongly differ across surfaces. Furthermore, we determine solvation energies of adsorbates to be non-transferable between metal surfaces, in contrast to standard practice.

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

Enhancing the connection between computation and experiments in electrocatalysis

Combining computational and experimental methods is a powerful approach to understand the variables that govern catalyst performance and ultimately design improved materials. However, the effectiveness of this approach rests on the strength of the relationships between calculated parameters and experimental measurements. These relationships are complicated by the intricacy and dynamic behaviour of catalytic active sites, and by the non-trivial relationship between calculated reaction energetics and observed rates. In this Perspective, we highlight opportunities to enhance the connection between computation and experiment in electrocatalysis. These include measuring the intrinsic kinetic behaviour of catalysts, creating precise models for the active site and its environment, and forming clear relationships between calculated reaction energetics and observed rates. As experimental and computational methods continue to become more powerful, clear connections between the two will maximize their utility to guide the design of efficient and selective electrocatalysts.

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

Stability and Activity of Cobalt Antimonate for Oxygen Reduction in Strong Acid

Guided by computational Pourbaix screening and high-throughput experiments aimed at the development of precious-metal-free fuel cells, we investigate rutile CoSb 2 O 6 as an electrocatalyst for oxygen reduction in 1 M sulfuric acid. Following 4 h of catalyst conditioning at 0.7 V vs RHE, operation at this potential for 20 h yielded an average current density of –0.17 mA cm –2 with corrosion at a rate of 0.04 nm hour –1 that is stoichiometric with catalyst composition. Surface Pourbaix analysis of the (111) surface identified partial H coverage under operating conditions. The Sb active site has an HO* binding free energy of 0.49 eV, which is near the peak of the kinetic 4e – ORR volcano for transition-metal oxides in acidic conditions. Furthermore, the experimental demonstration of operational stability and computational identification of a reaction pathway with favorable energetics place rutile CoSb 2 O 6 among the most promising precious-metal-free electrocatalysts for oxygen reduction in acidic media.

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

Neural Network Sampling of the Free Energy Landscape for Nitrogen Dissociation on Ruthenium

In heterogeneous catalysis, free energy profiles of reactions govern the mechanisms, rates, and equilibria. Energetics are conventionally computed using the harmonic approximation (HA), which requires determination of critical states a priori. Here, we use neural networks to efficiently sample and directly calculate the free energy surface (FES) of a prototypical heterogeneous catalysis reaction—the dissociation of molecular nitrogen on ruthenium—at density-functional-theory-level accuracy. We find that the vibrational entropy of surface atoms, often neglected in HA for transition metal catalysts, contributes significantly to the reaction barrier. The minimum free energy path for dissociation reveals an “on-top” adsorbed molecular state prior to the transition state. While a previously reported flat-lying molecular metastable state can be identified in the potential energy surface, it is absent in the FES at relevant reaction temperatures. Finally, these findings demonstrate the importance of identifying critical points self-consistently on the FES for reactions that involve considerable entropic effects.

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

Extension of Solvent-Assisted Linker Exchange to Supported Metal–Organic Framework Thin Films

Solvent-assisted linker exchange (SALE) is a common postsynthetic approach used to accommodate metal–organic frameworks (MOFs) for targeted applications. As this technique is relatively untested in immobilized MOF structures, we repeated the generalized SALE approach with ZIF-8 deposited directly onto solid substrates and assessed the results with a variety of techniques, including time-of-flight secondary-ion mass spectrometry (ToF-SIMS). Our methods confirm that the supported films maintain their key structural properties while allowing for linker substitution with foreign species. Alongside experiments, we used computational modeling to quantify the relative energetics of SALE reactions with a variety of linkers. These techniques together reveal that SALE is translatable to immobilized MOFs. Additionally, we recognize how various aspects of the postsynthetic approach can bear consequences for the material’s final micro- and macroscopic properties. Furthermore, the results uncover basic yet critical details needed to inform future efforts to engineer MOF films with functional properties.

Ligands↗
DOE OSTI2022

Computational Aspects of Single-Molecule Kinetics for Coupled Catalytic Cycles: A Spectral Analysis

Catalysis from single active sites is analyzed using methods developed from single molecule kinetics. Using a stochastic Markov state description, the observable properties of general catalytic networks of reactions are expressed using an eigenvalue decomposition of the transition matrix for the Markov process. By the use of a sensitivity analysis, the necessary eigenvalues and eigenvectors are related to the energies of controlling barriers and wells located along the reaction routes. A generalization of the energetic span theory allows the eigenvalues to be computed from several activation energies corresponding to distinct barrier-well pairings. The formalism is demonstrated accurately for model problems for a physically realistic mechanism for an alkene hydrogenation reaction on a single atom catalyst.

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

Extemporaneous Mechanochemistry: Shock-Wave-Induced Ultrafast Chemical Reactions Due to Intramolecular Strain Energy

We report regions of energy localization referred to as hotspots are known to govern shock initiation and the run-to-detonation in energetic materials. Mounting computational evidence points to accelerated chemistry in hotspots from large intramolecular strains induced via the interactions between the shock wave and microstructure. However, definite evidence mapping intramolecular strain to accelerated or altered chemical reactions has so far been elusive. From a large-scale reactive molecular dynamics simulation of the energetic material 1,3,5-triamino-2,4,6-trinitrobenzene, we map decomposition kinetics to molecular temperature and intramolecular strain energy prior to reaction. Both temperature and intramolecular strain are shown to accelerate chemical kinetics. A detailed analysis of the atomistic trajectory shows that intramolecular strain can induce a mechanochemical alteration of decomposition mechanisms. The results in this paper could inform continuum-level chemistry models to account for a wide range of mechanochemical effects.

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

Exploring the Effects of Node Topology, Connectivity, and Metal Identity on the Binding of Nerve Agents and Their Hydrolysis Products in Metal–Organic Frameworks

Recent studies have shown that metal–organic frameworks (MOFs) built from hexanuclear M(IV) oxide cluster nodes are effective catalysts for nerve agent hydrolysis, where the properties of the active sites on the nodes can strongly influence the reaction energetics. The connectivity and metal identity of these M6 nodes can be easily tuned, offering extensive opportunities for computational screening to predict promising new materials. Thus, we used density functional theory (DFT) to examine the effects of node topology, connectivity, and metal identity on the binding energies of multiple nerve agents and their corresponding hydrolysis products. By computing an optimization metric based on the relative binding strengths of key hydrolysis reaction species (water, agent, and bidentate-bound products), we predicted optimal M6 nodes for hydrolyzing specific nerve agent and simulant molecules, where our results are in qualitative agreement with observed experimental trends. This analysis highlighted the notion that no single metal or node topology is optimal for all possible organophosphates, suggesting that MOFs should be selected based on the agent of interest. Using the large amount of data generated from our DFT calculations, we then derived quantitative structure–activity relationship (QSAR) models to help explain the complex trends observed in the binding energies. Through linear regression, we identified the most important descriptors for describing the binding of nerve agents and their hydrolysis products to M6 nodes. These results suggested that both molecular and node properties, including both structural and chemical features, collectively contribute to the binding energetics. By performing a thorough statistical analysis, we showed that our QSAR models are capable of making quantitatively accurate binding energy predictions for nerve agents and their hydrolysis products in a wide variety of M(IV)-MOFs. In conclusion, the insights gained herein can be used to guide future experiments for the synthesis of MOFs with enhanced catalytic activity for organophosphate hydrolysis.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗