Energy Conversion and Entropy Production in Biased Random Walk Processes—From Discrete Modeling to the Continuous Limit

We considered discrete and continuous representations of a thermodynamic process in which a random walker (e.g., a molecular motor on a molecular track) uses periodically pumped energy (work) to pass N sites and move energetically downhill while dissipating heat. Interestingly, we found that, starting from a discrete model, the limit in which the motion becomes continuous in space and time (N→∞) is not unique and depends on what physical observables are assumed to be unchanged in the process. In particular, one may (as usually done) choose to keep the speed and diffusion coefficient fixed during this limiting process, in which case, the entropy production is affected. In addition, we also studied processes in which the entropy production is kept constant as N→∞ at the cost of a modified speed or diffusion coefficient. Furthermore, we also combined this dynamics with work against an opposing force, which made it possible to study the effect of discretization of the process on the thermodynamic efficiency of transferring the power input to the power output. Interestingly, we found that the efficiency was increased in the limit of N→∞. Finally, we investigated the same process when transitions between sites can only happen at finite time intervals and studied the impact of this time discretization on the thermodynamic variables as the continuous limit is approached.


I. INTRODUCTION
Non-equilibrium thermodynamics deals with general laws of a (driven) system transferring energy from one or more heat bath(s) to useful work.The second law, however, restricts this transformation as only part the input energy may be "accessible" as the entropy production, related to heat production, must not decrease [1][2][3].The "system" is usually described by many degrees of freedom, e.g., 10 23 gas particles with individual fluctuating trajectories, where the exchanged heat and the extracted work are determined from statistical averages over all these particles trajectories [4].For "small" systems, such as (bio)polymers, colloid particle, enzymes or molecular motors, the dynamics is described by only few degrees of freedom and the fluctuation of individual "state" trajectories becomes more prominent [3,5,6].Notably, it has been shown, e.g., see Refs.[7][8][9], that the resulting (relative) fluctuations of a statistical averaged observable, like the number of steps R of a molecular motor mimicked by a random walk along a track of molecules, can be related to the entropy production.A relative uncertainty or fluctuation ∆R 2 = 2D/v 2 of the motor steps along the molecular track with diffusion coefficient D and velocity v requires at least a entropy production rate σ of 2k B /∆R 2 .This leads to the inequality, known as the thermodynamic uncertainty relation (TUR) [7,8], Different approaches are used for describing small system dynamics including fluctuations and the resulting * khenning@sas.upenn.edu(stochastic) thermodynamic properties.Two of these approaches have been prominently elaborated in recent years [3]: (i) Dynamics on a discrete set of states: A system is described by its microstates whose dynamics is captured by a master equation.The rates to interchange between microstates is governed by local detailed balance relation determined by a thermal bath [3].
(ii) Dynamics on continuous trajectories: Consideration of individual continuous trajectories of a (driven) colloidal particle whose velocity is described by a Langevin equation where the probability distributions of a particle at position x with velocity v and diffusion D are determined by the Fokker-Planck equation.
Starting from discrete set of states, we can approximate the dynamics by a continuous description, such as the Fokker-Planck equation by standard procedures like the Kramers-Moyal expansion, when (infinite) many states may be visited during the time of observation, e.g., a chemical reaction network of multiple reaction steps [10][11][12][13][14]. Similarly, starting from a continuous description for the dynamics, a discrete representation is obtained by standard mathematical steps of replacing derivatives by finite difference quotients.While such mathematical transformations are expected to lead to equivalent descriptions of the underlying physics, we show below that not all physical observables can be kept invariant under such transformations.In particular, we show that the discrete to continuous limiting process is not unique and depends on which observables are chosen to be invariant under this process.Recent studies have already shown that the entropy production might differ under different coarse-graining schemes since under coarse-graining some "information" is lost while also the mathematical derivation of the differential entropy starting from the discrete Shannon entropy has revealed some discrepancy [15][16][17][18].Here we suggest a different point of view by showing that it is possible to impose an invariant entropy production (or a given heat exchanged with the thermal environment) when proceeding from the discrete to the continuous description of the system dynamics.Under such restriction some other system observable cannot be kept invariant.Alternatively, one may even ask if this observation may be translated into realistic system processes.Put differently, can we, by adding intermediate states between an initial and a target state, optimize the process of transferring input energy into useful work.
The aim of the present manuscript is to investigate the transition from a discrete to a continuous state-space system for an exactly solvable master equation by keeping distinct system observables constant while studying the impact on other observables.Explicitly, we consider a simple model: A Brownian particle moving on a downhill slope with an energy pumping step taking place at constant length intervals L that restores the energy of the particle into its initial value which represents a cycle.In the discrete representation, the particle is a random walker moving among N equally spaced sites per fixed length L. N → ∞ that represents the continuous limit for the cyclic dynamics.More details about the model are provided in Sec.(II).Other processes with period boundary conditions can be mapped into this form.Realization of such a process can be molecular motors which transform free energy liberated in a chemical reaction by a succession of steps on a track into mechanical work (motion) [5,8].Driven rotational Brownian particles through periodic potential wells for experimentally testing thermodynamics laws represent another example of many others [19].Given the specific chosen conditions, like constant speed v and the entropy production rate of the cycle, we can study the diffusion coefficient D in dependence of the number of sites N per cycle and the continuous limit when N → ∞.Interestingly, under an additional action of an opposing force during the process under study, we can also determine the impact of the number of sites per cycle on thermodynamic performance by transferring input energy into useful work.Furthermore, when keeping the energy drop per cycle and the speed constant, we interestingly find that in the limiting case N → ∞ the diffusion coefficient (equivalent to the variance in the site distribution) approaches the thermal Einstein relation when N → ∞ usually expected only for small velocities in the linear response limit of small driving.
The paper is structured as follows.In the first Section II, we describe the cyclic process as biased random between (energy) sites governed by a master equation restricted by periodic boundary conditions.We further discuss the implication for the entropy and heat production.In Sec.III, we study the entropy production rate and heat transfer rate into the environment of the cyclic process given the velocity v and diffusion constant D for different number of sites N in the dynamical (relaxation from initial site) and in steady state.We then deter-mine the performance in transferring input power to useful output power for the cycle process under an opposing force by varying N and, so, by going from the discrete to the continuous (N → ∞) state space.We study in Sec.IV the transition to N → ∞ given constant entropy production for the process, and, under either constant velocity v or diffusion coefficient D. The impact on the respective system variable when increasing N is discussed.Section V is devoted to examining the randomness parameter of the cycle, which is dependent on the number of sites N , while maintaining various physical observables as constants.Through this randomness factor, we can also determine the variance in cycle completion time and study its dependence on N .In Sec.VI, we investigate the process by the same biased random walk, but where the system evolves at finite time intervals, and study the impact of this time discretization on the entropy production and diffusion coefficient.Sec.VII concludes this work.

