Thermodynamic and Kinetic Bounds for Finite-frequency Fluctuation-Response

Abstract: Fluctuation-response relations encode fundamental constraints on nonequilibrium systems. While time-domain static response is bounded by activity and entropy production, finite-frequency extensions for time-dependent perturbations remain largely unexplored. Here, we derive frequency-domain fluctuation-response inequalities for steady-state Markov processes with time-dependent perturbations. For barrier and entropic perturbations, the spectral signal-to-noise ratio (SNR) is universally bounded by dynamical activity. Furthermore, for state-current observables, the SNR is bounded by the entropy production rate (EPR). We illustrate our results using the F1-ATPase model to infer EPR. These finite-frequency inequalities provide a practical route to infer dissipation from power spectra measurements.

• The paper derives frequency-domain fluctuation-response inequalities for steady-state Markov jump processes, using spectral-matrix positivity to bound response spectra by dynamical activity.
• The paper shows that state-current responses under barrier perturbations obey a thermodynamic bound proportional to half the entropy production rate, providing a nonequilibrium constraint that vanishes at equilibrium.
• The paper demonstrates with an F1-ATPase model that intermediate-frequency power spectra can produce tighter entropy-production estimates than zero-frequency, time-domain measurements.

Overview

This paper derives fluctuation-response inequalities (FRIs) in the frequency domain for steady-state Markov jump processes subject to time-dependent perturbations, and shows that these yield both kinetic bounds (in terms of dynamical activity) and thermodynamic bounds (in terms of entropy production rate, EPR) on finite-frequency signal-to-noise ratios (SNRs). The work extends a body of time-domain results—where linear-response SNRs were bounded by activity and dissipation—to experimentally relevant spectral quantities such as power spectra of currents. The practical payoff is a route to infer steady-state EPR directly from measured power spectra, illustrated on a coarse-grained F1-ATPase model.

Setup: Markov jump processes under time-dependent perturbations

The system is an $n$-state reversible Markov jump process with local detailed balance, so that edge currents $j_{ij} = r_{ij}\pi_j - r_{ji}\pi_i$, edge traffic $a_{ij} = r_{ij}\pi_j + r_{ji}\pi_i$, total dynamical activity $a = \sum_{i<j} a_{ij}$, and total EPR $$\dot\sigma = \sum_{i<j} j_{ij}\ln(r_{ij}\pi_j / r_{ji}\pi_i)$$ are all well defined. Transition rates are parameterized as

$$r_{ij} = \exp\left[b_{ij}(\zeta) + f_{ij}(\xi)/2\right],$$

with symmetric barrier parameters $b_{ij}$ (tunable by catalysts, magnetic fields, or nano-electronic techniques) and antisymmetric entropic-force parameters $f_{ij}$. Perturbations are weak, time-dependent modulations $$\lambda \mapsto \lambda + \varepsilon\phi_\lambda(t)$$ of any of these parameters. Observables considered are state-current type functionals combining dwelling-time increments and antisymmetrically weighted jump counts; their fluctuations are characterized by the spectral density $\mathcal{S}(\omega)$, the Fourier transform of the autocovariance of $j_{ij} = r_{ij}\pi_j - r_{ji}\pi_i$0.

Finite-frequency fluctuation-response inequalities

The derivation starts from the path-probability representation of the response kernel, $j_{ij} = r_{ij}\pi_j - r_{ji}\pi_i$1, where $j_{ij} = r_{ij}\pi_j - r_{ji}\pi_i$2 is the impulse observable obtained by differentiating the log path probability with respect to the perturbation pattern. Building the joint spectral density matrix of $j_{ij} = r_{ij}\pi_j - r_{ji}\pi_i$3 and invoking its Hermiticity and positive definiteness via the Schur complement recovers the finite-frequency FRI of Dechant et al., which the authors then specialize to single-edge perturbations. A key structural result is that for barrier perturbations the spectral matrix of the impulse observables is diagonal and frequency-independent, $j_{ij} = r_{ij}\pi_j - r_{ji}\pi_i$4, while for entropic perturbations it equals $j_{ij} = r_{ij}\pi_j - r_{ji}\pi_i$5. This yields the first main result:

$j_{ij} = r_{ij}\pi_j - r_{ji}\pi_i$6

For state-current observables, an exact ratio identity holds,

$j_{ij} = r_{ij}\pi_j - r_{ji}\pi_i$7

which produces alternative FRIs weighted by squared currents rather than activities. Numerical simulations on a fully connected three-state network with random rates confirm saturation efficiency $j_{ij} = r_{ij}\pi_j - r_{ji}\pi_i$8 across frequencies, with $j_{ij} = r_{ij}\pi_j - r_{ji}\pi_i$9 decaying to zero at high frequency and flattening at low frequency where $a_{ij} = r_{ij}\pi_j + r_{ji}\pi_i$0 converges to twice the diffusion coefficient $a_{ij} = r_{ij}\pi_j + r_{ji}\pi_i$1 of the accumulated observable.

