Nonlinear Response Relations and Fluctuation-Response Inequalities for Nonequilibrium Stochastic Systems
Pith reviewed 2026-05-18 13:34 UTC · model grok-4.3
The pith
Nonlinear responses of observables in nonequilibrium stochastic systems are expressed as covariances with a Bell-polynomial conjugate variable derived from stochastic entropy production.
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
Our formalism expresses nonlinear responses of observables in terms of the covariance between the observable and a nonlinear conjugate variable. The nonlinear conjugate variable is subject to the complete Bell polynomial form and is determined by the stochastic entropy production. In addition, the Fluctuation-Response Inequalities (FRIs) are also derived for nonlinear responses, unraveling the general trade-off relations between nonlinear response and systems' fluctuations far from equilibrium.
What carries the argument
The nonlinear conjugate variable in complete Bell polynomial form, determined by the stochastic entropy production, which allows expressing nonlinear responses via covariances.
If this is right
- The framework unifies linear and nonlinear response theories for nonequilibrium systems.
- It provides general trade-off relations between nonlinear responses and fluctuations.
- The theory applies to systems like the symmetric exclusion process, as verified numerically.
- It can guide the design of sensitive and adaptive networks in synthetic biology and chemistry.
Where Pith is reading between the lines
- If true, this could extend response theory to more complex perturbations in biological signaling pathways.
- The inequalities might inspire new bounds for information processing in nonequilibrium steady states.
- Numerical checks on other models beyond SEP could test broader applicability.
Load-bearing premise
The system follows Markovian stochastic dynamics so that path probabilities and entropy production can be used to build the conjugate variable with Bell polynomials.
What would settle it
A numerical simulation or experiment on a Markov process such as the symmetric exclusion process where the measured nonlinear response deviates from the predicted covariance with the Bell-polynomial conjugate variable would falsify the claim.
Figures
read the original abstract
Predicting how systems respond to external perturbations far from equilibrium remains a fundamental challenge across physics, chemistry, and biology. We present a unified response framework for stochastic Markov dynamics that integrates linear and nonlinear perturbations. Our formalism expresses nonlinear responses of observables in terms of the covariance between the observable and a nonlinear conjugate variable. The nonlinear conjugate variable is subject to the complete Bell polynomial form and is determined by the stochastic entropy production. In addition, the Fluctuation-Response Inequalities (FRIs) are also derived for nonlinear responses, unraveling the general trade-off relations between nonlinear response and systems' fluctuations far from equilibrium. The validity of our theory is verified by the numerical results from a symmetric exclusion process (SEP). By unifying and extending nonequilibrium linear response theories, our approach can provide principled design rules for sensitive, adaptive synthetic and biological networks.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The manuscript develops a unified framework for linear and nonlinear response relations in nonequilibrium Markovian stochastic systems. It expresses the nonlinear response of an observable as a covariance between the observable and a nonlinear conjugate variable constructed via complete Bell polynomials applied to the stochastic entropy production. The work also derives fluctuation-response inequalities (FRIs) that bound nonlinear responses in terms of fluctuations and verifies the formalism numerically on a symmetric exclusion process (SEP).
Significance. If the central claims hold with the stated generality, the paper would provide a valuable extension of stochastic thermodynamics and linear response theory to the nonlinear regime. The Bell-polynomial construction offers a compact, cumulant-based representation that could yield falsifiable predictions and design principles for adaptive networks. The derivation appears parameter-free and grounded in path probabilities rather than ad-hoc fitting.
major comments (2)
- [Main derivation / SEP section] The central derivation (presumably in the main text following the abstract) assumes that the ratio of path probabilities under nonlinear perturbation factors exactly through the cumulants of the entropy production in complete Bell-polynomial form. For general state-dependent nonlinear driving (e.g., rates nonlinear in the control parameter), the response functional generally contains irreducible multi-time correlations that are not captured by a single entropy-production functional. The SEP verification uses symmetric exclusion dynamics whose rates admit this factorization; it does not test the general case. A counter-example or explicit assumption on the form of the perturbation should be added.
- [FRI derivation] The abstract states that FRIs are derived for nonlinear responses, but without the explicit inequality statement, the bounding functional, or the equality case, it is unclear whether the trade-off is tight or merely a loose bound. If the FRI reduces to a known linear-response result in the appropriate limit, this should be shown explicitly.
minor comments (2)
- [Numerical results] The numerical verification on the SEP should include error bars, system-size scaling, and a clear statement of the observable and perturbation protocol used.
- [Notation] Notation for the nonlinear conjugate variable and the Bell polynomials should be introduced with a short table or explicit definition to improve readability.
Simulated Author's Rebuttal
We thank the referee for the careful reading and constructive comments on our manuscript. We address each major comment point by point below, indicating the revisions we will incorporate.
read point-by-point responses
-
Referee: [Main derivation / SEP section] The central derivation (presumably in the main text following the abstract) assumes that the ratio of path probabilities under nonlinear perturbation factors exactly through the cumulants of the entropy production in complete Bell-polynomial form. For general state-dependent nonlinear driving (e.g., rates nonlinear in the control parameter), the response functional generally contains irreducible multi-time correlations that are not captured by a single entropy-production functional. The SEP verification uses symmetric exclusion dynamics whose rates admit this factorization; it does not test the general case. A counter-example or explicit assumption on the form of the perturbation should be added.
Authors: We agree that the central derivation relies on the path-probability ratio being expressible via the cumulants of the stochastic entropy production, which permits the complete Bell-polynomial representation. This structure holds for the class of Markovian systems and perturbations considered in the work, where the nonlinear driving enters the transition rates such that the entropy production along trajectories serves as the conjugate variable without additional irreducible multi-time correlations. The symmetric exclusion process is chosen precisely because its rates admit this factorization, providing a clean numerical test. To address the referee’s concern, we will revise the manuscript to state explicitly the assumption on the perturbation form (global nonlinear driving that preserves the Markovian property and allows the entropy-production functional to capture the response). We will also add a short discussion of the limitations for more general state-dependent nonlinearities that may introduce extra correlations. No counter-example is supplied because the framework is intended to apply within the stated assumptions; we believe this clarification resolves the generality issue without altering the core results. revision: yes
-
Referee: [FRI derivation] The abstract states that FRIs are derived for nonlinear responses, but without the explicit inequality statement, the bounding functional, or the equality case, it is unclear whether the trade-off is tight or merely a loose bound. If the FRI reduces to a known linear-response result in the appropriate limit, this should be shown explicitly.
Authors: We thank the referee for highlighting this presentational gap. The fluctuation-response inequalities are derived in the main text after the nonlinear response relations, but the abstract and the surrounding discussion do not spell out the explicit form. In the revised manuscript we will (i) update the abstract to include a concise statement of the FRI, (ii) present the explicit inequality together with the bounding functional expressed in terms of fluctuations and the equality case, and (iii) demonstrate that the nonlinear FRI reduces to the standard linear fluctuation-dissipation relation in the small-perturbation limit, thereby recovering a known result and confirming consistency. revision: yes
Circularity Check
No significant circularity; derivation self-contained from Markov path measures and entropy production
full rationale
The paper constructs the nonlinear conjugate variable via complete Bell polynomials applied to the stochastic entropy production, which is defined directly from the ratio of forward and backward path probabilities under the Markov assumption. This step follows from the standard definition of entropy production in stochastic thermodynamics and does not presuppose or fit the target nonlinear response of the observable. The covariance expression for the response is then obtained by differentiating the generating function, without reducing to a fitted parameter or a self-citation chain. Numerical verification on the symmetric exclusion process serves as external validation rather than a definitional input. No load-bearing self-citation, ansatz smuggling, or renaming of known results occurs in the central derivation.
Axiom & Free-Parameter Ledger
axioms (1)
- domain assumption The underlying dynamics are continuous-time Markov processes on a discrete state space.
Forward citations
Cited by 2 Pith papers
-
Mutual Linearity in and out of Stationarity for Markov Jump Processes: A Trajectory-Based Approach
A trajectory-level derivation shows mutual linearity holds for non-stationary Markov jump processes and generalizes to other systems.
-
Mutual Linearity in Nonequilibrium Langevin Dynamics
Local perturbations in nonequilibrium Langevin dynamics induce linear relations between stationary densities and currents at different positions due to an underlying one-dimensional response structure.
Reference graph
Works this paper leans on
-
[1]
Linear Conjugate Variable In general, a control parameterλcan impact the value of one or many transition rates{R ij}. From here and below, we useB n,ij to represent the corresponding con- jugate variables when choosingλ= lnR ij. By doing so, the linear conjugate variable Λ becomes thedynamical discrepancydefined for the corresponding edgej→i: Λij[Xτ]≡N ij...
-
[2]
First, by summing over all trajec- tories on both sides of Eq
Nonlinear Conjugate Variables The nonlinear conjugate variableB n shares similar properties with Λ. First, by summing over all trajec- tories on both sides of Eq. (6), we find that⟨B n⟩= 0. Then, by taking trajectory-ensemble average on both sides of the recurrence relation Eq. (8), we immediately get⟨B n+1⟩=⟨ΛB n⟩+⟨∂ λBn⟩= 0, which leads to ⟨ΛBn⟩=−⟨∂ λBn...
work page 2025
-
[3]
Linear Conjugate V ariable Let us takeQ= Λ ij, and consider the response ∂lnR kl ⟨Λij⟩. The response is given by ∂lnR kl ⟨Λij⟩=⟨∂ lnR klΛij⟩+ Cov(Λ ij,Λ kl) = 0.(A1) The response is always zero since⟨Λ ij⟩= 0 for any pa- rameter. Ifi=kandj=l, then Λ ij =N ij −R ijTj depends on lnR ij. Therefore, we have Var[Λij] = Cov(Λij,Λ ij) =−⟨∂ lnR ijΛij⟩=⟨R ijTj⟩=⟨N...
-
[4]
Ifi̸=korj̸=l, thenB n,ij is independent of lnR kl
Non-linear Conjugate V ariable Let us takeQ=B n,ij, and consider the response ∂n lnR kl ⟨Bn,ij⟩. Ifi̸=korj̸=l, thenB n,ij is independent of lnR kl. In this case, we have ∂n lnR kl ⟨Bn,ij⟩= Cov(B n,ij, Bn,kl) = 0,(A7) where⟨B n,ij⟩= 0 for all parameters. Ifi=kandj=l, thenB n,ij depends on lnR ij. In this case, we need to use the chain rule expansion. For c...
-
[5]
Kubo, Statistical-mechanical theory of irreversible processes
R. Kubo, Statistical-mechanical theory of irreversible processes. i. general theory and simple applications to magnetic and conduction problems, Journal of the phys- ical society of Japan12, 570 (1957)
work page 1957
-
[6]
Kubo, The fluctuation-dissipation theorem, Reports on progress in physics29, 255 (1966)
R. Kubo, The fluctuation-dissipation theorem, Reports on progress in physics29, 255 (1966)
work page 1966
-
[7]
G. S. Agarwal, Fluctuation-dissipation theorems for systems in non-thermal equilibrium and applications, Zeitschrift f¨ ur Physik A Hadrons and nuclei252, 25 (1972)
work page 1972
- [8]
- [9]
-
[10]
H. Feng and J. Wang, Potential and flux decomposi- tion for dynamical systems and non-equilibrium ther- modynamics: Curvature, gauge field, and generalized fluctuation-dissipation theorem, The Journal of chemical physics135(2011)
work page 2011
-
[11]
U. M. B. Marconi, A. Puglisi, L. Rondoni, and A. Vulpi- ani, Fluctuation–dissipation: response theory in statisti- cal physics, Physics reports461, 111 (2008)
work page 2008
-
[12]
G. Diezemann, Nonlinear response theory for markov processes: Simple models for glassy relaxation, Physical Review E85, 051502 (2012)
work page 2012
-
[13]
M. Kryvohuz and S. Mukamel, Nonlinear response theory in chemical kinetics, The Journal of chemical physics140 (2014)
work page 2014
-
[14]
S. Mukamel, V. Khidekel, and V. Chernyak, Classical chaos and fluctuation-dissipation relations for nonlinear response, Physical Review E53, R1 (1996)
work page 1996
-
[15]
R. L. Peterson, Formal theory of nonlinear response, Re- views of modern physics39, 69 (1967)
work page 1967
-
[16]
U. Seifert and T. Speck, Fluctuation-dissipation theorem in nonequilibrium steady states, Europhysics Letters89, 10007 (2010)
work page 2010
- [17]
-
[18]
M. Baiesi and C. Maes, An update on the nonequilib- rium linear response, New Journal of Physics15, 013004 (2013)
work page 2013
-
[19]
Maes, Response theory: a trajectory-based approach, Frontiers in Physics8, 229 (2020)
C. Maes, Response theory: a trajectory-based approach, Frontiers in Physics8, 229 (2020)
work page 2020
-
[20]
J. Zheng and Z. Lu, Universal response inequalities be- yond steady states via trajectory information geometry, Physical Review E112, L012103 (2025)
work page 2025
-
[21]
K. Ptaszy´ nski, T. Aslyamov, and M. Esposito, Dissipa- tion bounds precision of current response to kinetic per- turbations, Physical Review Letters133, 227101 (2024)
work page 2024
-
[22]
T. Aslyamov, K. Ptaszy´ nski, and M. Esposito, Nonequi- librium fluctuation-response relations: From identities to bounds, Physical Review Letters134, 157101 (2025)
work page 2025
-
[23]
K. Ptaszynski, T. Aslyamov, and M. Esposito, Nonequi- librium fluctuation-response relations for state observ- ables, arXiv preprint arXiv:2412.10233 (2024)
-
[24]
K. Ptaszynski, T. Aslyamov, and M. Esposito, Nonequi- librium fluctuation-response relations for state-current correlations, arXiv preprint arXiv:2506.08877 (2025)
-
[25]
E. Kwon, H.-M. Chun, H. Park, and J. S. Lee, Fluctuation-response inequalities for kinetic and entropic perturbations, arXiv preprint arXiv:2411.18108 (2024)
-
[26]
Unified Linear Fluctuation-Response Theory Arbitrarily Far from Equilibrium
J. Zheng and Z. Lu, Spatial correlation unifies nonequi- librium response theory for arbitrary markov jump pro- cesses, arXiv preprint arXiv:2501.01050 (2025)
work page internal anchor Pith review Pith/arXiv arXiv 2025
- [27]
-
[28]
J. A. Owen, T. R. Gingrich, and J. M. Horowitz, Univer- sal thermodynamic bounds on nonequilibrium response with biochemical applications, Physical Review X10, 011066 (2020)
work page 2020
-
[29]
G. Fernandes Martins and J. M. Horowitz, Topologically constrained fluctuations and thermodynamics regulate nonequilibrium response, Physical Review E108, 044113 (2023)
work page 2023
-
[30]
H.-M. Chun and J. M. Horowitz, Trade-offs between number fluctuations and response in nonequilibrium chemical reaction networks, The Journal of Chemical Physics158(2023)
work page 2023
-
[31]
J. A. Owen and J. M. Horowitz, Size limits the sensitiv- ity of kinetic schemes, Nature communications14, 1280 (2023)
work page 2023
-
[32]
P. E. Harunari, S. Dal Cengio, V. Lecomte, and M. Polet- tini, Mutual linearity of nonequilibrium network currents, Physical Review Letters133, 047401 (2024)
work page 2024
-
[33]
T. Aslyamov and M. Esposito, General theory of static response for markov jump processes, Physical Review Letters133, 107103 (2024)
work page 2024
-
[34]
T. Aslyamov and M. Esposito, Nonequilibrium response for markov jump processes: exact results and tight bounds, Physical Review Letters132, 037101 (2024)
work page 2024
-
[35]
U. Basu, M. Kr¨ uger, A. Lazarescu, and C. Maes, Fre- netic aspects of second order response, Physical Chem- istry Chemical Physics17, 6653 (2015)
work page 2015
-
[36]
Y. Hasegawa and T. Van Vu, Uncertainty relations in stochastic processes: An information inequality ap- proach, Physical Review E99, 062126 (2019)
work page 2019
-
[37]
A. Dechant and S.-i. Sasa, Fluctuation–response in- equality out of equilibrium, Proceedings of the National Academy of Sciences117, 6430 (2020)
work page 2020
- [38]
-
[39]
L´ evy, Sur certains processus stochastiques homog` enes, Compositio mathematica7, 283 (1940)
P. L´ evy, Sur certains processus stochastiques homog` enes, Compositio mathematica7, 283 (1940)
work page 1940
-
[40]
C. Godreche and J. Luck, Statistics of the occupation time of renewal processes, Journal of Statistical Physics 104, 489 (2001)
work page 2001
-
[41]
S. N. Majumdar and A. Comtet, Local and occupation time of a particle diffusing in a random medium, Physical review letters89, 060601 (2002)
work page 2002
-
[42]
A. Lapolla and A. Godec, Manifestations of projection- induced memory: General theory and the tilted single file, Frontiers in Physics7, 182 (2019)
work page 2019
-
[43]
D. Hartich and A. Godec, Violation of local detailed bal- ance upon lumping despite a clear timescale separation, Physical Review Research5, L032017 (2023). 9
work page 2023
-
[44]
Ramaswamy, The mechanics and statistics of active matter, Annu
S. Ramaswamy, The mechanics and statistics of active matter, Annu. Rev. Condens. Matter Phys.1, 323 (2010)
work page 2010
-
[45]
S. Ramaswamy, Active matter, Journal of Statistical Me- chanics: Theory and Experiment2017, 054002 (2017)
work page 2017
-
[46]
M. C. Marchetti, J.-F. Joanny, S. Ramaswamy, T. B. Liverpool, J. Prost, M. Rao, and R. A. Simha, Hydrody- namics of soft active matter, Reviews of modern physics 85, 1143 (2013)
work page 2013
-
[47]
O’Connor,Time-correlated single photon counting (Academic press, 2012)
D. O’Connor,Time-correlated single photon counting (Academic press, 2012)
work page 2012
-
[48]
G. Margolin and E. Barkai, Nonergodicity of blinking nanocrystals and other l´ evy-walk processes, Physical re- view letters94, 080601 (2005)
work page 2005
- [49]
-
[50]
I. V. Gopich and A. Szabo, Theory of the energy transfer efficiency and fluorescence lifetime distribution in single- molecule fret, Proceedings of the National Academy of Sciences109, 7747 (2012)
work page 2012
-
[51]
W. Bialek and S. Setayeshgar, Physical limits to biochem- ical signaling, Proceedings of the National Academy of Sciences102, 10040 (2005)
work page 2005
-
[52]
R. G. Endres and N. S. Wingreen, Accuracy of direct gra- dient sensing by single cells, Proceedings of the National Academy of Sciences105, 15749 (2008)
work page 2008
-
[53]
T. Mora and N. S. Wingreen, Limits of sensing tempo- ral concentration changes by single cells, Physical review letters104, 248101 (2010)
work page 2010
-
[54]
C. C. Govern and P. R. ten Wolde, Fundamental limits on sensing chemical concentrations with linear biochemical networks, Physical review letters109, 218103 (2012)
work page 2012
-
[55]
C. C. Govern and P. R. ten Wolde, Energy dissipation and noise correlations in biochemical sensing, Physical review letters113, 258102 (2014)
work page 2014
-
[56]
S. E. Harvey, S. Lahiri, and S. Ganguli, Universal energy- accuracy tradeoffs in nonequilibrium cellular sensing, Physical Review E108, 014403 (2023)
work page 2023
- [57]
-
[58]
W. A. Catterall, Ion channel voltage sensors: structure, function, and pathophysiology, Neuron67, 915 (2010)
work page 2010
-
[59]
B. Roux, T. Allen, S. Berneche, and W. Im, Theoreti- cal and computational models of biological ion channels, Quarterly reviews of biophysics37, 15 (2004)
work page 2004
-
[60]
E. T. Bell, Exponential polynomials, Annals of Mathe- matics35, 258 (1934)
work page 1934
-
[61]
U. Seifert, Stochastic thermodynamics, fluctuation the- orems and molecular machines, Reports on progress in physics75, 126001 (2012)
work page 2012
-
[62]
L. Peliti and S. Pigolotti,Stochastic thermodynamics: an introduction(Princeton University Press, 2021)
work page 2021
-
[63]
T. Hatano and S.-i. Sasa, Steady-state thermodynam- ics of langevin systems, Physical review letters86, 3463 (2001)
work page 2001
-
[64]
C. Maes and K. Netoˇ cn` y, A nonequilibrium extension of the clausius heat theorem, Journal of Statistical Physics 154, 188 (2014)
work page 2014
-
[65]
A. Dechant, S.-i. Sasa, and S. Ito, Geometric decomposi- tion of entropy production into excess, housekeeping, and coupling parts, Physical Review E106, 024125 (2022)
work page 2022
-
[66]
L. Torbjørn Stutzer, C. Dieball, and A. Godec, Stochas- tic calculus for pathwise observables of markov-jump pro- cesses: Unification of diffusion and jump dynamics, arXiv e-prints , arXiv (2025)
work page 2025
-
[67]
G. M. Sch¨ utz, Exact solution of the master equation for the asymmetric exclusion process, Journal of statistical physics88, 427 (1997)
work page 1997
-
[68]
T. M. Liggett, Ergodic theorems for the asymmetric simple exclusion process, Transactions of the American Mathematical Society213, 237 (1975)
work page 1975
-
[69]
C. Giardina, J. Kurchan, and L. Peliti, Direct evaluation of large-deviation functions, Physical review letters96, 120603 (2006)
work page 2006
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.