Next Article in Journal
Information Hiding Based on Two-Level Mechanism and Look-Up Table Approach
Next Article in Special Issue
About Stability of Nonlinear Stochastic Differential Equations with State-Dependent Delay
Previous Article in Journal
Some Fuzzy Riemann–Liouville Fractional Integral Inequalities for Preinvex Fuzzy Interval-Valued Functions
Previous Article in Special Issue
Nonlinear Fokker-Planck Equation for an Overdamped System with Drag Depending on Direction
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Exact Time-Dependent Solutions and Information Geometry of a Rocking Ratchet

1
Fluid and Complex System Research Centre, Coventry University, Coventry CV1 2TT, UK
2
Department of Applied Mathematics, University of Leeds, Leeds LS2 9JT, UK
*
Author to whom correspondence should be addressed.
Symmetry 2022, 14(2), 314; https://doi.org/10.3390/sym14020314
Submission received: 30 December 2021 / Revised: 15 January 2022 / Accepted: 1 February 2022 / Published: 3 February 2022

Abstract

:
The noise-induced transport due to spatial symmetry-breaking is a key mechanism for the generation of a uni-directional motion by a Brownian motor. By utilising an asymmetric sawtooth periodic potential and three different types of periodic forcing G ( t ) (sinusoidal, square and sawtooth waves) with period T and amplitude A, we investigate the performance (energetics, mean current, Stokes efficiency) of a rocking ratchet in light of thermodynamic quantities (entropy production) and the path-dependent information geometric measures. For each G ( t ) , we calculate exact time-dependent probability density functions under different conditions by varying T, A and the strength of the stochastic noise D in an unprecedentedly wide range. Overall similar behaviours are found for different cases of G ( t ) . In particular, in all cases, the current, Stokes efficiency and the information rate normalised by A and D exhibit one or multiple local maxima and minima as A increases. However, the dependence of the current and Stokes efficiency on A can be quite different, while the behaviour of the information rate normalised by A and D tends to resemble that of the Stokes efficiency. In comparison, the irreversibility measured by a normalised entropy production is independent of A. The results indicate the utility of the information geometry as a proxy of a motor efficiency.

1. Introduction

Symmetry plays a key role in physics and other sciences. One important example is the so-called Brownian motor (e.g., [1,2,3,4,5,6]) whose function hinges on the very presence of symmetry-breaking. As its name indicates, it produces uni-directed motion and works even in the absence of any net macroscopic forces and potential gradients via a noise-induced transport [2]. Specifically, unlike man-made deterministic motors where noise has a negative effect on its performance, the Brownian motor works in a noisy environment far from equilibrium in the presence of spatial asymmetry; thermal fluctuations are preferentially rectified in one direction due to the asymmetry to allow them in the favoured direction while blocking those in the opposite direction [4].
It is a useful mathematical model of molecular motors [5] of the size O (1–100) nanometres in living organisms that play a vital role in organising and orchestrating various transport processes and movement in cells. Due to their small size, they produce kinetic energy which is comparable to thermal fluctuating energy and consequently have small inertia. Thus, their motion is approximated by an overdamped stochastic process. Important examples include myosin [7], responsible for muscle contraction in cardiac and skeleton muscles, or kinesin/dynein for pulling cargos (e.g., organelles). They have the capability of producing force directly rather than via an intermediate energy, by converting chemical energy, e.g., adenosine triphosphate (ATP), to kinetic energy.
One of the key questions is if there is an optimal level of fluctuations that maximises the motor performance. Previous works have suggested that the answer to this question is likely to depend on how a motor performance is defined, such as by the current, work or different types of efficiencies [1], and if an external force introduces an additional symmetry-breaking in time [8,9,10].
For instance, in the rocked thermal ratchet model where a spatially periodic sawtooth (asymmetric) potential V 0 ( x ) is rocked periodically in time by a symmetric square wave force G ( t ) , the optimal current was obtained analytically for a finite fluctuation level (temperature) D = D c under the adiabatic approximation [8] in the limit of a very slow time-variation. Using a similar adiabatic approximation, ref. [9] showed that the peak of efficiency is different from the peak of current and that the efficiency degrades as the fluctuation increases. However, the numerical simulations without using the adiabatic approximation [10] showed different results, the efficiency being optimised for a finite D = D c . On the other hand, even under the adiabatic approximation, the efficiency can be a peaked function of D if temporal asymmetry is introduced in the external force G ( t ) [11]. These are some examples demonstrating the importance of exploring a broad range of parameter values and calculating a fully time-dependent solution without making approximations such as slow or fast variation [12].
Given strong fluctuations, a Brownian motor constitutes an important example of non-equilibrium complex systems where traditional equilibrium thermodynamics or statistical physics do not hold. Therefore, a Brownian motor provides a useful framework in which to develop or consolidate/test newly emerging theories of non-equilibrium complex systems. In particular, far from equilibrium, statistical properties change with time and the time-evolution of a system does not obey time-symmetry (time-irreversibility) even when the external force is symmetric in time. As a measure of time-irreversibility, entropy production, fluctuation theorems (e.g., see [13,14,15] and references therein), etc., are investigated. Furthermore, since the proposal of Feymann’s ratchet, there has been growing interest in information theory [13,16,17,18] to understand or else optimise a Brownian motor [19,20,21,22,23].
In this paper, we propose to investigate the dynamics of a Brownian motor from the point of geometry and the distance (metric tensor) by utilising the information geometric theory—the application of the differential geometry to probability and statistics [18,24,25,26]. It is a powerful tool for elucidating the disparity between different probabilities as well as for linking complexity and geometry (e.g., see [18,24] and references therein). To capture a temporal variation, our focus will be on the path-dependent information geometric concept (information length and rate) that quantifies the time evolution of a system in terms of a dimensionless distance in a statistical space [18,27,28,29,30,31,32,33,34,35,36,37,38] or the change in information. Its key properties are summarised in Section 2.
The main aim of this paper is to investigate exact time-dependent solutions of a Brownian motor under different conditions by varying the values of parameters in an unprecedentedly wide range. By utilising an asymmetric sawtooth periodic potential and three different types of periodic forcing G ( t ) (sinusoidal, square and sawtooth waves) with period T and amplitude A, we investigate the performance (mean current, Stokes efficiency) of a rocking ratchet in light of thermodynamic quantities (energetics, entropy production, entropy flow) and the path-dependent information geometric measures. For each G ( t ) , we calculate exact time-dependent probability density functions (PDFs) by numerically solving the Fokker–Planck equation [39] under different conditions by varying T, A and the strength of the stochastic noise D. It is worth noting that this paper focuses on elucidating the path-dependent information geometry in the efficiency of the Brownian motor and comparing with some of the popular measures of irreversibility (e.g., entropy production).
The remainder of this paper is organised as follows. Section 2 provides our model and key thermodynamic measures, motor efficiency and the path-dependent informational geometry measures. Section 3 provides our numerical methods and key diagnostics. Section 4 and Section 5 provide results and discussions. We conclude in Section 6. We provide general thermodynamic relations in Appendix A and Appendix B to make the paper self-contained.

