Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “kinetic pathways”

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 OSTI2022

Deciphering the Mechanistic Role of Individual Oxide Phases and Their Combinations in Supported Mn–Na 2 WO 4 Catalysts for Oxidative Coupling of Methane

Oxidative coupling of methane (OCM) is an attractive direct route for upgrading methane to valuable chemicals. In this study, Temporal Analysis of Products (TAP) and steady state experiments are conducted to understand the role of individual oxide phases and their combinations in supported Mn-Na 2 WO 4 /SiO 2 catalysts for OCM. The results from TAP transient kinetic studies indicate that Mn plays an important role in promoting gas phase oxygen activation, while NaO x /SiO 2 and WO x /SiO 2 are relatively inert towards gas phase oxygen and methane activation. However, the supported catalyst combining Na and W in the form Na 2 WO 4 show enhanced gas phase oxygen activation exhibiting a much lower oxygen activation energy (148 kJ/mol) and enhanced activity toward methane activation as compared to the individual supported oxide catalysts. Addition of Mn to Na 2 WO 4 /SiO 2 further decreases the oxygen activation energy by 40 kJ/mol. Moreover, methane activation is also enhanced with CH 3 as the main intermediate but, with increasing Mn content, more CH 2 intermediates are observed. Different forms of oxygen (both dioxygen and atomic) are detected on the catalyst surface using isotopic pump/probe pulsing and their distribution is found to depend on the catalyst composition. An optimal Mn content in the Na 2 WO 4 /SiO 2 catalyst system is needed to enhance the amount of dioxide surface species (e.g., superoxide 16 O 2 - or peroxide 16 O 2 2- ) associated with the Na 2 WO 4 leading to high C 2 selectivity for OCM. When the Mn content is too high, the larger MnOx domains are shown to contribute to the formation of higher concertation of monoxide surface species that lead to nonselective OCM pathways. This insight from transient kinetic characterization using TAP combined with conventional steady state studies, provides a deeper understanding of the role of individual oxide phases and their combination on supported catalysts toward the formation of intermediate surface species and their impact on the OCM reaction mechanism. This knowledge is critical toward designing superior catalyst formulations for OCM.

10 SYNTHETIC FUELS↗
DOE OSTI2025

Epoxide Alcoholysis over M-BEA Zeolites: Effects of Alcohol Chain Length on Rates and Regioselectivities

The structures of nucleophilic reactants affect their coordination behavior among solvent molecules and kinetics of reactions with surface intermediates within the confines of fluidfilled pores of zeolites and other microporous materials. Consequently, rates and regioselectivities of diverse chemistries may depend sensitively on nucleophile identity in manners not observed for classic fluid phase reactions. Here, we examine the impact of varying the primary alcohol (ROH) chain length on the kinetics of 1,2-epoxybutane (C 4 H 8 O) ring-opening within Brønsted (Al-BEA) and Lewis acid (Zr-BEA) zeolites. Turnover rates increase by factors of ∼6 (Al-BEA) and 4-fold (Zr-BEA) between methanol and 1-hexanol, yet the reaction mechanisms remain comparable. Despite modest rate differences, apparent activation enthalpies calculated from rates and activities of solvated reactants decrease linearly by 12 (Al-BEA) to 33 kJ mol −1 (Zr- BEA) with increased proton affinity, which suggests bond formation energies for the nucleophile strongly influence rate increases. The molecular interpretation of these trends demonstrates, however, that the solvation of ring-opening transition states by zeolite pore structures and solvent molecules also governs rates. The impact of local solvating interactions appears most directly as changes in regioselectivities, which tend to enhance terminal alcohol formation with increasing ROH chain length. Regioselectivities largely do not vary with differences in fluid composition for a given ROH. The addition of H 2 O increases the number of hydrogen bonds among reactive species, and trends in regioselectivities imply that the decreased hydrogen bonding ability of longer chain ROH, and not the nucleophile strength or steric bulk, determines the regioselectivities of the resulting products. This work provides direct experimental evidence that nucleophilicity and hydrogen bonding influence reaction barriers and regioselectivities in zeolitecatalyzed epoxide ring-opening, offering pathways to better control reaction kinetics.

