Search NASA⌕ Search

SEARCH · Search NASA

Results for “Computational catalysis”

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 127 records · Page 7
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Application of machine learning to discover new intermetallic catalysts for the hydrogen evolution and the oxygen reduction reactions

The adsorption energies for hydrogen, oxygen, and hydroxyl were calculated by means of density functional theory on the lowest energy surface of 24 pure metals and 332 binary intermetallic compounds with stoichiometries AB, A 2 B, and A 3 B taking into account the effect of biaxial elastic strains. This information was used to train two random forest regression models, one for the hydrogen adsorption and another for the oxygen and hydroxyl adsorption, based on 9 descriptors that characterized the geometrical and chemical features of the adsorption site as well as the applied strain. All the descriptors for each compound in the models could be obtained from physico-chemical databases. The random forest models were used to predict the adsorption energy for hydrogen, oxygen, and hydroxyl of ≈2700 binary intermetallic compounds with stoichiometries AB, A 2 B, and A 3 B made of metallic elements, excluding those that were environmentally hazardous, radioactive, or toxic. This information was used to search for potential good catalysts for the HER and ORR from the criteria that their adsorption energy for H and O/OH, respectively, should be close to that of Pt. Further, this investigation shows that the suitably trained machine learning models can predict adsorption energies with an accuracy not far away from density functional theory calculations with minimum computational cost from descriptors that are readily available in physico-chemical databases for any compound. Moreover, the strategy presented in this paper can be easily extended to other compounds and catalytic reactions, and is expected to foster the use of ML methods in catalysis.

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

Polyketide synthase-based controlled synthesis of polycyclopropanated fuel molecules

Reducing carbon emissions from aviation and long-distance transportation sectors requires the development of sustainable biofuels with suitable energy density, freezing point, and other physical properties. We previously demonstrated biological production of high energy polycyclopropanated fatty acids (POP-FAs, class I) using an iterative polyketide synthase (iPKS) pathway in a Streptomyces host. Here, we used a computational model of fuel properties to identify chain length and cyclopropanation control as critical steps to engineer this iPKS for biofuel applications. We next explored the natural diversity of POP biosynthesis by investigating homologous pathways. Then, by in vivo gene exchange, we determined cyclopropanase (CP) catalysis to be key for POP-FA engineering. Leveraging both natural and engineered pathway product diversity, we demonstrate targeted production of improved POP-FAs, namely shortened POP-FAs with predicted superior freezing point properties for aviation, as well as fully cyclopropane-saturated POP-FAs which should have superior energy-density. These precise and controllable modifications to POP-FA structure open the door for bioproduction of designer POP fuels.

Yin, Kevin↗
DOE OSTI2025

High entropy powering green energy: hydrogen, batteries, electronics, and catalysis

A reformation in energy is underway to replace fossil fuels with renewable sources, driven by the development of new, robust, and multi-functional materials. High-entropy materials (HEMs) have emerged as promising candidates for various green energy applications, having unusual chemistries that give rise to remarkable functionalities. This review examines recent innovations in HEMs, focusing on hydrogen generation/storage, fuel cells, batteries, semiconductors/electronics, and catalysis—where HEMs have demonstrated the ability to outperform state-of-the-art materials. We present new master plots that illustrate the superior performance of HEMs compared to conventional systems for hydrogen generation/storage and heat-to-electricity conversion. We highlight the role of computational methods, such as density functional theory and machine learning, in accelerating the discovery and optimization of HEMs. The review also presents current challenges and proposes future directions for the field. We emphasize the need for continued integration of modeling, data, and experiments to investigate and leverage the underlying mechanisms of the HEMs that are powering progress in sustainable energy.

batteries↗
NASA NTRS1992

Stagnation-point heat-transfer rate predictions at aeroassist flight conditions

