Multicyclic Norias: A First-Transition Approach to Extreme Values of the Currents

For continuous-time Markov chains we prove that, depending on the notion of effective affinity F, the probability of an edge current to ever become negative is either 1 if \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$F< 0$$\end{document}F<0 else \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\sim \exp - F$$\end{document}∼exp-F. The result generalizes a “noria” formula to multicyclic networks. We give operational insights on the effective affinity and compare several estimators, arguing that stopping problems may be more accurate in assessing the nonequilibrium nature of a system according to a local observer. Finally we elaborate on the similarity with the Boltzmann formula. The results are based on a constructive first-transition approach.


Introduction
Let us first provide an illustrative simple example of some of our main results.Perform a random walk on the network where every time you cross one specific edge in one direction you get 1e, and when you cross it in the opposite direction you pay 1e.All other transitions also give and take credit, but in liras.Initially your pockets are empty of euros, and while you have a virtually infinite reservoir of liras your bank would not convert them into euros.The question we pose and answer in this paper is: assuming that you live forever, what is the probability f − that you eventually get broke (get to -1e, that you cannot actually afford to pay)?
In the special case where you are initially in ⋆ (so that at least you get a chance to not get immediately broke!), we find, on the assumption F > 0, with where each diagram implies multiplication by the corresponding rates of the Markov chain.If instead F < 0, the probability to go bankrupt is 1.
In the case where the diagonal transition has vanishing rates in both directions, in the last expression the diagrams containing diagonal terms disappear and we are left with the log-ratio of two cyclic contributions This object is known as cycle affinity, a measure of the probability of performing the cycle in one direction relative to the opposite.For unicyclic networks ("norias") Eq. ( 1) boils down to a remarkable formula first obtained by Bauer and Cornu [1] (which in fact is slightly more general in that the initial state can be chosen arbitarily).As reasonable, the cycle is completed more often in the favourable direction (A > 0) rather than the unfavourable one; then the Bauer-Cornu formula yields the probability of the rare event of the cycle to ever be completed more often in the unfavourable direction.Generalizations have been proven for the total entropy production (a weighted balance of euros and liras) using the tools of martingale theory [2][3][4] 1 , and discussed in the light of first-exit time problems [5].
Beyond the probabilistic interpretation, cycle affinities also afford two different thermodynamic interpretations, one global and one local, of which we give here an intuitive sketch (but see Sec. 3.2 for details), on the assumption that forward and backward rates of a transition x ↔ y have ratio exp δq xy /T xy with δq the heat (in our example: currency) exchange along that transition and T the temperature (in our example this could be the interest rate) of that transition, i.e. a measure of how inconvenient it is to perform that transition (to borrow money).Then the global one is as Carnot's entropy production along a cyclic process [6] where is a shorthand for the sum over cyclic transitions.This takes into account all contributions, e.g. from liras and from euros, even if for the problem at hand liras to not play any role.Thus this latter interpretation has little operational value.The local one is in terms the specific transaction in euros, its temperature T e , and the value of such temperature T ∅ at which forward and backward euro transactions happen with the same frequency (that is: the rest of the world does not care much of whatever you do with your euros).Going back to our main result Eq. ( 1) (which, we remind, holds for arbitrary non-unicyclic networks), here the effective affinity does not have such a simple global interpretation, but still retains the local interpretation.What is lost with respect to norias is the independence from the initial state: while it does not matter where you are initially along a cycle to complete that cycle, it does matter where you are in a network to perform an arbitrary composition of cycles that touch the initial state.We generalize Eq. ( 1) appropriately.
We then show by computational examples that first-passage and extreme value problems such as the one above might give a better estimate of the effective affinity than do fluctuation and fluctuation-dissipation relations at fixed stopping time.Finally, Eq. ( 1) is reminiscent of Boltzmann's formula, but turned upside down.We linger on this analogy towards the conclusions.