2. Model

To gain key insight, we consider a rocking ratchet model [8,9,10] governed by the following Langevin equation
d x d t = f + ξ = V x + ξ = x V 0 + G ( t ) + F + ξ .
Here, f = x V is the force due to the potential V ( x , t ) . The potential V ( x , t ) = V 0 ( x ) + V 1 ( x , t ) + V 2 ( x ) contains the three parts. V 0 ( x ) = V 0 ( x + L ) is a spatially periodic potential with the period L, V 1 ( t ) = x G ( t ) is a time-dependent (rocking) potential given by a function G ( t ) = G ( t + T ) that is periodic in time with period T, V 2 = x F is a potential due to a constant force F. We note that T in this paper denotes a temporal periodicity of G and not temperature.
ξ is assumed to have a zero mean ξ = 0 and the strength D with the following property
ξ ( t ) ξ ( t ) = 2 D δ ( t t ) ,
where the angular brackets denote the ensemble average over ξ . We note that we are using the notation D for the fluctuation in heat bath (temperature) instead of T (the periodicity of the force G ( t ) ).
The Fokker–Planck equation [39] for a time-dependent probability density function (PDF) p ( x , t ) corresponding to Equations (1) and (2) is
t p ( x , t ) = J x = x V x p ( x , t ) + D x p ( x , t )
where V = V 0 ( x ) x G ( t ) x F ; J = V x p ( x , t ) + D x p ( x , t ) = f p D x p = ( x V 0 + G ( t ) + F ) p D x p is the probability current.
Because of the spatial periodicity of V 0 ( x ) = V 0 ( x + L ) , we have p ( x , t ) = p ( x + L , t ) and J ( x , t ) = J ( x + L , t ) . So, we normalise the PDF p ( x , t ) over L as 0 L d x p ( x , t ) = 1 at any time t. The ensemble average of a variable, say, A ( x , t ) , is then expressed in terms of the PDF of A ( x , t ) as
A ( t ) = 0 L d x A ( x , t ) p ( x , t ) .
We use the double angular brackets to denote the average over the space and time, for instance
A = 1 T 0 T d t 0 L d x A ( x , t ) p ( x , t ) .
For the current J, Equations (4) and (5) translate to
J = 1 L 0 L d x J ( x , t ) ,
J = 1 L T 0 T d t 0 L d x J ( x , t ) .

2.1. Thermodynamic Relations

Thermodynamic laws concern the energy conservation and entropy relations and are provided in Appendix A and Appendix B for a general overdamped process. When V = V 0 ( x ) x G ( t ) x F in Equation (1) [ V 1 ( x , t ) = x G ( t ) , V 2 ( x ) = x F ], the energy conservation in Equation (A3) (with W ˙ = E ˙ i n ) can be expressed as
d d t U 0 + U 1 + U 2 = d d t V E ˙ i n Q ˙ .
Here, U 0 = 0 L d x p ( x , t ) V 0 , U 1 = 0 L d x p ( x , t ) ( x G ( t ) ) , and U 2 = 0 L d x p ( x , t ) x F , and E ˙ i n and Q ˙ are given by
E ˙ i n = 0 L d x ( t V ) p = t V 1 = 0 L d x x G ˙ ( t ) p = x G ˙ ,
Q ˙ = 0 L d x V ( t p ) = 0 L d x J f = f x ˙ = D S ˙ M .
Here, E ˙ i n is the rate of energy input (e.g., by chemical agency, ATP, etc.) and Q ˙ the heat flow from the system to the environment. The rate at which a system does work against the external force is defined by
E ˙ o = d U 2 d t = d d t 0 L p ( x , t ) V 2 = F 0 L d x ( t p ) x = F 0 L d x ( x J ) x = F 0 L d x J .
Using Equations (9)–(11) in Equation (8) gives us
d d t U 0 + U 1 + Q ˙ = E ˙ i n E ˙ o .
The total energy input over one period T is given by the time-integral of Equation (9) as
E i n = t t + T d t 1 0 L d x x G ˙ ( t 1 ) p ( x , t 1 ) = t t + T d t 1 0 L d x x G ( t 1 ) p ˙ ( x , t 1 ) = t t + T d t 1 0 L d x x G ( t 1 ) x J ( x , t 1 ) = t t + T d t 1 0 L d x G ( t 1 ) J ( x , t 1 ) .
Furthermore, over one period T, 0 T d t 1 d d t 1 ( U 0 ( t 1 ) + U 1 ( t 1 ) ) = 0 . Thus, Equation (12) is simplified as
Q = E i n E o .
Q in Equation (14) is related to the entropy flow S M from the system to the environment as Q = D S M . It is also related to the entropy production S T = S + S M 0 in Equation (A1) (see Appendix A for details), where S = 0 L d x p ( x , t ) ln p ( x , t ) is the differential entropy. Specifically, Equation (A2) reads
S ˙ T = d S T d t = 0 L d x 1 D p J 2 , S ˙ M = d S M d t = 0 L d x 1 D J f ,
where S ˙ T and S ˙ M denotes the entropy production rate and entropy flow rate, respectively. Equation (15) shows that S ˙ T 0 ; S ˙ T > 0 serves as a measure of irreversibility far from equilibrium. S ˙ M is positive when the entropy flows from the system to the environment. Due to the periodicity of p ( x , t ) = p ( x , t + T ) , the total change over the cycle in the differential entropy S = 0 T d t 1 S ˙ ( t 1 ) = 0 . Thus, we have S T = S + S M = S M = Q / D 0 where
S T = 0 T d t S ˙ T ( t ) , S M = 0 T d t S ˙ M ( t ) .

2.2. Efficiency

Using the relations in Section 2.1, we have the instantaneous efficiency ( η ) and cumulative efficiency ( η c ) defined over one cycle as follows:
η = E ˙ o E ˙ i n ,
η c = E o E i n = 1 Q E i n = 1 D S T E i n .
Equation (18) shows explicitly that the non-zero entropy production S T > 0 or heat flow Q reduces the efficiency.
To take into account the viscous work E ˙ s = γ x ˙ 2 = γ d x J 2 ( γ = 1 for Equation (1)), Equations (17) and (18) were generalised to include the Stokes efficiency [40,41] as follows
η = E ˙ o + E s ˙ E ˙ i n = E s ˙ E ˙ i n ,
η c = E s E i n ,
where E s = t t + T d t 1 E ˙ s ( t 1 ) .
We note that even when there is no external force F = 0 , the system can generate a non-zero current J > 0 . A stopping force F s < 0 is the force that is required to make the current zero J = 0 above/below which the sign of the current changes ( J > 0 for F > F s and J < 0 for F < F s ). To separate the additional effect of F 0 , we will focus on F = 0 in Section 3, Section 4 and Section 5. Moreover, since in our problem (e.g., for a sinusoidal G ( t ) ) E ˙ i n can take zero value at certain time with a singular η , we will calculate η c in Equation (20) in Section 3, Section 4 and Section 5.

