# Diffusion in periodic, correlated random forcing landscapes

David S Dean, Shamik Gupta, Gleb Oshanin, Alberto Rosso, Grégory Schehr

► To cite this version:

David S Dean, Shamik Gupta, Gleb Oshanin, Alberto Rosso, Grégory Schehr. Diffusion in periodic, correlated random forcing landscapes. Journal of Physics A: Mathematical and Theoretical, 2014, 47 (37), pp.372001. 10.1088/1751-8113/47/37/372001. hal-01069753

HAL Id: hal-01069753

https://hal.science/hal-01069753

Submitted on 23 Feb 2016

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>
# Diffusion in periodic, correlated random forcing landscapes

David S Dean<fnref n="1" />, Shamik Gupta<fnref n="2" />, Gleb Oshanin<fnref n="3" />,<fnref n="4" />,
Alberto Rosso<fnref n="2" /> and Grégory Schehr<fnref n="2" />

E mail: shamikg1@gmail.com

## Abstract

We study the dynamics of a Brownian particle in a strongly correlated quenched random potential defined as a periodically extended (with period $L$) finite trajectory of a fractional Brownian motion with arbitrary Hurst exponent $H \in (0, 1)$. While the periodicity ensures that the ultimate long time behavior is diffusive, the generalized Sinai potential considered here leads to a strong logarithmic confinement of particle trajectories at intermediate times. These two competing trends lead to dynamical frustration and result in a rich statistical behavior of the diffusion coefficient $D_L$: although one has the typical value $D_L^{typ} \sim \exp(-\beta L^H)$, we show via an exact analytical approach that the positive moments ($k > 0$) scale like $\langle D_L^k \rangle \sim \exp[-c'(k\beta L^H)^{1/(1+H)}]$, and the negative ones as $\langle D_L^{-k} \rangle \sim \exp(a'(k\beta L^H)^2)$, $c'$ and $a'$ being numerical constants and $\beta$ the inverse temperature. These results demonstrate that $D_L$ is strongly non-self-averaging. We further show that the probability distribution of $D_L$ has a log normal left tail and a highly singular, one sided log stable right tail reminiscent of a Lifshitz singularity.

Keywords: Brownian motion, quenched disorder, anomalous dynamics
PACS numbers: 05.40.–a, 02.50.–r, 05.10.Ln

<footnote n="1">¹ Université de Bordeaux and CNRS, Laboratoire d'Ondes et Matière d'Aquitaine (LOMA), UMR 5798, F 33400 Talence, France</footnote>

<footnote n="2">² Laboratoire de Physique Théorique et Modèles Statistiques (UMR CNRS 8626), Université de Paris Sud, Orsay Cedex, France</footnote>

<footnote n="3">³ Sorbonne Universités, UPMC Univ Paris 06, UMR 7600, LPTMC, F 75005 Paris, France</footnote>

<footnote n="4">⁴ CNRS, UMR 7600, Laboratoire de Physique Théorique de la Matière Condensée, F 75005 Paris, France</footnote>
<figure><img src="image_2.png" /><figcaption>Figure 1. Sketch of the potential V (x) for the ratchet (left) and the translationally invariant case (right). Here, $H = 1/3$.</figcaption></figure>

Transport in random media is extensively studied due to its practical and fundamental importance [1, 3]. In many cases, the dynamics is modelled as a Langevin process, with a drift generated by a quenched disordered potential. In theoretical analysis, the potential landscape is taken to be either infinitely extended or periodic in space. Stochastic dynamics in a periodic potential, both random and deterministic, is commonly encountered in many different contexts, including modulated structures [4], superionic conductors [5], colloids in light fields [6, 7], diffusion on regular [8, 11] and disordered [12, 16] solid surfaces, molecular motors on disordered tracks like DNA/RNA [17, 19], and motion in a tilted potential due to a random polymer [20].

