Search NASA⌕ Search

SEARCH · Search NASA

Results for “reaction 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 181 records · Page 10
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Intrinsic activation energies for ring contraction of allylic cations in zeolites

Cyclic carbocations are important intermediates in zeolite-catalyzed chemistries such as methanol-to-hydrocarbons conversion, naphthenes ring opening, or coke formation. While information about their thermodynamic stability exists, little is known about the kinetics of formation and transformation of cyclic carbocations in zeolites. To fill this knowledge gap, ring contraction of the 1,3,5,5-tetramethylcyclohexenyl cation (C 10 H 17 + ), a representative of 6-membered ring allylic cations, was investigated by in situ UV–vis and IR spectroscopy. Protonic forms of zeolites served as catalysts, at temperatures from 80 °C to 135 °C. Significant oligomerization and hydride transfer in BEA and FAU hampered kinetics analysis, whereas ring contraction dominated in the channels of MOR. The reactant cation, characterized by an electronic absorption at 314 nm and an allylic stretch at 1549 cm −1 , contracted to both a 1,3-alkyl-substituted cyclopentenyl cation (287 nm and 1506 cm −1 ) and a 1,2,3-alkyl-substituted cyclopentenyl cation (297 nm and 1489 cm −1 ). Collection of time-resolved IR spectra and fitting of the intensities with various kinetic models revealed a third transformation, which is expected from thermodynamics: the 1,3-alkyl-substituted cyclopentenyl cation isomerizes to the 1,2,3-alkyl-substituted cyclopentenyl cation. Series of IR spectra recorded at different temperatures delivered intrinsic activation enthalpies (entropies) in MOR of 67 ± 2 kJ mol −1 (−130 ± 6 J mol −1 K −1 ) and 90 ± 3 kJ mol −1 (−70 ± 8 J mol −1 K −1 ) for the contraction to 1,3- and 1,2,3-substituted species, and of 88 ± 5 kJ mol −1 (−84 ± 12 J mol −1 K −1 ) for the isomerization. The findings characterize one path – via contraction of larger rings – to different cyclopentenyl species in zeolites; and the associated, moderate activation energies suggest such transformations contribute to many complex hydrocarbon reaction networks.

Acid catalysis↗
DOE OSTI2022

Machine Learning Models for Binary Molecular Classification using VUV Absorption Spectra

Machine learning methods were combined with differential absorption spectroscopy measurements in the vacuum-ultraviolet region (5.167 – 9.920 eV) in order to develop predictive capabilities for inferring molecular structure from the spectra. Several types of species were analyzed and, for modeling purposes, were defined using a single classification: (1) alkane, (2) conjugation with oxygen (e.g. diacetyl, ethyl vinyl ether), (3) non-conjugated alkene (e.g. 1-butene, 1,4-cyclohexadiene), (4) oxygen-containing (e.g. 1-butanol, tetrahydrofuran), or (5) cyclic (e.g. cyclopentane, cyclohexanone). The latter molecular classification excluded cyclic ethers. Several modeling methods were employed in the analysis of 102 absorption spectra, 24 of which were measured for the first time. The primary objective was to identify suitable methods that enable accurate predictions of molecular structure classifications with minimized statistical uncertainties. Rather than identifying a single, unifying method to reliably predict molecular structure contributions to VUV absorption spectra, coordination is required among a particular method, the type of molecular structure detail (e.g. conjugation), and absorption region of interest. The latter is accomplished using a binning approach, wherein absorption regions of ~0.5 eV were utilized rather than the entire ~4.8 eV range. Photon energy binning enabled analysis of region-specific predictions of accuracy, precision, and recall. The outcome from the binning approach is that, rather than utilizing the entire spectrum, optimal determination of molecular structure using machine learning methods depends on the absorption region. Furthermore, the present work provides separate machine learning models for each molecular classification, which enables the identification of multi-functional species relevant to atmospheric chemistry and combustion chemistry, where isomer-resolved speciation is critical to understanding complex reaction networks.

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

Probing the limits of statistical neutron capture for the r process: Experimental constraints on 141 Cs nuclear level densities