The results are presented for the stagnation-point heat-transfer rates used in the design process of the Aeroassist Flight Experiment (AFE) vehicle over its entire aeropass trajectory. The prediction methods used in this investigation demonstrate the application of computational fluid dynamics (CFD) techniques to a wide range of flight conditions and their usefulness in a design process. The heating rates were computed by a viscous-shock-layer (VSL) code at the lower altitudes and by a Navier-Stokes (N-S) code for the higher altitude cases. For both methods, finite-rate chemically reacting gas was considered, and a temperature-dependent wall-catalysis model was used. The wall temperature for each case was assumed to be radiative equilibrium temperature, based on total heating. The radiative heating was estimated by using a correlation equation. Wall slip was included in the N-S calculation method, and this method implicitly accounts for shock slip. The N-S/VSL combination of projection methods was established by comparison with the published benchmark flow-field code LAURA results at lower altitudes, and the direct simulation Monte Carlo results at higher altitude cases. To obtain the design heating rate over the entire forward face of the vehicle, a boundary-layer method (BLIMP code) that employs reacting chemistry and surface catalysis was used. The ratio of the VSL or N-S method prediction to that obtained from the boundary-layer method code at the stagnation point is used to define an adjustment factor, which accounts for the errors involved in using the boundary-layer method.

Gupta, Roop N.↗
DOE OSTI2025

Multimodal Ligand Non-Innocence for Reduction Catalysis (Final Technical Report)

The goal of this project was to develop multimodal ligand non-innocence in late transition metal complexes of aminophenolate ligands and to apply this ligand-centered reactivity to small-molecule reduction reactions relevant to energy-conversion chemistry. Specific aims were, 1) Prepare new coordination complexes of aminophenolate ligands and characterize them using a battery of structural, spectroscopic, and computational methods. 2) Benchmark the thermodynamics and kinetics of e - , H + , H • , and H - transfer reactivity to establish the baseline reactivity of new metal complexes containing aminophenolate ligands. 3) Apply new aminophenolate complexes to catalytic reduction reactions.

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

Accelerating Embedding Potential Optimization by Reconstructing the Pseudo-Valence Electron Density

Density functional embedding theory (DFET) enables use of electronic structure methods with higher accuracy than density functional theory in a local region, with applications thus far ranging from (photo/electro)catalysis to reactions in solution. DFET partitions a large collection of atoms into smaller groups that interact via a shared embedding (interaction) potential V emb , determined via functional optimization. The optimized effective potential (OEP) process used to optimize V emb is time-consuming and becomes a computational bottleneck due to sharp, oscillating features of V emb near nuclei. Here, similar to pseudopotential theory, by reconstructing electron densities used in the OEP process from smoother pseudo-valence-only (PVO) electron densities as proxies for total densities of the full system and subsystems, we can retain accuracy in the embedded electronic structure calculations while potentially reducing the overhead of V emb construction, within the projector augmented-wave (PAW) formalism. We explore three different chemical reactions as exemplars to test PVO–DFET, namely, H 2 dissociative adsorption on a Cu(111) surface, H 2 O adsorption on a Pt(111) surface, and aqueous [Ca 2+ –SO 4 2– ] ion-pair formation. The PVO approximation works well for all three systems with minimal loss of accuracy (∼10–70 meV error relative to the original exact-derivative (ED) approach) while accelerating V emb generation for the Cu and Pt systems respectively by 20× and 5×. Given proper numerical convergence parameters, the spatial distributions of differences between PVO- and ED-based V emb outside the core regions are small, explaining the exceptional agreement between the two approaches. Finally, we anticipate that this more efficient PVO–DFET approximation will be useful whenever computation of V emb is much more expensive than subsequent embedded high-level electron correlation calculations.

approximation↗
DOE OSTI2024

Unveiling Spatial and Temporal Dynamics of Plasmon-Enhanced Localized Fields in Metallic Nanoframes through Ultrafast Electron Microscopy

Plasmonic nanomaterials, particularly noble metal nanoframes (NFs), are important for applications such as catalysis, biosensing, and energy harvesting due to their ability to enhance localized electric fields and atomic efficiency via localized surface plasmon resonance (LSPR). Yet the fundamental structure-function relationships and plasmonic dynamics of the NFS are difficult to study experimentally and thus far rely predominately on computational methodologies, limiting their utilization. This study leverages the capabilities of ultrafast electron microscopy (UEM), specifically photon-induced near-field electron microscopy (PINEM), to probe the light-matter interactions within plasmonic NF structures. Here, the effects of shape, size, and plasmonic coupling of Pt@Au core-shell NFs on spatial and temporal characteristics of plasmon-enhanced localized electric fields are explored. Importantly, time-resolved PINEM analysis reveals that the plasmonic fields around hexagonal NF prisms exhibit a spatially dependent excitation and decay rate, indicating a nuanced interplay between the spatial geometry of the NF and the temporal evolution of the localized electric field. These results and observations uncover nanophotonic energy transfer dynamics in NFs and highlight their potential for applications in biosensing and photocatalysis.

