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 271 records · Page 15
All sourcesNASA NTRSDOE OSTI
NASA NTRS1998

Computational Analysis of Shock Layer Emission Measurements in an Arc-Jet Facility

This paper reports computational analysis of radiation emission experiments in a high enthalpy arc-jet wind tunnel at NASA Ames Research Center. Recently, as part of ongoing arc-jet characterization work, spectroscopic radiation emission experiments have been conducted at the 20 MW NASA Ames arc-jet facility. The emission measurements were obtained from the arc-jet freestream and from a shock layer formed in front of flatfaced models. Analysis of these data is expected to provide valuable information about the thermodynamic state of the gas in the arc-jet freestream and in the shock layer as well as thermochemical equilibration processes behind the shock in arc-jet flows. Knowledge of the thermodynamic state of the gas in arc-jet test flows and especially within the shock layer is essential to interpret the heat transfer measurements such as in surface catalysis experiments. The present work is a continuation of previous work and focuses on analysis of the emission data obtained at relatively low-pressure conditions for which the arc-jet shock layer is expected to be in thermal and chemical nonequilibrium. Building blocks of the present computational analysis are: (1) simulation of nonequilibrium expanding flow in the converging-diverging conical nozzle and supersonic jet; (2) simulation of nonequilibrium shock layer formed in front of the flat-faced cylinder model; and (3) prediction of line-of-sight radiation from the computed flowfield. For computations of the nonequilibrium flow in the conical nozzle and shock layer, multi-temperature nonequilibrium codes with the axisymmetric formulation are used. For computations of line-of-sight radiation. a nonequilibrium radiation code (NEQAIR) is used to predict emission spectra from the computed flowfield. Computed line-of-sight averaged flow properties such as vibrational and rotational temperatures, species number densities within the shock layer will be compared with those deduced from the experimental spectra. Detailed comparisons of computational and experimental spectra will also be presented.

Gokcen, Tahir↗
DOE OSTI2020

Active Learning A Neural Network Model For Gold Clusters & Bulk From Sparse First Principles Training Data

Small metal clusters are of fundamental scientific interest and of tremendous significance in catalysis. These nanoscale clusters display diverse geometries and structural motifs depending on the cluster size; a knowledge of this size-dependent structural motifs and their dynamical evolution has been of longstanding interest. Given the high computational cost of first-principles calculations, molecular modeling and atomistic simulations such as molecular dynamics (MD) has proven to be an important complementary tool to aid this understanding. Classical MD typically employ predefined functional forms which limits their ability to capture such complex size-dependent structural and dynamical transformation. Neural Network (NN) based potentials represent flexible alternatives and in principle, well-trained NN potentials can provide high level of flexibility, transferability and accuracy on-par with the reference model used for training. A major challenge, however, is that NN models are interpolative and requires large quantities (similar to 10 4 or greater) of training data to ensure that the model adequately samples the energy landscape both near and far-from-equilibrium. A highly desirable goal is minimize the number of training data, especially if the underlying reference model is first-principles based and hence expensive. In this work, we introduce an active learning (AL) scheme that trains a NN model on-the-fly with minimal amount of first-principles based training data. Our AL workflow is initiated with a sparse training dataset (similar to 1 to 5 data points) and is updated on-the-fly via a Nested Ensemble Monte Carlo scheme that iteratively queries the energy landscape in regions of failure and updates the training pool to improve the network performance. Using a representative system of gold clusters, we demonstrate that our AL workflow can train a NN with similar to 500 total reference calculations. Using an extensive DFT test set of similar to 1100 configurations, we show that our AL-NN is able to accurately predict both the DFT energies and the forces for clusters of a myriad of different sizes. Our NN predictions are within 30 meV/atom and 40 meV/angstrom of the reference DFT calculations. Moreover, our AL-NN model also adequately captures the various size-dependent structural and dynamical properties of gold clusters in excellent agreement with DFT calculations and available experiments. We finally show that our AL-NN model also captures bulk properties reasonably well, even though they were not included in the training data.

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

