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 37 records · Page 2
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Frontiers in divalent $f$-block chemistry

Divalent f-block chemistry has undergone rapid expansion in recent years, driven by advances in the stabilization of low-valent lanthanide and actinide complexes. Although f-elements were historically accessed in the +3 or higher oxidation states, the isolation of +2 species has revealed unusual spectroscopic signatures, distinct bonding motifs, and multiple accessible electronic configurations that influence their chemical and physical properties. These developments have generated opportunities in areas including quantum information science, molecular magnetism, catalysis, and luminescence. Computational chemistry has played a pivotal role in interpreting the electronic structure and reactivity of these systems. However, accurately modeling divalent f-block complexes remains challenging because of strong electronic correlation, multiconfigurational character, and the presence of close-in-energy competing electronic states. As experimental capabilities and theoretical methodologies continue to advance, a comprehensive assessment of divalent chemistry is both timely and needed. In this Review, we integrate experimental and theoretical perspectives to provide a systematic analysis of divalent lanthanide and actinide complexes across the f-block series. We critically assess the role of ligands in determining competing ground-state configurations (f n d 1 vs. f n+1 ) resulting in their distinct spectroscopic, magnetic, and bonding properties. Finally, we discuss emerging strategies for predictive modeling and rational design of low-valent f-block molecular systems, with potential applications ranging from dinitrogen reduction to quantum technologies.

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

Tuning Optical and Electrical Properties of Vanadium Oxide with Topochemical Reduction and Substitutional Tin

Vanadium oxides are widely tunable materials, with many thermodynamically stable phases suitable for applications spanning catalysis to neuromorphic computing. The stability of vanadium in a range of oxidation states enables mixed-valence polymorphs of kinetically accessible metastable materials. Low-temperature synthetic routes to, and the properties of, these metastable materials are poorly understood and may unlock new optoelectronic and magnetic functionalities for expanded applications. In this work, we demonstrate topochemical reduction of α-V 2 O 5 to produce metastable vanadium oxide phases with tunable oxygen vacancies (>6%) and simultaneous substitutional tin incorporation (>3.5%). The chemistry is carried out at low temperature (65 °C) with solution-phase SnCl 2 , where Sn 2+ is oxidized to Sn 4+ as V 5+ sites are reduced to V 4+ during oxygen vacancy formation. Despite high oxygen vacancy and tin concentrations, the transformations are topochemical in that the symmetry of the parent crystal remains intact, although the unit cell expands. Band structure calculations show that these vacancies contribute electrons to the lattice, whereas substitutional tin contributes holes, yielding a compensation doping effect and control over the electronic properties. The SnCl 2 redox chemistry is effective on both solution-processed V 2 O 5 nanoparticle inks and mesoporous films cast from untreated inks, enabling versatile routes toward functional films with tunable optical and electronic properties. The electrical conductance rises concomitantly with the SnCl 2 concentration and treatment time, indicating a net increase in density of free electrons in the host lattice. This work provides a valuable demonstration of kinetic tailoring of electronic properties of vanadium–oxygen systems through top-down chemical manipulation from known thermodynamic phases.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Diffusion Model-Guided Inverse Design of Bimetallic Catalysts for Ammonia Decomposition

In the past decade, artificial intelligence and deep learning have played increasingly prominent roles in materials design and discovery. Among these, generative AI models, known for their ability to create unique and complex structures, have emerged as state-of-the-art tools for materials screening due to their high efficiency and low computational cost. In catalysis, one of the major challenges is identifying promising material candidates within an immense chemical space. This challenge can be addressed using generative approaches, such as diffusion-based inverse design models. In this study, we present a machine learning-guided workflow that employed a diffusion model for the inverse design of bimetallic alloy catalysts for low-carbon ammonia decomposition, a key reaction for ammonia emission control and sustainable hydrogen production. Catalyst candidates were evaluated using nitrogen adsorption energy as the key descriptor, inspired by multiscale modeling. The proposed workflow identified low-cost, environmentally friendly catalysts with excellent catalytic performance, which have been validated theoretically and experimentally. Our framework decoupled the generative and property-prediction components, enhancing both flexibility and accuracy in the catalytic material design process.

