Abstract

Complex chemical reaction networks, which underlie many industrial and biological processes, often exhibit non-monotonic changes in chemical species concentrations, typically described using nonlinear models. Such non-monotonic dynamics are in principle possible even in linear models if the matrices defining the models are non-normal, as characterized by a necessarily non-orthogonal set of eigenvectors. However, the extent to which non-normality is responsible for non-monotonic behavior remains an open question. Here, using a master equation to model the reaction dynamics, we derive a general condition for observing non-monotonic dynamics of individual species, establishing that non-normality promotes non-monotonicity but is not a requirement for it. In contrast, we show that non-normality is a requirement for non-monotonic dynamics to be observed in the R\'enyi entropy. Using hydrogen combustion as an example application, we demonstrate that non-monotonic dynamics under experimental conditions are supported by a linear chain of connected components, in contrast with the dominance of a single giant component observed in typical random reaction networks. The exact linearity of the master equation enables development of rigorous theory and simulations for dynamical networks of unprecedented size (approaching $10^5$ dynamical variables, even for a network of only 20 reactions and involving less than 100 atoms). Our conclusions are expected to hold for other combustion processes, and the general theory we develop is applicable to all chemical reaction networks, including biological ones.

Highlights

  • Matrix non-normality is perhaps best known for its role in a counterintuitive form of nonlinear instability [1]

  • How prevalent is non-normality in dynamical networks and how often does it lead to non-monotonic dynamics? While non-normality is known to be widespread among matrices encoding network structures [12], the question is open for matrices representing dynamical interactions, which have direct implications on non-monotonic dynamics

  • This equation establishes a connection between information theory and dynamical systems theory: the rate of change of the Rényi entropy is directly proportional to the Rayleigh quotient of S, which is of interest in non-normal growth analysis [2,3]

Read more

Summary

MASTER-EQUATION FORMULATION

Given a set of Ns chemical species and a set of Nr reactions involving them, we consider the dynamics of reactions in a mixture of these species. From a given initial state of this H2/O2 combustion process with specified numbers of molecules of each species, the number of states that are accessible through a sequence of reactions is finite, and we denote this number by N This is because the number of each atomic species in the gas is conserved during each reaction event. Since the transition rates Wi j span multiple orders of magnitude, the flow of probability under Eq (1) is typically limited to a small portion of the network This observation can be used up front to reduce the size of the network, and the computational burden, without significantly affecting the dynamics.

SPECTRAL ANALYSIS
NON-MONOTONIC ENTROPY DYNAMICS
CONDITION FOR NON-MONOTONIC
CONNECTED COMPONENT ANALYSIS
RANDOM NETWORKS
Findings
CONCLUSIONS
Full Text
Published version (Free)

Talk to us

Join us for a 30 min session where you can share your feedback and ask us any queries you have

Schedule a call