Pith. sign in

REVIEW 6 minor 30 references

Robustness in power law kinetic systems with reactant-determined interactions

T0 review · 0 major / 6 minor · reviewed 2026-08-14 · deepseek-v4-flash

Challenge this review Re-run · record.json Download PDF Read on arXiv ↗ Review JSON

Pith's one-line read A structural robustness theorem now covers power-law reaction kinetics

desk verdict A sound, honest extension of Shinar–Feinberg ACR to PL-RDK systems; the proof is correct, and the carbon-cycle application is clearly scoped to a power-law approximation. read the letter →

arxiv 1908.04497 v2 pith:4CXIU7ZB submitted 2019-08-13 math.DS

Noel T. Fortun , Angelyn R. Lao , Luis F. Razon , Eduardo R. Mendoza This is my paper · ORCID
classification math.DS MSC 37N2580A30
keywords absoluteconcentrationrobustnesschemicalreactionnetworkpowerlawkineticsreactant-determinedinteractionsdeficiency-onecarboncyclemodelShinar-Feinbergtheorem
verification ladder T0 review T1 audit T2 compute T3 formal

The pith

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The reading

This paper proves that absolute concentration robustness (ACR) — a species taking the same value at every positive steady state, independent of initial conditions — can be certified from network structure for a broad class of power-law kinetic systems, not only the mass-action systems covered by the original Shinar-Feinberg theorem. The relevant class is power-law kinetics with reactant-determined interactions (PL-RDK), where reactions that share a reactant complex share the same kinetic order vector. The main theorem states that in a deficiency-one reaction network admitting a positive equilibrium, if two nonterminal complexes have kinetic order vectors differing only in species $X_i$, then the system has ACR in $X_i$. This is a genuine extension because power-law kinetics models reactions in crowded intracellular or fractal-like media, where kinetic orders need not be the integer stoichiometric coefficients of mass action. The paper illustrates the result on a power-law approximation of a pre-industrial global carbon-cycle model, where the atmospheric carbon pool is absolutely robust when the human terrestrial carbon off-take term vanishes.

What carries the argument

The load-bearing object is the T-matrix, whose columns are the kinetic order vectors of the reactant complexes; PL-RDK is exactly the condition that these columns are well defined because reactions with the same reactant complex have identical kinetic orders. The proof leans on the Structure Theorem of the Laplacian Kernel, which supplies a basis of the Laplacian kernel supported on terminal strong linkage classes, and on the deficiency-one bound $\dim(\ker YA_\kappa)\le 1+t$. These ingredients yield Proposition 2's log-linear relation between any two positive equilibria for nonterminal complexes. Theorem 1 follows because nonterminal complexes are necessarily reactant complexes, so the relevant columns of $T$ exist, and a single-coordinate difference in those columns forces equality of the $\log$-coordinates.

What would settle it

The theorem would be falsified by a deficiency-one PL-RDK system satisfying its hypotheses but possessing two positive steady states with different $X_i$ concentrations. In the carbon setting, one concrete check is to solve the original ODE system (A.1) with $\alpha=0$ and verify whether every positive steady state satisfies $A_2=(k_2/k_1)^{1/(q_1-q_2)}$; a positive steady state with a different $A_2$ would show the approximation does not transfer ACR.

Watch

Extended reading notes

Core claim