Adsorption↗
DOE OSTI2026

Accelerating catalytic advancements through the precision of high-throughput experiments & calculations

The growing demand for energy-efficient processes to support a sustainable future drives the need for research to rapidly explore chemical and material space through accelerated catalyst discovery initiatives. Recent breakthroughs in high-throughput experimental and computational methods are transforming the catalysis field, surpassing traditional approaches to manipulating variables in catalytic processes. Key advancements in innovation include the integration of machine learning for efficient catalyst screening, high-throughput experimentation, data-driven methodologies employing comprehensive databases, and in situ and in operando techniques for realistic observations. This progress has undoubtedly been intertwined with a collaborative framework across disciplines, reshaping catalyst discovery methods in both industry and academia. This Opinion article presents a multifaceted perspective from coauthors with expertise spanning various stages of the Technology Readiness Level spectrum, highlighting both opportunities and persistent challenges in integrating computational and experimental approaches in catalysis. These challenges span from obtaining high-quality experimental data, scaling simulations to industrially relevant materials and process conditions to navigating the complexity and predictive accuracy of computational models.

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

Nanocluster Active Sites Formed on Heterogeneous Thermal Catalysts and Electrocatalysts by Operando Reactive Environments

Here, in this Viewpoint, we summarize our recent studies in this endeavor, ending with a forward-looking perspective into the remaining challenges and open questions for consideration by the catalysis community. Through a computationally efficient framework of systematically investigating metal atom ejection from and migration onto metal surfaces, leading to the formation of metal atom clusters, we show that the surface of a metal catalyst is dynamic and much less rigid than previously thought. We suggest that the operando formed metal atom clusters are novel active sites, which can dominate the activity of steps or kinks (let aside terraces) from where their constituent atoms were ejected. Through more readily responding to reaction conditions by forming catalytically active clusters, solid metal catalysts approach not only homogeneous catalyst motifs, thus establishing a bridge between homogeneous and heterogeneous catalysis, but also liquid metals which have been suggested to be better than traditional heterogeneous catalysts.

adsorption↗
DOE OSTI2026

Mechanistic Studies of an Iron-Catalyzed Intermolecular C–H Amination Reaction under Catalytic Conditions and Having a Large KIE

The conversion of C–H bonds into amines by nitrene insertion is an attractive transformation since it is both atom- and step-economical, and provides a direct route to functionalizing hydrocarbons. Using an iron catalyst [{( tBu pyrr) 2 pyr}Fe(OEt 2 )] (1-OEt 2 ) (( tBu pyrr) 2 pyr 2– = 3,5- t Bu 2 -bis(pyrrolyl)pyridine), we recently demonstrated the catalytic conversion of weak C–H bonds into secondary amines using aryl azides as the nitrene source [Zars, E.; Angew. Chem., Int. Ed. 2023, 62, e202311749]. Here, we describe detailed mechanistic studies of this intermolecular C–H amination reaction under catalytic conditions. We find by Variable Time Normalization Analysis (VTNA) that the conversion of xanthene (2-H 2 ) and 2,4,6-trimethyl-phenyl azide ( Me 3) catalyzed by 1-OEt 2 is an overall 3/2 order process, being 1 st order in 2-H 2 and half order in Me 3. A kinetic isotope effect study (KIE) using 2-d 2 results in a significant decrease in the rate (KIE = 61(15)), which clearly implicates the C–H insertion step as rate-determining. Furthermore, treatment of 1-OEt 2 with one equivalent of N 3 -2,6- i Pr 2 –C 6 H 3 yields the mixed-valence C–N coupled product [( tBu pyrr) 2 pyrFe-N═C(2,6 i Pr 2 –Ph)═N-(2,6 i Pr 2 –Ph))Fe tBu pyrrpyr(2-H-pyrr)] (5 iPr ). Quantum chemical calculations confirm the electronic structure of the mixed-valence dimer in 5 iPr and rationalize the Hammett correlation by a delicate balance in the dinuclearization of the catalytically active monomers. Calculations further indicate significant tunneling for the pivotal H atom abstraction by the iron-imidyl complex. Combining all these results allows us to propose a mechanism consisting of imido formation in equilibrium with a radical-coupled diiron system, followed by stepwise C–H insertion via a linear H atom abstraction transition state and subsequent radical rebound.