State-space processes
Consider an irreducible continuous-time Markov chain on finite state space X .We characterize it in terms of the probability p X t = {p X t (x), x ∈ X } of being at x at time t, which satisfies the continuous-time master equation starting from a given initial distribution p X 0 , with non-negative rates r(x|x ′ ) of jumping from x ′ to x.We also define the continuous-time (adjoint) generator R with matrix entries where δ is Kroenecker's delta and is the exit rate out of a state.The master equation reads in vector form d dt p X t = Rp X t and its stationary distribution solves Rp X ∞ = 0. From now on we do not specify the range of summation unless necessary.
We now focus on a pair of connected states, namely x, x ′ = 1, 2 without loss of generality.We assume that edge 1 ↔ 2 is not a bridge, that is, that its removal does not disconnect the system, and denote X ∅ (or simply ∅) a system where edge 1 ↔ 2 is removed.Transitions between these states are deemed to be visible to an external observer.Let ℓ ∈ L = {→ = 2 → 1, ← = 1 → 2} denote transitions between such states, to and from, and t(•), s(•) ∈ {1, 2} the source and target states of a transition, i.e. t(→) = s(←) = 1 and t(←) = s(→) = 2.
Letting n(ℓ) be the number of times transition ℓ occurs along a realization of the process, we define the visible activity and the cumulated current respectively as Notice that they typically grow linearly in time; thus we denote the mean stationary current (i.e.cumulated current per unit time) as Our final goal is to compute the probability f ± that the cumulated current c takes value ±1 at least once as the process unfolds from time t = 0 to time t → +∞.

Transition-state processes
Our strategy is to lift the description of the process from state space X to transition space L, following the treatment of Ref. [7].The central objects to be calculated are the trans-transition probabilities p(ℓ|ℓ ′ ) that the next observable transition is ℓ given that the previous was ℓ ′ .An intuitive way to go about this would by brute-force coarse-graining of a stochastic trajectory x, τ = x 0 , τ 0 → x 1 , τ 1 → . . .→ x k , τ k in state space, where x j are the visited states and τ j are the permanence times.Then p(ℓ|ℓ ′ ) can be computed by marginalization of the probability density function at fixed number of jumps k by integrating away the intermediate times and summing over all trajectories between the target state of ℓ ′ and the source state of ℓ but that otherwise do not contain observable transitions, and multiplying by the rate of this latter transition.This direct procedure is illustrated in the Appendix of Ref. [7].A more elegant line is the following.Notice that with probability one the visible activity takes any positive integer value and at any given time t does not depend on future information.Thus the time when the activity reaches a certain value n for the first time is a valid stopping time.Then by the strong Markov property [8] the event of being at x after n visible transitions is also a Markov process in state space.Let p L n (ℓ) be the probability that the n-th visible transition is ℓ.Notice that the probability that the next transition is ℓ given that the previous was ℓ ′ only depends on the target state of ℓ ′ .Thus we conclude that p L n satisfies a discrete-time Markov chain in transition space evolving from some initial probability p L 1 (ℓ) that the first transition is ℓ.The p(ℓ|ℓ ′ ) are the the so-called trans-transition probabilities; we arrange them in a trans-transition matrix P with entries P ℓ,ℓ ′ = p(ℓ|ℓ ′ ).
Both the initial transition probability and the trans-transition probabilities can be obtained from the initial state probability p X 0 and the transition rates r(x|x ′ ) by solving first-transition time problems.In particular the probability that, starting from x, ℓ is the first visible transition and that it occurs in the time interval [t, t + dt) is given by where S is the matrix obtained from R by setting to zero the off-diagonal entries corresponding to the visible transition, namely By integrating Eq. ( 14) from t = 0 to infinity and evaluating at x = t(ℓ ′ ) we find, for all ℓ, ℓ ′ ∈ L, the trans-transition probabilities and, for all ℓ ∈ L, the probability of the first transition It was proven [7,Appendix] that trans-transition probabilities and the initial transition probability are positive and normalized, as they should be: Explicitly, the trans-transition matrix is given by where, letting A \(x1,...,xn|x ′ 1 ,...,x ′ n ) be a matrix from which rows x 1 , . . ., x n and columns x ′ 1 , . . ., x ′ n are removed, we have A proof of these expressions is given in Appendix A.1.