Plasmonics↗
DOE OSTI2025

Off-Equilibrium Reactivity of Boron-Enriched Metal Diboride Surfaces in Electroreduction Conditions

Boron-based materials, featuring B-dependent reactivity and diverse phases, are emerging as promising catalyst systems. However, the catalytic mechanism on many borides remains poorly understood due to complex surface reconstructions under reaction conditions. Here, we investigate the MoB 2 surface in conditions of hydrogen evolution reaction in acidic media, using grand canonical global optimization, grand canonical density functional theory, ab initio molecular dynamics, free energy surface sampling, and an analytical model for electrochemical barrier evaluation. We propose a boron-enrichment strategy to tune the surface reactivity of the hexagonal face of MoB 2 . We reveal the dynamic nature of the B-enriched surface under H coverage and kinetic trapping of the system in the metastable regime with an extensive examination of the deactivation pathways. The metastable center B site on B-enriched surfaces, featuring buckled-up configuration and a usual relaxation effect, is found to be highly active toward HER via the Volmer–Heyrovsky mechanism. In conclusion, this work demonstrates how off-equilibrium behaviors can arise from the interplay between adsorbate coverage and surface reconstruction on a seemingly simple surface, and we present a theoretical framework and computational workflows to address these behaviors, along with other realistic complexities, in kinetics simulations.

Adsorption↗
DOE OSTI2024

Machine learning the electric field response of condensed phase systems using perturbed neural network potentials

Abstract The interaction of condensed phase systems with external electric fields is of major importance in a myriad of processes in nature and technology, ranging from the field-directed motion of cells (galvanotaxis), to geochemistry and the formation of ice phases on planets, to field-directed chemical catalysis and energy storage and conversion systems including supercapacitors, batteries and solar cells. Molecular simulation in the presence of electric fields would give important atomistic insight into these processes but applications of the most accurate methods such as ab-initio molecular dynamics (AIMD) are limited in scope by their computational expense. Here we introduce Perturbed Neural Network Potential Molecular Dynamics (PNNP MD) to push back the accessible time and length scales of such simulations. We demonstrate that important dielectric properties of liquid water including the field-induced relaxation dynamics, the dielectric constant and the field-dependent IR spectrum can be machine learned up to surprisingly high field strengths of about 0.2 V Å −1 without loss in accuracy when compared to ab-initio molecular dynamics. This is remarkable because, in contrast to most previous approaches, the two neural networks on which PNNP MD is based are exclusively trained on molecular configurations sampled from zero-field MD simulations, demonstrating that the networks not only interpolate but also reliably extrapolate the field response. PNNP MD is based on rigorous theory yet it is simple, general, modular, and systematically improvable allowing us to obtain atomistic insight into the interaction of a wide range of condensed phase systems with external electric fields.

Science & Technology - Other Topics↗
NASA NTRS1993

Analysis of hypersonic arcjet flow fields and surface heating of blunt bodies

A Gauss-Seidel implicit aerothermodynamic Navier-Stokes computational CFD code with thermochemical surface conditions (GIANTS) and with appropriate boundary conditions was developed for the analysis of hypersonic arcjet flows over large angle blunt bodies. Two-dimensional axisymmetric multicomponent full Navier-Stokes equations were solved, together with boundary conditions, to account for the surface slip, surface catalysis, and energy and mass conservation at the body surface. It is shown that the predictions for isothermal, noncatalytic, and no-slip surface conditions which are made using the GIANTS code with simple boundary conditions agree with those of the LAURA code, and the predicted normalized heating distributions over the forebody surface with boundary conditions for various cone angles are consistent with the trends of data from arcjet experiments.

Chen, Y.-K.↗
DOE OSTI2024

Mechanism of Ni-Catalyzed Photochemical Halogen Atom-Mediated C(sp 3 )–H Arylation

