Functional renormalization group and the field theory of disordered elastic systems

Entanglement Growth and Minimal Membranes in (d+1) Random Unitary Circuits

Understanding the nature of entanglement growth in many-body systems is one of the fundamental questions in quantum physics. Here, we study this problem by characterizing the entanglement fluctuations and distribution of a (d+1)-dimensional qubit lattice evolved under a random unitary circuit. Focusing on Clifford gates, we perform extensive numerical simulations of random circuits in 1≤d≤4 dimensions. Our findings demonstrate that properties of growth of bipartite entanglement entropy are characterized by the roughening exponents of a d-dimensional membrane in a (d+1)-dimensional elastic medium.

Experimental Test of Sinai's Model in DNA Unzipping

The experimental measurement of correlation functions and critical exponents in disordered systems is key to testing renormalization group (RG) predictions. We mechanically unzip single DNA hairpins with optical tweezers, an experimental realization of the diffusive motion of a particle in a one-dimensional random force field, known as the Sinai model. We measure the unzipping forces F_{w} as a function of the trap position w in equilibrium and calculate the force-force correlator Δ_{m}(w), its amplitude, and correlation length, finding agreement with theoretical predictions. We study the universal scaling properties since the effective trap stiffness m^{2} decreases upon unzipping. Fluctuations of the position of the base pair at the unzipping junction u scales as u∼m^{-ζ}, with a roughness exponent ζ=1.34±0.06, in agreement with the analytical prediction ζ=4/3. Our study provides a single-molecule test of the functional RG approach for disordered elastic systems in equilibrium.

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.

Dimensional reduction breakdown and correction to scaling in the random-field Ising model

We provide a theoretical analysis by means of the nonperturbative functional renormalization group (NP-FRG) of the corrections to scaling in the critical behavior of the random-field Ising model (RFIM) near the dimension d_{DR}≈5.1 that separates a region where the renormalized theory at the fixed point is supersymmetric and critical scaling satisfies the d→d-2 dimensional reduction property (d>d_{DR}) from a region where both supersymmetry and dimensional reduction break down at criticality (d<d_{DR}). We show that the NP-FRG results are in very good agreement with recent large-scale lattice simulations of the RFIM in d=5 and we detail the consequences for the leading correction-to-scaling exponent of the peculiar boundary-layer mechanism by which the dimensional-reduction fixed point disappears and the dimensional-reduction-broken fixed point emerges in d_{DR}.

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.

Bose-glass phase of a one-dimensional disordered Bose fluid: Metastable states, quantum tunneling, and droplets

We study a one-dimensional disordered Bose fluid using bosonization, the replica method, and a nonperturbative functional renormalization-group approach. We find that 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 k, quantum tunneling between the ground state and low-lying metastable states leads to a rounding of the cusp singularity into a quantum boundary layer (QBL). The width of the QBL depends on an effective Luttinger parameter K_{k}∼k^{θ} that vanishes with an exponent θ=z-1 related to the dynamical critical exponent z. The QBL encodes the existence of rare "superfluid" regions, controls the low-energy dynamics, and yields a (dissipative) conductivity vanishing as ω^{2} in the low-frequency limit. These results reveal the glassy properties (pinning, "shocks," or static avalanches) of the Bose-glass phase and can be understood within the "droplet" picture put forward for the description of glassy (classical) systems.

Manifolds in a 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≫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≥1 they admit a finite "massless" limit μ=0 which is used to provide an upper bound for the depinning threshold under an applied force.

Depinning Transition of Charge-Density Waves: Mapping onto O(n) Symmetric ϕ^{4} Theory with n→-2 and Loop-Erased Random Walks

Driven periodic elastic systems such as charge-density waves (CDWs) pinned by impurities show a nontrivial, glassy dynamical critical behavior. Their proper theoretical description requires the functional renormalization group. We show that their critical behavior close to the depinning transition is related to a much simpler model, O(n) symmetric ϕ^{4} theory in the unusual limit of n→-2. We demonstrate that both theories yield identical results to four-loop order and give both a perturbative and a nonperturbative proof of their equivalence. As we show, both theories can be used to describe loop-erased random walks (LERWs), the trace of a random walk where loops are erased as soon as they are formed. Remarkably, two famous models of non-self-intersecting random walks, self-avoiding walks and LERWs, can both be mapped onto ϕ^{4} theory, taken with formally n=0 and n→-2 components. This mapping allows us to compute the dynamic critical exponent of CDWs at the depinning transition and the fractal dimension of LERWs in d=3 with unprecedented accuracy, z(d=3)=1.6243±0.001, in excellent agreement with the estimate z=1.62400±0.00005 of 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.

