BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Gheorghe Craciun

Gheorghe Craciun appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
math.DSMolecular Networksmath.PRDiscrete Mathematicsmath.AGQuantitative MethodsSubcellular Processes

Catalog footprint

What is connected

17works
7topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

17 published item(s)

preprint2026arXiv

A flux-based approach for analyzing the disguised toric locus of reaction networks

Dynamical systems with polynomial right-hand sides are very important in various applications, e.g., in biochemistry and population dynamics. The mathematical study of these dynamical systems is challenging due to the possibility of multistability, oscillations, and chaotic dynamics. One important tool for this study is the concept of reaction systems, which are dynamical systems generated by reaction networks for some choices of parameter values. Among these, disguised toric systems are remarkably stable: they have a unique attracting fixed point, and cannot give rise to oscillations or chaotic dynamics. The computation of the set of parameter values for which a network gives rise to disguised toric systems (i.e., the disguised toric locus of the network) is an important but difficult task. We introduce new ideas based on network fluxes for studying the disguised toric locus. We prove that the disguised toric locus of any network $G$ is a contractible manifold with boundary, and introduce an associated graph $G^{\max}$ that characterizes its interior. These theoretical tools allow us, for the first time, to compute the full disguised toric locus for many networks of interest.

0 citations0 reviewsNo code linkedOpen access
preprint2026arXiv

Extinction in Reaction Network Models

In this paper, we study extinction in dynamical systems generated by reaction networks. We introduce two notions: weak extinction and strong extinction, and relate them to the structure of the underlying network through Lyapunov functions and LaSalle's invariance principle. In particular, for all deficiency-zero networks that are not weakly reversible, we provide a geometric construction of linear Lyapunov functions. Using these functions, we establish that if these networks have bounded invariant subspaces, then they must exhibit weak extinction within every such subspace. Also, for linear networks that are not weakly reversible, we show that every species outside a terminal strongly connected component undergoes strong extinction. Moreover, in order to further emphasize the difference between weak and strong extinction, we construct an example of a reaction system (based on the Ivanova network) that exhibits weak extinction for all the species, but does not exhibit strong extinction in any species.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

An algorithm for finding weakly reversible deficiency zero realizations of polynomial dynamical systems

Systems of differential equations with polynomial right-hand sides are very common in applications. On the other hand, their mathematical analysis is very challenging in general, due to the possibility of complex dynamics: multiple basins of attraction, oscillations, and even chaotic dynamics. Even if we restrict our attention to mass-action systems, all of these complex dynamical behaviours are still possible. On the other hand, if a polynomial dynamical system has a weakly reversible deficiency zero ($WR_0$) realization, then its dynamics is known to be remarkably simple: oscillations and chaotic dynamics are ruled out and, up to linear conservation laws, there exists a single positive steady state, which is asymptotically stable. Here we describe an algorithm for finding $WR_0$ realizations of polynomial dynamical systems, whenever such realizations exist.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Disguised toric dynamical systems

We study families of polynomial dynamical systems inspired by biochemical reaction networks. We focus on complex balanced mass-action systems, which have also been called toric. They are known or conjectured to enjoy very strong dynamical properties, such as existence and uniqueness of positive steady states, local and global stability, persistence, and permanence. We consider the class of disguised toric dynamical systems, which contains toric dynamical systems, and to which all dynamical properties mentioned above extend naturally. By means of (real) algebraic geometry we show that some reaction networks have an empty toric locus or a toric locus of Lebesgue measure zero in parameter space, while their disguised toric locus is of positive measure. We also propose some algorithms one can use to detect the disguised toric locus.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Multistationarity in cyclic sequestration-transmutation networks

We consider a natural class of reaction networks which consist of reactions where either two species can inactivate each other (i.e., sequestration), or some species can be transformed into another (i.e., transmutation), in a way that gives rise to a feedback cycle. We completely characterize the capacity of multistationarity of these networks. This is especially interesting because such networks provide simple examples of "atoms of multistationarity", i.e., minimal networks that can give rise to multiple positive steady states

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Reaction Network Motifs for Static and Dynamic Absolute Concentration Robustness

Networks with absolute concentration robustness (ACR) have the property that a translation of a coordinate hyperplane either contains all steady states (static ACR) or attracts all trajectories (dynamic ACR). The implication for the underlying biological system is robustness in the concentration of one of the species independent of the initial conditions as well as independent of the concentration of all other species. Identifying network conditions for dynamic ACR is a challenging problem. We lay the groundwork in this paper by studying small reaction networks, those with 2 reactions and 2 species. We give a complete classification by ACR properties of these minimal reaction networks. The dynamics is rich even within this simple setting. Insights obtained from this work will help illuminate the properties of more complex networks with dynamic ACR.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Single-Target Networks

We characterize the dynamics of all single-target networks under mass-action kinetics: either the system is (i) globally stable for all choice of rate constants (in fact, dynamically equivalent to a detailed-balanced system) or (ii) has no positive steady states for any choice of rate constants and all trajectories must converge to the boundary of the positive orthant or to infinity. Moreover, global stability occurs if and only if the target vertex of the network is in the relative interior of the convex hull of the source vertices.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Delay stability of reaction systems

Delay differential equations are used as a model when the effect of past states has to be taken into account. In this work we consider delay models of chemical reaction networks with mass action kinetics. We obtain a sufficient condition for absolute delay stability of equilibrium concentrations, i.e., local asymptotic stability independent of the delay parameters. Several interesting examples on sequestration networks with delays are presented.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Minimal invariant regions and minimal globally attracting regions for toric differential inclusions