Theorem 1 is the paper's central claim: for a deficiency-one chemical reaction network $(S,\mathcal{C},\mathcal{R})$ equipped with PL-RDK kinetics and a positive equilibrium, if two nonterminal complexes $y,y'$ have kinetic order vectors that differ only in species $X_i$, then every positive steady state assigns the same concentration to $X_i$. The proof derives, for any two positive equilibria $c^*$ and $c^{**}$, the log-linear equation $(T_{\cdot,y}-T_{\cdot,y'})\cdot \log(c^{**}/c^*)=0$, where $T$ is the reactant-complex kinetic order matrix. When the two columns differ only in coordinate $i$, this equation leaves only $\log c_i^{**}=\log c_i^*$, so the species is absolutely robust. The paper also exhibits the concrete equilibrium set of the approximated carbon-cycle scenario, in which $A_2=(k_2/k_1)^{1/(q_1-q_2)}$ is fixed while $A_1$ adjusts to conserve total carbon.

Load-bearing premise

The theorem assumes the network has deficiency one and admits a positive equilibrium; for the carbon illustration, the added load-bearing assumptions are that the power-law approximation faithfully represents the original nonlinear dynamics and that the human terrestrial off-take coefficient is exactly zero.

Editorial extensions

If this is right

  • Any PL-RDK system that meets the structural hypotheses inherits the Shinar-Feinberg guarantee: the distinguished species has one fixed concentration across all positive steady states, regardless of rate constants or initial conditions.
  • The carbon-cycle application identifies a pre-industrial scenario in which the approximated atmospheric carbon pool $A_2$ is absolutely robust, with the explicit fixed value $(k_2/k_1)^{1/(q_1-q_2)}$.
  • Because mass-action kinetics is a special case of PL-RDK, the theorem contains the original Shinar-Feinberg ACR theorem rather than merely mimicking it.
  • The hypothesis is readable off the network and kinetic-order matrix, so practitioners can certify ACR in models with fractional kinetic orders without solving the steady-state equations.

Reading between the lines

Editorial extensions of the paper, not claims the author makes directly.

  • The theorem itself is parameter-free, but the carbon conclusion is not: it relies on the power-law (GMA) approximation of the Anderies model and on setting the human terrestrial off-take coefficient to zero. A testable extension is to check whether the original logistic-exponential model also has a unique positive $A_2$ at equilibrium under $\alpha=0$.
  • Proposition 2's log-linear relation holds for any pair of nonterminal complexes, not just pairs differing in one species; this suggests a possible general robustness criterion or quantitative constraints on ratios of species concentrations, which the paper does not develop.
  • Because the condition is a comparison of two columns of $T$, robustness certification could be automated: scan a PL-RDK system's T-matrix for single-coordinate differences among nonterminal reactant complexes before doing any dynamics, which would make ACR screening practical for large networks.
Share X Bluesky LinkedIn Reddit HN

Signed reviews

No signed human review yet.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, and a circularity audit.

Referee Report

0 major / 6 minor

Summary. The paper extends the Shinar-Feinberg absolute concentration robustness (ACR) theorem from mass-action kinetics to power-law kinetics with reactant-determined interactions (PL-RDK). Theorem 1 states that for a deficiency-one CRN and a PL-RDK system admitting a positive equilibrium, if two nonterminal complexes have kinetic order vectors that differ only in species Xi, then the system has ACR in Xi. The proof is based on Proposition 2, which generalizes the original Shinar-Feinberg argument by replacing stoichiometric complexes with kinetic-order vectors encoded in the T-matrix. The paper applies the theorem to a GMA power-law approximation of the Anderies et al. pre-industrial carbon-cycle model, showing ACR in atmospheric carbon A2 when the human terrestrial carbon off-take coefficient α vanishes.

Significance. The main result is a genuine and useful generalization of a celebrated theorem. The proof is complete and transparent, relying only on standard CRNT results (the Structure Theorem of the Laplacian Kernel and the deficiency-bound Corollary 1). The carbon-cycle illustration demonstrates a plausible application, though the authors correctly limit their claim to the power-law approximation rather than the original nonlinear model. This is a solid contribution to chemical reaction network theory, of interest to researchers in CRNT and biochemical systems theory.

minor comments (6)
  1. [Section 1, last paragraph] The sentence "Particulary, for the pre-industrial scenario where there are anthropogenic causes that reduce the capacity of terrestrial carbon pool to store carbon, the power law approximation leads to an ACR-possessing PL-RDK system" is inconsistent with the derivation in Section 4 and Appendix A, where ACR is obtained when the human terrestrial carbon off-take coefficient α vanishes. Please rephrase to avoid implying that ACR arises in the presence of anthropogenic off-take.
  2. [Section 3, proof of Proposition 2, after Eq. (3.15)] The sentence "Observe that each vector bi, i = 0, 1,...,t , has its support entirely on terminal complexes" is incorrect for i=0, since b_0 = 1_C has support on all complexes. The subsequent argument only requires that b_i for i=1,...,t have terminal support; please correct the indexing.
  3. [Section 3, proof of Proposition 2, after Eq. (3.12)] The phrase "Therefore, c* and c** are positive equilibria ... if and only if Equations (3.11) and (3.12) hold" is slightly misleading, since (3.11) holds for c* by construction of κ. Consider rephrasing to state that (3.11) is equivalent to c* being an equilibrium and (3.12) is equivalent to c** being an equilibrium.
  4. [Section 5, item 2] The phrase "which accounts for the which accounts for human activities" is a typo; please remove the duplicated words.
  5. [Abstract and Section 1] There are minor typographical errors: "Particulary" should be "Particularly" in the abstract and Section 1, and "a phenomena" should be "a phenomenon."
  6. [Section 3, Proposition 3] The proof that "T = ŷ∘Y_res is also an isomorphism" is compressed; since Y_res may not be surjective onto R^m, it is more precise to say T is injective (hence the system is PL-RLK). Please elaborate or rephrase to avoid confusion.

Circularity Check

0 steps flagged · score 0.0 of 10

No significant circularity: Theorem 1 is proven from standard CRNT results, and the carbon-cycle application checks the theorem's structural condition rather than fitting the conclusion.

full rationale

The central claim is Theorem 1, which is derived from Proposition 2. Proposition 2 adapts the Shinar-Feinberg argument using the Structure Theorem of the Laplacian Kernel and Corollary 1 (dim Ker Y A_kappa <= delta + t). The proof does not assume ACR: it uses the fact that 1_C lies in Ker Y A_kappa while 1_C is not in Ker A_kappa to show equality in the dimension bound, then compares coefficients of nonterminal complexes in Eq. (3.15) to get Eq. (3.16), which becomes the T-matrix relation in Eq. (3.4). This is a genuine derivation, not the conclusion smuggled in by definition. The carbon-cycle illustration is an application, not a circular prediction. The condition p1 = p2 follows from setting the human off-take term alpha = 0, and q1 != q2 is computed from the logarithmic derivatives of the flux approximations. The ACR in A2 is then confirmed by an explicit equilibrium set in Section 4. The self-citations to [5] and [13] supply the total-CRN representation and the GMA approximation; these are explicit, checkable constructions and are not used to assert the ACR conclusion itself. No fitted parameter is renamed as a prediction, and no self-citation chain forces the result. The paper is self-contained against standard external CRNT results for its main theorem, so the appropriate finding is no circularity.

Assumptions & free parameters 3 free parameters · 5 assumptions · 0 invented entities

The theorem relies on standard CRNT results (STLK and the deficiency corollary) plus the structural deficiency-one assumption. The carbon application adds model-specific assumptions: the GMA power-law approximation faithfully represents the original carbon-cycle dynamics near the operating point, and the human terrestrial off-take term is zero (α=0). No new physical entities are introduced; all 'parameters' are either structural, computed elasticities, or standard model parameters.

free parameters (3)
  • Operating point (A1, A2, A3) = (0.69, 0.155, 0.155)
    Chosen near the computed steady state of the original model to construct the power-law approximation; the numerical kinetic orders depend on this choice, though the condition p1=p2 does not.
  • Kinetic orders p1, p2, q1, q2 = p1=p2=-68; q1=0.580148; q2=0.910864
    Computed as elasticities of the carbon fluxes (K1,K2) at the operating point. For ACR, only p1=p2 and q1≠q2 are needed; the fixed value of A2 in steady state depends on q1-q2.
  • Human terrestrial carbon off-take rate α = 0
    Set to zero to define the pre-industrial scenario; this ensures p1=p2 and is load-bearing for the application's conclusion.
assumptions (5)
  • standard math Structure Theorem of the Laplacian Kernel (STLK)
    Invoked as Proposition 1 to obtain a basis of Ker Aκ supported on terminal strong linkage classes; used in the proof of Proposition 2, Section 3.
  • standard math Corollary 1: dim(Ker Y Aκ) ≤ δ+t for a deficiency-δ network with t terminal strong linkage classes
    Used in the proof of Proposition 2 to bound the dimension of Ker Y Aκ and show {1_C,b1,...,bt} is a basis.
  • domain assumption Deficiency-one network structure
    Required by Theorem 1; the theorem does not address networks with deficiency ≠1.
  • domain assumption Validity of GMA power-law approximation of the carbon cycle model
    Section 4 and Appendix A: the ACR conclusion is about the approximated system (4.1), not the original model (A.1); the approximation is a local log-linearization.
  • domain assumption Existence of a positive equilibrium
    Assumed by Theorem 1; for the carbon application a steady state is computed numerically, and positivity requires the conserved total carbon A0 to exceed (1+1/β)(k2/k1)^(1/(q1-q2)).

how reviews work

0 comments
Cite this review

Pith. "Pith review of Robustness in power law kinetic systems with reactant-determined interactions." pith.science (2026). https://pith.science/paper/4CXIU7ZB

@misc{pith2026190804497,
  author       = {Pith},
  title        = {Pith review of: Robustness in power law kinetic systems with reactant-determined interactions},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/4CXIU7ZB}},
  note         = {Machine review of arXiv:1908.04497}
}
read the original abstract

