Source author record

Pierre Le Doussal

Pierre Le Doussal appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

99works
17topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

99 published item(s)

preprint2022arXiv

Clusters in an epidemic model with long-range dispersal

In presence of long range dispersal, epidemics spread in spatially disconnected regions known as clusters. Here, we characterize exactly their statistical properties in a solvable model, in both the supercritical (outbreak) and critical regimes. We identify two diverging length scales, corresponding to the bulk and the outskirt of the epidemic. We reveal a nontrivial critical exponent that governs the cluster number, the distribution of their sizes and of the distances between them. We also discuss applications to depinning avalanches with long range elasticity.

preprint2022arXiv

Counting statistics for non-interacting fermions in a rotating trap

We study the ground state of $N \gg 1$ noninteracting fermions in a two-dimensional harmonic trap rotating at angular frequency $Ω>0$. The support of the density of the Fermi gas is a disk of radius $R_e$. We calculate the variance of the number of fermions ${\cal N}_R$ inside a disk of radius $R$ centered at the origin for $R$ in the bulk of the Fermi gas. We find rich and interesting behaviours in two different scaling regimes: (i) $Ω/ ω<1 $ and (ii) $1 - Ω/ ω= O(1/N)$, where $ω$ is the angular frequency of the oscillator. In the first regime (i) we find that ${\rm Var}\,{\cal N}_{R}\simeq\left(A\log N+B\right)\sqrt{N}$ and we calculate $A$ and $B$ as functions of $R/R_e$, $Ω$ and $ω$. We also predict the higher cumulants of ${\cal N}_{R}$ and the bipartite entanglement entropy of the disk with the rest of the system. In the second regime (ii), the mean fermion density exhibits a staircase form, with discrete plateaus corresponding to filling $k$ successive Landau levels, as found in previous studies. Here, we show that ${\rm Var}\,{\cal N}_{R}$ is a discontinuous piecewise linear function of $\sim (R/R_e) \sqrt{N}$ within each plateau, with coefficients that we calculate exactly, and with steps whose precise shape we obtain for any $k$. We argue that a similar piecewise linear behavior extends to all the cumulants of ${\cal N}_{R}$ and to the entanglement entropy. We show that these results match smoothly at large $k$ with the above results for $Ω/ω=O(1)$. These findings are nicely confirmed by numerical simulations. Finally, we uncover a universal behavior of ${\rm Var}\,{\cal N}_{R}$ near the fermionic edge. We extend our results to a three-dimensional geometry, where an additional confining potential is applied in the $z$ direction.

preprint2022arXiv

Half-space stationary Kardar-Parisi-Zhang equation beyond the Brownian case

We study the Kardar-Parisi-Zhang equation on the half-line $x \geqslant 0$ with Neumann type boundary condition. Stationary measures of the KPZ dynamics were characterized in recent work: they depend on two parameters, the boundary parameter $u$ of the dynamics, and the drift $-v$ of the initial condition at infinity. We consider the fluctuations of the height field when the initial condition is given by one of these stationary processes. At large time $t$, it is natural to rescale parameters as $(u,v)=t^{-1/3}(a,b)$ to study the critical region. In the special case $a+b=0$, treated in previous works, the stationary process is simply Brownian. However, these Brownian stationary measures are particularly relevant in the bound phase ($a<0$) but not in the unbound phase. For instance, starting from the flat or droplet initial data, the height field near the boundary converges to the stationary process with $a>0$ and $b=0$, which is not Brownian. For $a+b\geqslant 0$, we determine exactly the large time distribution $F_{a,b}^{\rm stat}$ of the height function $h(0,t)$. As an application, we obtain the exact covariance of the height field in a half-line at two times $1\ll t_1\ll t_2$ starting from stationary initial data, as well as estimates, when starting from droplet initial data, in the limit $t_1/t_2\to 1$.

preprint2022arXiv

More on the Brownian force model: avalanche shapes, tip driven, higher $d$

The Brownian force model (BFM) is the mean-field model for the avalanches of an elastic interface slowly driven in a random medium. It describes the spatio-temporal statistics of the velocity field, and, to some extent is analytically tractable. We extend our previous studies to obtain several observables for the BFM with short range elasticity, related to the local jump sizes $S(x)$ and to the avalanche spatial extension in $d=1$, or the avalanche span in $d>1$. In $d=1$ we consider both driving (i) by an imposed force (ii) by an imposed displacement "at the tip" and obtain in each case the mean spatial shape $\langle S(x) \rangle$ at fixed extension, or at fixed seed to edge distance. We find that near an edge $x_e$, $S(x) \simeq σ|x-x_e|^3$ where $σ$ has a universal distribution that we obtain. We also obtain the spatiotemporal shape near the edge. In $d>1$ we obtain (i) the mean shape $\langle S(x_1,x_\perp) \rangle$ for a fixed span, which exhibits a non-trivial dependence in the transverse distance to the seed $x_\perp$ (ii) the mean shape around a point which has not moved, $\langle S(x) \rangle_{S(0)=0}$, which vanishes at the center as $|x|^{b_d}$ with non trivial exponents, $b_1=4$, $b_2=2 \sqrt{2}$ and $b_3=\frac{1}{2} (\sqrt{17}-1)$. We obtain the probability distributions in any $d$ for the maximal radius of an avalanche and the minimal distance of approach to a given point, as well as the probability of not hitting a cone in $d=2$. These results equivalently apply to the continuum limit of some branching diffusions, as detailed in a companion paper.

preprint2022arXiv

Quench dynamics of noninteracting fermions with a delta impurity

We study the out-of-equilibrium dynamics of noninteracting fermions in one dimension and in continuum space, in the presence of a delta impurity potential at the origin whose strength $g$ is varied at time $t=0$. The system is prepared in its ground state with $g=g_0=+\infty$, with two different densities and Fermi wave-vectors $k_L$ and $k_R$ on the two half-spaces $x>0$ and $x<0$ respectively. It then evolves for $t>0$ as an isolated system, with a finite impurity strength $g$. We compute exactly the time dependent density and current. For a fixed position $x$ and in the large time limit $t \to \infty$, the system reaches a non-equilibrium stationary state (NESS). We obtain analytically the correlation kernel, density, particle current, and energy current in the NESS, and characterize their relaxation, which is algebraic in time. In particular, in the NESS, we show that, away from the impurity, the particle density displays oscillations which are the non-equilibrium analog of the Friedel oscillations. In the regime of "rays", $x/t=ξ$ fixed with $x, t \to \infty$, we compute the same quantities and observe the emergence of two light cones, associated to the Fermi velocities $k_L$ and $k_R$ in the initial state. Interestingly, we find non trivial quantum correlations between two opposite rays with velocities $ξ$ and $-ξ$ which we compute explicitly. We extend to a continuum setting and to a correlated initial state the analytical methods developed in a recent work of Ljubotina, Sotiriadis and Prosen, in the context of a discrete fermionic chain with an impurity. We also generalize our results to an initial state at finite temperature, recovering, via explicit calculations, some predictions of conformal field theory in the low energy limit.

preprint2021arXiv

Hole probability for noninteracting fermions in a $d$-dimensional trap

The hole probability, i.e., the probability that a region is void of particles, is a benchmark of correlations in many body systems. We compute analytically this probability $P(R)$ for a spherical region of radius $R$ in the case of $N$ noninteracting fermions in their ground state in a $d$-dimensional trapping potential. Using a connection to the Laguerre-Wishart ensembles of random matrices, we show that, for large $N$ and in the bulk of the Fermi gas, $P(R)$ is described by a universal scaling function of $k_F R$, for which we obtain an exact formula ($k_F$ being the local Fermi wave-vector). It exhibits a super exponential tail $P(R)\propto e^{- κ_d (k_F R)^{d+1}}$ where $κ_d$ is a universal amplitude, in good agreement with existing numerical simulations. When $R$ is of the order of the radius of the Fermi gas, the hole probability is described by a large deviation form which is not universal and which we compute exactly for the harmonic potential. Similar results also hold in momentum space.

preprint2021arXiv

Inverse scattering solution of the weak noise theory of the Kardar-Parisi-Zhang equation with flat and Brownian initial conditions

We present the solution of the weak noise theory (WNT) for the Kardar-Parisi-Zhang equation in one dimension at short time for flat initial condition (IC). The non-linear hydrodynamic equations of the WNT are solved analytically through a connexion to the Zakharov-Shabat (ZS) system using its classical integrability. This approach is based on a recently developed Fredholm determinant framework previously applied to the droplet IC. The flat IC provides the case for a non-vanishing boundary condition of the ZS system and yields a richer solitonic structure comprising the appearance of multiple branches of the Lambert function. As a byproduct, we obtain the explicit solution of the WNT for the Brownian IC, which undergoes a dynamical phase transition. We elucidate its mechanism by showing that the related spontaneous breaking of the spatial symmetry arises from the interplay between two solitons with different rapidities.

preprint2021arXiv

Steady state of the KPZ equation on an interval and Liouville quantum mechanics

We obtain a simple formula for the stationary measure of the height field evolving according to the Kardar-Parisi-Zhang equation on the interval $[0,L]$ with general Neumann type boundary conditions and any interval size. This is achieved using the recent results of Corwin and Knizel (arXiv:2103.12253) together with Liouville quantum mechanics. Our formula allows to easily determine the stationary measure in various limits: KPZ fixed point on an interval, half-line KPZ equation, KPZ fixed point on a half-line, as well as the Edwards-Wilkinson equation on an interval.

preprint2020arXiv

Kernels for noninteracting fermions via a Green's function approach with applications to step potentials

The quantum correlations of $N$ noninteracting spinless fermions in their ground state can be expressed in terms of a two-point function called the kernel. Here we develop a general and compact method for computing the kernel in a general trapping potential in terms of the Green's function for the corresponding single particle Schrödinger equation. For smooth potentials the method allows a simple alternative derivation of the local density approximation for the density and of the sine kernel in the bulk part of the trap in the large $N$ limit. It also recovers the density and the kernel of the so-called {\em Airy gas} at the edge. This method allows to analyse the quantum correlations in the ground state when the potential has a singular part with a fast variation in space. For the square step barrier of height $V_0$, we derive explicit expressions for the density and for the kernel. For large Fermi energy $μ>V_0$ it describes the interpolation between two regions of different densities in a Fermi gas, each described by a different sine kernel. Of particular interest is the {\em critical point} of the square well potential when $μ=V_0$. In this critical case, while there is a macroscopic number of fermions in the lower part of the step potential, there is only a finite $O(1)$ number of fermions on the shoulder, and moreover this number is independent of $μ$. In particular, the density exhibits an algebraic decay $\sim 1/x^2$, where $x$ is the distance from the jump. Furthermore, we show that the critical behaviour around $μ= V_0$ exhibits universality with respect with the shape of the barrier. This is established (i) by an exact solution for a smooth barrier (the Woods-Saxon potential) and (ii) by establishing a general relation between the large distance behavior of the kernel and the scattering amplitudes of the single-particle wave-function.

preprint2020arXiv

Large deviations for the KPZ equation from the KP equation

Recently, Quastel and Remenik \cite{QRKP} [arXiv:1908.10353] found a remarkable relation between some solutions of the finite time Kardar-Parisi-Zhang (KPZ) equation and the Kadomtsev-Petviashvili (KP) equation. Using this relation we obtain the large deviations at large time and at short time for the KPZ equation with droplet initial conditions, and at short time with half-Brownian initial conditions. It is consistent with previous results and allows to obtain sub-leading corrections, as well as results at intermediate time. In addition, we find that the appropriate generating function associated to the full Brownian initial condition also satisfies the KP equation. Finally, generating functions for some linear statistics of the Airy point process are also found to satisfy the KP property, and consequences are discussed.

preprint2020arXiv

Manifolds in high dimensional random landscape: complexity of stationary points and depinning

