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 181 records · Page 10

Simplifying computational workflows with the Multiscale Atomic Zeolite Simulation Environment (MAZE)

Zeolites, an important class of 3-dimensional nanoporous materials, have been widely explored for a variety of applications including gas storage, separations, and catalysis. As the properties of these aluminosilicate materials depend on a number of factors (e.g., framework topology, Si/Al ratio, extra-framework cations etc.), detailed experiments (e.g., catalytic properties, adsorption capacities etc.) are often limited to only a handful of materials. Computational methods have played an important role in (1) providing molecular level insights to rationalize experimental observations, and (2) screening large libraries of zeolites to identify promising candidates for experimental synthesis and validation. Different levels of theory and computational chemistry codes are necessary to describe the range of relevant phenomena such as adsorption (e.g., grand canonical Monte Carlo), diffusion (e.g., molecular dynamics), and chemical reactions (e.g., density functional theory). Manipulation of atomic structures, handling of input files, and developing robust workflows becomes quite cumbersome. To mitigate these challenges, we describe the development of the Multiscale Atomic Zeolite Simulation Environment (MAZE) – a Python package that simplifies zeolite-specific calculation workflows by providing a user-friendly interface for systematically manipulating zeolite structures

97 MATHEMATICS AND COMPUTING↗

Computational procedure for three-dimensional boundary layers on aircraft and aerospace vehicles

An implicit numerical procedure is developed to solve the full three-dimensional laminar, boundary layer equations over realistic aircraft and spacecraft configurations. This procedure utilizes splined functions in the direction normal to the surface and finite differences in the cross flow direction. Calculation of off-pitch plane attachment and detachment lines can be calculated. The solution procedure permits the incorporation of a wide variety of turbulent models and boundary conditions including; entropy layer procedures, specified surface temperatures or heat flux distributions, surface catalysis, mass addition, and the coupling to both three-dimensional inviscid codes and surface ablation codes. In addition, reasonable computing times are obtained.

Kendall, R. M.↗

Simulating Catalysis with Realistic Pellet Geometries Using Mesoflow: A Case Study of Catalytic Propane Dehydrogenation

We present a case study of catalytic propane dehydrogenation with our open-source multiphysics solver, Mesoflow. The solver was developed to simulate reactive flow coupled to heterogeneous catalytic reactions and deactivation in the context of complex, mesoscale geometry. The method leverages cartesian block-structured adaptive mesh refinement to capture realistic catalyst microstructural features acquired directly from X-ray computed tomography data. A kinetic model for propane dehydrogenation and catalyst deactivation was developed based on temporal analysis of products (TAP) reactor experiments. The TAP reactor experiments allow for precise characterization of intrinsic kinetic reaction steps which are implemented into Mesoflow simulations to model the spatial and temporal evolution of reactants, products, and catalyst active sites. The short-term and long-term deactivation behavior is studied by using XCT data collected from fresh and aged catalyst pellets, which exhibit different microstructural features. This study employs time-splitting algorithms to connect disparate reaction and flow timescales, enabling the simulations to achieve realistic deactivation timescales on the order of minutes while the flow time-scales for small particles (100 microns) are several milliseconds. We also introduce a flexible automated python script that writes the necessary files to construct a Mesoflow simulation from user-created chemical mechanisms. We will also introduce a few new features that are added to Mesoflow such as higher order schemes, implicit chemistry integrators and the ability to run on AMD and NVIDIA graphics-processing-units.

AMReX↗

Numerical procedures for three-dimensional computational surface thermochemistry

Models and equations for surface thermochemistry and near-surface thermophysics of aerodynamically-heated thermal protection materials are reviewed, with particular emphasis on computational boundary conditions for surface mass and energy transfer. The surface energy and mass balances, coupled with an appropriate ablation or surface catalysis model, provide complete thermochemical boundary conditions for a true multidisciplinary solution of the fully coupled fluid-dynamics/solid mechanics problem. Practical approximate solutions can be obtained by using a detailed model with full thermophysics for either the solid or fluid phase amd a semianalytic method for the other half of the problem. A significant increase in the state-of-the-art in aerothermal computational fluid dynamics is possible by uniting CFD methodology with surface thermochemistry boundary conditions and the heat-balance-integral method.

