The goal of this paper is to gather and develop some necessary and sufficient criteria for injectivity and multistationarity in vector fields associated with a chemical reaction network under a variety of more or less general assumptions on the natur
e of the network and the reaction rates. The results are primarily linear algebraic or matrix-theoretic, with some graph-theoretic results also mentioned. Several results appear in, or are close to, results in the literature. Here, we emphasise the connections between the results, and where possible, present elementary proofs which rely solely on basic linear algebra and calculus. A number of examples are provided to illustrate the variety of subtly different conclusions which can be reached via different computations. In addition, many of the computations are implemented in a web-based open source platform, allowing the reader to test examples including and beyond those analysed in the paper.
We describe combinatorial approaches to the question of whether families of real matrices admit pairs of nonreal eigenvalues passing through the imaginary axis. When the matrices arise as Jacobian matrices in the study of dynamical systems, these con
ditions provide necessary conditions for Hopf bifurcations to occur in parameterised families of such systems. The techniques depend on the spectral properties of additive compound matrices: in particular, we associate with a product of matrices a signed, labelled digraph termed a DSR^[2] graph, which encodes information about the second additive compound of this product. A condition on the cycle structure of this digraph is shown to rule out the possibility of nonreal eigenvalues with positive real part. The techniques developed are applied to systems of interacting elements termed interaction networks, of which networks of chemical reactions are a special case.
In this paper we discuss the question of how to decide when a general chemical reaction system is incapable of admitting multiple equilibria, regardless of parameter values such as reaction rate constants, and regardless of the type of chemical kinet
ics, such as mass-action kinetics, Michaelis-Menten kinetics, etc. Our results relate previously described linear algebraic and graph-theoretic conditions for injectivity of chemical reaction systems. After developing a translation between the two formalisms, we show that a graph-theoretic test developed earlier in the context of systems with mass action kinetics, can be applied to reaction systems with arbitrary kinetics. The test, which is easy to implement algorithmically, and can often be decided without the need for any computation, rules out the possibility of multiple equilibria for the systems in question.
We extend previous work on injectivity in chemical reaction networks to general interaction networks. Matrix- and graph-theoretic conditions for injectivity of these systems are presented. A particular signed, directed, labelled, bipartite multigraph
, termed the ``DSR graph, is shown to be a useful representation of an interaction network when discussing questions of injectivity. A graph-theoretic condition, developed previously in the context of chemical reaction networks, is shown to be sufficient to guarantee injectivity for a large class of systems. The graph-theoretic condition is simple to state and often easy to check. Examples are presented to illustrate the wide applicability of the theory developed.