Robustness against the presence of environmental disruptions can be observed in many systems of chemical reaction network. However, identifying the underlying components of a system that give rise to robustness is often elusive. The influential work of Shinar and Feinberg established simple yet subtle network-based conditions for absolute concentration robustness (ACR), a phenomena in which a species in a mass-action system has the same concentration for any steady state the network may admit. In this contribution, we extend this result to embrace kinetic systems more general than mass-action systems, namely, power-law kinetic systems with reactant-determined interactions (denoted by "PL-RDK"). In PL-RDK, the kinetic order vectors (which we call "interactions") of reactions with the same reactant complex are identical. As illustration, we considered a scenario in the pre-industrial state of global carbon cycle. A power-law approximation of the dynamical system of this scenario is found to be dynamically equivalent to an ACR-possessing PL-RDK system.

Figures

Figures reproduced from arXiv: 1908.04497 by the authors.

Figure 4.1
Figure 4.1. Biochemical map of the pre-industrial carbon cycle model of Anderies et al. [1] [PITH_FULL_IMAGE:figures/full_fig_p009_4_1.png] view at source ↗

Discussion (0). Continue with ORCID to comment.

Reference graph

Works this paper leans on

30 extracted references · 30 canonical work pages

  1. [13]

    Journal of Mathematical Chemistry 56(10), 2929–2962 (2018)

    Fortun, N., Lao, A., Razon, L., Mendoza, E.: A deficiency-one algorithm for a power law kinetics with reactant determined interactions. Journal of Mathematical Chemistry 56(10), 2929–2962 (2018)

    work page 2018

  2. [1]

    Environmental Research Letters 8(4), 044–048 (2013)

    Anderies, J., Carpenter, S., Steffen, W., Rockstr¨ om, J.: The topology of non-linear global car- bon dynamics: from tipping points to planetary boundaries. Environmental Research Letters 8(4), 044–048 (2013)

    work page 2013

  3. [2]

    J R Soc Interface 11(93) (2014)

    Anderson, D.F., Enciso, G.A., Johnston, M.D.: Stochastic analysis of biochemical reaction networks with absolute concentration robustness. J R Soc Interface 11(93) (2014)

    work page 2014

  4. [3]

    Mathematical biosciences 283, 13–29 (2017)

    Arceo, C.P.P., Jose, E.C., Lao, A., Mendoza, E.R.: Reaction networks and kinetics of bio- chemical systems. Mathematical biosciences 283, 13–29 (2017)

    work page 2017

  5. [4]

    Journal of Mathematical Chemistry 56(2), 395–422 (2018)

    Arceo, C.P.P., Jose, E.C., Lao, A., Mendoza, E.R.: Reactant subspaces and kinetics of chemical reaction networks. Journal of Mathematical Chemistry 56(2), 395–422 (2018)

    work page 2018

  6. [5]

    Mathematical Biosciences 269, 135–52 (2015)

    Arceo, C.P.P., Jose, E.C., Mar´ ın-Sanguino, A., Mendoza, E.R.: Chemical reaction network approaches to biochemical systems theory. Mathematical Biosciences 269, 135–52 (2015)

    work page 2015

  7. [6]

    Mathematical Biosciences 215(1), 35 – 47 (2008)

    Zeljko Bajzer, Huzak, M., Neff, K.L., Prendergast, F.G.: Mathematical analysis of models for reaction kinetics in intracellular environments. Mathematical Biosciences 215(1), 35 – 47 (2008)

    work page 2008

  8. [7]

    Cell Biophysics 12, 237–253 (1988)

    Clarke, B.L.: Stoichiometric network analysis. Cell Biophysics 12, 237–253 (1988)

    work page 1988

