Cores in parameter-rich Reaction Networks

Parameter-rich Chemical Reaction Networks are built upon kinetic laws that have sufficient flexibility to allow the independent choice of fixed-point coordinates and the absolute numerical values of the non-zero entries of the Jacobian. Although this seems to be a tall order, Michaelis-Menten and generalized mass action kinetics are of this type. In this setting, the existence of choices of parameters that make a given inner fixed point Hurwitz unstable is determined entirely by the stoichiometric matrix. More precisely, the existence of certain quadratic sub-matrices is sufficient. Minimal submatrices of this type form unstable cores. Autocatalytic cores (in the sense of Nghe) turn out to be unstable cores that, in addition, have the special structure of Metzler matrices. Autocatalysis, therefore can always act as source of instability given suitable choices of parameters, but there are also other causes of dynamic instability, entirely unrelated to autocatalysis. The notion of cores is surprisingly versatile and can be extended to other forms of dynamical behavior that is determined by spectral properties. In particular, it pertains to the search for oscillatory behavior by virtue of guaranteeing the existence of Hopf bifurcation.

(Joint work with Nicola Vassena and Alex Blockhuis)