Statement and derivation of the main result
We can now formulate our problem of calculating the probability that the cumulated current ever hits value −1 (case +1 for later) as where f (n) − is the probability that the cumulated current c takes value −1 for the first time at the n−th visible transition.The first is just the probability that the transition occurs right-away: Notice instead that the cumulated current cannot be −1 after two visible transitions: f For the cumulated current to be −1 for the first time at the third visible transition, we need that the first visible transition is → and the second and third are ←, therefore: To go beyond, first notice that all probabilities at even n vanish.At odd n, we need to count all different paths of 2n + 1 steps that perform a ← transition leading from c = 0 to c = −1 for the first time as the last step, and multiply each path by the corresponding probability.Namely, we need to count all different sequences (ℓ 1 , ℓ 2 , . . ., ℓ 2n+1 ) such that: 2) at any intermediate step the number of ← is never greater than the number of → and at 2n the number of ← is exactly equal to the number of →; 3) they have a given amount k ≤ n of ←|→ trans-transitions (which also fixes the number of ←|←, →|→, and →|← trans-transitions).
In fact, this question maps to a well-known enumeration problem: if we replace → with a 45 • unit segment and ← with a −45 • unit segment, we need to count all "mountains" of length 2n that can be drawn without lifting the pencil and that have exactly k peaks (see Fig. 1).This problem is well-known to be solved by the Narayana numbers [9] Fig. 1 A mountain with n = 5 up and down slopes, 3 peaks and 2 valleys.
Therefore for n ≥ 1 the result to our problem is Notice that the prefactor in Eq. ( 26) accounts for the initial transition (which must be ←), for the last transition (which must be ← given that the previous was also ←), and for the fact that in such mountains valleys are one less than peaks.
We now use the fact that Narayana numbers admit the generating function [10] Then, letting we find that Now notice that, using normalization of the trans-transition probabilities Eq. ( 18), we have After some tedious but revealing calculation (see appendix A.2 for details) one obtains that the square root in Eq. ( 27) has real-valued solution in terms of the absolute value | • |.We then find the remarkably simple expression Plugging this latter into Eq.( 30) we find our central result where the two values are obtained respectively for p(→|→) < p(←|←) and for p(→|→) ≥ p(←|←).To express f − as a minimum between two values we used the fact that, because p(←|→)/p(→|←) = [1 − p(→|→)]/[1 − p(←|←)], the second value is monotonically increasing in p(←|←) and decreasing in p(→|→), and is only 1 for p(→|→) = p(←|←).Now notice that the stationary distribution in transition space (eigenvector of the trans-transition matrix relative to eigenvalue 1, P p L ∞ = p L ∞ ) is easily found to be p L ∞ (ℓ) ∝ p(ℓ|ℓ), where ℓ denotes the reverse transition of ℓ (i.e.→ = ←, ← = →).Therefore we can rewrite the above expression as Now consider the probability f + that the cumulated current ever reaches value +1.A quick review of the above derivation promptly leads to where the two values are taken respectively for p(→|→) < p(←|←) and for p(→|→) ≥ p(←|←).

