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

David Siegel

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

ResearcherUnclaimed source record
math.DSmath.CAMolecular Networksphysics.chem-ph

Catalog footprint

What is connected

9works
4topics
3close 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

9 published item(s)

preprint2012arXiv

Computing Weakly Reversible Linearly Conjugate Chemical Reaction Networks with Minimal Deficiency

Mass-action kinetics is frequently used in systems biology to model the behaviour of interacting chemical species. Many important dynamical properties are known to hold for such systems if they are weakly reversible and have a low deficiency. In particular, the Deficiency Zero and Deficiency One Theorems guarantee strong regularity with regards to the number and stability of positive equilibrium states. It is also known that chemical reaction networks with disparate reaction structure can exhibit the same qualitative dynamics. The theory of linear conjugacy encapsulates the cases where this relationship is captured by a linear transformation. In this paper, we propose a mixed-integer linear programming algorithm capable of determining weakly reversible reaction networks with a minimal deficiency which are linearly conjugate to a given reaction network.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

A Linear Programming Approach to Weak Reversibility and Linear Conjugacy of Chemical Reaction Networks

A numerically effective procedure for determining weakly reversible chemical reaction networks that are linearly conjugate to a known reaction network is proposed in this paper. The method is based on translating the structural and algebraic characteristics of weak reversibility to logical statements and solving the obtained set of linear (in)equalities in the framework of mixed integer linear programming. The unknowns in the problem are the reaction rate coefficients and the parameters of the linear conjugacy transformation. The efficacy of the approach is shown through numerical examples.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Asymptotic Behaviour Near a Nonlinear Sink

In this paper, we will develop an iterative procedure to determine the detailed asymptotic behaviour of solutions of a certain class of nonlinear vector differential equations which approach a nonlinear sink as time tends to infinity. This procedure is indifferent to resonance in the eigenvalues. Moreover, we will address the writing of one component of a solution in terms of the other in the case of a planar system. Examples will be given, notably the Michaelis-Menten mechanism of enzyme kinetics.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Dynamical Equivalence and Linear Conjugacy of Chemical Reaction Networks: New Results and Methods

In the first part of this paper, we propose new optimization-based methods for the computation of preferred (dense, sparse, reversible, detailed and complex balanced) linearly conjugate reaction network structures with mass action dynamics. The developed methods are extensions of previously published results on dynamically equivalent reaction networks and are based on mixed-integer linear programming. As related theoretical contributions we show that (i) dense linearly conjugate networks define a unique super-structure for any positive diagonal state transformation if the set of chemical complexes is given, and (ii) the existence of linearly conjugate detailed balanced and complex balanced networks do not depend on the selection of equilibrium points. In the second part of the paper it is shown that determining dynamically equivalent realizations to a network that is structurally fixed but parametrically not can also be written and solved as a mixed-integer linear programming problem. Several examples illustrate the presented computation methods.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Linear Conjugacy of Chemical Reaction Networks

Under suitable assumptions, the dynamic behaviour of a chemical reaction network is governed by an autonomous set of polynomial ordinary differential equations over continuous variables representing the concentrations of the reactant species. It is known that two networks may possess the same governing mass-action dynamics despite disparate network structure. To date, however, there has only been limited work exploiting this phenomenon even for the cases where one network possesses known dynamics while the other does not. In this paper, we bring these known results into a broader unified theory which we call conjugate chemical reaction network theory. We present a theorem which gives conditions under which two networks with different governing mass-action dynamics may exhibit the same qualitative dynamics and use it to extend the scope of the well-known theory of weakly reversible systems.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Weak Dynamic Non-Emptiability and Persistence of Chemical Kinetics Systems

A frequently desirable characteristic of chemical kinetics systems is that of persistence, the property that if all the species are initially present then none of them may tend toward extinction. It is known that solutions of deterministically modelled mass-action systems may only approach portions of the boundary of the positive orthant which correspond to semi-locking sets (alternatively called siphons). Consequently, most recent work on persistence of these systems has been focused on these sets. In this paper, we focus on a result which states that, for a conservative mass-action system, persistence holds if every critical semi-locking set is dynamically non-emptiable and the system contains no nested locking sets. We will generalize this result by introducing the notion of a weakly dynamically non-emptiable semi-locking set and making novel use of the well-known Farkas' Lemma. We will also connect this result to known results regarding complex balanced systems and systems with facets.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

A Stratum Approach to Global Stability of Complex Balanced Systems

It has long been known that complex balanced mass-action systems exhibit a restrictive form of behaviour known as locally stable dynamics. This means that within each compatibility class $\mathcal{C}_{\mathbf{x}_0}$---the forward invariant space where solutions lies---there is exactly one equilibrium concentration and that this concentration is locally asymptotically stable. It has also been conjectured that this stability extends globally to $\mathcal{C}_{\mathbf{x}_0}$. That is to say, all solutions originating in $\mathcal{C}_{\mathbf{x}_0}$ approach the unique positive equilibrium concentration rather than $\partial \mathcal{C}_{\mathbf{x}_0}$ or $\infty$. To date, however, no general proof of this conjecture has been found. In this paper, we approach the problem of global stability for complex balanced systems through the methodology of dividing the positive orthant into regions called strata. This methodology has been previously applied to detailed balanced systems---a proper subset of complex balanced systems---to show that, within a stratum, trajectories are repelled from any face of $\mathbb{R}_{\geq 0}^m$ adjacent to the stratum. Several known global stability results for detailed balanced systems are generalized to complex balanced systems.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Properties of the Lindemann Mechanism in Phase Space

We study the planar and scalar reductions of the nonlinear Lindemann mechanism of unimolecular decay. First, we establish that the origin, a degenerate critical point, is globally asymptotically stable. Second, we prove there is a unique scalar solution (the slow manifold) between the horizontal and vertical isoclines. Third, we determine the concavity of all scalar solutions in the nonnegative quadrant. Fourth, we establish that each scalar solution is a centre manifold at the origin given by a Taylor series. Moreover, we develop the leading-order behaviour of all planar solutions as time tends to infinity. Finally, we determine the asymptotic behaviour of the slow manifold at infinity by showing that it is a unique centre manifold for a fixed point at infinity.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Properties of the Michaelis-Menten Mechanism in Phase Space

We study the two-dimensional reduction of the Michaelis-Menten reaction of enzyme kinetics. First, we prove the existence and uniqueness of a slow manifold between the horizontal and vertical isoclines. Second, we determine the concavity of all solutions in the first quadrant. Third, we establish the asymptotic behaviour of all solutions near the origin, which generally is not given by a Taylor series. Finally, we determine the asymptotic behaviour of the slow manifold at infinity.

0 citations0 reviewsNo code linkedOpen access