# Negative mobility of a Brownian particle:

strong damping regime

###### Abstract

We study impact of inertia on directed transport of a Brownian particle under non-equilibrium conditions: the particle moves in a one-dimensional periodic and symmetric potential, is driven by both an unbiased time-periodic force and a constant force, and is coupled to a thermostat of temperature . Within selected parameter regimes this system exhibits negative mobility, which means that the particle moves in the direction opposite to the direction of the constant force. It is known that in such a setup the inertial term is essential for the emergence of negative mobility and it cannot be detected in the limiting case of overdamped dynamics. We analyse inertial effects and show that negative mobility can be observed even in the strong damping regime. We determine the optimal dimensionless mass for the presence of negative mobility and reveal three mechanisms standing behind this anomaly: deterministic chaotic, thermal noise induced and deterministic non-chaotic. The last origin has never been reported. It may provide guidance to the possibility of observation of negative mobility for strongly damped dynamics which is of fundamental importance from the point of view of biological systems, all of which in situ operate in fluctuating environments.

###### keywords:

Brownian motion, periodic symmetric systems, negative mobility, strong damping regime## 1 Introduction

When a system at thermal equilibrium is exposed to a weak external static force, its response is in the same direction as this of applied bias towards a new equilibrium. This restriction is no longer valid under nonequilibrium conditions when already an unperturbed system may exhibit a current due to the ratchet effect denisov2014 (). Another example is the seemingly paradoxical situation of the negative mobility phenomenon when the system response is opposite to the applied constant force hanggi2009 (). Such anomalous transport behaviour was predicted theoretically in 2007 in a system consisting of an inertial Brownian particle moving in a one-dimensional periodic symmetric potential machura2007 (). Within a year of this discovery, negative mobility was confirmed experimentally in the experiment involving determination of current-voltage characteristics of the microwaved-driven Josephson junction nagel2008 (). Yet further examples of this phenomenon has been described theoretically in companionship of coloured noise kostur2009 (), white Poissonian noise spiechowicz2013jstatmech (); spiechowicz2014pre (), dichotomous process spiechowicz2016jstatmech () and for Brownian motion with presence of time-delayed feedback hennig2009 (); du2011 (), non-uniform space-dependent damping mulhern2013 () and potential phase modulation dandogbessi2015 (). Other illustrations include a vibrational motor du2012 (), two coupled resistively shunted Josephson junctions januszewski2011 (); machura2012 (), active Janus particles in a corrugated channel ghosh2014 (), entropic electrokinetics malgaretti2014 () as well as nonlinear response of inertial tracers in steady laminar flows sarracino2016 ().

Modelling systems and understanding their generic properties discloses which components of the setup are crucial and which elements may be sub-relevant. For instance, transport in the micro-world is strongly influenced by fluctuations and random perturbations. In some systems, like biological cells bressloff2013 (), they can even play a dominant role and a typical situation is that motion of particles is strongly damped. This fact justifies the use of an overdamped dynamics for which the particle inertial term can be formally neglected in comparison to the dissipation term ( is the particle mass, is the friction coefficient and dot denotes a differentiation with respect to time ). Omission of the inertial term enormously simplifies the modelling and in many cases allows for an analytical solutions of the corresponding Fokker-Planck equation. However, properties and features which are allowed to occur in systems with inertia can completely disappear when the inertial term is put to zero. Certainly a more correct approach in such a situation is to include the inertial term and use a technique of mathematical sequences of smaller and smaller dimensionless mass. Our main objective is to investigate impact of inertia on negative mobility of a Brownian particle moving in one-dimensional periodic systems. It is known that in such setups the inertial term is one of the key ingredients for the occurrence of this form of anomalous transport speer2007 (); spiechowicz2016pla () and negative mobility is absent for the overdamped dynamics when . We address the question whether it is still possible to observe the negative mobility phenomenon in strongly dissipative systems. In doing so, we first formulate the model and introduce the quantities of interest. Then we investigate the general transport behaviour as a function of model parameters and detect the optimal dimensionless mass for the presence of negative mobility. In the next part we demonstrate three mechanisms responsible for the emergence of this anomalous transport phenomenon: deterministic chaotic, thermal noise induced and deterministic non-chaotic. Finally, we discuss impact of inertia on the directed long time particle velocity and provide some conclusions.