The effective affinity
Let us define This quantity has been given an operational thermodynamic interpretation in Refs.[11,12] as follows.By Eqs. ( 19) and (20) we have where in the second expression is the stationary probability of the system where transition 1 ↔ 2 is removed, i.e.R ∅ p ∅ ∞ = 0 (see Appendix A.2 for a direct proof; the distribution is unique by the assumption that edge 1 ↔ 2 is not a bridge).Notice that this is a stalling system, that is, one where (by non-existence of the transition!) the mean stationary current vanishes.Let us now parametrize rates according to the principle of local detailed balance [13,14] in terms of a energy increment δq xx ′ = −δq x ′ x and a temperature profile x describing the influence of a local bath's degrees of freedom.We assume that temperature T 12 is specific of transition 1 ↔ 2 (that is, its variation does not affect other rates).Then it was proven [11] that there exists a value T ∅ 12 for which the mean current stalls (but here the transition is possible!).Nevertheless, a simple argument shows that the stationary values p ∅ ∞ (2) and p ∅ ∞ (1) are the same as in the system where the transition is removed altogether (see Appendix A.3).We therefore have This latter local expression grants an operational procedure to measure F , on the assumption that δq 12 is measured or theoretically determined by a microphysical theory of the system describing energy levels, that T 12 is tunable, and that the mean current ⟨ ċ⟩ is observable.The procedure consists in tuning T 12 to the value T ∅ 12 for which the observable mean current vanishes.Then, if δq 12 is known, F is determined in terms of the inverse temperature difference.
As regards the global acceptation of affinity mentioned in the introduction, for systems containing a single oriented cycle C it is easily shown [15,16] that F = A is the cycle affinity, namely the ratio of the products of rates along the cycle, in opposite directions For vanishing A (Kolmogorov condition) one finds an equilibrium state with vanishing mean current.From the above relation one immediately finds for the equilibrium temperature the relation For generic multicyclic systems, this latter identification with a specific thermodynamic cycle is not possible.However, the cumulated current c = C c(C) can in fact be envisioned as the sum of the winding numbers over all cycles that include the visible transition (see Refs. [17,18] for some insights on such winding numbers).Notice that a stalling mean current does not imply global equilibrium, as these cycles may have circulation even if overall the visible mean current stalls.An explicit expression of F in terms of such cycles is where w(C) is some cycle weight, independent of the cycle's orientation [11].
Nevertheless, defining entropy production as the Kullback-Leibler distance of random processes from their time-reversed, it has been shown that F ⟨ ċ⟩ is indeed the entropy production estimated by an external observer who only has access to the sequence of visible transitions [15,16].

Special cases and the noria, and a generalization
We consider two special cases where our main results write in terms of the effective affinity.Here we resolve the explicit dependency of the stopping probability in terms of the probability p L 1 of the first transition, f ± = f ± p L 1 .Remember that such probability can eventually be computed from the initial probability in state space p X 0 via Eq.( 17).Finally we generalize the above results to the probabilty of hitting arbitrary low values.

Stationary case
In the first case we sample the initial transition from the stationary distribution.We easily find from Eqs. ( 35) and (36) From an operational point of view this is particularly simple because it only requires to wait long enough for the system to stationarize.Then p L ∞ can be computed explicitly from the time series of the transitions, by just counting the relative frequency of →'s and ←'s.

Cyclic case
In the second case, we prepare the system just after a visible transition is performed and then wait for the same transition to occur again, thus completing a cycle.Therefore for c = +1 we prepare the system at the tipping point of →, which gives p X 0 (x) = δ x,1 so that, after Eq. ( 17) is applied, p L 1 (ℓ) = p(ℓ|→).For c = −1 we prepare the system at the tipping point of ←, which gives p X 0 (x) = δ x,2 and p L 1 (ℓ) = p(ℓ|←).After some calculation trick such as This result is analogous to the one derived in Ref. [1] for unicyclic systems, with the exception that in the unicyclic case the choice of initial state (or, equivalently, the final transition) is not relevant, given that all states share the same cycle and therefore the explicit dependency on the initial state drops and the above result simplifies to where f ± [ • ] is just the probability that the cycle is ever completed in either direction, independently of the initial state.

Hitting −n
The above hitting result for the cumulated current to ever become −1 (for F > 0) lends itself to a simple generalization to the case of the cumulated current hitting value −n, for n ∈ N. Intuitively (given that denumerable + denumerable = denumerable) this is just given by reiterating the hitting problem (renewal property), with the initial condition stabilizing to the previous occurrence of ← just after the first occurrence.One immediately obtains where we rewrote p(→|←)/p(←|→) = exp F ↔ as the effective affinity of a system whose trans-transition matrix P ↔ has the columns swapped with respect to P ; interestingly this auxiliary dynamics also plays a role in formulating the transient fluctuation relation in Ref. [7], but its physical interpretation has still to be clarified.