Within the context of Ni photoredox catalysis, halogen atom photoelimination from Ni has emerged as a fruitful strategy for enabling hydrogen atom transfer (HAT)-mediated C(sp 3 )–H functionalization. Despite the numerous synthetic transformations invoking this paradigm, a unified mechanistic hypothesis that is consistent with experimental findings on the catalytic systems and accounts for halogen radical formation and facile C(sp 2 )–C(sp 3 ) bond formation remains elusive. We employ kinetic analysis, organometallic synthesis, and computational investigations to decipher the mechanism of a prototypical Ni-catalyzed photochemical C(sp 3 )–H arylation reaction. Our findings revise the previous mechanistic proposals, first by examining the relevance of SET and EnT processes from Ni intermediates relevant to the HAT-based arylation reaction. Our investigation highlights the ability for blue light to promote efficient Ni–C(sp 2 ) bond homolysis from cationic Ni III and C(sp 2 )–C(sp 3 ) reductive elimination from bipyridine Ni II complexes. However interesting, the rates and selectivities of these processes do not account for the productive catalytic pathway. Instead, our studies support a mechanism that involves halogen atom evolution from in situ generated Ni II dihalide intermediates, radical capture by a Ni II (aryl)(halide) resting state, and key C–C bond formation from Ni III . Oxidative addition to Ni I , as opposed to Ni 0 , and rapid Ni III /Ni I comproportionation play key roles in this process. Furthermore, the findings presented herein offer fundamental insight into the reactivity of Ni in the broader context of catalysis.

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

Insights from designing an artificial cascade catalysis system using principles from substrate channeling in enzymes

Generalizing the key requirements of highly-selective, multi-step chemical conversions involving spatially separated reaction centers remains one of the grand challenges of chemistry. Much work towards this effort has focused on decomposing multi-step conversions into their constituent reactions, whose intermediates are successively upgraded in a chemical cascade via diffusion from center to center. This approach for synthesizing more complex molecules takes its cues from biochemical networks, where near-unit conversion of even complex carbohydrates is achieved by upgrading chemical precursors via enzymatic cascades. In this computational study we examine a simple cascade involving coupled Ag and Cu catalysts that sequentially converts CO2 to CO and then CO2 and CO to reduced products, generically named CO2Product and COProduct. The system architecture is inspired by the phenomenon of biological substrate channeling, and components are examined to evaluate their effects on conversion efficiency in the cascade. Aside from a substrate channel linking two reaction centers, we find efficient cascades must also incorporate directional substrate diffusion, compartmentalization of the reaction centers, and proper timing of substrate arrival at the active center. We make explicit linkages between these requirements and chemical conversion in known biological systems, revealing additional control elements that could be incorporated.

CO 2 reduction↗
DOE OSTI2026

Summary Report from the 2025 Interfaces for Energy and the Environment Conference

The inaugural Interfaces for Energy and the Environment Conference (IEEC) took place on May 19-23, 2025, at Pacific Northwest National Laboratory (PNNL), Richland, Washington (USA). The aim of this first-of-its-kind interdisciplinary meeting was to provide a forum for participants to share the latest cutting edge experimental and computational advances in interfacial science across energy and environmental applications. The sessions below (elaborated further in the report summaries) highlighted fundamental and applied collaborative research aimed at understanding the interactions occurring at interfaces in aqueous environments, including, but not limited to, the fields of geochemistry, atmospheric chemistry, agriculture, environmental management, and catalysis. They were organized to stimulate and provide opportunities to create, renew, and deepen collaborations. The conference included activities such as oral and poster presentations, honoree mentoring session, and a team building exercise to support all career stages (detailed summaries of these activities are in the Appendices).

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

Dynamics of McMillan mappings III. Symmetric map with mixed nonlinearity

This article extends the study of the dynamical properties of the symmetric McMillan map, emphasizing its utility in understanding and modeling complex nonlinear systems. Although the map features six parameters, we demonstrate that only two are irreducible: the linearized rotation number at the fixed point and a nonlinear parameter representing the ratio of terms in the biquadratic invariant. Through a detailed analysis, we classify regimes of stable motion, provide exact solutions to the mapping equations, and derive a canonical set of action-angle variables, offering analytical expressions for the rotation number and nonlinear tune shift. We further establish connections between general standard-form mappings and the symmetric McMillan map, using the area-preserving Hénon map and accelerator lattices with thin sextupole magnet as representative case studies. Our results show that, despite being a second-order approximation, the symmetric McMillan map provides a highly accurate depiction of dynamics across a wide range of system parameters, demonstrating its practical relevance in both theoretical and applied contexts.

43 PARTICLE ACCELERATORS↗
DOE OSTI2025

Barrierless nucleation in glassy precursors drives zeolite formation

