pith. sign in

arxiv: 2607.02367 · v1 · pith:F2TCI7NZnew · submitted 2026-07-02 · 🧮 math.PR · math-ph· math.CA· math.MP

Flux solutions for stochastic chemical systems with sources and sinks

E. Franco , J. J. L. Vel\'azquez This is my paper

Pith reviewed 2026-07-03 06:30 UTC · model grok-4.3

classification 🧮 math.PR math-phmath.CAmath.MP
keywords stochastic chemical systemssources and sinksMarkov chainsnon-explosionstationary distributionsout-of-equilibrium statesmembrane transportcontinuous-time processes
0
0 comments X

The pith

Adding sources and sinks to chemical reaction networks yields non-explosive Markov chains that converge to unique out-of-equilibrium stationary distributions.

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

The paper studies stochastic chemical systems obtained by augmenting conservative reaction networks with source and sink terms. These modifications produce continuous-time Markov chains on the non-negative integer lattice that violate detailed balance and complex balance. The authors prove the chains are non-explosive and converge to a unique stationary distribution as time tends to infinity. The resulting stationary measures describe out-of-equilibrium states that support persistent fluxes. The framework is applied to compute molecular transport through membrane channels.

Core claim

By adding sources and sinks to conservative chemical systems the authors obtain well-defined continuous-time Markov chains whose transition rates meet the conditions for non-explosion. Standard theorems then imply convergence to a unique invariant probability measure. This measure is an out-of-equilibrium stationary solution because the source-sink terms sustain net fluxes through the system.

What carries the argument

The modified transition rates created by adjoining source and sink Poisson processes to a conservative reaction network, which preserve the rate conditions needed for non-explosion and ergodicity on the countable state space.

If this is right

  • Every such augmented system admits a unique invariant probability measure.
  • The long-time distribution of molecule counts is given by this measure independently of the initial state.
  • Net molecular fluxes through any reaction or channel equal the expectations computed under the stationary measure.
  • The construction applies to any finite set of reactions supplemented by constant-rate reservoirs.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

  • The same rate conditions could be checked numerically for larger networks before proving convergence analytically.
  • The stationary distributions might serve as a baseline for comparing stochastic and deterministic flux predictions in open systems.
  • Extensions that include spatial diffusion or multiple compartments would require only local rate bounds of the same type.

Load-bearing premise

The birth and death rates induced by the sources and sinks must obey the growth bounds required by the general non-explosion and convergence theorems for Markov chains on countable spaces.

What would settle it

An explicit choice of source and sink rates for which the total molecule count reaches infinity in finite time with positive probability would disprove non-explosion.

Figures

Figures reproduced from arXiv: 2607.02367 by E. Franco, J. J. L. Vel\'azquez.

Figure 1
Figure 1. Figure 1: The system can jumps from the state (n, N, nX) = (3, 6, 1, 0) to the state (n, N, nX) = (3, 5, 0, 1) due to the reaction ρ1. It can then jump to state (n, N, nX) = (4, 5, 0, 1) via the reaction ρ2. The state of the system is characterized by the triple (n, N, nX) where n, N ∈ N and nX ∈ N 2 is the number of channels at state 0 and at state 1. In particular, we assume that initially we have only one channel… view at source ↗
Figure 2
Figure 2. Figure 2: The channel at the 8 different states. The state of the channel in the figure in the third line and [PITH_FULL_IMAGE:figures/full_fig_p031_2.png] view at source ↗
read the original abstract

In this paper we study a class of stochastic chemical systems that, in general, do not satisfy the property of detailed balance nor the property of complex balance. These systems are obtained by adding sources and sinks to conservative chemical systems. This procedure is a way to define rigorously stochastic chemical systems in contact with reservoirs. We prove that these systems are non-explosive Markov chains and we prove that they converge to a steady state as time tends to infinity. The stationary solution are out of equilibrium solutions. We conclude the paper by applying our results in order to describe fluxes of molecules through some membrane channels.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

1 major / 2 minor

Summary. The manuscript studies stochastic chemical reaction networks obtained by adding sources and sinks to conservative CRNs. It claims to prove that the resulting systems are non-explosive continuous-time Markov chains on the non-negative integer lattice that converge to a unique stationary distribution as t → ∞, with these stationary distributions being out of equilibrium. The results are applied to model molecular fluxes through membrane channels.

Significance. If the proofs are complete and the rate conditions are verified, the work supplies a rigorous probabilistic treatment of open stochastic chemical systems in contact with reservoirs, extending conservative CRN theory to cases lacking detailed or complex balance. The membrane-channel application illustrates potential utility in biophysics.

major comments (1)
  1. [Abstract and Introduction] Abstract and Introduction: The claim that the modified rates 'satisfy the conditions needed' for non-explosion and ergodicity is asserted without exhibiting the explicit verification. For general source/sink intensities on ℕ^d (e.g., state-independent sources paired with linear sinks in multi-species systems), it is not shown that q(x) < ∞ for all x or that a Lyapunov drift condition QV(x) ≤ K − εV(x) holds outside a finite set, leaving the invocation of standard theorems (Norris, Reuter, or equivalent) unsecured. This is load-bearing for the central existence and convergence statements.