acidic zeolites↗
DOE OSTI2021

Microstructural evolution in Cu-Nb processed via friction consolidation

Immiscible alloys, whether in well-mixed or layered forms, are of increasing interest based on their novel structural and functional properties, such as enhanced thermal stability against grain growth or radiation-induced defect trapping at the interfaces. To address the need for new approaches to tailor microstructures, the microstructural development of an immiscible Cu-Nb alloy processed via friction consolidation of elemental powders is investigated. Friction consolidation is a solid-phase processing technique that imparts severe plastic strain into a deforming volume resulting in elevated temperatures below the melting temperature of the alloy. Two distinct processing pathways were chosen to understand the effect of thermomechanical conditions on the final microstructure. The microstructure was characterized using scanning electron microscopy, scanning transmission electron microscopy, and X-ray diffraction techniques. Path 1 exhibited larger strain, strain rate, and temperature as compared with path 2. In path 1, agglomerated Nb particles were present in the recrystallized ultrafine-grained Cu matrix, while in path 2 extremely fine and dispersed Nb particles were present in a highly deformed Cu matrix. In both pathways, supersaturation of Cu in Nb lattices was noted, but not vice-versa. The asymmetry in mixing is explained based on deformation-based, thermodynamic and kinetic factors. These findings provide a pathway for creation of novel tailored microstructures and improved properties in any number of binary immiscible alloy systems.

Cu-Nb alloys, Microstructural refinement, Friction↗
DOE OSTI2022

Mechanistic pathways and role of oxygen in oxidative coupling of methane derived from transient kinetic studies

Oxidative coupling of methane (OCM) is a promising industrial process to upgrade natural gas to high value chemicals. In this study, Temporal Analysis of Products (TAP) and steady-state experiments were conducted to distinguish how the composition of surface and gas phase oxygen influence mechanistic details of the selective conversion of CH 4 to C 2 H 4 over the Mn-Na 2 WO 4 /SiO 2 catalyst. The results from TAP studies indicate that methane activation on this catalyst proceeds predominantly via a short-lived, transient surface oxygen species and there is a competition for this species to form either CO or methyl radicals on the surface. This active species has a total lifetime of 3 s and is identified to have a dioxygen (e.g. O 2 2- or O 2 -) form. We show that the concentration of the transient surface oxygen species significantly impacts the OCM performance. Oxygen attributed to the catalyst lattice (in a singular form e.g., O - ), is found to activate methane to a lesser degree, but exclusively forms CO 2 . Evidence for these surface pathways for methyl radical, CO and CO 2 formation identified by TAP are also validated through steady-state experiments. Finally, by distinguishing different catalyst oxygen species and their role in selective/nonselective pathways, important screening criteria have been identified for the advancement of superior catalyst formulations.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Stimulus-Responsive Modulation of Solvation Environments in Solid Catalysts

Liquid environments play a crucial role in the biological processes occurring in living organisms as well as in many human-made processes involving electrochemistry, photo-, and thermocatalysis. In the majority of these systems, aqueous phases are ubiquitous due to water’s natural abundance. Water molecules, however, can exert large changes in the chemical environment of catalytically active sites, altering the reaction rates, selectivity, and catalyst stability. These solvation effects induced by water molecules near catalytic sites can drastically change the energy landscape and unlock unique reaction pathways with far more favorable kinetics. In nature, living organisms couple these complex interactions with detection, communication, and actuation mechanisms to induce self-regulatory behavior, ensuring stability of the system and thus long-term durability. Extrapolating this behavior to heterogeneous catalysis is desirable because the resulting “smart materials” can potentially unlock new chemical conversion processes with higher atom efficiency, rates, and stability. The combination of polymer chemistry and heterogeneous catalysis has introduced versatile approaches for creating materials that can respond to cues in the reaction medium that alter the accessibility, intrinsic activity, and selectivity of the catalyst. To achieve this, one could combine stimulus-responsive polymers, which undergo a large volumetric phase transition in response to an external stimulus, with a solid catalyst. This chemo-mechanical response has been employed to create a variety of nanoreactor vessels with stimulus-responsive character that turn on- and off- depending on the reaction conditions. In this Account, we focus on the impact of these polymer coatings on the solvation environment around the active site and the implications of these effects on the reaction energy landscape, molecular arrangement of the solvent, electric fields at the catalyst–liquid interface, binding energy, and mobility of surface reaction intermediates. These seemingly subtle changes in solvent molecules induced by the presence of polymers can have a tremendous impact on the development of bioinspired heterogeneous catalysts, reliable chemical clocks, micro/nanoreactors, and robots. The large library of polymer chemistries offers a plethora of combinations of stimulus-responsive mechanisms (e.g., temperature, pH, light, magnetic field, solvent composition), providing the possibility of creating homeostatic catalysts à la carte.

