pith. sign in

arxiv: 2604.18528 · v1 · submitted 2026-04-20 · ❄️ cond-mat.stat-mech

Emergent nonreciprocity in open thermodynamically-consistent chemical reaction networks

Daniel Evans , Yizhi Shen , Ahmad K. Omar This is my paper

Pith reviewed 2026-05-10 03:47 UTC · model grok-4.3

classification ❄️ cond-mat.stat-mech
keywords chemical reaction networksnonreciprocityOnsager reciprocitychemostatsnonequilibrium steady statesoscillatory instabilitiesreaction-diffusion systems
0
0 comments X

The pith

Topology of open chemical reaction networks induces nonreciprocity leading to oscillatory instabilities.

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

This paper establishes that the structure of open chemical reaction networks, when connected to chemostats, can break Onsager reciprocity and produce complex eigenvalues in the linearized dynamics around steady states. The result is oscillatory instabilities even though the networks remain thermodynamically consistent. The local equilibrium hypothesis ensures that the dissipative dynamics stay variational, so the system continues to minimize a free-energy-like quantity. Readers should care because this provides a mechanism for nonequilibrium patterns such as oscillations to emerge purely from network topology without additional assumptions about activity or driving.

Core claim

The topology of open, thermodynamically-consistent chemical reaction networks can result in oscillatory instabilities near nonequilibrium steady states. These instabilities arise from chemostat-induced breaking of Onsager reciprocity, while the local equilibrium hypothesis preserves the variational structure of the dissipative part of the dynamics. Numerical results confirm that such nonreciprocity in reaction-diffusion systems produces oscillatory dynamics that nevertheless minimize a free energy.

What carries the argument

Chemostat-induced breaking of Onsager reciprocity through network topology, which generates an asymmetric Jacobian while preserving variational dissipation under local equilibrium.

If this is right

  • Nonreciprocity can appear in thermodynamically consistent systems and generate oscillatory dynamics.
  • Reaction-diffusion systems with such networks exhibit oscillations that still minimize free energy.
  • Instabilities arise near nonequilibrium steady states due to topology alone.
  • The variational structure of dissipation is preserved despite broken reciprocity.

Where Pith is reading between the lines

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

  • Similar topological mechanisms could be used to design chemical oscillators in synthetic biology.
  • This suggests that nonreciprocity in other open systems like metabolic networks might lead to unexpected temporal behaviors.
  • Extensions to larger networks could reveal how specific topologies control the frequency of oscillations.

Load-bearing premise

The local equilibrium hypothesis continues to hold and maintain the variational character of dissipation even after chemostats break reciprocity in the network.

What would settle it

A counterexample where a specific open CRN topology with chemostats produces no oscillatory instability despite the predicted breaking of Onsager reciprocity, or where the dynamics fail to minimize the free energy during oscillations.

Figures

Figures reproduced from arXiv: 2604.18528 by Ahmad K. Omar, Daniel Evans, Yizhi Shen.