Fluctuation relations
In the unicyclic case, one easily finds the fluctuation relation In the multicylic case, from Eqs. (49, 48) we have This looks formally like a fluctuation relation, with a caveat: in fluctuation relations the probabilities being compared should be the same, while in this case they are different probabilities, as they are conditioned on two different initial distributions, viz.p(•|→) and p(•|←).This, as we will see, has consequences on the computational or experimental interpretation of data, given that one should prepare different experiments for forward and backward processes and post-select their outcome, which is not desirable.In the next section we comment further on this aspect, arguing that Eq. ( 53) may in fact be the best chance of an estimator of nonequilibrium despite approximations.Furthermore, in Ref. [7] it was proven (Eq.( 21)) that, by sampling the initial transition from distribution p L 1 (ℓ) ∝ p(ℓ|ℓ), the following fluctuation relation holds where we remind that c is given by Eq. ( 10) and p n (c) is the probability that the cumulated current is a certain value c ∈ Z after n visible transitions.One can then further derive the relation where N is any subset of N.This is reminiscent of Eq. ( 53), but notice that these latter are not independent probabilities.Finally, fluctuation relations for single edge currents at stopping times different than the total number of visible transitions (in particular at "clock time" t) do not generally hold -but in the unicyclic case -because the statistics of a specific current depends on all other currents flowing through the network.This is what makes relations such as Eqs.( 53) and (55) particularly appealing, as they are local and phenomenological, and do not depend on knowledge of the whole system.

Estimation of the effective affinity
Many of the above expressions can be used to build estimators of the effective affinity.We will focus on the ones coming from cyclic processes.
Consider M independent realizations of a trajectory performing N visible transitions: Define the cumulated current after the n-th visible transition It has empirical distribution and empirical mean and variance Define the empirical stopping times and the estimators of the stopping probabilities where the minimum is introduced to avoid possible divergences in the case f± = 0 (see also Eq. ( 25) in Ref. [19]).
Notice that due to the finite cutoff on the number of transitions, given Eq. ( 21) these latter are biased.In particular they systematically underestimate (on average) the true stopping probability due to the fact that all occurrences of c = ±1 after N visible transitions are discarded.
Assuming that we can ignore the initial conditions, we can invert Eqs. ( 48) and ( 49) to obtain an estimator of the effective affinity We can compare this to the estimator coming from the stopping fluctuation relation which is generally biased due to the different initial conditions in Eqs.(53).
We complement these stopping-problem estimators with an estimator coming from the theory of linear response out of stalling states [20] and with an estimator obtained from the standard entropy production expression as a Kullback-Leibler divergence (properly regularized to avoid taking log 0) This latter is well-known to be a biased estimator, and better practices in evaluating relative entropies correct these biases but also greatly increase the running time (see Supplementary Material in Ref. [15]).We do not concern ourselves with this issue here.In Fig. 2 we compare the behaviour of these estimators in a simple model.The linear regime estimator Flr performs better near the stalling condition F = 0, while it diverges significantly out of stalling.On the contrary, the cyclic estimator Fcy converges far from stalling, but it systematically suffers from the finite N cutoff.The entropy production estimator Fkl is also biased and noisy due to the tails of the cumulated current's distribution.The stopping fluctuation-relation estimator Ffr instead appears to not be affected by all these issues, despite the approximation due the bias due to the different initial state in Eq. (53).
The reason behind the left-right asymmetry in the above plot is due to the fact that for simplicity we decided to only perturb one rate R 1,2 = exp F and keep all others fixed.This choice is useful to show that entropy production estimators can lead to noise in the tails depending on time-scale separation between rates.Had we distributed the perturbation among R 1,2 and R 2,1 , we would have obtained a more symmetric plot.