Inverse methods for design of soft materials

Functional soft materials, comprising colloidal and molecular building blocks that self-organize into complex structures as a result of their tunable interactions, enable a wide array of technological applications. Inverse methods provide a systematic means for navigating their inherently high-dimensional design spaces to create materials with targeted properties. Furthermore, while multiple physically motivated inverse strategies have been successfully implemented in silico, their translation to guiding experimental materials discovery has thus far been limited to a handful of proof-of-concept studies. In this perspective, we discuss recent advances in inverse methods for design of soft materials that address two challenges: (1) methodological limitations that prevent such approaches from satisfying design constraints and (2) computational challenges that limit the size and complexity of systems that can be addressed. Strategies that leverage machine learning have proven particularly effective, including methods to discover order parameters that characterize complex structural motifs and schemes to efficiently compute macroscopic properties from the underlying structure. We also highlight promising opportunities to improve the experimental realizability of materials designed computationally, including discovery of materials with functionality at multiple thermodynamic states, design of externally directed assembly protocols that are simple to implement in experiments, and strategies to improve the accuracy and computational efficiency of experimentally relevant models.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Steric Mapping, Ligand Dynamics, and Cycloisomerization Catalysis with Redox Robust Mn I/0/-I Dicarbenes

Manganese is perhaps the most electronically versatile element, yet the redox properties, reactivity, and catalytic applications of low-valent Mn 0 /Mn −I complexes remain underexplored due to the propensity for Mn 0 to dimerize, quenching highenergy metalloradicals. We report a series of redox-active monometallic Mn I , Mn 0 and Mn −I complexes containing a BH 2 - bridged dicarbene, characterized using a suite of experimental and cutting-edge computational (DFT) methods. Slow electron transfer kinetics at Mn I/0 are observed, with computations and electrochemical simulations in excellent agreement with experimental values. Despite the lack of steric bulk at the BH 2 -bridged Mn 0 , the t Bu groups at the dicarbene provide adequate steric protection to prevent dimerization, with percent buried volume (% V bur ) serving as a valuable steric ranking tool. We also show that a %V bur > 83% prevents dimerization for a diverse array of Mn 0 complexes from the literature. Ligand sterics of BPh 2 -and BH 2 -bridged complexes dictate reaction outcomes when Mn I and Mn −I are exposed to nucleophiles and electrophiles, respectively, while Mn0 facilitates the radical cycloisomerization catalysis of 6-iodo-1- hexene at room temperature. Furthermore, this work underscores the importance of ligand sterics in rationalizing reactivity patterns at Mn and provides valuable insights for designing chelating ligands that can selectively leverage Mn I/0/‑I states in redox-mediated catalytic reactions.

Ligands↗
DOE OSTI2022

Accessing and Photo-Accelerating Low-Overpotential Pathways for CO 2 Reduction: A Bis-Carbene Ruthenium Terpyridine Catalyst

A ruthenium catalyst bearing a bidentate bis(carbene) ligand is prepared and studied as a catalyst for CO 2 electroreduction. The catalyst [Ru(tpy)(bis-mim)(MeCN)][PF 6 ] 2 (tpy) is 2,2′,:6′,2″-terpyridine; bis-mim is (methylenebis(N-methylimidazol-2-ylidene)) mediates reduction of CO 2 into CO with a turnover frequency of 630 s –1 and Faradaic efficiency (FE) of 30% at an overpotential of 730 mV. The strongly donating bis(carbene) ligand also enables access to a pathway operating at a lower overpotential of ca. 310 mV. While low-overpotential catalysis is slow in the dark (TOF = 0.01 s –1 ), visible light illumination increases the rate 10-fold (TOF = 0.11 s –1 ). Here, a full mechanistic picture is developed using kinetic analysis from cyclic voltammetry, spectroelectrochemistry, and computational methods, with the bis-mim ligand facilitating rapid CO 2 activation at low overpotentials. Comparisons with other ruthenium catalysts yield insight into the ability to tune the rate of chemical steps (e.g., ligand dissociation and CO 2 nucleophilic attack) and the overpotential by tailoring the primary coordination sphere while retaining the “redox-active” tpy ligand.