We obtain explicit expressions for the annealed complexities associated respectively with the total number of (i) stationary points and (ii) local minima of the energy landscape for an elastic manifold with internal dimension $d<4$ embedded in a random medium of dimension $N \gg 1$ and confined by a parabolic potential with the curvature parameter $μ$. These complexities are found to both vanish at the critical value $μ_c$ identified as the Larkin mass. For $μ<μ_c$ the system is in complex phase corresponding to the replica symmetry breaking in its $T=0$ thermodynamics. The complexities vanish respectively quadratically (stationary points) and cubically (minima) at $μ_c^-$. For $d\geq 1$ they admit a finite "massless" limit $μ=0$ which is used to provide an upper bound for the depinning threshold under an applied force.

preprint2020arXiv

Noninteracting trapped Fermions in double-well potentials: inverted parabola kernel

We study a system of $N$ noninteracting spinless fermions in a confining, double-well potential in one dimension. When the Fermi energy is close to the value of the potential at its local maximum we show that physical properties, such as the average density and the fermion position correlation functions, display a universal behavior that depends only on the local properties of the potential near its maximum. This behavior describes the merging of two Fermi gases, which are disjoint at sufficiently low Fermi energies. We describe this behavior in terms of a new correlation kernel that we compute analytically and we call it the inverted parabola kernel". As an application, we calculate the mean and variance of the number of particles in an interval of size $2L$ centered around the position of the local maximum, for sufficiently small $L$. Finally, we discuss the possibility of observing our results in experiments, as well as the extensions to nonzero temperature and to higher space dimensions.

preprint2020arXiv

Replica Bethe Ansatz solution to the Kardar-Parisi-Zhang equation on the half-line

We consider the Kardar-Parisi-Zhang (KPZ) equation for the stochastic growth of an interface of height $h(x,t)$ on the positive half line with boundary condition $\partial_x h(x,t)|_{x=0}=A$. It is equivalent to a continuum directed polymer (DP) in a random potential in half-space with a wall at $x=0$ either repulsive $A>0$, or attractive $A<0$. We provide an exact solution, using replica Bethe ansatz methods, to two problems which were recently proved to be equivalent [Parekh, arXiv:1901.09449]: the droplet initial condition for arbitrary $A \geqslant -1/2$, and the Brownian initial condition with a drift for $A=+\infty$ (infinite hard wall). We study the height at $x=0$ and obtain (i) at all time the Laplace transform of the distribution of its exponential (ii) at infinite time, its exact probability distribution function (PDF). These are expressed in two equivalent forms, either as a Fredholm Pfaffian with a matrix valued kernel, or as a Fredholm determinant with a scalar kernel. For droplet initial conditions and $A> - \frac{1}{2}$ the large time PDF is the GSE Tracy-Widom distribution. For $A= \frac{1}{2}$, the critical point at which the DP binds to the wall, we obtain the GOE Tracy-Widom distribution. In the critical region, $A+\frac{1}{2} = εt^{-1/3} \to 0$ with fixed $ε= \mathcal{O}(1)$, we obtain a transition kernel continuously depending on $ε$. Our work extends the results obtained previously for $A=+\infty$, $A=0$ and $A=- \frac{1}{2}$.

preprint2020arXiv

Spatial clustering of depinning avalanches in presence of long-range interactions

Disordered elastic interfaces display avalanche dynamics at the depinning transition. For short-range interactions, avalanches correspond to compact reorganizations of the interface well described by the depinning theory. For long-range elasticity, an avalanche is a collection of spatially disconnected clusters. In this paper we determine the scaling properties of the clusters and relate them to the roughness exponent of the interface. The key observation of our analysis is the identification of a Bienaym{é}-Galton-Watson process describing the statistics of the cluster number. Our work has a concrete importance for experimental applications where the cluster statistics is a key probe of avalanche dynamics.

preprint2020arXiv

Statistics of extremes in eigenvalue-counting staircases

We consider the number ${\cal N}_{θ_A}(θ)$ of eigenvalues $e^{i θ_j}$ of a random unitary matrix, drawn from CUE$_β(N)$, in the interval $θ_j \in [θ_A,θ]$. The deviations from its mean, ${\cal N}_{θ_A}(θ) - \mathbb{E}({\cal N}_{θ_A}(θ))$, form a random process as function of $θ$. We study the maximum of this process, by exploiting the mapping onto the statistical mechanics of log-correlated random landscapes. By using an extended Fisher-Hartwig conjecture for Toeplitz determinants, supplemented with the freezing duality conjecture for log-correlated fields, we obtain the cumulants of the distribution of that maximum for any $β>0$. It exhibits combined features of standard counting statistics of fermions (free for $β=2$ and with Sutherland-type interaction for $β\ne 2$) in an interval and extremal statistics of the fractional Brownian motion with Hurst index $H=0$. The $β=2$ results are expected to apply to the statistics of zeroes of the Riemann Zeta function

preprint2020arXiv

Stochastic growth in time dependent environments

We study the Kardar-Parisi-Zhang (KPZ) growth equation in one dimension with a noise variance $c(t)$ depending on time. We find that for $c(t)\propto t^{-α}$ there is a transition at $α=1/2$. When $α>1/2$, the solution saturates at large times towards a non-universal limiting distribution. When $α<1/2$ the fluctuation field is governed by scaling exponents depending on $α$ and the limiting statistics are similar to the case when $c(t)$ is constant. We investigate this problem using different methods: (1) Elementary changes of variables mapping the time dependent case to variants of the KPZ equation with constant variance of the noise but in a deformed potential (2) An exactly solvable discretization, the log-gamma polymer model (3) Numerical simulations.

preprint2020arXiv

Tilted elastic lines with columnar and point disorder, non-Hermitian quantum mechanics and spiked random matrices: pinning and localization

We revisit the problem of an elastic line (e.g. a vortex line in a superconductor) subject to both columnar disorder and point disorder in dimension $d=1+1$. Upon applying a transverse field, a delocalization transition is expected, beyond which the line is tilted macroscopically. We investigate this transition in the fixed tilt angle ensemble and within a one-way model where backward jumps are neglected. From recent results about directed polymers and their connections to random matrix theory, we find that for a single line and a single strong defect this transition in presence of point disorder coincides with the Baik-Ben Arous-Peche (BBP) transition for the appearance of outliers in the spectrum of a perturbed random matrix in the GUE. This transition is conveniently described in the polymer picture by a variational calculation. In the delocalized phase, the ground state energy exhibits Tracy-Widom fluctuations. In the localized phase we show, using the variational calculation, that the fluctuations of the occupation length along the columnar defect are described by $f_{KPZ}$, a distribution which appears ubiquitously in the Kardar-Parisi-Zhang universality class. We then consider a smooth density of columnar defect energies. Depending on how this density vanishes at its lower edge we find either (i) a delocalized phase only (ii) a localized phase with a delocalization transition. We analyze this transition which is an infinite-rank extension of the BBP transition. The fluctuations of the ground state energy of a single elastic line in the localized phase (for fixed columnar defect energies) are described by a Fredholm determinant based on a new kernel. The case of many columns and many non-intersecting lines, relevant for the study of the Bose glass phase, is also analyzed. The ground state energy is obtained using free probability and the Burgers equation.

preprint2020arXiv

Universal scaling of the velocity field in crack front propagation

The propagation of a crack front in disordered materials is jerky and characterized by bursts of activity, called avalanches. These phenomena are the manifestation of an out-of-equilibrium phase transition originated by the disorder. As a result avalanches display universal scalings which are however difficult to characterize in experiments at finite drive. Here we show that the correlation functions of the velocity field along the front allow to extract the critical exponents of the transition and to identify the universality class of the system. We employ these correlations to characterize the universal behavior of the transition in simulations and in an experiment of crack propagation. This analysis is robust, efficient and can be extended to all systems displaying avalanche dynamics.

preprint2020arXiv

Universal survival probability for a $d$-dimensional run-and-tumble particle

We consider an active run-and-tumble particle (RTP) in $d$ dimensions and compute exactly the probability $S(t)$ that the $x$-component of the position of the RTP does not change sign up to time $t$. When the tumblings occur at a constant rate, we show that $S(t)$ is independent of $d$ for any finite time $t$ (and not just for large $t$), as a consequence of the celebrated Sparre Andersen theorem for discrete-time random walks in one dimension. Moreover, we show that this universal result holds for a much wider class of RTP models in which the speed $v$ of the particle after each tumbling is random, drawn from an arbitrary probability distribution. We further demonstrate, as a consequence, the universality of the record statistics in the RTP problem.

preprint2020arXiv

Velocity and diffusion constant of an active particle in a one dimensional force field

We consider a run an tumble particle with two velocity states $\pm v_0$, in an inhomogeneous force field $f(x)$ in one dimension. We obtain exact formulae for its velocity $V_L$ and diffusion constant $D_L$ for arbitrary periodic $f(x)$ of period $L$. They involve the "active potential" which allows to define a global bias. Upon varying parameters, such as an external force $F$, the dynamics undergoes transitions from non-ergodic trapped states, to various moving states, some with non analyticities in the $V_L$ versus $F$ curve. A random landscape in the presence of a bias leads, for large $L$, to anomalous diffusion $x \sim t^μ$, $μ<1$, or to a phase with a finite velocity that we calculate.

preprint2019arXiv

Correlations between avalanches in the depinning dynamics of elastic interfaces

We study the correlations between avalanches in the depinning dynamics of elastic interfaces driven on a random substrate. In the mean field theory (the Brownian force model), it is known that the avalanches are uncorrelated. Here we obtain a simple field theory which describes the first deviations from this uncorrelated behavior in a $ε=d_c-d$ expansion below the upper critical dimension $d_c$ of the model. We apply it to calculate the correlations between (i) avalanche sizes (ii) avalanche dynamics in two successive avalanches, or more generally, in two avalanches separated by a uniform displacement $W$ of the interface. For (i) we obtain the correlations of the total sizes, of the local sizes and of the total sizes with given seeds (starting points). For (ii) we obtain the correlations of the velocities, of the durations, and of the avalanche shapes. In general we find that the avalanches are {\it anti-correlated}, the occurence of a larger avalanche making more likely the occurence of a smaller one, and vice-versa. Examining the universality of our results leads us to conjecture several new exact scaling relations for the critical exponents that characterize the different distributions of correlations. The avalanche size predictions are confronted to numerical simulations for a $d=1$ interface with short range elasticity. They are also compared to our recent related work on static avalanches (shocks). Finally we show that the naive extrapolation of our result into the thermally activated creep regime at finite temperature, predicts strong positive correlations between the forward motion events, as recently observed in numerical simulations.

preprint2019arXiv

Delta-Bose gas on a half-line and the KPZ equation: boundary bound states and unbinding transitions

We revisit the Lieb-Liniger model for $n$ bosons in one dimension with attractive delta interaction in a half-space $\mathbb{R}^+$ with diagonal boundary conditions. This model is integrable for arbitrary value of $b \in \mathbb{R}$, the interaction parameter with the boundary. We show that its spectrum exhibits a sequence of transitions, as $b$ is decreased from the hard-wall case $b=+\infty$, with successive appearance of boundary bound states (or boundary modes) which we fully characterize. We apply these results to study the Kardar-Parisi-Zhang equation for the growth of a one-dimensional interface of height $h(x,t)$, on the half-space with boundary condition $\partial_x h(x,t)|_{x=0}=b$ and droplet initial condition at the wall. We obtain explicit expressions, valid at all time $t$ and arbitrary $b$, for the integer exponential (one-point) moments of the KPZ height field $\bar{e^{n h(0,t)}}$. From these moments we extract the large time limit of the probability distribution function (PDF) of the scaled KPZ height function. It exhibits a phase transition, related to the unbinding to the wall of the equivalent directed polymer problem, with two phases: (i) unbound for $b>-\frac{1}{2}$ where the PDF is given by the GSE Tracy-Widom distribution (ii) bound for $b<-\frac{1}{2}$, where the PDF is a Gaussian. At the critical point $b=-\frac{1}{2}$, the PDF is given by the GOE Tracy-Widom distribution.

