The complete set of steady states for rate constants in the range 0?ki<∞ and concentrations in the range 0?Xi<∞ is given explicitly in parametric form for the general chemical reaction system. The only assumptions are that the stoichiometries are real numbers, and the reaction rates are proportional to functions of class Ck, k≳0; the functions are assumed to be positive in the interior of the domain. Hence, these results apply far more generally than just to chemical systems and should be valid for many ecological and economic models as well. The set of steady states in the interior of the domain is in general a simply connected differentiable manifold M of dimension n+r−d, where n = number of species, r = number of reactions, and d = rank ν (ν = stoichiometric matrix). The full set of steady states M* consists of those in the interior (M) plus a frequently very complicated set of steady states lying in the boundary. M* is the union of a set of differentiable manifolds but is not itself a differentiable manifold in general. If the reaction rate laws are all monomials (i.e., products of concentrations raised to some ’’order of kinetics’’), a working APL algorithm is given which provides enough information to construct explicitly the set of manifolds comprising M*. This algorithm is used to treat an early version of the Belousov–Zhabotinski reaction mechanism. The algorithm may be used for more general rate laws with minor modifications.
Read full abstract