The frequency independence of $a_{ij} = r_{ij}\pi_j + r_{ji}\pi_i$2 is the structural reason these bounds extend nontrivially from static to time-dependent perturbations: it follows from the delta-correlated nature of Markovian noise combined with steady-state time-translation symmetry.

Kinetic and thermodynamic bounds on spectral SNR

Applying the chain rule and Cauchy–Schwarz inequality to the FRIs yields the second main result. For perturbations mediated by global parameters $a_{ij} = r_{ij}\pi_j + r_{ji}\pi_i$3 (barriers) and $a_{ij} = r_{ij}\pi_j + r_{ji}\pi_i$4 (entropic forces), the spectral SNR satisfies kinetic bounds controlled solely by dynamical activity:

$a_{ij} = r_{ij}\pi_j + r_{ji}\pi_i$5

For barrier perturbations acting on state-current observables, the current-weighted FRI gives a thermodynamic bound involving the pseudo-EPR $a_{ij} = r_{ij}\pi_j + r_{ji}\pi_i$6, which is itself bounded above by the true EPR via the log-mean inequality:

$a_{ij} = r_{ij}\pi_j + r_{ji}\pi_i$7

This is the central thermodynamic statement of the paper: unlike the activity bound, which applies equally at and away from equilibrium, the EPR bound vanishes at equilibrium and therefore captures a genuine nonequilibrium signature accessible through finite-frequency response measurements. The appendix extends analogous bounds to overdamped Langevin dynamics with multiplicative noise (anti-Ito discretization for thermodynamic consistency); there, perturbing $a_{ij} = r_{ij}\pi_j + r_{ji}\pi_i$8 yields $a_{ij} = r_{ij}\pi_j + r_{ji}\pi_i$9, whose spatial integral is exactly the steady-state EPR, so the Langevin thermodynamic bound carries over verbatim to finite frequency.

Frequency-integrated inequalities

Since integrating $a = \sum_{i<j} a_{ij}$0 over all frequencies recovers the time-domain variance, $a = \sum_{i<j} a_{ij}$1, the same chain-rule arguments applied to the integrated SNR' (noise measured by $a = \sum_{i<j} a_{ij}$2, often easier to obtain experimentally than a full spectrum) give matching integrated bounds by activity and by $a = \sum_{i<j} a_{ij}$3. These imply a conservation-like trade-off: for fixed activity or dissipation, amplifying the response-to-noise ratio in one frequency band necessarily suppresses it elsewhere.

Application to F1-ATPase

The authors instantiate the thermodynamic bound on a minimal three-state continuous-time Markov model of the F1-ATPase rotary motor (empty, ATP-bound, ADP-bound states), using experimentally determined energy differences ($a = \sum_{i<j} a_{ij}$4, $a = \sum_{i<j} a_{ij}$5, $a = \sum_{i<j} a_{ij}$6), driving forces ($a = \sum_{i<j} a_{ij}$7, $a = \sum_{i<j} a_{ij}$8, $a = \sum_{i<j} a_{ij}$9), and attempt rates $$\dot\sigma = \sum_{i<j} j_{ij}\ln(r_{ij}\pi_j / r_{ji}\pi_i)$$0. Taking the counterclockwise net probability current as the observable and comparing heterogeneous, uniform, and local Jacobian patterns $$\dot\sigma = \sum_{i<j} j_{ij}\ln(r_{ij}\pi_j / r_{ji}\pi_i)$$1, they find that the perturbation pattern affects only the magnitude of the response, not the validity of the bound. Notably, at intermediate frequencies the entropy-production bound on the spectral SNR can be tighter than the corresponding zero-frequency (time-domain) bound, suggesting that frequency-domain measurements may yield sharper EPR estimates than previously reported time-domain relations.

Limitations and open questions

The derivation rests on three assumptions stated explicitly by the authors: steady-state time-translation invariance, the Markov property, and linearity of response. Generalization to nonstationary and non-Markovian systems remains open, as does extension to nonlinear finite-frequency response. The thermodynamic bound additionally requires state-current observables, since it relies on the response-ratio identity specific to that class; whether an EPR-type bound exists for general trajectory observables is not settled here. Finally, extension to macroscopic and open quantum systems—where time-domain analogues have recently appeared—is left as future work.

Conclusion

The paper establishes rigorous finite-frequency fluctuation-response inequalities for steady-state Markov jump processes and overdamped Langevin systems, showing that spectral SNRs under barrier and entropic perturbations are bounded by dynamical activity, and that barrier-perturbation responses of state-current observables are bounded by half the entropy production rate. The diagonal, frequency-independent structure of the impulse-observable spectrum under Markovian noise is what makes the generalization from static to time-dependent perturbations possible. The F1-ATPase demonstration indicates that power-spectra-based dissipation inference can outperform time-domain estimators at intermediate frequencies, providing a concrete experimental protocol grounded in stochastic thermodynamics.