## 2 Model

The model of a Brownian particle moving in a one-dimensional periodic landscape has been already well established in the literature risken (). It has been used to explore a wide range of phenomena including ratchet effects spiechowicz2014prb (); spiechowicz2015njp (), noise induced transport luczka1995 (); spiechowicz2015physscripta (), the negative mobility machura2007 (), the enhancement of transport spiechowicz2014pre () and diffusion spiechowicz2015chaos (), the anomalous diffusion spiechowicz2015pre (); spiechowicz2016scirep () and the non-monotonic temperature dependence of diffusion spiechowicz2016njp (); spiechowicz2017chaos (). Here, we consider exactly the same model as in machura2007 (): a classical inertial Brownian particle of mass , which moves in a spatially periodic potential of period and is subjected to both an unbiased time-periodic force of amplitude and angular frequency and an external static force . Dynamics of such a particle is described by the following Langevin equation machura2007 ()

(1) |

where prime denotes a differentiation with respect to the particle coordinate . Thermal fluctuations due to the coupling of the particle with the thermal bath of temperature are modelled by Gaussian white noise of zero mean and unity intensity, namely

(2) |

The noise intensity factor (where is the Boltzmann constant) follows from the fluctuation-dissipation theorem marconi2008 () and ensures the canonical equilibrium Gibbs state when and . The potential is assumed to be in a symmetric form with the period and the barrier height , namely,

(3) |

There exists a wealth of physical systems that can be described by the Langevin equation (1). An important cases that come to mind are the semiclassical dynamics of a phase difference across a resistively and capacitively shunted Josephson junction blackburn2016 () and a cold atom moving in an optical lattice denisov2014 (); lutz2013 (). Other examples include superionic conductors fulde1975 (), dipoles rotating in external field coffey2004 (), charge density waves gruner1981 () and adatoms on a periodic surface guantes2001 ().

### 2.1 Scaling and dimensionless Langevin equation

Since only relations between scales of length, time and energy are relevant for the observed phenomena, not their absolute values, we next formulate the above presented equation of motion in its dimensionless form. This can be achieved in several ways machura2008 (). Because investigation of impact of the particle inertia on the system dynamics is our main goal, in the present consideration we propose the use of the following scales as the characteristic units of length and time machura2008 ()

(4) |

Under such a procedure the Langevin equation (1) takes the dimensionless form spiechowicz2016njp ()

(5) |

In this scaling, the dimensionless mass is

(6) |

where the second characteristic time is and the dimensionless friction coefficient is . Other parameters are: , , . The rescaled thermal noise reads and assumes the same statistical properties as , namely and . The dimensionless noise intensity is the ratio of thermal energy and half of the activation energy the particle needs to overcome the nonrescaled potential barrier. The dimensionless potential possesses the period and the barrier height . From now on, we will use only the dimensionless variables and therefore, in order to simplify the notation, we will omit the hat notation in the above equation.

### 2.2 Quantities of interest

In the present study we are particularly interested in the impact of inertia on properties of directed transport of particles in the stationary state. In the dimensionless formulation (5) it can be realized by changing the dimensionless mass (6). The case corresponds to overdamped dynamics and the setting represents the strong damping regime, which means that . The characteristic time is obtained from a particular form of Eq. (1), i.e. and has the interpretation of the relaxation time of the velocity of the free Brownian particle. The parameter is extracted from the equation which can be viewed as the characteristic time to travel a distance from a maximum of the potential to its minimum in the overdamped case (it is not exactly this time which is infinite in the considered case but scales it). It is remarkable that parameters of the potential such as its barrier height and period are crucial for controlling the regimes of weak or strong damping. For instance, if and are fixed and the system is in a weak damping regime , the transition to the strong damping case can be achieved by lowering the barrier height and lengthening the period of . We have checked that for values and smaller the system (5) can be considered to be in the strong damping regime.