minor comments (2)
  1. The notation for the state space, transition rates, and source/sink terms should be introduced with a dedicated preliminary section before the main theorems.
  2. Clarify whether the stationary distributions are shown to be unique or merely that convergence occurs to some stationary measure.

Simulated Author's Rebuttal

1 responses · 0 unresolved

We thank the referee for their thorough review and valuable comments on our manuscript. We address the major concern regarding the explicit verification of conditions for non-explosion and ergodicity below. We believe the proofs in the paper are complete, but agree that additional clarity in the introduction would strengthen the presentation.

read point-by-point responses

  1. Referee: [Abstract and Introduction] The claim that the modified rates 'satisfy the conditions needed' for non-explosion and ergodicity is asserted without exhibiting the explicit verification. For general source/sink intensities on ℕ^d (e.g., state-independent sources paired with linear sinks in multi-species systems), it is not shown that q(x) < ∞ for all x or that a Lyapunov drift condition QV(x) ≤ K − εV(x) holds outside a finite set, leaving the invocation of standard theorems (Norris, Reuter, or equivalent) unsecured. This is load-bearing for the central existence and convergence statements.

    Authors: The manuscript provides complete proofs of non-explosion and convergence in Sections 3 and 4, where we verify that the generator satisfies the conditions of the cited theorems for the specific form of source and sink rates considered. These rates are not arbitrary but derived from adding sources and sinks to conservative CRNs, which ensures q(x) remains finite and allows construction of a suitable Lyapunov function V (typically a linear or quadratic function in the species counts) to obtain the required drift. We concede that the abstract and introduction do not exhibit this verification explicitly, which may leave the application of the theorems less transparent. We will revise the introduction to include a short paragraph summarizing the key steps: (i) finiteness of q(x) follows from the boundedness of reaction rates in each state due to the conservative part plus linear growth from sinks; (ii) the drift condition is established by choosing V(x) = sum x_i and computing QV(x) explicitly using the structure of the reactions. This revision will be made. Note that the results are for the class of systems described, not necessarily all possible general intensities, but the example given (state-independent sources with linear sinks) is covered by our assumptions. revision: yes

Circularity Check

0 steps flagged

No circularity: standard Markov chain existence and ergodicity proofs applied to modified rates

full rationale

The paper establishes non-explosion and convergence to stationary distributions for CTMCs obtained by adding sources/sinks to conservative chemical reaction networks. These are direct applications of standard theorems (e.g., on generators satisfying drift conditions or Reuter/Norris criteria) to the modified transition rates. No equations reduce to fitted parameters, no self-definitional loops, and no load-bearing self-citations that substitute for verification. The central claims are existence statements about the resulting Markov process, which remain independent of the paper's own inputs.

Axiom & Free-Parameter Ledger

0 free parameters · 1 axioms · 0 invented entities

Abstract-only review; the modeling step implicitly invokes standard continuous-time Markov chain theory on countable state spaces.

axioms (1)
  • domain assumption The chemical reaction network with sources and sinks can be modeled as a continuous-time Markov chain on the state space of molecule counts.
    Foundational modeling choice stated in the abstract.

pith-pipeline@v0.9.1-grok · 5629 in / 1136 out tokens · 28158 ms · 2026-07-03T06:30:49.829492+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Reference graph

Works this paper leans on