Toric differential inclusions occur as key dynamical systems in the context of the Global Attractor Conjecture. We introduce the notions of minimal invariant regions and minimal globally attracting regions for toric differential inclusions. We describe a procedure for constructing explicitly the minimal invariant and minimal globally attracting regions for two-dimensional toric differential inclusions. In particular, we obtain invariant regions and globally attracting regions for two-dimensional weakly reversible or endotactic dynamical systems (even if they have time-dependent parameters).

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

On classes of reaction networks and their associated polynomial dynamical systems

In the study of reaction networks and the polynomial dynamical systems that they generate, special classes of networks with important properties have been identified. These include reversible, weakly reversible}, and, more recently, endotactic networks. While some inclusions between these network types are clear, such as the fact that all reversible networks are weakly reversible, other relationships are more complicated. Adding to this complexity is the possibility that inclusions be at the level of the dynamical systems generated by the networks rather than at the level of the networks themselves. We completely characterize the inclusions between reversible, weakly reversible, endotactic, and strongly endotactic network, as well as other less well studied network types. In particular, we show that every strongly endotactic network in two dimensions can be generated by an extremally weakly reversible network. We also introduce a new class of source-only networks, which is a computationally convenient property for networks to have, and show how this class relates to the above mentioned network types.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Weakly reversible mass-action systems with infinitely many positive steady states

We show that weakly reversible mass-action systems can have a continuum of positive steady states, coming from the zeroes of a multivariate polynomial. Moreover, the same is true of systems whose underlying reaction network is reversible and has a single connected component. In our construction, we relate operations on the reaction network to the multivariate polynomial occurring as a common factor in the system of differential equations.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

An efficient characterization of complex-balanced, detailed-balanced, and weakly reversible systems

Very often, models in biology, chemistry, physics, and engineering are systems of polynomial or power-law ordinary differential equations, arising from a reaction network. Such dynamical systems can be generated by many different reaction networks. On the other hand, networks with special properties (such as reversibility or weak reversibility) are known or conjectured to give rise to dynamical systems that have special properties: existence of positive steady states, persistence, permanence, and (for well-chosen parameters) complex balancing or detailed balancing. These last two are related to thermodynamic equilibrium, and therefore the positive steady states are unique and stable. We describe a computationally efficient characterization of polynomial or power-law dynamical systems that can be obtained as complex-balanced, detailed-balanced, weakly reversible, and reversible mass-action systems.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Toric Differential Inclusions and a Proof of the Global Attractor Conjecture

The global attractor conjecture says that toric dynamical systems (i.e., a class of polynomial dynamical systems on the positive orthant) have a globally attracting point within each positive linear invariant subspace -- or, equivalently, complex balanced mass-action systems have a globally attracting point within each positive stoichiometric compatibility class. A proof of this conjecture implies that a large class of nonlinear dynamical systems on the positive orthant have very simple and stable dynamics. The conjecture originates from the 1972 breakthrough work by Fritz Horn and Roy Jackson, and was formulated in its current form by Horn in 1974. We introduce toric differential inclusions, and we show that each positive solution of a toric differential inclusion is contained in an invariant region that prevents it from approaching the origin. We use this result to prove the global attractor conjecture. In particular, it follows that all detailed balanced mass action systems and all deficiency zero weakly reversible networks have the global attractor property.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Lyapunov functions, stationary distributions, and non-equilibrium potential for chemical reaction networks

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.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Dynamical Properties of Discrete Reaction Networks

Reaction networks are commonly used to model the evolution of populations of species subject to transformations following an imposed stoichiometry. This paper focuses on the efficient characterisation of dynamical properties of Discrete Reaction Networks (DRNs). DRNs can be seen as modelling the underlying discrete nondeterministic transitions of stochastic models of reactions networks. In that sense, any proof of non-reachability in DRNs directly applies to any concrete stochastic models, independently of kinetics laws and constants. Moreover, if stochastic kinetic rates never vanish, reachability properties are equivalent in the two settings. The analysis of two global dynamical properties of DRNs is addressed: irreducibility, i.e., the ability to reach any discrete state from any other state; and recurrence, i.e., the ability to return to any initial state. Our results consider both the verification of such properties when species are present in a large copy number, and in the general case. The obtained necessary and sufficient conditions involve algebraic conditions on the network reactions which in most cases can be verified using linear programming. Finally, the relationship of DRN irreducibility and recurrence with dynamical properties of stochastic and continuous models of reaction networks is discussed.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Persistence and permanence of mass-action and power-law dynamical systems

Persistence and permanence are properties of dynamical systems that describe the long-term behavior of the solutions, and in particular specify whether positive solutions approach the boundary of the positive orthant. Mass-action systems (or more generally power-law systems) are very common in chemistry, biology, and engineering, and are often used to describe the dynamics in interaction networks. We prove that two-species mass-action systems derived from weakly reversible networks are both persistent and permanent, for any values of the reaction rate parameters. Moreover, we prove that a larger class of networks, called endotactic networks, also give rise to permanent systems, even if we allow the reaction rate parameters to vary in time. These results also apply to power-law systems and other nonlinear dynamical systems. In addition, ideas behind these results allow us to prove the Global Attractor Conjecture for three-species systems.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Product-form stationary distributions for deficiency zero chemical reaction networks

We consider stochastically modeled chemical reaction systems with mass-action kinetics and prove that a product-form stationary distribution exists for each closed, irreducible subset of the state space if an analogous deterministically modeled system with mass-action kinetics admits a complex balanced equilibrium. Feinberg's deficiency zero theorem then implies that such a distribution exists so long as the corresponding chemical network is weakly reversible and has a deficiency of zero. The main parameter of the stationary distribution for the stochastically modeled system is a complex balanced equilibrium value for the corresponding deterministically modeled system. We also generalize our main result to some non-mass-action kinetics.

0 citations0 reviewsNo code linkedOpen access