BZPEERReading, trust and collaboration for research
Menu

Discover

GraphFollow connections around a paper, topic or saved view.

Workspaces

HomePick up where your research flow left off.InboxHandle requests, replies and notifications from one place.CollectionsCurate reading lists, evidence sets and shareable dossiers.ResearchSee your private provenance graph, trail and imports.CollaborationMove from discovery into focused discussion rooms.PublishDraft, preview and publish full-text research articles.

Network

ExpertsFind specialists worth following or asking for help.CommunitiesJoin research circles organized around shared work.PartnersSee external organizations active around the same topics.

Opportunities

ListingsBrowse roles, calls and collaborations with clearer fit signals.

Account

Log inReturn to your reading and collaboration workspace.Create accountStart saving work, following people and joining teams.
Log inCreate account

Researcher profile

Bryan S. Hernandez

Bryan S. Hernandez contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate
math.DS

Trust snapshot

Quick read

Trust 21 - Emerging
7works
0followers
1topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

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

7 published item(s)

preprint2022arXiv

Positive equilibria of power law kinetics on networks with independent linkage classes

Studies about the set of positive equilibria ($E_+$) of kinetic systems have been focused on mass action, and not that much on power law kinetic (PLK) systems, even for PL-RDK systems (PLK systems where two reactions with identical reactant complexes have the same kinetic order vectors). For mass action, reactions with different reactants have different kinetic order rows. A PL-RDK system satisfying this property is called factor span surjective (PL-FSK). In this work, we show that a cycle terminal PL-FSK system with $E_+\ne \varnothing$ and has independent linkage classes (ILC) is a poly-PLP system, i.e., $E_+$ is the disjoint union of log-parametrized sets. The key insight for the extension is that factor span surjectivity induces an isomorphic digraph structure on the kinetic complexes. The result also completes, for ILC networks, the structural analysis of the original complex balanced generalized mass action systems (GMAS) by Müller and Regensburger. We also identify a large set of PL-RDK systems where non-emptiness of $E_+$ is a necessary and sufficient condition for non-emptiness of each set of positive equilibria for each linkage class. These results extend those of Boros on mass action systems with ILC. We conclude this paper with two applications of our results. Firstly, we consider absolute complex balancing (ACB), i.e., the property that each positive equilibrium is complex balanced, in poly-PLP systems. Finally, we use the new results to study absolute concentration robustness (ACR) in these systems. In particular, we obtain a species hyperplane containment criterion to determine ACR in the system species.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Weakly Reversible CF-Decompositions of Chemical Kinetic Systems

This paper studies chemical kinetic systems which decompose into weakly reversible complex factorizable (CF) systems. Among power law kinetic systems, CF systems (denoted as PL-RDK systems) are those where branching reactions of a reactant complex have identical rows in the kinetic order matrix. Mass action and generalized mass action systems (GMAS) are well-known examples. Schmitz's global carbon cycle model is a previously studied non-complex factorizable (NF) power law system (denoted as PL-NDK). We derive novel conditions for the existence of weakly reversible CF-decompositions and present an algorithm for verifying these conditions. We discuss methods for identifying independent decompositions, i.e., those where the stoichiometric subspaces of the subnetworks form a direct sum, as such decompositions relate positive equilibria sets of the subnetworks to that of the whole network. We then use the results to determine the positive equilibria sets of PL-NDK systems which admit an independent weakly reversible decomposition into PL-RDK systems of PLP type, i.e., the positive equilibria are log-parametrized, which is a broad generalization of a Deficiency Zero Theorem of Fortun et al. (2019).

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

A Computational Approach to Multistationarity in Poly-PL Kinetic Systems

One important question that interests those who work in chemical reaction network theory (CRNT) is this: Does the system obtained from a reaction network admit a positive equilibrium and if it does, can there be more than one within a stoichiometric class? The higher deficiency algorithm (HDA) of Ji and Feinberg provided a method of determining the multistationarity capacity of a CRN with mass action kinetics (MAK). An extension of this, called Multistationarity Algorithm (MSA), recently came into the scene tackling CRNs with power law kinetics (PLK), a kinetic system which is more general (having MAK systems as a special case). For this paper, we provide a computational approach to study the multistationarity feature of reaction networks endowed with kinetics which are non-negative linear combinations of power law functions called poly-PL kinetics (PYK). The idea is to use MSA and combine it with a transformation called STAR-MSC (i.e., $S$-invariant Termwise Addition of Reactions via Maximal Stoichiometric Coefficients) producing PLKs that are dynamically equivalent to PYKs. This leads us to being able to determinine the multistationarity capacity of a much larger class of kinetic systems. We show that if the transformed dynamically equivalent PLK system is multistationary for a stoichiometric class for a set of particular rate constants, then so is its original corresponding PYK system. Moreover, the monostationarity property of the transformed PLK system also implies the monostationarity property of the original PYK system.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

A Computational Approach to Multistationarity of Power-Law Kinetic Systems

