Introduction

Transport processes on the microscale can exhibit features rather different from those encountered in the macroscopic world. On such a scale they are strongly influenced by the presence of ubiquitous fluctuations and ambient noise. Rather than fighting the role of ambient fluctuations those may be put to work for selective manipulation of matter and diverse beneficial transport phenomena. Recent examples of such intriguing, noise-assisted phenomena include stochastic resonance1, enhancement of transport efficiency2,3,4, amplification of diffusion5,6, hyperdiffusion7 or the vast field of Brownian motor transport8. The latter constitutes an archetype scheme in which, even in the absence of an externally applied bias, directed motion may emerge by rectifying environmental noise via the mechanism of breaking of spatiotemporal symmetries of the system. This working principle may be seen as a key for understanding processes ranging from intracellular transport9, cancer cell metastasis10, transport of ions through nanopores11, cold atoms in optical lattices12,13,14 as well as vortices and Josephson phase in superconductors15,16,17,18,19,20. In view of the widespread applications of Brownian ratchets, directed motion controllability has become a focal point of research in nonequilibrium statistical physics which inspired a plethora of new microscale devices displaying unusual transport features21,22,23,24,25.

On the other hand, the phenomenon of diffusion is ever growing in attracting interest since Einstein’s and Smoluchowski’s groundbreaking studies26,27. This stems from the fact of the universal character and ubiquitous presence of diffusion in both regimes of classical transport and as well for quantum setups28,29. One century after those pioneering Einstein-Smoluchowski papers a new class of systems emerged that does not obey the common Gaussian law of large numbers30. Their behavior is dominated by large and rare fluctuations that are described by nonexponential decay laws, commonly referred to as Lévy statistics. Other characteristic features of these systems are ergodicity breaking, that is, non-equivalence between time averages and the corresponding ensemble averages31,32,33 as well as aging, i.e. a manifest dependence of physical observable on the time span between initialisation of the system and the start of the measurement34,35. Both ergodicity breaking and aging are in essence two sides of the same coin as they are intimately connected to non-stationarity or ultra-slow relaxation36,37. These systems exhibit various forms of diffusion anomalies38,39,40,41 which were also demonstrated in numerous experiments40,41. This kind of dynamics may not survive until the asymptotic long time regime, nonetheless, lately even its transient nature has been predicted theoretically and observed experimentally42,43,44.

Nevertheless, despite many years of intense and beneficial research in ratchet physics, a more detailed analysis of diffusion in such systems has been addressed only recently45. There, an archetypal model of the Brownian ratchet was studied revealing that the mean square deviation of the particle position from its average proceeds within three intermediate stages: initially starting as superdiffusion, followed next as subdiffusion and finally as normal diffusion in the asymptotic long time limit. These transient diffusion anomalies may last many orders longer than the characteristic time scales of the setup and their lifetime can be controlled by variation of the system temperature and its asymmetry. However, the mechanism behind this phenomenon as well as its controllability remains somewhat moot, thus being far from fully understood. The detected sequence of superdiffusion-subdiffusion-normal diffusion is not an exception and seems to be universal for some systems. A similar pattern has been observed also in other setups like Brownian particles in random potentials39,43,44, granular gases in a homogeneous cooling state46 and glass forming liquids47. While the above systems are significantly different, their diffusion behaviour is very similar. In this paper we base our study on a simple ratchet model to explain the underpinning mechanism responsible for such a sequence of diffusion anomalies. We analyze the extraordinary properties of systems exhibiting anomalous diffusion like ergodicity breaking and nonstationarity to establish an interrelation between directed ratchet transport and diffusion anomalies.

Model

We consider the generic model of a ratchet system which consists of (i) a classical inertial particle of mass M, (ii) moving in a deterministic asymmetric ratchet potential U(x), (iii) driven by an unbiased time-periodic force A cos(Ωt) of amplitude A and angular frequency Ω and (iv) subjected to thermal noise of temperature T. The corresponding Langevin equation reads45