azides↗
DOE OSTI2025

Nature of molybdenum carbide surfaces for catalytic hydrogen dissociation using machine-learned potentials: an ensemble-averaged perspective

Molybdenum carbides with an electronic structure similar to noble metals have gained attention as a promising low-cost catalyst for biomass valorization and the hydrogen evolution reaction. However, our fundamental understanding of the catalyst surface and how different phases of these catalysts behave at varying reaction conditions is limited to ground state density functional theory calculations as ab initio molecular dynamics (AIMD) is computationally prohibitive at relevant length and time scales. Here, in this work, we train a multi-atomic cluster expansion (MACE) machine-learned interatomic potentials (MLIP) to study hydrogen dissociation and dynamics over Mo, δ-MoC, α-Mo 2 C, and β-Mo 2 C surfaces at varying temperatures and hydrogen partial pressures. Our simulations identify unique and different molecular and atomic hydrogen adsorption sites on different surfaces that do not depend on the temperature. At low hydrogen pressures, the surface coverage is monolayer, which transitions to two-layer adsorption at higher pressures. We find that atomic hydrogen diffusion and recombinations are preferred over molybdenum atom hollow sites, while the diffusion over carbon-terminated facets was negligible, signifying particularly strong C–H interactions. In contrast, molecular hydrogen adsorption occurs mostly atop Mo or the bridging sites. At a comparable hydrogen loading, β-Mo 2 C (001) is the most active surface for hydrogen dissociation reaction. This work provides insights into the dynamic nature of the hydrogen dissociation chemistry and the diversity of hydrogen adsorption sites on molybdenum carbides.

08 HYDROGEN↗
NASA NTRS2001

Aerothermal Analysis of the Project Fire II Afterbody Flow

Computational fluid dynamics (CFD) is used to simulate the wake flow and afterbody heating of the Project Fire II ballistic reentry to Earth at 11.4 km/sec. Laminar results are obtained over a portion of the trajectory between the initial heat pulse and peak afterbody heating. Although non-catalytic forebody convective heating results are in excellent agreement with previous computations, initial predictions of afterbody heating were about a factor of two below the experimental values. Further analysis suggests that significant catalysis may be occurring on the afterbody heat shield. Computations including finite-rate catalysis on the afterbody surface are in good agreement with the data over the early portion of the trajectory, but are conservative near the peak afterbody heating point, especially on the rear portion of the conical frustum. Further analysis of the flight data from Fire II shows that peak afterbody heating occurs before peak forebody heating, a result that contradicts computations and flight data from other entry vehicles. This result suggests that another mechanism, possibly pyrolysis, may be occurring during the later portion of the trajectory, resulting in less total heat transfer than the current predictions.

Wright, Michael J.↗
NASA NTRS1987

Computational chemistry

With the advent of supercomputers, modern computational chemistry algorithms and codes, a powerful tool was created to help fill NASA's continuing need for information on the properties of matter in hostile or unusual environments. Computational resources provided under the National Aerodynamics Simulator (NAS) program were a cornerstone for recent advancements in this field. Properties of gases, materials, and their interactions can be determined from solutions of the governing equations. In the case of gases, for example, radiative transition probabilites per particle, bond-dissociation energies, and rates of simple chemical reactions can be determined computationally as reliably as from experiment. The data are proving to be quite valuable in providing inputs to real-gas flow simulation codes used to compute aerothermodynamic loads on NASA's aeroassist orbital transfer vehicles and a host of problems related to the National Aerospace Plane Program. Although more approximate, similar solutions can be obtained for ensembles of atoms simulating small particles of materials with and without the presence of gases. Computational chemistry has application in studying catalysis, properties of polymers, all of interest to various NASA missions, including those previously mentioned. In addition to discussing these applications of computational chemistry within NASA, the governing equations and the need for supercomputers for their solution is outlined.

Arnold, J. O.↗
DOE OSTI2025