Theoretical approaches often assume that the dynamics in a periodic potential reproduces the behavior in an infinitely extended potential. This is implemented by setting the period in the final result, e.g. for the velocity (if any) or the diffusion coefficient, to infinity [21, 22]. It is crucial to investigate how far such an assumption holds. Especially in the context of numerical simulations carried out for periodic systems, one may ask how reliably their results may be extrapolated to infinite systems.

In this work, we address these fundamental questions for a Langevin dynamics $x(t)$ in a periodic, quenched random potential $V(x(t))$ (with $V(x(t) + L) = V(x(t))$):

$$ \eta \frac{dx}{dt} = -\frac{dV(x)}{dx} + \xi(t), \qquad (1) $$

with $\eta$ the friction coefficient, $\xi(t)$ a Gaussian white noise with zero mean and correlations $\xi(t)\xi(t') = 2\eta T\delta(t - t')$, the overbar being an average over the noise, and the temperature $T$ is in units of the Boltzmann constant. We consider two cases:

* the ratchet case where $V(x)$ is a fractional Brownian motion (fBm) in time $x \in [0, L]$, with $V(0) = 0$ and $V(L)$ arbitrary. Thus, $V(x)$ is a Gaussian process with zero mean, $\langle V(x) \rangle = 0$, and variance

$$ \langle [V(x) - V(y)]^2 \rangle = \frac{V_0^2}{l^{2H}} |x - y|^{2H}; \quad x, y \in [0, L], \qquad (2) $$

where $H \in (0, 1)$ is the Hurst exponent, $V_0$ and $l$ define respectively the typical amplitude of $V(x)$ and its scale of variation over $x$. In (2) the angular brackets denote averaging with respect to $V(x)$. Figure 1 (left) shows a realization of $V(x)$, with a jump at $x=L$.

* the translationally invariant case, where $V(x)$ is a stationary Gaussian process, which at short length scales $|x - y| \ll L$ has the variance (2), and satisfies $V(0) = V(L)$, so that all points are statistically equivalent. The particle in this case diffuses on a ring. A realization of such a $V(x)$ is shown in figure 1 (right).

The dynamics (1) involves a combination of two paradigmatic situations: random motion in a periodic potential and the generalized Sinai dynamics in presence of a force $F(x) = -dV(x)/dx$ that is a time-independent stochastic variable with spatial correlations (except for $H = 1/2$ when $V(x)$ is the trajectory of a Brownian motion itself so that (1) is the periodic Sinai model [23]). While the latter produces an archetypal subdiffusion with
logarithmically-confined trajectories, the periodicity of the random potential enforces a long time diffusive behavior with a diffusion coefficient $D_L$. Here, we show that the trade off between these two competing trends results in a rich statistical behavior of $D_L$. In particular, $D_L$ is strongly non-self-averaging, with both negative and positive moments exhibiting an anomalous dependence on the temperature, period $L$ and the order of the moment, and being supported by atypical realizations of $V(x)$. For the ratchet case, we obtain exact analytical results, relying on exact bounds, for both positive and negative moments of $D_L$. We also discuss the full form of the probability distribution of $D_L$, and show that it is characterized by a log normal left tail and a highly singular log stable right tail, reminiscent of a Lifshitz singularity. We finally highlight the issue of sample to sample fluctuations of $D_L$. From standard scaling arguments and physical intuition, one expects that our exact results for the ratchet case also hold for the translationally invariant situation, which is harder to analyze analytically. This is confirmed below by thorough numerical simulations [24].

The dynamics (1) in an infinite system for arbitrary $H$, where $H > 1/2$ ($H < 1/2$) implies positively (negatively) correlated increments and superdiffusive (subdiffusive) $V(x)$, respectively, was discussed in [25] where it was shown that $\lim_{t \to \infty} \langle x^2(t) \rangle \sim \ln^{2/H}(t)$ (see also [26]). In contrast, in a periodic system, the long time motion is diffusive for any given realization of the potential $V(x)$, so that we have the diffusion coefficient $D_L \equiv \lim_{t \to \infty} x^2(t)/(2t)$, with $D_L$ given by [27, 31] (see also [7, 8, 12, 32])