preprint2019arXiv

Distribution of Brownian coincidences

We study the probability distribution, $P_N(T)$, of the coincidence time $T$, i.e. the total local time of all pairwise coincidences of $N$ independent Brownian walkers. We consider in details two geometries: Brownian motions all starting from $0$, and Brownian bridges. Using a Feynman-Kac representation for the moment generating function of this coincidence time, we map this problem onto some observables in three related models (i) the propagator of the Lieb Liniger model of quantum particles with pairwise delta function interactions (ii) the moments of the partition function of a directed polymer in a random medium (iii) the exponential moments of the solution of the Kardar-Parisi-Zhang equation. Using these mappings, we obtain closed formulae for the probability distribution of the coincidence time, its tails and some of its moments. Its asymptotics at large and small coincidence time are also obtained for arbitrary fixed endpoints. The universal large $T$ tail, $P_N(T) \sim \exp(- 3 T^2/(N^3-N))$ is obtained, and is independent of the geometry. We investigate the large deviations in the limit of a large number of walkers through a Coulomb gas approach. Some of our analytical results are compared with numerical simulations.

preprint2019arXiv

Entanglement entropy growth in stochastic conformal field theory and the KPZ class

We introduce a model of effective conformal quantum field theory in dimension $d=1+1$ coupled to stochastic noise, where Kardar-Parisi-Zhang (KPZ) class fluctuations can be observed. The analysis of the quantum dynamics of the scaling operators reduces to the study of random trajectories in a random environment, modeled by Brownian vector fields. We use recent results on random walks in random environments to calculate the time-dependent entanglement entropy of a subsystem interval, starting from a factorized state. We find that the fluctuations of the entropy in the large deviation regime are governed by the universal Tracy-Widom distribution. This enlarges the KPZ class, previously observed in random circuit models, to a family of interacting many body quantum systems.

preprint2019arXiv

Manifolds pinned by a high-dimensional random landscape: Hessian at the global energy minimum

We consider an elastic manifold of internal dimension $d$ and length $L$ pinned in a $N$ dimensional random potential and confined by an additional parabolic potential of curvature $μ$. We are interested in the mean spectral density $ρ(λ)$ of the Hessian matrix $K$ at the absolute minimum of the total energy. We use the replica approach to derive the system of equations for $ρ(λ)$ for a fixed $L^d$ in the $N \to \infty$ limit extending $d=0$ results of our previous work. A particular attention is devoted to analyzing the limit of extended lattice systems by letting $L\to \infty$. In all cases we show that for a confinement curvature $μ$ exceeding a critical value $μ_c$, the so-called "Larkin mass", the system is replica-symmetric and the Hessian spectrum is always gapped (from zero). The gap vanishes quadratically at $μ\to μ_c$. For $μ<μ_c$ the replica symmetry breaking (RSB) occurs and the Hessian spectrum is either gapped or extends down to zero, depending on whether RSB is 1-step or full. In the 1-RSB case the gap vanishes in all $d$ as $(μ_c-μ)^4$ near the transition. In the full RSB case the gap is identically zero. A set of specific landscapes realize the so-called "marginal cases" in $d=1,2$ which share both feature of the 1-step and the full RSB solution, and exhibit some scale invariance. We also obtain the average Green function associated to the Hessian and find that at the edge of the spectrum it decays exponentially in the distance within the internal space of the manifold with a length scale equal in all cases to the Larkin length introduced in the theory of pinning.

preprint2019arXiv

Probing the large deviations of the Kardar-Parisi-Zhang equation at short time with an importance sampling of directed polymers in random media

The one-point distribution of the height for the continuum Kardar-Parisi-Zhang (KPZ) equation is determined numerically using the mapping to the directed polymer in a random potential at high temperature. Using an importance sampling approach, the distribution is obtained over a large range of values, down to a probability density as small as $10^{-1000}$ in the tails. The short time behavior is investigated and compared with recent analytical predictions for the large-deviation forms of the probability of rare fluctuations, showing a spectacular agreement with the analytical expressions. The flat and stationary initial conditions are studied in the full space, together with the droplet initial condition in the half-space.

preprint2016arXiv

Distribution of joint local and total size and of extension for avalanches in the Brownian force model

The Brownian force model (BFM) is a mean-field model for the local velocities during avalanches in elastic interfaces of internal space dimension $d$, driven in a random medium. It is exactly solvable via a non-linear differential equation. We study avalanches following a kick, i.e. a step in the driving force. We first recall the calculation of the distributions of the global size (total swept area) and of the local jump size for an arbitrary kick amplitude. We extend this calculation to the joint density of local and global sizes within a single avalanche, in the limit of an infinitesimal kick. When the interface is driven by a single point we find new exponents $τ_0=5/3$ and $τ=7/4$, depending on whether the force or the displacement is imposed. We show that the extension of a single avalanche along one internal direction (i.e. the total length in $d=1$) is finite and we calculate its distribution, following either a local or a global kick. In all cases it exhibits a divergence $P(\ell) \sim \ell^{-3}$ at small $\ell$. Most of our results are tested in a numerical simulation in dimension $d=1$.

preprint2016arXiv

Exact solution for a random walk in a time-dependent 1D random environment: the point-to-point Beta polymer

We consider the Beta polymer, an exactly solvable model of directed polymer on the square lattice, introduced by Barraquand and Corwin. We study the statistical properties of its point to point partition sum. The problem is equivalent to a model of a random walk in a time-dependent (and in general biased) 1D random environment. In this formulation, we study the sample to sample fluctuations of the transition probability distribution function (PDF) of the random walk. Using the Bethe ansatz we obtain exact formulas for the integer moments, and Fredholm determinant formulas for the Laplace transform of the directed polymer partition sum/random walk transition probability. The asymptotic analysis of these formulas at large time $t$ is performed both (i) in a diffusive vicinity, $x \sim t^{1/2}$, of the optimal direction (in space-time) chosen by the random walk, where the fluctuations of the PDF are found to be Gamma distributed; (ii) in the large deviations regime, $x \sim t$, of the random walk, where the fluctuations of the logarithm of the PDF are found to grow with time as $t^{1/3}$ and to be distributed according to the Tracy-Widom GUE distribution. Our exact results complement those of BC for the cumulative distribution function of the random walk in regime (ii), and in regime (i) they unveil a novel fluctuation behavior. We also discuss the crossover regime between (i) and (ii), identified as $x \sim t^{3/4}$. Our results are confronted to extensive numerical simulations of the model.

preprint2016arXiv

Joint min-max distribution and Edwards-Anderson's order parameter of the circular $1 / f$-noise model

We calculate the joint min--max distribution and the Edwards-Anderson's order parameter for the circular model of $1 / f$-noise. Both quantities, as well as generalisations, are obtained exactly by combining the freezing-duality conjecture and Jack-polynomial techniques. Numerical checks come with significantly improved control of finite-size effects in the glassy phase, and the results convincingly validate the freezing-duality conjecture. Application to diffusive dynamics is discussed. We also provide a formula for the pre-factor ratio of the joint/marginal Carpentier-Le Doussal tail for minimum/maximum which applies to any logarithmic random energy model.

preprint2016arXiv

Large deviations for the height in 1D Kardar-Parisi-Zhang growth at late times

We study the atypically large deviations of the height $H \sim {\cal O}(t)$ at the origin at late times in $1+1$-dimensional growth models belonging to the Kardar-Parisi-Zhang (KPZ) universality class. We present exact results for the rate functions for the discrete single step growth model, as well as for the continuum KPZ equation in a droplet geometry. Based on our exact calculation of the rate functions we argue that models in the KPZ class undergo a third order phase transition from a strong coupling to a weak coupling phase, at late times.

preprint2016arXiv

Non-interacting fermions at finite temperature in a $d$-dimensional trap: universal correlations

We study a system of $N$ non-interacting spin-less fermions trapped in a confining potential, in arbitrary dimensions $d$ and arbitrary temperature $T$. The presence of the trap introduces an edge where the average density of fermions vanishes. Far from the edge, near the center of the trap (the so called "bulk regime"), physical properties of the fermions have traditionally been understood using the Local Density Approximation. However, this approximation drastically fails near the edge where the density vanishes. In this paper we show that, even near the edge, novel universal properties emerge, independently of the details of the confining potential. We show that for large $N$, these fermions in a confining trap, in arbitrary dimensions and at finite temperature, form a determinantal point process. As a result, any $n$-point correlation function can be expressed as an $n \times n$ determinant whose entry is called the kernel. Near the edge, we derive the large $N$ scaling form of the kernels. In $d=1$ and $T=0$, this reduces to the so called Airy kernel, that appears in the Gaussian Unitary Ensemble (GUE) of random matrix theory. In $d=1$ and $T>0$ we show a remarkable connection between our kernel and the one appearing in the $1+1$-dimensional Kardar-Parisi-Zhang equation at finite time. Consequently our result provides a finite $T$ generalization of the Tracy-Widom distribution, that describes the fluctuations of the rightmost fermion at $T=0$. In $d>1$ and $T \geq 0$, while the connection to GUE no longer holds, the process is still determinantal whose analysis provides a new class of kernels, generalizing the $1d$ Airy kernel at $T=0$ obtained in random matrix theory. Some of our finite temperature results should be testable in present-day cold atom experiments, most notably our detailed predictions for the temperature dependence of the fluctuations near the edge.

preprint2016arXiv

On integrable directed polymer models on the square lattice

In a recent work Povolotsky provided a three-parameter family of stochastic particle systems with zero-range interactions in one dimension which are integrable by coordinate Bethe ansatz. Using these results we obtain the corresponding condition for integrability of a class of directed polymer models with random weights on the square lattice. Analyzing the solutions we find, besides known cases, a new two-parameter family of integrable DP model, which we call the Inverse-Beta polymer, and provide its Bethe ansatz solution.

preprint2016arXiv

One step replica symmetry breaking and extreme order statistics of logarithmic REMs

Building upon the one-step replica symmetry breaking formalism, duly understood and ramified, we show that the sequence of ordered extreme values of a general class of Euclidean-space logarithmically correlated random energy models (logREMs) behave in the thermodynamic limit as a randomly shifted decorated exponential Poisson point process. The distribution of the random shift is determined solely by the large-distance ("infra-red", IR) limit of the model, and is equal to the free energy distribution at the critical temperature up to a translation. the decoration process is determined solely by the small-distance ("ultraviolet", UV) limit, in terms of the biased minimal process. Our approach provides connections of the replica framework to results in the probability literature and sheds further light on the freezing/duality conjecture which was the source of many previous results for log-REMs. In this way we derive the general and explicit formulae for the joint probability density of depths of the first and second minima (as well its higher-order generalizations) in terms of model-specific contributions from UV as well as IR limits. In particular, we show that the second min statistics is largely independent of details of UV data, whose influence is seen only through the mean value of the gap. For a given log-correlated field this parameter can be evaluated numerically, and we provide several numerical tests of our theory using the circular model of $1/f$-noise.

preprint2016arXiv

The crossing probability for directed polymers in random media

We study the probability that two directed polymers in the same random potential do not intersect. We use the replica method to map the problem onto the attractive Lieb-Liniger model with generalized statistics between particles. We obtain analytical expressions for the first few moments of this probability, and compare them to a numerical simulation of a discrete model at high-temperature. From these observations, several large time properties of the non-crossing probabilities are conjectured. Extensions of our formalism to more general observables are discussed.

preprint2016arXiv

Universal correlations between shocks in the ground state of elastic interfaces in disordered media

