Abstract
We study far from equilibrium transport of a periodically driven inertial Brownian particle moving in a periodic potential. As detected for a SQUID ratchet dynamics, the mean square deviation of the particle position from its average may involve three distinct intermediate, although extended diffusive regimes: initially as superdiffusion, followed by subdiffusion and finally, normal diffusion in the asymptotic long time limit. Even though these anomalies are transient effects, their lifetime can be many, many orders of magnitude longer than the characteristic time scale of the setup and turns out to be extraordinarily sensitive to the system parameters like temperature or the potential asymmetry. In the paper we reveal mechanisms of diffusion anomalies related to ergodicity of the system, symmetry breaking of the periodic potential and ultraslow relaxation of the particle velocity towards its steady state. Similar sequences of the diffusive behaviours could be detected in various systems including, among others, colloidal particles in random potentials, glass forming liquids and granular gases.
Similar content being viewed by others
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.
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).
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.
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.
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.
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).
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).
References
Gammaitoni, L., Hänggi, P., Jung, P. & Marchesoni, F. Stochastic resonance. Rev. Mod. Phys. 70, 223 (1998).
Spiechowicz, J., Hänggi, P. & Åuczka, J. Brownian motors in the microscale domain: Enhancement of efficiency by noise. Phys. Rev. E 90, 032104 (2014).
Spiechowicz, J. & Åuczka, J. Efficiency of the SQUID ratchet driven by external current. New J. Phys. 17, 023054 (2015).
Spiechowicz, J., Åuczka, J. & Machura, Å. Efficiency of transport in periodic potentials: dichotomous noise contra deterministic force, J. Stat. Mech. 054038 (2016).
Reimann, P., Van den Broeck, C., Linke, H., Hänggi, P., Rubi, J. M. & Perez-Madrid, A. Giant Acceleration of Free Diffusion by Use of Tilted Periodic Potentials. Phys. Rev. Lett. 87, 010602 (2001).
Spiechowicz, J. & Åuczka, J. Josephson phase diffusion in the superconducting quantum interference device ratchet. Chaos 25, 053110 (2015).
Siegle, P., Goychuk, I. & Hänggi, P. Origin of Hyperdiffusion in Generalized Brownian Motion. Phys. Rev. Lett. 105, 100602 (2010).
Hänggi, P. & Marchesoni, F. Artificial Brownian motors: Controlling transport on the nanoscale. Rev. Mod. Phys. 81, 387 (2009).
Bressloff, B. C. & Newby, J. M. Stochastic models of intracellular transport. Rev. Mod. Phys. 85, 135 (2013).
Mahmud, G. et al. Directing cell motions on micropatterned ratchets. Nat. Phys. 5, 606 (2009).
Karnik, R. et al. Rectification of Ionic Current in a Nanofluidic Diode. Nano Lett. 7, 547 (2007).
Schiavoni, M., Sanchez-Palencia, L., Renzoni, F. & Grynberg, G. Phase control of directed diffusion in a symmetric optical lattice. Phys. Rev. Lett. 90, 094101 (2003).
Gommers, R., Bergamini, S. & Renzoni, F. Dissipation-induced symmetry breaking in a driven optical lattice. 95, 073003 (2005).
Denisov, S., Flach, S. & Hänggi, P. Tunable transport with broken space-time symmetries. Phys. Rep. 538, 77 (2014).
Villegas, J. E. et al. A Superconducting Reversible Rectifier That Controls the Motion of Magnetic Flux Quanta. Science 302, 1188 (2003).
de Souza Silva, C. C., Van de Vondel, J., Morelle, M. & Moshchalkov, V. V. Controlled multiple reversals of a ratchet effect. Nature 440, 651 (2006).
Zapata, I., Bartussek, R., Sols, F. & Hänggi, P. Voltage Rectification by a SQUID Ratchet. Phys. Rev. Lett. 77, 2992 (1996).
Spiechowicz, J., Hänggi, P. & Åuczka, J. Josephson junction ratchet: The impact of finite capacitances. Phys. Rev. B 90, 054520 (2014).
Weiss, S. et al. Ratchet effect in dc SQUIDs. Europhys. Lett. 51, 499 (2000).
Sterk, A., Kleiner, R. & Koelle, D. Three-Junction SQUID Rocking Ratchet. Phys. Rev. Lett. 95, 177006 (2005).
Hernández, J. V., Kay, E. R. & Leigh, D. A. A reversible synthetic rotary molecular motor. Science 306, 1532 (2004).
Costache, M. V. & Valenzuela, S. O. Experimental spin ratchet. Science 330, 1645 (2010).
Serreli, V., Lee, C. F., Kay, E. R. & Leigh, D. A. A molecular information ratchet. Nature 445, 523 (2007).
Drexler, C. et al. Magnetic quantum ratchet effect in graphene. Nat. Nanotechnol. 8, 104 (2013).
Grossert, C., Leder, M., Denisov, S., Hänggi, P. & Weitz, M. Experimental control of transport resonances in a coherent quantum rocking ratchet. Nature Communications 7, 10440 (2016).
Einstein, A. Ãber die von der molekularkinetischen Theorie geforderte Bewegung von in ruhenden Flssigkeiten suspendierten Teilchen. Annalen der Physik 17, 4549 (1905).
Smoluchowski, M. Zur kinetischen Theorie der Brownschen Molekularbewegung und der Suspensionen. Annalen der Physik 21, 756 (1906).
Hänggi, P. & Marchesoni, F. 100 Years of Brownian Motion. Chaos 15, 026101 (2005).
Sokolov, I. M. & Klafter, J. From diffusion to anomalous diffusion: A century after Einsteinâs Brownian motion. Chaos 15, 026103 (2005).
Lutz, E. & Renzoni, F. Beyond Boltzmann-Gibbs statistical mechanics in optical lattices. Nat. Phys. 9, 615 (2013).
Jeon, J. H. et al. In vivo anomalous diffusion and weak ergodicity breaking of lipid granules. Phys. Rev. Lett. 106, 048103 (2011).
Burov, S., Jeon, J. H., Metzler, R. & Barkai, E. Single particle tracking in systems showing anomalous diffusion: the role of weak ergodicity breaking. Phys. Chem. Chem. Phys. 13, 1800 (2011).
Meroz, Y. & Sokolov, I. M. A toolbox for determining subdiffusive mechanisms. Phys. Rep. 573, 1 (2015).
Bouchaud, J. P. Weak ergodicity breaking and aging in disordered systems. J. Phys. I 2, 1705 (1992).
Burov, S., Metzler, R. & Barkai, E. Aging and nonergodicity beyond the Khinchin theorem. PNAS 107 13228 (2010).
Sokolov, I. M. Models of anomalous diffusion in crowded environments. Soft Matter. 8, 9043 (2012).
Barkai, E., Garini, Y. & Metzler, R. Strange Kinetics of single molecules in living cells. Phys. Today 65, 29 (2012).
Metzler, R. & Klafter, J. The random walkâs guide to anomalous diffusion: a fractional dynamics approach. Phys. Rep. 339, 1 (2000).
Goychuk, I. & Kharchenko, V. O. Anomalous Features of Diffusion in Corrugated Potentials with Spatial Correlations: Faster than Normal and Other Surprises. Phys. Rev. Lett. 113, 100601 (2014).
Metzler, R., Jeon, J. H., Cherstvy, A. G. & Barkai, E. Anomalous diffusion models and their properties: non-stationarity, non-ergodicity and ageing at the centenary of single particle tracking. Phys. Chem. Chem. Phys. 16, 24128 (2014).
Zaburdaev V., Denisov S. & Klafter J. Lévy walks. Rev. Mod. Phys. 87, 483 (2015).
Bronstein I. et al. Transient Anomalous Diffusion of Telomeres in the Nucleus of Mammalian Cells. Phys. Rev. Lett. 103, 018102 (2009).
Hanes, R. D. L. & Egelhaaf, S. U. Dynamics of individual colloidal particles in one-dimensional random potentials: a simulation study. J. Phys. Condens. Matter. 24, 464116 (2012).
Hanes, R. D. L. et al. Colloids in one dimensional random energy landscapes. Soft Matter. 8, 2714 (2012).
Spiechowicz, J. & Åuczka, J. Diffusion anomalies in ac-driven Brownian ratchets. Phys. Rev. E 91, 062104 (2015).
Bodrova, A., Dubey, A. K., Puri, S. & Brilliantov, N. Intermediate Regimes in Granular Brownian Motion: Superdiffusion and Subdiffusion. Phys. Rev. Lett. 109, 178001 (2012).
Haxton, T. K. & Liu, A. J. Kinetic heterogeneities at dynamical crossovers. EPL 90, 66004 (2010).
Jung, P., Kissner, J. G. & Hänggi, P. Regular and Chaotic Transport in Asymmetric Periodic Potentials: Inertia Ratchets. Phys. Rev. Lett. 76, 3436 (1996).
Machura, Å., Kostur, M. & Åuczka, J. Transport characteristics of molecular motors. Biosystems 94, 253 (2008).
Khoury, M., Lacasta, A. M., Sancho, J. M. & Lindenberg, K. Weak Disorder: Anomalous Transport and Diffusion Are Normal Yet Again. Phys. Rev. Lett. 106, 090602 (2011).
Jung, P. Periodically driven stochastic systems. Phys. Rep. 234, 175 (1993).
Hänggi, P., Talkner, P. & Borkovec, M. Reaction Rate Theory: Fifty Years After Kramers. Rev. Mod. Phys. 62, 251 (1990).
Reiner, M. The Deborah number. Phys. Today 17, 62 (1964).
Spiechowicz, J. & Åuczka, J. Addendum to âDiffusion anomalies in ac-driven Brownian ratchetsâ. arXiv:1506.00105v2 (2016).
Kessler, D. A. & Barkai, E. Infinite Covariant Density for Diffusion in Logarithmic Potentials and Optical Lattices. Phys. Rev. Lett. 105, 120602 (2010).
Zaburdaev, V., Denisov, S. & Hänggi, P. Space-Time Velocity Correlation Function for Random Walks. Phys. Rev. Lett. 110, 170604 (2013).
Rebenshtok, A., Denisov, S., Hänggi, P. & Barkai, E. Non-normalizable densities in strong anomalous diffusion: Beyond the central limit theorem. Phys. Rev. Lett. 112, 110601 (2014).
Sancho, J. M. et al. Diffusion on a Solid Surface: Anomalous is Normal. Phys. Rev. Lett. 92, 250601 (2004).
Sune, M., Sancho, J. M. & Lindenberg, K. Transport and diffusion of underdamped Brownian particles in random potentials. Eur. Phys. J. B 87, 201 (2014).
Spiechowicz, J., Kostur, M. & Machura, Å. GPU accelerated Monte Carlo simulation of Brownian motors dynamics with CUDA. Comp. Phys. Commun. 191, 140 (2015).
Acknowledgements
This work was supported in part by the MNiSW program via a Diamond Grant (J.S.).
Author information
Authors and Affiliations
Contributions
J.S. performed all numerical simulations in this project. All authors contributed extensively to the planning, discussion of the results and writing up of this work.
Ethics declarations
Competing interests
The authors declare no competing financial interests.
Rights and permissions
This work is licensed under a Creative Commons Attribution 4.0 International License. The images or other third party material in this article are included in the articleâs Creative Commons license, unless indicated otherwise in the credit line; if the material is not included under the Creative Commons license, users will need to obtain permission from the license holder to reproduce the material. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/
About this article
Cite this article
Spiechowicz, J., Åuczka, J. & Hänggi, P. Transient anomalous diffusion in periodic systems: ergodicity, symmetry breaking and velocity relaxation. Sci Rep 6, 30948 (2016). https://doi.org/10.1038/srep30948
Received:
Accepted:
Published:
DOI: https://doi.org/10.1038/srep30948
This article is cited by
-
Collective transient ratchet transport induced by many elastically interacting particles
Scientific Reports (2021)
-
Tunable particle separation via deterministic absolute negative mobility
Scientific Reports (2020)
-
Noise-enabled optical ratchets
Scientific Reports (2017)
-
Subdiffusion via dynamical localization induced by thermal equilibrium fluctuations
Scientific Reports (2017)