$$ \frac{D_L}{D^{(0)}} = \left( \int_0^L \frac{dx}{L} \int_0^L \frac{dy}{L} e^{\beta[V(x)-V(y)]} \right)^{-1}, \qquad (3) $$

where $\beta$ is the inverse temperature, and $D^{(0)} = T/\eta$. Clearly, $D_L$ is a random variable that fluctuates between realizations of $V(x)$, and has support on $[0, D^{(0)}]$. The inverse of $D_L$ may be regarded as a product of partition functions in potentials $V(x)$ and $-V(x)$, respectively. The Brownian version of this quantity finds applications in disordered systems and has been extensively studied, while our results for $H \neq 1/2$ apply to more general situations (note that the marginal case $H=0$, when $V(X)$ is log-correlated, was studied in [33]). The expression in equation (3) also describes the ground state energy in a toy model of localization, and its average value was studied in [34] for $H=1/2$.

Turning to the discussion of the behavior of $D_L$, we first reduce the number of parameters. In the following, we measure $L$ in units of $l$ [see equation (2)], absorb $V_0$ into $\beta$, and measure $D_L$ in units of $D^{(0)}$, so that $D_L$ has support on $[0, 1]^5$. Now, the typical behavior of $D_L$ is easy to estimate as $D_L^{\text{typ}} \propto L^2/\tau_{\text{typ}}$, where the dimensionless $L$ sets the scale of an inhomogeneous region, and $\tau_{\text{typ}}$ defines the typical (dimensionless) time a particle spends in this region. The random potential $V(x)$ being a fBm, the typical height of the potential barrier over a length $L$ scales like $L^H$, for both the ratchet and the translationally invariant case. Assuming Arrhenius type activation, one expects $\tau_{\text{typ}} \sim e^{\beta L^H}$, which implies

$$ D_L^{\text{typ}} \sim L^2 \exp(-\beta L^H). \qquad (4) $$

The average behavior of $D_L$ is a much more delicate question because, as can be seen from equation (3), computing the statistics of $D_L$ is a highly non-trivial task that involves the study of an exponential functional of fBm, for which standard methods like the Feynman Kac formula are of little use for $H \neq 1/2$.

Let us first summarize our main analytical results obtained for the ratchet case: We find that the average of the logarithm of $D_L$ is given, to leading order in $L$, by

<footnote n="5">⁵ $l$, $V_0$ and $D^{(0)}$ may be made explicit in the final result by the replacements $L \to L/l$, $\beta \to \beta V_0$ and $D \to D/D^{(0)}$.</footnote>
$$ \langle \ln D_L \rangle \propto -2 m \beta L^H, \qquad (5) $$

with $m = \langle \max_{s \in [0,1]} V(s) \rangle$. The result in equation (5) is consistent with the logarithmic growth of the disorder averaged mean square displacement in an infinite system: $\lim_{t \to \infty} \langle x^2(t) \rangle \sim \ln^{2/H}(t)$, and the argument leading to equation (4). Further on, we obtain sharp bounds for the positive moments ($k > 0$) of the random variable $D_L$:

$$ A_k(L) \le \langle D_L^k \rangle \le B_k(L), \qquad (6) $$

where, in the limit $L \to \infty$, the bounds satisfy

$$ \begin{aligned} A_k(L) &= \exp \left[ -(1+H)(2k\beta)^{\frac{1}{1+H}} \left(\frac{C}{H}L\right)^{\frac{H}{1+H}} \right], \\ B_k(L) &= \exp \left[ -(1+H)(k\beta)^{\frac{1}{1+H}} \left(\frac{c}{H}L\right)^{\frac{H}{1+H}} \right], \end{aligned} \qquad (7) $$

with constants $c$ and $C$, $0 < c \le C < \infty$, being independent of $L$, $k$ and $\beta$. Both bounds exhibit the same dependence on $L$, $k$, $\beta$, from which we infer that the exact asymptotic result has the same functional form. Finally, for the negative moments of $D_L$, we find