where the dot and the prime denote differentiation with respect to time t and the Brownian particle coordinate x, respectively. The parameter Γ stands for the friction coefficient and kB is the Boltzmann constant. Thermal equilibrium fluctuations are modeled by δ-correlated, Gaussian white noise ξ(t) of zero mean and unit intensity, i.e.,

The spatially periodic potential U(x) is assumed to be in a double-sine form48 of period 2πL and a barrier height ΔU, namely

It is reflection-symmetric whenever there exists a shift x0 such that U(x0 + x) = U(x0 − x) for any x. The relative phase φ between the two harmonics serves as a control parameter of the reflection-asymmetry of this potential. If φ ≠ 0 then generally its reflection symmetry is broken which we in turn classify as as a ratchet-device8, cf. Fig. 1. The studied model describes a wealth of physical systems, as cited in the introductory section.

Figure 1
figure 1

Panel (a): The potential Eq. (3) for ΔU = 1 and L = 2π depicted in the symmetric case φ = 0 in comparison with ratchet one for two values of the potential asymmetry parameter φ = 0.5π and φ = 0.65π. Panel (b): the conservative potential force −U′(x).

As only relations between scales of length, time and energy are relevant but not their absolute values we next formulate the above equations of motion in its dimensionless form. To do so, we first introduce the characteristic dimensionless scales for the system under consideration49

so that the dimensionless form of the Langevin dynamics (1) reads

Here, the dimensionless potential possesses the period 2π and the unit barrier height is ΔU = 1. Other parameters are: m = M/(Γτ0), a = (L/ΔU)A, ω = τ0Ω. The rescaled thermal noise reads and assumes the same statistical properties, namely and . The dimensionless noise intensity Q = kBT/ΔU is the ratio of thermal and the activation energy the particle needs to overcome the nonrescaled potential barrier. From now on we will use only the dimensionless variables and shall omit the hat in all quantities appearing in the Langevin equation (5). We stress that all forces on the right hand side of Eq. (5) are non-biased: the spatial period-average of U′(x), the temporal period average of acos(ωt) and the noise average of ξ(t) are zero.

Quantifiers of diffusive transport

There exist suitable quantifiers which characterize the Brownian particle diffusion and spread of its trajectories. The most common one is the mean square deviation (or variance) of the coordinate degree of freedom x(t), namely,

where the averaging is over all possible thermal noise realizations as well as over initial conditions for the position x(0) and the velocity . The latter is mandatory due to the fact that especially in the deterministic limit of vanishing thermal fluctuations intensity Q → 0 the dynamics may not be ergodic and corresponding results may be affected by a specific choice of those initial conditions. Even though the diffusive motion of x(t) deviates at intermediate times from being normal it is nevertheless appealing to introduce a time-dependent “diffusion coefficient” D(t), reading45,50

For the diffusion process, the asymptotic time evolution of the mean square deviation (MSD) 〈Δx2(t)〉 becomes an increasing function of elapsing time and typically grows according to a power law40,41

The exponent α specifies a type of anomalous diffusion. Normal diffusion is characterized for α = 1. The two distinct regimes of anomalous diffusion are quantified as40,41: subdiffusion if 0 < α < 1 and superdiffusion if α > 1. In the former, subdiffusive case, the MSD increases over time slower than normal while for superdiffusion it is growing faster than normal diffusion. The time dependent diffusion coefficient D(t) therefore allows to differentiate between these anomalous intermediate regimes of diffusion; namely superdiffusion occurs when D(t) increases, the case of decreasing D(t) corresponds to subdiffusion and for D(t) = const. normal diffusion is developed. Only when asymptotically α approaches unity the time-independent diffusion coefficient D is given by

To gain insight into the origin of the diverse diffusion phenomena in the system under consideration it helps to study the occurrence of non-vanishing directed transport, that is the noise and time averaged velocity v(t), defined by the relation