catalysts↗
DOE OSTI2021

Measurements of Intermediate Species in Fuel-Rich Oxidation of Ethylene, Toluene, and n -Decane

Polycyclic aromatic hydrocarbons (PAHs) are important precursors to the formation of soot and are also pollutant emissions from combustion devices, including internal combustion engines. To understand the formation of PAHs, the low-molecular-weight hydrocarbon chemistry leading to their formation needs to be understood. Toward this goal, fuel-rich oxidation of three different hydrocarbons, i.e., ethylene, toluene, and n-decane, has been investigated in an atmospheric pressure flow reactor at varying temperatures (1000–1350 K), equivalence ratios (φ = 3.0–12.0), and residence times of 0.25–1.5 s. The major C 1 to C 7 intermediates, such as methane, acetylene, ethylene, allene, propyne, propylene, vinylacetylene, 1,3-butadiene, cyclopentadiene, benzene, and toluene, were quantified using a gas chromatograph equipped with a flame ionization detector. The experimental tendency of intermediate species formation, e.g., the dependence of temperature, equivalence ratio, and residence time, was similar in ethylene oxidation and n-decane oxidation. The concentrations of intermediate species up to C 4 were higher in ethylene/n-decane oxidation than toluene oxidation, while a nearly equal or larger amount of C5–C7 species was produced in toluene oxidation. The experimental data were compared with modeling results using a detailed chemical kinetic mechanism. The calculated data using the kinetic model were in agreement with the experimental results. Additionally, a comprehensive kinetic analysis on the reaction pathways of each species was conducted to assess the differences in the oxidation chemistry with the change in the structure of hydrocarbons. In particular, in-depth analysis of benzene formation was performed, elucidating that two benzene formation pathways were important in ethylene and n-decane oxidation: (1) H elimination of the 2,4-cyclohexadienyl radical produced from the isomerization of the 2,4-cyclopentadienylmethyl radical, and (2) dehydrogenation of cyclohexadiene produced from reactions of vinyl radical with 1,3-butadiene. In toluene oxidation, it was found that benzene was primarily produced from toluene through the replacement of the methyl group with hydrogen.

30 DIRECT ENERGY CONVERSION↗
DOE OSTI2026

Cation Site Occupancy in Natural Ferroan Double Carbonates via Mössbauer and Fe K-Edge X-ray Absorption Spectroscopy

Double carbonates are minerals with a calcite type structure with alternating cation layers composed of Ca and Mg/Fe coordinated by carbonate groups. While the perfectly ordered AB-stacked crystal is the thermodynamically most stable configuration, natural mineral formation pathways can leave signatures through kinetically trapped disorder such as AB antisite cation substitutions. This study probes the degree of cation ordering in naturally occurring double carbonate samples.in particular, the distribution of Fe A and B crystallographic sites (A site is Ca layer and B site is Mg/Fe layer). Mössbauer, X-ray diffraction, X-ray absorption fine structure (both X-ray absorption near edge structure, XANES, and extended x-ray absorption fine structure, EXAFS), and energy dispersive X-ray spectroscopies were used to collect a comprehensive experimental data set, which we interpret using density functional theory to elucidate the structural effects of cation disorder. Our results show that the A (nominally Ca) site can host a relatively high Fe fraction. We discuss the implications in terms of mineral formation.

Boglaienko, Daria (ORCID:0000000160943877)↗
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 OSTI2023

Constraining the Reaction Rate of Criegee Intermediates with Carboxylic Acids during the Multiphase Ozonolysis of Aerosolized Alkenes