The r-process abundance peaks, particularly near mass number A ∼ 130, reflect underlying nuclear structure effects such as closed neutron shells, yet modeling the nucleosynthesis in this region remains hindered by uncertain neutron-capture rates. These rates are especially sensitive to nuclear level densities (NLDs) and γ-ray strength functions of neutron-rich nuclei, where experimental data are scarce. We present the first experimental constraint on the NLD of 141 Cs using the β-Oslo method, extending sensitivity to the neutron-rich regime near the N = 82 closed shell. Our data allow for critical calibration of microscopic NLD models and reveal that 141 Cs lies near the limit of statistical model applicability. Using this experimental input, we evaluate radiative neutron-capture rates across neighboring isotones using both Hauser–Feshbach (HF) and High Fidelity Resonance (HFR) models. Our results show order-of-magnitude rate increases for nuclei along the N = 86 line, signaling a transition to resonance-dominated capture in this region. These findings underscore the importance of constraining NLDs to improve r-process reaction network predictions, particularly in environments where the validity of statistical models breaks down.

Nuclear level density↗
DOE OSTI2020

Comparing the performance of density functionals in describing the adsorption of atoms and small molecules on Ni(111)

Reaction networks calculated using density functional theory (DFT) play a crucial role in understanding and predicting heterogeneous catalyst performance. However, results obtained using different density functionals can vary significantly. Accordingly, it is important to assess the accuracy of available functionals in capturing adsorption energetics derived from single-crystal adsorption microcalorimetry experiments. Here, we evaluate how different density functionals describe the molecular and dissociative adsorption of CH 3 I, the adsorption of CH 3 , I, and H, and the dissociative adsorption of CH 4 on Ni(111). First, we identify energetically-preferred adsorbed states for each adsorbate at three surface coverages and calculate the corresponding enthalpies of adsorption at 160 K. Then, we compare the results to experimental measurements published in the literature. For each adsorption process we find at least one functional that is quantitatively accurate within experimental error, but no functional is accurate (within experimental error) for all adsorption processes studied. Even when assuming an additional DFT error of ±20 kJ/mol on top of the experimental error, only the PBE-D3 and RPBE-D3 functionals are accurate for all considered adsorption systems. Overall, these results indicate that quantitative agreement between density functional theory calculations and experimental measurements is both system- and functional-dependent. Importantly, single-crystal adsorption microcalorimetry experiments will continue to play an important role for benchmarking density functionals.

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

A Unified Workflow for Sensitivity-Based Kinetic Analysis in Microkinetic Models

Degrees of rate control (DRC), apparent activation energies, and apparent reaction orders are established local sensitivity diagnostics for interpreting microkinetic models, but applying them routinely to large mechanisms often requires substantial reaction-specific bookkeeping, perturbation design, and postprocessing. Here, in this study, we present a unified derivative-based workflow that evaluates these quantities from a single compiled reaction-network model and target-rate definition. For any user-provided microkinetic model, the workflow compiles the mechanism into stoichiometrically consistent mass-action rate equations, solves the surface dynamics, and uses automatic differentiation to compute sensitivities with respect to rate constants, temperature, and gas partial pressures. By combining their calculations in the same framework, the workflow clearly demonstrates the relationships between different DRCs and the apparent activation energy. Using existing examples of propylene partial oxidation and methane oxidation on Pd(100), we verify expected transient redistribution of rate control, distinguish net Campbell DRCs from one-sided directional sensitivities, and show how apparent activation energy can be reconstructed either from one-sided DRCs or from state-based DRCs while critical mechanistic insights are obtained consistently. In the methane oxidation case, a pathway-subset test further illustrates how a simplified mechanism preserves key kinetic signatures of a full model, showing the potential of our user-friendly tool for model construction beyond kinetic analysis.

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 OSTI2024

Diffusion-Limited Kinetics in Reactive Systems

A proper representation of chemical kinetics is vital to understanding, modeling, and optimizing many important chemical processes. In liquid and surface phases, where diffusion is slow, the rate at which the reactants diffuse together limits the overall rate of many elementary reactions. Commonly, the textbook Smoluchowski theory is utilized to estimate effective rate coefficients in the liquid phase. On surfaces, modelers commonly resort to much more complex and expensive Kinetic Monte Carlo (KMC) simulations. Here, in this study, we extend the Smoluchowski model to allow the diffusing species to undergo chemical reactions and derive analytical formulas for the diffusion-limited rate coefficients for 3D, 2D, and 2D/3D interface cases. With these equations, we are able to demonstrate that when species react faster than they diffuse they can react orders of magnitude faster than predicted by Smoluchowski theory, through what we term “the reactive transport effect”. We validate the derived steady-state equations against particle Monte Carlo (PMC) simulations, KMC simulations, and non-steady-state solutions. Furthermore, using PMC and KMC simulations, we propose corrections that agree with all limits and the computed data for the 2D and 2D/3D interface steady-state equations, accounting for unique limitations in the associated derived equations. Additionally, we derive equations to handle couplings between diffusion-limited rate coefficients in reaction networks. We believe these equations should make it possible to run much more accurate mean-field simulations of liquids, surfaces, and liquid–surface interfaces accounting for diffusion limitations and the reactive transport effect.