Due to the presence of the external time-periodic driving , as well as the friction term , the particle velocity approaches a unique non-equilibrium asymptotic long time state, in which it is characterized by a temporally periodic probability density. This latter function has the same period as the driving jung1993 (). Therefore, the first statistical moment of the instantaneous particle velocity assumes for an asymptotic long time regime the form of a Fourier series over all possible harmonics jung1993 ()

(7) |

where is the directed (time independent) velocity, while denote time periodic higher harmonics of vanishing time-average over the fundamental period . The observable of foremost interest in this study is the directed transport component , which due to the mentioned particular decomposition can be obtained in the following way

(8) |

where indicates averaging over all realizations of thermal noise as well as over initial conditions for the position and the velocity . The latter is obligatory for the deterministic limit when dynamics may be non-ergodic and results can be affected by specific choice of initial conditions spiechowicz2016scirep ().

Due to the multidimensionality of the parameter space of the considered model, as well as its nonlinearity, the force-velocity curve is typically a nonlinear function of the applied bias . From the symmetries of the underlying Langevin equation (5) it follows that this observable is odd as a function of the external static force , i.e. and in consequence denisov2014 (). This is in clear contrast to the case of a ratchet mechanism, which exhibits the finite directed transport even at the vanishing static bias when hanggi2009 (). Since the observable of our interest is symmetric around , we limit our consideration to the positive bias . Then, for sufficiently small values of the external force the directed transport velocity is usually its increasing function. Such regimes correspond to the normal, expected transport behaviour. However, in the parameter space there are also regimes for which the particle moves on average in the direction opposite to the applied bias, i.e. for , exhibiting anomalous transport behaviour in the form of the negative mobility phenomenon machura2007 (); speer2007 (). It has been already shown that there are two fundamentally various mechanisms responsible for negative mobility in this setup, (i) generated by chaotic dynamics and (ii) induced by thermal equilibrium fluctuations machura2007 (). The latter situation is nevertheless rooted in the sophisticated evolution of the corresponding deterministic system described by Eq. (5) with . Its three-dimensional phase space is minimal for chaotic evolution, which is important for negative mobility to occur.

For the considered deterministic system with there are three Lyapunov exponents , and . It can be easily checked that the system is dissipative, i.e. the phase space volume is contracting during the time evolution. Therefore the sum of all Lyapunov exponents must be negative gollub ()

(9) |

One of the exponents, say and the other, say . If the system is chaotic must be positive indicating divergence of the trajectories. Therefore to detect chaotic behaviour of the system it is sufficient to calculate the maximal Lyapunov exponent and check whether it is larger than zero boffetta2002 ().

## 3 Numerical simulation

Unluckily, the Fokker-Planck equation corresponding to the Langevin equation (5) cannot be handled by any known analytical methods. For this reason, in order to analyse transport properties of the system, we carried out comprehensive numerical simulations. We integrated the Langevin equation (5) by employing a weak version of the stochastic second order predictor corrector algorithm with a time step typically set to about . We chose the initial coordinates and velocities equally distributed over the intervals and , respectively. The quantities of interest were ensemble averaged over different trajectories, which evolved over periods of the external harmonic driving. All numerical calculations were performed by the use of CUDA environment implemented on a modern desktop GPU. This gave us possibility to speed up the computations up to a factor of the order times as compared to a common present day CPU method. Details on this promising scheme can be found in Ref. spiechowicz2015cpc ().

Dynamics described by Eq. (5) is characterized by a 5-dimensional parameter space , the detailed exploration of which is a very challenging task even for our innovative computational method. However, we focus on the impact of the particle inertia on the anomalous transport processes occurring in this setup. This task is very tractable numerically with the currently available hardware. We start our analysis by looking at the deterministic system . We set the bias to a low value and check how the directed velocity depends on the remaining parameters. In doing so we performed scans of the following area at a resolution of 200 points per dimension to determine the general behaviour of the system. Our results reveal that negative mobility is not present for and . This is in agreement with the approximate solutions of Eq. (5). In the limit of low frequencies an adiabatic approximation is valid chi1990 (), while for high frequencies a solution can be formulated in terms of Bessel functions kautz1996. Moreover, there is no net transport for .

