Search NASA⌕ Search

SEARCH · Search NASA

Results for “kinetic network”

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

A kinetic model for simulating non-equilibrium mass transport in oxides applied to hematite growth under irradiation

We propose an approach to simulating the dynamic evolution and transport of charged point defects within and through the oxide scales that form during corrosion. The method follows the cluster dynamics formalism widely adopted for radiation damage in solids, which can apply in both quasi-static and far from equilibrium conditions, such as in radiation environments. By treating each charged defect as a cluster of an atomic defect and associated unit charges, the proposed model flexibly allows charge state transitions and shifts in the Fermi level by absorption and emission of charge carriers from point defects in a reaction network with rates constrained by local equilibrium. Applying this model to hematite predicts changes in self-diffusion and oxidation kinetics in irradiation environments, surprisingly showing reduced oxidation rates in many conditions, despite enhanced self-diffusion. In conclusion, the origin of this effect lies in a change in Fermi level induced by excess vacancies formed under irradiation, which in turn suppresses the transport of cation interstitials which facilitate hematite growth.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Integrative teaching of metabolic modeling and flux analysis with interactive python modules

Abstract The modeling of rates of biochemical reactions—fluxes—in metabolic networks is widely used for both basic biological research and biotechnological applications. A number of different modeling methods have been developed to estimate and predict fluxes, including kinetic and constraint‐based (Metabolic Flux Analysis and flux balance analysis) approaches. Although different resources exist for teaching these methods individually, to‐date no resources have been developed to teach these approaches in an integrative way that equips learners with an understanding of each modeling paradigm, how they relate to one another, and the information that can be gleaned from each. We have developed a series of modeling simulations in Python to teach kinetic modeling, metabolic control analysis, 13C‐metabolic flux analysis, and flux balance analysis. These simulations are presented in a series of interactive notebooks with guided lesson plans and associated lecture notes. Learners assimilate key principles using models of simple metabolic networks by running simulations, generating and using data, and making and validating predictions about the effects of modifying model parameters. We used these simulations as the hands‐on computer laboratory component of a four‐day metabolic modeling workshop and participant survey results showed improvements in learners' self‐assessed competence and confidence in understanding and applying metabolic modeling techniques after having attended the workshop. The resources provided can be incorporated in their entirety or individually into courses and workshops on bioengineering and metabolic modeling at the undergraduate, graduate, or postgraduate level.

Kaste, Joshua A. M.↗
DOE OSTI2025

Machine Learning-Assisted Recovery of Delicate Kinetic Information from Transient Reactor Experiments

Identifying active sites and their roles in chemical reaction steps remains a vital challenge in heterogeneous catalysis. Transient experiments offer a unique way to probe active sites and distinguish subtle kinetic features. Although physics-based analysis methods may be well-developed, they can be highly susceptible to experimental noise, and smoothing methods may erase or even distort important features; a smooth curve is not always the best curve. We demonstrate a new workflow for the direct interpretation of intrinsic kinetic information from exit flux curves measured in transient reactor experiments. This workflow contains three artificial neural networks (ANNs), including a noise reducer, a concentration predictor, and a rate predictor to analyze experimental data, followed by the virtual TAP (VTAP) physics-based reactor model and density functional theory (DFT) calculations of adsorption energies on specific sites. We use this workflow to analyze the data from experiments titrating Pt/Al 2 O 3 and Pt/SiO 2 catalysts with carbon monoxide (CO) in the temporal analysis of products (TAP) reactor. Our workflow separates the time-evolving chemical reaction and mass transfer information contained in the TAP pulse response. The existence of strong- and weak-binding sites on the Pt/Al 2 O 3 catalyst is observed in the catalyst titration experiment in the transient reactor. The structures of the strong- and weak-binding sites are then identified by using DFT calculations. We find that the Pt/SiO 2 catalyst has only strong-binding sites, which aligns with the inactive support effect of SiO 2 . We demonstrate how machine learning methods provide unique insights with high-resolution data analysis that cannot be achieved by using state-of-the-art physics-based methods.