Milos, Frank S.↗

Review of numerical procedures for computational surface thermochemistry

Models and equations for surface thermochemistry and near-surface thermophysics of aerodynamically heated thermal protection materials are reviewed, with particular emphasis on computational boundary conditions for surface mass and energy transfer. The surface energy and mass balances, coupled with an appropriate ablation or surface catalysis model, provide complete thermochemical boundary conditions for a true multidisciplinary solution of the fully coupled fluid-dynamics/solid mechanics problem. Practical approximate solutions can be obtained by using a detailed model with full thermophysics for either the solid or fluid phase and a semianalytic method for the other half of the problem. A significant increase in the state-of-the-art in aerothermal computational fluid dynamics is possible by uniting computational fluid dynamic (CFD) methodology with surface thermochemistry boundary conditions and the heat-balance-integral method.

SOLID MECHANICS↗

The 2020 plasma catalysis roadmap

Abstract Plasma catalysis is gaining increasing interest for various gas conversion applications, such as CO 2 conversion into value-added chemicals and fuels, CH 4 activation into hydrogen, higher hydrocarbons or oxygenates, and NH 3 synthesis. Other applications are already more established, such as for air pollution control, e.g. volatile organic compound remediation, particulate matter and NO x removal. In addition, plasma is also very promising for catalyst synthesis and treatment. Plasma catalysis clearly has benefits over ‘conventional’ catalysis, as outlined in the Introduction. However, a better insight into the underlying physical and chemical processes is crucial. This can be obtained by experiments applying diagnostics, studying both the chemical processes at the catalyst surface and the physicochemical mechanisms of plasma-catalyst interactions, as well as by computer modeling. The key challenge is to design cost-effective, highly active and stable catalysts tailored to the plasma environment. Therefore, insight from thermal catalysis as well as electro- and photocatalysis is crucial. All these aspects are covered in this Roadmap paper, written by specialists in their field, presenting the state-of-the-art, the current and future challenges, as well as the advances in science and technology needed to meet these challenges.

Bogaerts, Annemie (ORCID:0000000198756460)↗

Experimental and Computational Evaluation of Nicotinamide Cofactor Biomimetics

Oxidoreductase enzymes are widely used biocatalysts due to their high enantioselectivity and broad substrate compatibility in useful transformations. Many oxidoreductases require nicotinamide cofactors (i.e., NAD(P)H). To replace this costly natural cofactor, synthetic nicotinamide cofactor biomimetics (NCBs) offer different shapes, binding affinities, and reducing potentials that exceed the capabilities of wild-type NAD(P)H. However, the ill-defined structure–activity relationships (SARs) of various NCBs slow rationally guided innovation, such as customized reducing potentials. Here, we dissect two essential elements of NCB design, holding the nicotinamide invariant. First, the linker length between the nicotinamide and an unconjugated aromatic ring uncovered unexpected benefits to redox activity for two or three carbon linkers. Second, substitution on this unconjugated aryl group (Ring 2) might not be expected to affect activity. However, SAR trends demonstrate substantial benefits to reductive potential conferred by electron-donating functionalities on Ring 2. Furthermore, catalysis by two enzymes demonstrates enzyme-dependent tolerance or sensitivity to the NCB structures. Density functional theory (DFT) and computational modeling provide a theoretical framework to understand and build upon these observations. Ring 2 reaches up to the nicotinamide to stabilize its positive charge after oxidation through π–π stacking and charge transfer. Thus, the systematic examination of NCB’s stability, electrochemical redox potentials, and kinetics uncovers trends for the improved design of NCBs.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Practical Considerations for Understanding Surface Reaction Mechanisms Involved in Heterogeneous Catalysis

