Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userRobust Lyapunov Functions for Reaction Networks: An Uncertain System Framework
65
0
0.0
(
0
)
Ask ChatGPT about the research
No Arabic abstract
We present a framework to transform the problem of finding a Lyapunov function of a Chemical Reaction Network (CRN) in concentration coordinates with arbitrary monotone kinetics into finding a common Lyapunov function for a linear parameter varying system in reaction coordinates. Alternative formulations of the proposed Lyapunov function is presented also. This is applied to reinterpret previous results by the authors on Piecewise Linear in Rates Lyapunov functions, and to establish a link with contraction analysis. Persistence and uniqueness of equilibria are discussed also.
rate research
Read More
Piecewise-Linear in Rates (PWLR) Lyapunov functions are introduced for a class of Chemical Reaction Networks (CRNs). In addition to their simple structure, these functions are robust with respect to arbitrary monotone reaction rates, of which mass-action is a special case. The existence of such functions ensures the convergence of trajectories towards equilibria, and guarantee their asymptotic stability with respect to the corresponding stoichiometric compatibility class. We give the definition of these Lyapunov functions, prove their basic properties, and provide algorithms for constructing them. Examples are provided, relationship with consensus dynamics are discussed, and future directions are elaborated.
We consider the relationship between stationary distributions for stochastic models of reaction systems and Lyapunov functions for their deterministic counterparts. Specifically, we derive the well known Lyapunov function of reaction network theory as a scaling limit of the non-equilibrium potential of the stationary distribution of stochastically modeled complex balanced systems. We extend this result to general birth-death models and demonstrate via example that similar scaling limits can yield Lyapunov functions even for models that are not complex or detailed balanced, and may even have multiple equilibria.
In 1961, Renyi discovered a rich family of non-classical Lyapunov functions for kinetics of the Markov chains, or, what is the same, for the linear kinetic equations. This family was parameterised by convex functions on the positive semi-axis. After works of Csiszar and Morimoto, these functions became widely known as $f$-divergences or the Csiszar--Morimoto divergences. These Lyapunov functions are universal in the following sense: they depend only on the state of equilibrium, not on the kinetic parameters themselves. Despite many years of research, no such wide family of universal Lyapunov functions has been found for nonlinear reaction networks. For general non-linear networks with detailed or complex balance, the classical thermodynamics potentials remain the only universal Lyapunov functions. We constructed a rich family of new universal Lyapunov functions for {em any non-linear reaction network} with detailed or complex balance. These functions are parameterised by compact subsets of the projective space. They are universal in the same sense: they depend only on the state of equilibrium and on the network structure, but not on the kinetic parameters themselves. The main elements and operations in the construction of the new Lyapunov functions are partial equilibria of reactions and convex envelopes of families of functions.
Homochirality, i.e. the dominance across all living matter of one enantiomer over the other among chiral molecules, is thought to be a key step in the emergence of life. Building on ideas put forward by Frank and many others, we proposed recently one such mechanism in G. Laurent et al., PNAS (2021) based on the properties of large out of equilibrium chemical networks. We showed that in such networks, a phase transition towards an homochiral state is likely to occur as the number of chiral species in the system becomes large or as the amount of free energy injected into the system increases. This paper aims at clarifying some important points in that scenario, not covered by our previous work. We first analyze the various conventions used to measure chirality, introduce the notion of chiral symmetry of a network, and study its implications regarding the relative chiral signs adopted by different groups of molecules. We then propose a generalization of Franks model for large chemical networks, which we characterize completely using methods of random matrices. This analysis can be extended to sparse networks, which shows that the emergence of homochirality is a robust transition.
We introduce the notion of non-oscillation, propose a constructive method for its robust verification, and study its application to biological interaction networks (also known as, chemical reaction networks). We begin by revisiting Muldowneys result on non-existence of periodic solutions based on the study of the variational system of the second additive compound of the Jacobian of a nonlinear system. We show that exponential stability of the latter rules out limit cycles, quasi-periodic solutions, and broad classes of oscillatory behavior. We focus then on nonlinear equations arising in biological interaction networks with general kinetics, and we show that the dynamics of the aforementioned variational system can be embedded in a linear differential inclusion. We then propose algorithms for constructing piecewise linear Lyapunov functions to certify global robust non-oscillatory behavior. Finally, we apply our techniques to study several regulated enzymatic cycles where available methods are not able to provide any information about their qualitative global behavior.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