Application of machine learning interatomic potentials in heterogeneous catalysis

Heterogeneous catalysts are crucial in modern societies as they promote sustainability by enabling lower-energy pathways for various chemical reactions. While Density Functional Theory (DFT) computations can provide critical insights into how heterogeneous catalysts operate at the atomic level, they are limited by computational costs and unfavorable scaling with system size. Recently, machine learning interatomic potentials (MLIPs) have emerged as a promising alternative to DFT, offering near-DFT accuracy at significantly reduced cost. Here, in this perspective, we discuss the application of MLIPs in heterogeneous catalyst modeling as a surrogate for DFT. We detail how MLIPs have been applied in thermal catalysis to probe active sites, enable studying complex metallic and nanoporous catalysts, and investigate the reconstruction of catalytic surfaces. We review the use of MLIPs in electrocatalysis and photocatalysis, emphasizing their capabilities in studying transition metal oxide surfaces and solid–liquid interfaces. We also discuss the current limitations of MLIPs, particularly their challenges with transferability and description of non-local interactions. Finally, we conclude by identifying promising and underexplored domains in which MLIPs can further advance our understanding of heterogeneous catalysts.

Catalytic surfaces↗
DOE OSTI2025

Transferable dispersion-aware machine learning interatomic potentials for multilayer transition metal dichalcogenide heterostructures

Stacking atomically thin transition metal dichalcogenides (TMDs) into heterostructures enables exploration of exotic quantum phases, particularly through twist-angle-controlled moiré superlattices. These structures exhibit novel electronic and optical behaviors driven by atomic-scale structural reconstruction. However, studying such systems with DFT is computationally demanding due to their large unit cells and van der Waals (vdW) interactions between layers. To address this, we develop a transferable neural network potential (NNP) that includes long-range vdW corrections up to 12Å with minimal overhead. Trained on vdW-corrected DFT data for Mo- and W-based TMDs with S, Se, and Te, the NNP accurately models monolayers, bilayers, heterostructures, and their interaction with h-BN substrates. It reproduces equilibrium structures, energy landscapes, phonon dispersions, and matches experimental atomic reconstructions in twisted WS2 and MoS2/WSe2 systems. We demonstrate that our NNP achieves DFT-level accuracy and high computational efficiency, enabling large-scale simulations of TMD-based moiré superlattices both with and without substrates.

materials for energy and catalysis↗
DOE OSTI2026

Roadmap for transforming heterogeneous catalysis with artificial intelligence

Artificial intelligence (AI) is poised to transform heterogeneous catalysis, opening avenues for catalytic materials discovery. By uncovering intricate patterns in high-dimensional data, AI has been reshaping our pursuit of sustainable catalytic processes across the energy, environmental and chemical sectors. This promise, however, hinges on overcoming fundamental barriers, including limitations in data availability and quality, challenges in the generalizability and interpretability of data-augmented decisions, and the persistent gap between in silico predictions and experiments. Furthermore, we outline a forward-looking roadmap for deeply integrating AI into heterogeneous catalysis with an AI-ready data ecosystem, multimodal foundation models, and ultimately autonomous laboratories to accelerate the development of next-generation catalytic technologies via AI-empowered human–machine collaboration.

Computational methods↗
DOE OSTI2026

Tandem bulk oxygen diffusion and surface reactions in reducible metal oxides control redox cycle dynamics

The interplay between bulk oxygen diffusion and surface reactions in reducible metal oxides is key in heterogeneous catalysts, but direct measurements of oxygen mobility, transient kinetics, and in situ spectroscopies have been lacking. Here, we reveal complex dynamic behavior of ceria-zirconia by H 2 using transient kinetics via mass spectrometry and in situ Raman and near-ambient pressure x-ray photoelectron spectroscopies. Molecular dynamics simulations with a machine learning potential delineate competitive oxygen diffusion mechanisms, with an optimal mobility at intermediate reductions. We expose a compensation between vacancy availability and lattice distortion at intermediate to high reductions and Frenkel defects at low reductions, underscoring a potential deficiency of 16 O/ 18 O exchange experiments in deducing oxygen mobility. Vacancies in proximity require electron localization on Ce atoms further away. The continuous replenishment of surface oxygen results in a varying reduction rate, with H 2 dissociation being the rate-limiting step. Multiscale transient simulations, consistent with experiments, indicate catalysts of potentially spatially varying oxidation states. The approach is broadly applicable to reducible oxide materials.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Finite-element-based simulations of electrodes for CO 2 cascade reduction reactions