Adsorption↗
DOE OSTI2023

Understanding hydrogen electrocatalysis by probing the hydrogen-bond network of water at the electrified Pt–solution interface

Rational construction of the electrode–solution interface where electrochemical processes occur is of paramount importance in electrochemistry. Efforts to gain better control and understanding of the interface have been hindered by lack of probing methods. Here, in this study, we show that the hydrogen evolution and oxidation reactions (HER/HOR) catalysed by platinum in base can be promoted by introduction of N-methylimidazoles at the platinum–water interface. In situ spectroscopic characterization together with simulations indicate that the N-methylimidazoles facilitate diffusion of hydroxides across the interface by holding the second layer of water close to platinum surfaces, thereby promoting the HER/HOR. We thus propose that the HER/HOR kinetics of platinum in acid and base is governed by diffusion of protons and hydroxides, respectively, through the hydrogen-bond network of interfacial water by the Grotthuss mechanism. Moreover, we demonstrate a 40% performance improvement of an anion exchange membrane electrolyser by adding 1,2-dimethylimidazole into the alkali fed into its platinum cathode.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Biopolymer-Templated Titania Film Formation for Nanostructured Coatings Revealed by Machine Learning-Supported Time-Resolved Analysis

This study presents a machine learning approach to derive the film formation of biopolymer-templated titania nanostructures during spray deposition, in combination with in situ grazing-incidence small-angle X-ray scattering (GISAXS). A neural network trained on synthetic GISAXS data directly predicts domain-size distributions from experimental two-dimensional scattering patterns, capturing the full kinetics of nanostructure evolution with high temporal resolution. The predictions reveal hierarchical size distributions and periodic growth features, consistent with layer-by-layer spray deposition and validated by complementary scanning electron microscopy (SEM) imaging. Quantitative comparison with conventional parametric GISAXS fits shows good qualitative agreement, with systematic differences explained by domain-shape assumptions and resolved by applying a geometric scaling factor. Simulated SEM-like surfaces derived from neural network outputs reproduce the porous, foam-like nanoscale morphology observed experimentally, reinforcing the method’s credibility. This integrated approach enables real-time, nondestructive, statistically averaged monitoring of bulk nanostructure development in functional coatings, offering a scalable methodology to accelerate the characterization and process control of sustainably manufactured nanostructured titania films for energy-related applications such as photocatalysis and photovoltaics.

Heger, JulianEliah↗
DOE OSTI2020

Reaction Mechanism of Vapor-Phase Formic Acid Decomposition over Platinum Catalysts: DFT, Reaction Kinetics Experiments, and Microkinetic Modeling

A combination of periodic density functional theory (DFT, PW91-GGA) calculations, reaction kinetics experiments, and mean-field microkinetic modeling is used to derive insights on the reaction mechanism and determine the nature of the active site under reaction conditions for the vapor-phase decomposition of formic acid (FA, HCOOH) over Pt/C catalysts. Microkinetic models formulated using DFT energetics derived on the clean Pt(100) and Pt(111) required large parameter adjustments to reproduce the experimentally measured apparent activation energies and reaction orders. Further, these models predicted high surface coverage of adsorbed carbon monoxide (CO*), inconsistent with the environment of the active site in the DFT calculations on the clean surfaces. Consequently, we reperformed DFT calculations for the entire reaction network on partially CO*-covered (4/9 monolayer, ML) Pt(111) and Pt(100). The resultant microkinetic models, with thermochemistry and kinetics explicitly dependent on CO* coverage, were able to reproduce the experimentally determined activation energies and reaction orders, in addition to being self-consistent in CO* coverage. Our results suggest that Pt(100) is likely poisoned by CO* under typical reaction conditions and does not contribute significantly to the experimentally observed reactivity. Instead, we find that Pt(111) better represents the active site for FA decomposition reaction on Pt/C catalysts. The optimized model on 4/9 ML CO*-covered Pt(111) suggests that the reaction occurs via the carboxyl (COOH*) intermediate and that the spectator CO*-assisted pathways play a significant role under reaction conditions. Lastly, this study underscores the importance of spectator species on the energetics and the mechanism of a catalytic reaction and their key role in developing a model that better addresses the nature of the active site under realistic catalytic reaction conditions.

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

