Transient anomalous diffusion in periodic systems: ergodicity, symmetry breaking and velocity relaxation

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.

experiments 40,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 experimentally [42][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 recently 45 . 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 potentials 39,43,44 , granular gases in a homogeneous cooling state 46 and glass forming liquids 47 . 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 reads 45 ξ 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 k B 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 form 48 of period 2πL and a barrier height Δ U, namely It is reflection-symmetric whenever there exists a shift x 0 such that U(x 0 + x) = U(x 0 − 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-device 8 , 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 consideration 49  Scientific RepoRts | 6:30948 | DOI: 10.1038/srep30948 Here, the dimensionless potential 2 ) possesses the period 2π and the unit barrier height is Δ U = 1. Other parameters are: 0 and assumes the same statistical properties, namely ξ =t ( ) 0 and ξ ξ δ = −ˆˆˆˆt s t s ( ) ( ) ( ). The dimensionless noise intensity Q = k B T/Δ 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  x(0). 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), reading 45,50 For the diffusion process, the asymptotic time evolution of the mean square deviation (MSD) 〈 Δ x 2 (t)〉 becomes an increasing function of elapsing time and typically grows according to a power law 40,41 The exponent α specifies a type of anomalous diffusion. Normal diffusion is characterized for α = 1. The two distinct regimes of anomalous diffusion are quantified as 40,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  x t ( ) takes the form of a Fourier series over all possible higher harmonics 51 where v is the time-independent (dc) component while v nω (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 paper 45 we revealed the remarkable regime of thermal noise induced ratchet effect, (cf. Fig. 6  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 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 U 0 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 motion 48 . Roughly speaking, there are three classes of trajectories: ∼ .
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 2 2 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 breaking 33 . 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 breaking 33 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 attractors 52 . Therefore in a presence of noise the system is ergodic 45 . 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 (x 0 , v 0 ) for which in the deterministic case Q = 0 the particle moves along the trajectory ∼ − . x t t ( ) 0 4 . 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 ∼ .
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 ∼ − . x t t ( ) 0 4 ) 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 ≈ 10 7  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   Anomalous Diffusion: Impact of Potential Asymmetry. The sequence of superdiffusion-subdiffusionnormal 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 10 4 .

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 breaking 33 . This kind of behaviour can be quantified by the Deborah number De 53 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 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 − , U 0 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 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 =  x v (0) (0). 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, respectively 18 . 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 10 2 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 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. 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 regime [55][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 dynamics 34,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 distributions 30 , nor disorder 50,58,59 or many-body physics 41 .
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 lattices 14 and asymmetric SQUID devices 18,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 nanophysics 8 .

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  x(0). 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  x(0) to be equally distributed over the intervals [0,2π] and [− 2, 2], respectively. All our quantities of interest were averaged over 10 3 − 10 5 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.
Scientific RepoRts | 6:30948 | DOI: 10.1038/srep30948 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 10 3 times as compared to a common present-day CPU method 60 .