## 4 Results

### 4.1 General behaviour of the system

The study of various aspects of transport in the system (5) has been presented elsewhere machura2007 (); speer2007 (). Here, we focus our analysis primarily on the relationship between inertia and negative mobility. In Fig. 1 we present the asymptotic long time directed velocity depicted as a function of the amplitude and angular frequency of the external harmonic driving , for the strongly damped Brownian particle with the dimensionless mass and different values of the external bias . Surprisingly, despite the fact that the inertia is one order smaller than the dimensionless friction coefficient , there are regions in the parameter space of the system where negative mobility occurs. They form a band-like structure. The stripes of negative velocity are interspersed with the stripes of positive velocity and the difference in its magnitude in the neighbouring regions can be significant. This suggest that the system may be very sensitive to a small change of values of parameters. It can be observed that larger values of lead to reduction of the negative mobility areas towards the lower and . At the same time the bands of negative velocity become wider and more intense. In contrast, the regions of positive mobility increase population and supersede their negative counterparts. Overall, as the bias is increased the band structure seems to zoom out until the directed velocity points to the direction given by and becomes almost constant in the whole map. In Fig. 2 we depict the same characteristic but for different values of thermal noise intensity with fixed and . One can expect that thermal noise perturbs deterministic dynamics. We observe that larger areas of negative mobility are relatively stable with respect to increasing temperature, while the smaller ones disappear more quickly. Thermal noise first blurs the band-like structure of negative mobility areas, erasing the finer details of the regions visible in the deterministic case . It seems to be obvious since thermal noise enables random transitions between deterministically coexisting basins of attraction. For high enough temperatures, negative mobility disappears completely. A careful inspection of Fig. 2 reveals that there are regions in the parameter space where the directed velocity is zero or positive in the deterministic case, but becomes negative upon the introduction of noise. This fact suggests that thermal fluctuations may induce negative mobility or reverse its sign even for the strongly damped Brownian particle. Finally, in Fig. 3 we present the directed velocity versus the amplitude and the angular frequency for different values of the mass with and . The stripes of negative mobility move towards lower values of and higher values of as is increased. The band-like structure changes its inclination and the regions becomes more horizontally oriented. Moreover, for larger masses some new negative mobility bands appear, while at the same time the negative velocity tends to disappear in other regions. This effect suggest that there should exist an optimal mass for which the occurrence of negative mobility is mostly pronounced.

### 4.2 Optimal mass for the presence of negative mobility

Since the impact of the mass seems to be non-trivial, an interesting task is to find the value of for which the presence of negative mobility is the most common. Our numerical scans of the parameter space allowed us to determine this value. The result is shown in Fig. 4. In panel (a) we depict fraction of the negative mobility area in the analysed parameter space for the deterministic system and different values of the bias . A perhaps surprising finding is that for the small values of the static load the optimal mass for the presence of negative mobility is . It is significantly less than the magnitude of dimensionless friction coefficient, which for the employed scaling is equal to unity . This fact indicates that the friction plays prevalent role for the emergence of negative mobility. Moreover, even for the very strongly damped Brownian particle the area of negative mobility is non-zero and relatively large. In both limiting regimes of the overdamped and the underdamped motion there are no regions of negative mobility. Ipso facto we confirmed numerically the no-go theorem formulated in Ref. speer2007 (). This conclusion is valid also for larger values of the bias , however, then the optimal mass for the presence of negative mobility becomes shifted towards higher, e.g. for . Moreover, as the static load is increased the overall occurrence of negative mobility in the analysed parameter space is decreased and the presented curves come to be more flattened. In panel (b) we show the same characteristic but for the fixed bias and different values of thermal noise intensity . As temperature grows the regions of negative mobility in the parameter space tend to contract which is illustrated in the figure. Apart from this fact for stronger thermal noise, c.f. the case , the curve becomes noticeably more bimodal. This observation is most likely due to parameter regimes for which negative mobility is induced by thermal fluctuations.