Encapsulation of Nanoparticle Organic Hybrid Materials within Electrospun Hydrophobic Polymer/Ceramic Fibers for Enhanced CO 2 Capture

Liquid-like nanoparticle organic hybrid materials (NOHMs) consisting of a silica core with ionically grafted branched polyethyleneimine chains (referred to as NIPEI) are encapsulated within submicron-scale polyacrylonitrile (PAN)/polymer-derived-ceramic electrospun fibers. The addition of a room-temperature curable, liquid-phase organopolysilazane (OPSZ) ceramic precursor to the PAN/NOHM solution enhances the internal dispersion of NOHMs and forms a thin ceramic sheath layer on the fiber exterior, shielding the hydrophilic NIPEI to produce near-superhydrophobic non-woven fiber mats with contact angles exceeding 140°. 60:40 loadings of NOHMs to PAN/OPSZ can be reliably achieved with low OPSZ percentages required, and up to 4:1 NOHM:polymer loadings are possible before losing hydrophobicity. These fibers demonstrate up to ≈2 mmol CO 2 g -1 fiber capture capacities in a pure CO 2 atmosphere through the nonwoven fibrous networks and the permeability of the OPSZ shell. In conclusion, the hybrid fibers additionally show enhanced capture kinetics under pure CO 2 and 400 ppm CO 2 conditions, indicating their promising application as a direct air capture platform.

36 MATERIALS SCIENCE↗
DOE OSTI2023

The Identity and Chemistry of C 7 H 7 Radicals Observed during Soot Formation

Here we used aerosol mass spectrometry coupled with tunable synchrotron photoionization to measure radical and closed-shell species associated with particle formation in premixed flames and during pyrolysis of butane, ethylene, and methane. We analyzed photoionization (PI) spectra for the C 7 H 7 radical to identify the isomers present during particle formation. For the combustion and pyrolysis of all three fuels, the PI spectra can be fit reasonably well with contributions from four radical isomers: benzyl, tropyl, vinylcyclopentadienyl, and o-tolyl. Although there are significant experimental uncertainties in the isomeric speciation of C7H 7 , the results clearly demonstrate that the isomeric composition of C 7 H 7 strongly depends on the combustion or pyrolysis conditions and the fuel or precursors. Fits to the PI spectra using reference curves for these isomers suggest that all of these isomers may contribute to m/z 91 in butane and methane flames, but only benzyl and vinylcyclopentadienyl contribute to the C 7 H 7 isomer signal in the ethylene flame. Only tropyl and benzyl appear to play a role during pyrolytic particle formation from ethylene, and only tropyl, vinylcyclopentadienyl, and o-tolyl appear to participate during particle formation from butane pyrolysis. There also seems to be a contribution from an isomer with an ionization energy below 7.5 eV for the flames but not for the pyrolysis conditions. Kinetic models with updated and new reactions and rate coefficients for the C 7 H 7 reaction network predict benzyl, tropyl, vinylcyclopentadienyl, and o-tolyl to be the primary C 7 H 7 isomers and predict negligible contributions from other C 7 H 7 isomers. These updated models provide better agreement with the measurements than the original versions of the models but, nonetheless, underpredict the relative concentrations of tropyl, vinylcyclopentadienyl, and o-tolyl in both flames and pyrolysis and overpredict benzyl in pyrolysis. Our results suggest that there are additional important formation pathways for the vinylcyclopentadienyl, tropyl, and o-tolyl radicals and/or loss pathways for the benzyl radical that are currently unaccounted for in the present models.

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

Cleavage of C-O and C-C Bonds in Lignin-Derived Compounds to Produce Aromatics Using Molybdenum-Containing MFI Zeolites