The multielectron reduction of CO 2 to liquid fuels could be a path to scalable energy storage, but reaching this goal requires major advances in catalysis and systems engineering. Cascade catalysis, which couples sequential reactions without isolating intermediates, has emerged as a promising route to enhance selectivity and efficiency in CO 2 reduction (CO 2 R). In this review, we examine how finite-element-based simulations of continuum model [finite element method (FEM)] approaches are being used to analyze and guide CO 2 R cascade systems. We first outline the fundamentals of cascade catalysis and recent advances in catalytic materials (metallic, molecular, and hybrid architectures). We then focus on FEM developments at the electrode and device scales, emphasizing how these models capture transport phenomena, local microenvironments, and geometry-dependent effects. To clarify design principles, we present case studies of cascade electrodes organized in systems without and with integrated semiconductors. We further emphasize the integration of FEM with multiscale frameworks (density functional theory, molecular dynamics, kinetic Monte Carlo) and its role in bridging atomic-level insights with device-level performance. Finally, we identify current limitations and future prospects, including improved boundary conditions, coupling with operando experiments, and machine learning-accelerated model development. Together, these insights provide design principles for next-generation CO 2 R cascade systems for efficient solar fuel production.

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

Static Subspace Approximation for Random Phase Approximation Correlation Energies: Implementation and Performance

Developing theoretical understanding of complex reactions and processes at interfaces requires using methods that go beyond semilocal density functional theory to accurately describe the interactions between solvent, reactants and substrates. Methods based on many-body perturbation theory, such as the random phase approximation (RPA), have previously been limited due to their computational complexity. However, this is now a surmountable barrier due to the advances in computational power available, in particular through modern GPU-based supercomputers. In this work, we describe the implementation of RPA calculations within BerkeleyGW and show its favorable computational performance on large complex systems relevant for catalysis and electrochemistry applications. Our implementation builds off of the static subspace approximation which, by employing a compressed representation of the frequency dependent polarizability, enables the evaluation of the RPA correlation energy with significant acceleration and systematically controllable accuracy. We find that the computational cost of calculating the RPA correlation energy scales only linearly with system size for systems containing up to 50 thousand bands, and is expected to scale quadratically thereafter. We also show excellent strong scaling results across several supercomputers, demonstrating the performance and portability of this implementation.

algorithmic development↗
DOE OSTI2025

Thermodynamic Stability and Site‐Specific Distribution of Graphitic and Pyridinic Nitrogen in Graphene Moiré on Ru(0001)

Abstract Graphene‐like materials are of interest for large‐scale hydrogen storage applications due to their lightweight, durable, and scalable properties. Nitrogen‐doping minimizes kinetic limitations in diffusion and recombination on surfaces, however, the role of graphitic nitrogen (GN) and pyridinic nitrogen (PN) is not well understood. Nitrogen‐doped graphene is synthesized on Ru(0001) using chemical vapor deposition (CVD) of pyridine and ion irradiation. Scanning tunneling microscopy (STM), x‐ray photoelectron spectroscopy (XPS), and density functional theory (DFT) are used to identify the structure, location, and thermodynamic stability of nitrogen species within the graphene moiré. CVD of pyridine results in a low nitrogen concentration (<0.1at%), while the post‐growth nitrogen ion irradiation allows us to increase the concentration further. The concentration of GN and PN is controlled by varying the ion dose and annealing temperature. Comparison of measured and simulated STM images of GN and PN yield an excellent agreement, allowing us to confidently establish that GN is preferentially located near the center of the Atop region, while PN is located in the valley region of the graphene moiré. This report explicitly confirms the site assignments and provides a foundation for the site synthesis and analysis of structural and electronic properties that drive the reactivity of N‐doped graphene.