2.3. Information Rate and Length

Information rate Γ and length L are the information geometric measures that quantify how information unfolds in time-varying stochastic processes [18,27,28,29,30,31,32,33,34,35,36,37,38]. In a nutshell, Γ is proportional to the square root of infinitesimal (symmetric) relative entropy (Kullback–Leibler divergence) while L measures the total change in the information along its evolution path. For example, for a PDF of one variable x evolving in time t,
Γ 2 ( t ) = d x p ( x , t ) ( t ln p ( x , t ) ) 2 = 4 d x ( t q ( t ) ) 2 ,
L ( t ) = 0 t d t 1 Γ ( t 1 ) .
We note that Γ 2 in terms of q = p is well defined for p = 0 . When the parameters of p ( x , t ) are known as λ i ’s ( i = 1 , 2 , ), Γ 2 can be expressed in terms of the metric tensor g i j = d x p ( x , t ) ( λ i ln p ( x , t ) ) ( λ j ln p ( x , t ) ) as Γ 2 = i , j ( t λ i ) g i j ( t λ j ) .
Γ 1 has the dimensions of time and is linked to the smallest timescale of fluctuations [38]. As a non-decreasing function of time, L ( t ) is a dimensionless, path-dependent distance measuring the deviation from the initial state in terms of the total number of statistically different states that a system passes through. For a Gaussian PDF with a constant variance, Δ L = 1 when p ( x , t ) moves from p ( x , 0 ) by one standard deviation since the latter provides the uncertainty in measuring the PDF position. Even when p ( x , 0 ) = p ( x , T ) , Γ ( t ) is sensitive to temporal changes at 0 < t < T and take different values depending on the path. Γ ( t ) and L ( t ) are invariant under (time-independent) change of variables (unlike entropy), enabling us to compare the evolution of different variables.
We measure the change in information associated with the dynamics of a Brownian ratchet by calculating
L = 0 T d t 1 Γ ( t 1 ) ,
over the time interval of one period. Equation (23) represents the total number of different statistical states that a ratchet passes through in one period.

3. Numerical Models and Methods

3.1. Potential V 0 ( x ) and Rocking Force G ( t )