Acquiring useful knowledge about the active site(s) of a catalyst, nature of reactant–catalyst interactions, nature of reactive intermediates, rate-determining step, reaction rate orders that affect various process parameters, and reaction mechanism as a whole is exceedingly challenging. This is especially true in the case of heterogeneous catalysts due to the complexity of the nature of surface active sites and their nonstatic behavior. Here, we present our perspective on differentiating between various surface reaction mechanisms in light of pioneering studies by leaders in the field, with the aim of clarifying some of the confusion associated with these complex mechanisms, especially the Eley–Rideal mechanism. Using bibliometric analysis, we identify and discuss the following four reactions that most commonly invoke the Eley– Rideal mechanism: H 2 activation, CO oxidation, esterification of alcohols by acids, and selective catalytic reduction (SCR) of NO x with NH 3 . Our analysis of studies utilizing well-suited experimental and computational methodologies for differentiating surface reaction mechanisms suggests that the above-mentioned four reactions do not occur via the Eley–Rideal mechanism. Instead, each reaction occurs via the Langmuir–Hinshelwood mechanism with nonidealities present. Lastly, we highlight practical considerations regarding select experimental (characterization methods and differential kinetics) and computational modeling that we believe can provide useful insights to accurately discern between the various possible reaction mechanisms in heterogeneous catalysis.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Incorporating Coverage-Dependent Reaction Barriers into First-Principles-Based Microkinetic Models: Approaches and Challenges

Mean-field microkinetic models (MKMs) are appealing for their relatively facile construction, computational tractability, and high-throughput catalyst screening capabilities. As such, they will continue to be a valuable tool for materials design in heterogeneous catalysis even as the field aims to describe more complex systems. Numerous prior reports have provided the groundwork for constructing first-principles-based MKMs, including the analysis of strategies for incorporating lateral interactions into thermodynamic parameters (e.g., adsorption energies). Yet, there remains a need for concerted dialogue on methods for calculating and incorporating coverage-dependent kinetic parameters into MKMs. In this Perspective, we assess strategies for doing so, including the corresponding key physical implications and computational challenges. Here, we emphasize that decoupling thermodynamic and kinetic parameters within MKMs can violate thermodynamic consistency and risk unphysical solutions. For some reactions and catalyst materials, scaling relationships can predict coverage-dependent activation energies, but there are several exceptions evident in the literature, indicating that this approach is not universally applicable and that the field could benefit from research aimed at elucidating the limitations. Conducting high-coverage transition state searches is a rigorous but computationally costly strategy, and the effects of various methods for mitigating this cost on resulting energetics have yet to be broadly explored and validated. The goal of this Perspective is to generate discussion on and inspire focused research into the physical relevance of approaches for describing coverage-dependent reaction barriers in MKMs, including the development of computationally tractable methodologies, to advance the applicability of MKMs across diverse reaction chemistries and conditions.

36 MATERIALS SCIENCE↗

Growth in heterogeneous, evolving macromolecular networks: toward functional, biomimetic material

This report summarizes research carried out by the Balazs and Matyjaszewski Labs at Pitt and CMU, respectively, during 2025-26, supported by the DOE grant. It has produced the following 3 research articles and 2 review publications that acknowledged grant ER45998: 1. Computational Modeling of Hyperbranched Polymers 2. Synthesis of Structurally Tailored Networks 3. Sustainable and Oxygen-Tolerant Catalysis 4. A review paper on Current Status and Outlook for ATRP 5. A review paper on Future Directions for Atom Transfer Radical Polymerization

Balazs, Anna (ORCID:0000000255552692)↗

Combining Reactive Quantum-Mechanical Molecular-Dynamics Simulations with Mutagenesis, Crystallography, and Enzyme Kinetics to Reveal Plausible Steps of Isocyanide Hydratase Catalysis