Disorder-Driven Quantum Transition in Relativistic Semimetals: Functional Renormalization via the Porous Medium Equation

In the presence of randomness, a relativistic semimetal undergoes a quantum transition towards a diffusive phase. A standard approach relates this transition to the U(N) Gross-Neveu model in the limit of N→0. We show that the corresponding fixed point is infinitely unstable, demonstrating the necessity to include fluctuations beyond the usual Gaussian approximation. We develop a functional renormalization group method amenable to include these effects and show that the disorder distribution renormalizes following the so-called porous medium equation. We find that the transition is controlled by a nonanalytic fixed point drastically different from that of the U(N) Gross-Neveu model. Our approach provides a unique mechanism of spontaneous generation of a finite density of states and also characterizes the scaling behavior of the broad distribution of fluctuations close to the transition. It can be applied to other problems where nonanalytic effects may play a role, such as the Anderson localization transition.

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.

Different universality classes at the yielding transition of amorphous systems

We study the yielding transition of a two-dimensional amorphous system under shear by using a mesoscopic elasto-plastic model. The model combines a full (tensorial) description of the elastic interactions in the system and the possibility of structural reaccommodations that are responsible for the plastic behavior. The possible structural reaccommodations are encoded in the form of a "plastic disorder" potential, which is chosen independently at each position of the sample to account for local heterogeneities. We observe that the stress must exceed a critical value σ_{c} in order for the system to yield. In addition, when the system yields a flow curve (relating stress σ and strain rate γ[over ̇]) of the form γ[over ̇]∼(σ-σ_{c})^{β} is obtained. Remarkably, we observe the value of β to depend on some details of the plastic disorder potential. For smooth potentials a value of β≃2.0 is obtained, whereas for potentials obtained as a concatenation of smooth pieces a value β≃1.5 is observed in the simulations. This indicates a dependence of critical behavior on details of the plastic behavior. In addition, by integrating out nonessential, harmonic degrees of freedom, we derive a simplified scalar version of the model that represents a collection of interacting Prandtl-Tomlinson particles. A mean-field treatment of this interaction reproduces the difference of β exponents for the two classes of plastic disorder potentials and provides values of β that compare favorably with those found in the full simulations.

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.

Spinodals with Disorder: From Avalanches in Random Magnets to Glassy Dynamics

We revisit the phenomenon of spinodals in the presence of quenched disorder and develop a complete theory for it. We focus on the spinodal of an Ising model in a quenched random field (RFIM), which has applications in many areas from materials to social science. By working at zero temperature in the quasistatically driven RFIM, thermal fluctuations are eliminated and one can give a rigorous content to the notion of spinodal. We show that the latter is due to the depinning and the subsequent expansion of rare droplets. We work out the associated critical behavior, which, in any finite dimension, is very different from the mean-field one: the characteristic length diverges exponentially and the thermodynamic quantities display very mild nonanalyticities much like in a Griffith phenomenon. From the recently established connection between the spinodal of the RFIM and glassy dynamics, our results also allow us to conclusively assess the physical content and the status of the dynamical transition predicted by the mean-field theory of glass-forming liquids.

Ground-state statistics of directed polymers with heavy-tailed disorder

In this mostly numerical study, we reconsider 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 Kardar-Parisi-Zhang scaling behavior of the directed polymer, and new exponents prevail. After confirming again that the Flory argument accurately predicts these exponents 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.

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.

Finite-scale singularity in the renormalization group flow of a reaction-diffusion system

We study the nonequilibrium critical behavior of the pair contact process with diffusion (PCPD) by means of nonperturbative functional renormalization group techniques. We show that usual perturbation theory fails because the effective potential develops a nonanalyticity at a finite length scale: Perturbatively forbidden terms are dynamically generated and the flow can be continued once they are taken into account. Our results suggest that the critical behavior of PCPD can be either in the directed percolation or in a different (conjugated) universality class.

Avalanche dynamics of elastic interfaces