As noted in Section 1, we consider a sawtooth potential V 0 ( x ) , as shown in the first panel of Figure 1. Note first how V 0 ( x ) is periodic in x, with period 1. Any other period L could always be rescaled to 1 by introducing a new variable x ˜ = x / L , so there is no loss of generality in only considering L = 1 . V 0 then has maxima at x = 0 mod ( 1 ) and minima at x = 0.2 mod ( 1 ) . The particular value x = 0.2 as the location of the minimum is not critical, as long as it is not at 0.5, which would yield a potential that was symmetric under x x . In contrast, this potential V 0 ( x ) is clearly asymmetric under x x , and this asymmetry underlies the entire rocking ratchet problem.
In particular, consider what happens if V 0 ( x ) is rocked back and forth, by the addition of the term V 1 ( t ) = x G ( t ) . Suppose that the maximum positive and negative amplitudes of G ( t ) are ± A , for some constant A. The second panel in Figure 1 shows the results for A = 0.5 : the original V 0 ( x ) is simply rocked back and forth. Now, consider what would happen to any particles situated in the V 0 potential wells. For this particular rocking amplitude of 0.5 , nothing would happen yet (at least in the absence of any noise), because the local minima always remain local minima throughout the entire rocking process. However, by considering the slopes of V 0 ( x ) on its two segments, that is, 5 on x [ 0 , 0.2 ] and 1.25 on x [ 0.2 , 1 ] , we see that if A had been chosen in the range 1.25 < A < 5 , any particles in the potential wells would always spill over to the right at that point in the rocking cycle but would never spill back to the left even at the opposite point in the cycle. This is the essence of the ratchet mechanism; the original asymmetry in V 0 ( x ) causes this spill-over effect to be uni-directional, resulting in a directed particle current. The interesting dynamics about this problem are then all about what happens in the presence of noise, how do the resulting current and other diagnostics vary with the amplitude and period of the rocking, etc.
There are various possible choices for the time-dependence G ( t ) of the rocking. One simple choice is sinusoidal,
G ( t ) = A sin ( ω t ) ,
where ω = 2 π / T relates the frequency ω to the period T. Another choice is a square wave, with
G ( t ) = { A 0 t < T / 2 mod ( T ) , A T / 2 t < T mod ( T ) .
A third choice is a sawtooth wave, with
G ( t ) = A ( 1 2 t / T ) 0 t < T mod ( T ) ,
consisting of a linear variation from A at t = 0 to A at t = T , followed by a discontinuous jump back to A for the next period, etc.
Note also that positive versus negative A in Equation (26) are not equivalent—unlike in Equation (24) or (25), where changing the sign of A merely shifts G ( t ) by half a period in time. In contrast, in Equation (26) positive and negative A are different, in the following sense: taking A positive starts with the upward sloping potential in Figure 1b, gradually tilts that to the downward sloping potential, then discontinuously jumps back to the upward sloping potential, etc. Because the original V 0 ( x ) already introduced the ± x asymmetry, this is different from what a negative A would do, namely to have the gradual tilting going from downward to upward sloping.
Since we prefer to always think of the amplitudes A as being positive though, we modify our definition of the sawtooth wave to be
G ± ( t ) = ± A ( 1 2 t / T ) 0 t < T mod ( T ) ,
so that now we have two different sawtooth waves, but each one separately can be restricted to considering only A > 0 , just as for the sinusoidal and square wave choices. We will see then that the two choices G ± yield results that are very similar in some diagnostic quantities but different in others.

3.2. Numerical Methods

To numerically solve the Fokker–Planck Equation (3), with a spatially periodic potential, we use the following
p ( x , t ) = a 0 ( t ) + k = 1 K a k ( t ) cos ( 2 π k x ) + b k ( t ) sin ( 2 π k x ) .
See also [42] for the theory underlying Fourier series in numerical analysis. The coefficients a k and b k are then time-stepped using second-order Runge–Kutta [43]. The term x ( V x p ) is separated out into the relevant Fourier components using a fast Fourier transform. Since our chosen potential V 0 ( x ) has discontinuous derivatives, special attention is required to ensure that results are adequately resolved.
Resolutions in the range K = 2 12 2 20 were used and time-steps in the range 10 4 10 6 . One very basic test of the numerical implementation is simply to monitor the coefficient a 0 : this is time-stepped along with the others but must, in fact, remain constant if the total probability p d x = a 0 is to remain constant. It was found that if the initial condition is set to have a 0 = 1 , then this was indeed preserved throughout the subsequent evolution. Note that there are also other numerical methods of solving Fokker–Planck equations that are very different from the direct time-stepping implemented here, instead using a variety of neural network approaches (e.g., [44,45,46]).
Since we are not interested in initial transients, but only in the final solutions p ( x , t ) that are periodic in both x (with period 1) and t (with period T), we start with the very simple initial condition a 0 = 1 and all other Fourier coefficients set to zero, and integrate for however many periods T are required to settle in to a solution that is periodic in time. For some parameter choices this happened after just a few periods, but for other parameter choices many dozens of periods were required, significantly increasing the total computational effort.
As noted above, throughout this entire evolution the solutions always maintained a 0 = 1 , but this is not a sufficiently stringent check to convince ourselves that the numerics are adequately dealing with the rather challenging potential V 0 ( x ) . To see how the code is coping with this, Figure 2 shows some sample results for the sinusoidal choice Equation (24) with A = 6 and two noise levels D as indicated in the figure. Note that, in particular, according to our arguments above, this A value is sufficiently large that spill-over should occur in both directions at different points in the rocking cycle. Sure enough, at different times the solutions are concentrated at locations far away from the x = 0.2 potential well of just V 0 ( x ) alone.
If we focus on the PDFs at t = 0 mod ( T ) though, they do indeed form a peak around x = 0.2 , with the width also becoming narrower as D is decreased. The first panel in Figure 3 examines this in more detail, showing how the width scales as D and the height as 1 / D . Note also how the solutions appear to be forming true cusps at x = 0.2 , with p continuous but x p discontinuous at x = 0.2 . This is exactly what one would expect, based on the discontinuity in d d x V 0 ( x ) . The second panel in Figure 3 shows the corresponding Fourier power spectra a k 2 + b k 2 . Note how they are flat up to a certain k, and then start dropping off. Decreasing D by an order of magnitude increases the transition value of k by an order of magnitude, corresponding to the width scaling as D. Finally, once the wave-numbers k are sufficiently large to be well in the drop-off regime, the power scaling is as k 4 . From the general theory of Fourier series, this is exactly what one expects from a function that is continuous but with discontinuities in its first derivative. We see therefore that this potential V 0 ( x ) is challenging, especially at very small D, but that if one is prepared to set the Fourier series resolution K sufficiently high, any arbitrarily small D can still be achieved and yields a well-resolved power spectrum.
In the following section, we will present systematic scans over a range of A values for D as small as 10 3 , which is already small enough for clear trends to emerge. We will focus in particular on the diagnostic quantities current (Equation (7)), efficiency (Equation (20)), L (Equation (23)) and S T (Equation (16)) integrated over the period T and then mapped out as functions of A, as well as some representative cases of current (Equation (6)), Γ (Equation (21)), S ˙ (Equation (A1)) and S ˙ T (Equation (15)) as functions of time over a period.

4. Results for G ( t ) = A sin ( ω t )

We start with the sinusoidal rocking Equation (24). We investigated a variety of periods T, ultimately focusing on the four values T = 0.4 , 0.7 , 1 and 2 to consider in detail. For each of these, we scanned over the range A [ 0 , 20 ] , for the three noise values D = 10 1 , 10 2 and 10 3 .
Figure 4 shows the current J from Equation (7). One feature that immediately stands out is how J tends to zero if A is too small, consistent with our discussion above regarding the dynamics of the ratchet mechanism, and why there is a non-zero current at all. It is also noticeable that J is much less for T = 0.4 than for the larger values. This is also as one might expect: if the system is rocked back and forth too rapidly, there simply is not enough time for much spill-over to occur, even if A is large enough that it otherwise would. For T even smaller than 0.4, the current also becomes even smaller.
Regarding the variation with D, we note how the various peaks become increasingly sharp and distinct as D is reduced. The limit D 0 would consist of a set of peaks that vary discontinuously with A, being either zero or 1 / n T , for n = 1 , 2 , (e.g., see [47]). To understand such a complicated pattern, and why the current does not vary monotonically with A, we recall that once A > 5 for our particular V 0 here, the spill-over can also go to the left, the opposite of the ‘desired’ direction. When J abruptly drops to zero are then values of A where this negative spill-over occurred, whereas J jumps back to some non-zero value for those values of A where the positive spill-over extends one further well of V 0 than before.
Figure 5 shows the Stokes efficiency η c from Equation (20). If A is too small it is relatively small, but for larger A it can approach 1, even for T = 0.4 where the current was small. As with the current, there are various up-and-down jumps that become sharper and more distinct for smaller D. Comparing Figure 4 and Figure 5, one can already see that the jumps in J and η c occur at similar values of A; we will consider this phase relationship between these quantities in more detail below (Figure 8).
Figure 6 shows the information length L from Equation (23). The two features that immediately stand out are that L generally increases with increasing A, and decreasing D. This suggests to collapse the results by instead plotting L · D 1 / 2 / A . As shown in Figure 7, this does indeed remove the overall increase with A, allowing one to focus more on the jumps that are similar to those seen before in the current and efficiency. Figure 7 also suggests that in the D 0 limit, this combination might become independent of D, including the details of the various jumps. However, D would have to be reduced quite a bit further to fully clarify this question. The overall trend though that L is broadly proportional to D 1 / 2 is very likely robust. Roughly speaking, this is (i) because L is a dimensionless distance that measures the change in the mean location of a PDF (mean value) with respect to the measurement error given by the standard deviation (width of a PDF) as long as the mean (A) is not too small compared with the standard deviation and (ii) the change in the mean and standard deviation increase with A and D 1 / 2 , respectively.
The three quantities presented so far, J , η c and L D 1 / 2 / A , all exhibit similar up and down jumps as A is increased. To clarify the phase relationships between these quantities, Figure 8 super-imposes the three, for the D = 10 3 case where the jumps become clearest. The range in A is also restricted to A [ 3 , 15 ] to focus on the first few jumps. Only T = 0.7 and 1 are shown here, but 0.4 and 2 exhibited similar behaviour. In particular, we see how the jumps in all three quantities are indeed very closely correlated, with maxima or minima in η c and L D 1 / 2 / A tending to line up with those A values where J is changing most rapidly. This suggests that our (normalised) information geometric diagnostics L D 1 / 2 / A approximates the Stokes efficiency better than the current. We however notice that the phase relations seem to be less clear for larger A 10 .
Based on Figure 8, the particular values A = 6 , 7 , 8 for T = 0.7 represent one of the transitions from J being large at A = 6 to essentially zero at A = 8 , and correspondingly A = 7 is a minimum for L . The first two panels in Figure 9 then show the variation in time of the quantities J and Γ . Note how J is positive for almost the entire cycle for A = 6 , which is why J is so large at this A. For A = 7 and 8, J is if anything slightly greater for the first half of the cycle, but for the second half it is strongly negative. This is precisely the point noted above, that for sufficiently large A the spill-over can also go in the ‘wrong’ direction, reducing the average current J to essentially zero at A = 8 .
In comparison, Γ measures the rate at which the change in a PDF occurs, taking a large value when a PDF changes rapidly. Thus, a spill-over regardless of its direction can cause a sudden change in a PDF shape. This is seen in Figure 9 where Γ tends to be large at points in the cycle where the absolute value of J is also large. That is, Γ is insensitive to the sign of any spill-over, but it does measure that spill-over is occurring. Furthermore, the first peak of J occurs at different times t 0.2 , 0.25 for A = 6 , 7 , 8 , while the peak of Γ occurs at a similar time t 0.12 , earlier than that of J . This suggests that Γ detects a rapid change in a PDF prior to the appearance of J peak.
Panels (c,d) in Figure 9 show the variation in time of the quantities S ˙ and S ˙ T = S ˙ + S ˙ M . A few points to note here are: First, S ˙ can be both positive and negative, see Equation (A1), and indeed integrates to zero due to the temporal periodicity S ( t ) = S ( t + T ) , as noted in Section 2.2. In contrast, S ˙ T is strictly non-negative, as it must be according to Equations (15) and (16). Next, we can note that a typical magnitude of S ˙ is roughly 50 times smaller than a typical magnitude of S ˙ T , O ( 100 ) versus O ( 50 , 000 ) . Finally, if we compare S ˙ with J in the panel directly above it, we notice that the jumps in S ˙ almost invariably line up with certain features in J , for each of the three values A = 6 , 7 , 8 . It seems therefore that S ˙ is also sensitive to the spill-over dynamics that cause variations in J . This is because the spill-over inevitably causes the change in a PDF such as broadening or narrowing of a PDF width and thus entropy.
Figure 10 shows the period-integrated quantity D · S T = D · S M , for the same parameter values as before in Figure 4, Figure 5, Figure 6 and Figure 7. There are again three curves in each panel, with the same colour-coding as before, but we see that the three curves almost completely overlap virtually everywhere. That is, the scaling S T D 1 holds so accurately that D S T becomes almost completely identical for the three different D values. The other point to note is how the curves in all four panels are following parabolic profiles for all but the very smallest values of A. That is, once A > 5 or so, S T A 2 becomes an almost exact fit.

5. Results for Square and Sawtooth G ( t )

Similar scans were performed for the square wave Equation (25) and sawtooth waves Equation (27). We will only present T = 1 here; other values were also considered and exhibited similar trends. Figure 11 again shows the current J as a function of A, for the same three D values as before. We note how the maximum value of J for the square wave is twice what it was for the sinusoidal G ( t ) and also for the sawtooth waves. This simply reflects the fact that the square wave G ( t ) = ± A has twice the rms value as the sinusoidal G ( t ) . Another interesting point is that J for the square wave does not seem to drop back to zero between the various peaks, unlike for the sinusoidal or sawtooth waves. Finally, note how the two sawtooth waves G ± yield virtually identical results for this particular diagnostic quantity J .
Figure 12 shows the Stokes efficiency η c . The results are again broadly similar to the sinusoidal results in Figure 5, with efficiencies approaching one in all cases. Again, just as before the particularly small currents for the sinusoidal G ( t ) with T = 0.4 had no particular effect on the efficiency, here also the particularly large currents for the square wave have no discernible effect on its efficiency. Another point to note here is that the results for the two sawtooth waves are still similar, but noticeably different, unlike in Figure 11 where the currents were virtually identical.
Finally, for the information length L , we dispense with the plots of L itself and proceed directly to plots of L D 1 / 2 / A , as shown in Figure 13. The overall pattern is again similar to what it was for the sinusoidal G ( t ) . In particular, there is again at least the suggestion that this combination might be tending to a limit independent of D, but considerably smaller values would be required to confirm this. The last point to note is that now the results for the two sawtooth waves G ± are significantly different, more even than previously for the Stokes efficiency. Information length is evidently a very delicate diagnostic quantity that picks out differences in PDFs that other diagnostics are far less sensitive to.

6. Conclusions

We investigated exact time-dependent solutions of a Brownian motor under different conditions to elucidate the role of symmetry-breaking and information geometry. Specifically, we utilised an asymmetric sawtooth periodic potential and three different types of periodic forcing G ( t ) given by sinusoidal, square and sawtooth waves with period T and amplitude A. For each G ( t ) , we calculated exact time-dependent probability density functions (PDFs) by numerically solving the Fokker–Planck equation [39] under different conditions by varying T, A and the strength of the stochastic noise D in an unprecedentedly wide range. We showed that a non-differentiable potential led to a non-differentiable PDF and ensured that our solutions were well resolved by including a sufficiently large number of Fourier modes. From our time-dependent PDF, we performed a systematic investigation of comparing mean current, Stokes efficiency, entropy production and the path-dependent information geometric measures (information rate and information length). Interestingly, regardless of different forms of G ( t ) , we found overall similar behaviours. This highlights the essential role of symmetry-breaking and robustness of the results that are less sensitive to the detailed temporal evolution of a motor (driven by different forcings G ( t ) ).
In particular, current, Stokes efficiency and the information rate normalised by A and D exhibit one or multiple local maxima and minima as A increases. However, the dependence of current and Stokes efficiency on A can be quite different, while the information rate normalised by A and D tends to resemble that of the Stokes efficiency. In comparison, the irreversibility measured by a normalised entropy production is independent of A. The results thus suggest that the entropy production is not a good measure of a motor efficiency. Instead, our information geometry provides a useful proxy of the Stokes efficiency. Whether this result holds in general when there is a non-zero constant external force F in our model Equation (1) will be investigated in future. Future work will also address different types of Brownian motors such as the (two-state) flashing ratchet [48] driven by the dichotomous noise in addition to a short-correlated noise or Levy-type noise. It will also be of interest to compare our path-dependent information geometric theory with other measures of irreversibility such as the Komogorov eigenvalue or non-equilibrium Lyapunov function [49].

Author Contributions

Original research idea, E.-j.K.; investigation, E.-j.K. and R.H.; software, R.H.; visualisation, R.H.; writing—original draft, E.-j.K. and R.H. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

Not applicable.

Data Availability Statement

Data is contained within the article. Further details, e.g., regarding the numerical code, are available from the authors on reasonable request.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A. Entropy Production Rate and Flow

Entropy measures disorder or the lack of information [17]. For a continuous PDF p ( x , t ) , a differential entropy S is defined by S ( t ) = d x p ( x , t ) ln ( p ( x , t ) ) . By using Equation (8) [ p ˙ = x J , J = f p D x p , f = x V ], the periodic boundary condition over x = [ 0 , L ] and the integration by parts, we express the time derivative of entropy S ˙ = d S ( x , t ) d t as
S ˙ = 0 L d x p ˙ ln p = 0 L d x [ x J ln p ] = 0 L d x [ J x ln p ] = S ˙ T S ˙ M .
Here,
S ˙ T = d S T d t = 0 L d x 1 D p J 2 , S ˙ M = d S M d t = 0 L d x 1 D J f .
S ˙ T 0 represents entropy production rate (which is non-negative as can be seen from Equation (A2)) and is a measure of irreversibility far from equilibrium. S M is the entropy flux to the environment, e.g., being positive when the entropy flows from the system to the environment. In an equilibrium reversible process S ˙ T = 0 , Equation (A1) establishes an equilibrium thermodynamic relation between S ˙ and Q ˙ = D S ˙ M (see Appendix B) as S ˙ = S ˙ M = Q / D . On the other hand, if S ˙ = 0 , S ˙ T = S ˙ M .

Appendix B. Energetics: First Thermodynamic Law

We summarise the energetics that links the entropy relation in Appendix A to the thermodynamics. For the average potential energy U = V , we have
d U d t = d d t V W ˙ Q ˙ .
Here, W ˙ represents the power (the rate of the work) due to the explicit time-dependence of V; Q ˙ is the dissipated heat defined by
W ˙ = 0 L d x ( t V ) p = t V ,
Q ˙ = 0 L d x V ( t p ) = 0 L d x J f = f x ˙ = D S ˙ M .
Equation (A3) gives our non-equilibrium thermodynamic relation U = W Q that the energy function increases due to W ˙ while decreasing for the total dissipated heat to the environment Q ˙ . Using Equations (A3), (A5), and (A1) in a non-equilibrium (information) free energy F ( t ) = U ( t ) D S ( t ) [50] gives us
F ˙ = U ˙ D S ˙ = W ˙ D S ˙ T .
Using that S ˙ T 0 , we then rewrite Equation (A6) as
D S ˙ T = W ˙ F ˙ W ˙ D 0
where W ˙ D is the dissipated power which is lost to the environment. The time-integral of Equation (A7) for time [ t 0 , t ] gives us the well-known result that the average work performed on the system is bounded below by the free energy difference as W Δ F = W D 0 . For molecular motors, the energy input E ˙ i n is used for W ˙ .

References

  1. Parrondo, J.M.R.; De Cisneros, B.J.; Brito, R. Energetics of Brownian motors: A review. Appl. Phys. A 2002, 75, 179–191. [Google Scholar] [CrossRef]
  2. Reimann, P. Brownian motors: Noisy transport far from equilibrium. Phys. Rep. 2002, 361, 57. [Google Scholar] [CrossRef] [Green Version]
  3. Hänggi, P.; Marchesoni, F. Artificial Brownian motors: Controlling transport on the nanoscale. Rev. Mod. Phys. 2009, 81, 387–442. [Google Scholar] [CrossRef] [Green Version]
  4. Xing, J.; Wang, H.; Oster, G. From continuum Fokker-Planck models to discrete kinetic models. Biophys. J. 2005, 89, 1551. [Google Scholar] [CrossRef] [Green Version]
  5. Julicher, F.; Ajdari, A.; Prost, P. Modeling molecular motors. Rev. Mod. Phys. 1997, 69, 1269. [Google Scholar] [CrossRef] [Green Version]
  6. Astumian, R.D.; Derényi, I. Fluctuation driven transport and models of molecular motors and pumps. Eur. Biophys. J. 1998, 27, 474–489. [Google Scholar] [CrossRef] [PubMed]
  7. Warshaw, D.W. Throttling back the heart’s molecular motor. Science 2016, 351, 556. [Google Scholar] [CrossRef] [Green Version]
  8. Magnasco, M.O. Forced thermal ratchet. Phys. Rev. Lett. 1993, 71, 1477–1481. [Google Scholar] [CrossRef]
  9. Kamegawa, H.; Hondou, T.; Takagi, F. Energetics of forced thermal ratchet. Phys. Rev. Lett. 1998, 80, 5251. [Google Scholar] [CrossRef] [Green Version]
  10. Sumithra, K.; Sintes, T. Efficiency optimization in forced ratchets due to thermal fluctuations. Physica A 2001, 297, 1–12. [Google Scholar] [CrossRef] [Green Version]
  11. Ai, B.-Q.; Wang, X.-J.; Liu, G.-T.; Wen, D.-H.; Xie, H.-Z.; Chen, W.; Liu, L.-G. Efficiency optimization in a correlation ratchet with asymmetric unbiased fluctuations. Phys. Rev. E 2003, 68, 061105. [Google Scholar]
  12. Ai, B.-Q.; Wang, L.; Liu, L.-G. Flashing motor at high transition rate. Chaos Solitons Frac. 2003, 34, 1265–1271. [Google Scholar] [CrossRef]
  13. Murashita, Y. Absolute Irreversibility in Information Thermodynamics. Master’s Thesis, Department of Physics, University of Tokyo, Tokyo, Japan, 2015. Available online: https://arxiv.org/pdf/1506.04470.pdf (accessed on 2 May 2021).
  14. Seifert, U. Stochastic thermodynamics, fluctuation theorems and molecular machines. Rep. Prog. Phys. 2012, 75, 126001. [Google Scholar] [CrossRef] [Green Version]
  15. Gupta, D. Exact distribution for work and stochastic efficiency of an isothermal machine. J. Stat. Mech. 2018, 2018, 073201. [Google Scholar] [CrossRef] [Green Version]
  16. Haken, H. Information and Self-Organization: A Macroscopic Approach to Complex Systems, 3rd ed.; Springer: Berlin/Heidelberg, Germany, 2006. [Google Scholar]
  17. Frieden, B.R. Science from Fisher Information; Cambridge University Press: Cambridge, UK, 2004. [Google Scholar]
  18. Kim, E. Information geometry, fluctuations, non-equilibrium thermodynamics, and geodesics in complex systems. Entropy 2021, 23, 1393. [Google Scholar] [CrossRef]
  19. Landauer, R. Irreversibility and heat generation in the computing process. IBM J. Res. Dev. 1961, 5, 183–191. [Google Scholar] [CrossRef]
  20. Bérut, A.; Arakelyan, A.; Petrosyan, A.; Ciliberto, S.; Dillenschneider, R.; Lutz, E. Experimental verification of Landauer’s principle linking information and thermodynamics. Nature 2012, 483, 187. [Google Scholar] [CrossRef]
  21. Leff, H.S.; Rex, A.F. Maxwell’s Demon: Entropy, Information, Computing; Princeton University Press: Princeton, NJ, USA, 1990. [Google Scholar]
  22. Toyabe, S.; Sagawa, T.; Ueda, M.; Muneyuki, E.; Sano, M. Experimental demonstration of information-to-energy conversion and validation of the generalized Jarzynski equality. Nat. Phys. 2010, 6, 988–992. [Google Scholar] [CrossRef] [Green Version]
  23. Lee, D.Y.; Um, J.; Paneru, G. An experimentally-achieved information-driven Brownian motor shows maximum power at the relaxation time. Sci. Rep. 2018, 8, 12121. [Google Scholar] [CrossRef]
  24. Kowalski, A.M.; Martin, M.T.; Plastino, A.; Rosso, O.A.; Casas, M. Distances in probability space and the statistical complexity setup. Entropy 2011, 13, 1055. [Google Scholar] [CrossRef] [Green Version]
  25. Parr, T.; Da Costa, L.; Friston, K.J. Markov blankets, information geometry and stochastic thermodynamics. Phil. Trans. Roy. Soc. A 2019, 378, 20190159. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  26. Oizumi, M.; Tsuchiya, N.; Amari, S. Unified framework for information integration based on information geometry. Proc. Nat. Acad. Sci. USA 2016, 113, 14817. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  27. Guel-Cortez, A.J.; Kim, E. Information geometric theory in the prediction of abrupt changes in system dynamics. Entropy 2021, 23, 694. [Google Scholar] [CrossRef]
  28. Kim, E.; Hollerbach, R. Geometric structure and information change in phase transitions. Phys. Rev. E 2017, 95, 062107. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  29. Kim, E.; Jacquet, Q.; Hollerbach, R. Information geometry in a reduced model of self-organised shear flows without the uniform coloured noise approximation. J. Stat. Mech. 2019, 2019, 023204. [Google Scholar] [CrossRef] [Green Version]
  30. Kim, E.; Heseltine, J.; Liu, H. Information length as a useful index to understand variability in the global circulation. Mathematics 2020, 8, 299. [Google Scholar] [CrossRef] [Green Version]
  31. Kim, E.; Lee, U.; Heseltine, J.; Hollerbach, R. Geometric structure and geodesic in a solvable model of nonequilibrium process. Phys. Rev. E 2016, 93, 062127. [Google Scholar] [CrossRef] [Green Version]
  32. Heseltine, J.; Kim, E. Comparing information metrics for a coupled Ornstein-Uhlenbeck process. Entropy 2019, 21, 775. [Google Scholar] [CrossRef] [Green Version]
  33. Kim, E.; Hollerbach, R. Time-dependent probability density functions and information geometry of the low-to-high confinement transition in fusion plasma. Phys. Rev. Res. 2020, 2, 023077. [Google Scholar] [CrossRef] [Green Version]
  34. Hollerbach, R.; Kim, E.; Schmitz, L. Time-dependent probability density functions and information diagnostics in forward and backward processes in a stochastic prey-predator model of fusion plasmas. Phys. Plasmas 2020, 27, 102301. [Google Scholar] [CrossRef]
  35. Guel-Cortez, A.J.; Kim, E. Information length analysis of linear autonomous stochastic processes. Entropy 2020, 22, 1265. [Google Scholar] [CrossRef] [PubMed]
  36. Kim, E. Investigating information geometry in classical and quantum systems through information length. Entropy 2018, 20, 574. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  37. Kim, E. Information geometry and non-equilibrium thermodynamic relations in the over-damped stochastic processes. J. Stat. Mech. 2021, 2021, 093406. [Google Scholar] [CrossRef]
  38. Nicholson, S.B.; García-Pintos, L.P.; del Campo, A.; Green, J.R. Time-information uncertainty relations in thermodynamics. Nat. Phys. 2020, 16, 1211–1215. [Google Scholar] [CrossRef]
  39. Risken, H. The Fokker-Planck Equation: Methods of Solution and Applications; Springer: Berlin/Heidelberg, Germany, 1996. [Google Scholar]
  40. Derényi, I.; Bier, M.; Astumian, R.D. Generalized efficiency and its application to microscopic engines. Phys. Rev. Lett. 1999, 83, 903. [Google Scholar] [CrossRef] [Green Version]
  41. Mamata, S.; Jayannavar, A.M. Stokes efficiency and its stochastic properties. Phys. A 2017, 465, 40–48. [Google Scholar]
  42. Boyd, J.P. Chebyshev and Fourier Spectral Methods, 2nd ed.; Dover: Mineola, NY, USA, 2013. [Google Scholar]
  43. Gottlieb, S.; Ketcheson, D.; Shu, C.-W. Strong Stability Preserving Runge-Kutta and Multistep Time Discretizations; World Scientific: Singapore, 2011. [Google Scholar]
  44. Xu, Y.; Zhang, H.; Li, Y.; Zhou, K.; Liu, Q.; Kurths, J. Solving Fokker-Planck equation using deep learning. Chaos 2020, 30, 013133. [Google Scholar] [CrossRef]
  45. Chen, X.; Yang, L.; Duan, J.; Karniadakis, G.E. Solving inverse stochastic problems from discrete particle observations using the Fokker-Planck equation and physics-informed neural networks. SIAM J. Sci. Comput. 2021, 43, B811–B830. [Google Scholar] [CrossRef]
  46. Tang, K.; Wan, X.; Liao, Q. Adaptive deep density approximation for Fokker-Planck equations. arXiv 2021, arXiv:2103.11181. [Google Scholar]
  47. Bartussek, R.; Hänggi, H.; Kissner, J.G. Periodically rocked thermal ratchets. Europhys. Lett. 1994, 28, 459. [Google Scholar] [CrossRef]
  48. Jülicher, F.; Prost, J. Cooperative molecular motors. Phys. Rev. Lett. 1995, 75, 2618. [Google Scholar] [CrossRef] [PubMed]
  49. Caceres, M.O.; Lobos, A.M. Theory of eigenvalues for periodic non-stationary Markov processes: The Kolmogorov operator and its applications. J. Phys. A 2006, 39, 1547–1564. [Google Scholar] [CrossRef] [Green Version]
  50. Deffner, S.; Lutz, E. Information free energy for nonequilibrium states. arXiv 2012, arXiv:1201.3888. [Google Scholar]
Figure 1. (a) The potential V 0 ( x ) . Note how V 0 is periodic, with period 1, and has minima at x = 0.2 mod ( 1 ) . (b) V 0 is now rocked back and forth according to V 0 ( x ) ± 0.5 x , with the ± indicated by red/blue, respectively.
Figure 1. (a) The potential V 0 ( x ) . Note how V 0 is periodic, with period 1, and has minima at x = 0.2 mod ( 1 ) . (b) V 0 is now rocked back and forth according to V 0 ( x ) ± 0.5 x , with the ± indicated by red/blue, respectively.
Symmetry 14 00314 g001
Figure 2. The PDFs p ( x , t ) corresponding to G ( t ) = A sin ( ω t ) , with A = 6 and T = 0.7 . D = 10 1 and 10 2 in the top and bottom panels, as labelled. The times shown are t = n · T / 4 mod ( T ) , with n = 0 , 1 , 2 , 3 colour-coded as black, red, blue, green, respectively.
Figure 2. The PDFs p ( x , t ) corresponding to G ( t ) = A sin ( ω t ) , with A = 6 and T = 0.7 . D = 10 1 and 10 2 in the top and bottom panels, as labelled. The times shown are t = n · T / 4 mod ( T ) , with n = 0 , 1 , 2 , 3 colour-coded as black, red, blue, green, respectively.
Symmetry 14 00314 g002
Figure 3. (a) p ( x ) · D as a function of ( x 0.2 ) / D . Scaling p ( x ) and ( x 0.2 ) in these ways confirms the scalings of the height and the width. As in Figure 2, A = 6 and T = 0.7 , and t = 0 mod ( T ) , corresponding to the black curves in Figure 2. (b) The Fourier power spectra a k 2 + b k 2 of the solutions in (a). For both (a,b), D = 10 1 to 10 4 are colour-coded as blue, black, red, green, respectively.
Figure 3. (a) p ( x ) · D as a function of ( x 0.2 ) / D . Scaling p ( x ) and ( x 0.2 ) in these ways confirms the scalings of the height and the width. As in Figure 2, A = 6 and T = 0.7 , and t = 0 mod ( T ) , corresponding to the black curves in Figure 2. (b) The Fourier power spectra a k 2 + b k 2 of the solutions in (a). For both (a,b), D = 10 1 to 10 4 are colour-coded as blue, black, red, green, respectively.
Symmetry 14 00314 g003
Figure 4. J as functions of A. As indicated in the panels, T = 0.4 , 0.7 , 1 and 2. D = 10 1 , 10 2 and 10 3 are colour-coded as blue, black, red, respectively.
Figure 4. J as functions of A. As indicated in the panels, T = 0.4 , 0.7 , 1 and 2. D = 10 1 , 10 2 and 10 3 are colour-coded as blue, black, red, respectively.
Symmetry 14 00314 g004
Figure 5. The Stokes efficiency η c as functions of A. As in Figure 4, the four panels are at T = 0.4 , 0.7 , 1 and 2, as indicated, and D = 10 1 , 10 2 and 10 3 are colour-coded as blue, black, red, respectively.
Figure 5. The Stokes efficiency η c as functions of A. As in Figure 4, the four panels are at T = 0.4 , 0.7 , 1 and 2, as indicated, and D = 10 1 , 10 2 and 10 3 are colour-coded as blue, black, red, respectively.
Symmetry 14 00314 g005
Figure 6. L as functions of A. The periods T and the colour-coding are as in Figure 4 and Figure 5.
Figure 6. L as functions of A. The periods T and the colour-coding are as in Figure 4 and Figure 5.
Symmetry 14 00314 g006
Figure 7. The same data as in Figure 6, but now plotted as L D 1 / 2 / A versus A.
Figure 7. The same data as in Figure 6, but now plotted as L D 1 / 2 / A versus A.
Symmetry 14 00314 g007
Figure 8. The three quantities J in red, η c in black and L D 1 / 2 / A in blue, as functions of A. T = 0.7 and 1 in the (left) and (right) panels as labelled, and all results at D = 10 3 .
Figure 8. The three quantities J in red, η c in black and L D 1 / 2 / A in blue, as functions of A. T = 0.7 and 1 in the (left) and (right) panels as labelled, and all results at D = 10 3 .
Symmetry 14 00314 g008
Figure 9. (a) J as a function of time, for T = 0.7 , D = 10 3 and A = 6 (blue), A = 7 (red) and A = 8 (black). (b) Γ as a function of time, for the same parameters and colour-coding as in (a). (c) S ˙ as a function of time. (d) S ˙ T = S ˙ + S ˙ M as a function of time.
Figure 9. (a) J as a function of time, for T = 0.7 , D = 10 3 and A = 6 (blue), A = 7 (red) and A = 8 (black). (b) Γ as a function of time, for the same parameters and colour-coding as in (a). (c) S ˙ as a function of time. (d) S ˙ T = S ˙ + S ˙ M as a function of time.
Symmetry 14 00314 g009
Figure 10. D · S T as functions of A. The four panels correspond to T = 0.4 , 0.7 , 1 and 2, as before in Figure 4, Figure 5, Figure 6 and Figure 7. D = 10 1 , 10 2 and 10 3 are colour-coded as blue, black, red, respectively, again as before, although here the curves overlap almost everywhere.
Figure 10. D · S T as functions of A. The four panels correspond to T = 0.4 , 0.7 , 1 and 2, as before in Figure 4, Figure 5, Figure 6 and Figure 7. D = 10 1 , 10 2 and 10 3 are colour-coded as blue, black, red, respectively, again as before, although here the curves overlap almost everywhere.
Symmetry 14 00314 g010
Figure 11. J as functions of A. (a) is the square wave Equation (25). (b) is the sawtooth wave G + , (c) is the sawtooth wave G , as given by Equation (27). T = 1 , and D = 10 1 , 10 2 , 10 3 colour-coded as blue, black, red.
Figure 11. J as functions of A. (a) is the square wave Equation (25). (b) is the sawtooth wave G + , (c) is the sawtooth wave G , as given by Equation (27). T = 1 , and D = 10 1 , 10 2 , 10 3 colour-coded as blue, black, red.
Symmetry 14 00314 g011
Figure 12. The Stokes efficiency η c as functions of A. (a) is the square wave Equation (25). (b) is the sawtooth wave G + , (c) is the sawtooth wave G , as given by Equation (27). T = 1 , and D = 10 1 , 10 2 , 10 3 colour-coded as blue, black, red.
Figure 12. The Stokes efficiency η c as functions of A. (a) is the square wave Equation (25). (b) is the sawtooth wave G + , (c) is the sawtooth wave G , as given by Equation (27). T = 1 , and D = 10 1 , 10 2 , 10 3 colour-coded as blue, black, red.
Symmetry 14 00314 g012
Figure 13. L D 1 / 2 / A as functions of A. (a) is the square wave Equation (25). (b) is the sawtooth wave G + , (c) is the sawtooth wave G , as given by Equation (27). T = 1 , and D = 10 1 , 10 2 , 10 3 colour-coded as blue, black, red.
Figure 13. L D 1 / 2 / A as functions of A. (a) is the square wave Equation (25). (b) is the sawtooth wave G + , (c) is the sawtooth wave G , as given by Equation (27). T = 1 , and D = 10 1 , 10 2 , 10 3 colour-coded as blue, black, red.
Symmetry 14 00314 g013
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Kim, E.-j.; Hollerbach, R. Exact Time-Dependent Solutions and Information Geometry of a Rocking Ratchet. Symmetry 2022, 14, 314. https://doi.org/10.3390/sym14020314

AMA Style

Kim E-j, Hollerbach R. Exact Time-Dependent Solutions and Information Geometry of a Rocking Ratchet. Symmetry. 2022; 14(2):314. https://doi.org/10.3390/sym14020314

Chicago/Turabian Style

Kim, Eun-jin, and Rainer Hollerbach. 2022. "Exact Time-Dependent Solutions and Information Geometry of a Rocking Ratchet" Symmetry 14, no. 2: 314. https://doi.org/10.3390/sym14020314

APA Style

Kim, E. -j., & Hollerbach, R. (2022). Exact Time-Dependent Solutions and Information Geometry of a Rocking Ratchet. Symmetry, 14(2), 314. https://doi.org/10.3390/sym14020314

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop