Search NASA⌕ Search

SEARCH · Search NASA

Results for “reaction rate equations”

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 91 records · Page 5
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Modeling delayed thermal runaway in nitric acid-soaked cat litter mixed with radioactive waste

Thermal ignition of radioactive waste within a 55-gallon drum was simulated by using a pressure-dependent waste decomposition model (Hobbs et al. in Process Saf Environ Prot https://doi.org/10.1016/j.psep.2022.09.047, 2022) calibrated with data from full-scale drum experiments (Parker et al. in The thermolytic response of a surrogate RNS waste mixture at the drum scale. Los Alamos National Laboratory Report LA-UR-16-21760, 2016) and validated with experiments from multiple laboratories (Hobbs et al. in Thermal analysis of aged nitric acid-soaked kitty litter in TRU waste drums-23370.WM2023 Conference, Phoenix, AZ, 2023). The acceleration of nitric acid chemistry reacting with an organic cat litter leading to thermal ignition was likely triggered by a restricted vent in the drum. Here, we address whether the form of the rate equation in (Hobbs et al. in Process Saf Environ Prot https://doi.org/10.1016/j.psep.2022.09.047, 2022) is sufficient to extrapolate thermal ignition within aged drums of similar content that have been stored in Texas for over nine years by investigating four different reaction rate forms for waste decomposition. A critical reaction rate reduction analysis is performed on each of these models to determine if delayed thermal runaway within vented aged waste is possible after nine years. We found that a pressure-dependent first-order rate expression not only predicted the accidental ignition of the waste drum, but the form also matches multiple experiments from different laboratories. Even though the waste composition decreases over time, the model predicts that acceleration leading to thermal runaway is possible if the waste is confined, even after 9 years. In conclusion, waste containing oxidizers such as nitric acid should not be mixed with organic adsorbents, especially if the waste is confined.

cookoff↗
DOE OSTI2023

Theoretical Kinetics Predictions for Reactions on the NH 2 O Potential Energy Surface

Recent modeling studies of ammonia oxidation, which are motivated by the prospective role of ammonia as a zero-carbon fuel, have indicated significant discrepancies among the existing literature mechanisms. In this study, high-level theoretical kinetics predictions have been obtained for reactions on the NH 2 O potential energy surface, including the NH 2 + O, HNO + H, and NH + OH reactions. These reactions have previously been highlighted as important reactions in NH 3 oxidation with high sensitivity and high uncertainty. The potential energy surface is explored with coupled cluster calculations, including large basis sets and high-level corrections to yield high-accuracy (~0.2 kcal/mol 2σ uncertainty) estimates of the stationary point energies. Variational transition state theory is used to predict the microcanonical rate constants, which are then incorporated in master equation treatments of the temperature- and pressure-dependent kinetics. For radical–radical channels, the microcanonical rates are obtained from variable reaction coordinate transition state theory implementing directly evaluated multireference electronic energies. The analysis yields predictions for the total rate constants as well as the branching ratios. We find that the NO + H 2 channel contributes 10% of the total NH 2 + O flux at combustion temperatures, although this channel is not included in modern NH 3 oxidation mechanisms. Modeling is used to illustrate the ramifications of these rate predictions on the kinetics of NH 3 oxidation and NO x formation. Finally, the present results for NH 2 + O are important for predicting the chain branching and formation of NO in the oxidation of NH 3 and thermal DeNO x .

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

Theoretical kinetics predictions for NH 2 + HO 2

Recent modeling studies of NH 3 oxidation, which are motivated by the prospective role of ammonia as a zero-carbon fuel, have indicated significant discrepancies between existing literature mechanisms. In this study high level theoretical kinetics predictions have been obtained for the reaction of NH 2 with HO 2 , which has previously been highlighted as an important reaction with high sensitivity and high uncertainty. The potential energy surface is explored with coupled cluster calculations including large basis sets and high-level corrections to yield high accuracy (~0.2 kcal/mol) estimates of the stationary point energies. Variational transition state theory is used to predict the microcanonical rate constants, which are then incorporated in master equation treatments of the temperature and pressure dependent kinetics. For the radical-radical channels, the microcanonical rates are obtained from variable reaction coordinate transition state theory implementing directly evaluated multireference electronic energies. The analysis yields predictions for the total rate constant as well as the branching to the NH 3 + O 2 , H 2 NO + OH, and HNO + H 2 O channels. Rate constants are also reported for the H 2 NO + OH reaction as they arise naturally from the analysis. The rate constant and branching fraction determined in this work for the NH 2 + HO 2 reaction deviate significantly from values used in most previous modeling studies. The fact that the main product channel is chain terminating, rather than propagating, has strong implications for modeling NH 3 ignition and oxidation, in particular at intermediate temperatures and elevated pressure.

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

Mechanism, Thermochemistry, and Kinetics for the CH + N 2 Reaction Leading to Prompt NO Formation in Combustion

Here, the reaction of CH + N 2 forming H + NCN is a remarkable example of activation of the nitrogen triple bond and is an important source of prompt NO in combustion. The reaction pathway is complex and proceeds through two competing mechanisms: a cyclic addition channel initiated by c-HC(NN) and a chain-addition channel initiated by HCNN, both of which eventually form HNCN prior to dissociation to H + NCN. This work reinvestigates this reaction with composite coupled cluster protocols, including a novel spin-flip equation of motion coupled cluster scheme, combined with pragmatic two-dimensional master equation simulations of the resulting rate coefficients. These improved calculations predict the CH + N 2 rate coefficient between the two more recent previous theoretical results and reduce the uncertainties of the best theoretical models of this reaction to less than a factor of 1.3. Additionally, we provide a closer theoretical investigation of the simultaneous dependence of the CH + N 2 rate coefficient on pressure and temperature, and affirm that collisionally stabilized HNCN, another potential source of prompt NO, emerges as an appreciable product of this reaction under conditions relevant to automotive internal combustion engines and aircraft gas turbine engines.

Nguyen, Thanh Lam [Univ. of Florida, Gainesville, ↗
DOE OSTI2023

Mass measurement of 27 P to constrain type-I x-ray burst models and validate the isobaric multiplet mass equation for the $A = 27, T = \frac{3}{2}$ isospin quartet

Background: Light curves are the primary observable of type-I x-ray bursts. Computational x-ray burst models must match simulations to observed light curves. Most of the error in simulated curves comes from uncertainties in rp process reaction rates, which can be reduced via precision mass measurements of neutron-deficient isotopes in the rp process path. Purpose: Perform a precise atomic mass measurement of 27 P. Use this new measurement to calculate rp process reaction rates and input these rates into an x-ray burst model to reduce simulated light curve uncertainty. Use the mass measurement of 27 P to validate the isobaric multiplet mass equation (IMME) for the A = 27 T = 3/2 isospin quartet which 27 P belongs to. Method: High-precision Penning trap mass spectrometry utilizing the time-of-flight ion cyclotron resonance technique was used to determine the atomic mass of 27 P. The MESA code (Modules for Experiments in Stellar Astrophysics) was then used to simulate x-ray bursts using a one-dimensional multizone model to produce updated light curves. Results: The mass excess of 27 P was measured to be –670.7(6) keV, a 14-fold precision increase over the mass reported in the 2020 Atomic Mass Evaluation (AME2020). The 26 Si(p, $γ$) 27 P– 27 P($γ$, p) 26 Si rate equilibrium has been determined to a higher precision based on the precision mass measurement of 27 P. x-ray burst light curves were produced with the MESA code using the new reaction rates. Changes in the mass of 27 P seem to have minimal effect on light curves, even in burster systems tailored to maximize impact. Conclusion: The mass of 27 P does not play a significant role in x-ray burst light curves. It is important to understand that more advanced models do not just provide more precise results, but often qualitatively different ones. This result brings us a step closer to being able to extract stellar parameters from individual x-ray burst observations. In addition, the IMME has been validated for the A = 27, T = 3/2 quartet. The normal quadratic form of the IMME using the latest data yields a reduced $χ$ 2 of 2.9. Finally, the cubic term required to generate an exact fit to the latest data matches theoretical attempts to predict this term.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

Inferring small-sample-reactivity worth from oscillation experiment with auto-corrected local flux

To estimate reactivity changes as a function of time from the signal of detectors count rates, one usually uses point kinetic equations. However, when detectors are located inside the core, it is often necessary to compute a correction factor to account for the flux-shape deformation (i.e. the change of detector efficiency that may have been induced by neutron absorbers movements). Various simulation methods have been developed depending on transients observed. They range from the resolution of static equations up to full-time-dependent neutron transport in 3D. However, those methods require high computational-cost calculations of local reaction rates that rely on nuclear data evaluations. Therefore, a 'measured' reactivity is never determined through pure measurements. Taking advantage of the known oscillation movement of a small sample at the center of the MINERVE reactor, the innovative algorithm presented in this paper derives the correction factor needed to get no time drift on reactivity steps. Derived correction factors are in good agreement with those obtained independently with TRIPOLI-4 through sensitivity analysis. Low sensitivity to the sample composition and moderate sensitivity to the detectors position reinforce the reliability of those results. This application paves the way to almost pure nuclear data reactivity measurements. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2024

Theory Manual for the Reaction Summary Module- SUMMAR

The post-processing module of DIF3D called SUMMAR was created to produce commonly used reactor engineering outputs. It is a optional output option for DIF3D and provides alternative summarizes of power, flux, and neutron balance over the user specified regions and areas. It can also provide isotope-wise microscopic and macroscopic reaction rates along with the effective one-group microscopic cross sections needed for solving the Bateman equations. This report serves as both the theory manual, user guide, and software verification of the SUMMAR module. The methodology and equations for all SUMMAR calculated parameters are presented here and discussed. Then the verification work, based upon simple DIF3D models, is presented to validate the derived equations. The SUMMAR inputs and outputs are detailed and examples are given to demonstrate it.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2021

Substitution Reactions in the Pyrolysis of Acetone Revealed through a Modeling, Experiment, Theory Paradigm

The development of high-fidelity mechanisms for chemically reactive systems is a challenging process that requires the compilation of rate descriptions for a large and somewhat ill-defined set of reactions. The present unified combination of modeling, experiment, and theory provides a paradigm for improving such mechanism development efforts. Here we combine broadband rotational spectroscopy with detailed chemical modeling based on rate constants obtained from automated ab initio transition state theory-based master equation calculations and high-level thermochemical parametrizations. Broadband rotational spectroscopy offers quantitative and isomer-specific detection by which branching ratios of polar reaction products may be obtained. Using this technique, we observe and characterize products arising from H atom substitution reactions in the flash pyrolysis of acetone (CH 3 C(O)CH 3 ) at a nominal temperature of 1800 K. The major product observed is ketene (CH 2 CO). Minor products identified include acetaldehyde (CH 3 CHO), propyne (CH 3 CCH), propene (CH 2 CHCH 3 ), and water (HDO). Literature mechanisms for the pyrolysis of acetone do not adequately describe the minor products. The inclusion of a variety of substitution reactions, with rate constants and thermochemistry obtained from automated ab initio kinetics predictions and Active Thermochemical Tables analyses, demonstrates an important role for such processes. Furthermore, the pathway to acetaldehyde is shown to be a direct result of substitution of acetone's methyl group by a free H atom, while propene formation arises from OH substitution in the enol form of acetone by a free H atom.

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

Fast estimation of reaction rates in spherical and non-spherical porous catalysts

Here we present a methodology for modeling multi-step reaction rates in porous catalyst particles for use in CFD-DEM and two fluid models. Single-step effectiveness factors based on a Thiele modulus, while useful, cannot accurately capture the cascading reaction systems common in high temperature vapor-phase chemical reactors like fluidized catalytic cracking units and catalytic biomass fast pyrolysis systems. Instead, multi-step effectiveness vectors derived from steady-state solutions to the governing reaction-diffusion equations are needed. Solutions for various catalyst shapes are presented, including spheres, cylinders, and prisms. Computational challenges inherent in repeated evaluation of reaction rates with diffusion limitations are discussed, and an efficient implementation based on pre-computed lookup tables is proposed and demonstrated on a simulation of a fluidized bed reactor. Open-source code is provided for the compilation of reaction rate tables for use in ODE, DEM, and two-fluid models.

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

Polysulfone shock compressed above the decomposition threshold: Velocimetry and modeling of two-wave structures

Polysulfone was shock loaded to pressures of 14.6–26.2 GPa in a series of gas gun-driven plate-impact experiments measuring material response with embedded electro-magnetic particle velocity gauges and optical velocimetry. The embedded electro-magnetic particle velocity gauges did not show a distinct two-wave structure but did show rounding that suggested a reaction but not a distinct separation of the reactants and product waves. In contrast, the transmission experiments fielded with optical velocimetry, with product pressures ranging from 21.0 to 26.2 GPa, showed well-defined two-wave structures due to shock-driven chemical decomposition of the polymer to products at a higher density. Distinct two-wave structures have not previously been published in literature during polymer compression; here, we observed these two-wave structures at both the polymer/lithium fluoride and polymer/polymethylmethacrylate interfaces. Hydrodynamic simulations of the experiments were performed using a newly constructed SESAME equation of state (EOS) for the polymer, a thermochemical EOS for the decomposition products, and an Arrhenius reaction rate model for the kinetics of the reaction. Finally, simulation results also demonstrated two-wave structures but were unable to quantitatively reproduce either the embedded gauge or optical velocimetry data.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

AWSD reactive flow model for PBX 9404

An Arrhenius–Wescott–Stewart–Davis (AWSD) reactive flow model for high explosive PBX 9404 is developed. We specifically calibrate an AWSD model for PBX 9404 by fitting equations of state for reactants and detonation products to the results of thermochemical calculations and to experimental data from multiple sources. The calibrated equations of state are then coupled with an Arrhenius rate law based on shock temperature that describes the reaction progress during PBX 9404 detonation. The parameter values in the rate law are calibrated to experimental gas-gun data and diameter effect results. The results of the calibrated AWSD model are in strong agreement with available experimental data for PBX 9404. A similar level of agreement between predicted and experimental results is observed when the calibrated model is validated on data that were not used in the model parameterization procedure. Our results illustrate that the AWSD model is capable of accurately describing the many important properties and observables in the reactive burn of PBX 9404. Because of the historical significance of PBX 9404 in high explosives research and its current use in aging studies, this work provides an important model of a legacy material, which can be used to make comparisons to new high explosive formulations.

45 MILITARY TECHNOLOGY, WEAPONRY, AND NATIONAL DEF↗
DOE OSTI2022

Dissociation-induced depletion of high-energy reactant molecules as a mechanism for pressure-dependent rate constants for bimolecular reactions

In 1922, Lindemann proposed the now-well-known mechanism for pressure-dependent rate constants for unimolecular reactions: reactant molecules with sufficiently high energies dissociate more quickly than collisions can reestablish the Boltzmann distribution of the internal energies of the molecule during its dissociation at low pressures – yielding pressure-dependent rate constants for unimolecular reactions due to the preferential depletion of the high energy states capable of dissociation. In the last century, incredible progress has been made in achieving a far greater understanding of and quantitative predictions for unimolecular and association reactions. In the modern era, pressure-dependent phenomenological rate constants are now nearly universally used to describe the rates of unimolecular and associative reactions in phenomenological kinetic modeling. However, there is a second, more indirect, implication of Lindemann’s mechanism that relates to how these dissociation-induced non-equilibrium distributions impact bimolecular reactions, including non-associative bimolecular reactions – which are generally not considered to have pressure-dependent rate constants. Yet, as we show herein, the same high energy states depleted due to dissociation would otherwise react most rapidly in high-activation-energy bimolecular reactions – yielding a mechanism for pressure-dependent rate constants for bimolecular reactions (including non-associative reactions). Here, we present results from a case study for CH 2 O dissociation, isomerization, and bimolecular reaction with O2 to explore this question. Results from our master equation calculations indicate that the effect of dissociation-induced non-equilibrium distributions on bimolecular reactions can be substantial – even when chemical timescales are well separated from internal energy relaxational timescales (i.e. when the traditional rate constant description would be thought to apply). This effect is found to be more pronounced – and more complex – for bimolecular reactions involving molecular entities whose chemical timescales are merged with the internal energy relaxational timescales. Lastly, we present some ideas for discussion regarding what should be considered as “chemical species” in phenomenological kinetic models.

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

On the prediction of pressure effects for the combination kinetics of two alkyl radicals with the geometric mean rule

Here, alkyl combination reactions, along with their reverse alkane unimolecular decomposition reactions, are widely present in pyrolysis and combustion systems. These reactions occur mainly at temperatures at which their rate coefficients depend on pressure, and only a few studies have explored this dependence for large size alkyls. The accurate calculation of these kinetic data is limited by the size of the alkyl fragments involved, since the theoretical methods used are computationally expensive. Therefore, little reliable data is available for the modeling of real fuels alkane decomposition kinetics. In this paper, we propose an alternative approach that allows an accurate estimation of pressure-dependent kinetic rates for the combination of long-chain alkyl radicals, based on tabulated reaction models. The rates of combination of linear alkyls radicals forming ethane to pentane were calculated using the CASPT2-F12/cc-pVDZ-F12 level within VRC-TST theory combined with master-equation simulations. The computed high-pressure limit rate coefficients showed that n-propyl could be a minimum-size representative structure for larger n-alkyl chain, but this approach fails to capture the pressure dependence of the reaction that remains size-dependent for larger alkyls. The Geometric Mean Rule (GMR), a simple rule that links the self and cross-reactions of three alkyl radicals, is used to address this issue. Our theoretical calculations show for the first time that the GMR applies not only at the high-pressure limit, but also for every pressure explored. The GMR is shown to be accurate, with a maximum 22 % deviation in the falloff region for the estimation of rates for small alkyl combinations. Using this method, we propose a tabulated reaction rate rule for the combination of alkyl radicals up to C5+C5. The use of computed rate coefficients with the GMR approximation opens ways to develop pressure-dependent reaction rate rules for large n-alkyl combinations, based on accurate and computationally expensive theoretical calculations on smaller alkyl fragments.

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

An exploratory steady–state redox model of photosynthetic linear electron transport for use in complete modelling of photosynthesis for broad applications

A photochemical model of photosynthetic electron transport (PET) is needed to integrate photophysics, photochemistry, and biochemistry to determine redox conditions of electron carriers and enzymes for plant stress assessment and mechanistically link sun-induced chlorophyll fluorescence to carbon assimilation for remotely sensing photosynthesis. Towards this goal, we derived photochemical equations governing the states and redox reactions of complexes and electron carriers along the PET chain. These equations allow the redox conditions of the mobile plastoquinone pool and the cytochrome b 6 f complex (Cyt) to be inferred with typical fluorometry. The equations agreed well with fluorometry measurements from diverse C 3 /C 4 species across environments in the relationship between the PET rate and fraction of open photosystem II reaction centres. We found the oxidation of plastoquinol by Cyt is the bottleneck of PET, and genetically improving the oxidation of plastoquinol by Cyt may enhance the efficiency of PET and photosynthesis across species. Redox reactions and photochemical and biochemical interactions are highly redundant in their complex controls of PET. Although individual reaction rate constants cannot be resolved, they appear in parameter groups which can be collectively inferred with fluorometry measurements for broad applications. The new photochemical model developed enables advances in different fronts of photosynthesis research.

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

Reaction–drift–diffusion models from master equations: application to material defects

We present a general method to produce well-conditioned continuum reaction–drift–diffusion equations directly from master equations on a discrete, periodic state space. We assume the underlying data to be kinetic Monte Carlo models (i.e. continuous-time Markov chains) produced from atomic sampling of point defects in locally periodic environments, such as perfect lattices, ordered surface structures or dislocation cores, possibly under the influence of a slowly varying external field. Our approach also applies to any discrete, periodic Markov chain. Here, the analysis identifies a previously omitted non-equilibrium drift term, present even in the absence of external forces, which can compete in magnitude with the reaction rates, thus being essential to correctly capture the kinetics. To remove fast modes which hinder time integration, we use a generalized Bloch relation to efficiently calculate the eigenspectrum of the master equation. A well conditioned continuum equation then emerges by searching for spectral gaps in the long wavelength limit, using an established kinetic clustering algorithm to define a proper reduced, Markovian state space.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Approach to Evaluating Reorganization Energies of Interfacial Electrochemical Reactions

Reaction rate coefficients for electron-transfer processes at the electrode–electrolyte interface are commonly estimated by using the Butler–Volmer equation, but their values are inaccurate beyond a few tenths of volts of overpotential. The Marcus–Hush–Chidsey (MHC) formalism yields correct asymptotic behavior of the rate coefficients vs applied overpotential but has complex dependencies on the redox system’s intrinsic parameters, which can be difficult to model or measure. In this work, we bridge the two kinetics formalisms to estimate the reorganization energy, one of the important parameters for the MHC formalism, and investigate its dependence on other intrinsic parameters such as activation barriers, electronic coupling strength, and the density of states of the electrode surface. We examine the sensitivity of the reorganization energy to these parameters, establish some general relationships for accurately predicting rate coefficients using the MHC formalism over a wide range of applied overpotentials, and compare this approach to calculating MHC rate constants with other empirical approaches for the mechanisms of CO 2 reduction on different metal electrode surfaces.

Butler−Volmer↗
DOE OSTI2021

Dynamic Catalysis Fundamentals: I. Fast calculation of limit cycles in dynamic catalysis

Dynamic catalysis—the forced oscillation of catalytic reaction coordinate potential energy surfaces (PES)—has recently emerged as a promising method for the acceleration of heterogeneously-catalyzed reactions. Theoretical study of enhancement of rates and supra-equilibrium product yield via dynamic catalysis has, to-date, been severely limited by onerous computational demands of forward integration of stiff, coupled ordinary differential equations (ODEs) that are necessary to quantitatively describe periodic cycling between PESs. Here, we establish a new approach that reduces, by ≳108×, the computational cost of finding the time-averaged rate at dynamic steady state (i.e. the limit cycle for linear and nonlinear systems of kinetic equations). Our developments are motivated by and conceived from physical and mathematical insight drawn from examination of a simple, didactic case study for which closed-form solutions of rate enhancement are derived in explicit terms of periods of oscillation and elementary step rate constants. Generalization of such closed-form solutions to more complex catalytic systems is achieved by introducing a periodic boundary condition requiring the dynamic steady state solution to have the same periodicity as the kinetic oscillations and solving the corresponding differential equations by linear algebra or Newton-Raphson-based approaches. The methodology is well-suited to extension to non-linear systems for which we detail the potential for multiple solutions or solutions with different periodicities. For linear and non-linear systems alike, the acute decrement in computational expense enables rapid optimization of oscillation waveforms and, consequently, accelerates understanding of the key catalyst properties that enable maximization of reaction rates, conversions, and selectivities during dynamic catalysis.

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

A numerical approach for detecting switch-like bistability in mass action chemical reaction networks with conservation laws

Theoretical analysis of signaling pathways can provide a substantial amount of insight into their function. One particular area of research considers signaling pathways capable of assuming two or more stable states given the same amount of signaling ligand. This phenomenon of bistability can give rise to switch-like behavior, a mechanism that governs cellular decision making. Investigation of whether or not a signaling pathway can confer bistability and switch-like behavior, without knowledge of specific kinetic rate constant values, is a mathematically challenging problem. Recently a technique based on optimization has been introduced, which is capable of finding example parameter values that confer switch-like behavior for a given pathway. Although this approach has made it possible to analyze moderately sized pathways, it is limited to reaction networks that presume a uniterminal structure. It is this limited structure we address by developing a general technique that applies to any mass action reaction network with conservation laws. In this paper we developed a generalized method for detecting switch-like bistable behavior in any mass action reaction network with conservation laws. The method involves (1) construction of a constrained optimization problem using the determinant of the Jacobian of the underlying rate equations, (2) minimization of the objective function to search for conditions resulting in a zero eigenvalue, (3) computation of a confidence level that describes if the global minimum has been found and (4) evaluation of optimization values, using either numerical continuation or directly simulating the ODE system, to verify that a bistability region exists. The generalized method has been tested on three motifs known to be capable of bistability. We have developed a variation of an optimization-based method for the discovery of bistability, which is not limited to uniterminal chemical reaction networks. Successful completion of the method provides an S-shaped bifurcation diagram, which indicates that the network acts as a bistable switch for the given optimization parameters.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