Researcher profile

Gheorghe Craciun

Gheorghe Craciun contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
13works
0followers
5topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

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

Published work

13 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.

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.

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.

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.

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

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.

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.

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.

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).

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.

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.

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.

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.