Johnson, Matthew Sean↗
DOE OSTI2021

Moving from Fundamental Knowledge of Kinetics and Mechanisms on Surfaces to Prediction of Catalyst Performance in Reactors

Catalysis is defined by kinetics, but after over a century of catalyst development our kinetic and mechanistic knowledge governing real-world catalysis is still severely lacking. This Perspective considers how we can use the precise knowledge available from both theoretical and fundamental experimental techniques to understand and even predict catalytic performance under operational conditions. Here, we describe advances that link fundamental and "real-world" measurements, crucial elements required to successfully recreate a reaction network, and avenues which appear to offer viable routes to success.

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

Formic Acid Electrooxidation on Pt or Pd Monolayer on Transition-Metal Single Crystals: A First-Principles Structure Sensitivity Analysis

Here, we present a density functional theory analysis of trends for the electrooxidation of formic acid (FAO) on a single Pt or Pd monolayer supported on the close-packed (111) facet of transition metals (Pt*/M or Pd*/M): Au, Ag, Cu, Pt, Pd, Ir, and Rh, the close-packed (0001) facet of Os, Ru, and Re, and the open (100) facet of Au, Ag, Pt, Pd, Ir, and Rh. We show that the deposition of Pt or Pd pseudomorphic monolayers on these single crystals modifies the electronic structure of the Pt or Pd monolayer. Specifically, we found a direct correlation between the d-band center of the Pt and Pd monolayer and the free energy of adsorbed CO*, the latter being a reactivity descriptor for FAO. Together with the free energy of adsorbed OH* as a second reactivity descriptor, we depict the thermochemistry of the reaction network as phase diagrams showing calculated free energies across regions of rate-determining steps. We found that FAO is structure-sensitive on most surfaces studied. Pt*/Au(111) is predicted to be the most active among all Pt*/M(111/0001) surfaces studied, despite binding CO*, the strongest among the close-packed facets. This is the case because of its superiority in activating water to OH*, thus removing CO* at lower potentials than other surfaces. On similar grounds, Pt*/Pd(100), Pd*/Re(0001), Pd*/Au(111), Pd*/Ag(111), and Pd*/Pt(111) are predicted to show higher FAO activity than the corresponding monometallic Pt and Pd surfaces.

(100) facets↗
DOE OSTI2025

Low Temperature CO 2 Hydrogenation on Unsupported Mo 2 C Catalysts

CO 2 hydrogenation to methanol, a key reaction for decarbonizing the fuel and chemical industries, requires catalyst formulations that hydrogenate CO 2 selectively to methanol at temperatures where methanol conversion is not significantly equilibrium limited (<423 K). Herein we report continuous CO 2 hydrogenation at low temperatures (348-408 K, H 2 /CO 2 = 0.1-50, 5-35 bar) with high selectivity to methanol (up to ca. 80%) over unsupported β-Mo 2 C catalysts. Active site density quantification via titration with trifluoroacetic acid at reaction temperatures enables an assessment of site-specific rates. Methanation and reverse water gas shift (RWGS) occur concurrently with methanol synthesis during CO 2 hydrogenation over Mo 2 C. Reaction pathway analysis, product cofeeds, and reversibility formalisms show that all products form through primary reaction pathways from CO 2 , but secondary reactions of CO contribute significantly to rates of methanation. Dependences of forward rates on reactant and product concentration determined by independently varying the CO 2 , H 2 , CO, H 2 O, CH 3 OH, and CH 4 pressure in conjunction with reversibility formalisms reveal that all products form through H-assisted CO 2 activation and involve partially hydrogenated CO 2 -derived intermediates. Here, these inferences were verified by quantitative agreement between measured site-time yields and site-time yields predicted by closed form kinetic rate expressions in an integral reactor model over widely varying conditions (85-2000 kPa H 2 , 80-1500 kPa CO 2 , 0-45 kPa H 2 O, 0-21 kPa CO, 0-25 kPa CH 3 OH, 0-75 kPa CH 4 , 5-87 mol Mo s s mol CO 2 -1 ). Coverages calculated based on the kinetic model reveal that the Mo 2 C surface is covered with bidentate CO- and CO 2 -derived intermediates of the stoichiometry H 2 CO 2 and H 2 CO, indicating that H 2 and CO x do not compete for surface occupancy but instead adsorb cooperatively to form partially hydrogenated intermediates. Hydrogenation of the CO-derived H 2 CO** intermediate favors methanation, while hydrogenation of CO 2 -derived H 2 CO 2 ** favors methanol synthesis. Together, these findings demonstrate the ability of unsupported Mo 2 C to catalyze the hydrogenation of CO 2 to methanol at low temperatures and provide insight into the reaction network and mechanisms involved in its formation.

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 OSTI2021