28 extracted references · 28 canonical work pages

  1. [1]

    Alberts, D

    B. Alberts, D. Bray, J. H. Wilson, J. Lewis, M. Raff, K. Roberts, J. D. Watson, et al.,Molecular biology of the cell 7th edition, vol. 2, 2022

    work page 2022

  2. [2]

    D. F. Anderson,Boundedness of trajectories for weakly reversible, single linkage class reaction sys- tems, Journal of Mathematical Chemistry, 49 (2011), pp. 2275–2290

    work page 2011

  3. [3]

    D. F. Anderson, D. Cappelletti, M. Koyama, and T. G. Kurtz,Non-explosivity of stochastically modeled reaction networks that are complex balanced, Bulletin of Mathematical Biology, 80 (2018), pp. 2561–2579

    work page 2018

  4. [4]

    D. F. Anderson, G. Craciun, and T. G. Kurtz,Product-form stationary distributions for deficiency zero chemical reaction networks, Bulletin of Mathematical Biology, 72 (2010), pp. 1947–1970

    work page 2010

  5. [5]

    D. F. Anderson and T. G. Kurtz,Continuous time Markov chain models for chemical reaction networks, in Design and analysis of biomolecular circuits: engineering approaches to systems and synthetic biology, Springer, 2011, pp. 3–42

    work page 2011

  6. [6]

    674, Springer, 2015

    ,Stochastic analysis of biochemical systems, vol. 674, Springer, 2015

    work page 2015

  7. [7]

    Avanzini, M

    F. Avanzini, M. Bilancioni, V. Cavina, S. Dal Cengio, M. Esposito, G. Falasco, D. Forastiere, J. N. Freitas, A. Garilli, P. E. Harunari, et al.,Methods and conversations in (post) modern thermody- namics, SciPost Physics Lecture Notes, (2024), p. 080

    work page 2024

  8. [8]

    Avanzini, M

    F. Avanzini, M. Esposito, and E. Penocchio,Coarse graining photo-isomerization reactions: Ther- modynamic consistency and implications for molecular ratchets, The Journal of Chemical Physics, 163 (2025)

    work page 2025

  9. [9]

    K. Ball, T. G. Kurtz, L. Popovic, and G. Rempala,Asymptotic analysis of multiscale approximations to reaction networks, The Annals of Applied Probability, 716 (2006), pp. 1925–1961

    work page 2006

  10. [10]

    Boyd,Detailed balance in chemical kinetics as a consequence of microscopic reversibility, The Journal of Chemical Physics, 60 (1974), pp

    R. Boyd,Detailed balance in chemical kinetics as a consequence of microscopic reversibility, The Journal of Chemical Physics, 60 (1974), pp. 1214–1222

    work page 1974

  11. [11]

    Cappelletti and C

    D. Cappelletti and C. Wiuf,Product-form Poisson-like distributions and complex balanced reaction systems, SIAM Journal on Applied Mathematics, 76 (2016), pp. 411–432

    work page 2016

  12. [12]

    Catacuzzeno, M

    L. Catacuzzeno, M. V. Leonardi, F. Franciolini, C. Domene, A. Michelucci, and S. Furini,Building predictive markov models of ion channel permeation from molecular dynamics simulations, Biophys- ical Journal, 123 (2024), pp. 3832–3843

    work page 2024

  13. [13]

    L. C. Evans,Partial differential equations, vol. 19, American Mathematical Society, 2022

    work page 2022

  14. [14]

    Feinberg,Complex balancing in general kinetic systems, Archive for Rational Mechanics and Analysis, 49 (1972), pp

    M. Feinberg,Complex balancing in general kinetic systems, Archive for Rational Mechanics and Analysis, 49 (1972), pp. 187–194

    work page 1972

  15. [15]

    ,Foundations of Chemical Reaction Network Theory, Springer International Publishing, 2019

    work page 2019

  16. [16]

    M. A. Ferreira, E. Franco, and J. J. Vel´ azquez,On the self-similar behavior of coagulation systems with injection, Annales de l’Institut Henri Poincar´ e C, 40 (2022), pp. 803–861

    work page 2022

  17. [17]

    Franco and J

    E. Franco and J. Vel´ azquez,The detailed balance property and chemical systems out of equilibrium, Bulletin of Mathematical Biology, 87 (2025), pp. 108–159

    work page 2025

  18. [18]

    Horn and R

    F. Horn and R. Jackson,General mass action kinetics, Archive for Rational Mechanics and Analysis, 47 (1972), pp. 81–116

    work page 1972

  19. [19]

    T. G. Kurtz,Strong approximation theorems for density dependent Markov chains, Stochastic Pro- cesses and their Applications, 6 (1978), pp. 223–240

    work page 1978

  20. [20]

    T. M. Liggett,Interacting particle systems, vol. 2, Springer, 1985

    work page 1985

  21. [21]

    Maas and A

    J. Maas and A. Mielke,Modeling of chemical reaction systems with detailed balance using gradient structures, Journal of Statistical Physics, 181 (2020), p. 2257. 34

    work page 2020

  22. [22]

    M. C. Mackey,Ion transport through biological membranes: an integrated theoretical approach, Springer Science & Business Media, 2013

    work page 2013

  23. [23]

    J. R. Norris,Markov chains, no. 2, Cambridge University Press, 1998

    work page 1998

  24. [24]

    Polettini and M

    M. Polettini and M. Esposito,Irreversible thermodynamics of open chemical networks. I. Emergent cycles and broken conservation laws, The Journal of Chemical Physics, 141 (2014)

    work page 2014

  25. [25]

    Rao and M

    R. Rao and M. Esposito,Conservation laws and work fluctuation relations in chemical reaction networks, The Journal of Chemical Physics, 149 (2018)

    work page 2018

  26. [26]

    Remlein, M

    B. Remlein, M. Esposito, and F. Avanzini,What is a chemostat? Insights from hybrid dynamics and stochastic thermodynamics, The Journal of Chemical Physics, 162 (2025)

    work page 2025

  27. [27]

    R. T. Rockafellar,Convex analysis, vol. 28, Princeton University Press, 1997

    work page 1997

  28. [28]

    Wegscheider, ¨Uber simultane gleichgewichte und die beziehungen zwischen thermodynamik und reaktionskinetik homogener systeme, Zeitschrift f¨ ur physikalische Chemie, 39 (1902), pp

    R. Wegscheider, ¨Uber simultane gleichgewichte und die beziehungen zwischen thermodynamik und reaktionskinetik homogener systeme, Zeitschrift f¨ ur physikalische Chemie, 39 (1902), pp. 257–303. 35

    work page 1902