A complete understanding of enzyme mechanisms requires atomistic details of chemical reactions. Quantum-based molecular dynamics simulations (QMD) are a potential source of this information, but trade-offs between accuracy and computational cost have limited their use. We previously developed extended Lagrangian Born–Oppenheimer molecular dynamics (XL-BOMD) methods that leverage a negligible compromise in accuracy to substantially decrease the cost of QMD simulations. Here, we develop a reactive QMD approach using the latest XL-BOMD formulation, which enables efficient simulations of highly reactive systems, and use it to investigate mechanisms of intermediate formation in isocyanide hydratase (ICH) catalysis. In QMD simulations, molecular analogs of ICH active site residues reacted with para-nitrophenyl isocyanide, forming a thioimidate. Analysis of simulated atomic configurational and charge dynamics revealed a pathway where protonation of the isocyanide carbon occurs prior to thioimidate formation and suggested a possible role of Asp17 as a proton donor in the early phase of ICH catalysis. To test whether the pathway seen using the reactive QMD approach might be relevant to ICH catalysis, we performed X-ray crystallography and pre-steady-state enzyme kinetics studies of wild-type and D17N mutant ICH. Both the structure and kinetics are sensitive to the D17N mutation in a manner that is consistent with the order of the reaction steps seen in the simulations. Mobile protons play essential roles in many enzymes, yet they are difficult to observe experimentally, making the ordering of proton-dependent steps ambiguous in many enzyme mechanisms. The ability to directly simulate model reactions for the design of experiments that provide information about enzyme mechanisms involving mobile protons demonstrates the significance of our reactive QMD approach and motivates further biological applications.

36 MATERIALS SCIENCE↗

Simulation and Optimization of Volatile Fatty Acid Upgrading Strategies for Sustainable Transportation Fuel Production

The multicomponent nature of and variability of biomass makes chemical upgrading via a single process stream to a single end use infeasible for most feedstocks. A more promising approach is to identify upgrading strategies (encompassing bioprocessing, catalysis, and separations) which valorize varied biomass fractions to distinct products to which each is best suited. A robust procedure to carry out this goal should incorporate comparison of various upgrading procedures as well as suitability of products to a slate of identified end uses. We applied this strategy to one biomass-derived feedstock, volatile fatty acids (VFAs) derived from wet waste arrested anaerobic digestion, by developing a computer program, VFA Upgrading to Liquid Transportation fUels Refinery Estimation (VULTURE) which evaluates VFA catalytic upgrading to liquid transportation fuels. VULTURE considers multiple separations, catalysis (ketonization, hydrogenation), and fuel application options, generating hundreds of candidate scenarios for a given VFA stream, then selects several promising strategies that optimize bio-content of products with properties best suited for target fuel types. We find that VFAs are upgraded most efficiently when separate light alcohol (C3-6) and heavy hydrocarbon or alcohol (C7-13) fractions are targeted to create gasoline and heavy-duty (diesel or jet) fuels or fuel blends. Surrogate property testing of VULTURE-recommended fuels reveals that most predictive models employed are robust, while rigorous process simulation shows that the simple unit operation assumptions used in VULTURE are largely accurate, especially for heavy-duty fuel synthesis. Techno-economic and life-cycle analyses show that VFA-derived fuels are profitable and have dramatically (=57%) lower carbon intensities than fossil analogs.

BIOMASS FUELS↗

Computing inelastic neutron scattering spectra from molecular dynamics trajectories

Abstract Inelastic neutron scattering (INS) provides a weighted density of phonon modes. Currently, INS spectra can only be interpreted for perfectly crystalline materials because of high computational cost for electronic simulations. INS has the potential to provide detailed morphological information if sufficiently large volumes and appropriate structural variety are simulated. Here, we propose a method that allows direct comparison between INS data with molecular dynamics simulations, a simulation method that is frequently used to simulate semicrystalline/amorphous materials. We illustrate the technique by analyzing spectra of a well-studied conjugated polymer, poly(3-hexylthiophene-2,5-diyl) (P3HT) and conclude that our technique provides improved volume and structural variety, but that the classical force field requires improvement before the morphology can be accurately interpreted.

36 MATERIALS SCIENCE↗

Spontaneous dynamical disordering of borophenes in MgB2 and related metal borides