The ground state of an elastic interface in a disordered medium undergoes collective jumps upon variation of external parameters. These mesoscopic jumps are called shocks, or static avalanches. Submitting the interface to a parabolic potential centered at $w$, we study the avalanches which occur as $w$ is varied. We are interested in the correlations between the avalanche sizes $S_1$ and $S_2$ occurring at positions $w_1$ and $w_2$. Using the Functional Renormalization Group (FRG), we show that correlations exist for realistic interface models below their upper critical dimension. Notably, the connected moment $ \langle S_1 S_2 \rangle^c$ is up to a prefactor exactly the renormalized disorder correlator, itself a function of $|w_2-w_1|$. The latter is the universal function at the center of the FRG; hence correlations between shocks are universal as well. All moments and the full joint probability distribution are computed to first non-trivial order in an $ε$-expansion below the upper critical dimension. To quantify the local nature of the coupling between avalanches, we calculate the correlations of their local jumps. We finally test our predictions against simulations of a particle in random-bond and random-force disorder, with surprisingly good agreement.

preprint2016arXiv

Universal ground state properties of free fermions in a $d$-dimensional trap

The ground state properties of $N$ spinless free fermions in a $d$-dimensional confining potential are studied. We find that any $n$-point correlation function has a simple determinantal structure that allows us to compute several properties exactly for large $N$. We show that the average density has a finite support with an edge, and near this edge the density exhibits a universal (valid for a wide class of potentials) scaling behavior for large $N$. The associated edge scaling function is computed exactly and generalizes to any $d$ the edge electron gas result of Kohn and Mattsson in $d=3$ [Phys. Rev. Lett. 81, 3487 (1998)]. In addition, we calculate the kernel (that characterizes any $n$-point correlation function) for large $N$ and show that, when appropriately scaled, it depends only on dimension $d$, but has otherwise universal scaling forms, at the edges. The edge kernel, for higher $d$, generalizes the Airy kernel in one dimension, well known from random matrix theory.

preprint2016arXiv

Universality in the mean spatial shape of avalanches

Quantifying the universality of avalanche observables beyond critical exponents is of current great interest in theory and experiments. Here, we improve the characterization of the spatio-temporal process inside avalanches in the universality class of the depinning of elastic interfaces in random media. Surprisingly, at variance with the temporal shape, the spatial shape of avalanches has not yet been predicted. In part this is due to a lack of an analytically tractable definition: how should the shapes be centered? Here we introduce such a definition, accessible in experiments, and study the mean spatial shape of avalanches at fixed size centered around their starting point (seed). We calculate the associated universal scaling functions, both in a mean-field model and beyond. Notably, they are predicted to exhibit a cusp singularity near the seed. The results are in good agreement with a numerical simulation of an elastic line.

preprint2015arXiv

Crossing probability for directed polymers in random media: exact tail of the distribution

We study the probability $p \equiv p_η(t)$ that two directed polymers in a given random potential $η$ and with fixed and nearby endpoints, do not cross until time $t$. This probability is itself a random variable (over samples $η$) which, as we show, acquires a very broad probability distribution at large time. In particular the moments of $p$ are found to be dominated by atypical samples where $p$ is of order unity. Building on a formula established by us in a previous work using nested Bethe Ansatz and Macdonald process methods, we obtain analytically the leading large time behavior of {\it all moments} $\overline{p^m}\simeq γ_m/t$. From this, we extract the exact tail $\sim ρ(p)/t$ of the probability distribution of the non-crossing probability at large time. The exact formula is compared to numerical simulations, with excellent agreement.

preprint2015arXiv

Experimental evidence for three universality classes for reaction fronts in disordered flows

Self-sustained reaction fronts in a disordered medium subject to an external flow display self-affine roughening, pinning and depinning transitions. We measure spatial and temporal fluctuations of the front in $1+1$ dimensions, controlled by a single parameter, the mean flow velocity. Three distinct universality classes are observed, consistent with the Kardar-Parisi-Zhang (KPZ) class for fast advancing or receding fronts, the quenched KPZ class (positive-qKPZ) when the mean flow approximately cancels the reaction rate, and the negative-qKPZ class for slowly receding fronts. Both quenched KPZ classes exhibit distinct depinning transitions, in agreement with the theory.

preprint2015arXiv

Finite temperature free fermions and the Kardar-Parisi-Zhang equation at finite time

We consider the system of $N$ one-dimensional free fermions confined by a harmonic well $V(x) = mω^2 {x^2}/{2}$ at finite inverse temperature $β= 1/T$. The average density of fermions $ρ_N(x,T)$ at position $x$ is derived. For $N \gg 1$ and $β\sim {\cal O}(1/N)$, $ρ_N(x,T)$ is given by a scaling function interpolating between a Gaussian at high temperature, for $β\ll 1/N$, and the Wigner semi-circle law at low temperature, for $β\gg N^{-1}$. In the latter regime, we unveil a scaling limit, for $β{\hbar ω}= b N^{-1/3}$, where the fluctuations close to the edge of the support, at $x \sim \pm \sqrt{2\hbar N/(mω)}$, are described by a limiting kernel $K^{\rm ff}_b(s,s')$ that depends continuously on $b$ and is a generalization of the Airy kernel, found in the Gaussian Unitary Ensemble of random matrices. Remarkably, exactly the same kernel $K^{\rm ff}_b(s,s')$ arises in the exact solution of the Kardar-Parisi-Zhang (KPZ) equation in 1+1 dimensions at finite time $t$, with the correspondence $t= b^3$.

preprint2015arXiv

Spatial shape of avalanches in the Brownian force model

We study the Brownian force model (BFM), a solvable model of avalanche statistics for an interface, in a general discrete setting. The BFM describes the overdamped motion of elastically coupled particles driven by a parabolic well in independent Brownian force landscapes. Avalanches are defined as the collective jump of the particles in response to an arbitrary monotonous change in the well position (i.e. in the applied force). We derive an exact formula for the joint probability distribution of these jumps. From it we obtain the joint density of local avalanche sizes for stationary driving in the quasi-static limit near the depinning threshold. A saddle-point analysis predicts the spatial shape of avalanches in the limit of large aspect ratios for the continuum version of the model. We then study fluctuations around this saddle point, and obtain the leading corrections to the mean shape, the fluctuations around the mean shape and the shape asymmetry, for finite aspect ratios. Our results are finally confronted to numerical simulations.

preprint2014arXiv

An exact mapping of the stochastic field theory for Manna sandpiles to interfaces in random media

We show that the stochastic field theory for directed percolation in presence of an additional conservation law (the C-DP class) can be mapped exactly to the continuum theory for the depinning of an elastic interface in short-range correlated quenched disorder. On one line of parameters commonly studied, this mapping leads to the simplest overdamped dynamics. Away from this line, an additional memory term arises in the interface dynamics; we argue that it does not change the universality class. Since C-DP is believed to describe the Manna class of self-organized criticality, this shows that Manna stochastic sandpiles and disordered elastic interfaces (i.e. the quenched Edwards-Wilkinson model) share the same universal large-scale behavior.

preprint2014arXiv

Avalanche shape and exponents beyond mean-field theory

Elastic systems, such as magnetic domain walls, density waves, contact lines, and cracks, are all pinned by substrate disorder. When driven, they move via successive jumps called avalanches, with power law distributions of size, duration and velocity. Their exponents, and the shape of an avalanche, defined as its mean velocity as function of time, have recently been studied. They are known approximatively from experiments and simulations, and were predicted from mean-field models, such as the Brownian force model (BFM), where each point of the elastic interface sees a force field which itself is a random walk. As we showed in EPL 97 (2012) 46004, the BFM is the starting point for an $ε= d_{\rm c}-d$ expansion around the upper critical dimension, with $d_{\rm c}=4$ for short-ranged elasticity, and $d_{\rm c}=2$ for long-ranged elasticity. Here we calculate analytically the ${\cal O}(ε)$, i.e. 1-loop, correction to the avalanche shape at fixed duration $T$, for both types of elasticity. The exact expression is well approximated by $\left< \dot u(t=x T)\right>_T\simeq [ Tx(1-x)]^{γ-1} \exp\left( {\cal A}\left[\frac12-x\right]\right)$, $0<x<1$. The asymmetry ${\cal A}\approx - 0.336 (1-d/d_{\rm c})$ is negative for $d$ close to $d_{\rm c}$, skewing the avalanche towards its end, as observed in numerical simulations in $d=2$ and $3$. The exponent $γ=(d+ζ)/z$ is given by the two independent exponents at depinning, the roughness $ζ$ and the dynamical exponent $z$. We propose a general procedure to predict other avalanche exponents in terms of $ζ$ and $z$. We finally introduce and calculate the shape at fixed avalanche size, not yet measured in experiments or simulations.

preprint2014arXiv

Crossover from droplet to flat initial conditions in the KPZ equation from the replica Bethe ansatz

We show how our previous result based on the replica Bethe ansatz for the Kardar Parisi Zhang (KPZ) equation with the "half-flat" initial condition leads to the Airy$_2$ to Airy$_1$ (i.e. GUE to GOE) universal crossover one-point height distribution in the limit of large time. Equivalently, we obtain the distribution of the free energy of a long directed polymer (DP) in a random potential with one fixed endpoint and the other one on a half-line. We then generalize to a DP when each endpoint is free on its own half-line. It amounts, in the limit of large time, to obtain the distribution of the maximum of the transition process Airy$_{2\to 1}$ (minus a half-parabola) on a half line.

preprint2014arXiv

From the sine-Gordon field theory to the Kardar-Parisi-Zhang growth equation

We unveil a remarkable connection between the sine-Gordon quantum field theory and the Kardar-Parisi-Zhang (KPZ) growth equation. We find that the non-relativistic limit of the two point correlation function of the sine-Gordon theory is related to the generating function of the height distribution of the KPZ field with droplet initial conditions, i.e. the directed polymer free energy with two endpoints fixed. As shown recently, the latter can be expressed as a Fredholm determinant which in the large time separation limit converges to the GUE Tracy-Widom cumulative distribution. Possible applications and extensions are discussed.

preprint2014arXiv

Interaction quench in a Lieb-Liniger model and the KPZ equation with flat initial conditions

Recent exact solutions of the 1D Kardar-Parisi-Zhang equation make use of the 1D integrable Lieb-Liniger model of interacting bosons. For flat initial conditions, it requires the knowledge of the overlap between the uniform state and arbitrary exact Bethe eigenstates. The same quantity is also central in the study of the quantum quench from a 1D non-interacting Bose-Einstein condensate upon turning interactions. We compare recent advances in both domains, i.e. our previous exact solution, and a new conjecture by De Nardis et al.. This leads to new exact results and conjectures for both the quantum quench and the KPZ problem.

preprint2014arXiv

Log-Gamma directed polymer with fixed endpoints via the replica Bethe Ansatz

We study the model of a discrete directed polymer (DP) on the square lattice with homogeneous inverse gamma distribution of site random Boltzmann weights, introduced by Seppalainen. The integer moments of the partition sum, $\overline{Z^n}$, are studied using a transfer matrix formulation, which appears as a generalization of the Lieb-Liniger quantum mechanics of bosons to discrete time and space. In the present case of the inverse gamma distribution the model is integrable in terms of a coordinate Bethe Ansatz, as discovered by Brunet. Using the Brunet-Bethe eigenstates we obtain an exact expression for the integer moments of $\overline{Z^n}$ for polymers of arbitrary lengths and fixed endpoint positions. Although these moments do not exist for all integer n, we are nevertheless able to construct a generating function which reproduces all existing integer moments, and which takes the form of a Fredholm determinant (FD). This suggests an analytic continuation via a Mellin-Barnes transform and we thereby propose a FD ansatz representation for the probability distribution function (PDF) of $Z$ and its Laplace transform. In the limit of very long DP, this ansatz yields that the distribution of the free energy converges to the GUE Tracy-Widom distribution up to a non-trivial average and variance that we calculate. Our asymptotic predictions coincide with a result by Borodin et al. based on a formula obtained by Seppalainen using the gRSK correspondence. In addition we obtain the dependence on the endpoint position and the exact elastic coefficient at large time. We argue the equivalence between our formula and the one of Borodin et al. As we discuss, this open the way to explore the connections between quantum integrability and tropical geometry.