Figure 1
Figure 1. Figure 1: FIG. 1. Schematic of the steady state in an open system whose [PITH_FULL_IMAGE:figures/full_fig_p003_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: FIG. 2. Time evolution of (a) the dimensionless relative densities and [PITH_FULL_IMAGE:figures/full_fig_p004_2.png] view at source ↗
read the original abstract

Nonreciprocity, a hallmark of nonequilibrium systems, can generate dynamics not possible near thermodynamic equilibrium, including oscillatory and rotating patterns. The onset of temporal oscillations is often evident in linearized dynamics, where nonreciprocity appears as complex eigenvalues of an asymmetric Jacobian. Here, we show that the topology of open, thermodynamically-consistent chemical reaction networks can result in oscillatory instabilities near nonequilibrium steady states. These instabilities arise from chemostat-induced breaking of Onsager reciprocity, while the local equilibrium hypothesis preserves the variational structure of the dissipative part of the dynamics. Numerical results confirm that such nonreciprocity in reaction-diffusion systems produces oscillatory dynamics that nevertheless minimize a free energy.

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

2 major / 2 minor

Summary. The paper claims that the topology of open, thermodynamically consistent chemical reaction networks permits chemostat-induced breaking of Onsager reciprocity, which produces an asymmetric Jacobian and oscillatory instabilities near nonequilibrium steady states. The local-equilibrium hypothesis is invoked to ensure that the dissipative sector remains a gradient flow of a free-energy functional, so that the resulting oscillations still minimize free energy. Numerical simulations of reaction-diffusion systems are presented as confirmation.

Significance. If the asserted separation between nonreciprocal circulatory terms and a variational dissipative sector holds, the work supplies a concrete topological route to nonreciprocity inside thermodynamically consistent CRNs. This could explain or predict oscillatory behavior in open chemical and biological systems while preserving a free-energy minimization principle for the dissipative part. The numerical evidence for oscillations that nevertheless descend a free-energy landscape is a concrete, falsifiable prediction.

major comments (2)
  1. [§3.2] §3.2 (linearized dynamics around the NESS): the argument that chemostat boundary conditions break reciprocity only in the circulatory sector while leaving the friction matrix symmetric relies on an unexpanded form of the mass-action rates. An explicit expansion of the effective Jacobian under fixed external concentrations is needed to confirm that off-diagonal nonreciprocity does not enter the dissipative block.
  2. [Eq. (12)] Eq. (12) (variational structure): the claim that the dissipative dynamics remain a pure gradient flow under the local-equilibrium hypothesis is asserted after chemostatting; the proof that the effective mobility matrix stays symmetric (or symmetrizable) once external concentrations are fixed is not shown explicitly and is load-bearing for the central claim.
minor comments (2)
  1. [Figure 2] Figure 2 caption: the color scale for the free-energy landscape is not labeled, making it difficult to verify that trajectories descend the functional.
  2. [Introduction] The relation between the network topology and the specific chemostat placements that induce complex eigenvalues should be stated more explicitly in the introduction for readers unfamiliar with CRN stoichiometry.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for their careful reading of the manuscript and for highlighting these important points on the explicit structure of the linearized dynamics and the variational properties. We address each major comment below and will incorporate the requested clarifications in the revised version.

read point-by-point responses

  1. Referee: [§3.2] §3.2 (linearized dynamics around the NESS): the argument that chemostat boundary conditions break reciprocity only in the circulatory sector while leaving the friction matrix symmetric relies on an unexpanded form of the mass-action rates. An explicit expansion of the effective Jacobian under fixed external concentrations is needed to confirm that off-diagonal nonreciprocity does not enter the dissipative block.

    Authors: We agree that an explicit expansion strengthens the argument. In the revised manuscript we will linearize the mass-action rates around the NESS with fixed chemostat concentrations, writing the effective Jacobian as J = J_diss + J_circ. The dissipative block J_diss derives from the symmetric Onsager matrix contracted with the stoichiometric structure and remains symmetric; the chemostat terms appear only in the antisymmetric circulatory sector. This expansion will be added to §3.2 together with the resulting eigenvalue analysis. revision: yes

  2. Referee: [Eq. (12)] Eq. (12) (variational structure): the claim that the dissipative dynamics remain a pure gradient flow under the local-equilibrium hypothesis is asserted after chemostatting; the proof that the effective mobility matrix stays symmetric (or symmetrizable) once external concentrations are fixed is not shown explicitly and is load-bearing for the central claim.

    Authors: We acknowledge that the symmetry of the effective mobility matrix after chemostatting was stated rather than derived in detail. We will add an explicit derivation (either in the main text or as an appendix) showing that, under the local-equilibrium hypothesis, the mobility matrix L remains symmetric because the chemostat concentrations enter only through the thermodynamic forces, not through the kinetic coefficients. Consequently the dissipative sector continues to be a gradient flow of the free-energy functional, preserving Eq. (12). revision: yes

Circularity Check

0 steps flagged

No circularity: derivation grounded in topology and standard local-equilibrium assumption

full rationale

The paper asserts that chemostat-induced breaking of Onsager reciprocity produces an asymmetric Jacobian while the dissipative sector remains a gradient flow, with the separation justified by the local-equilibrium hypothesis. This is a standard modeling assumption in nonequilibrium thermodynamics rather than a self-definition or fitted input. No equations reduce the target result to a prior fit or self-citation chain by construction, and the numerical confirmation is presented as independent verification. The derivation chain is therefore self-contained against external thermodynamic benchmarks.

Axiom & Free-Parameter Ledger

0 free parameters · 2 axioms · 0 invented entities

The central claim rests on thermodynamic consistency of the networks and the local equilibrium hypothesis; no free parameters or invented entities are mentioned in the abstract.

axioms (2)
  • domain assumption Local equilibrium hypothesis
    Preserves the variational structure of the dissipative part of the dynamics.
  • domain assumption Thermodynamic consistency of open chemical reaction networks
    Ensures the system respects nonequilibrium thermodynamics when driven by chemostats.

pith-pipeline@v0.9.0 · 5411 in / 1304 out tokens · 38976 ms · 2026-05-10T03:47:51.445091+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

15 extracted references · 15 canonical work pages

  1. [1]

    & Esposito, M.J

    Polettini, M. & Esposito, M.J. Chem. Phys.141(2014)

    work page 2014

  2. [2]

    & Esposito, M.Phys

    Rao, R. & Esposito, M.Phys. Rev. X6, 041064 (2016). 7

    work page 2016

  3. [3]

    Avanzini, F., Aslyamov, T., Étienne Fodor & Esposito, M.J. Chem. Phys.161(2024)

    work page 2024

  4. [4]

    & Nelson, D

    Hatano, N. & Nelson, D. R.Phys. Rev. Lett.77, 570–573 (1996)

    work page 1996

  5. [5]

    & Nelson, D

    Hatano, N. & Nelson, D. R.Phys. Rev. B56, 8651–8673 (1997)

    work page 1997

  6. [6]

    & Chen, Y .-F.Advances in Physics: X7, 2109431 (2022)

    Zhang, X., Zhang, T., Lu, M.-H. & Chen, Y .-F.Advances in Physics: X7, 2109431 (2022)

    work page 2022

  7. [7]

    & Vaikuntanathan, S.Nat

    Murugan, A. & Vaikuntanathan, S.Nat. Commun.8, 13881 (2017)

    work page 2017

  8. [8]

    & Golestanian, R.Phys

    Tang, E., Agudo-Canalejo, J. & Golestanian, R.Phys. Rev. X11, 031015 (2021)

    work page 2021

  9. [9]

    Ezawa, M.Phys. Rev. B105, 125421 (2022)

    work page 2022

  10. [10]

    & Sato, M.Annual Review of Condensed Matter Physics14, 83–107 (2023)

    Okuma, N. & Sato, M.Annual Review of Condensed Matter Physics14, 83–107 (2023)

    work page 2023

  11. [11]

    & Tang, E.Phys

    Nelson, A. & Tang, E.Phys. Rev. B110, 155116 (2024)

    work page 2024

  12. [12]

    & Tang, E.Nat

    Zheng, C. & Tang, E.Nat. Commun.15, 6453 (2024)

    work page 2024

  13. [13]

    Nature627, 528–533 (2024)

    Veenstra, J.et al. Nature627, 528–533 (2024)

    work page 2024

  14. [14]

    Belyansky, R., Weis, C., Hanai, R., Littlewood, P. B. & Clerk, A. A.Phys. Rev. Lett.135, 123401 (2025)

    work page 2025

  15. [15]

    Poggioli, A. R. & Limmer, D. T.Phys. Rev. Lett.130, 158201 (2023). 8

    work page 2023