II. CYCLIC PROCESS
We consider a cyclic process, where the n = 1, . . ., N sites are aligned on a circle such that each site has two neighboring sites with the periodic boundary condition N +n = n with N ≥ 3. The process dynamics is captured by a biased random walk with forward and backward transfer rates α and β, which are equal at each site, such that the master equation for the probability distribution of the n-site of the cycle reads Ṗ (n, t) = αP (n − 1, t) + βP (n + 1, t) − (α + β)P (n, t). ( Starting form a well defined site at t = 0 one finds in the limit of long times t → ∞ the (nonequilibrium) steady state distribution for the n site to be P (n, t → ∞) = 1/N (see Appendix C).When writing the master equation (2) in the form Ṗ = M P we see that the matrix M is irreducible with the dominant eigenvalue λ = 0 while all other eigenvalues have a strictly negative real part according to the Perron-Frobenius theorem signaling the existence of a stable steady state vector [20,21].While our cyclic process can take place in the state space of any given system, it is convenient to consider the equivalent random walk with forward and backward rates α and β on a cycle of constant circumference L = 2πR.In this equivalent random walk problem the velocity v and diffusion constant D at steady state are defined by [14] v and where ∆x = 2πR/N is the equidistant step size with R and N being the radius of the cycle and the number of total sites N of the cycle respectively.Evidently, the time for completing a full cycle is τ = N (α − β) −1 .The transition between the neighboring sites, n → n ± 1, are considered as autonomous Markov jump processes where each site has its energy E(n).Thermodynamic consistency is introduced by the local detailed balance condition where, for simplicity, we assume isothermal conditions with T on all sites, and where ∆E = E(n)− E(n+ 1) > 0 is taken the same for all nearest neighbor sites.Note that we count the heat exchange with the bath as Q ≡ −∆E, i.e., if Q < 0 the amount of heat is taken from the system to the bath while the environment provides heat to the system, if Q > 0. This implies that during a cycle, the amount of heat ∆Q ≡ −N ∆E = k B T N ln α β is dissipated, so E(N + 1) ≡ E(1) − N ∆E may be compared to a downhill process of energy loss N ∆E.To remain consistent with the periodic boundary conditions we further assume that between sites N and N + 1 = 1 an upward energy jump occurs, in which the same amount of work, W = N ∆E, is returned to the system.
Next we calculate the time-dependent change of the Boltzmann-Gibbs entropy, [2,22] for the system using the master equation (2).We obtain Note that in Eq. ( 6) we could replace ln P (n, t) by ln[P (n, t)/P (m, t)] because the sum that multiply ln P (m, t) vanishes.The result (8) can be recast in the form Ṡ(t) = Ṡe + σ(t), [23,24] where, using n P (n, t) = 1, Next we show that first term Ṡe (Eq.( 9)) is the entropy flow into the environment while the second term is entropy production rate σ(t) (Eq.( 10)).
Consider first Eq.( 9).Because the rates α and β were assumed not depend on the site identity, Ṡe is timeindependent.Ṡe can be written as where J ≡ α − β is the cumulative flux: The sum over nearest neighbor site-pair fluxes (see first expression in Eq. ( 9)).Q is the heat exchanged during a single nearest neighbor transfer event with the environment according to the local detailed balance relation in Eq. ( 5).The product JQ in Eq. ( 11) is the heat flux into the thermal environment of temperature T per cycle.For the given model this heat exchange is time-independent.
According to the Clausius principle the change in system entropy is bounded by the (negative) heat amount exchanged with the environment Ṡ ≥ JQ T where equality is reached for reversible processes [1].Motivated by this inequality, one defines the total entropy production rate by σ(t) = Ṡ − JQ T ≥ 0. Indeed, the second term σ(t) (Eq.( 10)) meets the two important properties: (i) It is nonegative because the firs term in Eq. ( 8) can be recast into (x−y) ln(x/y) ≥ 0 and (ii) it vanishes for thermal equilibrium, when microscopic reversibility or detailed balance condition, αP (m, t) = βP (n, t), is obeyed and no entropy is produced.At (nonequilibrium) steady state Eq. ( 6) is zero and the entropy production equals the negative entropy flow into the environment σ = − Ṡe , Ref. [24,25].
In the following, we investigate the entropy production rate σ(t) and the physical measurable entropy or heat flow Ṡe into the environment for prototype N -site cyclic process given different conditions.In particular we are interested in the limit N → ∞ to investigate the transition from the discrete to the continuous state-space.

III. N-SITE CYCLIC PROCESS UNDER CONSTANT VELOCITY v AND DIFFUSION CONSTANT D
We first consider the dependence of the number of sites N per cycle on the system that is done under the conditions that (a) the speed v (Eq.( 3)) and (b) the diffusion constant D (Eq. ( 4)) are kept constant.The first condition requires that the time, τ = N (α − β) −1 , for a full completion of the cycle remains constant.Together with the second condition, the forward and backward rate must depend on N .We get To ensure the positivity of the rate (13), β ≥ 0, it immediately follows that v ≤ N D/Rπ.Note that under the conditions restricted by Eqs. ( 12) and ( 13), the detailed balance relation ∆E = k B T ln α β (Eq.( 5)) will be a nonlinear function of N .Consider the entropy production rate σ(t) (Eq.( 10)) given a well define initial cycle site n = 1 at t = 0, so P (n = 1, t = 0) = 1.The entropy production rate (the entropy produced per unit time) decreases over time before reaching a (N -dependent) steady state value, see Fig. 1.This observation can be understood as follows.Initially, the system starts at a given site and evolves in its cyclic dynamics according to the master equation (2).At later time t the system will be found at a given site n with probability P (n, t) (see Appendices A-C).The loss of the initially "knowledge" about the exact system site n increases the entropy of the system.As the probability distribution relaxes to its steady state distribution P n,SS ≡ P (n, t → ∞) = 1/N (see Appendix C) the entropy production rate, σ, decreases to its steady state value Note that entropy is produced at constant rate when running the cyclic process under steady state conditions.Interestingly, as seen in Fig. 1, the more sites N are included in the cycle of finite length the faster the entropy production rate decreases to its steady state value.At steady state the entropy production equals the negative (in this model time-independent) entropy flow into the environment σSS = − Ṡe .Consider the N -dependence of the the steady state entropy production rate σSS for an N -site cycle, Eq. ( 14).The dependence is shown in Fig. 2. For N → ∞ the steady state entropy production rate ( 14) is reduced to k B v 2 /D.This minimum entropy production rate in this limit can be understood as the forward and backward rates (Eqs.( 12) and ( 13)) will become more alike.Therefore, consecutive transitions between sites n and n + 1 become more time-symmetric and less entropy per step, ∆E/T = ln [α(N )/β(N )]/T , is produced.
Interestingly, for N → ∞, the TUR relation, the minimal required entropy production rate σ for a given relative fluctuation 2D/v 2 (as introduced in Eq. ( 1) for a molecular motor moving along a molecular track) states equality, σ2D/v 2 = 2k B .In this limit, the dynamics of the cyclic process is comparable to a continuous Brownian diffusive motion of a particle of constant speed v and diffusion D described by a Fokker-Planck equation [26].
Assume that a forward (downhill) step on the cycle takes place against a constant applied force f .Then, the local detail balance relation (5) must be redefined by ∆ Ẽ ≡ E(n) − E(n + 1) − f ∆x = ∆E − f ∆x > 0, as part of the energy per step ∆E is transferred to work f ∆x, where ∆x is the step size.We defined the system's heat exchange per step with the bath as Q ≡ −∆ Ẽ, while Q < 0 is the amount of heat taken from the system to the bath.At steady state the heat flow into the bath is − Q = T σSS , see discussion in Sec.II.The power output per cycle when running against the force at steady state is In Eq. ( 15) αP (n) − βP (n + 1) is the probability flux between neighboring sites and where n P (n) = 1.By the first law of thermodynamics the total supplied power must be P − Q such that we can define the thermodynamic efficiency as [8] η(N ) where σSS (N ), Eq. ( 14), depends on the total number of sites N per cycle.The cycle can be compared to a process of going down a slope against a constant force and with friction.The friction force is usually taken as F f r = γv such that the related heat dissipated in the cycle per unit time is for the present cycle process using Eq. ( 14).Note that, as expected, the friction coefficient γ goes to its linear response value for a small velocity, γ(v → 0) → k B T /D, independent of N .Surprisingly, as depicted in the inset of Fig. 3, γ takes the same value also for finite v in the limit of N → ∞.In this limit the system assumes some features similar to equilibrium even though the flux is finite.Not only in the limit v → 0, but also in the limiting process N → ∞ can be compared to thermodynamic cycles where the system changes adiabatically slowly to be always in thermal equilibrium throughout the process.Equivalently, in the limit N → ∞ we can write the thermodynamic efficiency (16) as [27] Fig. 3 portrays the efficiencies for different N .As expected, the efficiency increases with N as less heat − Q = T σSS will be produced per cycle given constant power output P .Interestingly, the slower the velocity v is chosen the more one can reach maximal efficiency.The efficiency is bound from above, η ≤ 1, while equality is reached for v = 0 or f → ∞, D → ∞, see Eq. ( 17).The last conditions, however, do not produce useful output power, as the cycle will stop.

IV. CYCLIC PROCESS WITH CONSTANT VELOCITY v OR DIFFUSION CONSTANT D AND CONSTANT ENERGY DROP PER CYCLE
As stated in the introduction, the process of going to the continuous description for the state-space dynamics (N → ∞) is not unique.We therefore use the same methodology as above to describe the unicyclic process as 1-D random walk but where we require now the total entropy produced per cycle is constant (given by the constant energy drop per cycle) and either (A) the velocity v or (B) the diffusion coefficient D are constant.We assume that all energy invested into the system is dissipated as heat to the environment, so W = σT .The steady state entropy production σ is where N is the total number of sites.At steady state the entropy production equals the heat going into the environment, see Sec.II.(A) Given the constant velocity v, Eq. ( 3), the forward and backward rates are where R is the radius of the cycle (see Sec. II).(B) Given the constant diffusion coefficient D, Eq. ( 4), the forward and backward rates are Consider first case (A).Eq. ( 4) with Eqs. ( 19) and ( 20) lead to which is shown in Fig. 4. The diffusion coefficient decreases with increasing N .As the diffusion coefficient is the variance of the site distribution on our equivalent cycle, see Eq. ( 4), the related fluctuations in localization of a site is reduced during a cycle with increase of N .Assuming that the energy falls uniformly along the cycle, so that ∆E = T σ/N , we find that in the regime of linear response v → 0 and in the limit N → ∞ the diffusion coefficient is captured by the (Einstein) relation D = v/F k B T = µk B T with the mobility µ = v/F and the related force F = T σ/2πR in analogue to the friction coefficient (as discussed in Sec.III) [28,29].Next, consider case (B).Eq. (3) with Eqs. ( 21) and ( 22) leads to the velocity on N The velocity v increases with N and so the time of a full completion of the system cycle τ = 2πR/v is reduced and minimizes for N → ∞, see Fig. 5.

V. RANDOMNESS PARAMETER AND VARIANCE IN CYCLE COMPLETION TIME
In the previous Sections III and IV, we have analyzed, for a biased random process, the dependence of different physical observables on the number intermediate sites N taken to complete a given cycle of operation, as a way to demonstrate the non-uniqueness of going to the continuous limit (N → ∞) of this process.Another perspective of this problem has been studied in Refs.[30][31][32][33], where the values of physical observables associated with enzyme catalyzed cycles were used to set bounds on the number of intermediate cycle steps [34].It was pointed out that in addition to the average speed v and diffusion coefficient D, their ratio provides important information on the observed walk statistics [30][31][32][33].Explicitly, we characterize the random process in terms of their forward and backward rates, α and β, respectively, and an equidistant step length ∆x.If the process starts at n = 1 and a random site n is reached at time t, then the process will be on average at site n = vt/∆x, whilst the random diffusive process produces a variance in the site by δn 2 ≡ n 2 − n 2 = 2Dt/∆x 2 , see Eqs. ( 3) and ( 4) in the limit t → ∞.These two quantities can be combined into a randomness parameter, which for the given step size ∆x reads [30][31][32][33] r ≡ where v and D are defined by Eqs. ( 3) and ( 4), respectively.Alternatively, we may consider the random passage time τ at which, starting from n = 1, the walk reached for the first time the site N , namely a distance N ∆x from the starting point.For walks of uniform step length and finite bias, it has been shown [30,31] that for large enough N , the randomness parameter can be expressed in terms of the first two moments of the passage time distribution where τ is the average time for a cycle completion and δτ 2 its variance.Note that for many enzyme reaction cycle, the backward reaction rates are often sufficiently low as to be negligible.In such cases the pathway consists of a sequence of N forward reactions only and the randomness parameter (r = N −1 min ) can be used to estimate the minimal number of kinetic sites that compose the underlying kinetic model [31,33].In general, when considering forward and backward steps and using the average cycle completion time τ , we can calculate the variance in cycle completion time to We can now apply the results of Secs.III and IV to examine the behavior of these observables in our different limiting cases.In our case N corresponds to the number of sites per cycle and, consequently, τ is the time for the process to complete the cycle.Increasing N is done with eventually approaching a continuous description, so the cycle length N ∆x = 2πR is kept fixed.
Consider first the condition of a constant velocity and diffusion constant by increasing the number of sites N , see Sec.III.Using Eqs. ( 12) and ( 13) in ( 25), the randomness parameter is and is linear in N .
Next, for the condition of constant entropy production σ per cycle under either constant velocity v or constant diffusion D, see Sec.IV, we find in both cases, using the respective rates in Eqs. ( 19) and (20) or Eqs. ( 21) and ( 22), the randomness parameter to be The randomness parameter ( 29) is determined by the number of sites N and the thermodynamic entropy production or heat dissipation into the environment which equals, when neglecting the movement against an external force, the energy drop per cycle, see Sec.II.Note that r tends to infinity in the limit σ → 0 and in the continuous limit N → ∞ since both limits reflect the equilibrium situation where forward and backward rates will be alike.
In the limit σ → ∞, given finite N , we find r = 1, which is expected for the so-called "Poisson" motion since the infinite energy drop per cycle leads to an unidirectional motion [30].
With the randomness parameter at hand we can now study the variance in the cycle completion time (Eq.( 27)).It has been shown that for a biased random walk for large N that the average completion time τ = N (α − β) −1 [35].Consider first the cases (A) of keeping the velocity and diffusion coefficient constant as well as (B) keeping the entropy production per cycle and velocity constant.We find for (A) given Eq. ( 28), together with Eqs. ( 12) and ( 13) in Eq. ( 27) where τ = N (α − β) −1 the variance in cycle completion time to be  24).We show δτ 2 for three different choices of σ.
and, equivalent for case (B) by using Eq. ( 29), together with Eqs. ( 19) and ( 20) in Eq. ( 27) In both Eqs. ( 30) and ( 31) the variance in cycle completion time increases monotonously with site number N .This reflects the fact that with increasing number of sites N per cycle the inter-site rates become more alike which increases the overall "randomness", and, thus, δτ 2 for the total cycle completion.
In contrast when N is changed while keeping a constant diffusion coefficient and entropy production, Eq. ( 26), together with ( 24) and (29) lead to Interestingly, δτ 2 goes through a minimum with increasing site number N , see Fig. 6.It should be kept in mind however that Eq. ( 26) and, consequently, Eq. ( 27), was derived under the assumption N is large so that this observation should not be regarded conclusive.