preprint2014arXiv

Random matrices and entanglement entropy of trapped Fermi gases

We exploit and clarify the use of random matrix theory for the calculation of the entanglement entropy of free Fermi gases. We apply this method to obtain analytic predictions for Renyi entanglement entropies of a one-dimensional gas trapped by a harmonic potential in all the relevant scaling regimes. We confirm our findings with accurate numerical calculations obtained by means of an ingenious discretisation of the reduced correlation matrix.

preprint2014arXiv

Revisiting Directed Polymers with heavy-tailed disorder

In this mostly numerical study, we revisit the statistical properties of the ground state of a directed polymer in a $d=1+1$ "hilly" disorder landscape, i.e. when the quenched disorder has power-law tails. When disorder is Gaussian, the polymer minimizes its total energy through a collective optimization, where the energy of each visited site only weakly contributes to the total. Conversely, a hilly landscape forces the polymer to distort and explore a larger portion of space to reach some particularly deep energy sites. As soon as the fifth moment of the disorder diverges, this mechanism radically changes the standard "KPZ" scaling behaviour of the directed polymer, and new exponents prevail. After confirming again that the Flory argument accurately predicts these exponent in the tail-dominated phase, we investigate several other statistical features of the ground state that shed light on this unusual transition and on the accuracy of the Flory argument. We underline the theoretical challenge posed by this situation, which paradoxically becomes even more acute above the upper critical dimension.

preprint2014arXiv

Superfluid/Bose-glass transition in one dimension

We consider a one-dimensional system of interacting bosons in a random potential. At zero temperature, it can be either in the superfluid or in the insulating phase. We study the transition at weak disorder and moderate interaction. Using a systematic approach, we derive the renormalization group equations at two-loop order and discuss the phase diagram. We find the universal form of the correlation functions at the transitions and compute the logarithmic corrections to the main universal power-law behavior. In order to mimic large density fluctuations on a single site, we study a simplified model of disordered two-leg bosonic ladders with correlated disorder across the rung. Contrarily to the single-chain case, the latter system exhibits a transition between a superfluid and a localized phase where the exponents of the correlation functions at the transition do not take universal values.

preprint2013arXiv

Avalanche dynamics of elastic interfaces