$$ \langle D_L^{-k} \rangle \sim \exp \left( \frac{a^2 k^2 \beta^2 L^{2H}}{4} \right), \qquad (8) $$

where $a$ is a constant (independent of $L$, $k$ and $\beta$). Our exact results, equations (6) and (8), thus show that the positive and negative moments are dominated by atypical realizations of $V(x)$, in contrast to $\langle \ln D_L \rangle$, see equation (5).

We now turn to a derivation of our results. Using equation (3), the logarithm of $D_L$ can be formally written as $\ln D_L = \ln J_+(L) + \ln J_-(L) + 2 \ln L$, where $J_\pm(L)$ are stationary currents through a finite sample of length $L$ with potentials $\pm V(x)$ [26]: $J_\pm(L) = [\int_0^L dx \exp[\pm \beta V(x)]]^{-1}$. Statistical properties of these currents for the Sinai problem ($H = 1/2$) are known [35, 40]. Using the results of [41], and noting that with $\langle V(x) \rangle = 0$, $J_+(L)$ and $J_-(L)$ have equal moments, we have for arbitrary $H$ and to leading order in $L$, $\langle \ln J_+(L) \rangle = \langle \ln J_-(L) \rangle \propto -m\beta L^H$, which yields equation (5).

The proof of the result in equation (6) is based on a Theorem due to Monrad and Rootzén [42] on the probability that a fBm $V(x)$, with $V(0) = 0$, remains within a strip of width $\epsilon$ for the time $x \in [0, L]$. Defining $M_L = \max_{0 \le x \le L} |V(x)|$, the Monrad Rootzén theorem, in our notation, states that $P(M_L \le \epsilon)$ satisfies

$$ \exp \left( -CL \epsilon^{-\frac{1}{H}} \right) \le P(M_L \le \epsilon) \le \exp \left( -cL \epsilon^{-\frac{1}{H}} \right), \qquad (9) $$

for $0 < \epsilon \le L^H$, $c$ and $C$ being $L$-independent constants [see equation (6)].

Consider the lower bound in equation (6). Suppose we average the positive definite quantity $D_L^k$ by considering instead of the entire set $\Omega$ of all possible paths $V(x)$ only a subset $\Omega' \subset \Omega$ of paths such that $M_L \le \epsilon$. This gives the lower bound $\langle D_L^k \rangle_\Omega \ge \langle D_L^k \rangle_{\Omega'} P(M_L \le \epsilon)$. However, for paths in $\Omega'$, we have $\exp(\beta[V(x) - V(y)]) \le \exp(2\beta\epsilon)$, and hence, $(\int_0^L \frac{dx}{L} \int_0^L \frac{dy}{L} e^{\beta[V(x) - V(y)]})^{-1} \ge e^{-2\beta\epsilon}$. Therefore, we obtain $\langle D_L^k \rangle_\Omega \ge e^{-2k\beta\epsilon} P(M_L \le \epsilon)$. Making the inequality more stringent by choosing the lower bound in equation (9), we get
$$ \langle D_L^k \rangle_\Omega \ge e^{-2k\beta\epsilon} \exp(-CL\epsilon^{-\frac{1}{H}}), \qquad (10) $$

which holds for any $\epsilon$ with $0 < \epsilon \le L^H$.

The function on the right-hand side (rhs) of equation (10) is a non-monotonic function of $\epsilon$, attaining its maximum at $\epsilon = \epsilon_{opt} = (CL/k2\beta H)^{H/(1+H)}$. Clearly, the best lower bound corresponds to the choice $\epsilon = \epsilon_{opt}$, leading to the lower bound in equation (6). Note that to satisfy the conditions of validity of the Monrad Rootzén theorem, we require that $\epsilon_{opt} \le L^H$, that is, $2k\beta L \ge C/H$, which is easily realized for sufficiently large $L$. The derivation of the lower bound is an example of the Lifshitz optimal fluctuation method [43], which has been used to bound the survival probability of particles diffusing in the presence of randomly scattered immobile traps (see, e.g., [44]).

We now discuss the derivation of the upper bound. To this end, we discretize $x$, and write the rhs of equation (3) as

$$ \int_0^L \frac{dx}{L} \int_0^L \frac{dy}{L} e^{\beta[V(x)-V(y)]} \sim \sum_{i,j=1}^N e^{\beta[V(j)-V(i)]}. \qquad (11) $$

Given that the fBm starts at $V(0) = 0$, at least one term in the double sum on the rhs takes the value $\exp(\beta M_L)$, corresponding to the point $x = 0$ and the point where $|V(x)|$ attains its maximal value $M_L$. Now, as all the other terms are positive, we have the bound $D_L^{-1} \ge \exp(\beta M_L)$, and, thus,

$$ \langle D_L^k \rangle \le \langle e^{-k\beta M_L} \rangle = k\beta \int_0^\infty d\epsilon P(M_L \le \epsilon) e^{-k\beta\epsilon}. \qquad (12) $$

The integral in the rhs is dominated, for large $L$, by the small $\epsilon$ region, where we can thus use the upper bound in (9). Performing the remaining integral over $\epsilon$ by the saddle-point method and omitting the pre-exponential terms lead to the upper bound in equation (6). The result in (6) has several striking features. Namely, the function

$$ \mu(k, \beta, L) = -\ln \langle D_L^k \rangle, \qquad (13) $$

(as compared to its typical counterpart given by equation (4) as $\mu_{typ}(k, \beta, L) \sim k\beta L^H$), (a) grows sub-linearly with $k$ (multifractality), (b) is a non-analytic, sublinear function of $\beta$, which implies a rather unusual sub-Arrhenius dependence of the positive moments on the temperature, and (c) exhibits a slower anomalous growth with $L$ as $\sim L^{\frac{H}{1+H}}$. This means that the disorder-averaged $D_L$ is generically larger than the one expected on the basis of typical realizations of the disorder. In turn, this implies that the behavior of the positive moments of $D_L$ is supported by atypical realizations of disorder, reminiscent of the so-called Lifshitz singularities, as discussed above. In conclusion, one cannot infer the dynamical behavior in an infinite system from the positive moments of $D_L$. This is surprising at first glance, as $\langle x^2(t) \rangle$ is linearly proportional to $\langle D_L \rangle$, and shows that the limits $t \to \infty$ and $L \to \infty$ do not commute in this system.

The behavior of the negative moments $\langle D_L^{-k} \rangle$ with $k = 1, 2, ...$ is determined by essentially the same approach as above. Note that both the lower and the upper bound on $D_L^{-k}$ are made tighter for a given realization of $V(x)$ by using $D_L^{-k} \sim \exp(k\beta S)$, where $S$ is the span of $V(x)$ (the difference between the maximum and minimum) on the interval $[0, L]$. Therefore, in contrast to the positive moments, the negative moments are supported by realizations of $V(x)$ with a large span. Using the result that for large $S$, $P(M_L = S) \sim \exp(-S^2/aL^2^H)$, with $a$ a constant, integration of equation (12) gives the result announced in equation (8), which displays a super-Arrhenius dependence on the temperature,
<figure><img src="image_3.png" /><figcaption>**Figure 2.** (a) $-\ln[\langle D_L \rangle/L^2]$ versus $L$ for three values of $H$ corresponding to diffusive, subdiffusive, and superdiffusive $V(x)$. The inset shows $-\langle \ln[D_L/L^2] \rangle$ as a function of $L$ for three values of $H$, see equation (5). (b) $\mu(k, \beta, L)$ in equation (13) as a function of $k$ for $L=512$ and three values of $H$. In all cases, the symbols denote simulation results for the translationally invariant $V(x)$, while the slopes of the solid lines correspond to the results derived for the ratchet case.</figcaption></figure>

a superlinear dependence on $k$, and a strong dependence on $L$. A similar result was obtained earlier in [45]. We note that, as for the positive moments, one cannot deduce the behavior in an infinite system from that of negative moments of $D_L$ in a periodic system, as the latter is supported by atypical realizations of $V(x)$ that have anomalously large span scaling as $S \sim L^{2H}$, while the typical behavior is $S_{typ} \sim L^H$.

Based on our results for the moments, we now obtain the probability distribution $P(D_L)$. As already explained, the behavior of the negative moments is supported by anomalously stretched trajectories of $V(x)$ for which the value of $D_L$ is small. One may thus expect in view of the form of the moments in equation (8) that for small $D_L$, $P(D_L)$ is log-normal:

$$P(D_L) \sim \frac{1}{a\beta L^H D_L} \exp\left(-\frac{\ln^2(D_L)}{a^2\beta^2 L^{2H}}\right). \quad (14)$$

To analyze the behavior of $P(D_L)$ for $D_L$ close to 1, we recall the formal definition of the one-sided Lévy distribution $\mathcal{L}_\nu(z)$, $0 \le z < \infty$, of order $\nu$ (see, e.g., [46]):

$$\int_0^\infty dz e^{-pz} \mathcal{L}_\nu(z) = e^{-p\nu}. \quad (15)$$

The asymptotic behavior of $\mathcal{L}_\nu(z)$ is well-known [46], and, in particular, one has $\mathcal{L}_\nu(z) \sim z^{-\sigma} \exp(-b/z^\tau)$ for $z \to 0$, where $b$ is a computable constant, $\sigma = (2 - \nu)/(2(1 - \nu))$, and $\tau = \nu/(1 - \nu)$. It is important to note that this precise asymptotic form is responsible for the stretched-exponential behavior in equation (15), which is immediately verified by substituting the form in equation (15), and performing the integration by the saddle-point method. Moreover, one realizes by making in equation (15) a change of the integration variable $z = \ln(1/D_L)/\beta L^H$, choosing $\nu = 1/(1 + H)$, and setting $p = k/\beta L^H$ that equation (15) becomes identical to the result in equation (6), up to numerical factors. It follows that for $D_L$ close to 1 (i.e., $z$ close to 0), the distribution function behaves as

$$P(D_L) \sim \frac{1}{\beta L^H D_L} \mathcal{L}_{1+H} \left[ \ln (D_L^{-1}) / \beta L^H \right]. \quad (16)$$

Using the asymptotic $\mathcal{L}_\nu$ given above, we get that $P(D_L)$ is highly singular near the right edge, $P(D_L) \sim \exp[-b\beta^{1/H}L/(1 - D_L)^{1/H}]$, similar to the Lifshitz singularity.
<figure><img src="image_4.png" /><figcaption>**Figure 3.** (a) Distribution $P(\ln D_L)$ for $L=1024$ and $H = 1/3$. The red line denotes numerical results, while analytical predictions for the right and left tails behaving as $\exp(-a(H)/\ln^{1/H}D_L^{-1})$ and $\exp(-b(H)\ln^2(D_L))$, respectively, with $a(H)$ and $b(H)$ being constants, are shown by black lines. (b) Numerical results for $P(D)$ for different $L$ and $H = 1/3$.</figcaption></figure>

We now consider the translationally invariant case. Here, we expect our above analysis, in particular, result (6) to hold, up to possible numerical factors. To demonstrate this, we now present results of extensive numerical simulations: figure 2 for $-\ln((D_L)/L^2)$, $-\langle\ln(D_L/L^2)\rangle$, and $\mu(k, \beta, L)$, and figure 3(a) for $P(\ln D_L)$ indeed show a very good agreement that supports our expectations.