Criegee intermediates (CIs) are central to the ozonolysis of unsaturated organic compounds, controlling the formation of small fragmentation reaction products and higher molecular weight oligomeric compounds through competing unimolecular and bimolecular pathways. In particular, the reaction of thermalized CI with carboxylic acids (RCOOH) to produce oligomeric α-acyloxyalkyl hydroperoxides (AAHPs) is of interest to the chemical transformation of organic compounds in the atmosphere. In the gas phase, the CI + RCOOH reaction proceeds at supercollisional rates, consistent with a barrierless reaction mechanism. However, the rate of this reaction in condensed phases remains uncertain. Here, we report the results of aerosol flow tube studies of the ozonolysis of submicron organic aerosol containing mixtures of an internal n-alkene (Z-9-tricosene, tri) and a saturated acid (2-hexyldecanoic acid, HDA). The decay of the reactants and formation of reaction products are monitored using atmospheric pressure chemical ionization mass spectrometry (APCI–MS). A 6-fold decrease in the reactive uptake of ozone is observed with increasing acid concentration, consistent with a previous study using alcohol additives. The decay of HDA with ozone exposure shows that HDA scavenges CIs at rates comparable to competing unimolecular pathways involving CIs. The decay kinetics of HDA are analyzed using a simple kinetic model in order to constrain the ratio of the CI + RCOOH and unimolecular CI rate constants. Empirical fitting to this model, combined with theoretical estimates of the unimolecular loss rate, yields a CI + RCOOH rate constant of 1.85 ± 0.27 × 10 –19 cm 3 molec –1 s –1 , which is six orders of magnitude slower than the expected diffusion limit in a liquid organic matrix. Stochastic kinetic simulations are conducted with different values of the CI + RCOOH rate constant to validate this result and show that the diffusion-limited rate is indeed much too rapid to reproduce the experimental results.

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

Stable excited dication: trapping on the S 1 state of formaldehyde dication after strong field ionization

Combined theoretical and experimental work examines the dynamics of dication formaldehyde produced by strong field ionization. Trajectory surface hopping dynamics on the first several singlet electronic states of the formaldehyde dication are used to examine the relaxation pathways and dissociation channels, while kinetic energy distributions after strong field ionization of formaldehyde and deuterated formaldehyde are used to confirm the theoretical predictions. Here, we find that the first excited state of the formaldehyde dication is stable, neither decays to the ground state nor dissociates, even though the ground state and higher lying states are directly dissociative. The stability of the first excited state is explained by its symmetry which does not allow for radiative or nonradiative transitions to the ground state and by large barriers to dissociate on the excited state surface.

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

Reaction dynamics of S( 3 P) with 1,3-butadiene and isoprene: crossed-beam scattering, low-temperature flow experiments, and high-level electronic structure calculations

Sulfur atoms serve as key players in diverse chemical processes, from astrochemistry at very low temperature to combustion at high temperature. Building upon our prior findings, showing cyclization to thiophenes following the reaction of ground-state sulfur atoms with dienes, we here extend this investigation to include many additional reaction products, guided by detailed theoretical predictions. The outcomes highlight the complex formation of products during intersystem crossing (ISC) to the singlet surfaces. Here, we employed crossed-beam velocity map imaging and high-level ab initio methods to explore the reaction of S( 3 P) with 1,3-butadiene and isoprene under single-collision conditions and in low-temperature flows. For the butadiene reaction, our experimental results show the formation of thiophene via H 2 loss, a 2H-thiophenyl radical through H loss, and thioketene through ethene loss at a slightly higher collision energy compared to previous observations. Complementary Chirped-Pulse Fourier-Transform mmWave spectroscopy (CP-FTmmW) measurements in a uniform flow confirmed the formation of thioketene in the reaction at 20 K. For the isoprene reaction, we observed analogous products along with the 2H-thiophenyl radical arising from methyl loss and C 3 H 4 S (loss of ethene or H 2 + acetylene). CP-FTmmW detected the formation of thioformaldehyde via loss of 1,3-butadiene, again in the 20 K flow. Coupled-cluster calculations on the pathways found by the automated kinetic workflow code KinBot support these findings and indicate ISC to the singlet surface, leading to the generation of various long-lived intermediates, including 5-membered heterocycles.

