Measuring functional renormalization group fixed-point functions for pinned manifolds

Theory and experiments for disordered elastic manifolds, depinning, avalanches, and sandpiles

Domain walls in magnets, vortex lattices in superconductors, contact lines at depinning, and many other systems can be modeled as an elastic system subject to quenched disorder. The ensuing field theory possesses a well-controlled perturbative expansion around its upper critical dimension. Contrary to standard field theory, the renormalization group (RG) flow involves a function, the disorder correlator Δ(w), and is therefore termed the functional RG. Δ(w) is a physical observable, the auto-correlation function of the center of mass of the elastic manifold. In this review, we give a pedagogical introduction into its phenomenology and techniques. This allows us to treat both equilibrium (statics), and depinning (dynamics). Building on these techniques, avalanche observables are accessible: distributions of size, duration, and velocity, as well as the spatial and temporal shape. Various equivalences between disordered elastic manifolds, and sandpile models exist: an elastic string driven at a point and the Oslo model; disordered elastic manifolds and Manna sandpiles; charge density waves and Abelian sandpiles or loop-erased random walks. Each of the mappings between these systems requires specific techniques, which we develop, including modeling of discrete stochastic systems via coarse-grained stochastic equations of motion, super-symmetry techniques, and cellular automata. Stronger than quadratic nearest-neighbor interactions lead to directed percolation, and non-linear surface growth with additional Kardar-Parisi-Zhang (KPZ) terms. On the other hand, KPZ without disorder can be mapped back to disordered elastic manifolds, either on the directed polymer for its steady state, or a single particle for its decay. Other topics covered are the relation between functional RG and replica symmetry breaking, and random-field magnets. Emphasis is given to numerical and experimental tests of the theory.

Mean-field theories for depinning and their experimental signatures

Mean-field theory is an approximation replacing an extended system by a few variables. For depinning of elastic manifolds, these are the position u of its center of mass and the statistics of the forces F(u). There are two proposals how to model the latter: as a random walk (ABBM model), or as uncorrelated forces at integer u (discretized particle model, DPM). While for many experiments the ABBM model (in the literature misleadingly equated with mean-field theory) makes quantitatively correct predictions for the distributions of velocities, or avalanche size and duration, the microscopic disorder force-force correlations cannot grow linearly, and thus unboundedly as a random walk, with distance. Even the effective (renormalized) disorder forces which do so at small distances are bounded at large distances. To describe both regimes, we model forces as an Ornstein-Uhlenbeck process. The latter has the statistics of a random walk at small scales, and is uncorrelated at large scales. By connecting to results in both limits, we solve the model largely analytically, allowing us to describe in all regimes the distributions of velocity, avalanche size, and duration. To establish experimental signatures of this transition, we study the response function, and the correlation function of position u, velocity u[over ̇], and forces F under slow driving with velocity v>0. While at v=0 force or position correlations have a cusp at the origin and then decay at least exponentially fast to zero, this cusp is rounded at a finite driving velocity. We give a detailed analytic analysis for this rounding by velocity, which allows us, given experimental data, to extract the timescale of the response function, and to reconstruct the force-force correlator at v=0. The latter is the central object of the field theory, and as such contains detailed information about the universality class in question. We test our predictions by careful numerical simulations extending over up to ten orders in magnitude.

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 anticorrelated, the occurrence of a larger avalanche making more likely the occurrence of a smaller one, and vice versa. Examining the universality of our results leads us to conjecture several 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.

Glassy properties of the Bose-glass phase of a one-dimensional disordered Bose fluid

We study a one-dimensional disordered Bose fluid using bosonization, the replica method, and a nonperturbative functional renormalization-group approach. The Bose-glass phase is described by a fully attractive strong-disorder fixed point characterized by a singular disorder correlator whose functional dependence assumes a cuspy form that is related to the existence of metastable states. At nonzero momentum scale, quantum tunneling between these metastable states leads to a rounding of the nonanalyticity in a quantum boundary layer that encodes the existence of rare superfluid regions responsible for the ω^{2} behavior of the (dissipative) conductivity in the low-frequency limit. These results can be understood within the "droplet" picture put forward for the description of glassy (classical) systems.

Topological Order Generated by a Random Field in a 2D Exchange Model

We study a 2D exchange model with a weak static random field on lattices containing over 10^{8} spins. Ferromagnetic correlations persist on the Imry-Ma scale inversely proportional to the random-field strength and decay exponentially at greater distances. We find that the average energy of the correlated area is close to the ground-state energy of a Skyrmion, while the topological charge of the area is close to ±1. The correlation function of the topological charge density changes sign at a distance determined by the ferromagnetic correlation length, while its Fourier transform exhibits a maximum. These findings suggest that static randomness transforms a 2D ferromagnetic state into a Skyrmion-anti-Skyrmion glass.