where T = 2π/ω is a period of the time-periodic force. In the asymptotic long time limit this quantity becomes time-independent, while the noise-averaged variant alone assumes a time-periodic function of the asymptotic phase-space probability. Put differently, in such a limit, the mean velocity takes the form of a Fourier series over all possible higher harmonics51

where v is the time-independent (dc) component while vnω(t) denote time-periodic higher harmonic functions of zero average over the fundamental period T = 2π/ω of the driving. The dc component v is obtained as

Due to the presence of the external driving the Brownian particle is taken far away from thermal equilibrium and a time-dependent nonequilibrium state is reached in the asymptotic long time regime. Since all forces in the right hand side of Eq. (5) are non-biased, a necessary condition for the occurrence of directed transport v ≠ 0 is the breaking of the reflection symmetry of the potential U(x)8.

Anomalous Diffusion: The Role of Finite Temperature

Even with the help of our innovative computational simulations, see the section Methods, the system described by Eq. (5) is too complex to analyze numerically in a systematic manner. It possess a five-dimensional parameter space {m, a, ω, φ, Q}. However, in our recent paper45 we revealed the remarkable regime of thermal noise induced ratchet effect, (cf. Fig. 6 in ref. 3): At zero temperature Q = 0 the averaged velocity v ≈ 0. If Q starts to increase, the velocity v also increases attaining the maximal value v ≈ 0.4 for Q ≈ 0.0004 and next it decreases to zero as Q grows. It is so for φ = π/2 in the potential (3) implying its most asymmetrical ratchet form. In this regime, evolution of the mean square deviation can be divided into three time domains as depicted in Fig. 2(a) (see also Fig. 3 in ref. 45): the early period of superdiffusion τ1 (where D(t) is an increasing function of time), the intermediate interval τ2 where subdiffusion is developed (there D(t) is a decreasing function of time) and approaches the asymptotic long time regime where normal diffusion occurs (i.e. D(t) ≈ constant). In Fig. 2(c) we show how the crossover time τ1 depends on temperature. The lifetime τ1 of the superdiffusion regime is extremely long in the low temperature limit, c.f. the curve corresponding to Q = 0.0001 in Fig. 2(a). However, persistent superdiffusion occurs only in the deterministic case when formally Q = 0. The deflection of this stage of diffusion can expressively be noted as temperature increases. For sufficiently high temperature the motion is initially superdiffusive and next normal diffusion occurs, see the case of Q = 0.005 in Fig. 2(a).

Figure 2
figure 2

Control of anomalous diffusion regimes by temperature Q ∝ T in the ratchet potential Eq. (3) with φ = π/2.

Panel (a): the diffusion coefficient D(t). Panel (b): relaxation of the velocity v(t) to its nonequilibrium stationary state. Panel (c): the crossover time τ1 separating superdiffusion and subdiffusion stages of the diffusion process and the velocity relaxation time τv versus system temperature. Panel (d): basins of attraction for the asymptotic long time particle velocity v. Red colour indicates the running states with the positive velocity v = 0.4, blue corresponds to its negative counterpart v = −0.4 and green marks the locked states v ≈ 0. Parameters are: m = 6, a = 1.899, ω = 0.403 and Q = 0.0004.

Strong ergodicity breaking