### 4.3 The mechanisms of negative mobility

To gain further insight into the nature of negative mobility in this system, as the next step we identify, exemplify and analyse the mechanisms standing behind emergence of this phenomenon. In Fig. 5 we present the regime of parameters for which negative mobility occurs on grounds which are rooted solely in the complex, strongly damped and deterministic chaotic dynamics. Panel (a) depicts the directed velocity versus the external static bias . It is a nonlinear function without any obvious relation to the magnitude of the force . Clearly, there are two windows of the latter parameter for which negative mobility is observed. The first starts at and ends at . The second is present for the bias larger than approximately . In panel (b) we present the directed velocity of the Brownian particle as a function of thermal noise intensity for the fixed bias which corresponds to the sharp minimum of depicted in the panel (a). In the limiting case of the very low temperatures the measured directed velocity is less than zero indicating that negative mobility has its origin in the complex deterministic dynamics. For increasing temperature the directed velocity grows as well up to the critical thermal noise intensity , for which the Brownian particle response reverses its sign . In the high temperature limit (not depicted) all forces in the right hand side of Eq. (5) become negligible in comparison to thermal noise and thus the directed velocity vanishes completely . Panel (c) of the same figure presents the bifurcation diagram of the latter quantity illustrated as the function of the external bias for the deterministic system with . Each blue dot represents an attractor for the asymptotic long time directed velocity . For almost all values in the considered range of the bias there is the continuum of the directed velocity solutions. This fact suggests that the system is predominantly chaotic in this interval. We confirm this hypothesis in panel (d) where we depict the maximal Lyapunov exponent for the deterministic system described by Eq. (5) with versus the biasing force . Accordingly, this quantity is positive in almost entire considered interval of the parameter . In particular, it is so for the values of corresponding to negative mobility. Therefore, we conclude that in the presented parameter regime this phenomenon is induced solely by the chaotic deterministic dynamics of the system given by Eq. (5). Such a mechanism has been already reported in literature machura2007 (); speer2007 (), however, here we prove that it may operate also for the strongly damped Brownian particle .

The second mechanism of the emergence of negative mobility is exemplified in Fig. 6. In panel (a) we present the directed velocity of the strongly damped Brownian particle versus the external static bias for thermal fluctuations intensity . In this case a very small amount of noise yields negative mobility in the linear response regime, i.e. for small forces . For larger values of the bias the directed velocity is positive and the particle moves in the direction pointed by the static force . Interestingly, there is an optimal value of for which negative mobility is most pronounced. To gain further insight into the origin of the discussed anomaly in the presented regime, in the neighbouring panel (b) we study the directed velocity as a function of thermal noise intensity . Contrary to the previously presented case, here at low temperature the Brownian particle velocity is positive. The above described negative mobility manifests itself only in finite interval of temperature . Further increase of thermal noise intensity leads to disappearance of this phenomenon. Although a solely noise induced negative mobility can occur only under impact of thermal fluctuations, the underlying relevant mechanism is strongly influenced by the deterministic dynamics as it was already shown in Ref. machura2007 ().