Slowly driven elastic interfaces, such as domain walls in dirty magnets, contact lines wetting a nonhomogeneous substrate, or cracks in brittle disordered material proceed via intermittent motion, called avalanches. Here 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 [Alessandro, Beatrice, Bertotti, and Montorsi (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 (nonlinear) instanton equation. We focus on the limit of slow uniform driving. This tree approximation is the mean-field theory (MFT) for realistic interfaces in short-ranged disorder, up to the renormalization of two parameters at d=d(uc). We calculate a number of observables of direct experimental interest: Both for the center of mass, and for a given Fourier mode q, we obtain various correlations and probability distribution functions (PDF's) of the velocity inside an avalanche, as well as the avalanche shape and its fluctuations (second shape). Within MFT we find that velocity correlations at nonzero q are asymmetric under time reversal. Next we calculate, beyond MFT, i.e., including loop corrections, the one-time PDF of the center-of-mass velocity u[over ·] for dimension d<d(uc). The singularity at small velocity P(u[over ·])~1/u[over ·](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.

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.

Random-field and random-anisotropy O(N) spin systems with a free surface

We study the surface scaling behavior of a semi-infinite d-dimensional O(N) spin system in the presence of a quenched random field and random anisotropy disorders. It is known that above the lower critical dimension d(LC) = 4 the infinite models undergo a paramagnetic-ferromagnetic transition for N > N(c) (N(c) = 2.835 for the random field and N(c) =9.441 for random anisotropy). For N < N(c) and d < d(LC) there exists a quasi-long-range-order phase with a zero order parameter and a power-law decay of spin correlations. Using a functional renormalization group, we derive the surface scaling laws that describe the ordinary surface transition for d > d(LC) and the long-range behavior of spin correlations near the surface in the quasi-long-range-order phase for d < d(LC). The corresponding surface exponents are calculated to one-loop order. The obtained results can be applied to the surface scaling of periodic elastic systems in disordered media, amorphous magnets, and (3)He-A in aerogel.

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.

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 nonlinear equilibrium profile, which is carried out exactly. We then make detailed predictions for the avalanche-size distribution at quasistatic 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.

Spontaneous versus explicit replica symmetry breaking in the theory of disordered systems

We investigate the relation between spontaneous and explicit replica symmetry breaking in the theory of disordered systems. On general ground, we prove the equivalence between the replicon operator associated with the stability of the replica-symmetric solution in the standard replica scheme and the operator signaling a breakdown of the solution with analytic field dependence in a scheme in which replica symmetry is explicitly broken by applied sources. This opens the possibility to study, via the recently developed functional renormalization group, unresolved questions related to spontaneous replica symmetry breaking and spin-glass behavior in finite-dimensional disordered systems.

Perturbation theory for ac-driven interfaces in random media

We study D-dimensional elastic manifolds driven by ac forces in a disordered environment using a perturbation expansion in the disorder strength and the mean-field approximation. We find that for D ≤ 4 , perturbation theory produces nonregular terms that grow unboundedly in time. The origin of these nonregular terms is explained. By using a graphical representation we argue that the perturbation expansion is regular to all orders for D > 4. Moreover, for the corresponding mean-field problem we prove that ill-behaved diagrams can be resummed in a way that their unbounded parts mutually cancel. Our analytical results are supported by numerical investigations. Furthermore, we conjecture the scaling of the Fourier coefficients of the mean velocity with the amplitude of the driving force h .

Depinning transition in the failure of inhomogeneous brittle materials

The dynamics of cracks propagating in elastic inhomogeneous materials is investigated experimentally. The variations of the average crack velocity with the external driving force are measured for a brittle rock and shown to display two distinct regimes: an exponential law characteristic of subcritical propagation at a low driving force and a power law above a critical threshold. This behavior can be explained quantitatively by extending linear elastic fracture mechanics to disordered systems. In this description, the motion of a crack is analogous to the one of an elastic line driven in a random medium, and critical failure occurs when the external force is sufficiently large to depin the crack front from the heterogeneities of the material.

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 .

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 .

Creep of a fracture line in paper peeling

The slow motion of a crack line is studied via an experiment in which sheets of paper are split into two halves in a "peel-in-nip" (PIN) geometry under a constant load, in creep. The velocity-force relation is exponential. The dynamics of the fracture line exhibits intermittency, or avalanches, which are studied using acoustic emission. The energy statistics is a power law, with the exponent beta ~ 1.8 +/- 0.1. Both the waiting times between subsequent events and the displacement of the fracture line imply complicated stick-slip dynamics. We discuss the correspondence to tensile PIN tests and other similar experiments on in-plane fracture and the theory of creep for elastic manifolds.

Measuring functional renormalization group fixed-point functions for pinned manifolds

Exact numerical minimization of interface energies is used to test the functional renormalization group 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 the 1-loop calculation are compared to 2-loop corrections and chaos is measured.

Statics and dynamics of elastic manifolds in media with long-range correlated disorder

We study the statics and dynamics of an elastic manifold in a disordered medium with quenched defects correlated as approximately r{-a} for large separation r. We derive the functional renormalization-group equations to one-loop order, which allow us to describe the universal properties of the system in equilibrium and at the depinning transition. Using a double epsilon=4-d and delta=4-a expansion we compute the fixed points characterizing different universality classes and analyze their regions of stability. The long-range disorder-correlator remains analytic but generates short-range disorder whose correlator exhibits the usual cusp. The critical exponents and universal amplitudes are computed to first order in epsilon and delta at the fixed points. At depinning, a velocity-versus-force exponent beta larger than unity can occur. We discuss possible realizations using extended defects.

Universal distribution of threshold forces at the depinning transition

We study the distribution of threshold forces at the depinning transition for an elastic system of finite size, driven by an external force in a disordered medium at zero temperature. Using the functional renormalization group technique, we compute the distribution of pinning forces in the quasistatic limit. This distribution is universal up to two parameters, the average critical force and its width. We discuss possible definitions for threshold forces in finite-size samples. We show how our results compare to the distribution of the latter computed recently within a numerical simulation of the so-called critical configuration.

Chaos and residual correlations in pinned disordered systems

We study, using functional renormalization, two copies of an elastic system pinned by mutually correlated random potentials. Short scale decorrelation depends on a nontrivial boundary layer regime with (possibly multiple) chaos exponents. Large scale mutual displacement correlations behave as [x - x'](2zeta-mu), mu proportional to the difference between Flory (or mean field) and exact roughness exponents zeta. For short range disorder mu>0 and small; e.g., for random bond interfaces mu=5zeta-epsilon, epsilon=4-d, and mu=epsilon{[(2pi)(2)/36]-1} for the one component Bragg glass. Random field (i.e., long range) disorder exhibits finite residual correlations (no chaos mu=0) described by new functional renormalization fixed points. Temperature and dynamic chaos (depinning) are discussed.

Random-field spin models beyond 1 loop: a mechanism for decreasing the lower critical dimension

The functional renormalization group for the random-field and random-anisotropy O(N) sigma models is studied to 2 loop. The ferromagnetic-disordered (F-D) transition fixed point is found to next order in d = 4 + epsilon for N > N(c) (N(c) = 2.834 740 8 for random field, N(c) = 9.441 21 for random anisotropy). For N < N(c) the lower critical dimension d = d(lc) plunges below d(lc) = 4: we find two fixed points, one describing the quasiordered phase, the other is novel and describes the F-D transition. d(lc) can be obtained in an (N(c)-N) expansion. The theory is also analyzed at large N and a glassy regime is found.

Unified picture of ferromagnetism, quasi-long-range order, and criticality in random-field models

By applying the recently developed nonperturbative functional renormalization group (FRG) approach, we study the interplay between ferromagnetism, quasi-long-range order (QLRO), and criticality in the d-dimensional random-field O(N) model in the whole (N, d) diagram. Even though the "dimensional reduction" property breaks down below some critical line, the topology of the phase diagram is found similar to that of the pure O(N) model, with, however, no equivalent of the Kosterlitz-Thouless transition. In addition, we obtain that QLRO, namely, a topologically ordered "Bragg glass" phase, is absent in the 3-dimensional random-field XY model. The nonperturbative results are supplemented by a perturbative FRG analysis to two loops around d = 4.

Two-loop functional renormalization for elastic manifolds pinned by disorder in N dimensions

We study elastic manifolds in an N -dimensional random potential using a functional renormalization group. We extend to N>1 our previous construction of a field theory renormalizable to two loops. For isotropic disorder with O (N) symmetry we obtain the fixed point and roughness exponent to next order in epsilon=4-d , where d is the internal dimension of the manifold. Extrapolation to the directed polymer limit d=1 allows some handle on the strong coupling phase of the equivalent N -dimensional Kardar-Parisi-Zhang growth equation, and eventually suggests an upper critical dimension d(u) approximately 2.5.