To explain the above diffusion anomalies let us now study the deterministic dynamics Q = 0 and the corresponding structure of basins of attraction for the asymptotic long time velocity v. The result is shown in Fig. 2(d). There exist only three attractors: the set U+ which consists of all running states with positive velocity v ≈ 0.4 (marked by red colour), the set U− of states running with negative velocity v ≈ −0.4 (marked by blue colour) and the set U0 of locked states v ≈ 0 when the motion is bounded to a finite number of the potential wells (green colour). This simple structure is crucial for the occurrence of initial superdiffusive (ballistic) stage of motion48. Roughly speaking, there are three classes of trajectories: , and . Moreover, the running states of either positive or negative velocity occupy almost even area of the initial phase space of the system. This fact causes that the contribution of the averaged trajectory to the mean-square deviation 〈Δx2(t)〉 = 〈x2(t)〉 − 〈x(t)〉2 is negligible small 〈x(t)〉 ≈ 0 in comparison to its second moment . As a consequence also for the asymptotic long time limit t → ∞. Therefore then the superdiffusive regime is persistent and its lifetime τ1 → ∞ when Q → 0. In fact, it is an example of ergodicity breaking33. The phase space of the system consists of three (v ± 0.4, v ≈ 0) disjoint sets which are mutually inaccessible. In other words a phase point starting e.g. in the blue region in Fig. 2(d) and evolving according to Eq. (5) will stay in this region forever and cannot penetrate the other two sets. As a consequence, statistical properties of the system deduced from time average along one sufficiently long trajectory are not the same as determined by the ensemble averaging. Such a decomposability of the phase space is called strong ergodicity breaking33 in the language accepted in the community of nonequilibrium statistical physicists. This is essentially a refinement of original Boltzmann’s definition and in other words means that a trajectory associated with time evolution of an initial condition cannot pass arbitrarily close to any point of a phase-space, put simply initial conditions are never fully forgotten - a property which finds its analogue in any definition of ergodicity.

Superdiffusion versus thermal fluctuations

We expect that thermal noise with the intensity Q ≠ 0 perturbs deterministic dynamics and equilibrium fluctuations enable stochastic escape events connecting coexisting deterministic disjoint attractors52. Therefore in a presence of noise the system is ergodic45. In Fig. 3 we depict 1024 sample Brownian particle trajectories for different thermal fluctuations intensity Q ∝ T. Let us consider for definiteness an initial condition (x0, v0) for which in the deterministic case Q = 0 the particle moves along the trajectory . When Q > 0, from time to time the particle is kicked by thermal fluctuations and escapes from its initial trajectory jumping on the other like or . After some time, the particle can again break away the ongoing trajectory and jump on the different one. This process reduces a number of trajectories which contribute to superdiffusion (here ) and destroys the initial simple structure of only three attractors. It can be directly seen in Fig. 4 where we present the probability distributions P(v) of the asymptotic long time velocity v for different values of temperature Q ∝ T. In the deterministic limit of vanishing thermal noise intensity Q = 0 there are two almost equally pronounced δ-peaks representing the running solutions of either positive v = 0.4 or negative v = −0.4 direction. The deviation of this structure of states is distinctly noticeable as temperature increases. Thermal noise first blurs δ-peaks observed in the deterministic case and then eliminates both the negative running as well as locked states in favor of the positive ones v ≈ 0.4. Further increase in temperature causes additional smearing of the probability distribution P(v). Clearly, in presence of thermal noise the system forgets about its initial conditions so that ergodicity is restored.

Figure 3
figure 3

Various typical sample trajectories of the Brownian particle dynamics.

These are depicted for different temperature of the system Q ∝ T. Panel (a–c) corresponds to Q = 0.0001, Q = 0.0004 and Q = 0.005, respectively. Other parameters are the same as in Fig. 2.

Figure 4
figure 4

The probability distributions (histograms) P(v) of the individual asymptotic long time particle velocity v are presented for different values of temperature Q ∝ T.

Panel (a): Q = 0.0001, (b): Q = 0.0004 and (c): Q = 0.005. Other parameters are the same as in Fig. 2.

Overall, if temperature grows the jumps between different types of solutions are more and more frequent (as in Fig. 3(c)) and the mean time to destroy the deterministic structure of attractors becomes shorter. As a consequence, the lifetime τ1 of superdiffusion is also shortened. It is now obvious that τ1 is a decreasing function of temperature Q ∝ T in accordance with results presented in Fig. 2(a,c).

Superdiffusion versus velocity relaxation