Finally, in Fig. 7 we present a case of negative mobility induced by the deterministic non-chaotic dynamics of the system. Panel (a) illustrates the directed velocity as a function of the external bias for the deterministic case and in the strong damping regime . For very small values of the force the directed velocity oscillates around zero. When the bias is of moderate magnitude there is a window for which the particle response is opposite to the applied constant perturbation, so we detect there negative mobility. Further increase of the external force rapidly reverses the particle current and causes its monotonic growth. In panel (b) we study the impact of thermal fluctuations on in the parameter regime with corresponding to the minimal plateau depicted in the plot (a). Indeed, for the deterministic limit of the dynamics the directed velocity is negative suggesting that the observed phenomenon of negative mobility has the deterministic origin. An increase of thermal noise intensity ceases this effect. For this parameter regime there is a surprisingly simple structure of attractors for the directed velocity which is visualized in the panel (c) in the deterministic case with . In the considered interval of the external force there are two asymptotically stable solutions corresponding to . Notably, in the bias window where negative mobility is observed only the attractor survives. This unexpected simplicity of solutions suggests that in the considered parameter regime the deterministic dynamics is non-chaotic, nonetheless still exhibits negative mobility. Our finding is confirmed in panel (d) where we depict the maximal Lyapunov exponent versus the bias for the system with . An interesting observation is that the dynamics is generally chaotic () when two attractors coexist and it is non chaotic with in the window where negative mobility emerges. To the best of the authors knowledge such a mechanism has never been reported. It may potentially open a possibility of observation of the negative mobility phenomenon for the discontinuous or non-ergodic one-dimensional nonequilibrium overdamped dynamics, corresponding to the formal substitution in Eq. (5) speer2007 ().

### 4.4 Impact of inertia

To conclude this section we present in Fig. 8 the representative dependence of the directed velocity on the inertia for the deterministic and noisy system. In panel (a) the amplitude and frequency of the periodic driving are the same as in the previous figures 5-7. Here, the most pronounced negative mobility is observed for the mass . For specifically tailored parameter sets this phenomenon could be detected for even smaller mass. We exemplify this situation in panel (b) where it is observed for which is indeed the regime of very strong damping.

We observe that the system response is very sensitive to even smallest changes of the particle mass . Moreover, there are multiple reversals of the sign of the directed velocity which are characteristic for a massive setup driven by the external harmonic force bartussek1994 (); jung1996 (). This finding can be utilized to particle sorting eichhorn2010 (). For instance, one can see from the above figure that particles with different masses can easily be guided into opposite direction by a suitable choice of the system parameters. In addition, we want to point out that the occurrence of the three presented mechanisms of negative mobility are controlled by the magnitude of the particle inertia. Depending on its value the deterministic chaotic, thermal noise induced and the deterministic non-chaotic anomalous transport can be observed. It is worth to explicitly note that a tiny change from to transforms the nature of negative mobility effect from deterministic non-chaotic to thermal noise induced, c.f. Fig. 6 and Fig. 7. Therefore, the outlined mechanisms of the occurrence of negative mobility are also very sensitive to alteration of the particle inertia.

## 5 Conclusion

In this work we investigated the impact of inertia on transport properties of a Brownian particle moving in a periodic symmetric structure, which in addition is exposed to a harmonic ac driving as well as a constant bias. The parameter space is 5-dimensional and its complete numerical exploration is far beyond the scope of this work. We first analysed the general behaviour of the directed velocity as a function of the amplitude and the angular frequency of the driving, for selected values of the remaining system parameters: the particle inertia , the external bias and thermal noise intensity . These results reveal especially that the negative mobility phenomenon emerges also for the strongly damped motion of a Brownian particle when the dissipation dominates over inertia. Our scans of the parameter space allowed us to determine the optimal mass for the presence of the negative mobility phenomenon. By observing that the fraction of the negative mobility area in the analysed parameter space disappears for both the limiting cases of the overdamped and underdamped motion, we confirmed with precise numerics the no-go theorem formulated in Ref. speer2007 (). We gained further insights into the origin of negative mobility in this system by revealing three classes of mechanisms responsible for this anomalous transport process. It can (i) be caused by the complex deterministic chaotic dynamics of the system or (ii) induced by the thermal noise, or (iii) associated with the deterministic, yet non-chaotic system evolution. In particular, according to the best of authors knowledge, the latter origin has never been reported before. It may provide guidance to the possibility of observation of the negative mobility phenomenon for the discontinuous or non-ergodic one-dimensional nonequilibrium overdamped dynamics when the particle inertia is fixed to zero . This case is of fundamental importance from the point of view of biological systems, all of which in situ operate in strongly fluctuating environments. Finally, we depict the illustrative impact of the particle inertia on its transport properties to unravel its spectacular sensitiveness to variation of this parameter. We detect the phenomenon of multiple velocity reversals which may constitute a cornerstone of particle sorting. Moreover, a small change in the particle inertia can radically alter the mechanism of negative mobility.

