Equilibrium-like behavior in far-from-equilibrium chemical reaction networks

In an equilibrium chemical reaction mixture, the number of molecules present obeys a Poisson distribution. We ask when the same is true of the steady state of a nonequilibrium reaction network and obtain an essentially complete answer. In particular, we show that networks with certain topological features must have a Poisson distribution, whatever the reaction rates. Such driven systems also obey an analog of the fluctuation-dissipation theorem. Our results may be relevant to biological systems and to the larger question of how equilibrium concepts might apply to nonequilibrium systems.

💡 Research Summary

The paper addresses a fundamental question: under what circumstances does a nonequilibrium chemical reaction network exhibit the same Poisson statistics for molecular copy numbers that are characteristic of equilibrium mixtures? By combining graph‑theoretic analysis of reaction networks with the master‑equation formalism, the authors derive an essentially complete answer.

First, the authors formalize a reaction network as a directed graph whose vertices are complexes (linear combinations of species) and whose edges are elementary reactions. Assuming mass‑action kinetics, the stochastic dynamics are governed by a master equation for the probability distribution (P(\mathbf{n},t)) over the vector of molecule numbers (\mathbf{n}). In equilibrium, detailed balance guarantees that each species follows an independent Poisson distribution, but this property is not automatically inherited by driven steady states.

The central result is that topology alone can enforce Poisson statistics, irrespective of the numerical values of the rate constants. Specifically, if the network is (i) weakly reversible—every reaction lies on a directed cycle, so each connected component is strongly linked—and (ii) has deficiency zero (δ = |C| – ℓ – s = 0, where |C| is the number of complexes, ℓ the number of linkage classes, and s the dimension of the stoichiometric subspace), then the network is guaranteed to be complex‑balanced. Complex balance means that for each complex the total inflow rate equals the total outflow rate in the steady state. Under these conditions the stationary solution of the master equation is exactly a multivariate Poisson distribution:

\