Abstract Layered boron compounds have attracted significant interest in applications from energy storage to electronic materials to device applications, owing in part to a diversity of surface properties tied to specific arrangements of boron atoms. Here we report the energy landscape for surface atomic configurations of MgB 2 by combining first-principles calculations, global optimization, material synthesis and characterization. We demonstrate that contrary to previous assumptions, multiple disordered reconstructions are thermodynamically preferred and kinetically accessible within exposed B surfaces in MgB 2 and other layered metal diborides at low boron chemical potentials. Such a dynamic environment and intrinsic disordering of the B surface atoms present new opportunities to realize a diverse set of 2D boron structures. We validated the predicted surface disorder by characterizing exfoliated boron-terminated MgB 2 nanosheets. We further discuss application-relevant implications, with a particular view towards understanding the impact of boron surface heterogeneity on hydrogen storage performance.

36 MATERIALS SCIENCE↗

Model-based iterative reconstruction with adaptive regularization for artifact reduction in electron tomography

Obtaining high-quality 3D reconstructions from electron tomography of crystalline particles embedded in lighter support elements is crucial for various material systems such as catalysts for fuel cell applications. However, significant challenges arise due to the limited tilt range, sparse and low signal-to-noise ratio of the measurements. In addition, small metal particles can cause strong streaking and shading artifacts in the 3D reconstructions when using conventional reconstruction algorithms due to the presence of Bragg diffraction and the large scattering cross-section difference between the materials of the particles and the background support regions. These artifacts lead to errors in the downstream characterization affecting extraction of critical features such as the size of the metal particles, their distribution and the volume of the lighter support regions. In this paper, we present a two-stage algorithm based on metal artifact reduction, utilizing model-based iterative reconstruction methods with adaptive adjustment of regularization parameters. Our approach yields high-quality 3D reconstructions compared to traditional algorithms, accurately capturing both the metal particles as well as the background support. We demonstrate the effectiveness of our algorithm through simulated and experimental bright-field electron tomography data, showing significant improvements in reconstruction quality compared to traditional methods.

97 MATHEMATICS AND COMPUTING↗

Hydrocarbon, Oxidation, Dehydrogenation and Coupling Over Model Metal Oxide Surfaces

Final report for a 24.5 year single investigator project aimed at understanding structure/function relationships in adsorption and reaction on metal oxide surfaces for understanding heterogeneous catalysis. An experimental surface science approach was taken using single crystal surfaces as model catalysts and supplemented by density functional theory (DFT) calculations. Materials studied experimentally and computationally include the base metal oxide SnO 2 and transition metal oxides α-Cr 2 O 3 , α-Fe 2 O 3 and MnO.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Global Optimization of Chemical Cluster Structures: Methods, Applications, and Challenges

Chemical clusters are relevant to many applications in catalysis, separations, materials, and energy sciences. Experimentally, the structure of clusters is difficult to determine, but it is very important in understanding their chemistry and properties. Computational methods can be used to examine cluster structure, however finding the most stable structure is not simple, particularly as the cluster size increases. Global optimization techniques have long been used to tackle the problem of the most stable structure, but such approaches would have to look for a global minimum, while sampling local minima over the whole potential energy surface as well. In this review, the state-of-the-art theory of global optimization theory is summarized. First, the definition, significance, relation to experiments, and a brief history of global optimization is presented. We then discuss, in more detail, three versatile global optimization methods: the basin hopping, the artificial bee colony algorithm, and the genetic algorithm. We close with some representative application examples of global optimization of clusters since 2016 and the challenges, open questions and opportunities in this field.

Global optimization, Chemical clusters, Artificial↗

Code interoperability extends the scope of quantum simulations

The functionality of many materials is critically dependent on the integration of dissimilar components and on the interfaces that arise between them. The description of such heterogeneous components requires the development and deployment of first principles methods, coupled to appropriate dynamical descriptions of matter and advanced sampling techniques, in order to capture all the relevant length and time scales of importance to the materials’ performance. It is thus essential to build simple, streamlined computational schemes for the prediction and design of multiple properties of broad classes of materials, by developing interoperable codes which can be efficiently coupled to each other to perform complex tasks. We discuss the use of interoperable codes to simulate the structural and spectroscopic characterization of materials, including chemical reactions for catalysis, the description of defects for quantum information science, and heat and charge transport.

36 MATERIALS SCIENCE↗