Spatial shape of avalanches

In disordered elastic systems, driven by displacing a parabolic confining potential adiabatically slowly, all advance of the system is in bursts, termed avalanches. Avalanches have a finite extension in time, which is much smaller than the waiting time between them. Avalanches also have a finite extension ℓ in space, i.e., only a part of the interface of size ℓ moves during an avalanche. Here we study their spatial shape 〈S(x)〉_{ℓ} given ℓ, as well as its fluctuations encoded in the second cumulant 〈S^{2}(x)〉_{ℓ}^{c}. We establish scaling relations governing the behavior close to the boundary. We then give analytic results for the Brownian force model, in which the microscopic disorder for each degree of freedom is a random walk. Finally, we confirm these results with numerical simulations. To do this properly we elucidate the influence of discretization effects, which also confirms the assumptions entering into the scaling ansatz. This allows us to reach the scaling limit already for avalanches of moderate size. We find excellent agreement for the universal shape and its fluctuations, including all amplitudes.

Crack propagation through disordered materials as a depinning transition: A critical test of the theory

The dynamics of a planar crack propagating within a brittle disordered material is investigated numerically. The fracture front evolution is described as the depinning of an elastic line in a random field of toughness. The relevance of this approach is critically tested through the comparison of the roughness front properties, the statistics of avalanches, and the local crack velocity distribution with experimental results. Our simulations capture the main features of the fracture front evolution as measured experimentally. However, some experimental observations such as the velocity distribution are not consistent with the behavior of an elastic line close to the depinning transition. This discrepancy suggests the presence of another failure mechanism not included in our model of brittle failure.

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 〈S_{1}S_{2}〉^{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 nontrivial 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.

Dynamical selection of critical exponents

In renormalized field theories there are in general one or few fixed points that are accessible by the renormalization-group flow. They can be identified from the fixed-point equations. Exceptionally, an infinite family of fixed points exists, parameterized by a scaling exponent ζ, itself a function of a nonrenormalizing parameter. Here we report a different scenario with an infinite family of fixed points of which seemingly only one is chosen by the renormalization-group flow. This dynamical selection takes place in systems with an attractive interaction V(ϕ), as in standard ϕ^{4} theory, but where the potential V at large ϕ goes to zero, as, e.g., the attraction by a defect.

Depinning of stiff directed lines in random media

Driven elastic manifolds in random media exhibit a depinning transition to a state with nonvanishing velocity at a critical driving force. We study the depinning of stiff directed lines, which are governed by a bending rigidity rather than line tension. Their equation of motion is the (quenched) Herring-Mullins equation, which also describes surface growth governed by surface diffusion. Stiff directed lines are particularly interesting as there is a localization transition in the static problem at a finite temperature and the commonly exploited time ordering of states by means of Middleton's theorems [Phys. Rev. Lett. 68, 670 (1992)] is not applicable. We employ analytical arguments and numerical simulations to determine the critical exponents and compare our findings with previous works and functional renormalization group results, which we extend to the different line elasticity. We see evidence for two distinct correlation length exponents.

Avalanches and dimensional reduction breakdown in the critical behavior of disordered systems

We investigate the connection between a formal property of the critical behavior of several disordered systems, known as "dimensional reduction," and the presence in these systems at zero temperature of collective events known as "avalanches." Avalanches generically produce nonanalyticities in the functional dependence of the cumulants of the renormalized disorder. We show that this leads to a breakdown of the dimensional reduction predictions if and only if the fractal dimension characterizing the scaling properties of the avalanches is exactly equal to the difference between the dimension of space and the scaling dimension of the primary field. This is proven by combining scaling theory and the functional renormalization group. We therefore clarify the puzzle of why dimensional reduction remains valid in random field systems above a nontrivial dimension (but fails below), always applies to the statistics of branched polymer, and is always wrong in elastic models of interfaces in a random environment.

First-principles derivation of static avalanche-size distributions

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 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 Alessandro-Beatrice-Bertotti-Montorsi (ABBM) model in the nonequilibrium context of depinning. This model can itself be treated exactly in any dimension and its shock statistics is that of a Lévy 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.

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 Nguyen-Spivak-Shklovskii model of hopping conductivity in disordered insulators and the Chalker-Coddington model for the (spin) quantum Hall effect. We also propose a direct realization in an experiment with cold atoms. The model has three distinct phases: (I) a high-temperature or weak disorder phase, (II) a pinned phase for strong amplitude disorder, and (III) a 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 zeros 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.

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 =4-d expansion from a tree and one-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 dynamic avalanches are mentioned.

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) approximately S;{-3/2}e;{-S4S_{m}} , where S_{m} is a large-scale cutoff, in some cases calculable. Resumming all one-loop contributions, we find P(S) approximately S;{-tau}exp(C(SS_{m});{1/2}-B/4(S/S_{m});{delta}) , where B , C , delta , and tau are obtained to first order in =4-d . Our result is consistent to O() with the relation tau=tau_{zeta}:=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 disorders. 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(=2-d) . We discuss consequences for avalanches at depinning and for sandpile models, relations to Burgers turbulence and the possibility that the relation tau=tau_{zeta} be violated to higher loop order. Finally, we show that the avalanche-size distribution on a hyperplane of codimension one is in mean field (valid close to and above d=4 ) given by P(S) approximately K_{13}(S)S , where K is the Bessel- K function, thus tau_{hyperplane}=4/3 .