VI. N-SITE PROCESS WITH TIME STEP DISCRETIZATION
Consider now the same biased random walk process in which a cycle of length 2πR is traversed in N steps, so that N ∆x = 2πR, but where the system is restricted to move (by intersite distance ∆x) only at finite time intervals ∆t.Indeed, small systems which are periodically driven can be thought of as discrete-time processes, see [36].As shown below, ∆x and ∆t, are not independent of each other but some freedom exists in their choices.N ∆E is the energy drop per such cycle (see Sec. II, recall that ∆E determines the detailed balance ratio of the forward and backward rates according to Eq. ( 5)).The probability to be at site n, namely at position x = n∆x on the cycle at time t = M ∆t, is governed by the Makrov chain [37,38] + (1 − α∆t − β∆t)P (x, t).
Here, α∆t and β∆t are the probabilities (both assumed linear in ∆t) to move a step forward and backward, respectively.Note that (α + β)∆t ≤ 1 has to be imposed in Eq. ( 33) to ensure positivity.As before, the process has periodic bounderies so that, after the final site n = N has reached, it restarts at the beginning n = 1 and its original energy is restored by some external work reservoir between sites N and N + 1 = 1 (see discussion in Sec.II).We use this model to study the effect of time discretization on dynamical properties of the process.To calculate the velocity v and diffusion coefficient D, we determine the generating function P (s, t) = x s x P (x, t) where the moments can be calculated by x m = (s∂/∂s) m P (s, t)| s=1 .For the initial condition P (x = 0, t = 0) = 1 (so that P (s, t = 0) = 1) we find the generating function to be = [α∆ts ∆x + β∆ts −∆x + (1 − α∆t − β∆t)] M .
The velocity and diffusion coefficient are determined as follows.Starting at x = 0, we find in the long time limit (at steady state) t → ∞ (see details in Appendix D) The velocity v, Eq. ( 35), is the same as in the continuous time case (Eq.( 3)), whereas the diffusion coefficient in Eq. ( 36) is smaller by v 2 ∆t in comparison to the continuous time case of Eq. (4).Refs.[39,40] associate the bigger variance in the continuous-time master equation with higher fluctuations in the total number of hops observed in a given time interval.Next, consider the process as N increases.As in Secs.III and IV, we may consider an increase of N while keeping v and D constant or while keeping only one of them together with N ∆E constant.As examples of the effect of moving in discrete time steps, we study the cases (A) constant v and D and (B) constant v and N ∆E.
(A) Keeping v and D constant, we scale the rates again analog to Eqs. (19) and (20) with the total site number N (given by the intersite distance ∆x = 2πR/N ) Note that the modification of the rates α and β depends on ∆t.As before, this rescaling implies also a change in ∆E (see discussion in Sec.III) so that the detailed balance relation is maintained.The condition (α+β)∆t ≤ 1 in Eq. ( 33) together with Eqs. ( 37) and ( 38) restricts the choices for ∆t given ∆x = 2πR/N to which implies that ∆t and ∆x cannot be assigned independently of each other.In the limit N → ∞ (∆x → 0), this inequality (Eq.( 39)) becomes 0 ≤ ∆t ≤ ∆x 2 /2D.Consider next the entropy production for this discrete hopping process.The average entropy change per step is where the two terms represent the entropy change in a forward and backward step multiplied by the probabilities that the respective step occurs.The rate of entropy change at steady state (entropy production rate) is given by Interestingly, comparing the resulting expression (41) to its analog (14) for the continuous master equation, we obtain a similar result, but with an additional term v 2 ∆t/2 added to the diffusion constant.The additional term effectively modifies the TUR relation (Eq.( 1)) as the relative uncertainty 2D/v 2 changes.A similar observation was made in Ref. [41].The resulting entropy production rate in Eq. ( 41) for given site number N per cycle is reduced if we choose a finite ∆t (given the restriction on choices of ∆t by Eq. ( 39)), see Fig. 7.This might be understood since during a given time interval the variance in position x on the cycle is reduced (Eq.( 36)) by allowing intersite hops only in intervals ∆t.In the continuous limit (N → ∞, ∆x → 0 and ∆t → 0), however, Eq. ( 41) Here, v (chosen to be v = 3D 2πR ) needs to be within the bound implied by β ≥ 0 of Eq. ( 38) for N ≥ 3. (i) The blue circled line represents the entropy production rate in the continuous time limit ∆t → 0. (ii) The black squared line is the entropy production rate for a discrete time process using the maximal ∆t allowed by Eq. ( 39), ∆t = which is the same result obtained in this limit in Sec.III (see Eq. ( 14)).(B) Consider next the dependence on N under the condition of a constant energy drop N ∆E per cycle and constant velocity v.This is the analogue consideration as in Sec.IV, but where the intersite hops are only allowed at time intervals ∆t (that are restricted by the given N according to Eq. (46) below).We assume that all energy invested into the system is dissipated as heat into the environment, so the steady state entropy production per cycle is where N is the total number per cycle.Given the constant velocity, Eq. ( 35), the forward and backward rates are equivalent to Eqs. ( 19) and ( 20) The restriction (α + β)∆t ≤ 1 in Eq. ( 33) together with Eqs. ( 44) and (45), limits the choices for ∆t given ∆x = 2πR/N to In the limit N → ∞ (∆x → 0 and ∆t → 0), the inequality in Eq. ( 46) becomes 0 ≤ ∆t ≤ ∆xσ/(vN k B ).For a given N and ∆x = 2πR/N , the time step ∆t needs to satisfy the inequality (46).Here we take ∆t with 0 ≤ a ≤ 1 and use Eq. ( 36) to get Interestingly, the dependence on time discretization translates here to a dependence of D on the choice of a, see Fig. 8.Given N for different finite time ∆t (scaled between 0 ≤ a ≤ 1) the diffusion coefficient, and consequently, the variance in x = n∆x, are strongly reduced for increasing ∆t (less fluctuations in the total number of transitions for finite ∆t).As expected, in the limit N → ∞ and ∆t → 0, the diffusion constant takes the form D = v/F k B T = µk B T with the mobility µ = v/F and the corresponding force F = T σ/2πR, [28,29], see discussion in Sec.IV.Therefore, given the chosen ∆t and its above discussed effect on D, the latter increases or decreases with N to the final value D = µk B T as depicted in Fig. 8. (where 0 ≤ a ≤ 1) and depict D(a, N ) with a = 0 (blue curve), a = 0.6 (black curve) and a = 1 (green curve).We choose the produced entropy per cycle to be σ/kB = 10 for N ≥ 3.
To summarize this section, when describing the dynamics of a (cyclic) process in discrete time intervals, thermo-dynamical properties of the process, e.g., the entropy production or diffusion coefficient, are affected by this time discretization.The discretization in time, however, cannot be chosen arbitrarily but must obey the bounds given by the state-space discretization of the process.If the process dynamics is going to a continuous state-space dynamics (N → ∞), the time evolution needs to be described by intervals ∆t → 0, i.e., equivalently by a continuous time scale and all effects vanishes.