Now, in order to quantify the mean time to destroy the deterministic structure of attractors we study a relaxation process of the velocity degree of freedom. The result is presented in Fig. 2(b) where time evolution of the velocity v(t) defined in Eq. (10) is depicted. In the limiting regime of weak thermal noise the velocity quickly relaxes to zero v(t) ≈ 0, see e.g. the case Q = 7·10−5. It means that at least during the data acquisition time in the numerical experiment equilibrium fluctuations are not able to wipe out the deterministic structure of basins of attraction so that the system still remembers three solutions and the mean velocity is zero. However, since in a presence of thermal noise the attractors coexisting in the phase space are connected via stochastic escape events we expect that even if temperature is low a relevant trajectory will eventually fully sample the available state space of the system. For moderate and strong thermal fluctuations the system forgets about its initial conditions since the velocity relaxes to the nonzero steady state v = 0, c.f. the case Q = 0.005 in Fig. 2(b). Its time significantly influences the duration of superdiffusion τ1 as we demonstrate in Fig. 2(c). A rough synchronization of these two curves is observed there. It is remarkable that for the considered case of φ = 0.5π these two characteristic times increase over three orders of magnitude when the intensity Q of thermal fluctuations varies in the small interval [0.0002, 0.01]. The bottom bound of this window is chosen so that thermal noise is still able to destroy the deterministic structure of coexisting regular attractors in the data acquisition time .

Anomalous Diffusion: Impact of Potential Asymmetry

The sequence of superdiffusion-subdiffusion-normal diffusion is detected only when reflection symmetry of the potential U(x) is broken, i.e. when it has a ratchet form. We already demonstrated that the initial superdiffusive stage of motion stems from the deterministic structure of the coexisting counter-propagating regular directed transporting attractors, see also the discussion on the power exponent α in ref. 45. Figure 5(a) presents the impact of the potential asymmetry parameter φ on the diffusion coefficient D(t) in the vicinity of φ = π/2 where the crossover times τ1 and τ2 are extremely sensitive to its changes. Panel (b) of the same figure illustrates time evolution of the velocity v(t) for selected values of φ. Again, astonishingly, this quantity relaxes to its nonequilibrium steady state with the relaxation time τv strongly dependent on the magnitude of the potential asymmetry φ. Similarly as in the case of temperature influence, one can detect an interrelation between the superdiffusion crossover time τ1 and the velocity relaxation time τv so that if the potential asymmetry parameter φ increases both τ1 and τv grow in the synchronized way. This feature is better visualized in the panel (c) where we contrast the relaxation time τv together with the superdiffusion lifetime τ1 for the fixed temperature Q = 0.0004. Both of them vary in the similar way and are extraordinarily sensitive to alteration of the potential asymmetry φ: Small changes of the order 10−2 in φ are accompanied by the giant increase of τ1 and τv of the order 104.

Figure 5
figure 5

Impact of the asymmetry parameter φ of the potential U(x) on the diffusion process.

Panel (a): sensitivity of the diffusion coefficient D(t) in the vicinity of φ = π/2. Panel (b): relaxation of the velocity v(t) to its nonequilibrium steady state. Panel (c): the crossover time τ1 separating superdiffusion and subdiffusion stages of the diffusion process and the velocity relaxation time τv versus potential asymmetry parameter φ. Panel (d): the crossover time τ1 versus φ depicted for different values of temperature Q ∝ T. In the inset τ1 is shown for Q = 0.002 in larger interval of φ capturing its maximum. Other parameters are the same as in Fig. 2.

Weak ergodicity breaking