Lang, Jinxin↗
DOE OSTI2023

High-pressure CO 2 dissociation with nanosecond pulsed discharges

Herein we investigate the conversion of CO 2 into CO and O 2 with nanosecond repetitively pulsed (NRP) discharges in a high-pressure batch reactor. Stable discharges are obtained at up to 12 bar. By-products are measured with gas chromatography. The energy efficiency is determined for a range of processing times, pulse energy, and fill pressures. It is only weakly sensitive to the plasma operating parameters, i.e the extent of CO 2 conversion is almost linearly-dependent on the specific energy invested. A conversion rate as high as 14% is achieved with an energy efficiency of 23%. For long processing times, saturation in the yield and a drop in efficiency are observed, due to the increasing role of three-body recombination reactions, as described by zero-dimensional detailed kinetic modeling. The modeling reveals the presence of three-stage kinetics between NRP pulses, controlled by electron-impact CO 2 dissociation, vibrational relaxation, and neutral elementary kinetics. Transport effects are shown to be important for CO 2 conversion at high pressures. For fill pressures beyond 10 bar, CO 2 may locally transit into supercritical states. The supercritical plasma kinetics may bypass atomic oxygen pathways and directly convert CO 2 into O 2 . This work provides a detailed analysis of plasma-based high-pressure CO 2 conversion, which is of great relevance to future large-scale sustainable carbon capture, utilization, and storage.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Controlled Reaction Dynamics of Binary Hexafluorides to Explore New Chemical Signatures

Conclusions • Successfully halted the UF6 hydrolysis reaction and trapped intermediate complexes • Simple alcohol instead of water leads to shifting of bands in UF6 likely from change in matrix effects • MoF6 reaction shows some potential for methyl group being present in structure of intermediates • Future work with bulkier R-groups will hopefully reveal further information as to the reaction pathway and intermediate structure via kinetics and peak shifting

