Micro-kinetic modeling of temporal analysis of products data using kinetics-informed neural networks
Kinetics-informed neural networks improve fit quality for multi-pulse and noisy temporal analysis of products datasets.
SEARCH · Search NASA
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.
Kinetics-informed neural networks improve fit quality for multi-pulse and noisy temporal analysis of products datasets.
The Temporal Analysis of Products (TAP) pulse response technique provides valuable insights into catalytic function and reaction kinetics. However, complex fragmentation patterns in the TAP mass spectrometry signals can complicate precise quantification, particularly when analyzing transient gas flux data typical of TAP experiments. This work demonstrates a standard defragmentation method that deconvolves transient TAP signals while maintaining the temporal resolution of the experiment. First, the integrals of calibration gas fluxes are used to determine the fingerprint fragmentation pattern and construct a fragmentation matrix. This matrix is then used to defragment experimental flux data at each recorded time point via a non-negative least squares regression. The effectiveness of this method is demonstrated using virtual data and control experiments with a TAP reactor system. The defragmentation is then applied to the more complex propane dehydrogenation reaction on a chromia/alumina catalyst, which can contain up to ten significant gas species in the reactor outlet. Initial propane pulsing reveals an induction period during which propane is fully oxidized to CO2, followed by partial reduction to CO. Afterwards, there is a transition in chemistries towards coking and propylene production. Our example illustrates a practical method for the accurate determination of the time-dependent reactant/product concentrations and rates for a thorough analysis of the propane dehydrogenation kinetics. This approach can be broadly applied to any transient mass spectrometry experiment for a better understanding of catalyst-reaction dynamics.
This perspective reviews the current understanding of the Oxidative Coupling of Methane (OCM) reaction over the supported Mn-Na 2 WO 4 /SiO 2 catalyst, with a focus on recent insights gained from state-of-the-art in-situ and operando spectroscopic characterization and chemical probe experiments under controlled environments. The supported Mn-Na 2 WO 4 /SiO 2 catalyst exhibits dynamic structural changes during the OCM reaction, involving multiple reactive lattice and adsorbed oxygen species, each associated with different oxide phases. These oxygen species play distinct roles in various steps of the OCM mechanism. The catalytic active sites for activation of CH 4 are associated with isolated surface Na-WO x sites on the SiO 2 support and the role of surface MnO x sites on SiO 2 is to oxidatively dehydrogenate C 2 H 6 to C 2 H 4 . Furthermore, this paper provides a detailed discussion of these roles and also introduces new experimental data from Temporal Analysis of Products (TAP) studies to clarify the ongoing debate in the literature regarding the contributions of lattice versus adsorbed oxygen species in OCM reaction product formation. Additionally, recommendations are offered for optimizing the performance of supported Mn-Na 2 WO 4 /SiO 2 catalysts to enhance CH 4 activation and C 2 product selectivity.
Selective oxidation of propene to acrolein over industrial multicomponent bismuth molybdate (BMO) catalysts significantly depends on reaction conditions that include operating parameters and catalyst state. Here, this work investigates selective propene oxidation in the intrinsic kinetic paradigm of Temporal Analysis of Products (TAP) reactor by systematically varying catalyst redox state, temperature, and oxygen-to-propene feed ratio. A 93 % propene conversion with an acrolein yield of 80 % is achieved at elevated temperatures (450 °C) on oxidized catalysts under oxygen-rich conditions (O 2 :C 3 H 6 = 10). However, these conditions diminish acrolein-to-CO 2 selectivity due to enhanced total oxidation to CO 2 . In contrast, a reduced catalyst state, moderate temperature (350 °C), and lower oxygen feed (O 2 :C 3 H 6 = 1) nearly doubles the acrolein to-CO 2 selectivity, albeit at a lower acrolein yield (33 %). Transient kinetic studies together with a kinetic model that simplify the major product formation pathways in lumped non-elementary forms reveal that the availability of surface oxygen species plays a pivotal role in governing reaction pathways. Additionally, density functional theory (DFT) calculations on pure BMO catalysts inform the role of surface redox states on propene and oxygen activation barriers. Lattice oxygen at acrolein-selective sites drives both acrolein and CO 2 formation, while adsorbed oxygen at activation sites favors unselective CO 2 generation. This work establishes a critical relationship between transient product selectivity and surface oxygen availability, which is strongly influenced by catalyst redox state, feed ratio, and reaction temperature. These insights underscore the importance of dynamic reactor operation strategies and offer a foundation for designing next-generation propene oxidation processes with tunable acrolein selectivity.
Here, this paper introduces the use of statistical distributions based on transport differential equations for clear distinction of transport modes within transient kinetic experiments. More specifically, novel techniques are developed for the transient data obtained through the Temporal Analysis of Products (TAP) reactor and are applicable to experiments where pulse response into a gas flow is used. The methodology allows distinguishing between two domains of diffusion transport in heterogeneous catalytic systems, i.e., Knudsen and non-Knudsen diffusion, using statistical fingerprints, and finding the transition domain. Two distribution parameters were obtained that directly result in coefficients that correspond to the concentration and the rate of transport. Using a linear relationship between the rate and concentration coefficients, Knudsen diffusion is revealed when the rate of transport is constant and non-Knudsen diffusion is confirmed when the rate of transport coefficient is a function of the concentration coefficient. As a result, accurate transport information can be extracted from experimental data even in the presence of comprising instrument drift or noise particularly when analyzing higher pressure pulse responses with complex transport. This enables more direct investigation of experiments influenced by gas-phase reactions.
This software implements a complete preprocessing pipeline for transient mass spectrometry (MS) data collected during TAP (Temporal Analysis of Products) experiments. It is designed to extract chemically meaningful fluxes from overlapping ion signals by applying a calibrated defragmentation matrix and solving the resulting linear system using non-negative least squares (NNLS) regression. The core script, preprocess_mass_spec.py, performs the following operations: Gain correction: Applies amplifier gain scalars derived from inert-packed calibration pulses to normalize signal intensities across AMUs and acquisition settings. Background subtraction: Removes experiment baselines using user-defined time windows, ensuring compatibility with slow-diffusing species and preventing negative values that would interfere with NNLS. Options to subtract before and after defragmentation. Defragmentation: Constructs a fragmentation matrix A from zeroth moments of calibration pulses (equal molar gas:inert mixtures) and solves Ax=b at each time point, where b is the raw MS signal and x is the estimated species flux. The matrix is normalized to inert signals and accounts for instrument-specific fragmentation behavior. Pulse-mode handling: Supports both averaged and individual pulse modes, enabling statistical treatment of fluxes and calculation of standard deviations. Integration and output: Computes zeroth moments (integrated fluxes) and exports time-resolved and integrated data in CSV format, suitable for downstream kinetic modeling. The software is validated using both virtual TAP simulations (VTAP) and experimental data from propane dehydrogenation (PDH) on CrOx/Al2O3 catalysts. It preserves temporal resolution by applying NNLS point-by-point across the pulse duration (typically 6,000+ time slices per pulse), leveraging the linear superposition principle to reconstruct full flux profiles. The defragmented outputs are compatible with kinetic extraction methods such as the G and Y procedures, which are used to derive rate–concentration relationships from TAP data. The details of these validations are discussed in detail in the supporting manuscript and supporting information. Example data and output files are also included. The methodology is robust to experimental noise and drift, with calibration protocols that account for pulse size effects, MS aging, and inert gas normalization. The software is modular, reproducible, and tailored for high-throughput TAP-MS workflows in catalysis research.
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.
Direct study of changes in catalyst structure-kinetic properties during redox transitions using a new time-resolved technique The Temporal Analysis of Products (TAP) pulse response methodology1 is a transient technique offering the time scale needed to deconvolve many reaction steps from the complex network typical to industrial catalytic processes. The traditional TAP measurement observes the gas phase dynamic response but, hereto now, has been devoid of any direct measurement of change in the catalyst itself. The development of an operando technique that couples gas phase transient kinetics to dynamic metal centers and surface species is presented. Significance Unification of the TAP methodology with time-resolved spectroscopic measurement can offer unprecedented insight into the complex kinetic phenomena regulated by the solid catalyst.
The Temporal Analysis of Products (TAP) pulse response methodology is a transient technique that provides the time resolution needed to deconvolve reaction steps from the complex networks typical in industrial catalytic processes. Traditionally, TAP measurements observe gas phase dynamics at the reactor exit but lack direct measurements of changes in the catalyst itself. Recently, a new operando technique was developed that couples gas phase transients to dynamic changes in metal centers with the precise TAP methodology for nanomole titration. Using an industrial CrOx/Al2O3 catalyst used for propane dehydrogenation, we demonstrate the capabilities of this unique device to reveal key catalytic processes: 1) total propane oxidation not associated with chromia centers, 2) reduction of Cr6+ to Cr3+ correlated with selective product formation, and 3) subsequent carbon accumulation. By utilizing incremental pulsing in a diffusion-only transport regime, the spectrokinetic device allows us to resolve detailed changes in catalyst structure, composition, and kinetic function that are otherwise indistinguishable in conventional operando devices. The unification of the TAP methodology with time-resolved spectroscopic measurements offers new and unique insights into the complex kinetic phenomena regulated by solid catalyst surfaces.
The CO 2 methanation reaction (aka the Sabatier reaction) is one of the important carbon utilization technologies to reduce atmospheric greenhouse gases. In this study, kinetic modeling of the Sabatier reaction proceeding via the CO route has been presented by considering weak and strong chemisorption concepts rooted in electronic theory. Based on the single-route mechanism of the Sabatier reaction, the dependence of the reaction rate on the position of the catalyst Fermi level with respect to the antibonding states of the reagent (CO 2 ) and the product (CH 4 ), is demonstrated with the classic Sabatier volcano-shape. Furthermore, the model was applied to explain experimental results of reaction rates for Ni- and Ru-based catalytic systems and to understand origin of enhancement of turnover frequencies established for dynamic catalysis.
A brief literature review of recent experimental and computational efforts on the use of high entropy alloys (HEAs) as catalysts is presented, while also sharing some perspectives and future insights. To fully broach the vast compositional possibilities of HEA materials, integrating computational modeling with high-throughput experimental synthesis and validation is necessary to accelerate their design and development. Once identified, specific HEAs can be a class of materials for the next generation of catalysts when addressing global challenges related to energy independence, commodity chemical production, environmental remediation, abating emissions, and modernized domestic supply chain resilience.
Transient kinetic data contain a wealth of information about intrinsic features of a catalyst as well as the reaction mechanism. Currently, high volume transient data is underutilized, and data science methods could both increase the value of information that can be extracted from this data, integrate experimental with theoretical data sources, and accelerate the pace of catalyst technology advancement. Transient kinetic characterizations with simple probe molecules exhibiting reversible adsorption, irreversible adsorption and bulk-surface diffusion are presented as training components for similar experiments with more complex surface reactions. In conclusion, by increasing the availability and accessibility of transient kinetic data through details of its structure and acquisition, we aim to decrease the barrier for data scientists to apply machine learning methods to this valuable data source.
ABSTRACT Monoclonal antibodies (mAbs) represent a majority of biotherapeutics in the market today. These glycoproteins undergo posttranslational modifications, such as N‐linked glycosylation, that influence the structural & functional characteristics of the antibody. Glycosylation is a heterogenous posttranslational modification that may influence therapeutic glycoprotein stability and clinical efficacy, which is why it is often considered a critical quality attribute (CQA) of the mAb product. While much is known about the glycosylation pathways of Chinese Hamster Ovary (CHO) cells and how cell culture chemical modifiers may influence the N‐glycosylation profile of the final product, this knowledge is often based on the final cumulative glycan profile at the end of the batch process. Building a temporal understanding of N‐glycosylation and how mAb glycoform composition responds to real‐time changes in the biomanufacturing process will help build integrated process models that may allow for glycosylation control to produce a more homogenous product. Here, we look at the effect of specific nutrient feed media additives (e.g., galactose, manganese) and feeding times on the N‐glycosylation pathway to modulate N‐glycosylation of a Herceptin biosimilar mAb (i.e., Trastuzumab). We deploy the N‐GLYcanyzer process analytical technology (PAT) to monitor glycoforms in near real‐time for bench‐scale bioprocesses operated in both fed‐batch and perfusion modes to build an understanding of how temporal changes in mAb N‐glycosylation are dependent on specific media additives. We find that Trastuzumab terminal galactosylation is sensitive to media feeding times and intracellular nucleotide sugar pools. Temporal analysis reveals an increased desirable production of single and double galactose‐occupied glycoforms over time under glucose‐starved fed‐batch cultures. Comparable galactosylation profiles were also observed between fed‐batch (nutrient‐limited) and perfusion (non‐nutrient‐limited) bioprocess conditions. In summary, our results demonstrate the utility of real‐time monitoring of mAb glycoforms and feeding critical cell culture nutrients under fed‐batch and perfusion bioprocessing conditions to produce higher‐quality biologics.
Catalyst size, morphology, and crystal structure play crucial roles in determining the activity and selectivity of electrochemical CO 2 reduction reactions, which are known to change during the reaction process. A comprehensive understanding of how, when, and why these parameters evolve under operational conditions is essential for developing stable, efficient, and selective catalysts. In this study, we reveal that formate, one of the reaction products, contributes to the degradation of copper catalysts through a ligand-assisted dissolution mechanism. Utilizing in situ electrochemical atomic force microscopy and ex-situ scanning and transmission electron microscopies, we observed a significant reduction in the size of copper nanoparticles, which decreased from over 30 nm to less than 10 nm in diameter within 60 min of CO 2 RR. The temporal production of formate correlated with the particle size changes. Furthermore, analysis of the electrolyte using inductively coupled plasma optical emission spectroscopy confirmed the dissolution of copper nanoparticles. Control experiments involving various reaction products (H 2 , CO, and HCOO – ) demonstrated that formate significantly promotes copper dissolution, thereby highlighting its role in the ligand-assisted dissolution mechanism of copper electrocatalysts. In conclusion, our findings provide critical insights into copper catalyst behavior during electrochemical CO 2 reduction, facilitating the design of more resilient and effective electrocatalysts.
Abstract. This study explores the causal relationships between catchment water availability, vapor pressure deficit, and gross primary productivity (GPP) across 341 catchments in the contiguous US. Seasonal climatic, hydrological, and vegetation characteristics were represented using the Horton index, ecological aridity index, evaporative fraction index, and carbon uptake efficiency. Statistical methods, including circularity statistics, correlation analysis, and causality tests, were employed to determine the complex interactions between catchment wetness, atmospheric dryness, and vegetation carbon uptake. The results revealed a maximum lag of 2 months in the intra-annual variability of catchment water supply–productivity and atmospheric water demand–productivity relationships, with hysteresis patterns varying with the catchment's hydrological characteristics. In catchments not permanently under water-limited or energy-limited conditions, vegetation experiences hydrological stress during the peak growing period, coinciding with the highest gross primary productivity and carbon uptake efficiency being out of phase with the Horton index and in phase with the evaporative fraction index. Causality analysis highlights strong temporal continuity in GPP seasonal characteristics, with a cause–effect relationship between catchment water supply, atmospheric demand, and vegetation productivity spanning a maximum of 2 months. These findings underscore the need for a comprehensive functional framework that integrates catchment water supply, atmospheric demand, and vegetation productivity to enhance our understanding and predictive capabilities with regard to ecosystem responses to climate change.
Abstract Background The industrial feasibility of photosynthetic bioproduction using cyanobacterial platforms remains challenging due to insufficient yields, particularly due to competition between product formation and cellular carbon demands across different temporal phases of growth. This study investigates how circadian clock regulation impacts carbon partitioning between storage, growth, and product synthesis in Synechococcus elongatus PCC 7942, and provides insights that suggest potential strategies for enhanced bioproduction. Results After entrainment to light-dark cycles, PCC 7942 cultures transitioned to constant light revealed distinct temporal patterns in sucrose production, exhibiting three-fold higher productivity during subjective night compared to subjective day despite moderate down-regulation of genes from the photosynthetic apparatus. This enhanced productivity coincided with reduced glycogen accumulation and halted cell division at subjective night time, suggesting temporal separation of competing processes. Transcriptome analysis revealed coordinated circadian clock-driven adjustment of the cell cycle and rewiring of energy and carbon metabolism, with over 300 genes showing differential expression across four time points. The subjective night was characterized by altered expression of cell division-related genes and reduced expression of genes involved in glycogen synthesis, while showing upregulation of glycogen degradation pathways, alternative electron flow components, the pentose phosphate pathway, and oxidative decarboxylation of pyruvate. These molecular changes created favorable conditions for product formation through enhanced availability of major sucrose precursors (glucose-1-phosphate and fructose-6-phosphate) and maintained redox balance through multiple mechanisms. Conclusions Our analysis of circadian regulatory rewiring of carbon metabolism and redox balancing suggests two potential approaches that could be developed for improving cyanobacterial bioproduction: leveraging natural circadian rhythms for optimizing cultivation conditions and timing of pathway induction, and engineering strains that mimic circadian-driven metabolic shifts through controlled carbon flux redistribution and redox rebalancing. While these strategies remain to be tested, they could theoretically improve the efficiency of photosynthetic bioproduction by enabling better temporal separation between cell growth, carbon storage accumulation, and product synthesis phases.
The purpose of this experiment was to evaluate how circadian clock regulation impacts carbon partitioning between storage, growth, and product synthesis in Synechococcus elongatus PCC 7942 in providing insights to strategies for enhanced bioproduction. Culture samples were collected at 0, 0.5, 1, 2, 4, 6, and 8 hours for extracellular sucrose analysis. Circadian metabolomics data was acquired using a Agilent single quadrupole gas chromatography-mass spectrometer and processed using Agilent Mass Hunter for targeted sucrose quantification. Metabolomic analysis of PCC 7942 light-dark cycle cultures transitioned to constant light revealed distinct temporal patterns in sucrose production. Processed metabolomic datasets are openly accessible from the PNNL DataHub project dataset download page and contain secondary processed GC-MS results files and supporting metadata materials linked to relevant source code information supporting data transparency and reuse.
Oil and natural gas (ONG) extraction emits volatile organic compounds (VOCs). Certain VOCs are identified as hazardous air pollutants (HAPS) while others contribute to ozone formation. This study examines the impact of ONG operations on VOC levels during the development of multi-well ONG pads in suburban Broomfield, Colorado. From October 2018 to December 2020, weekly VOC measurements were taken at 18 sites across the area. These included spots near well pads, in adjacent neighborhoods, and at a background site, covering various stages of well pad development including drilling, hydraulic fracturing, flowback, and production. Analysis using Positive Matrix Factorization (PMF) identified six factors, including combustion, background/biogenic sources, light and complex alkanes, drilling activities, and ONG acetylene. Factors linked to local ONG activities exhibited clear temporal and spatial correlations with Broomfield well development. Benzene source analysis revealed distinct contribution gradients, with ONG-related sources notably influencing areas near the well pads, particularly in pre-production. ONG-related weekly benzene contributions varied from 9% to 63% at a community background site and 18% to 89% in a neighborhood close to a well pad.