CO2 reduction↗
NASA NTRS1984

Fundamental organometallic reactions: Applications on the CYBER 205

Two of the most challenging problems of Organometallic chemistry (loosely defined) are pollution control with the large space velocities needed and nitrogen fixation, a process so capably done by nature and so relatively poorly done by man (industry). For a computational chemist these problems are on the fringe of what is possible with conventional computers (large models needed and accurate energetics required). A summary of the algorithmic modification needed to address these problems on a vector processor such as the CYBER 205 and a sketch of findings to date on deNOx catalysis and nitrogen fixation are presented.

Rappe, A. K.↗
DOE OSTI2024

Actuating Liquid Crystals Rapidly and Reversibly by Using Chemical Catalysis

Abstract Microtubules and catalytic motor proteins underlie the microscale actuation of living materials, and they have been used in reconstituted systems to harness chemical energy to drive new states of organization of soft matter (e.g., liquid crystals (LCs)). Such materials, however, are fragile and challenging to translate to technological contexts. Rapid (sub‐second) and reversible changes in the orientations of LCs at room temperature using reactions between gaseous hydrogen and oxygen that are catalyzed by Pd/Au surfaces are reported. Surface chemical analysis and computational chemistry studies confirm that dissociative adsorption of H 2 on the Pd/Au films reduces preadsorbed O and generates 1 ML of adsorbed H, driving nitrile‐containing LCs from a perpendicular to a planar orientation. Subsequent exposure to O 2 leads to oxidation of the adsorbed H, reformation of adsorbed O on the Pd/Au surface, and a return of the LC to its initial orientation. The roles of surface composition and reaction kinetics in determining the LC dynamics are described along with a proof‐of‐concept demonstration of microactuation of beads. These results provide fresh ideas for utilizing chemical energy and catalysis to reversibly actuate functional LCs on the microscale.

Chemistry↗
DOE OSTI2026

Exploration of the Electronic and Catalytic Properties of [Co 5 MS 8 (PEt 3 ) 5 ] 1+ Nanoclusters: A Computational Study

Recent studies have demonstrated the relative stability of undercoordinated hexanuclear cobalt sulfide nanoclusters (NCs) with different charge states. Considering that these small metal NCs have atomically precise structures and high reactivity due to the open shell of the transition metals, and provide selectivity toward ligand loss, they are a vital model for catalysis. In this paper, the electronic structures of these NCs are investigated. These NCs are then used as the reference state to analyze the catalytic properties with respect to hydrogen evolution reaction (HER) and CO 2 reduction (CO 2 R). Further, to understand the effect of heteroatom incorporation, the geometry and reactivity of ten different metal dopants are analyzed. This work shows that the type of metal incorporation greatly affects the electronic structure and formation energies for ligand binding and catalysis. Particularly, the d-orbital occupancy in the cobalt atoms remains largely unchanged, while the heteroatom greatly influences the reactivity of the undercoordinated NCs. Most notably, this work highlights that transition metals in [Co 5 MS 8 (PEt 3 ) 5 ] 1+ NCs would competitively prefer electrochemical adsorption of H over COOH, while the main group metals prefer COOH adsorption.

catalysis↗
DOE OSTI2021

Sustainable Oxide Electrocatalyst for Hydrogen- and Oxygen-Evolution Reactions