Our results can readily be experimentally tested with a single Josephson junction device or cold atoms moving in an dissipative optical lattice. Suitable parameter values, for which the above effects are predicted to occur, are accessible experimentally.

## 6 Acknowledgement

The work was supported by the NCN grant 2015/19/B/ST2/02856 (J. S. & J. Ł.).

## References

- (1) Denisov S, Flach S, Hänggi P. Tunable transport with broken space-time symmetries. Phys. Rep. 2014;538:77-120
- (2) Hänggi P and Marchesoni F. Artificial Brownian motors: Controlling transport on the nanoscale. Rev. Mod. Phys. 2009;81:387
- (3) Machura Ł et al. Absolute negative mobility induced by thermal equilibrium fluctuations. Phys. Rev. Lett. 2007;98:40601
- (4) Nagel J et al. Observation of negative Absolute Resistance in a Josephson Junction. Phys. Rev. Lett. 2008;100:217001
- (5) Kostur M, Łuczka J and Hänggi P. Negative mobility induced by colored thermal fluctuations. Phys. Rev. E 2009;80:051121
- (6) Spiechowicz J, Łuczka J and Hänggi P. Absolute negative mobility induced by white Poissonian noise. J. Stat. Mech. 2013;P02044
- (7) Spiechowicz J, Hänggi P and Łuczka J. Brownian motors in the microscale domain: Enhancement of efficiency by noise. Phys. Rev. E 2014;90:032104
- (8) Spiechowicz J, Łuczka J and Machura Ł. Efficiency of transport in periodic potentials: dichotomous noise contra deterministic force. J. Stat. Mech. 2016;054038
- (9) Hennig D. Current control in a tilted washboard potential via time-delayed feedback. Phys. Rev. E 2009;79:041114
- (10) Du L and Mei D. Time delay control of absolute negative mobility and multiple current reversals in an inertial Brownian motor. J. Stat. Mech. 2011;P11016
- (11) Mulhern C. Persistence of uphill anomalous transport in inhomogeneous media. Phys. Rev. E 2013;88:022906
- (12) Dandogbessi B, Kenfack A. Absolute negative mobility induced by potential phase modulation. Phys. Rev. E 2015;92:062903
- (13) Du L and Mei D. Absolute negative mobility in a vibrational motor. Phys. Rev. E 2012;85:011148
- (14) Januszewski M and Łuczka J. Indirect control of transport and interaction-induced negative mobility in an overdamped system of two coupled particles. Phys. Rev. E 2011;83:051117
- (15) Machura Ł et al. Two coupled Josephson junctions: dc voltage controlled by biharmonic current. J. Phys. Condens. Matt. 2012;24:085702
- (16) Ghosh PK, Hänggi P, Marchesoni F and Nori F. Giant negative mobility of Janus particles in a corrugated channel. Phys. Rev. E 2014;89:062115
- (17) Malgaretti P, Pagonabarraga I and Rubi JM. Entropic Electrokinetics: Recirculation, Particle Separation, and Negative Mobility. Phys. Rev. Lett. 2014;113:128301
- (18) Sarracino A, Cecconi F, Puglisi A and Vulpiani A. Nonlinear Response of Inertial Tracers in Steady Laminar Flows: Differential and Absolute Negative Mobility. Phys. Rev. Lett. 2016;117:174501
- (19) Bressloff PC and Newby JM. Stochastic models of intracellular transport. Rev. Mod. Phys. 2013;85:135
- (20) Speer D, Eichhorn R and Reimann P. Transient chaos induces anomalous transport properties of an underdamped Brownian particle. Phys. Rev. E 2007;76:051110
- (21) Spiechowicz J, Kostur M and Łuczka J. Comment on ”Absolute negative mobility in a one-dimensional overdamped system”. Phys. Lett. A 2016;380:1499
- (22) Risken H. The Fokker-Planck Equation. Berlin: Springer-Verlag; 1984
- (23) Łuczka J, Bartussek R and Hänggi P. White noise induced transport in periodic structures. Europhys. Lett. 1995;31:431-436
- (24) Spiechowicz J and Łuczka J. Poissonian noise assisted transport in periodic systems. 2015;165:014015
- (25) Spiechowicz J, Hänggi P and Łuczka J. Josephson junction ratchet: The impact of finite capacitances. Phys. Rev. B 2014;90:054520
- (26) Spiechowicz J and Łuczka J. Efficiency of a SQUID ratchet driven by external current. New J. Phys. 2015;17:023504
- (27) Spiechowicz J and Łuczka J. Josephson phase diffusion in the SQUID ratchet. Chaos 2015;25:053110
- (28) Spiechowicz J and Łuczka J. Diffusion anomalies in ac-driven Brownian ratchets. Phys. Rev. E 2015;91:062104
- (29) Spiechowicz J, Łuczka J and Hänggi P. Transient anomalous diffusion in periodic systems: ergodicity, symmetry breaking and velocity relaxation. Sci. Rep. 2016;6:30948
- (30) Spiechowicz J, Talkner P, Hänggi P and Łuczka J. Non-monotonic temperature dependence of chaos-assisted diffusion in driven periodic systems. New J. Phys. 2016;18:123029
- (31) Spiechowicz J, Kostur M and Łuczka J. Brownian ratchets: How stronger thermal noise can reduce diffusion. Chaos 2017;27:023111
- (32) Marconi UMB, Puglisi A, Rondoni L and Vulpiani A. Fluctuationâdissipation: response theory in statistical physics. Phys. Rep. 2008;461:111-195
- (33) Blackburn JA, Cirillo M and Gronbech-Jensen N. A survey of classical and quantum interpretations of experiments on Josephson junctions at very low temperatures. Phys. Rep. 2016;611:1
- (34) Lutz E and Renzoni F. Beyond Boltzmann-Gibbs statistical mechanics in optical lattices. Nat. Phys. 2013;9:615
- (35) Fulde P et al. Problem of Brownian motion in a periodic potential. Phys. Rev. Lett. 1975;35:1776
- (36) Coffey W, Kalmykov Y and Waldron J. The Langevin equation. Singapore: World Scientific; 2004
- (37) Gruner G, Zawadowski A and Chaikin P. Nonlinear conductivity and noise due to charge-density-wave depinning in . Phys. Rev. Lett. 1981;46:511
- (38) Guantes R, Vega J and Miret-Artes S. Chaos and anomalous diffusion of adatoms on solid surfaces. Phys. Rev. B 2001;64:245415
- (39) Machura Ł, Kostur M and Łuczka J. Transport Characteristics of Molecular Motors. Biosystems 2008;05:033
- (40) Jung P. Periodically driven stochastic systems. Phys. Rep. 1993;234:175
- (41) Baker GL and Gollub JP. Chaotic Dynamics: An Introduction. Cambridge: Cambridge University Press; 1996
- (42) Boffetta G, Cencini M, Falcioni M and Vulpiani A. Predictability: a way to characterize complexity. Phys. Rep. 2002;356:367
- (43) Spiechowicz J, Kostur M and Machura Ł. GPU accelerated Monte Carlo simulation of Brownian motors dynamics with CUDA. Comp. Phys. Commun. 2015;191:140
- (44) Chi CC and Vanneste C. Onset of chaos and dc current-voltage characteristics of rf-driven Josephson junctions in the low-frequency regime. Phys. Rev. B 1990;42:9875
- (45) Bartussek R, Hänggi P and Kissner JG. Periodically rocked thermal ratchets. Europhys. Lett. 1994;28:459
- (46) Jung P, Kissner JG and Hänggi P. Regular and chaotic transport in asymmetric periodic potentials: inertia ratchets. Phys. Rev. Lett. 1996;76:3436
- (47) Eichhorn R, Regtmeier J, Anselmetti D and Reimann P. Negative mobility and sorting of colloidal particles. Soft Matter 2010;6:1858