VII. CONCLUSION
In this paper we have investigated a cyclic and unithermal thermodynamic process by a model of a biased random walk between N sites on a cycle of a given length by an exact solvable master equation.We note that many dynamical site (state)-space processes with periodic boundary conditions may be equivalently mapped to this process.The limit N → ∞ corresponds to the continuous limit that is usually captured by a Fokker-Planck equation.This limit is taken by keeping low order moments of system observables, i.e., the velocity v and diffusion coefficient D constant.We have shown that the entropy produced, or, equivalently, the energy drop per cycle is reduced when going to this continuous description.This has direct consequences for the efficiency of transferring input power into useful output power when an opposing force acts on the cyclic process.In particular, more power can be extracted from the process with increasing number of sites N per cycle length.
An important outcome of our analysis is that the procedure in going to the continuous description of the process is not unique and depends on the physical observables that are assumed to be invariant under this limiting process.In addition to taking the limit N → ∞ while keeping v and D constant we have also analyzed this limiting process while keeping v or D and N ∆E constant.Interestingly, considering the limiting process under constant v and N ∆E, we have shown that diffusion coefficient D for finite cycle velocity v in the limit N → ∞ takes the same value as in linear response v → 0 limit.Additionally, when analyzing the cycle randomness statistics and, in particular, the variance in the cycle completion time δτ 2 .We find that with increasing N , δτ 2 increases signaling the increasing randomness in the cycle.Finally, we have studied the dependence on N in the case where the transitions between sites are only allowed at fixed time intervals ∆t.We find that not only the entropy production rate per cycle (when v and D are kept constant), but also the diffusion coefficient D (when v and N ∆E are kept constant) are strongly affected by the way time discretization is introduced for a given N .
In conclusion, on can use the total site number N as a control parameter to design "usefull" physical and thermodynamic (cycle) processes by keeping desirable observables constant and affecting others.It may provide valuable insights into the engineering of small (molecular) machines capable of performing specific tasks with high efficiency and precision.Further investigations of the discrete to continuous transition in state-space and its potential impact on information-to-work conversion are kept future research.