The time evolution of the diffusion coefficient D(t) illustrated in Figs 2(a) and 5(a) is very similar. However, if temperature tends to zero Q ∝ T → 0 the superdiffusion crossover time monotonically increases to infinity τ1 → ∞ and ergodicity is broken. On the other hand, when the asymmetry parameter φ increases the crossover time τ1 monotonically grows reaching the maximal value and next it starts to decrease, see the inset in Fig. 5(d). In other words when Q ≠ 0 the lifetime τ1 does not tend to infinity for any value of the potential asymmetry parameter φ. It means that ergodicity cannot be strictly broken by change of φ although τ1 may be extremely long. In such a case the whole phase space is still accessible because of the presence of thermal fluctuations, however, the time after it is fully sampled might be enormously long. From an experimental point of view the problem of ergodicity is a matter of time scale. In practice, there is not much difference whether the system is nonergodic or ergodic but exhibiting an extremely slow relaxation. If this effect is pronounced one often speaks about weak ergodicity breaking33. This kind of behaviour can be quantified by the Deborah number De53

which is a ratio of a relaxation time τ of a given observable and the time of observation . If it diverges De → ∞ then the system gives the impression of being “frozen” and the question whether it may eventually sample the full state space becomes irrelevant. This can happen not only because is short but also because τ is extremely long. In our case the system behaves as weakly nonergodic when the superdiffusion lifetime τ1 or the velocity relaxation time τv is sufficiently large so that the condition is satisfied.

Superdiffusion versus potential asymmetry

To gain insight into the origin of the extremal sensitivity of the superdiffusion lifetime τ1 to changes of the potential asymmetry parameter φ in the vicinity of φ = π/2 we now study variation of the structure of basins of attraction with respect to φ. In Fig. 6 we present them in the deterministic limit Q = 0 and for different values of the parameter φ. Again, there are three classes of trajectories v = 0.4, v = −0.4 and v ≈ 0 with three corresponding basins of attraction U+, U−, U0 marked by red, blue and green colour, respectively. For φ = 0.488π the set U+ is much larger than U− and in consequence the superdiffusion lifetime τ1 is short. Surprisingly, we note that a tiny increase of the order 10−2 in the asymmetry parameter φ is accompanied by a rapid growth of the set U− of the initials conditions for which the asymptotic long time velocity v is negative v = −0.4. Then the size of the sets U+ and U− becomes comparable, there are more trajectories with the negative velocity and the spread of trajectories is large. Additionally, the noisy system needs more time to destroy this evident structure of basins of attraction via the mechanism of kicking out the particle from trajectories with the negative velocity v ≈ −0.4 onto those with positive ones v ≈ 0.4. Certainly, it is one of the reasons why the superdiffusion lifetime τ1 increases. However, we observe in the inset of Fig. 5(d) that the superdiffusion crossover time τ1 can still grow although the structure of basins of attraction for the asymptotic velocity v is not changing significantly any more, c.f. Figs 2(d) and 6(d). We note that for φ = π/2, asymmetry of the potential is maximal, the depth of the potential wells is maximal, mean first passage time T− to overcome the potential barrier to the left is much longer that mean first passage times T+ to overcome the potential barrier to the right (not depicted). Further increase of φ makes the potential less asymmetric, its depth decreases, both T− and T+ decrease and become comparable for φ ≈ 0.6. Therefore its easier for thermal fluctuations to kick the particle onto the trajectory running in negative direction with the asymptotic velocity v < 0. We support this point of view by resource to simulation of sample trajectories for different values of the asymmetry parameter φ. The result is shown in Fig. 7. Indeed, we observe there that as φ increases (i) progressively more particles move with the negative velocity and (ii) they start doing it earlier. These features are also visualized in Fig. 8 where we present the probability distributions P(v) of the asymptotic long time particle velocity for different values of the potential asymmetry parameter φ. For fixed temperature Q ≠ 0 an increase in φ causes elimination of the positive running states in favour of the negative ones v ≈ −0.4. In turn, this mechanism is reversed with respect to that which occurred in the case of temperature changes, cf. Fig. 4.

Figure 6
figure 6

Basins of attraction for the asymptotic long time particle velocity v are presented for different values of the potential asymmetry parameter φ.