Dorris, Austin L. [Savannah River National Laborat↗
DOE OSTI2022

A Note on the Effects of Linear Topology Preservation in Monte Carlo Simulations of Knotted Proteins

Monte Carlo simulations are a powerful technique and are widely used in different fields. When applied to complex molecular systems with long chains, such as those in synthetic polymers and proteins, they have the advantage of providing a fast and computationally efficient way to sample equilibrium ensembles and calculate thermodynamic and structural properties under desired conditions. Conformational Monte Carlo techniques employ a move set to perform the transitions in the simulation Markov chain. While accepted conformations must preserve the sequential bonding of the protein chain model and excluded volume among its units, the moves themselves may take the chain across itself. We call this a break in linear topology preservation. In this manuscript, we show, using simple protein models, that there is no difference in equilibrium properties calculated with a move set that preserves linear topology and one that does not. However, for complex structures, such as those of deeply knotted proteins, the preservation of linear topology provides correct equilibrium results but only after long relaxation. In any case, to analyze folding pathways, knotting mechanisms and folding kinetics, the preservation of linear topology may be an unavoidable requirement.

Especial, João N. C. (ORCID:0000000254765370)↗
NASA NTRS1991

Estimation of the initial equilibrium constants in the formation of tetragonal lysozyme nuclei

Results are presented from a study of the equilibria, kinetic rates, and the aggregation pathway which leads from a lysozyme monomer crystal to a tetragonal crystal, using dialyzed and recrystallized commercial hen eggwhite lysozyme. Relative light scattering intensity measurements were used to estimate the initial equilibrium constants for undersaturated lysozyme solutions in the tetragonal regime. The K1 value was estimated to be (1-3) x 10 exp 4 L/mol. Estimates of subsequent equilibrium constants depend on the crystal aggregation model chosen or determined. Experimental data suggest that tetragonal lysozyme crystal grows by addition of aggregates preformed in the bulk solution, rather than by monomer addition.

Pusey, Marc L.↗
NASA NTRS2003

Sulfur Oxidation and Contrail Precursor Chemistry

Sulfuric acid (H2SO4), formed in commercial aircraft operations via fuel-S (goes to) SO2 (goes to) SO3 (goes to) H2SO4 plays an important role in the formation of contrails. It is believed that the first step occurs inside the combustor, the second step in the engine exit nozzle, and the third step in the exhaust plume. Thus, measurements of the sulfur oxidation rates are critical to the understanding of contrail formation. Field measurements of contrails formed behind commercial aircraft indicate that significantly greater conversion of fuel-bound sulfur to sulfate aerosol occurs than can be explained by our current knowledge of contrail physics and chemistry. The conversion of sulfur from S(IV) to S(VI) oxidation state, required for sulfate aerosol formation, is thermodynamically favored for the conditions that exist within jet engines but is kinetically disfavored. The principal reaction pathway is O+SO2+M (goes to) SO3+M. The rates of this reaction have never been measured in the temperature and pressure regimes available to aircraft operation. In the first year (FY02) of this project, we performed a series of experiments to elucidate the rate information for the O+SO2+M (goes to) SO3+M reaction. The work performed is described following the proposed work plan. Because we used the H2/O2 system for an O-atom source and rate coefficients were obtained via computer simulation, construction of a reaction mechanism and either recalculation or estimation of thermodynamic properties of H(x)SO(y) species are described first.

DeWitt, Kenneth J.↗
DOE OSTI2022

Vertex protein PduN tunes encapsulated pathway performance by dictating bacterial metabolosome morphology

Engineering subcellular organization in microbes shows great promise in addressing bottlenecks in metabolic engineering efforts; however, rules guiding selection of an organization strategy or platform are lacking. Here, we study compartment morphology as a factor in mediating encapsulated pathway performance. Using the 1,2-propanediol utilization microcompartment (Pdu MCP) system from Salmonella enterica serovar Typhimurium LT2, we find that we can shift the morphology of this protein nanoreactor from polyhedral to tubular by removing vertex protein PduN. Analysis of the metabolic function between these Pdu microtubes (MTs) shows that they provide a diffusional barrier capable of shielding the cytosol from a toxic pathway intermediate, similar to native MCPs. However, kinetic modeling suggests that the different surface area to volume ratios of MCP and MT structures alters encapsulated pathway performance. Finally, we report a microscopy-based assay that permits rapid assessment of Pdu MT formation to enable future engineering efforts on these structures.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2020

CRNT4SBML: a Python package for the detection of bistability in biochemical reaction networks

Motivation: Signaling pathways capable of switching between two states are ubiquitous within living organisms. They provide the cells with the means to produce reversible or irreversible decisions. Switchlike behavior of biological systems is realized through biochemical reaction networks capable of having two or more distinct steady states which are dependent on initial conditions. Investigation of whether a certain signaling pathway can confer bistability involves a substantial amount of hypothesis testing. The cost of direct experimental testing can be prohibitive. Therefore, constraining the hypothesis space is highly bene?cial. One such methodology is based on Chemical Reaction Network Theory which uses computational techniques to rule out pathways that are not capable of bistability regardless of kinetic constant values and molecule concentrations. Although useful, these methods are complicated from both pureandcomputationalmathematicsperspectives.Thus,theiradoptionisverylimitedamongstbiologists. Results: We brought Chemical Reaction Network Theory approaches closer to experimental biologists by automating all the necessary steps in CRNT4SMBL. The input is based on SBML format which is the community standard for biological pathway communication. The tool parses SBML and derives C-graph representations of the biological pathway with mass action kinetics. Next steps involve an ef?cient search for potential saddle-node bifurcation points using an optimization technique. This type of bifurcation is important as it has the potential of acting as a switching point between two steady states. Finally, if any bifurcation points are present, numerical continuation analysis extends the equilibria branches for generating the diagram. Presence of an S-shaped bifurcation diagram indicates that the pathway acts as a bistable switch for the given optimization parameters. Availability: CRNT4SBML is available via the Python Package Index. The documentation can be found at https://crnt4sbml.readthedocs.io. CRNT4SBML is licensed under the Apache Software License 2.0. Contact: vladislav.petyuk@pnnl.gov

Reyes, Brandon C.↗