Brownian particles driven by spatially periodic noise

Abstract We discuss the dynamics of a Brownian particle under the influence of a spatially periodic noise strength in one dimension using analytical theory and computer simulations. In the absence of a deterministic force, the Langevin equation can be integrated formally exactly. We determine the short- and long-time behaviour of the mean displacement (MD) and mean-squared displacement (MSD). In particular, we find a very slow dynamics for the mean displacement, scaling as \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$t^{-1/2}$$\end{document}t-1/2 with time t. Placed under an additional external periodic force near the critical tilt value we compute the stationary current obtained from the corresponding Fokker–Planck equation and identify an essential singularity if the minimum of the noise strength is zero. Finally, in order to further elucidate the effect of the random periodic driving on the diffusion process, we introduce a phase factor in the spatial noise with respect to the external periodic force and identify the value of the phase shift for which the random force exerts its strongest effect on the long-time drift velocity and diffusion coefficient Graphical Abstract


Introduction
Dating back to the important paper by Einstein in the annus mirabilis 1905 [1], the dynamics of Brownian particles has been in the focus of statistical physics for more than 100 years now [2].The constant interest in Brownian particles is basically inspired by two facts: First, their stochastic description requires fundamental principles such as the Langevin or Smoluchowski picture such that they serve as paradigmatic models which can be made systematically more complex.Second, there is a variety of excellent realizations of Brownian particles including mesoscopic colloidal particles in suspension [3], random walkers in the macroscopic world (such as [4]) and in the microscopic biological context [5], and even elements of the stock exchange market [6].This facilitates a direct comparison of the stochastic averages between the stochastic modelling and real experimental data.
In its simplest one-dimensional form, the most basic model Langevin equation for a particle trajectory x(t) a e-mail: breoni@hhu.de as a function of time t is ẋ(t) = √ Dη(t) in which η(t) is white noise with zero mean and variance η(t)η(t ) = δ(t − t ) and D > 0 is the diffusion constant.Here, ... denotes a noise average.With the initial position x(t = 0) = x 0 , the mean displacement vanishes due to symmetry, x(t)−x 0 = 0, and the mean-squared displacement is purely diffusive, (x(t) − x 0 ) 2 = 2Dt.Clearly, this basic equation can be extended towards more complicated situations including an additional static external force, time-dependent external forcing, higher spatial dimensions, and many interacting particles, see [7][8][9][10] for some reviews.One particularly interesting way to extend the equation is to generalize it to a situation of multiplicative noise, where the noise strength is a positive function D(x, t).While the case where D is only an explicit function of time t is well studied, for example in the context of Brownian ratchets [11][12][13][14][15] and heat engines [16][17][18][19][20], in this work we focus on the case where we have a spatially dependent noise strength [21][22][23][24] modelled by a positive function D(x), i.e. a space-dependent diffusion coefficient, such that the most basic model for such processes is given by the Langevin equation ẋ(t) = D(x(t))η(t). (1) The special case of multiplicative noise where ẋ(t) = −κx(t)η(t) with positive κ [25], which is somehow related to this model, documents already that the spatial dependence of the noise gives rise to fundamentally new mathematical concepts also known as the Itô-Stratonovich problem [26].The mathematical difficulties associated with the formal treatment of Eq.( 1) are subject to intense discussion, see, e.g., the recent work by Leibovich and Barkai for the specific choice of D(x) as a power-law [27] and numerous other studies [24,[28][29][30][31][32][33][34][35][36][37][38][39][40].
In this paper we consider a variant of this model in the context of the discussion of particle motion in tilted potentials.There is a large literature on this topic, see [41][42][43][44][45][46][47][48][49][50][51][52][53][54].Following the original suggestion by Büttiker [21] and Landauer [22] the spatially-varying thermal noise source can be combined with a ratchet potential, as e.g. recently discussed by [55].Our model considers overdamped Brownian particles subject to an oscillating tilted potential and a space-dependent periodic noise amplitude with the same wave vector k as the force; furthermore, we will ultimately also allow a shifted phase φ in the random force.In its general form, the model is given by the Langevin equation where V (x) := −F 0 (x + sin(kx)/k) is the potential, T (x) := T 0 (1 + ν cos(kx + φ)) 2 is the space-dependent noise strength, γ is the friction coefficient, F 0 is the tilting force, T 0 is a reference temperature, η(t) is a white noise, as introduced before, and and ν are dimensionless parameters.The critical tilt in this model arises when = 1.In order keep the noise strength differentiable everywhere and its phase in a fixed frame we consider 0 ≤ ν ≤ 1.The period of both the force and the noise will be L = 2π/k.We remark that the case ν = 1 plays a special role insofar as there are special positions at which the noise is zero.In absence of forces, the particle will therefore never cross these positions but stay confined within a periodicity length L.
Our goal in this paper is to describe the particle dynamics as functions of , ν and φ, either in the vicinity of the critical tilt, or in the absence of the deterministic force, F 0 = 0, i.e. in the purely spatial random noise case.Among our main results are the very slow dynamics in the relaxation of the mean displacement (MD) and mean-squared displacement (MSD) for long times in the F 0 = 0 case and an essential singularity in the stationary current for F 0 = 0 and ν 1.In the case of the full model, we build upon the results of [21] by also considering extreme temperature oscillations where the noise strength vanishes (ν = 1) and adding an external driving force, while we expand on [41] by finding a theoretical approximation for both the long-time drift v L and diffusion constant D L and the phase value φ for which we have the largest increase of v L and D L for = 1 and ν = 0. Our results have been obtained both with numerical and analytical methods.
The paper is organized as follows: in the beginning we focus on the free case, for which we study the shortand long-time behaviour of MD and MSD, then we proceed with the full model, including the tilted potential, for which we study the stationary distribution and the dependence of long time diffusion and drift on φ and ν.Finally, we summarize the results obtained and discuss possible experimental realizations of the model.

Free particle case
In the case of a vanishing external force (F 0 = 0), the Langevin equation ( 2) now reads as where we set φ = 0 without loss of generality.We decided to approach this problem using the Stratonovich interpretation.For a given representation of the noise, this equation can be solved by direct integration in the particular case of periodic boundary conditions (PBC) in which we identify x(t) ± L with x(t).The PBC correspond to a ring-like geometry of the one-dimensional system.
and the limit of this solution for ν → 1 is We remark here that in the case with no boundaries, i.e. when we let the particle diffuse through the whole x axis, the analysis is harder and we were not able to find an analytical expression except for the special case ν = 1.In this limit PBC and the no boundaries case are identical as the particle can never trespass the points where the noise is zero.
Equations ( 4) and ( 5) can be used to express noiseaverages of any power of displacement.For an arbitrary moment M n (t) := (x(t) − x 0 ) n we obtain for ν = 1 and for ν = 1.Since we are going to focus on the mean displacement x(t) − x 0 and the mean-squared displacement (x(t) − x 0 ) 2 , we write the expressions for these two moments (n = 1, 2) explicitly: for ν = 1 and for ν = 1.

Short-time behavior
We can use equations (8)(9)(10)(11) to extract the short-time behavior of the MD and MSD.Expanding the integrand in powers of t using a Taylor series and integrating the terms separately we obtain for the MD: and for the MSD In the special limit ν = 1 we also add the second order correction as: Clearly, the first-order correction of ( 14) and ( 15) coincides with equations ( 12) and ( 13) in the limit ν → 1.
Moreover for ν = 0 we recover the white noise case solved by Einstein [1].
0.4 0.2 0.0 0.2 0.4 Fig. 1 Effective potential of the mean displacement V M , obtained from the short-time drift of the mean displacement, and space-dependent noise T (x) for ν = 0.5 and φ = 0 as functions of space x.While the averaged MD tends to the minima of V M , where the noise strength T (x) is largest, individual trajectories spend most of their time around the maxima of V M .
We now define an effective potential of the mean displacement such that a particle subject to this potential and constant white noise will experience the same average drift as a particle in a space-dependent noise landscape.In other words, following the spirit of the mapping proposed by Büttiker [21], the effective force resulting from this potential can be viewed as a substitute source for the drift when only white noise is considered.Hence we define this force F M (x) up to a friction coefficient prefactor γ as the first coefficient of the short-time expansion of the MD as follows The effective potential of the mean displacement is then defined by (18) This potential is shown in Fig. 1.Even though this potential is defined just by the short-time expansion of the MD, it is still significant for any finite time, as the particle is overdamped and feels at every time a shorttime drift depending only on its position.As a result, the MD of a particle subject to this potential and white noise can be perfectly mapped to the MD of a free particle with space-dependent noise.
While the average mean displacement behaves according to V M , moving over time towards the regions where V M is smaller and the noise strength is larger, we want to stress that individual trajectories will not accumulate in the minima of V M but will instead freely move over all the domain, spending most of their time in the maxima of V M instead.This is because when particles reach such low-noise regions they take a longer time escaping, as their fluctuations there are severely reduced.

Dynamics for finite and long times
Now we explore the behavior of the MD and MSD for finite and long times.First we present an asymptotic analysis for the special case ν = 1.Then we use a numerical solution of the integrals in ( 8) and ( 9) as well as computer simulations of the original Langevin equation to obtain data for finite times and arbitrary ν.

Asymptotic analysis for ν = 1 for long times
Here we present an asymptotic analysis for the MD and MSD by starting from Eq.( 10) and using the asymptotic approximation for large θ.We now expand arctan and insert this expansion in Eq.( 10) to obtain which yields As a result, the leading asymptotic behavior of x(t) is determined by the first term involving a scaling behavior of the MD in 1/ √ t.This is remarkably slow compared to typical behavior of a Brownian particle in a harmonic potential or of active Brownian motion where the MD reaches its asymptotic value exponentially in time [57-59] thus constituting an example of a very slow relaxation as induced by space-dependent noise.
Likewise an asymptotic analysis for ν = 1 yields for the long-time limit of the MSD which represents the degree of smearing of the particle distribution for long times.We want to remark that the MSD calculated from a distribution with periodic boundary conditions does not describe the effective diffusion coefficient D L in periodic systems with no boundaries, in contrast to the MD which can actually be calculated from the distribution with periodic boundary conditions even for open systems.

Computer simulations
We performed direct Brownian dynamics computer simulations of the original Langevin equations with a finite time step ∆t to obtain numerically results for the MD and MSD at any times.In order to properly simulate a system with space-dependent noise, we used the order O (∆t) Milstein scheme [60] with a time step of ∆t = 10 −3 τ , where τ := γL 2 k B T0 is a typical Brownian time scale of the system.For each simulation set we fixed the initial position x 0 within the first period [−L/2, L/2] and averaged typically over 200 trajectories of length 500τ .

MD and MSD for finite times
Data for the mean displacement and the mean position as a function of time are obtained by a numerical evaluation of the integral in Eq. ( 10) and by computer simulation.For ν = 1 results are presented in Fig. 2 together with the corresponding short-time and long-time asymptotics ( 14) and (22).The displacement starts linear in time t and saturates for long times.The mean position approaches zero slowly as a power law in time proportional to t −1/2 .For large times the statistical error in the simulation data is significant but nevertheless these data are compatible with the scaling prediction of the theory.
In order to understand the very slow behavior of the MD we note that while the MD tends to zero, i.e. to the point with largest noise, this is just an effect of averaging over particles spending most of their time at the points with the smallest noise on both sides of the x-axis: x −L/2 and x L/2.This particular mechanism explains why the MD approaches its final value so slowly, as the particles have to hop from one side to the other to symmetrize their distribution.In Fig. 3a this is clearly documented in the time evolution of the particle distribution function p(x, t), which gives the probability to find a particle after a time t at position x provided it started at time t = 0 at position x 0 .The system evolves from a single-peaked distribution around x 0 to a double-peaked distribution in ±L/2.Near the two points x = ±L/2 of zero noise the peaks are getting sharper as t → ∞ approaching to δ-peaks such that lim t→∞ p(x, t) = (δ(x − L/2) + δ(x + L/2))/2.The intuitive reason for this is that once a particle adsorbs at the points x = ±L/2 of zero noise it will never return to the region where the noise is finite.This peculiar behavior is clearly delineated from the relaxation in a symmetric double-well potential with white noise of strength T 0 .In order to reveal this, we have performed simulations for a Brownian particle in the double-well potential with two equal minima We set A := 48k B T 0 /L 4 and B := L 2 /2 in order to have the two wells in ±L/2 such that the energy barrier between the two minima is 3k B T 0 .Our simulation for (b) Fig. 3 (a) Probability density function p(x, t) for the particle position at different times t, with ν = 1 and x 0 = −0.2L.
Here we averaged over 10000 different trajectories of length 10τ .(b) Hopping rate Γ h between the two peaks in the particle distribution as a function of time t for a space-dependent noise with points of vanishing noise (ν = 1), and for a double-well potential with white noise.Here we have chosen x 0 = −0.2L.We also show the stationary state theoretical value of Γ h for the double well potential, defined as the inverse of the mean first passage time t e , derived in Eq. (26).
this white-noise reference case show that both the MD and the MSD decay exponentially in time t rather than with 1/ √ t, and hence much faster than for our case of space-dependent noise.We also defined a particle hopping rate Γ h between the two peaks of the distribution as where N h (t) is the number of times a particle hops from one peak to the other in the time interval [t/2, t].Note that the relevant time window in which hopping is considered is chosen to be proportional in time in order to improve the statistics.We have a hop whenever the particle trespasses the x = L/4 or x = −L/4 thresholds and previously was respectively in the left or right peak.
In fact, as we show in Fig. 3b, for the double-well potential, the hopping rate Γ h (t) converges to a constant for long times.This rate is maintaining the equilibrium state with a symmetrized occupation around the two minima.The rate saturates for t → ∞ to a value very close to the inverse of the mean first passage time (see for example [61]) in the double-well potential t e [62], which in our case is given by: Conversely, for our case of space-dependent noise, the hopping rate keeps decreasing as a function of time again with an inverse power law t −1/2 .This reflects the fact that the peaks of the space-dependent noise distribution keep growing indefinitely as the particles get in average closer to the points of zero noise.Now in Fig. 4 we explore the MD for the case ν = 1 where the particle crosses the position of minimal noise.
Here the boundary conditions do matter and we distinguish between no boundaries (Fig. 4a) with infinitely many oscillations and periodic boundary conditions of a ring-like geometry (Fig. 4b).While the short-time behavior is linear in time for both kind of boundary conditions, the MD saturates for long times to a finite value depending on ν and x 0 for the no boundaries case.This finite value is −x 0 for periodic boundary conditions since in this case the mean position will always end at zero due to symmetry.The asymptotic approach to zero is exponential in time as in the case of the double-well potential with noise as the particle stays mobile even when approaching the position where the noise is minimal.This is in marked contrast to the limit of ν = 1 where the particle gets immobilized at the boundaries.Now we turn to the MSD, first for the special case ν = 1 shown in Fig. 5a where boundary conditions do not matter.The MSD starts linear in time and then saturates to its long-time limit C := x 2 0 + (π/k) 2 .Its asymptotic approach to this saturation value is revealed by plotting the MSD shifted by C which decays to zero for large times, see Fig. 5b.Similar to the MD for ν = 1, we find that the asymptotic behavior is compatible with a 1/ √ t scaling.
In Fig. 6 we show the MSD for ν = 1 for both types of boundary conditions.In absence of boundary conditions (see Fig. 6a) the long-time behavior is linear in time ≈ 2D L t involving a long-time diffusion coefficient D L .Clearly the latter depends on ν but not on the initial position x 0 .This dependence is depicted in the inset of Fig. 6a.We found the empirical expression D L (ν) = D 0 (1 − ν 2 ) with D 0 = k B T 0 /γ to be a very good fit to the data.This can be regarded as a parabolic fit which fulfills the inflection symmetry in ν and the constraint D L (ν = 1) = 0.The same behavior was re-cently found in a similar system [63].
Finally, to better clarify the behaviors of the MD and MSD for ν = 1, we plot the dynamical exponents (Fig. 7) that define the scaling regimes for the MD (β, β ) and MSD (α, α ) close to their short-time and longtime limits, respectively: Both the MD and MSD for short times are linear, while for long times the scaling of the MD converges clearly to -0.5, that corresponds to 1/ √ t.Within the time window explored the MSD has not yet saturated to an ultimate dynamical exponent for long times.The asymptotics shown is compatible with a final scaling exponent of −1/2 although the approach to this final exponent is much slower for the MSD than for the MD where the saturation is clearly visible.
We remark that an algebraic asymptotic approach in the MSD was also found for equilibrium Brownian dynamics of repulsive interacting particles.Here the time-derivative of the time-dependent diffusion coefficient MSD/t scales as t −d/2 in d spatial dimensions [64-67] but the physical origins of the algebraic scaling laws are different.

Tilted potential
In this section we leave the situation in which the Brownian particle is a free particle only driven by spatiallydependent noise.We now consider the full model, including the deterministic tilted potential.We first look at the situation near the critical value of the amplitude = 1, where the tilted potential develops a plateau.The situation addressed in shown in Fig. 8.

The stationary current
Being weakly confined to a region of the deterministic potential in which the dynamics can be considered 'slow', a quasi-stationary distribution can be defined [49].The Fokker-Planck equation corresponding to the  Langevin equation, Eq. ( 1) in Stratonovich interpretation reads as ) with a(x) the force and b(x) the noise amplitude, a(x) := F 0 (1 + cos(kx)) , b(x) := 2γk B T (x) .(29) Following the discussion in [49], the dynamics near the critical tilt value for ≥ 1 is characterized by a stationary current given by the one-time integrated FPequation Defining (b(x)/2)p s (x) = p s (x) we can rewrite the last expression as with The equation can be solved with the Ansatz p s (x) = u(x) • v(x) which reduces the problem to two readily integrable first-order ordinary differential equations for u(x) and v(x).One obtains the final expression x 0 = −0.2L in theory and simulation as functions of time t.As we have already seen in Fig. 2 for the MD and Fig. 4 for the MSD, both quantities grow initially linearly in time and decay to their final limit with 1/ √ t for the MD and slower than 1/ √ t for the MSD. in which the current J s can be obtained from the normalization integral In the following we take for simplicity (setting all other constants to one) Setting b(z) := exp (−F (z)) = exp(− ln(1 + ν cos(z)), and expanding both b(z) and R(z) in Taylor series around the center of the flat region near z = L/2, the stationary current J s is given by in which the symbol ... indicates the Taylor-expanded functions, and The integration of R(z) yields a cubic polynomial, but due to cancellations the resulting expression in the exponential is Gaussian in x and cubic in y.The Gaussian integral in y can be calculated exactly, while the remaining expression in y needs to be evaluated numerically for each value of and ν.

Phase difference between noise and potential
For a tilted potential, we now explore the effect of a non-zero phase φ = 0 on the long-time behavior of the particle for different values of ν by using computer simulation.
As shown in [41], the long-time drift velocity and diffusion coefficients (v L and D L respectively) can be analytically calculated for the case ν = 0, where we set V (x) as potential.
Here the question is how the mismatch of the periodic noise and external forcing affects the long-time behavior of the particle.Intuitively one would expect that overcoming an energetic barrier is best if the maximum of the noise occurs where the external force is opposing most.Then the noise would help to bring the particle over the energetic barrier.The position where the force is opposing most is clearly given for x = L/2 + nL, where n is an integer.Then it is expected that mobility gets a maximum if the phase shift is φ = π.This is indeed what we confirm by simulation.We chose k B T 0 = 0.01F 0 L and = 1.3.The potential barrier ∆E is given by yielding ∆E 0.04F 0 L > 0.01F 0 L for = 1.3.Given these parameters, we simulated the system for different values of φ and ν and results are summarized in Fig. 9. Since to the best of our knowledge there is no easy generalization of the results in [41] for a spacedependent temperature, we have compared the simulation data with a mapping on the analytical results for v (0) L and D (0) L [41] which were obtained for a spatially constant temperature.Since the crucial position to hop over the barrier is at x = L/2 + nL where the opposing force is maximal, this represents the kinetic bottleneck for the dynamical process.Therefore it is tempting to compare our simulation results with the analytical ones where this local noise strength T (x = L/2) is inserted as a homogeneous temperature.We remark that this temperature T (x = L/2) depends both on the oscillation strength ν and the phase shift φ of T with respect to the potential.This mapping theory should work best if the particle spends most of its time close to the point x = L/2.In fact, Fig. 9 reveals that this simple mapping theory describes the simulation data well even for large ν.As a function of the phase mismatch φ, both D L and v L are enhanced when φ is between about 3  5 π and 7 5 π.Clearly around the value φ = π we find the maximal enhancement of both v L and D L .In the complementary case, the noise strength T (x) has its minimum closer to the crucial region where the opposing force is maximal, and as a result the drift velocity and diffusion are severely reduced.For ν = 1 they are even brought exactly to zero when |φ| < arcsec( ), since the particle is stuck and there is no systematic external force to drift over the positions of vanishing noise.

Conclusions and outlook
In conclusion we have presented a detailed study of a model for a Brownian particle moving in a one-dimensional environment with a space-periodic noise and under an external potential with a tilt near its critical value.In the free case we calculated the exact solution of the associated Langevin equation, and further explicitly obtained short-and long-time approximations of the MD and MSD.These results allow us to characterize the slow decay of these quantities at long times.Interesting relaxation dynamics occurs around points of vanishing noise which establish centers of growing peaks in the particle distribution, as particles slow down significantly in the neighborhoods of these points.Introducing the tilted periodic potential we first determined the stationary current for the quasi-stationary state, which for ≥ 1 displays an essential singularity for the maximal strength of the noise oscillations, ν.Finally, we determined numerically the effects of a space-periodic noise on the long-time diffusion and drift as functions of the phase difference between noise and potential φ and the strength of the noise oscillations ν, finding the largest enhancements to take place for a phase of φ = π and the maximal possible noise oscillations for ν = 1.
Our one-dimensional model with both periodic boundary conditions or no boundaries can be realized by a colloidal particle confined in a ring or a linear channel respectively by e.g.optical forces [48,[68][69][70].The space-dependent noise can be added by various means.First, one can change locally the solvent temperature.This realization has a limited applicability, since the state of the solvent can be changed drastically upon such a temperature variation.However, there are more general and more important realizations for our model.First of all, the viscosity or the friction coefficient can directly be changed without changing the ambient temperature.The solvent viscosity, for instance, can be tuned over orders of magnitude by imposed patterned substrates interacting with the solvent or even by varying the size of the colloids without changing the solvent phase.Second, space-dependent noise can stem from active internal fluctuations [72, 73] different from thermal fluctuations and can be embodied into an effective noise strength that can largely be tuned by activity [74][75][76][77][78]. Optical gradients can be used to steer activity as a function of the position, as realized and discussed in [63,[79][80][81].Another possibility is to tune the noise amplitude of skyrmions, which have a similar equation of motion [82].Last but not least, the noise can be mimicked in valuable model systems by applying randomized kicks of an external field to the particle.For example, the noise strength can largely be tuned externally without changing the solvent at all by tuning the rotational diffusion constant of the colloids [83,84].In fact, the effective diffusion constant of an active particle depends on its rotational diffusion constant, and in the limit of short persistence lengths one can indirectly tune the translational diffusion by tuning the rotational one.

Fig. 2 1 t 1 / 2
Fig.2Absolute values of the mean displacement x(t) − x 0 (a) and the mean position with a minus sign to ensure positivity − x(t) (b) for ν = 1 and x 0 = −0.2Las a function of time t.The numerical evaluation of the integral in Eq. (10) (theory) and its asymptotic short-and long-time expansions (14) and (21) are shown together with simulation data.The MD increases linear in time t for short times, while the decay to its limit scales in a very slow way with O t −1/2 .

Fig. 4
Fig. 4 Mean displacement x(t) − x 0 (a,b) as a function of time t for three values of ν = 0.25, 0.5, 0.75 and x 0 = −0.2Lfor no boundaries (a) and for periodic boundary conditions (b).

1 Fig. 5
Fig.5Absolute values of the mean-squared displacement (MSD) (x(t) − x 0 ) 2 (a) and shifted MSD−C, where C is the long time limit of the MSD (b) as a function of time t for ν = 1 and x 0 = −0.2L.The numerical evaluation the integral in Eq. (11) and its asymptotic short-and long-time expansions are shown together with simulation data.

1 Fig. 6 Fig. 7
Fig.6(a) Mean-squared displacement (MSD) (x(t) − x 0 ) 2 as a function of time t for ν = 1 and x 0 = −0.2Lboth for no boundaries (a) and for periodic boundary conditions (b).The inset shows the long-time diffusion constant D L as a function of ν for the no boundaries case.

Fig. 8
Fig. 8 Potential (a) and corresponding force (b) near the 'flat' regime = 1 as a function of the spatial coordinate x.

Fig. 9
Fig.9Long time drift (a) and diffusion (b) for = 1.3 as functions of ν and φ, simulations results and theory.For φ between about3  5 π and 7 5 π both D L and v L are enhanced, having a maximum in φ = π, while otherwise they are reduced.
Zvyagolskaya O, Bleil S, Eichhorn R, Bechinger C and Reimann P 2008 Physical Review E 77 041107 46.Evstigneev M, von Gehlen S and Reimann P 2009 Physical Review E 79 011116 47.Cheng L and Yip N K 2015 Physica D: Nonlinear Phenomena 297 1-32 48.Juniper M P N, Straube A V, Besseling R, Aarts D G A L and Dullens R P A 2015 Nature Communications 6 7187 49.Guérin T and Dean D S 2017 Physical Review E 95 012109 50.Bai Z W and Zhang W 2018 Chemical Physics 500 62-66 51.Bia las K and Spiechowicz J 2021 Chaos: An Interdisciplinary Journal of Nonlinear Science 31 123107 52.Zarlenga D G, Larrondo H A, Arizmendi C M and Family F 2007 Physical Review E 75 051101 53.Zarlenga D G, Larrondo H A, Arizmendi C M and Family F 2009 Physical Review E 80 011127 54.Zarlenga D G, Frontini G L, Family F and Arizmendi C M 2019 Physica A: Statistical Mechanics and its Applications 523 172-179 55.Mazzitello K I, Iguain J L, Jiang Y, Family F and Arizmendi C M 2019 Journal of Physics: Conference Series 1290 012022 56.Chien-Lih H 2005 The Mathematical Gazette 89 469-470 57.Howse J R, Jones R A L, Ryan A J, Gough T, Vafabakhsh R and Golestanian R 2007 Physical Review Letters 99 048102 58.ten Hagen B, van Teeffelen S and Löwen H 2011 Journal of Physics: Condensed Matter 23 194119 59.Sprenger A R, Jahanshahi S, Ivlev A V and Löwen H 2021 Physical Review E 103 042601 60.Mil'shtejn G N 1975 Theory of Probability & Its Applications 19 557-562 61.Malgaretti P, Puertas A M and Pagonabarraga I 2022 Journal of Colloid and Interface Science 608 2694-2702 62. Caprini L, Cecconi F and Marini Bettolo Marconi U 2021 The Journal of Chemical Physics 155 234902 63.Caprini L, Marconi U M B, Wittmann R and Löwen H 2022 Soft Matter 64.Cichocki B and Felderhof B U 1991 Physical Review A 44 6551-6558 65.Ackerson B J and Fleishman L 1982 The Journal of Chemical Physics 76 2675-2679 66. Löwen H 1992 Journal of Physics: Condensed Matter 4 10105-10116 67.Kollmann M 2003 Physical Review Letters 90 180602 68.Lutz C, Kollmann M and Bechinger C 2004 Physical Review Letters 93 026001 69.Lutz C, Kollmann M, Leiderer P and Bechinger C 2004 Journal of Physics: Condensed Matter 16 S4075-S4083 70.Juniper M P N, Zimmermann U, Straube A V, Besseling R, Aarts D G A L, Löwen H and Dullens R P A 2017 New Journal of Physics 19 013010 71.Berndt I, Pedersen J S and Richtering W 2005 Journal of the American Chemical Society 127 9372-9373 72.Frangipane G, Dell'Arciprete D, Petracchini S, Maggi C, Saglimbeni F, Bianchi S, Vizsnyiczai G, Bernardini M L and Di Leonardo R 2018 eLife 7 e36608 73.Arlt J, Martinez V A, Dawson A, Pilizota T and Poon W C K 2018 Nature Communications 9 768 74.Szamel G 2014 Physical Review E 90 012111 75.Wittmann R, Maggi C, Sharma A, Scacchi A, Brader J M and Marconi U M B 2017 Journal of Statistical Mechanics: Theory and Experiment 2017 113207 76.Caprini L and Marconi U M B 2018 Soft Matter 14 9044-9054 77.Dabelow L, Bo S and Eichhorn R 2019 Physical Review X 9 021009 78.Caprini L, Marconi U M B, Puglisi A and Vulpiani A 2019 Journal of Statistical Mechanics: Theory and Experiment 2019 053203 79.Lozano C, ten Hagen B, Löwen H and Bechinger C 2016 Nature Communications 7 12828 80. Lozano C, Liebchen B, ten Hagen B, Bechinger C and Löwen H 2019 Soft Matter 15 5185-5192