In this work, an electrocatalyst for water splitting based on earth-abundant metals is reported. This perovskite-oxide catalyst, CaSrFe 0.75 Co 0.75 Mn 0.5 O 6-δ (CSFCM), is examined using both experimental and computational methods. It demonstrates a combination of properties, which include (a) very high activity for the oxygen-evolution reaction with an overpotential of eta = 0.19 V at 10 mA/cm 2 , (b) high stability over 1000 cycles of catalysis, (c) the ability to catalyze the hydrogen-evolution reaction effectively in both acidic and basic conditions, and (d) catalytic activity as a single-phase bulk material without the need for any additional processing, multicomponent composite preparation, or nanofabrication. Therefore, the catalytic activity of CSFCM is intrinsic, making it a good benchmark compound for future studies of electrocatalytic parameters. This work also highlights the impact of systematic structural design on electrocatalytic activity. Results from density functional theory calculations indicate that in addition to an optimal e g occupancy of similar to 1, an additional descriptor, i.e., maximizing the number of free e g carriers, correlates with the electrocatalytic activity.

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

Stability of heterogeneous single-atom catalysts: a scaling law mapping thermodynamics to kinetics

Heterogeneous single-atom catalysts (SACs) hold the promise of combining high catalytic performance with maximum utilization of often precious metals. We extend the current thermodynamic view of SAC stability in terms of the binding energy (E bind ) of single-metal atoms on a support to a kinetic (transport) one by considering the activation barrier for metal atom diffusion. A rapid computational screening approach allows predicting diffusion barriers for metal–support pairs based on E bind of a metal atom to the support and the cohesive energy of the bulk metal (E c ). Metal–support combinations relevant to contemporary catalysis are explored by density functional theory. Assisted by machine-learning methods, we find that the diffusion activation barrier correlates with (E bind ) 2 /E c in the physical descriptor space. This diffusion scaling-law provides a simple model for screening thermodynamics to kinetics of metal adatom on a support.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Degrees of rate control and AutoDiff-driven direct sensitivity analysis in heterogeneous catalysis

Despite the wide application and benefits of the degree of rate control (DRC) analysis, several details remain argued, particularly about the conservation of DRCs at transient (TR) and steady-state (SS) conditions, especially for complex reaction networks. This work argues that previous proofs about the conservation properties of DRCs have been incomplete, and we provide new mathematical proofs at TR and SS conditions. In addition, we use both analytical (automatic differentiation) and numerical (finite difference) approaches to compute DRCs for the case study of ethane hydrogenolysis (EH) over Pt(111). This work confirms that at both TR and SS conditions, the sum of all DRCs, i.e., sum of the degrees of kinetic (DKRC) and thermodynamic rate control (DTRC), is conserved at zero. At SS conditions, the sum of DKRC is conserved at 1 while the sum of DTRC is conserved at −1. In corroboration of previous works, we show that the DTRC for any adsorbate at SS is equal to the product of the species coverage and a constant. In contrast, at TR conditions, the individual sums of both DTRC and DKRC are not conserved and can be any real number, with potential implications for the novel field of dynamic catalysis. Finally, we show that the conventional finite difference (FD) approach, only useful at SS, is prone to inaccuracy and very sensitive to the value of the differential change applied. The optimal differential value also varies significantly with system and rate definition. Consequently, we describe and illustrate in this work the application of the automatic differentiation (AD) approach for the more accurate determination of DRCs at both TR and SS conditions.

Automatic differentiation↗
DOE OSTI2025

Prediction of hydration energies of adsorbates at Pt(111) and liquid water interfaces using machine learning