Lignin, the most abundant source of renewable arenes, is a viable feedstock for the production of aromatic compounds. However, the prevalence of resilient C-C bonded oligomeric fragments in lignin-derived streams can compromise monomer yields during reductive catalytic fractionation (RCF). To address this issue, we developed a bifunctional molybdenum-containing MFI (Mo/H-MFI) zeolite catalyst capable of cleaving both C-O and C-C bonds in lignin-derived molecules to produce aromatic monomers. Using propylguaiacol as a model compound, we demonstrated the importance of proximity between metallic molybdenum carbide sites and the Bronsted acid sites in the zeolite in achieving high carbon yields (~80%) of benzene, toluene, propylbenzene, and phenol while maintaining catalyst stability (>98% stable conversion for 20 h). A reaction network involving both C-O and C-C bond cleavage pathways was proposed based on kinetic studies using key intermediates as feeds. Finally, we successfully depolymerized partially deoxygenated lignin oil obtained from the RCF of poplar using a continuous, two-pass catalytic process. This work highlights the potential of the bifunctional Mo/H-MFI catalyst in upgrading complex lignin feedstocks and provides a methodological approach for converting lignin-derived compounds into platform aromatic chemicals.

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

Analytical ab initio hessian from a deep learning potential for transition state optimization

Identifying transition states—saddle points on the potential energy surface connecting reactant and product minima—is central to predicting kinetic barriers and understanding chemical reaction mechanisms. In this work, we train a fully differentiable equivariant neural network potential, NewtonNet, on thousands of organic reactions and derive the analytical Hessians. By reducing the computational cost by several orders of magnitude relative to the density functional theory (DFT) ab initio source, we can afford to use the learned Hessians at every step for the saddle point optimizations. We show that the full machine learned (ML) Hessian robustly finds the transition states of 240 unseen organic reactions, even when the quality of the initial guess structures are degraded, while reducing the number of optimization steps to convergence by 2–3× compared to the quasi-Newton DFT and ML methods. All data generation, NewtonNet model, and ML transition state finding methods are available in an automated workflow.

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

Gradient-based optimization of complex nanoparticle heterostructures enabled by deep learning on heterogeneous graphs

Applications of deep learning (DL) to design nanomaterials are hampered by a lack of suitable data representations and training data. Here, in this study, we report efforts to overcome these limitations and leverage DL to optimize the nonlinear optical properties of core–shell upconverting nanoparticles (UCNPs). UCNPs, which have applications in fields such as biosensing, super-resolution microscopy and three-dimensional printing, can emit visible and ultraviolet light from near-infrared excitations. We report a large-scale dataset of UCNP emission spectra based on accurate but expensive kinetic Monte Carlo simulations (N > 6,000) and use these data to train a heterogeneous graph neural network using a physically motivated representation of UCNP nanostructure. Applying gradient-based optimization on the trained graph neural network, we identify structures with 6.5× higher predicted emission under 800-nm illumination than any UCNP in our training set. Our work reveals design principles for UCNP heterostructures and presents a roadmap for DL-based inverse design of nanomaterials.