Accurate Thermochemistry of Complex Lignin Structures via Density Functional Theory, Group Additivity, and Machine Learning

A molecular-level understanding of lignin structures and bond dissociation energies could facilitate depolymerization technologies. Still, this information is currently limited due to the lack of databases and the simplification of surrogate models. Here, substitution effects on seven common linkages in lignin polymers are systematically investigated. An automated reaction network generator is employed to create a database of structures. A new group additivity (GA) model based on principal component analysis (PCA) descriptors is introduced and trained on gas-phase density functional theory data of 4100 species at the M06-2X/6-311++G(d,p) level. Hydrogen bonds, local steric, and nonaromatic ring contributions are also incorporated. Lastly, we improve the accuracy of the group additivity model to reach the G4 theory by computing a data set of 770 species at this level and using a data fusion approach.

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

From Nonclassical to Classical: Crystallization Seeds Reshape Nucleation Mechanisms

Crystalline seeds are widely employed in crystallization to accelerate nucleation and control product polymorphs; yet, their impact on nucleation mechanisms remains poorly understood. While homogeneous nucleation of crystals from solution often proceeds through nonclassical pathways involving amorphous intermediates, it is unclear how seeds that promote heterogeneous nucleation reshape these mechanisms and govern polymorph selection. Here, in this study, we provide the first direct evidence that crystalline seeds can bypass the need for amorphous intermediates as nucleation sites, converting nonclassical nucleation mechanisms into classical, monomer-by-monomer crystallization pathways. Using molecular dynamics simulations of zeolite synthesis, we uncover a complex reaction network of competing nucleation processes mediated by intermediate interfacial polymorphs. The interplay between thermodynamic stability and kinetic favorability of these interfacial polymorphs dictates nucleation outcomes, creating a dynamic balance between the interfacial polymorph stability and crystallization rates. Furthermore, we show that the synthesis environment-whether monomers or aggregates serve as reactants-profoundly impacts these pathways. At moderate supersaturation, seeds eliminate amorphous intermediates and promote classical nucleation, whereas high supersaturation or aggregate-based reactants favor nonclassical pathways, even in the presence of seeds. These findings establish a general framework for understanding how seeds govern crystallization mechanisms, with broad implications for controlling nucleation kinetics, polymorph selection, and material properties. While focused on zeolites, this work reveals insights that may be applicable to biominerals, pharmaceuticals, functional materials, and catalysts, providing a basis for engineering crystallization pathways in diverse applications.