Panel (a): φ = 0.488π, (b): φ = 0.490π, φ = 0.494π and (d): φ = 0.51π. Red color indicates the running states with the positive velocity v = 0.4, blue corresponds to its negative counterpart v = −0.4 and green color marks the locked states v ≈ 0. Other parameters are the same as in Fig. 2, except thermal noise intensity which is fixed to zero Q = 0.

Figure 7
figure 7

Characteristic sample trajectories of the Brownian particle depicted for different potential asymmetry parameters φ.

Panel (a–c) corresponds to φ = 0.5π, φ = 0.505π and φ = 0.51π, respectively. Other parameters are the same as in Fig. 2 and Q = 0.0004.

Figure 8
figure 8

Probability distributions (histograms) P(v) of the individual asymptotic long time particle velocity v are presented for different values of the potential asymmetry parameter φ.

Panel (a): φ = 0.5π, (b): φ = 0.505π and (c): φ = 0.51π. Other parameters are the same as in Fig. 2 and Q = 0.0004.

Role of Initial Velocity Distribution

Last but not least, we analyze impact of an initial probability distribution for the Brownian particle velocity P(v) on the observed diffusion anomalies. Transient phenomena often depend on a state in which a system is initially prepared therefore it is a vital question to pose. In our case this state is described by probability distributions for the particle coordinate x(0) and velocity . We focus only on the latter since the former can hardly be manipulated in many experimentally accessible physical systems that can be described by the studied model. For example, in an asymmetric SQUID device composed of three resistively and capacitively Josephson junctions the Brownian particle coordinate x(t) and velocity v(t) translates to the Josephson phase Ψ(t) and voltage drop V(t) across the setup, respectively18. The first quantity cannot be directly manipulated whereas the second may be easily altered by applying an external electric field. For this reason we now compare the diffusion processes occurring in the model Eq. (5) starting from three different types of an initial condition for the velocity v(0): a Dirac delta P(v) = δ(v), a random variable distributed uniformly P(v) = U(−2, 2) and allocated normally with a zero mean and a standard deviation equals to one P(v) = N(0, 1). The result is shown in Fig. 9. In panel (a) we study time evolution of the diffusion coefficient D(t) for selected temperature Q = 0.001. The discrepancies between these distributions die out after approximately 102 characteristic units of time. The system relaxes to normal diffusion in a universal pattern which is independent on the form of the initial velocity distribution. It can be divided into three time domains: the early period of superdiffusion, the intermediate interval where subdiffusion is developed and approaches the asymptotic long time regime where normal diffusion occurs. Therefore at least for the studied regime the observed behaviour is omnipresent. It is confirmed in the remaining panels of this figure when we illustrate a duration of the period of superdiffusion respectively as a function of temperature Q ∝ T and the potential asymmetry parameter φ. In particular, the discussed extremal sensitivity to changes of the latter is robust with respect of variation of the initial velocity distribution P(v).

Figure 9
figure 9

Impact of the initial probability distribution P(v) for the Brownian particle velocity on the diffusion anomalies observed in the system.

Panel (a): the time dependent diffusion coefficient D(t) for temperature Q = 0.001. Panel (b,c): the crossover time τ1 separating superdiffusion and subdiffusion motion depicted versus temperature Q ∝ T and the potential asymmetry parameter φ, respectively. For the other parameters read Fig. 2.

Discussion

With this work we investigated diffusion anomalies occurring in the archetype model of an inertial Brownian ratchet. In particular, we established the clear connection between relaxation of the average velocity v(t) and the transient anomalous diffusion quantified by the “diffusion coefficient” D(t). The key ingredients for the emergence of the observed sequence superdiffusion-subdiffusion-normal diffusion are following: (i) the system is driven far away from thermal equilibrium into a time dependent nonequilibrium state, (ii) reflection symmetry of the potential is violated so that the directed transport may arise, i.e. v ≠ 0, (iii) ergodicity of the corresponding deterministic system is strongly broken due to the coexistence of two counter-propagating regular attractors which almost evenly occupy the accessible phase space. As we demonstrated upon an introduction of thermal noise the ratchet system will be typically weakly nonergodic with an extremely slow relaxation of the degrees of freedom leading to a whole range of anomalous diffusion phenomena including the superdiffusive and subdiffusive motion.