Show all 30 references
  1. [8]

    In: Stadtman, E.R., Chock, P.B

    Clegg, J.S.: Cellular infrastructure and metabolic organization. In: Stadtman, E.R., Chock, P.B. (eds.) From Metabolite, to Metabolism, to Metabolon. Current Topics in Cellular Regu- lation, vol. 33, pp. 3–14. Academic Press (1992)

    1992

  2. [9]

    Journal of Mathematical Chemistry 56(2), 336–357 (2018)

    Cortez, M.J., Nazareno, A., Mendoza, E.: A computational approach to linear conjugacy in a class of power law kinetic systems. Journal of Mathematical Chemistry 56(2), 336–357 (2018)

    2018

  3. [10]

    Dexter, J.P., Dasgupta, T., Gunawardena, J.: Invariants reveal multiple forms of robustness in bifunctional enzyme systems. Integr. Biol. 7, 883–894 (2015) 11

    2015

  4. [11]

    Notes of lectures given at the mathe- matics research center of the University of Wisconsin (1979)

    Feinberg, M.: Lectures on chemical reaction networks. Notes of lectures given at the mathe- matics research center of the University of Wisconsin (1979)

    1979

  5. [12]

    Feinberg, M.: The existence and uniqueness of steady states for a class of chemical reaction networks. Arch. Ration. Mech. Anal. 132, 311–370 (1995)

    1995

  6. [14]

    Horn, F., Jackson, R.: General mass action kinetics. Arch. Rational Mech. Anal 47, 187–194 (1972)

    1972

  7. [15]

    Nature Reviews Genetics 5(11), 826–837 (2004)

    Kitano, H.: Biological robustness. Nature Reviews Genetics 5(11), 826–837 (2004)

    2004

  8. [16]

    Journal of Statistical Physics 42, 185–200 (1986)

    Kopelman, R.: Rate processes on fractals: theory, simulations, and experiments. Journal of Statistical Physics 42, 185–200 (1986)

    1986

  9. [17]

    Science 241 4873, 1620–1626 (1988)

    Kopelman, R.: Fractal reaction kinetics. Science 241 4873, 1620–1626 (1988)

    1988

  10. [18]

    Israel Journal of Chemistry 31(2), 147–157 (1991)

    Kopelman, R., Koo, Y.: Reaction kinetics in restricted spaces. Israel Journal of Chemistry 31(2), 147–157 (1991)

    1991

  11. [19]

    Progress in Biophysics and Molecular Biology 75(1), 1 – 17 (2001)

    Kuthan, H.: Self-organisation and orderly processes by individual protein complexes in the bacterial cell. Progress in Biophysics and Molecular Biology 75(1), 1 – 17 (2001)

    2001

  12. [20]

    SIAM Journal of Applied Mathematics 72, 1926–1947 (2012)

    M¨ uller, S., Regensburger, G.: Generalized mass action systems: Complex balancing equilibria and sign vectors of the stoichiometric and kinetic-order subspaces. SIAM Journal of Applied Mathematics 72, 1926–1947 (2012)

    2012

  13. [21]

    In: Proceedings of the International Workshop on Computer Algebra in Scientific Computing(CASC) (2014)

    M¨ uller, S., Regensburger, G.: Generalized mass-action systems and positive solutions of poly- nomial equations with real and symbolic exponents (invited talk). In: Proceedings of the International Workshop on Computer Algebra in Scientific Computing(CASC) (2014)

    2014

  14. [22]

    The Journal of Physical Chemistry 92(6), 1538–1541 (1988)

    Newhouse, J.S., Kopelman, R.: Steady-state chemical kinetics on surface clusters and islands: segregation of reactants. The Journal of Physical Chemistry 92(6), 1538–1541 (1988)

    1988

  15. [23]

    some mathematical properties of the rate law for the component enzymatic reactions

    Savageau, M.A.: Biochemical systems analysis: I. some mathematical properties of the rate law for the component enzymatic reactions. American Journal of Science25(3), 365–369 (1969)

    1969

  16. [24]

    Biosystems 47(1), 9 – 36 (1998)

    Savageau, M.A.: Development of fractal kinetic theory for enzyme-catalysed reactions and implications for the design of biochemical pathways. Biosystems 47(1), 9 – 36 (1998)

    1998

  17. [25]

    Progress in biophysics and molecular biology 85 2–3 , 235–260 (2004)

    Schnell, S., Turner, T.E.: Reaction kinetics in intracellular environments with macromolecular crowding: simulations and rate laws. Progress in biophysics and molecular biology 85 2–3 , 235–260 (2004)

    2004

  18. [26]

    Science 327(5971), 1389–1391 (2010)

    Shinar, G., Feinberg, M.: Structural sources of robustness in biochemical reaction networks. Science 327(5971), 1389–1391 (2010)

    2010

  19. [27]

    Journal of Mathematical Chemistry 56(2), 358–394 (2018) 12

    Talabis, D.A.S.J., Arceo, C.P.P., Mendoza, E.R.: Positive equilibria of a class of power law kinetics. Journal of Mathematical Chemistry 56(2), 358–394 (2018) 12

    2018

  20. [28]

    Cambridge University Press (2000)

    Voit, E.: Computational analysis of biochemical systems: a practical guide for biochemists and molecular biologists. Cambridge University Press (2000)

    2000

  21. [29]

    ISRN Biomathematics 2013, 1–53 (2013)

    Voit, E.: Biochemical systems theory: A review. ISRN Biomathematics 2013, 1–53 (2013)

    2013

  22. [30]

    Wiuf, C., Feliu, E.: Power-law kinetics and determinant criteria for the preclusion of mul- tistationarity in networks of interacting species. SIAM J. Applied Dynamical Systems 12, 1685–1721 (2013) A Pre-industrial Carbon Cycle Model of Anderies et al. The complete set of ODEs...

    2013

Pith tools

Reviewed August 14, 2026 · model on record in the stance chip above.