Aqueous phase heterogeneous catalysis is important to various industrial processes, including biomass conversion, Fischer–Tropsch synthesis, and electrocatalysis. Accurate calculation of solvation thermodynamic properties is essential for modeling the performance of catalysts for these processes. Explicit solvation methods employing multiscale modeling, e.g., involving density functional theory and molecular dynamics have emerged for this purpose. Although accurate, these methods are computationally intensive. This study introduces machine learning (ML) models to predict solvation thermodynamics for adsorbates on a Pt(111) surface, aiming to enhance computational efficiency without compromising accuracy. In particular, ML models are developed using a combination of molecular descriptors and fingerprints and trained on previously published water–adsorbate interaction energies, energies of solvation, and free energies of solvation of adsorbates bound to Pt(111). These models achieve root mean square error values of 0.09 eV for interaction energies, 0.04 eV for energies of solvation, and 0.06 eV for free energies of solvation, demonstrating accuracy within the standard error of multiscale modeling. Feature importance analysis reveals that hydrogen bonding, van der Waals interactions, and solvent density, together with the properties of the adsorbate, are critical factors influencing solvation thermodynamics. Furthermore, these findings suggest that ML models can provide rapid and reliable predictions of solvation properties. This approach not only reduces computational costs but also offers insights into the solvation characteristics of adsorbates at Pt(111)–water interfaces.

Adsorption↗
DOE OSTI2026

Metalloborophenes: Structural Diversity and Emerging Properties of Metal–Boron Two‐Dimensional Frameworks

Metalloborophenes, an emerging subclass of 2D materials, have attracted growing attention owing to their exceptional structural diversity and highly tunable electronic and magnetic properties. Constructed from vacancy‐rich borophene frameworks stabilized by electron donation from incorporated metal atoms, metalloborophenes merge the chemical versatility of boron with the functional richness of metal dopants. The first experimental realization of Cu–borophene nanoribbons in 2024 marked a pivotal advance, confirming long‐standing theoretical predictions and revitalizing interest in this new frontier of boron‐based 2D chemistry. Despite this progress, most studies to date remain conceptual and theoretical, constrained by challenges in scalable synthesis, dopant precision, and substrate control. Computational investigations have revealed a broad landscape of stable metalloborophene structures, exhibiting metallic, semiconducting, and magnetic behavior across diverse dopant families, including alkali, alkaline‐earth, transition, and lanthanide elements. These tunable characteristics open promising avenues for applications in spintronics, catalysis, and hydrogen storage. This review provides a comprehensive overview of metalloborophenes, emphasizing the interplay between structure, stability, and functionality, and outlining future directions toward bridging predictive modeling with experimental realization of this rapidly evolving class of 2D materials.

2D materials↗
DOE OSTI2021

Identification of Known and Novel Monomers for Poly(hydroxyurethanes) from Biobased Materials

Sourcing polymers from biobased materials is desirable for a transition to a more sustainable world. However, finding new monomers sourced from biobased molecules requires extensive experimental effort, due to the fact that most biomolecules lack the functional groups or properties needed to serve as a monomer in a particular polymeric system without chemical modification. In this work, we demonstrate a computational screening method to obtain monomers of a non-isocyanate polyurethane (NIPU), polyhydroxyurethane (PHU), from biomolecules. Our method generates candidate monomers from a list of 15 selected biobased starting molecules by performing chemical transformations that are standard to chemical catalysis and organic synthesis on these starting molecules and then screening the results for molecules with moieties that can be used as monomers for a PHU. Product molecules from two generations of reaction, a total number of 20,057, were analyzed automatically to identify molecules that have two or more epoxy or cyclic carbonate functional groups as candidates for further derivatization. These 21 candidate molecules were then subjected to another two generations of reaction, producing a set of 17,913 additional molecules. Further scrutiny of this set for desired functional groups characterizing PHUs and querying against the PubChem database was carried out, and 15 known molecules that are commercially available were identified, with the desirable features that they are predicted to be derived from biobased starting molecules in four or fewer steps. Moreover, molecule space generated from this process was primarily comprised of molecules that were not known to PubChem, indicating the framework also has the ability to generate novel candidate monomers that originate from biobased starting points.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Modelling drug adsorption in metal–organic frameworks: the role of solvent

