{"id":"01c2b15e-efb1-4fe1-b790-00100b5a6871","arxiv_id":"2601.01273","paper_version":3,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":1,"one_line_summary":"For slowly driven reversible Markov chains, the friction metric equals the graph-Laplacian pseudoinverse (resistance distance and commute-time geometry), and the thermodynamic distance is an equilibrium-restricted discrete Wasserstein transport cost.","lead":"Slowly driven random systems pay extra energy to shuffle probability around; this paper shows that dissipation cost is exactly the electrical resistance — equivalently the random-walk commute time — of a network built from the system's rates. The result links stochastic thermodynamics, circuit theory, and optimal transport, and turns hard friction-metric calculations into simple circuit algebra.","discovery_kind":"unification","skeptic_critique":{"model":"deepseek-v4-flash","headline":"NESS generalization in App. A rests on a false response identity: Eq. (A2) fails for a standard 3-state cycle with fixed affinities, so the claimed survival of Eq. (12) beyond detailed balance is unsupported.","rationale":"The paper's core discrete, detailed-balance theorem (Eq. 12 and the circuit/commute-time equivalences) is well supported, and the reader's verification is credible. The single most load-bearing concern is the NESS extension, because the introduction explicitly promises that the core geometric structures survive non-equilibrium steady states, and Appendix A is the only support. Eq. (A2) is not a minor qualifier: it is the identity that turns the actual NESS response into the detailed-balance form. The explicit 3-state counterexample shows the identity fails under standard local detailed balance with fixed affinities. This does not overturn the main discrete result, so the reader's CONDITIONAL verdict remains appropriate; the paper should restrict the NESS claims or correct the derivation. I did not select the continuous-space divergence as primary because it does not touch the discrete central claim, and I did not select the OT flux-function conditionality because it is already acknowledged as a formal restriction in App. B. The NESS issue is the one that silently invalidates an explicit generalization claim in the introduction.","tokens_in":14295,"tokens_out":34964,"duration_ms":294891,"concrete_test":"Recompute the stationary distribution and its response for the 3-state cycle above: solve Wπ=0 and the linearized equation W δπ = −(dW/dε)π with 1^T δπ=0 at ε=0. Compare δπ with ππ^T−Dπ. The analytic result deviates by roughly 0.065 to 0.13 in the components, so Eq. (A2) does not hold. This directly tests the key premise of Appendix A.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central discrete reversible result, Eq. (12), appears sound; the reader's verification holds. The load-bearing problem is the claimed NESS generalization in Appendix A. Eq. (A2) asserts that under conservative driving with fixed nonconservative forces, dπ/d(βV) = ππ^T − Dπ, i.e., the same stationary-distribution response as under detailed balance. This is false in general. Consider a 3-state cycle with local detailed balance and fixed affinity a=1 on each edge: rates k01=e^{1−ε}, k10=1, k12=1, k21=e^{−1−ε}, k20=1, k02=e^{−1}, with energies V=(0,ε,0). At ε=0 the stationary distribution is π≈(0.2447,0.3985,0.3568). Differentiating Wπ=0 w.r.t. ε gives dπ/dε≈(0.1628,−0.1748,0.0120), whereas Eq. (A2) predicts (0.0975,−0.2397,0.1421). The deviations are large and not a matter of convention. Consequently Eq. (A3) and the claimed survival of the βg∼−1/2C equivalence for arbitrary NESS do not follow. The introduction's statement that 'the core geometric structures derived here survive relaxation to NESS' is therefore unsupported. This is a concrete mathematical failure, not a disagreement with consensus.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper establishes an equivalence, for slowly driven continuous-time Markov chains in the linear-response regime, between the thermodynamic friction metric on the probability simplex and three graph-theoretic objects: the Moore–Penrose pseudoinverse of the weighted graph Laplacian, the effective resistance (via a resistor-network mapping), and the symmetrized mean first-passage (commute) time. The central chain is Eq. (12): βg ~ L^+ ~ −(1/2)R_eff ~ −(1/2)C. The paper also derives closed-form friction metrics for linear and cyclic graphs, interprets linear-response dissipation as Joule heating, connects the thermodynamic distance to a restricted discrete L^2-Wasserstein cost, and claims extensions to nonequilibrium steady states (Appendix A) and continuous state spaces (Section VII and Appendix C).","tokens_in":14510,"tokens_out":13070,"duration_ms":123500,"significance":"The discrete reversible result is valuable and appears sound: it unifies linear-response thermodynamics with classical random-walk and circuit theory, and it provides practical computational tools (circuit reduction) for friction metrics. The derivation builds on classical identities (deviation matrix vs. MFPTs, resistance distance vs. Laplacian pseudoinverse) and the reader's hand-checked verification of the 2- and 3-state examples supports the core equivalence. The resistor-network picture, including the exact formulas for linear and cyclic topologies, is a useful contribution. However, the paper's claims beyond the detailed-balance discrete setting are not supported: the NESS generalization rests on a false response identity, and the continuous-space commute-time kernel is not finite in d≥2. These overclaims need to be corrected before the paper can be recommended for publication.","major_comments":[{"comment":"The NESS generalization rests on Eq. (A2), ∂π/∂(βV) = ππ^T − Dπ, cited to Ref. [41] for fixed nonconservative forces. This identity is not valid for general irreducible rate matrices satisfying local detailed balance with fixed affinities. For example, take the 3-state cycle with k01=e^{1−ε}, k10=1, k12=1, k21=e^{−1−ε}, k20=1, k02=e^{−1} and energies V=(0,ε,0). Differentiating Wπ=0 at ε=0 gives dπ/dε≈(0.163,−0.175,0.012), whereas Eq. (A2) predicts (0.098,−0.240,0.142). Consequently Eqs. (A3)–(A4) and the claimed survival of Eq. (12) for NESS do not follow. The NESS statements in Section III and the Introduction must be removed or rederived under assumptions that actually guarantee the response identity.","section":"Appendix A, Eq. (A2)"},{"comment":"The continuous-space commute-time kernel C(x,y) as defined is not finite for reversible diffusions in d≥2. The heat-kernel short-time asymptotics give p_t(x|x)−p_t(x|y) ∼ (4πDt)^{-d/2} for t→0 when x≠y, so the t-integral diverges at t=0. Exponential relaxation of the generator (the 'standard assumptions' cited) controls only the large-t tail, not the short-time singularity. Thus the statement that 'the integral in (32) remains finite' is incorrect, and the proof of βg ∼ −(1/2)C for continuous spaces is not well-defined. The authors should restrict the continuous claim to d=1 (where the expression matches the classical MFPT) or introduce a regularized kernel and clearly state the resulting approximation.","section":"Sec. VII / App. C, Eqs. (32)/(C1)"}],"minor_comments":[{"comment":"The infimum in Eq. (15) is written only with the constraint ˙π_s = L_s φ_s. Please specify explicitly that the infimum is over equilibrium paths π_s connecting π_0 to π_1, i.e., π_0 and π_1 fixed and π_s ∈ Δ_n for all s, so the 'equilibrium-path-restricted' nature of the OT cost is unambiguous.","section":"Section IV, Eq. (15)"},{"comment":"Please check the physical dimensions in Eq. (18). The commute-time distance C(x,y) has units of time (the MFPTs are times), so ε^2 k_B T C(x,y) does not have units of energy unless an implicit rate/frequency scale is specified. If rescaled variables or dimensionless rates are used, state this directly.","section":"Section V A, Eq. (18)"},{"comment":"The equivalence symbol used in Eq. (12) and elsewhere is compact but can confuse: R_eff is defined as a two-point function, while the equivalence is on the tangent space of the simplex. For readability, explicitly define the quadratic-form equivalence for R_eff and C, e.g., ˙π^T βg ˙π = −(1/2) Σ_{x,y} ˙π(x) R_eff(x,y) ˙π(y), before using the shorthand.","section":"General notation"}],"recommendation":"major_revision","confidential_remarks":"The core discrete reversible equivalence (Eq. 12) is solid and publishable after revision. However, the paper currently overclaims in two areas that are explicitly advertised: the NESS survival (Appendix A) and the continuous-space generalization (Section VII). The false response identity in Eq. (A2) is a concrete mathematical error, and the divergence of the continuous commute-time kernel in d≥2 is a standard heat-kernel fact. I would suggest the editor require the authors to either remove these extensions or substantially rework them, and to recalibrate the abstract and introduction accordingly. The paper would then be a solid contribution to stochastic thermodynamics and graph theory."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The central result here is the equivalence chain (12) for reversible continuous-time Markov chains: the linear-response friction metric, the graph Laplacian pseudoinverse, effective resistance, and the symmetrized MFPT matrix all coincide as quadratic forms on the simplex. I checked the 2- and 3-state path examples myself; they close exactly, including the holding-time factors. The algebra is short but the synthesis is real: tying thermodynamic geometry to classical random-walk and circuit theory is a genuinely useful unification, and it gives practical tools — circuit reduction for friction metrics, and clean closed forms for linear and cyclic graphs. The Joule-heating picture and the bottleneck interpretation are helpful. The paper deserves serious referee time.\n\nSoft spots, in proportion. First, the continuous-space extension (Sec. VII and App. C) is not correct as stated for d≥2. The claim that the integral in (32) remains finite is false: for a diffusion in R^d, the heat kernel has a short-time divergence, p_t(x|x) − p_t(x|y) ~ t^{-d/2}, so the integral diverges for d≥2 unless a cutoff is introduced. The metric equivalence holds only in d=1, or with explicit regularization. That is a fixable qualifier, not a fatal blow.\n\nSecond, and more serious, the NESS generalization in Appendix A rests on Eq. (A2), which asserts the equilibrium response relation dπ/d(βV) = ππᵀ − D_π holds even with fixed nonconservative forces. That is false in general. I ran the 3-state counterexample from the stress-test: with fixed edge affinities and a simple energy perturbation, the actual stationary-distribution response differs from (A2) substantially. Consequently Eqs. (A3)–(A4) and the claimed survival of βg ~ −C/2 for arbitrary NESS do not follow. The introduction's statement that core structures survive NESS is unsupported. This is the load-bearing flaw in the appendix, but it does not affect the detailed-balance derivation that is the paper's main contribution.\n\nThe OT identification (15) is also stated more broadly than warranted: it holds only when the flux function is a generalized mean under the chosen rate law, as the appendix makes clear but the main text glosses over. That is a presentation issue.\n\nBottom line: for reversible processes, the main equivalence is proven, original, and useful. The NESS claims and the continuous-space claims need correction or removal before publication. A serious referee should engage; with those sections fixed, this would be a solid contribution.\n\nRecommendation: send to peer review, but be prepared for major revision on the NESS and continuous extensions.","headline":"The core reversible-chain result is sound and worth publishing after fixes; the NESS generalization and the continuous-space extension both overreach.","tokens_in":15146,"tokens_out":1657,"would_cite":true,"duration_ms":19604,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"For slowly driven Markov chains, the thermodynamic friction metric is the same quadratic form as the graph Laplacian pseudoinverse, effective resistance, and mean commute time; and the thermodynamic distance is a discrete optimal-transport","keywords":["thermodynamic geometry","friction metric","Markov chains","linear response","graph Laplacian","resistance distance","commute time","optimal transport"],"falsifier":"Numerically solve the master equation for a three-state Markov chain driven through a cycle by time-varying energies, compute the exact mean excess work, and check whether the leading term as protocol duration τ goes to infinity matches the circuit prediction βg ≃ L+; the equivalence should hold to first order in 1/τ, with deviations appearing only at higher order. For the non-equilibrium steady-state extension, drive a three-state cycle with a fixed non-conservative force and compare the ratio of linear-response excess work to the detailed-balanced case against the predicted factor α = (a0+a1","tokens_in":14000,"feed_emoji":"⚡","tokens_out":7052,"duration_ms":75916,"temperature":0.7,"pith_summary":"The paper claims that the thermodynamic friction metric governing energy dissipation in slowly driven continuous-time Markov chains is, as a quadratic form, exactly the pseudoinverse of the weighted graph Laplacian. Because that Laplacian also determines effective electrical resistance and mean commute (round-trip) times, linear-response dissipation in any such network is equivalent to Joule heating in a resistor circuit and to a Euclidean embedding in which distance between states is commute time. The same framework identifies the linear-response thermodynamic distance as a discrete L2-Wasserstein optimal-transport cost evaluated along paths of equilibrium distributions. If these equivalences hold, a wide class of thermodynamic metric calculations reduce to circuit algebra, and dissipation acquires a concrete physical picture as the cost of routing probability through dynamical bottlenecks.","feed_headline":"Markov-chain friction is graph resistance and commute time","feed_subtitle":"Linear-response dissipation in slowly driven networks is Joule heating, so friction metrics become simple circuit algebra.","key_machinery":"The object that carries the whole equivalence is the weighted graph Laplacian L = −W Dπ, whose off-diagonal entries are the negative equilibrium fluxes w(x|y)π(y) and whose diagonal entries are the total outgoing fluxes. Its pseudoinverse L+ acts as the inverse on the subspace of probability-conserving directions; the projectors I − π1ᵀ strip the non-physical null direction and make βg, L+, −½ Reff, and −½ C numerically interchangeable as metrics on the simplex. The same Laplacian supplies the continuity equation π̇ = Lφ used in the optimal-transport representation, with φ the velocity (node) potential and edge currents given by Ohm's law.","core_discovery":"The paper's central claim is that, for a reversible continuous-time Markov chain driven slowly by conservative changes to state energies, the thermodynamic friction metric on the space of equilibrium distributions is the same quadratic form as the pseudoinverse of the weighted graph Laplacian, the electrical effective resistance, and the symmetrized mean first-passage (commute) time: βg ≃ L+ ≃ −½ Reff ≃ −½ C on the probability simplex. The argument uses the fact that the edge weights of the Laplacian are precisely the equilibrium fluxes, and that the generalized inverse of the rate matrix can be expressed through mean first-passage times, so all four geometries share one object. The paper al","pith_inferences":["The correspondence suggests a practical route the paper does not develop: because commute times and effective resistances can be estimated from trajectory data and spectral methods, thermodynamic friction metrics could be inferred directly from simulations or experiments on molecular and biological networks without reconstructing the full energy landscape.","If the discrete optimal-transport structure extends beyond the equilibrium-restricted setting, minimal-dissipation protocols on finite graphs may be obtainable as geodesics in a discrete Wasserstein geometry, potentially leading to finite-graph analogues of counterdiabatic driving that go beyond linear response.","The resistor-network picture implies that entropic bottlenecks—sparse connectivity—impose a dissipation cost that no conservative control can eliminate; this could be tested by comparing the excess work of two graphs that have identical equilibrium energies and populations but differ only in how edges are wired.","The NESS result that stationary currents reduce dissipation regardless of orientation is suggestive for autonomous molecular machines: it hints that background nonequilibrium flows can assist externally driven probability transport at no additional linear-response cost, a prediction that could be checked in a three-state cycle experiment."],"forward_implications":["Friction metrics for arbitrary Markov graphs can be computed with standard circuit rules (series and parallel reduction, Kron reduction) instead of matrix inversions, and adding an edge always lowers the linear-response dissipation by a classical monotonicity property of effective resistance.","The linear-response work cost of moving a small probability mass between two states is kBT times the squared commute-time distance between them, so bottlenecks appear as large embedded distances that are genuinely costly to traverse.","The squared thermodynamic distance between equilibrium distributions is a discrete optimal-transport cost evaluated along equilibrium paths, linking stochastic thermodynamics to discrete Wasserstein geometry and providing variational principles for minimal-dissipation protocols.","For non-equilibrium steady states with fixed non-conservative forces, the equivalence survives with the metric given by the symmetric part of the generalized Laplacian pseudoinverse, and stationary currents reduce the linear-response excess work compared with a detailed-balanced system with the same edge traffic."],"fun_headline_variants":["Friction metric equals graph resistance and commute time","Markov-chain dissipation is Joule heating on graphs","One geometry links friction, resistance, and commute time","Graph resistance is the friction cost of probability transport","Friction = resistance = commute time on graphs"],"cache_read_input_tokens":2304,"weakest_assumption_plain":"The load-bearing premise is the linear-response approximation: the protocol must be slow enough that the lag between the actual and instantaneous equilibrium distributions is small and proportional to the driving speed, so that the excess work is exactly the quadratic form assumed; if this fails for fast or strong driving, the metric equivalence and the circuit and optimal-transport identities no longer describe the dissipation.","fun_headline_variants_meta":{"raw":{"variants":["Friction metric equals graph resistance and commute time","Markov-chain dissipation is Joule heating on graphs","One geometry links friction, resistance, and commute time","Graph resistance is the friction cost of probability transport","Friction = resistance = commute time on graphs"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.00082,"raw_usage":{"total_tokens":3385,"prompt_tokens":661,"completion_tokens":2724,"prompt_tokens_details":{"cached_tokens":256},"prompt_cache_hit_tokens":256,"prompt_cache_miss_tokens":405,"completion_tokens_details":{"reasoning_tokens":2652}},"tokens_in":405,"tokens_out":2724,"duration_ms":18577,"temperature":1.0,"reasoning_tokens":2652,"cache_read_input_tokens":256,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-03T12:56:27.178610+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Numerically solve the master equation for a three-state Markov chain driven through a cycle by time-varying energies, compute the exact mean excess work, and check whether the leading term as protocol duration τ goes to infinity matches the circuit prediction βg ≃ L+; the equivalence should hold to first order in 1/τ, with deviations appearing only at higher order. For the non-equilibrium steady-state extension, drive a three-state cycle with a fixed non-conservative force and compare the ratio of linear-response excess work to the detailed-balanced case against the predicted factor α = (a0+a1","supporting_citations":[],"review_version":1}