A nonequilibrium Boltzmann formula?
Our Eq. ( 1) can be seen as a "nonequilibrium Boltzmann formula" given its similarities with S = log W connecting entropy S and probability 1/W (W being the volume of state space), elaborated by Boltzmann and refined by Planck (Boltzmann's constant set to unity).But with some precautions.
Einstein wrote about the Boltzmann formula: «To be able to calculate W , one needs a complete theory of the system under consideration.If considered from a phenomenological point of view [this] equation appears devoid of content».Einstein then inverted the equation to make it a rule for inferring probabilities from measured entropy differences between equilibrium stateswhich better capture the dynamical nature of processes -and used it to perfection Smoluchowski's theory of critical opalescence [21].However, at least since Kant, philosophers warn us that observations are not independent of conceptions, and therefore deduction from measurements needs theory (the fluctuations of what?), and theory needs the human touch.Still now we don't know which came first, whether the chickens of gases and thermal machines or the eggs of thermodynamics and statistical mechanics [22].At equilibrium the situation is aggravated by the fact that the construction of thermodynamic potentials requires many arbitrary choices by the observer [23], while the pursue of objectiveness requires a description of processes in terms of invariant quantities.
Far from equilibrium, flows of heat to and from the environment are not quantified by differences of a state function, but by "inexact differences".By the so-called principle of local detailed balance ratios of probabilities of forwardto-backward processes have been connected to so-called affinities that quantify the entropy production along cyclic processes, and which are invariant upon the redefinition of the fundamental degrees of freedom [23].However, until recently it has proven difficult to directly connect probabilities and meaningful physical quantities.In fact, despite some claims, there are no predictive variational principles far from equilibrium [24,25].This latter issue is solved by our equations, which allow to connect a statistical property and a physical object in a more direct way than did, for example, fluctuation relations.However, the former concern still plagues our result.What comes first: the egg of F , or the chicken of f − ?Only circumstances can tell.

Relation to a companion publication
The present manuscript is strictly related to a companion work [26] by the same Authors that addresses similar questions.Let us clarify in which ways.
Equation ( 51) is strictly related to Eq. ( 14) in Ref. [26].There the normalized probability p −n of the cumulated current taking minimum value −n is addressed, while in our case f −n allows that, after hitting value −n, the cumulate current may take even more negative values.Therefore we have, intuitively, that this latter is the cumulative distribution of the former Given that p −k is normalized, this identification allows to estimate the escape probability that the cumulated current never actually attains a negative value as p 0 = 1 − f −1 , which in view of Eqs.(34), the explicit expression for the trans-transition probabilibites Eqs.(19) and (20), and the explicit expression for the probability of the first transition Eq. ( 17) allows to express p 0 in terms of the (distribution of) the initial state (see below the explicit expression).
The other main difference between the two works is methodological.Here we follow a constructive but specific approach based on first-transition time techniques and combinatorics, while Ref. [26] is rooted in the more general theory of martingales.In particular in Ref. [26] it is shown that, upon a proper choice of initial state, exp −F c is a martingale, and in particular its expected value ⟨exp −F c⟩ is constant in time.Doob's optional stopping theorem then states that this time can be any proper stopping time.By choosing the moment when the cumulated current hits the boundary values n + > 0 or n − < 0 for the first time, and given that c starts from value 0, one obtains where f (n+) n− is the probability of hitting n − whilst not hitting n + .The nonequilibrium Boltzmann formula follows by taking n − = −1 and n + → ∞, in which limit f Interestingly, similar formulas were derived in an optimization context in Ref. [27].
Finally, here is a short dictionary of equivalent terms and concepts in the two papers: transition rates r(x|x ′ ) here are k u→v there; the observed edge 1 ↔ 2 is y ↔ x; "cumulated" currents c are "integrated" currents J; for stationary probabilities we have ∞ instead of "ss"; the effective affinity F is a * ; the extremum probabability p −n [p L 1 ], given Eq. ( 17), is p J inf x→y (−ℓ|X(0) = x 0 ); the escape probability p 0 is where S is the matrix with entries S u,v = k v→u − δ u,v w∈X ; w̸ =u k u→w if (u, v) ̸ = (x, y), (y, x), else S x,y = S y,x = 0, and we used the explicit expression of the effective affinity Eq. (37), that now translates into When x 0 = x we find where this latter passage follows from the algebraic manipulations in Appendix A.1.We thus recover Eq. ( 16) in the companion paper.We checked computationally the more general equivalence (implied by the theory) of Eq. (69) with Eq. (80) in the companion paper, but a direct proof has remained elusive.

Conclusions
Both martingale and first-transition methods are having a revival in connection to thermodynamic considerations [2-4, 7, 15, 16, 28], and they may lead to independent generalizations and applications of our results.In both approaches, the main open question is the generalization to an arbitrary subset of currents -neither the full entropy production nor a single edge current.
As regards the first-transition approach followed here, as soon as one steps out of the single-edge case the Markov property of the process in transition space is lost.Here the combinatorial approach may allow some exploration.
Since any Radon-Nicodym derivative of two probability distributions over realizations of the process is a martingale, martingales can be used to generalise the results in this paper.From this approach it would seem that one can generate an arbirary number of first-hitting results by building ad hoc auxiliary dynamics.However, the physical interpretation of this class of results may not be clear: it is crucial in our approach that the effective affinity has a clear operational interpretation.In particular, if one could tune its value by just "turning a knob", then the effective affinity is just the difference of that knob's value (in proper physical units) where one wants to perform the experiment and the value of that knob at which the observable current vanishes on average.This local operational interpretation dispenses one to compute the effective affinity from knowledge of all the inner details of the fundamental thermodynamic cycles that influence that particular current.
On a more speculative side, notice that in our derivation we made an arbitrary restriction of the solution of the Narayana generating function, based on the assumption that we expect probabilities to be real-valued.It may be interesting to explore the meaning of the complex-valued solution.
Notoriously, Boltzmann's epitaph is his formula.But it took a whole community (including Einstein, Planck etc.) to digest it.So who's formula is it?
First consider det S \(1|1) and det S \(2|2) .Since the removal of the first line and column, and of the second line and column, both take away the entries (1, 2) and (2, 1) which are the only ones that differ among S and R, these determinants are identical to det R \(1|1) and det R \(2|2) .Therefore from Eq. ( 16) we find where in the second passage we used the well-known fact that for any stochastic rate matrix R the cofactors (−1) x+x ′ det R \(x|x ′ ) are independent of x ′ (this follows for example from Eq. (85) in appendix A.3), and in the third we used the definition in Eqs.(19) and (20).A similar formula is found for p(→|←).
As regards det S \(1|2) (respectively, det S \(2|1) ), in this case the matrix resulting from the removal of the first row and second column only differs from R \(1|2) by entry R 2,1 .Using the Laplace cofactor expansion for determinants we thus obtain and given the definitions in Eq. ( 20) similar formulas as Eq.(73) follow for p(→|→) and p(←|←).
Finally, consider any matrix A, and decompose its determinant in terms that are respectively linear in A 1,2 , linear in A 2,1 or that contain where in the first passage we used Eq. ( 75) and in the second we used the explicit expressions for the paramters.In view of Eqs.(20) this yields the transtransition probabilities in Eq. (19).□ A.2 From the Narayana generating function to trans-transition probabilities where we used the fact that the mean current vanishes by assumption, and on the right-hand side we recognized the stationary equation R ∅ p ∅ ∞ = 0. □.

Fig. 2
Fig. 2 For a fully-connected four-state model with all unit rates except for R 1,2 = exp F and initial state x = 1 (that is p X 0 (x) = δ 1,x ): (dashed) the effective affinity F ; (continuous) estimator Ffr of log f + /f − ; (crossed) estimator Fcy of sign (f + − f − ) log fσ; (bullets) the linear regime estimator Flr ; (triangles) the entropy production estimator Fkl .The ultimate stopping time was set to N = 20 and the number of samples to M = 10000.