This paper presents a computational solution to determine if a chemical reaction network endowed with power-law kinetics (PLK system) has the capacity for multistationarity, i.e., whether there exist positive rate constants such that the corresponding differential equations admit multiple positive steady states within a stoichiometric class. The approach, which is called the "Multistationarity Algorithm for PLK systems" (MSA), combines (i) the extension of the "higher deficiency algorithm" of Ji and Feinberg for mass action to PLK systems with reactant-determined interactions, and (ii) a method that transforms any PLK system to a dynamically equivalent one with reactant-determined interactions. Using this algorithm, we obtain two new results: the monostationarity of a popular model of anaerobic yeast fermentation pathway, and the multistationarity of a global carbon cycle model with climate engineering, both in the generalized mass action format of biochemical systems theory. We also provide examples of the broader scope of our approach for deficiency one PLK systems in comparison to the extension of Feinberg's "deficiency one algorithm" to such systems.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Fundamental Decompositions and Multistationarity of Power-Law Kinetic Systems

The fundamental decomposition of a chemical reaction network (also called its "$\mathscr{F}$-decomposition") is the set of subnetworks generated by the partition of its set of reactions into the "fundamental classes" introduced by Ji and Feinberg in 2011 as the basis of their "higher deficiency algorithm" for mass action systems. The first part of this paper studies the properties of the $\mathscr{F}$-decomposition, in particular, its independence (i.e., the network's stoichiometric subspace is the direct sum of the subnetworks' stoichiometric subspaces) and its incidence-independence (i.e., the image of the network's incidence map is the direct sum of the incidence maps' images of the subnetworks). We derive necessary and sufficient conditions for these properties and identify network classes where the $\mathscr{F}$-decomposition coincides with other known decompositions. The second part of the paper applies the above-mentioned results to improve the Multistationarity Algorithm for power-law kinetic systems (MSA), a general computational approach that we introduced in previous work. We show that for systems with non-reactant determined interactions but with an independent $\mathscr{F}$-decomposition, the transformation to a dynamically equivalent system with reactant-determined interactions -- required in the original MSA -- is not necessary. We illustrate this improvement with the subnetwork of Schmitz's carbon cycle model recently analyzed by Fortun et al.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

On the Independence of Fundamental Decompositions of Power-Law Kinetic Systems

The fundamental decomposition of a chemical reaction network (CRN) is induced by partitioning the reaction set into "fundamental classes". It was the basis of the Higher Deficiency Algorithm for mass action systems of Ji and Feinberg, and the Multistationarity Algorithm for power-law kinetic systems of Hernandez et al. In addition to our previous work, we provide important properties of the independence (i.e., the network's stoichiometric subspace is the direct sum of the subnetworks' stoichiometric subspaces) and the incidence-independence (i.e., the image of the network's incidence map is the direct sum of the incidence maps' images of the subnetworks) of these decompositions. Feinberg established the essential relationship between independent decompositions and the set of positive equilibria of a network, which we call the Feinberg Decomposition Theorem (FDT). Moreover, Farinas et al. recently documented its version for incidence-independence. Fundamental decomposition divides the network into subnetworks of deficiency either 0 or 1 only. Hence, available results for lower deficiency networks, such as the Deficiency Zero Theorem (DZT), can be used. These justify the study of independent fundamental decompositions. A MATLAB program which (i) computes the subnetworks of a CRN under the fundamental decomposition and (ii) is useful for determining whether the decomposition is independent and incidence-independent is also created. Finally, we provide the following solution for determining multistationarity of CRNs with the following steps: (1) the use of the program, (2) the application of available results for CRNs with deficiency 0 or 1 (e.g., DZT), and (3) the use of FDT. We illustrate the solution by showing that the generalization of a subnetwork of Schmitz's carbon cycle model by Hernandez et al., endowed with mass action kinetics, has no capacity for multistationarity.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Positive Equilibria of Hill-Type Kinetic Systems

This work introduces a novel approach to study properties of positive equilibria of a chemical reaction network $\mathscr{N}$ endowed with Hill-type kinetics $K$, called a Hill-type kinetic (HTK) system $\left(\mathscr{N},K\right)$, including their multiplicity and concentration robustness in a species. We associate a unique positive linear combination of power-law kinetic systems called poly-PL kinetic (PYK) system $\left( {\mathscr{N},{K_\text{PY}}} \right)$ to the given HTK system. The associated system has the key property that its equilibria sets coincide with those of the Hill-type system, i.e., ${E_ + }\left( {\mathscr{N},K} \right) = {E_ + }\left( {\mathscr{N},{K_\text{PY}}} \right)$ and ${Z_ + }\left( {\mathscr{N},K} \right) = {Z_ + }\left( {\mathscr{N},{K_\text{PY}}} \right)$. This allows us to identify two novel subsets of the Hill-type kinetics, called PL-equilibrated and PL-complex balanced kinetics, to which recent results on absolute concentration robustness (ACR) of species and complex balancing at positive equilibria of power-law (PL) kinetic systems can be applied. Our main results also include the Shinar-Feinberg ACR Theorem for PL-equilibrated HT-RDK systems (i.e., subset of complex factorizable HTK systems), which establishes a foundation for the analysis of ACR in HTK systems, and the extension of the results of Müller and Regensburger on generalized mass action systems to PL-complex balanced HT-RDK systems. In addition, we derive the theory of balanced concentration robustness (BCR) in an analogous manner to ACR for PL-equilibrated systems. Finally, we provide further extensions of our results to a more general class of kinetics, which includes quotients of poly-PL functions.

0 citations0 reviewsNo code linkedOpen access