Solvent plays a key role in biological functions, catalysis, and drug delivery. Metal–organic frameworks (MOFs) due to their tunable functionalities, porosities and surface areas have been recently used as drug delivery vehicles. To investigate the effect of solvent on drug adsorption in MOFs, we have performed integrated computational and experimental studies in selected biocompatible MOFs, specifically, UiO-AZB, HKUST-1 (or CuBTC) and NH 2 -MIL-53(Al). The adsorption of three drugs, namely, 5-fluorouracil (5-FU), ibuprofen (IBU), and hydroxyurea (HU) were performed in the presence and absence of the ethanol. Our computational predictions, at 1 atmospheric pressure, showed a reasonable agreement with experimental studies performed in the presence of ethanol. We find that in the presence of ethanol the drug molecules were adsorbed at the interface of solvent and MOFs. Moreover, the computationally calculated adsorption isotherms suggested that the drug adsorption was driven by electrostatic interactions at lower pressures (<10 –4 Pa). Our computational predictions in the absence of ethanol were higher compared to those in the presence of ethanol. The MOF–adsorbate interaction (U HA ) energy decreased with decrease in the size of a drug molecule in all three MOFs at all simulated pressures. At high pressure the interaction energy increases with increase in the MOFs pore size as the number of molecules adsorbed increases. Thus, our research shows the important role played by solvent in drug adsorption and suggests that it is critical to consider solvent while performing computational studies.

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

SPARC-X: Quantum simulations at extreme scale - reactive dynamics from first principles

We have developed the massively parallel electronic structure code SPARC-X: a computational framework for performing Kohn-Sham Density Functional Theory (DFT) calculations that can scale linearly with the number of atoms in the system, while being able to leverage petascale and emerging exascale parallel computers to study chemical phenomena at unprecedented length and time scales. SPARC-X exploits a recent breakthrough in electronic structure methodologies: systematically improvable, strictly local, orthonormal, discontinuous real-space bases that efficiently and systematically capture the local chemistry of the system. With further adaptation using new machine-learning techniques and the use of the massively parallel Spectral Quadrature (SQ) electronic structure method, the algorithmic complexity and prefactor associated with DFT calculations involving semilocal as well as hybrid functionals are dramatically reduced. Using petascale computational resources, SPARC-X enables quantum mechanical simulations at length and time scales previously accessible only by empirical approaches, e.g., 1,000,000 atoms for a few picoseconds using semilocal functionals or 1,000 atoms for a few picoseconds using hybrid functionals. Using exascale resources, the sizes and times targeted are two orders of magnitude larger. Such a capability has applications in a wide variety of chemical sciences, including reactive interfaces where large length- and/or long time-scales are needed and traditional force fields fail. This is particularly important in dynamic catalysis, where bond breaking and formation must be understood in detail. We developed, tested, and employed the SPARC-X framework to understand the photocatalytic properties of TiO 2 nanoparticles, revealing finite size effects that cannot be captured with standard model systems or functionals. This integrated development and application strategy ensures that SPARC-X remains a robust, efficient, and scalable software package for quantum simulations on current petascale and emerging exascale computing resources.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Remote stereocontrol with azaarenes via enzymatic hydrogen atom transfer

Strategies for achieving asymmetric catalysis with azaarenes have traditionally fallen short of accomplishing remote stereocontrol, which would greatly enhance accessibility to distinct azaarenes with remote chiral centres. The primary obstacle to achieving superior enantioselectivity for remote stereocontrol has been the inherent rigidity of the azaarene ring structure. Here we introduce an ene-reductase system capable of modulating the enantioselectivity of remote carbon-centred radicals on azaarenes through a mechanism of chiral hydrogen atom transfer. This photoenzymatic process effectively directs prochiral radical centres located more than six chemical bonds, or over 6 Å, from the nitrogen atom in azaarenes, thereby enabling the production of a broad array of azaarenes possessing a remote γ-stereocentre. Finally, results from our integrated computational and experimental investigations underscore that the hydrogen bonding and steric effects of key amino acid residues are important for achieving such high stereoselectivities.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