Cusps, self-organization, and absorbing states

Elastic interfaces embedded in (quenched) random media exhibit metastability and stick-slip dynamics. These nontrivial dynamical features have been shown to be associated with cusp singularities of the coarse-grained disorder correlator. Here we show that annealed systems with many absorbing states and a conservation law but no quenched disorder exhibit identical cusps. On the other hand, similar nonconserved systems in the directed percolation class are also shown to exhibit cusps but of a different type. These results are obtained both by a recent method to explicitly measure disorder correlators and by defining an alternative new protocol inspired by self-organized criticality, which opens the door to easily accessible experimental realizations.

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

We review how the renormalized force correlator Delta(micro) , the function computed in the functional renormalization-group (RG) field theory, can be measured directly in numerics and experiments on the dynamics of elastic manifolds in the presence of pinning disorder. We show how this function can be computed analytically for a particle dragged through a one-dimensional random-force landscape. The limit of small velocity allows one 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 Fréchet. For each class we obtain analytically the universal function Delta(micro) , 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 Alessandro, Beatrice, Bertotti, and Montorsi (ABBM) model, avalanche distributions and Delta(micro) can be computed for any velocity. For two-dimensional disorder, we perform large-scale numerical simulations to calculate the renormalized force correlator tensor Delta_{ij}(micro[over ]) , 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 .

Universality class of replica symmetry breaking, scaling behavior, and the low-temperature fixed-point order function of the Sherrington-Kirkpatrick model

A scaling theory of replica symmetry breaking (RSB) in the Sherrington-Kirkpatrick (SK) model is presented in the framework of critical phenomena for the scaling regime of large RSB orders kappa , small temperatures T , and small (homogeneous) magnetic fields H . We employ the pseudodynamical picture [R. Oppermann, M. J. Schmidt, and D. Sherrington, Phys. Rev. Lett. 98, 127201 (2007)], where two critical points CP1 and CP2 are associated with the order function's pseudodynamical limits lim_{a-->infinity}q(a)=1 and lim_{a-->0}q(a)=0 at (T=0 , H=0 , 1kappa=0) . CP1 - and CP2 -dominated contributions to the free energy functional F[q(a)] require an unconventional scaling hypothesis. We determine the scaling contributions in accordance with detailed numerical self-consistent solutions for up to 200 orders of RSB. Power laws, scaling functions, and crossover lines are obtained. CP1 -dominated behavior is found for the nonequilibrium susceptibility, which decays like chi_{1}=kappa;{-53}f_{1}(Tkappa;{-53}) , for the entropy, which obeys S(T=0) approximately chi_{1};{2} , and for the subclass of diverging parameters a_{i}=kappa;{53}f_{a_{i}}(Tkappa;{-53}) [describing Parisi box sizes m_{i}(T) identical witha_{i}(T)T ], with f_{1}(zeta) approximately zeta and f_{a_{i}}(zeta) approximately 1zeta for zeta-->infinity , while f(0) is finite. CP2 -dominated behavior, controlled by the magnetic field H while temperature is irrelevant, is retrieved in the plateau height (or width) of the order function q(a) according to q_{pl}(H)=kappa;{-1}f_{pl}(H;{23}kappa;{-1}) with f_{pl}mid R:(zeta)mid R:_{zeta-->infinity} approximately zeta and f_{pl}(0) finite. Divergent characteristic RSB orders kappa_{CP1}(T) approximately T;{-35} and kappa_{CP2}(H) approximately H;{-23} , respectively, describe the crossover from mean field SK- to RSB-critical behavior with rational-valued exponents extracted with high precision from our RSB data. The order function q(a) is obtained as a fixed-point function q(a) of RSB flow, in agreement with integrated fixed-point energy and susceptibility distributions.