Zeolites are crystalline, microporous silicates widely used in catalysis and separations, yet the molecular mechanisms of their formation remain unresolved. Experiments indicate that hydrothermal synthesis of silica zeolites from clear solution proceeds through amorphous nanoaggregates that gradually develop zeolite order in an apparently continuous amorphous to crystal transformation. Here, we combine molecular simulations with advanced algorithms that identify zeolite order and computer vision to elucidate the pathway from clear solution to zeolite nanocrystal. We show that at conditions of hydrothermal synthesis of silica zeolites, the transformation of precursor aggregate into zeolite is not limited by nucleation barriers but by the slow dynamics of reorganization in the glassy precursor matrix. The negligible nucleation barriers result in spinodal-like crystallization that leads to a gradual formation of a mosaic of small crystallites that explain the seemingly continuous character of zeolite crystallization and the catalytic activity of X-ray amorphous, protozeolites and embryonic zeolites. We find that zeolite-like porosity and short-range order emerge early within glassy precursors, well before crystallinity is detected in transmission electron microscopy (TEM) images or X-ray diffraction. The nanoaggregate’s temperature-size phase diagram reveals a convergence of the zeolite–amorphous equilibrium and maximum crystallization rate at ~3 nm diameter nanoparticle diameters and ~200 °C. This convergence signals the termination of the first-order amorphous-to-zeolite transition. Our results provide a unifying framework for understanding nucleation of silica zeolites from solution and suggest that barrierless nucleation may govern the formation of other nanoparticle systems, including minerals and oxides synthesized far below their bulk melting points.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Electron transfer in polysaccharide monooxygenase catalysis

Polysaccharide monooxygenase (PMO) catalysis involves the chemically difficult hydroxylation of unactivated C–H bonds in carbohydrates. The reaction requires reducing equivalents and will utilize either oxygen or hydrogen peroxide as a cosubstrate. Two key mechanistic questions are addressed here: 1) How does the enzyme regulate the timely and tightly controlled electron delivery to the mononuclear copper active site, especially when bound substrate occludes the active site? and 2) How does this electron delivery differ when utilizing oxygen or hydrogen peroxide as a cosubstrate? Using a computational approach, potential paths of electron transfer (ET) to the active site copper ion were identified in a representative AA9 family PMO from Myceliophthora thermophila ( Mt PMO9E). When Y62, a buried residue 12 Å from the active site, is mutated to F, lower activity is observed with O 2 . However, a WT-level activity is observed with H 2 O 2 as a cosubstrate indicating an important role in ET for O 2 activation. To better understand the structural effects of mutations to Y62 and axial copper ligand Y168, crystal structures were solved of the wild type Mt PMO9E and the variants Y62W, Y62F, and Y168F. A bioinformatic analysis revealed that position 62 is conserved as either Y or W in the AA9 family. The Mt PMO9E Y62W variant has restored activity with O 2 . Overall, the use of redox-active residues to supply electrons for the reaction with O 2 appears to be widespread in the AA9 family. Furthermore, the results provide a molecular framework to understand catalysis with O 2 versus H 2 O 2 .

Sayler, Richard I. (ORCID:000000017252707X)↗
DOE OSTI2025

Chemistry of Materials Underpinning Photoelectrochemical Solar Fuel Production

Since its inception, photoelectrochemistry has sought to power the generation of fuels, particularly hydrogen, using energy from sunlight. Efficient and durable photoelectrodes, however, remain elusive. Here we review the current state of the art, focusing our discussion on advances in photoelectrodes made in the past decade. We open by briefly discussing fundamental photoelectrochemical concepts and implications for photoelectrode function. We next review a broad range of semiconductor photoelectrodes broken down by material class (oxides, nitrides, chalcogenides, and mature photovoltaic semiconductors), identifying intrinsic properties and discussing their influence on performance. We then identify innovative in situ and operando techniques to directly probe the photoelectrode|electrolyte interface, enabling direct assessment of structure–property relationships for catalytic surfaces in active reaction environments. We close by considering more complex photoelectrochemical fuel-forming reactions (carbon dioxide and nitrogen reduction, as well as alternative oxidation reactions), where product selectivity imposes additional criteria on electrochemical driving force and photoelectrode architecture. By contextualizing recent literature within a fundamental framework, we seek to provide direction for continued progress toward achieving efficient and stable fuel-forming photoelectrodes.

08 HYDROGEN↗