Sivonxay, Eric [Lawrence Berkeley National Laborat↗
DOE OSTI2022

Exact constraints and appropriate norms in machine-learned exchange-correlation functionals

Machine learning techniques have received growing attention as an alternative strategy for developing general-purpose density functional approximations, augmenting the historically successful approach of human-designed functionals derived to obey mathematical constraints known for the exact exchange-correlation functional. More recently, efforts have been made to reconcile the two techniques, integrating machine learning and exact-constraint satisfaction. We continue this integrated approach, designing a deep neural network that exploits the exact constraint and appropriate norm philosophy to de-orbitalize the strongly constrained and appropriately normed (SCAN) functional. The deep neural network is trained to replicate the SCAN functional from only electron density and local derivative information, avoiding the use of the orbital-dependent kinetic energy density. The performance and transferability of the machine-learned functional are demonstrated for molecular and periodic systems.

Artificial neural networks↗
DOE OSTI2025

Challenges in integrating dissolved organic matter chemodiversity into kinetic models of soil respiration

The chemodiversity of dissolved organic matter (DOM) in soil has been proposed to influence the microbial metabolism and fate of belowground organic carbon (C). However, integrating DOM chemistry into soil C cycle models to improve predictions of C stocks and fluxes—beyond simply considering DOM pool size—remains a challenge. While recent research suggests that incorporating DOM chemodiversity into models can improve predictions of microbial respiration, there is still a lack of mechanistic understanding describing how DOM chemodiversity affects microbial metabolism and soil respiration. Here, we evaluated whether DOM chemodiversity was a determinant of soil respiration using paired measurements of high-resolution DOM chemistry, obtained from Fourier transform ion cyclotron resonance mass spectrometry (FTICR-MS), and potential soil respiration rates from across the United States (U.S.), all data provided by the Molecular Observation Network. Our objectives were to (1) assess statistical relationships between DOM chemodiversity and microbial respiration, and (2) evaluate the ability of kinetic models to leverage DOM chemistry to explain empirical relationships found in statistical models. Statistical regressions revealed that DOM chemodiversity (alpha diversity) was nonlinearly related to potential soil respiration rates, both independently and through its interactions with DOM and total C concentrations. In soils with relatively high DOM but low total C concentrations, potential soil respiration rates were negatively correlated with DOM alpha diversity, whereas in soils with relatively low DOM and high total C concentrations showed the opposite trend. However, when metabolic transition theory kinetic models were modified to include chemodiversity, their performance was comparable to traditional Monod kinetics approaches, which simulate respiration rates as a function of DOM concentration. The inability to account for nonlinearities in DOM chemodiversity–respiration relationships highlight an opportunity to advance substrate uptake kinetics by establishing causal links between DOM chemodiversity, microbial metabolism trade-offs, and potential interactions under varied environmental conditions.

Bioenergetic model↗
DOE OSTI2024

A Simple Comparison of Biochemical Systems Theory and Metabolic Control Analysis

This paper explores some basic concepts of Biochemical Systems Theory (BST) and Metabolic Control Analysis (MCA), two frameworks developed to understand the behavior of biochemical networks. Initially introduced by Savageau, BST focuses on system stability and employs power laws in modeling biochemical systems. On the other hand, MCA, pioneered by authors such as Kacser and Burns and Heinrich and Rapoport, emphasizes linearization of the governing equations and describes relationships (known as theorems) between different measures. Despite apparent differences, both frameworks are shown to be equivalent in many respects. Through a simple example of a linear chain, the paper demonstrates how BST and MCA yield identical results when analyzing steady-state behavior and logarithmic gains within biochemical pathways. This comparative analysis highlights the interchangeability of concepts such as kinetic orders, elasticities and other logarithmic gains.

FOS: Biological sciences↗
DOE OSTI2025

Functionalized Porous Polymer Networks as High-Performance PFAS Adsorbents

Toxic per- and polyfluoroalkyl substances (PFAS) are now found in nearly every water source on the planet. Exposure to these molecules can have negative health consequences, but the low concentration of PFAS relative to other solutes in water makes their removal challenging. Adsorbents offer a promising treatment route, but often exhibit low selectivities and removal capacities, as well as slow kinetics. The performance in these metrics can be improved by chemically optimizing PFAS binding sites and maximizing PFAS-adsorbent interactions. To explore how to achieve this, a porous polymer network solid (PPN-6, also known as PAF-1) was postsynthetically modified with various chemical moieties capable of leveraging unique combinations of electrostatic, hydrogen-bonding, hydrophobic, and fluorophilic interactions with PFAS molecules. Batch adsorption experiments and computational studies revealed that electrostatic and hydrogen-bonding interactions drive short-chain PFAS adsorption, while hydrophobic and fluorophilic interactions improve long-chain PFAS adsorption. In complex water matrices, a combination of electrostatic and fluorophilic interactions led to the greatest total PFAS removal. The best-performing material, functionalized with a fluorinated alkylammonium (PPN-6-FNDMB), selectively adsorbs PFAS with high capacity (up to 4.0 mmol/g) and rapid kinetics (equilibrium reached in <30 s). Furthermore, PPN-6-FNDMB outperforms several commercial adsorbents, achieving near-complete removal of 21 different PFAS from a groundwater sample collected at a US Air Force base. The PFAS could subsequently be desorbed from PPN-6-FNDMB, concentrating them by a factor of over 50 times. The recycled PPN-6-FNDMB could then be reused with minimal losses in long-chain PFAS adsorption capacity over four cycles.

Pezoulas, Ethan R↗
DOE OSTI2026

Data for KETCHUP: Parameterizing of Large-Scale Kinetic Models Using Multiple Datasets with Different Reference States

Repository for Kinetic Estimation Tool Capturing Heterogeneous Datasets Using Pyomo (KETCHUP), a flexible parameter estimation tool that leverages a primal-dual interior-point algorithm to solve a nonlinear programming (NLP) problem that identifies a set of parameters capable of recapitulating the steady-state fluxes and concentrations in wild-type and perturbed metabolic networks. KETCHUP can use K-FIT [2] input files. Example K-FIT input files are located in the K-FIT repository at https://github.com/maranasgroup/K-FIT.

Metabolomics↗
DOE OSTI2025

SympGNNs: Symplectic Graph Neural Networks for identifying high-dimensional Hamiltonian systems and node classification

Existing neural network models to learn Hamiltonian systems, such as SympNets, although accurate in low-dimensions, struggle to learn the correct dynamics for high-dimensional many-body systems. Herein, we introduce Symplectic Graph Neural Networks (SympGNNs) that can effectively handle system identification in high-dimensional Hamiltonian systems, as well as node classification. SympGNNs combine symplectic maps with permutation equivariance, a property of graph neural networks. Specifically, we propose two variants of SympGNNs: (i) G-SympGNN and (ii) LA-SympGNN, arising from different parameterizations of the kinetic and potential energy. We demonstrate the capabilities of SympGNN on two physical examples: a 40-particle coupled Harmonic oscillator, and a 2000-particle molecular dynamics simulation in a two-dimensional Lennard-Jones potential. Furthermore, we demonstrate the performance of SympGNN in the node classification task, achieving accuracy comparable to the state-of-the-art. Finally, we also empirically show that SympGNN can overcome the oversmoothing and heterophily problems, two key challenges in the field of graph neural networks.

Deep learning↗
DOE OSTI2023

Theoretical assessments of CO 2 activation and hydrogenation pathways on transition-metal surfaces

Carbon dioxide (CO 2 ) hydrogenation on transition-metal active sites offers a promising carbon utilization route toward mitigating greenhouse gas emissions. C 1 products are often formed in parallel during CO 2 hydrogenation, prompting investigations into the intrinsic properties of transition metals that drive activity and product selectivity. Here, in this work, close-packed surfaces of a selection of transition-metal catalysts (Ni, Co, Rh, Ru, Pd, and Pt) were studied with density functional theory (DFT) calculations to understand their fundamental reactivities for CO 2 transformation reactions. Results indicate that CO 2 conversion proceeds through CO* formation and hydrogenation to form C 1 products (* denotes an adsorbed species). Ni, Co, Rh, and Ru favor CO/CH 4 formation, while Pd and Pt favor CO/CH 3 OH formation. The ability of a metal to dissociate C-O bonds drives selectivity between CH 4 and CH 3 OH, while competition between CO* desorption and surface hydrogenation describes CO selectivities. The C-O bond dissociation steps often impose the highest barrier along CH 4 formation reaction profiles, suggesting their kinetic relevance for CH 4 formation rates. The provided DFT-derived data sets detail a comprehensive reaction network of elementary steps relevant to C 1 chemistries, ultimately offering a benchmark for insights into design strategies for materials that exploit transition-metal active sites in carbon capture or utilization processes.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