Recognition: unknown
Spectral Fluctuation-Dissipation-Response Inequalities
Pith reviewed 2026-05-09 23:40 UTC · model grok-4.3
The pith
Finite-state Markov jump processes obey spectral inequalities that bound how far the fluctuation-dissipation theorem can break in driven steady states.
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
By comparing the causal susceptibility to its passive equilibrium reference, frequency-resolved and frequency-integrated inequalities are established that bound their mismatch in terms of the steady-state entropy production rate, probe variance, short-time perturbation diffusion, and reversible relaxation timescales. The bounds exactly recover the standard fluctuation-dissipation theorem at equilibrium and apply directly to measurable causal susceptibilities.
What carries the argument
The passive equilibrium reference constructed from the steady-state probabilities, which acts as the baseline against which the causal susceptibility mismatch is bounded.
If this is right
- The inequalities supply experimentally testable upper limits on fluctuation-dissipation theorem violations in any driven steady state.
- At equilibrium the bounds collapse exactly to the classical fluctuation-dissipation theorem.
- The frequency-integrated version supplies a global constraint linking total mismatch to entropy production and relaxation times.
- The bounds involve only steady-state probabilities and measurable response functions, so they can be checked without reconstructing the full transition matrix.
Where Pith is reading between the lines
- The same comparison technique might be applied to continuous-state diffusion processes or to quantum master equations to obtain analogous spectral bounds.
- In biological or chemical reaction networks the inequalities could be inverted to estimate entropy production from measured response functions alone.
- Frequency resolution of the bounds may allow experimenters to identify which dynamical timescales carry the dominant dissipative cost.
Load-bearing premise
The passive equilibrium reference correctly isolates the reversible part of the dynamics for the observable being measured.
What would settle it
Measure the causal susceptibility spectrum and the steady-state entropy production rate in the same driven Markov jump process, then check whether the observed mismatch stays inside the predicted bound for all frequencies.
Figures
read the original abstract
We derive spectral fluctuation--dissipation--response inequalities for finite-state Markov jump processes. By comparing the causal susceptibility to its passive equilibrium reference, we establish frequency-resolved and frequency-integrated inequalities that bound their mismatch in terms of the steady-state entropy production rate, probe variance, short-time perturbation diffusion, and reversible relaxation timescales. Our bounds exactly recover the standard fluctuation--dissipation theorem at equilibrium and apply directly to measurable causal susceptibilities, providing experimentally testable thermodynamic limits on FDT breakdown in driven steady states.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The manuscript derives spectral fluctuation-dissipation-response inequalities for finite-state Markov jump processes. By comparing the causal susceptibility to its passive equilibrium reference (constructed from the steady-state distribution p_ss with detailed balance imposed), it establishes frequency-resolved and frequency-integrated inequalities bounding the mismatch in terms of the steady-state entropy production rate, probe variance, short-time perturbation diffusion, and reversible relaxation timescales. The bounds recover the standard fluctuation-dissipation theorem exactly at equilibrium and are asserted to apply directly to measurable causal susceptibilities, yielding experimentally testable thermodynamic limits on FDT breakdown in driven steady states.
Significance. If the central derivation holds, the work would supply concrete, frequency-dependent bounds linking linear response violations to entropy production and other dynamical quantities in non-equilibrium Markovian systems. This could be useful for interpreting experiments on driven steady states and for quantifying how far the FDT can be violated in terms of measurable irreversibility. The exact recovery of the equilibrium case and the focus on finite-state processes (allowing potential numerical checks) are strengths. The emphasis on directly measurable susceptibilities increases potential impact in statistical mechanics and non-equilibrium thermodynamics.
major comments (1)
- [Derivation of passive equilibrium reference] The construction of the passive equilibrium reference (detailed in the derivation following the abstract and in the section introducing the reference susceptibility): the claim that imposing detailed balance on the same p_ss exactly isolates the reversible component of the linear response for arbitrary observables is load-bearing for all subsequent bounds. If the chosen observable couples to irreversible modes or if the reference does not match the reversible part of the generator, the mismatch bound becomes uncontrolled. The manuscript must provide an explicit proof or counterexample verification that this identification holds rigorously for the observables considered; without it the inequalities rest on an unverified step.
minor comments (2)
- [Abstract] The abstract and introduction would benefit from a brief statement of the precise class of observables to which the inequalities apply (e.g., whether they must be functions of the state or can be more general).
- [Introduction] Notation for the causal susceptibility and the passive reference should be introduced with a clear equation early in the text to avoid ambiguity when comparing them.
Simulated Author's Rebuttal
We thank the referee for their positive evaluation of the manuscript's potential impact and for the careful identification of the key foundational step in our derivation. We address the single major comment below and have revised the manuscript accordingly to strengthen the presentation.
read point-by-point responses
-
Referee: [Derivation of passive equilibrium reference] The construction of the passive equilibrium reference (detailed in the derivation following the abstract and in the section introducing the reference susceptibility): the claim that imposing detailed balance on the same p_ss exactly isolates the reversible component of the linear response for arbitrary observables is load-bearing for all subsequent bounds. If the chosen observable couples to irreversible modes or if the reference does not match the reversible part of the generator, the mismatch bound becomes uncontrolled. The manuscript must provide an explicit proof or counterexample verification that this identification holds rigorously for the observables considered; without it the inequalities rest on an unverified step.
Authors: We agree that an explicit verification of this step strengthens the manuscript. The passive reference is constructed by symmetrizing the rate matrix W to W_rev while preserving the identical steady-state p_ss, which enforces detailed balance and removes the irreversible probability currents. Because the linear-response operator acts on the observable through the same perturbation, the reference susceptibility χ_ref(ω) exactly captures the reversible contribution to the full causal susceptibility χ(ω). The difference χ(ω) − χ_ref(ω) is then bounded by the entropy-production term and the reversible relaxation spectrum, as already indicated in the frequency-resolved inequality. To make this rigorous, we have added a new Appendix A containing (i) the generator decomposition W = W_rev + W_irr, (ii) the resulting decomposition of the response function for arbitrary observables, and (iii) a short numerical check on a driven two-state system where the observable couples to both reversible and irreversible modes; the bound remains valid and saturates at equilibrium. These additions directly address the referee’s concern without altering the main results. revision: yes
Circularity Check
No significant circularity; derivation is self-contained
full rationale
The abstract describes deriving spectral inequalities by comparing causal susceptibility to a passive equilibrium reference constructed from steady-state probabilities, with bounds expressed in terms of entropy production, probe variance, diffusion, and relaxation timescales. This recovers the equilibrium FDT when entropy production is zero. No equations are provided that define the reference using the mismatch quantity itself or that reduce the bound to a tautology by construction. The approach follows standard non-equilibrium thermodynamic constructions that isolate reversible dynamics via detailed balance on p_ss, without self-citation chains or fitted inputs renamed as predictions. The central claim remains independent of its outputs.
Axiom & Free-Parameter Ledger
Forward citations
Cited by 4 Pith papers
-
Dynamical Fluctuation-Response Relations
Exact dynamical fluctuation-response relations are derived that split the finite-time covariance of time-integrated observables into initial variability and an integral of response kernels for nonautonomous Markov jum...
-
Finite-frequency fluctuation-response bounds for open quantum systems
A finite-frequency fluctuation-response inequality bounds the measured lock-in response-to-noise matrix by the output-field quantum Fisher information rate for Markovian open quantum systems.
-
Dynamical Fluctuation-Response Relations
Exact dynamical fluctuation-response relations are derived for time-integrated observables of nonautonomous Markov jump processes, with covariance split into initial variability and response kernel integral.
-
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]
Kubo, Reports on Progress in Physics29, 255 (1966)
R. Kubo, Reports on Progress in Physics29, 255 (1966)
1966
-
[2]
Baiesi, C
M. Baiesi, C. Maes, and B. Wynants, Physical Review Letters103, 010602 (2009)
2009
-
[3]
Seifert, Reports on Progress in Physics75, 126001 (2012)
U. Seifert, Reports on Progress in Physics75, 126001 (2012)
2012
-
[4]
T. Speck and U. Seifert, Europhysics Letters74, 391 (2006), arXiv:cond-mat/0511696 [cond-mat.stat-mech]
-
[5]
Prost, J.-F
J. Prost, J.-F. Joanny, and J. M. R. Parrondo, Physical Review Letters103, 090601 (2009)
2009
-
[6]
Harada and S.-i
T. Harada and S.-i. Sasa, Physical Review Letters95, 130602 (2005)
2005
-
[7]
Wang, Physical Review E97, 052125 (2018)
S.-W. Wang, Physical Review E97, 052125 (2018)
2018
-
[8]
A. C. Barato and U. Seifert, Physical Review Letters 114, 158101 (2015)
2015
-
[9]
T. R. Gingrich, J. M. Horowitz, N. Perunov, and J. L. England, Physical review letters116, 120601 (2016)
2016
-
[10]
Dechant and S.-i
A. Dechant and S.-i. Sasa, Proceedings of the National Academy of Sciences117, 6430 (2020)
2020
-
[11]
J. A. Owen, T. R. Gingrich, and J. M. Horowitz, Phys. Rev. X10, 011066 (2020)
2020
-
[12]
J. A. Owen and J. M. Horowitz, Nat. Commun.14, 1280 (2023)
2023
-
[13]
Fernandes Martins and J
G. Fernandes Martins and J. M. Horowitz, Phys. Rev. E 108, 044113 (2023)
2023
-
[14]
Aslyamov and M
T. Aslyamov and M. Esposito, Phys. Rev. Lett.132, 037101 (2024)
2024
-
[15]
Aslyamov and M
T. Aslyamov and M. Esposito, Phys. Rev. Lett.133, 107103 (2024)
2024
-
[16]
Zheng and Z
J. Zheng and Z. Lu, Phys. Rev. E112, L012103 (2025)
2025
-
[17]
Liu and J
K. Liu and J. Gu, Communications Physics8, 62 (2025)
2025
-
[18]
Kwon, H.-M
E. Kwon, H.-M. Chun, H. Park, and J. S. Lee, Physical Review Letters135, 097101 (2025)
2025
-
[19]
Van Vu, PRX Quantum6, 010343 (2025)
T. Van Vu, PRX Quantum6, 010343 (2025)
2025
- [20]
- [21]
-
[22]
Thermodynamic constraints on the power spectral density in and out of equilibrium,
A. Dechant, “Thermodynamic constraints on the power spectral density in and out of equilibrium,” (2023), arXiv:2306.00417 [cond-mat.stat-mech]
-
[23]
A. Dechant, Phys. Rev. Lett. (2026), 10.1103/3hs9-dz3d, accepted paper, accepted March 24, 2026
-
[24]
J. Zheng and Z. Lu, “Thermodynamic and ki- netic bounds for finite-frequency fluctuation-response,” (2026), arXiv:2602.18631 [cond-mat.stat-mech]
-
[25]
Aslyamov, K
T. Aslyamov, K. Ptaszy´ nski, and M. Esposito, Phys. Rev. Lett.136, 067102 (2026)
2026
-
[26]
Maes and B
C. Maes and B. Wynants, Markov Processes and Related Fields16, 45 (2010)
2010
-
[27]
Shiraishi, K
N. Shiraishi, K. Funo, and K. Saito, Phys. Rev. Lett. 121, 070601 (2018)
2018
-
[28]
Falasco and M
G. Falasco and M. Esposito, Physical Review Letters 125, 120604 (2020)
2020
-
[29]
N. Ohga, S. Ito, and A. Kolchinsky, Physical Review Letters131, 077101 (2023)
2023
-
[30]
Liang and S
S. Liang and S. Pigolotti, Physical Review E108, L062101 (2023)
2023
-
[31]
Shiraishi, Physical Review E108, L042103 (2023)
N. Shiraishi, Physical Review E108, L042103 (2023)
2023
-
[32]
Van Vu, V
T. Van Vu, V. T. Vo, and K. Saito, Physical Review Research6, 013273 (2024)
2024
-
[33]
Gu, Physical Review E109, L042101 (2024)
J. Gu, Physical Review E109, L042101 (2024)
2024
-
[34]
J. R. Gomez-Solano, A. Petrosyan, S. Ciliberto, R. Chetrite, and K. Gawedzki, Physical Review Letters 103, 040601 (2009)
2009
-
[35]
P. Jop, A. Petrosyan, and S. Ciliberto, EPL (Euro- physics Letters)81, 50005 (2008)
2008
-
[36]
J. V. Koski, T. Sagawa, O.-P. Saira, Y. Yoon, A. Kutvo- nen, P. Solinas, M. M¨ ott¨ onen, T. Ala-Nissila, and J. P. Pekola, Nature Physics9, 644 (2013)
2013
-
[37]
Toyabe, T
S. Toyabe, T. Okamoto, T. Watanabe-Nakayama, H. Taketani, S. Kudo, and E. Muneyuki, Physical Re- view Letters104, 198103 (2010)
2010
-
[38]
E. Oikonomou, Y. Juli, R. R. Kolan, L. Kern, T. Gru- ber, C. Alzheimer, P. Krauss, A. Maier, and T. Huth, Communications Chemistry (2024), 10.1038/s42004- 024-01369-y, article via DOI. 6
-
[39]
R. Roy, S. Hohng, and T. Ha, Nature Methods5, 507 (2008)
2008
-
[40]
Goldbeter and D
A. Goldbeter and D. E. J. Koshland, Proc. Natl. Acad. Sci. USA78, 6840 (1981)
1981
-
[41]
J. R. Szemere, H. G. Rotstein, and A. C. Ventura, npj Systems Biology and Applications7, 32 (2021)
2021
-
[42]
Bakry, I
D. Bakry, I. Gentil, and M. Ledoux,Analysis and Geometry of Markov Diffusion Operators, Grundlehren der mathematischen Wissenschaften, Vol. 348 (Springer, Cham, 2014)
2014
-
[43]
G. A. Pavliotis,Stochastic Processes and Applications: Diffusion Processes, the Fokker–Planck and Langevin Equations, Texts in Applied Mathematics, Vol. 60 (Springer, New York, 2014)
2014
-
[44]
Villani,Hypocoercivity, Memoirs of the American Mathematical Society, Vol
C. Villani,Hypocoercivity, Memoirs of the American Mathematical Society, Vol. 202 (American Mathemati- cal Society, Providence, RI, 2009)
2009
-
[45]
Kawai, J
R. Kawai, J. M. R. Parrondo, and C. van den Broeck, Physical Review Letters98, 080602 (2007)
2007
-
[46]
Rold´ an and J
´E. Rold´ an and J. M. R. Parrondo, Physical Review E85, 031129 (2012)
2012
-
[47]
G. S. Agarwal, Zeitschrift f¨ ur Physik252, 25 (1972)
1972
-
[48]
Seifert and T
U. Seifert and T. Speck, EPL (Europhysics Letters)89, 10007 (2010)
2010
-
[49]
Shiraishi, Journal of Statistical Physics185, 19 (2021)
N. Shiraishi, Journal of Statistical Physics185, 19 (2021)
2021
-
[50]
Dechant, J
A. Dechant, J. Phys. A: Math. Theor.55, 094001. END MA TTER Recall the stationary fluxesa ij =π jWij, net cur- rentsJ ij =a ij −a ji, Jji =−J ij, symmetric activi- tiesA ij =a ij +a ji, and thermodynamic forcesF ij = lna ij/aji, Fji =−F ij. The steady–state entropy produc- tion rate isσ= 1 2 P i,j JijFij ≥0 with equality if and only if detailed balance ho...
-
[51]
SinceW asym is antisymmetric in theπ-inner product and annihilates constants, one has ⟨v⟩π = 2ζβ⟨1, W asymb⟩π =−2ζβ⟨W asym1, b⟩π = 0
BASIC REPRESENT A TION AND CENTEREDNESS The mismatch is itself a causal susceptibility. SinceW asym is antisymmetric in theπ-inner product and annihilates constants, one has ⟨v⟩π = 2ζβ⟨1, W asymb⟩π =−2ζβ⟨W asym1, b⟩π = 0. Thus bothrandvbelong to the centered subspace H0 ≡ f: Ω→R ⟨f⟩ π = 0 . LetW † ⊥ denote the restriction ofW † toH 0. Its spectrum lies st...
-
[52]
In particular, Crv ∈L 1([0,∞))∩L 2([0,∞))
ANAL YTICITY, SYMMETR Y, AND V ANISHING CRITERIA Since the state space is finite,C rv(t) is a finite linear combination of exponentially decaying terms multiplied by finite polynomials int. In particular, Crv ∈L 1([0,∞))∩L 2([0,∞)). Therefore ∆χ(ω) is analytic forℑω >0 and obeys the standard Kramers–Kronig relations for causal response functions [1], ℜ∆χ(...
-
[53]
Therefore ∆χadmits a convergent Taylor expansion nearω= 0, ∆χ(ω) = ∞X n=0 (iω)n n! Mn.(14) 3 The first terms are ∆χ(ω) =M 0 +iωM 1 − ω2 2 M2 +O(ω 3)
LOW-FREQUENCY EXP ANSION BecauseC rv(t) decays exponentially, all time moments exist: Mn ≡ Z ∞ 0 dt tnCrv(t)<∞, n= 0,1,2, . . . . Therefore ∆χadmits a convergent Taylor expansion nearω= 0, ∆χ(ω) = ∞X n=0 (iω)n n! Mn.(14) 3 The first terms are ∆χ(ω) =M 0 +iωM 1 − ω2 2 M2 +O(ω 3). Hence ℜ∆χ(ω) =M 0 − ω2 2 M2 +O(ω 4),ℑ∆χ(ω) =ωM 1 +O(ω 3). The zero-frequency ...
-
[54]
The first three terms are ∆χ(ω) = i⟨rv⟩π ω − ⟨r, W †v⟩π ω2 − i⟨r,(W †)2v⟩π ω3 +O(|ω| −4)
HIGH-FREQUENCY ASYMPTOTICS Repeated integration by parts gives the asymptotic expansion ∆χ(ω) = m−1X n=0 ⟨r,(W †)nv⟩π (−iω)n+1 +O(|ω| −m−1),|ω| → ∞. The first three terms are ∆χ(ω) = i⟨rv⟩π ω − ⟨r, W †v⟩π ω2 − i⟨r,(W †)2v⟩π ω3 +O(|ω| −4). In particular, |∆χ(ω)|=O(|ω| −1). Thus the mismatch necessarily decays at high frequency. The leading 1/ωtail is deter...
-
[55]
NORM IDENTITIES AND GENERIC BOUNDS Let K(t)≡Θ(t)C rv(t), so that ∆χ(ω) is the Fourier transform ofK. Parseval then yields the exact identity 1 2π Z ∞ −∞ dω|∆χ(ω)| 2 = Z ∞ 0 dt|C rv(t)|2.(15) Equation (15) makes clear that the integrated mismatch is the squaredL 2 norm of the causal violation kernel. Before invoking thermodynamic closures, one already has ...
-
[56]
Consequently, |∆χ(ω)|2 =O(σ), 1 2π Z dω|∆χ(ω)| 2 =O(σ)
WEAK-DRIVING SCALING If all thermodynamic forces are scaled as Fij =ϵf ij, ϵ≪1, at fixed symmetric activitiesA ij, then Jij =A ij tanh(Fij/2) =O(ϵ), σ= 1 2 X i,j JijFij =O(ϵ 2), and by (13), v=O(ϵ),∆χ(ω) =O(ϵ). Consequently, |∆χ(ω)|2 =O(σ), 1 2π Z dω|∆χ(ω)| 2 =O(σ). Thus the inequalities are asymptotically sharp in scaling near equilibrium. The current–tr...
-
[57]
For a wavenumber q= 2πm N , m= 1,
RING CASE: EXACT MODE ANAL YSIS AND NEAR SA TURA TION Consider anN-state ring, with indices understood moduloN, and uniform clockwise and counterclockwise rates Wi+1,i =k e F/2 , W i−1,i =k e −F/2 .(19) The stationary distribution is uniform, πi = 1 N . For a wavenumber q= 2πm N , m= 1, . . . , N−1 2 , introduce the normalized real Fourier pair cq(i) = √ ...
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.