The detailed balance property and chemical systems out of equilibrium

E. Franco; J. J. L. Vel\'azquez

REVIEW 2 major objections 3 minor 24 references

The detailed balance property and chemical systems out of equilibrium

T0 review · 2 major / 3 minor · reviewed 2026-08-08 · deepseek-v4-flash

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

Pith's one-line read Every reversible chemical network can be embedded in a closed, detailed-balanced system built by adding new substances and reactions.

desk verdict A genuinely new embedding theorem for mass-action networks with a sound proof; the main caveat is interpretational, not mathematical. read the letter →

arxiv 2502.04754 v1 pith:25UPLJW2 submitted 2025-02-07 math.CA math-phmath.MP

E. Franco , J. J. L. Velázquez This is my paper · ORCID

classification math.CAmath-phmath.MP MSC 34A3437N2580A3092C42

keywords chemicalreactionnetworksdetailedbalancemassactionkineticsnon-equilibriumsystemskineticwithfluxesclosedcompletionreductionofWegscheidercriterion

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

Many biochemical models violate detailed balance, the condition that each reaction is balanced by its reverse at equilibrium. This paper argues that such non-equilibrium systems are still compatible with a thermodynamically consistent starting point: any reversible mass-action network can be obtained by freezing some concentrations in a larger closed network that satisfies detailed balance. The central result, Theorem 6.3, constructs a closed completion for every bidirectional kinetic system, adding substances and reactions so that the completion has no cycles and hence satisfies the circuit condition automatically. The reduction of that completion, with the added substances held at concentration 1, reproduces the original kinetics exactly. This matters because it makes the failure of detailed balance an effective, rather than fundamental, property of open biochemical systems.

What carries the argument

The load-bearing object is the reduction map $K[n_U]$, which extracts an effective rate for each projected reaction by summing original mass-action rates evaluated with the frozen concentrations $n_U$ held fixed, together with the complementary notion of a closed completion of the original network. Detailed balance is checked through the circuit condition: along every cycle the product of forward-to-reverse rate ratios must equal $1$. The completion construction makes all cycles disappear, so the circuit condition holds automatically, while setting the frozen concentrations of added substances to $1$ makes the reduction return the original rates. Proposition 5.1 supplies the companion statement that freezing concentrations at equilibrium values $e^{-E}$ preserves detailed balance in the reduced system.

What would settle it

Take the four-cycle $(1) \leftrightarrows (2) \leftrightarrows (3) \leftrightarrows (4) \leftrightarrows (1)$ with rate constants whose product around the cycle is not $1$, declare all four reactions constrained, and attempt to build a closed admissible completion: Proposition 6.8 says none exists, and verifying that no such completion can be written down would test the theorem's boundary. A complementary check is to simulate the completion constructed in Theorem 6.3 with added concentrations held at $1$ and confirm that the original ODEs are recovered exactly.

Watch

Extended reading notes

Core claim

On its own terms, the paper establishes that every bidirectional mass-action kinetic system $\Omega, R, K$ is admissible: there exists a closed kinetic system $\Omega_c, R_c, K_c$, detailed-balanced, conservative, and with no sources or sinks, whose reduction by freezing the added substances at concentration $1$ is exactly the original system. The proof is constructive: sources and sinks are paired with dummy substances; substances outside all conservation laws are given mirror copies to make the network conservative; and cycles are removed one by one by adding reactions and two new substances per cycle. A network without cycles satisfies the circuit condition, also called the Wegscheider criterion, for any choice of rates, so the completed system is automatically detailed-balanced. The same analysis shows that if certain reactions are constrained, a closed completion exists provided the constrained reactions do not lie in cycles; if an entire cycle whose rates violate the circuit condition is constrained, no closed admissible completion exists, as stated in Proposition 6.8.

Load-bearing premise

The theorem assumes the exchange with the environment is infinitely fast, so frozen concentrations are exactly constant; if that exchange is finite, the reduced system is only an approximation and the exact embedding statement no longer describes the open dynamics. It also assumes all reactions are reversible and obey mass-action kinetics.

Editorial extensions

If this is right