Gedara, Buddhika S. A. [Physical and Computational↗
NASA NTRS1995

Analysis of Aerothermal Heating Over Reusable Single-Stage to Orbit Vehicles

Using the General Aerodynamic Simulation Program (GASP) with finite-rate surface catalysis and radiative equilibrium surface temperature conditions, a computational study has been conducted to predict three-dimensional aero-convective heating over the proposed reusable single-stage to orbit vehicles. The solutions for three basic classes of usable single-stage-to-orbit vehicles, including wing body, lifting body and vertical lander, are presented and compared to provide information for vehicle design. The free stream conditions used in this study occur at an altitude of 57.5 km, ambient density of 4.129e-4 kg/cubic m, and velocity of 4,849 m/s. The convective heating over the vehicle surface in both laminar and turbulent flow regimes were computed. The Baldwin-Lomax algebraic model with a modified van Driest damping factor was applied in the turbulent flow calculations. A bifurcation model which guarantees the conservation of mass diffusive fluxes was adopted to account for the multi-species mass diffusion. The effects of turbulence and surface catalysis on the surface heating are discussed in detail The solutions indicate the degree of oxygen dissociation over the leeward side is sensitive to the shape of vehicle body. Predicted convective heating over the entire surface of lifting body is strongly affected by the presence of surface catalysis. For wing body and vertical lander, the surface catalysis effect does not appear to be important in the far downstream region, and turbulence is the dominant effect on the heating prediction in this region.

Chen, Y. K.↗
NASA NTRS2004

Montmorillonite, oligonucleotides, RNA and origin of life

Na-montmorillonite prepared from Volclay by the titration method facilitates the self-condensation of ImpA, the 5'-phosphorimidazolide derivative of adenosine. As was shown by AE-HPLC analysis and selective enzymatic hydrolysis of products, oligo(A)s formed in this reaction are 10 monomer units long and contain 67% 3',5'-phosphodiester bonds (Ferris and Ertem, 1992a). Under the same reaction conditions, 5'-phosphorimidazolide derivatives of cytidine, uridine and guanosine also undergo self-condensation producing oligomers containing up to 12-14 monomer units for oligo(C)s to 6 monomer units for oligo(G)s. In oligo(C)s and oligo(U)s, 75-80% of the monomers are linked by 2',5'-phosphodiester bonds. Hexamer and higher oligomers isolated from synthetic oligo(C)s formed by montmorillonite catalysis, which contain both 3',5'- and 2',5'-linkages, serve as catalysts for the non-enzymatic template directed synthesis of oligo(G)s from activated monomer 2-MeImpG, guanosine 5'-phospho-2-methylimidazolide (Ertem and Ferris, 1996). Pentamer and higher oligomers containing exclusively 2',5'-linkages, which were isolated from the synthetic oligo(C)s, also serve as templates and produce oligo(G)s with both 2',5'- and 3',5'-phosphodiester bonds. Kinetic studies on montmorillonite catalyzed elongation rates of oligomers using the computer program SIMFIT demonstrated that the rate constants for the formation of oligo(A)s increased in the order of 2-mer < 3-mer < 4-mer ... < 7-mer (Kawamura and Ferris, 1994). A decameric primer, dA(pdA)8pA bound to montmorillonite was elongated to contain up to 50 monomer units by daily addition of activated monomer ImpA to the reaction mixture (Ferris, Hill and Orgel, 1996). Analysis of dimer fractions formed in the montmorillonite catalyzed reaction of binary and quaternary mixtures of ImpA, ImpC, 2-MeImpG and ImpU suggested that only a limited number of oligomers could have formed on the primitive Earth rather than equal amounts of all possible isomers (Ertem and Ferris, 2000). Formation of phosphodiester bonds between mononucleotides by montmorillonite catalysis is a fascinating discovery, and a significant step forward in efforts to find out how the first RNA-like oligomers might have formed in the course of chemical evolution. However, as has been pointed out in several publications, these systems should be regarded as models rather than a literal representation of prebiotic chemistry (Orgel, 1998; Joyce and Orgel, 1999; Schwartz, 1999).

RNA/chemical synthesis↗