We want to stress that the mechanism of two counter-propagating attractors which leads to transient superdiffusion can also be detected in symmetric systems where the directed transport is strictly ruled out by symmetries, i.e. v = 0. Moreover, transient subdiffusion may also be found in these systems as well, see Addendum in ref. 54. After a sufficiently long time such anomalies die out and the diffusion is always normal. On the other hand, the sequence superdiffusion-subdiffusion-normal diffusion is detected in the presented parameter regime only for the ratchet system.

The discussed diffusional anomalies with slow relaxation may give the wrong impression that the anomalous diffusion is observed even in the asymptotic long time limit. The latter can happen only in systems whose dynamics is dominated by large and rare fluctuations that are characterized by broad distributions with power-law tails. A hallmark of these non-Gaussian distributions is the divergence of their second and/or first moment. These systems do not obey the law of large numbers and the strong convergence to the Gaussian according to the central limit theorem is broken. Consequently, diffusion anomalies may be observed even in the asymptotic long time regime55,56,57. Despite the fact that diffusion anomalies observed in the studied model have only transient nature their duration can be fine-tuned by adjustment of the setup parameters so that they last many orders longer than characteristic time scales of the system. Therefore from an experimental point of view they may be safely treated as persistent effects.

Conclusions

In summary, we presented the detailed qualitative theory illuminating the coexistence of the directed transport as well as the anomalous diffusion processes in a generic system of an inertial Brownian ratchet. We explained the underlying physical mechanism standing behind the emergence of diffusion anomalies and the parameter-dependent control of their extended regimes. Extremely slow relaxation and ergodicity breaking typically describe glassy dynamics34,35. Our work shows that these features can be observed even within a more straightforward, one dimensional classical Markovian dynamics with Brownian motion of its inherent long-time Gaussian nature, i.e. without the need to introduce heavy-tailed distributions30, nor disorder50,58,59 or many-body physics41.

The appealing strength and beauty of Brownian motion with its intrinsic Gaussian noise propagator lies in its universality and therefore our findings can be straightforwardly corroborated experimentally with a wealth of physical systems outlined in the introductory part of the article. One of the most promising setups for this purpose are optical lattices14 and asymmetric SQUID devices18,20. It is because of their high tunability: the period, the amplitude and the symmetry of the corresponding potential may be modified in a controlled way.

In view of the widespread applications of Brownian motor setups and ratchet devices our research may carry potential impact for further development of a working principle of a nanomotor operating on smallest scales, such as occurring in diverse areas of nanophysics8.

Methods

The Fokker-Planck equation corresponding to the analyzed Langevin Eq. (5) understandably cannot be solved by use of analytical means. Therefore, in order to obtain the relevant characteristics we have to carry comprehensive numerical simulations. We did so by employing a weak version of the stochastic second-order predictor-corrector algorithm with a time-step typically set to about (10−3 − 10−2) × T. Because Eq. (5) is a second-order differential equation, we need to specify two initial conditions, namely x(0) and . For some regimes the system dynamics does exhibit nonergodic behavior; thus, in order to avoid a dependence of the presented results on a specific choice of the initial conditions, unless stated otherwise, we have chosen x(0) and to be equally distributed over the intervals [0,2π] and [−2, 2], respectively. All our quantities of interest were averaged over 103 − 105 sample trajectories. All numerical calculations have been performed by use of a CUDA environment as implemented on a modern desktop GPU. This procedure did allow for a speedup of a factor of the order 103 times as compared to a common present-day CPU method60.

Additional Information

How to cite this article: Spiechowicz, J. et al. Transient anomalous diffusion in periodic systems: ergodicity, symmetry breaking and velocity relaxation. Sci. Rep. 6, 30948; doi: 10.1038/srep30948 (2016).