Any open bidirectional mass-action model, however far from detailed balance, has a closed, thermodynamically consistent embedding; lack of detailed balance can therefore be treated as an effective phenomenon caused by frozen concentrations.
If a cycle with nonzero energy imbalance contains only constrained reactions, the network has no closed completion, so such networks cannot be embedded without altering the constrained reactions.
Reduced systems whose frozen concentrations sit at equilibrium values inherit detailed balance and, per stoichiometric compatibility class, a unique globally stable steady state, as follows from Propositions 5.1 and 3.8.
For kinetic systems with fluxes at equilibrium frozen values, the free energy obeys $\partial_t F = -D_R + J_{\mathrm{ext}}$ and solutions relax to $e^{-E}$, with finite integrated external fluxes, as stated in Theorem 7.1.
The cycle-free completion can be chosen so that every original reaction is a one-to-one reduction and every added reaction has nonzero projection onto the original substances, making the reduction faithful.

Reading between the lines

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

If finite exchange rates are physically relevant, a natural next step is to study the completed system for large but finite $\alpha$; the exact embedding would then be the singular limit $\alpha \to \infty$, and finite-$\alpha$ corrections could be quantified.
The constructive proof suggests a quantitative measure of how far a network is from equilibrium: the minimal number of added substances and reactions needed for a closed completion, or the minimal number of reactions per cycle that must be modified. The paper does not pursue this measure.
Because the paper indicates that robust adaptation is impossible for closed chemical systems, a testable consequence is that any robustly adapting biochemical network must have a completion carrying nonzero external fluxes; this is an inference, not a claim made in the paper.
The completion with frozen concentrations fixed at $1$ effectively models the environment as infinite reservoirs, so an alternative reading is that every open mass-action network is the shadow of a conservative closed network, which could help in thermodynamic costing of biochemical functions.

Share X Bluesky LinkedIn Reddit HN

Editorial analysis

A structured set of objections, weighed in public.

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

Desk Editor's Note

read the letter

The headline is simple: this paper proves a new universal embedding theorem for mass-action networks, and the proof is sound. Any bidirectional mass-action kinetic system can be realized as a reduction of a closed, detailed-balanced system obtained by adding substances and freezing them at concentration 1. That is a real result, and I don't see a gap in the construction.

What is new is Theorem 6.3 together with the robust inheritance characterization in Section 5.2. The theorem says that every bidirectional kinetic system admits a closed completion; the earlier literature has related ideas, but not this universal construction. The paper also gives necessary and sufficient topological conditions for the reduced system to inherit detailed balance in a stable way: one-to-one reduced reactions, no zero-V reductions, and p(C)=C_V. That is a clean package, and the proofs are largely self-contained. The authors use Feinberg's book and the Wegscheider criterion as background, not as a substitute for the main argument. No data fitting, no back-computed constants; the circularity burden is minimal.

The soft spots are minor. The proof of Theorem 6.3 has the typo Cc != {0} where Proposition 6.2 proved Cc = {0}; that should be corrected. A few other small typos and wording slips appear in Section 5.2 and in the proof of Theorem 5.4, but none affect the substance. The stress-test note about the cycle-removal step in Proposition 6.2 overstating the dimension drop does not hold up on reading: when W2 = ζ - W1, on the cycle space one has W1 + W2 = ζ = 0, so the two added rows have rank one restricted to the cycles; the dimension drops by exactly one, as the paper says. The finite-α caveat is real for physical modeling, but the paper explicitly defines reduction through the α→∞ limit and does not claim finite-rate equivalence. That is an interpretational limitation, not a flaw in the completion construction.

For whom: this is for people working in chemical reaction network theory and mathematical biology, especially anyone who wants to justify out-of-equilibrium models as reductions of closed detailed-balanced systems. The flux-based measure in Section 7 is a nice extra.

Recommendation: send it to peer review. It deserves a serious referee. The referee should ask for typo fixes, a brief comment on the finite-α regime, and perhaps a short discussion connecting the topological conditions to known circuit/cycle conditions. The core result is sound and new.

Referee Report

2 major / 3 minor

Summary. The paper develops a reduction formalism for mass-action chemical reaction networks in which a subset of concentrations is frozen, and asks when such reduced systems inherit the detailed balance property from the original system. It proves that freezing at equilibrium values preserves detailed balance (Proposition 5.1), characterizes robust preservation in terms of cycle data (Theorem 5.4 and Propositions 5.6, 5.8, 5.10), and then proves the main completion theorem (Theorem 6.3): every bidirectional mass-action kinetic system is the reduction of a closed, conservative, detailed-balanced system with additional substances frozen at concentration 1. The paper also treats constrained completions and defines kinetic systems with fluxes.

Significance. If the main theorem is correct, the paper establishes a strong and useful embedding statement: arbitrary reversible open networks, even those that violate detailed balance, can be realized as fast-exchange reductions of thermodynamically consistent closed networks. This gives a rigorous framework for quantifying how far a biochemical model is from equilibrium. The proofs are self-contained, the constructions are explicit, and no fitted parameters or back-computed constants appear; the finite-α caveat in the reduction formalism is clearly stated as a modeling limitation rather than a mathematical claim. The central theorem is internally consistent, but one algebraic step in the proof of Proposition 6.2 needs correction as written.

major comments (2)

[Section 6.1, Proposition 6.2, Step 3] The displayed conservation law in Step 3 is incorrect as written. For an old reaction column R, with W1 + W2 = ζ = e_i^T Re and m_c = m_e - (1/2)m_e(i)e_i, the contribution of the old substances is -1/2 m_e(i) ζ(R), while the two new substances contribute a ζ(R) if each has coefficient a. To cancel, one needs a = 1/2 m_e(i), not 1/4 m_e(i). Thus the vector should be (m_c, 1/2 m_e(i), 1/2 m_e(i)) rather than (m_c, 1/4 m_e(i), 1/4 m_e(i)). Since conservativeness of the completed network is load-bearing for Theorem 6.3, this correction is necessary; after the correction the argument goes through, choosing m_e positive as guaranteed by conservativeness.
[Section 6.2, Proposition 6.7] The proof asserts that the row ζ used in Step 3 of Proposition 6.2 has ζ(i) = 0 for every constrained reaction R_i, but Step 3 only selects a row with two nonzero entries on the cycle being eliminated; it does not control the entries on reactions outside that cycle. Since Ra-admissibility requires W1 and W2 to vanish on constrained reactions, an additional linear-algebra argument is needed to show that a suitable row-space vector (or a different splitting) can be chosen with this vanishing property. As written, Proposition 6.7 is not fully proved.

minor comments (3)

[Theorem 6.3, proof] The sentence "hence Cc ≠ {0}" contradicts Proposition 6.2 and should read "hence Cc = {0}".
[Theorem 5.4, proof] In the final displayed computation, the phrase "since c ∈ CV" appears in a sum over the full reaction set R; the cycle membership that gives ∑ c(j)R_j(s) = 0 is the lifted cycle c ∈ C, so the notation should be corrected to avoid confusion.
[Section 1.1] The notation R* = [0,∞) and R+ = (0,∞) is nonstandard and should be flagged more prominently, since many readers will expect the opposite convention.

Circularity Check

0 steps flagged · score 0.0 of 10

No significant circularity: Theorem 6.3 is an explicit construction and the self-citations are not load-bearing.

full rationale

The central claim, Theorem 6.3 (informally Theorem 1.3), states that every bidirectional mass-action kinetic system admits a closed completion, and it is established by an explicit three-step construction rather than by importing the conclusion. In Proposition 6.2, Step 1 removes sources and sinks by adding two substances per bad reaction; Step 2 restores conservativeness by appending the coefficients (π_B R, -π_B R) for the non-conserved part B; Step 3 removes cycles by appending two reaction-index row functionals W1, W2 per basis cycle, so that the completed network has no cycles. Because the completed network has no cycles, the circuit condition of Lemma 3.2 holds for every rate function, and the identity K[n_U] = K at n_U = 1 follows directly from the reduction formula (4.5)-(4.6). Thus the reduced kinetics reproduce the original system by direct construction, with no fitted parameter, no back-computed constant, and no prediction that is equivalent to its input by definition. The self-citations to companion papers [11,12] are announcements of future or separate work and are not used in the proofs of the completion theorem; the detailed-balance facts cited from Feinberg [7,8] are standard background, not the source of the target result. The paper itself identifies the modeling limitation that reduction is defined via the α → ∞ limit of equation (1.1), so finite-rate exchange with the environment is not covered; this is a scope limitation, not a circular step. Minor typographical issues, such as the sign in the statement of Theorem 6.3 where "Cc ≠ {0}" should read "Cc = {0}" from context, do not affect the independence of the construction and are not circularity.

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

The theorems are derived from standard CRN theory plus explicit constructions. No parameters are fitted to data. The only new entities are the auxiliary substances of the completion, which are constructive devices rather than empirical postulates.

assumptions (5)

domain assumption Mass action kinetics
All kinetic systems in the paper use mass action rate laws (Section 2.2, equation (2.7)); the proofs of the detailed balance characterizations rely on this form.
standard math Detailed balance is equivalent to the circuit condition (Wegscheider criterion)
Lemma 3.2 cites Feinberg [8] and Feinberg [7] for this equivalence; it is background, not derived in the paper.
domain assumption Frozen concentrations are exactly constant (alpha -> infinity limit of equation (1.1))
Section 4 states the reduction is derived in the limit of infinitely fast exchange with the environment; the reduced ODEs are exact only in this limit.
domain assumption The original network is bidirectional
All theorems are stated for bidirectional networks (Definitions 2.1 and 2.11); the completion construction requires a reverse reaction for every reaction.
standard math Krein-Milman theorem and convex analysis results
Used in Lemma 2.9 to show that the extreme rays of the conservation-law cone form a basis of the space of conservation laws.

invented entities (1)

Auxiliary substances added in completions (Ωc \ Ω)
purpose: Allow the original open system to be represented as a reduction of a closed, detailed-balanced system by freezing these substances at concentration 1.
These are purely mathematical auxiliaries introduced by the construction in Proposition 6.2; they carry no independent falsifiable prediction and may not correspond to physically realizable chemicals.

how reviews work

0 comments

Cite this review

Pith. "Pith review of The detailed balance property and chemical systems out of equilibrium." pith.science (2026). https://pith.science/paper/25UPLJW2

@misc{pith2026250204754,
  author       = {Pith},
  title        = {Pith review of: The detailed balance property and chemical systems out of equilibrium},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/25UPLJW2}},
  note         = {Machine review of arXiv:2502.04754}
}

read the original abstract

The detailed balance property is a fundamental property that must be satisfied in all the macroscopic systems with a well defined temperature at each point. On the other hand, many biochemical networks work in non-equilibrium conditions and they can be effectively modelled using sets of equations in which the detailed balance condition fails. In this paper we study a class of "out of equilibrium" chemical networks that can be obtained freezing the concentration of some substances in chemical networks for which the detailed balance property holds. In particular, we prove that any chemical system with bidirectional chemical reactions can be extended to a system having additional substances and for which the detailed balance property holds.

Figures

Figures reproduced from arXiv: 2502.04754 by the authors.

Discussion (0). Continue with ORCID to comment.

Reference graph

Works this paper leans on

24 extracted references · 24 canonical work pages

[1]

Alon, An introduction to systems biology: design principles of biological circuits , Chapman and Hall/CRC, 2019

U. Alon, An introduction to systems biology: design principles of biological circuits , Chapman and Hall/CRC, 2019. 33

work page 2019
[2]

Barkai and S

N. Barkai and S. Leibler, Robustness in simple biochemical networks , Nature, 387 (1997), pp. 913– 917

work page 1997
[3]

Boltzmann, Lectures on Gas Theory [1896-8] , University of California Press, 1964

L. Boltzmann, Lectures on Gas Theory [1896-8] , University of California Press, 1964

work page 1964
[4]

Craciun and M

G. Craciun and M. Feinberg, Multiple equilibria in complex chemical reaction networks: I. The injectivity property, SIAM Journal on Applied Mathematics, 65 (2005), pp. 1526–1546

work page 2005
[5]

The species-reaction graph , SIAM Journal on Applied Mathematics, 66 (2006), pp

, Multiple equilibria in complex chemical reaction networks: II. The species-reaction graph , SIAM Journal on Applied Mathematics, 66 (2006), pp. 1321–1338

work page 2006
[6]

1859–1877

, Multiple equilibria in complex chemical reaction networks: semiopen mass action systems , SIAM Journal on Applied Mathematics, 70 (2010), pp. 1859–1877

work page 2010
[7]

Feinberg, Necessary and sufficient conditions for detailed balancing in mass action systems of arbitrary complexity, Chemical Engineering Science, 44 (1989), pp

M. Feinberg, Necessary and sufficient conditions for detailed balancing in mass action systems of arbitrary complexity, Chemical Engineering Science, 44 (1989), pp. 1819–1827

work page 1989
[8]

, Foundations of chemical reaction network theory , (2019)

work page 2019

Show all 24 references

[9]

J. E. Ferrell, Perfect and near-perfect adaptation in cell signaling , Cell systems, 2 (2016), pp. 62–67

2016
[10]

Flatt, D

S. Flatt, D. M. Busiello, S. Zamuner, and P. De Los Rios, ABC transporters are billion-year-old Maxwell Demons , Communications Physics, 6 (2023), p. 205

2023
[11]

Franco and J

E. Franco and J. J. Vel´ azquez,The adaptation property in non-equilibrium chemical systems , (2025)

2025
[12]

, A stochastic version of the Hopfield-Ninio kinetic proofreading model , SIAM Journal on Ap- plied Mathematics, (2025)

2025
[13]

Goldbeter and D

A. Goldbeter and D. E. Koshland Jr, An amplified sensitivity arising from covalent modification in biological systems., Proceedings of the National Academy of Sciences, 78 (1981), pp. 6840–6844

1981
[14]

Gorban and G

A. Gorban and G. Yablonsky, Extended detailed balance for systems with irreversible reactions , Chemical Engineering Science, 66 (2011), pp. 5388–5399

2011
[15]

G¨ urb¨ uz and A

B. G¨ urb¨ uz and A. D. Rendall,Analysis of a model of the Calvin cycle with diffusion of ATP , Discrete and Continuous Dynamical Systems-B, 27 (2022), pp. 5161–5177

2022
[16]

J. J. Hopfield, Kinetic proofreading: a new mechanism for reducing errors in biosynthetic processes requiring high specificity, Proceedings of the National Academy of Sciences, 71 (1974), pp. 4135–4139

1974
[17]

Milo and R

R. Milo and R. Phillips, Cell biology by the numbers , Garland Science, 2015

2015
[18]

Phillips, J

R. Phillips, J. Kondev, J. Theriot, and H. Garcia, Physical biology of the cell , Garland Science, 2012

2012
[19]

A. D. Rendall, A Calvin bestiary , in Patterns of Dynamics: Berlin, Springer, 2017, pp. 318–337

2017
[20]

A. D. Rendall and J. J. Vel´ azquez, Dynamical properties of models for the Calvin cycle , Journal of Dynamics and Differential Equations, 26 (2014), pp. 673–705

2014
[21]

R. T. Rockafellar, Convex analysis , (1970)

1970
[22]

Schnakenberg, Network theory of microscopic and macroscopic behavior of master equation sys- tems, Reviews of Modern physics, 48 (1976), p

J. Schnakenberg, Network theory of microscopic and macroscopic behavior of master equation sys- tems, Reviews of Modern physics, 48 (1976), p. 571

1976
[23]

L. A. Segel and M. Slemrod, The quasi-steady-state assumption: a case study in perturbation , SIAM Review, 31 (1989), pp. 446–477

1989
[24]

R. Wegscheider, ¨Uber simultane Gleichgewichte und die Beziehungen zwischen Thermodynamik und Reaktionskinetik homogener Systeme , Zeitschrift f¨ ur physikalische Chemie, 39 (1902), pp. 257–303. 34

1902

Pith tools

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

[1] [1]

Alon, An introduction to systems biology: design principles of biological circuits , Chapman and Hall/CRC, 2019

U. Alon, An introduction to systems biology: design principles of biological circuits , Chapman and Hall/CRC, 2019. 33

work page 2019

[2] [2]

Barkai and S

N. Barkai and S. Leibler, Robustness in simple biochemical networks , Nature, 387 (1997), pp. 913– 917

work page 1997

[3] [3]

Boltzmann, Lectures on Gas Theory [1896-8] , University of California Press, 1964

L. Boltzmann, Lectures on Gas Theory [1896-8] , University of California Press, 1964

work page 1964

[4] [4]

Craciun and M

G. Craciun and M. Feinberg, Multiple equilibria in complex chemical reaction networks: I. The injectivity property, SIAM Journal on Applied Mathematics, 65 (2005), pp. 1526–1546

work page 2005

[5] [5]

The species-reaction graph , SIAM Journal on Applied Mathematics, 66 (2006), pp

, Multiple equilibria in complex chemical reaction networks: II. The species-reaction graph , SIAM Journal on Applied Mathematics, 66 (2006), pp. 1321–1338

work page 2006

[6] [6]

1859–1877

, Multiple equilibria in complex chemical reaction networks: semiopen mass action systems , SIAM Journal on Applied Mathematics, 70 (2010), pp. 1859–1877

work page 2010

[7] [7]

Feinberg, Necessary and sufficient conditions for detailed balancing in mass action systems of arbitrary complexity, Chemical Engineering Science, 44 (1989), pp

M. Feinberg, Necessary and sufficient conditions for detailed balancing in mass action systems of arbitrary complexity, Chemical Engineering Science, 44 (1989), pp. 1819–1827

work page 1989

[8] [8]

, Foundations of chemical reaction network theory , (2019)

work page 2019