To close, we ask: if we have two different realizations of $V(x)$, and correspondingly, two different values, $D_L$ and $D'_L$, of the diffusion coefficient, how likely are these values equal? We introduce a random variable $D \equiv \frac{D_L}{D_L+D'_L}$, $D \in [0, 1]$, and analyze its distribution $P(D)$ via numerical simulations. Clearly, $D = 1/2$ maximizing $P(D)$ implies that the two values of $D_L$ are most likely very close to one another. Variables such as $D$ play a key role in various scale-independent hypothesis testing procedures, in classical problems in statistics, in signal processing (see, e.g., [47]), and in the analysis of chaotic scattering in few-channel systems [48]. Such variables are used to characterize the effective width of narrow distributions possessing moments of arbitrary order [49, 50].

In figure 3(b), we present numerical results for $P(D)$ for different values of $L$ and $H = 1/3$, for the translationally invariant case. We observe an interesting phenomenon of a change in the form of the distribution as $L$ is increased. For relatively small $L$, the distribution is bell-shaped and centered at $D = 1/2$. However, on increasing $L$, $P(D)$ broadens, becomes almost flat at a certain critical $L$ (whose value depends on $H$), and then changes its shape so that $D = 1/2$ minimizes the distribution. This implies that for sufficiently large $L$, two values of the diffusion coefficients obtained for two different realizations of $V(x)$ are most likely very different, and the event $D_L = D'_L$ is the least probable. Note that a similar dependence in the distribution of the prefactor in the Sinai law on the strength of disorder was recently observed in [51].

## Acknowledgments

GS is supported by the ANR grant 2011-BS04 013-01 WALKMAT. This project is partially supported by the Indo-French Centre for the Promotion of Advanced Research under Project 4604 3. SG, AR and GS thank the Galileo Galilei Institute for Theoretical Physics, Florence, Italy for the hospitality and the INFN for partial support during the completion of this work.
References

[1] Bouchaud J P, Comtet A, Georges A and le Doussal P 1990 *Ann. Phys.* **201** 285

[2] Bouchaud J P and Georges A 1990 *Phys. Rep.* **195** 127

[3] Oshanin G, Burlatsky S F, Moreau M and Gaveau B 1993 *Chem. Phys.* **177** 803

[4] Schneider T, Politi A and Sörensen M P 1988 *Phys. Rev. A* **37** 948

[5] Dietrich W and Peschel I 1977 *Z. Phys. B* **27** 177

[6] Evers F et al 2013 *Eur. Phys. J. Spec. Top.* **222** 2995

[7] Dean D S and Touya C 2008 *J. Phys. A: Math. Theor.* **41** 335002

[8] Reimann P, van den Broeck C, Linke H, Hänggi P, Rubi J M and Pérez Madrid A 2001 *Phys. Rev. Lett.* **87** 010602

Reimann P, van den Broeck C, Linke H, Hänggi P, Rubi J M and Pérez Madrid A 2002 *Phys. Rev. E* **65** 031104

[9] Sancho J M, Lacasta A M, Lindenberg K, Sokolov I M and Romero A H 2004 *Phys. Rev. Lett.* **92** 250601

[10] Lindenberg K, Sancho J M, Lacasta A M and Sokolov I M 2007 *Phys. Rev. Lett.* **98** 020602

[11] Lindner B, Kostur M and Schimansky Geier L 2001 *Fluct. Noise Lett.* **1** R25

[12] Lindenberg K, Sancho J M, Khoury M and Lacasta A M 2012 *Fluct. Noise Lett.* **11** 1240004

[13] Khoury M, Gleeson J P, Sancho J M, Lacasta A M and Lindenberg K 2009 *Phys. Rev. E* **80** 021123

[14] Reimann P and Eichhorn R 2008 *Phys. Rev. Lett.* **101** 180601

[15] Khoury M, Lacasta A M, Sancho J M and Lindenberg K 2011 *Phys. Rev. Lett.* **106** 090602

[16] Simon M S, Sancho J M and Lindenberg K 2013 *Phys. Rev. E* **88** 062105

[17] Kafri Y and Nelson D R 2005 *J. Phys.: Condens. Matter* **17** S3871