Chu-Jon, Carlos [Univ. of Utah, Salt Lake City, UT↗
DOE OSTI2025

Representing Soil Microbial Dynamics and Organo‐Mineral Interactions in the E3SM Land Model (ELM‐ReSOM)

Explicit representation of soil microbial processes and interactions with biotic and abiotic processes in Earth System Models (ESMs) remains limited, despite their importance in biogeochemical cycles. To address this gap, which hinders prediction of global biogeochemial cycling and responses to atmospheric conditions, we integrated a microbe- and mineral-surface-explicit model, the Reaction-network-based model of soil organic matter and Microbes (ReSOM), into the Energy Exascale ESM (E3SM) land model (ELM). Here, we describe ELM-ReSOM and show a case study at a conifer forest in California. ELM-ReSOM accurately simulated surface CO 2 fluxes and SOM stocks, demonstrating improved representations of microbial and mineral interactions compared to the default ELM. We examined ELM-ReSOM sensitivity to microbial traits, enzyme properties, and organo-mineral interactions. Microbial traits such as the maximum mortality rate, transporter-density scaling factor, and maximum monomer assimilation rate were strong controllers of heterotrophic respiration, while these microbial traits and enzyme-related properties collectively influenced SOM stocks. Mineral surfaces primarily affected SOM stocks by adsorbing enzymes, thereby limiting depolymerization. Synergies among processes led to stronger impacts of parameters when evaluated together versus separately (i.e., most parameters had greater indirect than direct effects). For example, due to interactions of microbial necromass with mineral surface adsorption, the indirect effect of the maximum microbial mortality rate was 33% larger than its direct effect on SOM stock. Thus, microbial and enzyme dynamics and their interactions with mineral surfaces play critical roles in SOM cycling. Tackling the challenges of microbe-explicit models will advance understanding and modeling of SOM dynamics.

Tao, Jing [Lawrence Berkeley National Laboratory (↗
DOE OSTI2024

Identifying Bayesian optimal experiments for uncertain biochemical pathway models

Abstract Pharmacodynamic (PD) models are mathematical models of cellular reaction networks that include drug mechanisms of action. These models are useful for studying predictive therapeutic outcomes of novel drug therapies in silico. However, PD models are known to possess significant uncertainty with respect to constituent parameter data, leading to uncertainty in the model predictions. Furthermore, experimental data to calibrate these models is often limited or unavailable for novel pathways. In this study, we present a Bayesian optimal experimental design approach for improving PD model prediction accuracy. We then apply our method using simulated experimental data to account for uncertainty in hypothetical laboratory measurements. This leads to a probabilistic prediction of drug performance and a quantitative measure of which prospective laboratory experiment will optimally reduce prediction uncertainty in the PD model. The methods proposed here provide a way forward for uncertainty quantification and guided experimental design for models of novel biological pathways.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Digital Twin for Chemical Science: a case study on water interactions on the Ag(111) surface

Directly visualizing chemical trajectories offers insights into catalysis, gas-phase reactions and photoinduced dynamics. Tracking the transformation of chemical species is best achieved by coupling theory and experiment. Here we developed Digital Twin for Chemical Science (DTCS) v.01, which integrates theory, experiment and their bidirectional feedback loops into a unified platform for chemical characterization. DTCS addresses a core question: given a set of experimental conditions, what is the expected outcome and why? It consists of a forward solver that takes a chemical reaction network and predicts spectra under experimental conditions, and an inverse solver that infers kinetics from measured spectra. We applied DTCS to ambient-pressure X-ray photoelectron spectroscopy measurements of the Ag–H2O interface as an example. This approach enables real-time knowledge extraction and guides experiments until a stopping condition is met based on accuracy and degeneracy. As a step toward autonomous chemical characterization, DTCS provides mechanistic knowledge in a verified, standardized manner.

Chemistry↗
DOE OSTI2026

Energy-efficient scientific computing using chemical reservoirs

The rapid growth of computing demands driven by scientific computing, data analytics, and artificial intelligence (AI) advancements has exposed the limitations of traditional digital processing systems. These systems are nearing physical energy barriers, making significant gains in energy efficiency increasingly unattainable. As we advance toward post-exascale computing, disruptive approaches are critical to overcoming these limitations. Among emerging analog solutions, biochemical computing offers a transformative path for achieving orders-of-magnitude improvements in energy efficiency. By leveraging the natural optimization capabilities of chemical reaction networks (CRNs), biochemical systems have the potential to meet high-performance computing needs through natural scalability. However, numerous challenges remain, including theoretical limitations in mapping computational problems to CRNs and practical barriers in implementing biochemical computing devices. In this paper, we present a framework for chemical computation using biochemical systems and introduce key components of our approach for energy-efficient scientific computing. We showcase the feasibility of this framework by solving a system of ordinary differential equations by emulating a chemical reservoir device, demonstrating its potential for addressing modern computing challenges. This work lays a foundational step toward harnessing the computational power of chemistry to design energy-efficient, scalable, high-performance next-generation computing systems.

Johnson, Connah G. M. [Pacific Northwest National ↗