Non-normality and non-monotonic dynamics in complex reaction networks

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ényi 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.

arXiv commentSoftware implementing our methods is available as a Github repository at https://github.com/znicolaou/ratematrix and an animated version of Fig. 1 is available at https://northwestern.box.com/s/otn3m2cov9gi3enht3r8jh5kjo9qnv6d

References in corpus (15)

• Statistical mechanics of complex networks
• Epidemic spreading in scale-free networks
• Random graphs with arbitrary degree distributions and their applications
• Equilibrium free energy differences from nonequilibrium measurements: a master equation approach
• Giant strongly connected component of directed networks
• Synchronization is optimal in non-diagonalizable networks
• Maximum Performance at Minimum Cost in Network Synchronization
• The limitations of model-based experimental design and parameter estimation in sloppy systems
• Nonnormal amplification in random balanced neuronal networks
• Generalized Master Equations for Non-Poisson Dynamics on Networks
• Robustness of Optimal Synchronization in Real Networks
• Giant amplification of noise in fluctuation-induced pattern formation
• Topological resilience in non-normal networked systems
• Universal transient behavior in large dynamical systems on networks
• Resilience for stochastic systems interacting via a quasi-degenerate network

Cited by in corpus (8)

• Non-reciprocal phase transitions
• Network structural origin of instabilities in large complex systems
• Synchronization dynamics in non-normal networks: the trade-off for optimality
• Universal dynamic scaling in chemical reactions at and away from equilibrium
• Prevalence of multistability and nonstationarity in driven chemical networks
• Revisiting weak values through non-normality
• What is a chemostat? Insights from hybrid dynamics and stochastic thermodynamics
• Nonmodal amplitude equations