Slowly driven elastic interfaces, such as domain walls in dirty magnets, contact lines, or cracks proceed via intermittent motion, called avalanches. We develop a field-theoretic treatment to calculate, from first principles, the space-time statistics of instantaneous velocities within an avalanche. For elastic interfaces at (or above) their (internal) upper critical dimension d >= d_uc (d_uc = 2, 4 respectively for long-ranged and short-ranged elasticity) we show that the field theory for the center of mass reduces to the motion of a point particle in a random-force landscape, which is itself a random walk (ABBM model). Furthermore, the full spatial dependence of the velocity correlations is described by the Brownian-force model (BFM) where each point of the interface sees an independent Brownian-force landscape. Both ABBM and BFM can be solved exactly in any dimension d (for monotonous driving) by summing tree graphs, equivalent to solving a (non-linear) instanton equation. This tree approximation is the mean-field theory (MFT) for realistic interfaces in short-ranged disorder. Both for the center of mass, and for a given Fourier mode q, we obtain probability distribution functions (PDF's) of the velocity, as well as the avalanche shape and its fluctuations (second shape). Within MFT we find that velocity correlations at non-zero q are asymmetric under time reversal. Next we calculate, beyond MFT, i.e. including loop corrections, the 1-time PDF of the center-of-mass velocity du/dt for dimension d< d_uc. The singularity at small velocity P(du/dt) ~ 1/(du/dt)^a is substantially reduced from a=1 (MFT) to a = 1 - 2/9 (4-d) + ... (short-ranged elasticity) and a = 1 - 4/9 (2-d) + ... (long-ranged elasticity). We show how the dynamical theory recovers the avalanche-size distribution, and how the instanton relates to the response to an infinitesimal step in the force.

preprint2013arXiv

Collective excitations in a large-d model for graphene

We consider a model of Dirac fermions coupled to flexural phonons to describe a graphene sheet fluctuating in dimension $2+d$. We derive the self-consistent screening equations for the quantum problem, exact in the limit of large $d$. We first treat the membrane alone, and work out the quantum to classical, and harmonic to anharmonic crossover. For the coupled electron-membrane problem we calculate the dressed two-particle propagators of the elastic and electron interactions and find that it exhibits a collective mode which becomes unstable at some wave-vector $q_{\rm c}$ for large enough coupling $g$. The saddle point analysis, exact at large $d$, indicates that this instability corresponds to spontaneous and simultaneous appearance of gaussian curvature and electron puddles. The relevance to ripples in graphene is discussed.

preprint2013arXiv

Exact form of the exponential correlation function in the glassy super-rough phase

We consider the random-phase sine-Gordon model in two dimensions. It describes two-dimensional elastic systems with random periodic disorder, such as pinned flux-line arrays, random field XY models, and surfaces of disordered crystals. The model exhibits a super-rough glass phase at low temperature $T<T_{c}$ with relative displacements growing with distance $r$ as $\bar{\langle [θ(r)-θ(0)]^2\rangle} \simeq A(τ) \ln^2 (r/a)$, where $A(τ) = 2 τ^2- 2 τ^3 +\mathcal{O}(τ^4)$ near the transition and $τ=1-T/T_{c}$. We calculate all higher cumulants and show that they grow as $\bar{\langle[θ(r)-θ(0)]^{2n}\rangle}_c \simeq [2 (-1)^{n+1} (2n)! ζ(2n-1) τ^2 + \mathcal{O}(τ^3) ] \ln(r/a)$, $n \geq 2$, where $ζ$ is the Riemann zeta function. By summation, we obtain the decay of the exponential correlation function as $\bar{\langle e^{iq\left[θ(r)-θ(0)\right]}\rangle} \simeq (a/r)^{η(q)} \exp\boldsymbol(-\frac{1}{2}\mathcal{A}(q)\ln^2(r/a)\boldsymbol)$ where $η(q)$ and ${\cal A}(q)$ are obtained for arbitrary $q \leq 1$ to leading order in $τ$. The anomalous exponent is $η(q) = c q^2 - τ^2 q^2 [2γ_E+ψ(q)+ψ(-q)]$ in terms of the digamma function $ψ$, where $c$ is non-universal and $γ_E$ is the Euler constant. The correlation function shows a faster decay at $q=1$, corresponding to fermion operators in the dual picture, which should be visible in Bragg scattering experiments.

preprint2013arXiv

Non-Gaussian effects and multifractality in the Bragg glass

We study, beyond the Gaussian approximation, the decay of the translational order correlation function for a d-dimensional scalar periodic elastic system in a disordered environment. We develop a method based on functional determinants, equivalent to summing an infinite set of diagrams. We obtain, in dimension d=4-epsilon, the even n-th cumulant of relative displacements as <[u(r)-u(0)]^n>^c = A_n ln r, with A_n = -(ε/3)^n Γ(n-1/2) ζ(2n-3)/π^(1/2), as well as the multifractal dimension x_q of the exponential field e^{q u(r)}. As a corollary, we obtain an analytic expression for a class of n-loop integrals in d=4, which appear in the perturbative determination of Konishi amplitudes, also accessible via AdS/CFT using integrability.

preprint2013arXiv

Statistics of Avalanches with Relaxation, and Barkhausen Noise: A Solvable Model

We study a generalization of the Alessandro-Beatrice-Bertotti-Montorsi (ABBM) model of a particle in a Brownian force landscape, including retardation effects. We show that under monotonous driving the particle moves forward at all times, as it does in absence of retardation (Middleton's theorem). This remarkable property allows us to develop an analytical treatment. The model with an exponentially decaying memory kernel is realized in Barkhausen experiments with eddy-current relaxation, and has previously been shown numerically to account for the experimentally observed asymmetry of Barkhausen-pulse shapes. We elucidate another qualitatively new feature: the breakup of each avalanche of the standard ABBM model into a cluster of sub-avalanches, sharply delimited for slow relaxation under quasi-static driving. These conditions are typical for earthquake dynamics. With relaxation and aftershock clustering, the present model includes important ingredients for an effective description of earthquakes. We analyze quantitatively the limits of slow and fast relaxation for stationary driving with velocity v>0. The v-dependent power-law exponent for small velocities, and the critical driving velocity at which the particle velocity never vanishes, are modified. We also analyze non-stationary avalanches following a step in the driving magnetic field. Analytically, we obtain the mean avalanche shape at fixed size, the duration distribution of the first sub-avalanche, and the time dependence of the mean velocity. We propose to study these observables in experiments, allowing to directly measure the shape of the memory kernel, and to trace eddy current relaxation in Barkhausen noise.

preprint2013arXiv

Statistics of shocks in a toy model with heavy tails

We study the energy minimization for a particle in a quadratic well in presence of short-ranged heavy-tailed disorder, as a toy model for an elastic manifold. The discrete model is shown to be described in the scaling limit by a continuum Poisson process model which captures the three universality classes. This model is solved in general, and we give, in the present case (Frechet class), detailed results for the distribution of the minimum energy and position, and the distribution of the sizes of the shocks (i.e. switches in the ground state) which arise as the position of the well is varied. All these distributions are found to exhibit heavy tails with modified exponents. These results lead to an "exotic regime" in Burgers turbulence decaying from a heavy-tailed initial condition.

preprint2013arXiv

Topology trivialization and large deviations for the minimum in the simplest random optimization

Finding the global minimum of a cost function given by the sum of a quadratic and a linear form in N real variables over (N-1)- dimensional sphere is one of the simplest, yet paradigmatic problems in Optimization Theory known as the "trust region subproblem" or "constraint least square problem". When both terms in the cost function are random this amounts to studying the ground state energy of the simplest spherical spin glass in a random magnetic field. We first identify and study two distinct large-N scaling regimes in which the linear term (magnetic field) leads to a gradual topology trivialization, i.e. reduction in the total number N_{tot} of critical (stationary) points in the cost function landscape. In the first regime N_{tot} remains of the order $N$ and the cost function (energy) has generically two almost degenerate minima with the Tracy-Widom (TW) statistics. In the second regime the number of critical points is of the order of unity with a finite probability for a single minimum. In that case the mean total number of extrema (minima and maxima) of the cost function is given by the Laplace transform of the TW density, and the distribution of the global minimum energy is expected to take a universal scaling form generalizing the TW law. Though the full form of that distribution is not yet known to us, one of its far tails can be inferred from the large deviation theory for the global minimum. In the rest of the paper we show how to use the replica method to obtain the probability density of the minimum energy in the large-deviation approximation by finding both the rate function and the leading pre-exponential factor.

preprint2012arXiv

Counting function fluctuations and extreme value threshold in multifractal patterns: the case study of an ideal $1/f$ noise

To understand the sample-to-sample fluctuations in disorder-generated multifractal patterns we investigate analytically as well as numerically the statistics of high values of the simplest model - the ideal periodic $1/f$ Gaussian noise. By employing the thermodynamic formalism we predict the characteristic scale and the precise scaling form of the distribution of number of points above a given level. We demonstrate that the powerlaw forward tail of the probability density, with exponent controlled by the level, results in an important difference between the mean and the typical values of the counting function. This can be further used to determine the typical threshold $x_m$ of extreme values in the pattern which turns out to be given by $x_m^{(typ)}=2-c\ln{\ln{M}}/\ln{M}$ with $c=3/2$. Such observation provides a rather compelling explanation of the mechanism behind universality of $c$. Revealed mechanisms are conjectured to retain their qualitative validity for a broad class of disorder-generated multifractal fields. In particular, we predict that the typical value of the maximum $p_{max}$ of intensity is to be given by $-\ln{p_{max}} = α_{-}\ln{M} + \frac{3}{2f'(α_{-})}\ln{\ln{M}} + O(1)$, where $f(α)$ is the corresponding singularity spectrum vanishing at $α=α_{-}>0$. For the $1/f$ noise we also derive exact as well as well-controlled approximate formulas for the mean and the variance of the counting function without recourse to the thermodynamic formalism.

preprint2012arXiv

Directed polymer near a hard wall and KPZ equation in the half-space

We study the directed polymer with fixed endpoints near an absorbing wall, in the continuum and in presence of disorder, equivalent to the KPZ equation on the half space with droplet initial conditions. From a Bethe Ansatz solution of the equivalent attractive boson model we obtain the exact expression for the free energy distribution at all times. It converges at large time to the Tracy Widom distribution $F_4$ of the Gaussian Symplectic Ensemble (GSE). We compare our results with numerical simulations of the lattice directed polymer, both at zero and high temperature.

preprint2012arXiv

Distribution of velocities and acceleration for a particle in Brownian correlated disorder: inertial case

We study the motion of an elastic object driven in a disordered environment in presence of both dissipation and inertia. We consider random forces with the statistics of random walks and reduce the problem to a single degree of freedom. It is the extension of the mean field ABBM model in presence of an inertial mass m. While the ABBM model can be solved exactly, its extension to inertia exhibits complicated history dependence due to oscillations and backward motion. The characteristic scales for avalanche motion are studied from numerics and qualitative arguments. To make analytical progress we consider two variants which coincide with the original model whenever the particle moves only forward. Using a combination of analytical and numerical methods together with simulations, we characterize the distributions of instantaneous acceleration and velocity, and compare them in these three models. We show that for large driving velocity, all three models share the same large-deviation function for positive velocities, which is obtained analytically for small and large m, as well as for m =6/25. The effect of small additional thermal and quantum fluctuations can be treated within an approximate method.

preprint2012arXiv

Functional renormalization-group approach to decaying turbulence

We reconsider the functional renormalization-group (FRG) approach to decaying Burgers turbulence, and extend it to decaying Navier-Stokes and Surface-Quasi-Geostrophic turbulence. The method is based on a renormalized small-time expansion, equivalent to a loop expansion, and naturally produces a dissipative anomaly and a cascade after a finite time. We explicitly calculate and analyze the one-loop FRG equations in the zero-viscosity limit as a function of the dimension. For Burgers they reproduce the FRG equation obtained in the context of random manifolds, extending previous results of one of us. Breakdown of energy conservation due to shocks and the appearance of a direct energy cascade corresponds to failure of dimensional reduction in the context of disordered systems. For Navier-Stokes in three dimensions, the velocity-velocity correlation function acquires a linear dependence on the distance, zeta_2=1, in the inertial range, instead of Kolmogorov's zeta_2=2/3; however the possibility remains for corrections at two- or higher-loop order. In two dimensions, we obtain a numerical solution which conserves energy and exhibits an inverse cascade, with explicit analytical results both for large and small distances, in agreement with the scaling proposed by Batchelor. In large dimensions, the one-loop FRG equation for Navier-Stokes converges to that of Burgers.

preprint2012arXiv

Nonstationary dynamics of the Alessandro-Beatrice-Bertotti-Montorsi model

We obtain an exact solution for the motion of a particle driven by a spring in a Brownian random-force landscape, the Alessandro-Beatrice-Bertotti-Montorsi (ABBM) model. Many experiments on quasi-static driving of elastic interfaces (Barkhausen noise in magnets, earthquake statistics, shear dynamics of granular matter) exhibit the same universal behavior as this model. It also appears as a limit in the field theory of elastic manifolds. Here we discuss predictions of the ABBM model for monotonous, but otherwise arbitrary, time-dependent driving. Our main result is an explicit formula for the generating functional of particle velocities and positions. We apply this to derive the particle-velocity distribution following a quench in the driving velocity. We also obtain the joint avalanche size and duration distribution and the mean avalanche shape following a jump in the position of the confining spring. Such non-stationary driving is easy to realize in experiments, and provides a way to test the ABBM model beyond the stationary, quasi-static regime. We study extensions to two elastically coupled layers, and to an elastic interface of internal dimension d, in the Brownian force landscape. The effective action of the field theory is equal to the action, up to 1-loop corrections obtained exactly from a functional determinant. This provides a connection to renormalization-group methods.

preprint2012arXiv

Phase transition of interacting disordered bosons in one dimension

Interacting bosons generically form a superfluid state. In the presence of disorder it can get converted into a compressible Bose glass state. Here we study such transition in one dimension at moderate interaction using bosonization and renormalization group techniques. We derive the two-loop scaling equations and discuss the phase diagram. We find that the correlation functions at the transition are characterized by universal exponents in a finite region around the fixed point.

preprint2012arXiv

Rings and Coulomb boxes in dissipative environments

We study a particle on a ring in presence of a dissipative Caldeira-Leggett environment and derive its response to a DC field. We show how this non-equilibrium response is related to a flux averaged equilibrium response. We find, through a 2-loop renormalization group analysis, that a large dissipation parameter ηflows to a fixed point η^R=\hbar/2π. We also reexamine the mapping of this problem to that of the Coulomb box and show that the relaxation resistance, of recent interest, is quantized for large η. For finite η>η^R we find that a certain average of the relaxation resistance is quantized. We propose a Coulomb box experiment to measure a quantized noise.

preprint2012arXiv

Short time growth of a KPZ interface with flat initial conditions

The short time behavior of the 1+1 dimensional KPZ growth equation with a flat initial condition is obtained from the exact expressions of the moments of the partition function of a directed polymer with one endpoint free and the other fixed. From these expressions, the short time expansions of the lowest cumulants of the KPZ height field are exactly derived. The results for these two classes of cumulants are checked in high precision lattice numerical simulations. The short time limit considered here is relevant for the study of the interface growth in the large diffusivity/weak noise limit, and describes the universal crossover between the Edwards-Wilkinson and KPZ universality classes for an initially flat interface.

preprint2012arXiv

Super-Rough Glassy Phase of the Random Field XY Model in Two Dimensions

We study both analytically, using the renormalization group (RG) to two loop order, and numerically, using an exact polynomial algorithm, the disorder-induced glass phase of the two-dimensional XY model with quenched random symmetry-breaking fields and without vortices. In the super-rough glassy phase, i.e. below the critical temperature $T_c$, the disorder and thermally averaged correlation function $B(r)$ of the phase field $θ(x)$, $B(r) = \bar{<[θ(x) - θ(x+ r) ]^2>}$ behaves, for $r \gg a$, as $B(r) \simeq A(τ) \ln^2 (r/a)$ where $r = |r|$ and $a$ is a microscopic length scale. We derive the RG equations up to cubic order in $τ= (T_c-T)/T_c$ and predict the universal amplitude ${A}(τ) = 2τ^2-2τ^3 + {\cal O}(τ^4)$. The universality of $A(τ)$ results from nontrivial cancellations between nonuniversal constants of RG equations. Using an exact polynomial algorithm on an equivalent dimer version of the model we compute ${A}(τ)$ numerically and obtain a remarkable agreement with our analytical prediction, up to $τ\approx 0.5$.

preprint2012arXiv

Super-rough phase of the random-phase sine-Gordon model: Two-loop results

We consider the two-dimensional random-phase sine-Gordon and study the vicinity of its glass transition temperature $T_c$, in an expansion in small $τ=(T_c-T)/T_c$, where $T$ denotes the temperature. We derive renormalization group equations in cubic order in the anharmonicity, and show that they contain two universal invariants. Using them we obtain that the correlation function in the super-rough phase for temperature $T<T_c$ behaves at large distances as $\bar{<[θ(x)-θ(0)]^2>} = \mathcal{A}\ln^2(|x|/a) + \mathcal{O}[\ln(|x|/a)]$, where the amplitude $\mathcal{A}$ is a universal function of temperature $\mathcal{A}=2τ^2-2τ^3+\mathcal{O}(τ^4)$. This result differs at two-loop order, i.e., $\mathcal{O}(τ^3)$, from the prediction based on results from the "nearly conformal" field theory of a related fermion model. We also obtain the correction-to-scaling exponent.

preprint2012arXiv

The KPZ equation with flat initial condition and the directed polymer with one free end

We study the directed polymer (DP) of length $t$ in a random potential in dimension 1+1 in the continuum limit, with one end fixed and one end free. This maps onto the Kardar-Parisi-Zhang growth equation in time $t$, with flat initial conditions. We use the Bethe Ansatz solution for the replicated problem which is an attractive bosonic model. The problem is more difficult than the previous solution of the fixed endpoint problem as it requires regularization of the spatial integrals over the Bethe eigenfunctions. We use either a large fixed system length or a small finite slope KPZ initial conditions (wedge). The latter allows to take properly into account non-trivial contributions, which appear as deformed strings in the former. By considering a half-space model in a proper limit we obtain an expression for the generating function of all positive integer moments $\bar{Z^n}$ of the directed polymer partition function. We obtain the generating function of the moments of the DP partition sum as a Fredholm Pfaffian. At large time, this Fredholm Pfaffian, valid for all time $t$, exhibits convergence of the free energy (i.e. KPZ height) distribution to the GOE Tracy Widom distribution

preprint2011arXiv

An exact solution for the KPZ equation with flat initial conditions

We provide the first exact calculation of the height distribution at arbitrary time $t$ of the continuum KPZ growth equation in one dimension with flat initial conditions. We use the mapping onto a directed polymer (DP) with one end fixed, one free, and the Bethe Ansatz for the replicated attractive boson model. We obtain the generating function of the moments of the DP partition sum as a Fredholm Pfaffian. Our formula, valid for all times, exhibits convergence of the free energy (i.e. KPZ height) distribution to the GOE Tracy Widom distribution at large time.

preprint2011arXiv

Distribution of velocities in an avalanche

For a driven elastic object near depinning, we derive from first principles the distribution of instantaneous velocities in an avalanche. We prove that above the upper critical dimension, d >= d_uc, the n-times distribution of the center-of-mass velocity is equivalent to the prediction from the ABBM stochastic equation. Our method allows to compute space and time dependence from an instanton equation. We extend the calculation beyond mean field, to lowest order in epsilon=d_uc-d.

preprint2011arXiv

Equilibrium avalanches in spin glasses

We study the distribution of equilibrium avalanches (shocks) in Ising spin glasses which occur at zero temperature upon small changes in the magnetic field. For the infinite-range Sherrington-Kirkpatrick model we present a detailed derivation of the density rho(Delta M) of the magnetization jumps Delta M. It is obtained by introducing a multi-component generalization of the Parisi-Duplantier equation, which allows us to compute all cumulants of the magnetization. We find that rho(Delta M) ~ (Delta M)^(-tau) with an avalanche exponent tau=1 for the SK model, originating from the marginal stability (criticality) of the model. It holds for jumps of size 1 << Delta M < N^(1/2) being provoked by changes of the external field by delta H = O(N^[-1/2]) where N is the total number of spins. Our general formula also suggests that the density of overlap q between initial and final state in an avalanche is rho(q) ~ 1/(1-q). These results show interesting similarities with numerical simulations for the out-of-equilibrium dynamics of the SK model. For finite-range models, using droplet arguments, we obtain the prediction tau= (d_f + theta)/d_m, where d_f,d_m and theta are the fractal dimension, magnetization exponent and energy exponent of a droplet, respectively. This formula is expected to apply to other glassy disordered systems, such as the random-field model and pinned interfaces. We make suggestions for further numerical investigations, as well as experimental studies of the Barkhausen noise in spin glasses.

preprint2011arXiv

First-principle derivation of static avalanche-size distribution

We study the energy minimization problem for an elastic interface in a random potential plus a quadratic well. As the position of the well is varied, the ground state undergoes jumps, called shocks or static avalanches. We introduce an efficient and systematic method to compute the statistics of avalanche sizes and manifold displacements. The tree-level calculation, i.e. mean-field limit, is obtained by solving a saddle-point equation. Graphically, it can be interpreted as a the sum of all tree graphs. The 1-loop corrections are computed using results from the functional renormalization group. At the upper critical dimension the shock statistics is described by the Brownian Force model (BFM), the static version of the so-called ABBM model in the non-equilibrium context of depinning. This model can itself be treated exactly in any dimension and its shock statistics is that of a Levy process. Contact is made with classical results in probability theory on the Burgers equation with Brownian initial conditions. In particular we obtain a functional extension of an evolution equation introduced by Carraro and Duchon, which recursively constructs the tree diagrams in the field theory.

preprint2011arXiv

Interference in disordered systems: A particle in a complex random landscape

We consider a particle in one dimension submitted to amplitude and phase disorder. It can be mapped onto the complex Burgers equation, and provides a toy model for problems with interplay of interferences and disorder, such as the NSS model of hopping conductivity in disordered insulators and the Chalker-Coddington model for the (spin) quantum Hall effect. The model has three distinct phases: (I) a {\em high-temperature} or weak disorder phase, (II) a {\em pinned} phase for strong amplitude disorder, and (III) a {\em diffusive} phase for strong phase disorder, but weak amplitude disorder. We compute analytically the renormalized disorder correlator, equivalent to the Burgers velocity-velocity correlator at long times. In phase III, it assumes a universal form. For strong phase disorder, interference leads to a logarithmic singularity, related to zeroes of the partition sum, or poles of the complex Burgers velocity field. These results are valuable in the search for the adequate field theory for higher-dimensional systems.

preprint2011arXiv

Rings and boxes in dissipative environments

We study a particle on a ring in presence of a dissipative Caldeira-Leggett environment and derive its response to a DC field. We find, through a 2-loop renormalization group analysis, that a large dissipation parameter $η$ flows to a fixed point $η_R=η_c=\hbar/2π$. We also reexamine the mapping of this problem to that of the Coulomb box and show that the relaxation resistance, of recent interest, is quantized for large $η$. For finite $η>η_c$ we find that a certain average of the relaxation resistance is quantized. We propose a box experiment to measure a quantized noise.

preprint2011arXiv

Shock statistics in higher-dimensional Burgers turbulence

We conjecture the exact shock statistics in the inviscid decaying Burgers equation in D>1 dimensions, with a special class of correlated initial velocities, which reduce to Brownian for D=1. The prediction is based on a field-theory argument, and receives support from our numerical calculations. We find that, along any given direction, shocks sizes and locations are uncorrelated.

preprint2010arXiv

Free-energy distribution of the directed polymer at high temperature

We study the directed polymer of length $t$ in a random potential with fixed endpoints in dimension 1+1 in the continuum and on the square lattice, by analytical and numerical methods. The universal regime of high temperature $T$ is described, upon scaling 'time' $t \sim T^5/κ$ and space $x = T^3/κ$ (with $κ=T$ for the discrete model) by a continuum model with $δ$-function disorder correlation. Using the Bethe Ansatz solution for the attractive boson problem, we obtain all positive integer moments of the partition function. The lowest cumulants of the free energy are predicted at small time and found in agreement with numerics. We then obtain the exact expression at any time for the generating function of the free energy distribution, in terms of a Fredholm determinant. At large time we find that it crosses over to the Tracy Widom distribution (TW) which describes the fixed $T$ infinite $t$ limit. The exact free energy distribution is obtained for any time and compared with very recent results on growth and exclusion models.

preprint2010arXiv

Freezing Transition in Decaying Burgers Turbulence and Random Matrix Dualities

We reveal a phase transition with decreasing viscosity $ν$ at ν=ν_c>0 in one-dimensional decaying Burgers turbulence with a power-law correlated random profile of Gaussian-distributed initial velocities <v(x,0)v(x',0)>\sim|x-x'|^{-2}. The low-viscosity phase exhibits non-Gaussian one-point probability density of velocities, continuously dependent on ν, reflecting a spontaneous one step replica symmetry breaking (RSB) in the associated statistical mechanics problem. We obtain the low orders cumulants analytically. Our results, which are checked numerically, are based on combining insights in the mechanism of the freezing transition in random logarithmic potentials with an extension of duality relations discovered recently in Random Matrix Theory. They are essentially non mean-field in nature as also demonstrated by the shock size distribution computed numerically and different from the short range correlated Kida model, itself well described by a mean field one step RSB ansatz. We also provide some insights for the finite viscosity behaviour of velocities in the latter model.

preprint2010arXiv

The Aharonov-Bohm effect in presence of dissipative environments

We study a particle on a ring in presence of various dissipative environments. We develop and solve a variational scheme assuming low frequency dominance. Our solution produces a renormalization group (RG) transformation to all orders in the inverse dissipation strength, and in particular reproduces known two loop results. Our RG leads to a weak dissipation parameter, for which a weak coupling expansion for the position correlation function shows a 1/t^2 decay in imaginary time.

preprint2010arXiv

The universal high temperature regime of pinned elastic objects

We study the high temperature regime within the glass phase of an elastic object with short ranged disorder. In that regime we argue that the scaling functions of any observable are described by a continuum model with a $δ$-correlated disorder and that they are universal up to only two parameters that can be explicitly computed. This is shown numerically on the roughness of directed polymer models and by dimensional and functional renormalization group arguments. We discuss experimental consequences such as non-monotonous behaviour with temperature.

preprint2009arXiv

Avalanche-size distribution at the depinning transition: A numerical test of the theory

We calculate numerically the sizes S of jumps (avalanches) between successively pinned configurations of an elastic line (d=1) or interface (d=2), pulled by a spring of (small) strength m^2 in a random-field landscape. We obtain strong evidence that the size distribution, away from the small-scale cutoff, takes the form P(S) = p(S/S_m) <S>/S_m^2, where S_m:=<S^2>/(2<S>), proportional to m^(-d-zeta), is the scale of avalanches, and zeta the roughness exponent at the depinning transition. Measurement of the scaling function f(s) := s^tau p(s) is compared with the predictions from a recent Functional RG (FRG) calculation, both at mean-field and one-loop level. The avalanche-size exponent tau is found in good agreement with the conjecture tau = 2- 2/(d+zeta), recently confirmed to one loop via the FRG. The function f(s) exhibits a shoulder and a stretched exponential decay at large s, with ln f(s) proportional to - s^delta, and delta approximately 7/6 in d=1. The function f(s), universal ratios of moments, and the generating function <exp(lambda s)> are found in excellent agreement with the one-loop FRG predictions. The distribution of local avalanche sizes, i.e. of the jumps of a subspace of the manifold of dimension d', is also computed and compared to our FRG predictions, and to the conjecture tau' =2- 2/(d'+zeta).

preprint2009arXiv

Elasticity of a contact-line and avalanche-size distribution at depinning

Motivated by recent experiments, we extend the Joanny-deGennes calculation of the elasticity of a contact line to an arbitrary contact angle and an arbitrary plate inclination in presence of gravity. This requires a diagonalization of the elastic modes around the non-linear equilibrium profile, which is carried out exactly. We then make detailed predictions for the avalanche-size distribution at quasi-static depinning: we study how the universal (i.e. short-scale independent) rescaled size distribution and the ratio of moments of local to global avalanches depend on the precise form of the elastic kernel.

preprint2009arXiv

Extreme value statistics from the Real Space Renormalization Group: Brownian Motion, Bessel Processes and Continuous Time Random Walks

We use the Real Space Renormalization Group (RSRG) method to study extreme value statistics for a variety of Brownian motions, free or constrained such as the Brownian bridge, excursion, meander and reflected bridge, recovering some standard results, and extending others. We apply the same method to compute the distribution of extrema of Bessel processes. We briefly show how the continuous time random walk (CTRW) corresponds to a non standard fixed point of the RSRG transformation.

preprint2009arXiv

Height fluctuations of a contact line: a direct measurement of the renormalized disorder correlator

We have measured the center-of-mass fluctuations of the height of a contact line at depinning for two different systems: liquid hydrogen on a rough cesium substrate and isopropanol on a silicon wafer grafted with silanized patches. The contact line is subject to a confining quadratic well, provided by gravity. From the second cumulant of the height fluctuations, we measure the renormalized disorder correlator Delta(u), predicted by the Functional RG theory to attain a fixed point, as soon as the capillary length is large compared to the Larkin length set by the microscopic disorder. The experiments are consistent with the asymptotic form for Delta(u) predicted by Functional RG, including a linear cusp at u=0. The observed small deviations could be used as a probe of the underlying physical processes. The third moment, as well as avalanche-size distributions are measured and compared to predictions from Functional RG.

preprint2009arXiv

Sinai model in presence of dilute absorbers

We study the Sinai model for the diffusion of a particle in a one dimension random potential in presence of a small concentration $ρ$ of perfect absorbers using the asymptotically exact real space renormalization method. We compute the survival probability, the averaged diffusion front and return probability, the two particle meeting probability, the distribution of total distance traveled before absorption and the averaged Green's function of the associated Schrodinger operator. Our work confirms some recent results of Texier and Hagendorf obtained by Dyson-Schmidt methods, and extends them to other observables and in presence of a drift. In particular the power law density of states is found to hold in all cases. Irrespective of the drift, the asymptotic rescaled diffusion front of surviving particles is found to be a symmetric step distribution, uniform for $|x| < {1/2} ξ(t)$, where $ξ(t)$ is a new, survival length scale ($ξ(t)=T \ln t/\sqrtρ$ in the absence of drift). Survival outside this sharp region is found to decay with a larger exponent, continuously varying with the rescaled distance $x/ξ(t)$. A simple physical picture based on a saddle point is given, and universality is discussed.

preprint2009arXiv

Statistical Mechanics of Logarithmic REM: Duality, Freezing and Extreme Value Statistics of $1/f$ Noises generated by Gaussian Free Fields

We compute the distribution of the partition functions for a class of one-dimensional Random Energy Models (REM) with logarithmically correlated random potential, above and at the glass transition temperature. The random potential sequences represent various versions of the 1/f noise generated by sampling the two-dimensional Gaussian Free Field (2dGFF) along various planar curves. Our method extends the recent analysis of Fyodorov Bouchaud from the circular case to an interval and is based on an analytical continuation of the Selberg integral. In particular, we unveil a {\it duality relation} satisfied by the suitable generating function of free energy cumulants in the high-temperature phase. It reinforces the freezing scenario hypothesis for that generating function, from which we derive the distribution of extrema for the 2dGFF on the $[0,1]$ interval. We provide numerical checks of the circular and the interval case and discuss universality and various extensions. Relevance to the distribution of length of a segment in Liouville quantum gravity is noted.

preprint2009arXiv

Winding of planar gaussian processes

We consider a smooth, rotationally invariant, centered gaussian process in the plane, with arbitrary correlation matrix $C_{t t'}$. We study the winding angle $ϕ_t$ around its center. We obtain a closed formula for the variance of the winding angle as a function of the matrix $C_{tt'}$. For most stationary processes $C_{tt'}=C(t-t')$ the winding angle exhibits diffusion at large time with diffusion coefficient $D = \int_0^\infty ds C'(s)^2/(C(0)^2-C(s)^2)$. Correlations of $\exp(i n ϕ_t)$ with integer $n$, the distribution of the angular velocity $\dot ϕ_t$, and the variance of the algebraic area are also obtained. For smooth processes with stationary increments (random walks) the variance of the winding angle grows as ${1/2} (\ln t)^2$, with proper generalizations to the various classes of fractional Brownian motion. These results are tested numerically. Non integer $n$ is studied numerically.

preprint2008arXiv

Critical interfaces in the random-bond Potts model

We study geometrical properties of interfaces in the random-temperature q-states Potts model as an example of a conformal field theory weakly perturbed by quenched disorder. Using conformal perturbation theory in q-2 we compute the fractal dimension of Fortuin Kasteleyn domain walls. We also compute it numerically both via the Wolff cluster algorithm for q=3 and via transfer-matrix evaluations. We obtain numerical results for the fractal dimension of spin cluster interfaces for q=3. These are found numerically consistent with the duality kappa(spin) * kappa(FK)= 16 as expressed in putative SLE parameters.

preprint2008arXiv

Driven particle in a random landscape: disorder correlator, avalanche distribution and extreme value statistics of records

We review how the renormalized force correlator Delta(u), the function computed in the functional RG field theory, can be measured directly in numerics and experiments on the dynamics of elastic manifolds in presence of pinning disorder. We show how this function can be computed analytically for a particle dragged through a 1-dimensional random-force landscape. The limit of small velocity allows to access the critical behavior at the depinning transition. For uncorrelated forces one finds three universality classes, corresponding to the three extreme value statistics, Gumbel, Weibull, and Frechet. For each class we obtain analytically the universal function Delta(u), the corrections to the critical force, and the joint probability distribution of avalanche sizes s and waiting times w. We find P(s)=P(w) for all three cases. All results are checked numerically. For a Brownian force landscape, known as the ABBM model, avalanche distributions and Delta(u) can be computed for any velocity. For 2-dimensional disorder, we perform large-scale numerical simulations to calculate the renormalized force correlator tensor Delta_{ij}(u), and to extract the anisotropic scaling exponents zeta_x > zeta_y. We also show how the Middleton theorem is violated. Our results are relevant for the record statistics of random sequences with linear trends, as encountered e.g. in some models of global warming. We give the joint distribution of the time s between two successive records and their difference in value w.

preprint2008arXiv

Fluctuation force exerted by a planar self-avoiding polymer

Using results from Schramm Loewner evolution (SLE), we give the expression of the fluctuation-induced force exerted by a polymer on a small impenetrable disk, in various 2-dimensional domain geometries. We generalize to two polymers and examine whether the fluctuation force can trap the object into a stable equilibrium. We compute the force exerted on objects at the domain boundary, and the force mediated by the polymer between such objects. The results can straightforwardly be extended to any SLE interface, including Ising, percolation, and loop-erased random walks. Some are relevant for extremal value statistics.

preprint2008arXiv

Size distributions of shocks and static avalanches from the Functional Renormalization Group

Interfaces pinned by quenched disorder are often used to model jerky self-organized critical motion. We study static avalanches, or shocks, defined here as jumps between distinct global minima upon changing an external field. We show how the full statistics of these jumps is encoded in the functional-renormalization-group fixed-point functions. This allows us to obtain the size distribution P(S) of static avalanches in an expansion in the internal dimension d of the interface. Near and above d=4 this yields the mean-field distribution P(S) ~ S^(-3/2) exp(-S/[4 S_m]) where S_m is a large-scale cutoff, in some cases calculable. Resumming all 1-loop contributions, we find P(S) ~ S^(-tau) exp(C (S/S_m)^(1/2) -B/4 (S/S_m)^delta) where B, C, delta, tau are obtained to first order in epsilon=4-d. Our result is consistent to O(epsilon) with the relation tau = 2-2/(d+zeta), where zeta is the static roughness exponent, often conjectured to hold at depinning. Our calculation applies to all static universality classes, including random-bond, random-field and random-periodic disorder. Extended to long-range elastic systems, it yields a different size distribution for the case of contact-line elasticity, with an exponent compatible with tau=2-1/(d+zeta) to O(epsilon=2-d). We discuss consequences for avalanches at depinning and for sandpile models, relations to Burgers turbulence and the possibility that the above relations for tau be violated to higher loop order. Finally, we show that the avalanche-size distribution on a hyper-plane of co-dimension one is in mean-field (valid close to and above d=4) given by P(S) ~ K_{1/3}(S)/S, where K is the Bessel-K function, thus tau=4/3 for the hyper plane.

preprint2008arXiv

Statistics of static avalanches in a random pinning landscape

We study the minimum-energy configuration of a d-dimensional elastic interface in a random potential tied to a harmonic spring. As a function of the spring position, the center of mass of the interface changes in discrete jumps, also called shocks or "static avalanches''. We obtain analytically the distribution of avalanche sizes and its cumulants within an epsilon=4-d expansion from a tree and 1-loop resummation, using functional renormalization. This is compared with exact numerical minimizations of interface energies for random field disorder in d=2,3. Connections to the Burgers equation and to dynamic avalanches are discussed.

preprint2006arXiv

Measuring functional renormalization group fixed-point functions for pinned manifolds

Exact numerical minimization of interface energies is used to test the functional renormalization group (FRG) analysis for interfaces pinned by quenched disorder. The fixed-point function R(u) (the correlator of the coarse-grained disorder) is computed. In dimensions D=d+1, a linear cusp in R''(u) is confirmed for random bond (d=1,2,3), random field (d=0,2,3), and periodic (d=2,3) disorders. The functional shocks that lead to this cusp are seen. Small, but significant, deviations from 1-loop FRG results are compared to 2-loop corrections. The cross-correlation for two copies of disorder is compared with a recent FRG study of chaos.

preprint2006arXiv

Stability of Random-Field and Random-Anisotropy Fixed Points at large N

In this note, we clarify the stability of the large-N functional RG fixed points of the order/disorder transition in the random-field (RF) and random-anisotropy (RA) O(N) models. We carefully distinguish between infinite N, and large but finite N. For infinite N, the Schwarz-Soffer inequality does not give a useful bound, and all fixed points found in cond-mat/0510344 (Phys. Rev. Lett. 96, 197202 (2006)) correspond to physical disorder. For large but finite N (i.e. to first order in 1/N) the non-analytic RF fixed point becomes unstable, and the disorder flows to an analytic fixed point characterized by dimensional reduction. However, for random anisotropy the fixed point remains non-analytic (i.e. exhibits a cusp) and is stable in the 1/N expansion, while the corresponding dimensional-reduction fixed point is unstable. In this case the Schwarz-Soffer inequality does not constrain the 2-point spin correlation. We compute the critical exponents of this new fixed point in a series in 1/N and to 2-loop order.

preprint2004arXiv

Can non-linear elasticity explain contact-line roughness at depinning?

We examine whether cubic non-linearities, allowed by symmetry in the elastic energy of a contact line, may result in a different universality class at depinning. Standard linear elasticity predicts a roughness exponent zeta=1/3 (one loop), zeta=0.39 (numerics) while experiments give zeta of about 0.5. Within functional RG we find that a non-local KPZ-type term is generated at depinning and grows under coarse graining. A fixed point with zeta=0.45 (one loop) is identified, showing that large enough cubic terms increase the roughness. This fixed point is unstable, revealing a rough strong-coupling phase. Experimental study of contact angles theta near pi/2, where cubic terms vanish in the energy, is suggested.

preprint2004arXiv

Derivation of the Functional Renormalization Group Beta-Function at order 1/N for Manifolds Pinned by Disorder

In an earlier publication, we have introduced a method to obtain, at large N, the effective action for d-dimensional manifolds in a N-dimensional disordered environment. This allowed to obtain the Functional Renormalization Group (FRG) equation for N=infinity and was shown to reproduce, with no need for ultrametric replica symmetry breaking, the predictions of the Mezard-Parisi solution. Here we compute the corrections at order 1/N. We introduce two novel complementary methods, a diagrammatic and an algebraic one, to perform the complicated resummation of an infinite number of loops, and derive the beta-function of the theory to order 1/N. We present both the effective action and the corresponding functional renormalization group equations. The aim is to explain the conceptual basis and give a detailed account of the novel aspects of such calculations. The analysis of the FRG flow, comparison with other studies, and applications, e.g. to the strong-coupling phase of the Kardar-Parisi-Zhang equation are examined in a subsequent publication.

preprint2000arXiv

Disorder Induced Transitions in Layered Coulomb Gases and Superconductors

A 3D layered system of charges with logarithmic interaction parallel to the layers and random dipoles is studied via a novel variational method and an energy rationale which reproduce the known phase diagram for a single layer. Increasing interlayer coupling leads to successive transitions in which charge rods correlated in N>1 neighboring layers are nucleated by weaker disorder. For layered superconductors in the limit of only magnetic interlayer coupling, the method predicts and locates a disorder-induced defect-unbinding transition in the flux lattice. While N=1 charges dominate there, N>1 disorder induced defect rods are predicted for multi-layer superconductors.

preprint2000arXiv

Glass transition of a particle in a random potential, front selection in nonlinear RG and entropic phenomena in Liouville and SinhGordon models

We study via RG, numerics, exact bounds and qualitative arguments the equilibrium Gibbs measure of a particle in a $d$-dimensional gaussian random potential with {\it translationally invariant logarithmic} spatial correlations. We show that for any $d \ge 1$ it exhibits a transition at $T=T_c>0$. The low temperature glass phase has a non trivial structure, being dominated by {\it a few} distant states (with replica symmetry breaking phenomenology). In finite dimension this transition exists only in this ``marginal glass'' case (energy fluctuation exponent $θ=0$) and disappears if correlations grow faster (single ground state dominance $θ>0$) or slower (high temperature phase). The associated extremal statistics problem for correlated energy landscapes exhibits universal features which we describe using a non linear (KPP) RG equation. These include the tails of the distribution of the minimal energy (or free energy) and the finite size corrections which are universal. The glass transition is closely related to Derrida's random energy models. In $d=2$ the glass transition of the particle exhibits interesting similarities with the weak to strong coupling transition in Liouville ($c=1$ barrier) and with a transition that we conjecture for the sinh-Gordon model. Glassy freezing of the particle is associated with the generation under RG of new local operators and of non-smooth configurations in Liouville. Applications to Dirac fermions in random magnetic fields at criticality reveals a peculiar ``quasi-localized'' regime (the glass phase for the particle) where eigenfunctions are concentrated over {\it a finite number} of distant regions, and allows to recover the multifractal spectrum.

preprint1996arXiv

Dynamics of particles and manifolds in a quenched random force field

We study the dynamics of a directed manifold of internal dimension D in a d-dimensional random force field. We obtain an exact solution for $d \to \infty$ and a Hartree approximation for finite d. They yield a Flory-like roughness exponent $ζ$ and a non trivial anomalous diffusion exponent $ν$ continuously dependent on the ratio $g_{T}/g_{L}$ of divergence-free ($g_{T}$) to potential ($g_{L}$) disorder strength. For the particle (D=0) our results agree with previous order $ε^2$ RG calculations. The time-translational invariant dynamics for $g_{T} >0$ smoothly crosses over to the previously studied ultrametric aging solution in the potential case.

preprint1996arXiv

Glass phase of two-dimensional triangular elastic lattices with disorder

We study two dimensional triangular elastic lattices in a background of point disorder, excluding dislocations (tethered network). Using both (replica symmetric) static and (equilibrium) dynamic renormalization group for the corresponding $N=2$ component model, we find a transition to a glass phase for $T < T_g$, described by a plane of perturbative fixed points. The growth of displacements is found to be asymptotically isotropic with $u_T^2 \sim u_L^2 \sim A_1 \ln^2 r$, with universal subdominant anisotropy $u_T^2 - u_L^2 \sim A_2 \ln r$. where $A_1$ and $A_2$ depend continuously on temperature and the Poisson ratio $σ$. We also obtain the continuously varying dynamical exponent $z$. For the Cardy-Ostlund $N=1$ model, a particular case of the above model, we point out a discrepancy in the value of $A_1$ with other published results in the litterature. We find that our result reconciles the order of magnitude of the RG predictions with the most recent numerical simulations.