FIG. 1 .
FIG.1.Entropy production rate σ(t) against time plotted for different N -site cycle by keeping the velocity v and diffusion coefficient D constant.The velocity v (here chosen to be v = 3D 2πR ) needs to be within the bounds implied by β ≥ 0 of Eq.(13) for N ≥ 3. The steady state value v 2 /D (black curve) is reached for N → ∞.

FIG. 2 .
FIG. 2.Entropy production rate at steady state σ(t → ∞) ≡ σSS[= − Ṡe] against the number of sites N per cycle by keeping the velocity v and diffusion coefficient D constant.The velocity v (here chosen to be v = 3D 2πR ) needs to be within the bounds implied by β ≥ 0 of Eq. (13) for N ≥ 3. The red line is the entropy production rate in the limit N → ∞ and takes the value σSS(N → ∞) = v 2 /D.

FIG. 3 .
FIG.3.Efficiency at steady state η plotted against the number of sites N per cycle by keeping the velocity v and diffusion constant D constant.We depict η for v1 = 0.5v, v2 = 0.8v and v3 = v, where we have chosen v = 3D πR to be within the bounds implied by β ≥ 0 of Eq. (13) for N ≥ 3. We have chosen the force f ≡ vT σSS(N → ∞) = 3kBT /πR.Inset:Damping constant γ ≡ γ(v)D/(kBT ) = T D σSS/(v2kBT ) in dependence of the velocity for different N .

FIG. 4 .
FIG.4.Diffusion constant D plotted against the number N of sites per cycle by keeping the entropy σ per cycle and velocity v constant for N ≥ 3. We choose the produced entropy per cycle to be σ/kB = 10.The red line is the value D(N → ∞) → v2πRkB/σ in the limit of N → ∞.

FIG. 5 .
FIG.5.Velocity v plotted against the number N of sites per cycle by keeping the entropy σ per cycle and the diffusion coefficient D constant for N ≥ 3. We choose the produced entropy per cycle to be σ/kB = 5.The red line is the value v(N → ∞) → Dσ/πRkB in the limit of N → ∞.

FIG. 6 .
FIG.6.Variance in cycle completion time δτ 2 plotted against the number N of sites per cycle by keeping the entropy σ per cycle and the diffusion coefficient D constant for N ≥ 3 while the velocity v(N ) results from Eq. (24).We show δτ 2 for three different choices of σ.
FIG. 7.Entropy production rate at steady state, σSS, against the number per cycle N by keeping the velocity v and diffusion coefficient D constant.