REVIEW 4 minor 1 cited by
Affine relationships between steady currents
T0 review · 0 major / 4 minor · reviewed 2026-08-11 · deepseek-v4-flash
Pith's one-line read Steady currents in finite Markov networks obey an exact affine relation when rates are perturbed, and the response coefficients are differences of mean first-passage times in the unperturbed process.
desk verdict A clean, short proof of the affine current relation with a genuinely new first-passage-time interpretation; the defects are textual, not mathematical. read the letter →
The pith
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
The reading
What carries the argument
The central object is the mean first-passage time $\tau_0(z,x)$ of the unperturbed Markov jump process, solved from the Poisson-type equation $L_0 \tau_0(z,x)=\delta_{xz}/\rho_0(x)-1$ for $z\neq x$, with $\tau_0(x,x)=0$. Combining that equation with stationarity of the perturbed process yields Eq. (III.1), showing the stationary distribution change $\rho(x)-\rho_0(x)$ is $\rho_0(x)$ times a sum of mean-first-passage-time differences weighted by perturbed currents. Feeding that into the definition of the steady current gives the affine relation and identifies the susceptibility $\lambda_0$ as mean-first-passage-time differences, or equivalently as a difference of excess-current quasipotentials satisfying a Poisson equation.
What would settle it
Take a small finite network, solve the stationary master equation before and after perturbing one edge exactly by matrix inversion, compute the right-hand side of Eq. (IV.2) from mean first-passage times of the unperturbed process, and check equality entry by entry; any mismatch beyond numerical precision would refute the formula. Alternatively, choose a perturbation that removes a bridge: the mean first-passage time between the two resulting components becomes infinite, so the stated formula cannot be evaluated, pinpointing the irreducible assumption.
Extended reading notes
Core claim
On a finite multigraph with transition rates $k_i(x,y)$, perturb some channels by $\Delta_i(x,y)$ while requiring irreducibility. For any channel whose own rates are unchanged, the paper proves Eq. (IV.2): $j_i(x,y)=j^0_i(x,y)+\sum_{z,z'}\lambda^0_i(zz',(x,y))\,j_\Delta(z,z')$, with $\lambda^0_i(zz',(x,y)) = \frac{1}{2}[\rho_0(x) k^0_i(x,y)(\tau_0(z,x)-\tau_0(z',x)) - \rho_0(y) k^0_i(y,x)(\tau_0(z,y)-\tau_0(z',y))]$. Here $\tau_0$ is the mean first-passage time in the unperturbed process, and $j_\Delta(z,z')$ is the current carried by the rate changes. The formula is exact for finite irreducible networks and does not rely on being close to equilibrium. The same object can be written as a difference of quasipotentials, and the appendix gives an equivalent Kirchhoff spanning-tree expression.
Load-bearing premise
The result assumes a finite state space and that the perturbed network remains a single connected component, so every state is still reachable from every other; if a perturbation disconnects the graph, mean first-passage times diverge and the formula is not defined.
Editorial extensions
If this is right
- The current on any unperturbed edge is an affine function of the perturbed-edge currents, with coefficients fixed by the original process alone.
- The coefficient $\lambda_0$ is a difference of mean first-passage times, so it can be interpreted as the response of the current to moving the initial condition from one vertex to the other.
- The relation is nonperturbative: no small-driving or linear-response approximation is needed.
- For several simultaneous perturbations the affine relation holds but is not additive; the susceptibility in that case is computed in the process where all perturbed edges are removed.
Reading between the lines
- Editorial: Because the coefficients live entirely in the unperturbed process, the relation can be turned into an estimation scheme: sample mean first-passage times from a single simulation of the reference dynamics and predict currents for any single-edge perturbation without re-solving the full network.
- Editorial: The nonadditivity in Section V implies that tuning several rates simultaneously is a genuinely joint problem; a control strategy built by adding single-edge responses can be off by an amount controlled by the cross terms $\lambda_0$ computed with both edges removed.
- Editorial: The irreducible assumption excludes bridge removal; a natural extension would treat absorbing or split networks, where the diverging mean first-passage times are replaced by escape-time or boundary-condition quantities, potentially covering open chemical reaction networks with input and output reservoirs.
Signed reviews
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The paper derives an exact, nonperturbative affine relation between steady currents in a finite-state continuous-time Markov jump process when transition rates are perturbed. Starting from the Poisson equation for mean first-passage times and an adjoint identity, Eq. (III.1) expresses the change in the stationary distribution in terms of the perturbing currents and mean first-passage time differences in the unperturbed process. Substituting this into the expression for an unperturbed channel current gives the main result (IV.2): the current is affine in the perturbing currents, with coefficients determined solely by the unperturbed process. The authors then specialize to the case of one added edge, give a quasipotential representation (V.5), extend the result to multiple perturbed edges, and illustrate the non-additivity of the multi-edge relation. An appendix supplies a Kirchhoff-tree representation of the susceptibility.
Significance. If correct, the result is a clean and exact generalization of the mutual linearity discovered by Harunari et al., with a shorter derivation and a new first-passage-time interpretation. The derivation is self-contained and requires no fitting, simulation, or back-reference to the target result; the extension to multiple simultaneous perturbations and the explicit demonstration that the affine relation is not additive are genuine additions. The graphical Kirchhoff representation in Appendix A is also useful for identifying which network features control current response. The main claims are mathematically sound within the stated irreducible-network setting, and the paper properly acknowledges the prior results in [1] and [9].
minor comments (4)
- [Section V, Eq. (V.9)] The displayed rate assignments are inconsistent: the edge k46 is assigned the values 2, 1, and 1 in the three lines of Eq. (V.9), and k65 and k54 are each assigned both 2 and 1. As a result, the numerical values in Eqs. (V.11)-(V.13), including j(4,6)=1/9, cannot be reproduced from the stated rates. Please correct the rate list and re-run the example so that the non-additivity illustration is verifiable.
- [Section V, Eq. (V.8)] The chain of equalities in Eq. (V.8) proves V0_xy(a)-V0_xy(b) = -lambda0(ab,(x,y)), not +lambda0 as printed; the final relation (V.5), which uses V0_xy(b)-V0_xy(a), is nevertheless correct. Please fix the sign in the displayed equality and in the sentence that follows it.
- [Section IV-V] There is a factor-of-two convention that should be stated explicitly: Eq. (IV.1) defines lambda0_i(zz',(x,y)) for an ordered pair (z,z'), while Eq. (V.3) defines the coefficient for an unordered added edge as twice the ordered-pair value, because both (a,b) and (b,a) contribute to the sum in (IV.2). Without an explicit remark, a reader applying (IV.2) to a single added edge may miss the factor of two.
- [Section V, Eq. (V.5)] There is a minor typo in Eq. (V.5): the argument should read (V0_xy(b)-V0_xy(a)) j(a,b), with a comma rather than a period in 'j(a.b)'.
Circularity Check
No significant circularity: Eq. (IV.2) is derived from standard Markov-chain identities and is not equivalent to its inputs by construction.
full rationale
The paper's derivation is self-contained. Section III proves the stationary-distribution change formula (III.1) from the stationarity condition L†ρ = 0, the identity L†0(ρ - ρ0) = Σ jΔ, and the standard Poisson equation (III.5) for mean first-passage times. Section IV then substitutes (III.1) into the expression for the current over an unperturbed channel, yielding the affine relation (IV.2) with the susceptibility λ0 expressed entirely in the unperturbed process. No parameter is fitted to data, and no target result is assumed in the derivation. The quantities jΔ are the actual perturbed currents, not fitted proxies; Eq. (IV.2) is an exact interrelation between currents rather than a closed prediction of jΔ from independent inputs, so it does not reduce to its own definition. The acknowledgments that Eq. (III.1) coincides with relation (44) of Harvey et al. and that the original mutual linearity was discovered by Harunari et al. are explicit statements of prior context, not load-bearing self-citations. The one self-citation, Ref. [10], supports the standard first-passage Poisson equation and is an independent, textbook-level result. The stated assumption that irreducibility is preserved is a scoping condition, not a hidden definitional shortcut. The numerical example V.1 contains rate-assignment inconsistencies and a sign slip in the displayed Eq. (V.8), but these are presentation errors that do not affect the derivation of (IV.2) and are not circularity. Overall, the central claim is internally derived and no circular step can be exhibited.
Assumptions & free parameters
assumptions (5)
- domain assumption The system is a continuous-time Markov jump process on a finite irreducible multigraph with unique stationary distribution.
- domain assumption Perturbations Δ(x,y) keep the vertex set fixed and preserve irreducibility.
- standard math Mean first-passage times satisfy L0 τ0(z,x) = δ_{xz}/ρ0(x)-1 for z≠x, with τ0(x,x)=0.
- standard math The adjoint identity Σ_z L0 f(z) g(z) = Σ_z f(z) L†0 g(z) holds.
- standard math Kirchhoff's spanning-tree formula gives stationary distribution and currents.
Cite this review
Pith. "Pith review of Affine relationships between steady currents." pith.science (2026). https://pith.science/paper/MSVLYW3C
@misc{pith2026241205019,
author = {Pith},
title = {Pith review of: Affine relationships between steady currents},
year = {2026},
howpublished = {\url{https://pith.science/paper/MSVLYW3C}},
note = {Machine review of arXiv:2412.05019}
}
read the original abstract
Perturbing transition rates in a steady nonequilibrium system, e.g. modelled by a Markov jump process, causes a change in the local currents. Their susceptibility is usually expressed via Green-Kubo relations or their nonequilibrium extensions. However, we may also wish to directly express the mutual relation between currents. Such a nonperturbative interrelation was discovered by P.E. Harunari et al. in [1] by applying algebraic graph theory showing the mutual linearity of currents over different edges in a graph. We give a novel and shorter derivation of that current relationship where we express the current-current susceptibility as a difference in mean first-passage times. It allows an extension to multiple currents, which remains affine but the relation is not additive.
Figures
Forward citations
Cited by 1 Pith paper
-
Exact First-Passage Time Response Theory from Steady-State Response
Mean first-passage time response to arbitrary single-rate perturbations is expressed exactly through unperturbed MFPTs and steady-state probabilities, via a fast-reset correspondence.
Reference graph
Works this paper leans on
-
[1]
P.E. Harunari, S. Dal Cengio, V. Lecomte, and M. Polettini. Mutual linearity of nonequilib- rium network currents. Phys. Rev. Lett., 133:047401, 2024
work page 2024
- [9]
-
[2]
M. S. Green. Markoff random processes and the statistical mechanics of time-dependent phenomena. ii. irreversible processes in fluids. J. Chem. Phys., 22(3):398–413, 1954
work page 1954
-
[3]
R. Kubo. Statistical-mechanical theory of irreversible processes. i. general theory and simple applications to magnetic and conduction problems. J. Phys. Soc. Jpn., 12(6):570–586, 1957
work page 1957
-
[4]
L. Onsager. Reciprocal relations in irreversible processes. i. Phys. Rev., 37(4):405–426, 1931
work page 1931
-
[5]
C. Maes. Response theory: a trajectory-based approach. Front. Phys., 2020
work page 2020
-
[6]
P. Gaspard. The Statistical Mechanics of Irreversible Phenomena. Cambridge University Press, 2022
work page 2022
-
[7]
Coplanarity of rooted spanning-tree vectors
M. Polettini, P. E. Harunari, S. Dal Cengio, and V. Lecomte. Coplanarity of rooted spanning- tree vectors. arXiv 2407.16093 [math.CO], 2024. 14
work page Pith review arXiv 2024
Show all 11 references
-
[8]
S. Redner. A Guide to First-Passage Processes. Cambridge University Press, 2001
2001
-
[10]
Khodabandehlou, C
F. Khodabandehlou, C. Maes, and K. Netoˆ cn´ y. On the Poisson equation for nonreversible Markov jump processes. J. Math. Phys., 65(4), April 2024
2024
-
[11]
Khodabandehlou, C
F. Khodabandehlou, C. Maes, and K. Netoˇ cn´ y. Trees and forests for nonequilibrium purposes: An introduction to graphical representations. J. Stat. Phys, 189:41, 2022
2022
Reviewed August 11, 2026 · model on record in the stance chip above.
Discussion (0). Continue with ORCID to comment.