[18] Kafri Y, Lubensky D K and Nelson D R 2005 *Phys. Rev. E* **71** 041906

[19] Kafri Y, Lubensky D K and Nelson D R 2004 *Biophys. J.* **86** 3373 91

[20] Salgado Garcia R and Maldonado C 2013 *Phys. Rev. E* **88** 062143
Salgado Garcia R arXiv:1404.2852

[21] Derrida B and 1983 *J. Stat. Phys.* **31** 433

[22] Dean D S, Drummond I T and Horgan R R 1997 *J. Phys. A: Math. Gen.* **30** 385

[23] Sinai Ya G 1982 *Theor. Probab. Appl.* **27** 256

[24] Santachiara R, Rosso A and Krauth W 2007 *J. Stat. Mech.* P02009

[25] Marinari E, Parisi G, Ruelle D and Windey P 1983 *Phys. Rev. Lett.* **50** 1223

[26] Oshanin G, Rosso A and Schehr G 2013 *Phys. Rev. Lett.* **110** 100602

[27] Lifson S and Jackson J L 1962 *J. Chem. Phys.* **36** 2410

[28] Festa R and Galleani d'Agliano E 1978 *Physica A* **90** 229

[29] Ferrari P A, Goldstein S and Lebowitz J L 1985 *Statistical Physics and Dynamical Systems* ed J Fritz, A Jaffe and D Szasz (Boston, MA: Birkhauser)

[30] Golden K, Goldstein S and Lebowitz J L 1985 *Phys. Rev. Lett.* **55** 2629

[31] Zwanzig R 1988 *Proc. Natl. Acad. Sci. USA* **85** 2029

[32] Baiesi M, Maes C and Wynants B 2011 *Proc. R. Soc. Lond. A* **467** 2792

[33] Castillo H E and Le Doussal P 2001 *Phys. Rev. Lett.* **86** 4859

[34] Monthus C, Oshanin G, Comtet A and Burlatsky S F 1996 *Phys. Rev. E* **54** 231

[35] Burlatsky S F, Oshanin G, Mogutov A and Moreau M 1992 *Phys. Rev. A* **45** R6955

[36] Oshanin G, Mogutov A and Moreau M 1993 *J. Stat. Phys.* **73** 379

[37] Monthus C and Comtet A 1994 *J. Phys. I France* **4** 635

[38] Comtet A, Monthus C and Yor M 1998 *J. Appl. Probab.* **35** 255

[39] Oshanin G and Redner S 2009 *Eur. Phys. Lett.* **85** 10008

[40] Oshanin G and Schehr G 2012 *Quant. Finance* **12** 1325

[41] Molchan G M 1999 *Commun. Math. Phys.* **205** 97

[42] Monrad D and Rootzén H 1995 *Probab. Theory Relat. Fields* **101** 173
Li W V, Linde W and Acad C R 1998 *Sci. Paris* **326** 1329

[43] Lifshitz I M 1963 *Sov. Phys. JETP* **17** 1159
Lifshitz I M 1965 *Sov. Phys. Usp.* **7** 549

[44] Grassberger P and Procaccia I 1983 *J. Chem. Phys.* **77** 6281

[45] Comtet A and Dean D S 1998 *J. Phys. A: Math. Gen.* **31** 8595

[46] Schehr G and le Doussal P 2010 *J. Stat. Mech.* P01009

[47] Dharmawansa P, McKay M R and Cheng Y 2013 *SIAM J. Matrix Anal. Appl.* **34** 257

[48] Mejia Monasterio C, Oshanin G and Schehr G 2011 *Phys. Rev. E* **84** 035203

[49] Mejia Monasterio C, Oshanin G and Schehr G 2011 *J. Stat. Mech.* P06022

[50] Mattos T G, Mejia Monasterio C, Metzler R and Oshanin G 2012 *Phys. Rev. E* **86** 031143

[51] Boyer D, Dean D S, Mejia Monasterio C and Oshanin G 2012 *Phys. Rev. E* **85** 031136