# Full text of "Higher correlations, universal distributions and finite size scaling in the field theory of depinning"

## See other formats

Higher correlations, universal distributions and finite size scaling in the field theory of depinning Pierre Le Doussal^ and Kay Jorg Wiese^ m o o (N m (N > in o m o O o CNRS-Laboratoire de Physique Theorique de I 'Ecole Normale Superieure, 24 rue Lhomond, 75005 Paris, France. ^ KITP University of California at Santa Barbara, Santa Barbara, CA 93106-4030, USA (Dated: February 2, 2008) Recently we constructed a renormalizable field theory up to two loops for the quasi-static depinning of elastic manifolds in a disordered environment. Here we explore further properties of the theory. We show how higher correlation functions of the displacement field can be computed. Drastic simplifications occur, unveiling much simpler diagrammatic rules than anticipated. This is applied to the universal scaled width-distribution. The expansion in d = 4 — e predicts that the scaled distribution coincides to the lowest orders with the one for a Gaussian theory with propagator G{q) = l/g'*^^^, ^ being the roughness exponent. The deviations from this Gaussian result are small and involve higher correlation functions, which are computed here for different boundary conditions. Other universal quantities are defined and evaluated: We perform a general analysis of the stability of the fixed point. We find that the correction-to-scaling exponent is uj = —e and not — e/3 as used in the analysis of some simulations. A more detailed study of the upper critical dimension is given, where the roughness of interfaces grows as a power of a logarithm instead of a pure power. I. INTRODUCTION Understanding the behavior of an elastic interface in a ran- dom potential is important for many experimental systems and still offers a considerable theoretical challenge [1^]. It is ex- pected that below the upper critical dimension due the inter- face is pinned by arbitrarily weak disorder, into some rough configurations and that at zero temperature it can acquire a non-zero velocity under an applied force / only if / is larger than the depinning threshold fc- A functional renormalization group (FRG) method predicts that due — 4 for the statics [5], and for the simplest universality class, the so called isotropic depinning [6-8]. There has been recent progress towards a precise descrip- tion of the depinning transition. From the theory side, the FRG for single component manifolds, originally studied to one loop in an expansion in e = due — d, has now been ex- tended to a field theory shown to be renormalizable to two loops. Renormalizable, we recall, means it has a well defined continuum limit, which is independent of all microscopic de- tails, and thus ensures universality of large scale observables. Presumably there exists a fully renormalizable theory to all orders, with full predictive power [9-11]. From the side of numerics a novel powerful algorithm allows to obtain the con- figurations at (or just below) depinning with much improved accuracy [12-14]. A reasonable agreement between the two methods was found in a measurement of the roughness expo- nent (, especially the clear conclusion that ( > e/3 contrarily to a previous conjecture [7, 8] (C = e/3) based on the 1-loop analysis. The field theory of depinning in its present form is uncon- ventional in that one must work with a non-analytic action. This peculiar feature is a deep part of the physics of the prob- lem and necessary to avoid the so called dimensional reduc- tion. It makes the perturbation theory superficially "ambigu- ous". A non-trivial step taken in [9-11] to define the theory at T = as the limit u — > of the moving phase, was to assume that the interface-position is monotonic in time. This removes the ambiguity and, remarkably, leads to a renormaliz- able theory, to at least two loops [9-11]. This is supported by the "non-crossing theorems" which apply to single component depinning and, remarkably, is the same property allowing to show ergodicity and to construct an efficient algorithm to find the exact critical configuration at depinning [12, 15]. The ori- gin of recent progresses in both numerics and field theory are thus related. Clearly one would like to test this novel field the- ory by calculating more universal measurable quantities and study its properties. In this paper we further explore the field theory constructed in [9-11]. We study displacement correlations of more than two points. We find that these correlations are static. Al- though physically natural, if one wants quasi-static depinning to make sense, this manifests itself through rather non-trivial massive cancellations in the time dependence of multi-point diagrams. We elucidate these cancellations and obtain as a consequence for a large class of diagrams much simpler di- agrammatic rules than previously anticipated. Basically, all time integrals become almost trivial, resulting in a theory with "quasi static" diagrams. We then apply these properties to the calculation of universal observables. One natural universal quantity is the so-called width distribution of the interface. Interestingly, to the two lowest leading orders in e — 4 — d, the distribution coincides with the one for a Gaussian theory with the full non-trivial propagator G{q) = 1/q'^^'^'^, C be- ing the depinning exponent. This is also the subject of a re- lated publication [16], where the distribution is also measured numerically. Here we give a detailed presentation and also compute the higher connected cumulants of the displacement field, i.e. deviations from the Gaussian. Some of these results are quoted in [16]. In a second part we study the theory at the upper critical dimension. The motivation is that no exact result is avail- able to confirm that due = 4 (the only exactly solved limit corresponding to fully connected models [17-19]). Thus the 2 question of what is the upper critical dimension d^c is still de- bated, even though the field theory of depinning [6-11, 20] clearly predicts due = 4. Also, in the other class of depin- ning transitions, the so called anisotropic depinning class with KPZ nonlinearities, there is not even a convincing prediction for due [21-23] and recent numerical studies have reopened the debate [24]. Recently it has become possible to study nu- merically depinning and statics in high dimensional spaces for reasonable system sizes with better precision, allowing for the hope to settle the issue of the upper critical dimension in the near future [12-14, 25]. It is thus important to give precise predictions for the behavior predicted by the FRG, in order to compare with numerics. Finally, we also clarify the issue of finite size scaling. In a previous work, used in several simulations, the value w = — e/3 was used for the finite size scaling exponent [26, 27]. We find that the correct value is instead w = —e. This may prove useful in numerical studies [35]. The paper is organized as follows. In Section 11 we define the model, briefly review the FRG method and field theory and define the main observable of interest here, the width dis- tribution. In Section III we compute the Laplace transform of the width distribution in perturbation theory and find that to lowest order in e it coincides with a Gaussian Approxima- tion. This approximation is introduced and further studied. Some results on Laplace Inversion are given. In Section IV we go beyond the Gaussian Approximation and compute higher connected cumulants. The detailed calculation of the fourth cumulant (4-point connected correlation function of the dis- placement field) is given, and the cancellations that occur in the field theory are studied. In Section VI we discuss the up- per critical dimension and in Section VII the finite size scal- ing. The effect of various boundary conditions is studied in Appendix A. II. MODEL AND OBSERVABLES A. Model We study the over-damped dynamics described by the equa- tion of motion vdtUxt = cVluxt + F{x, Ua;t) + f (2.1) with friction rj. Long range elasticity relevant for solid friction at the upper critical dimension d = 2, can be studied replac- ing cq"^ c\q\. In presence of an applied force / the center of mass velocity is f = L~'^ dtUxt- The pinning force is x) = —duV{u, x) and thus the second cumulant of the force is F(x, u)F{x', u') = A{u - u')6'^{x - x') , (2.2) such that A(u) = —R"{u) in the bare model, where R{u) is the correlator of the random potential. Random bond disor- der is modeled by a short range function R{u), random field (RF) disorder of amplitude a by R{v) ^ — cjul at large u and CDW disorder by a periodic function R{u). B. Review of FRG and field theory Let us briefly review the field theoretic approach, more details can be found in [10]. The dynamical action (MSR) averaged over disorder is given by e"*^ with S[u,u]= Uxt{dt-yDuxt-7i I UxtA{Uxt-Uxt')Uxt' Jxt ^ Jxtt' (2.3) Here and below we denote ■= J d'^x, in Fourier J^. := / and = / dt. The FRG shows that the full func- tion A(u) becomes relevant below d = due = 4 and a flow equation for its scale dependence has been derived to one and two loops, in an expansion in d = 4 — e. In Ref. [10] this was derived by adding a small mass m as an infrared cutoff and computing the flow of disorder, defined from the effective action r[u, ii] of the theory, as m decreases towards zero. As in [10] we will denote by Ao(m) the bare disorder correlator, i.e. the one appearing in the action S in (2.3), and by A(u) the renormaUzed one, appearing in T[u, u] which has a simi- lar expression as (2.3). The rescaled disorder is then defined by A(u) = ^rn^"^^A(um'^) , eh (2.4) where Ii = + m^)^^ is the 1-loop integral. It was then shown in [9, 10] that (2.3) leads to a functional renormalization group equation -m9„A(u) = (e - 2C)A(«) + CuK'{u) 1 {Mu)-m)\ {K{u)-m)K\uf (2.5) up to O(A^) terms. This equation implies that there are only two main universality classes at depinning, a single RF fixed point for interfaces and a periodic one for CDW type disor- der [9, 10]. Both ^ and the fixed point function A*{u) were determined to order O(e^) for these classes [9, 10]. The important feature of the field theory of depinning is that A(?i) has a cusp-like non-analyticity at u = 0. As was shown in [9, 10] calculations in the non-analytic theory (e.g. yielding (2.5)) are meaningfully performed using the expansion: A(?i) = A(0) + A(0+)|u| + ^A{0+)u^ (2.6) Performing Wick averages yields the usual diagrams, except that their actual values involve averages of e.g. sign func- tions of the fields. Replacing everywhere sgn(ut — Uf) sgn{t — t') is justified for single component quasi-static de- pinning (i.e. in the limit of vanishing velocity v = 0+). This yields diagrams with sometimes complicated internal time and momentum dependences. We find however that in some cases massive cancellations occur despite the complications due to the time dependence between various diagrams, con- tributing to the same observable. 3 C. Universal distributions and observables To motivate the present study let us consider one specific ex- ample of a universal observable, the width distribution of the configuration at depiiming (the so-called critical configura- tion). The width of a configuration is defined in a given disor- der realization as w uf (2.7) where u= J^Ux is the center of mass and L"^ the volume of the system. The basic observation is that the sample to sample probability distribution P{w'^) of is expected to be universal, with a single scale set by the disorder averaged second cumulant w"^, i.e: 1 P{w') = =f = w (2.8) f{z) is a universal function. This holds for thermal averages in a number of finite temperature problems of pure systems [28, 29]. Here we show that it also holds for depinning at T = and compute the distribution, first within a simple Gaussian approximation and second within the e-expansion. In the process we study higher point correlation functions in the field theory of depinning, define specific universal ratios of these, describing deviations from Gaussian behavior and compute them. Before turning to actual calculations let us first summarize the general spirit of the method and discuss the question of the universality of such observables. The hallmark of a renor- malizable theory is that if one expresses the correlation func- tions in an expansion in the renormalized disorder A, then the resulting expressions are UV finite, equivalently they have a well defined continuum limit, independent of short scale details. On a technical level, this can be achieved by com- puting correlations in standard perturbation theory to a given order in powers of Aq, and then using the relation between renormalized disorder A and bare one Aq to the same order, or equivalently through the definition of appropriate counter- terms [36]. Here, we restrict ourselves to calculations at dom- inant order in e and thus using either A or Aq makes no dif- ference. Beyond the Larkin scale, however, these are non- analytic functions, which is crucial. In the limit of large scales or large system sizes, the fixed point form reached by the rescaled A implies that the resulting observable, e.g. the width distribution, is universal. Universal means that these quantities do not depend on the short scale details. However they do depend on the details of the large scale infrared (IR) cutoff, i.e. of the type of chosen boundary conditions. Here we focus on periodic boundary conditions [34], also of interest for numerical simulations [16], although we sometimes give results for other types, for instance for the massive IR cutoff described in the previous paragraph. Since the FRG method developed in [9, 10] and summa- rized above uses a mass as IR-cutoff and defines disorder ver- tices at zero momentum, one should be careful in calculations with e.g. periodic boundary conditions. Since we only com- pute observables either to dominant order in e, or within a 1- loop approximation, it is easy to make the necessary replace- ments, as will be indicated below. For instance the 1-loop FRG equation remains identical to the two first fines of (2.5), the only changes being that 1 . — m9m A has to be replaced by A. 2. m ^ 1/ L'm the definition of the rescaled disorder 3. the 1-loop integral I\ = ^ entering into the k (fe^+m^) definition of the rescaled disorder has to be replaced by its homologue for periodic boundary conditions [37]: used below. (2.9) III. WIDTH DISTRIBUTION: PERTURBATION THEORY AND GAUSSIAN APPROXIMATION Let us start by giving the simplest approximation for this dis- tribution. It can be derived in two ways: (i) perturbation the- ory in the renormalized theory to lowest order in e (ii) a simple Gaussian approximation. In the end this will motivate going further, i.e. studying deviations from the Gaussian approxima- tion. A. Perturbation tlieory We now study perturbation theory. To compute the width dis- tribution using the dynamical field theoretic method [9, 10] one can start from the Laplace transform W{\) (3.1) with w"^ = Ylixi^x — uf. Here and below we omit the global multiplicative factor L~'^ in the definition of w"^, since in the end we will always normalize the distribution P{w^) by fix- ing its first moment to unity (in (3.1) it can be absorbed by a rescaling of A). Expanding in powers of the correlator of the pirming force A(u) (to lowest order this is equivalent to Aq see above), one finds that to leading order In W{X) is the sum of all connected 1-loop diagrams, as represented in Fig. 1. The loop with n disorder vertices and n insertions of w"^ is 2n E -2AA(0) \" {q^r ) (3.2) Here and below the sums over q thus runs over a d- dimensional hypercubic lattice with spacing and the 0- mode is excluded, as appropriate for periodic boundary con- ditions. If one were now to resum (3.2) over n one would find: W{X) = l[il + 2XGiq)y -1/2 (3.3) 4 r.h.s. of (3.5) generalizes into — □ — e^f FIG. 1: Examples of contributions to (3.1), terms , (w^)"^ and (w^)^ (bottom), together with the vertices for disorder (top left), (top center), and response function (top right). with G{q) = GLarkinC?) = A(0)/g^, a Larkin model type result if interpreted as naive perturbation theory (i.e. if A(0) was interpreted as the bare original disorder rather than the renormalized one). The correct procedure implies that A(0) is the running renormalized disorder A(0) A((0) — (e/i)-ie(2«-<^)'A*(0) where A*(0) is the (non-universal) value of the fixed point [9, 10]. Although I = ln(L) for the zero momentum disorder vertex, one notes that a momentum q flows in each vertex and one should take care of this by setting I ln(l/g). This yields finally (3.3) with Gig) = C/q- (3.4) where the value of C is non-universal and fixed by w"^ . As wiU become clear in the following section, the appropriate choice for G{q) is the 2-point finite-size scaling-function G^iq) = C/q''-+'^ig[qL) with g{0) = 1. The difference between the two above-mentioned choices for / simply amounts to the two different limits of small, or large qL. However to lowest order in e = 4— d they are identical. (For a calculation of the scahng function to next order in e see Appendix J in Ref. [10].) B. Gaussian Approximation and beyond A more general approach consistent with the previous calcu- lation is the following. We first note that the above result (3.3) would be exact if the distribution of the displacement fields u were Gaussian. It can thus be caUed the Gaussian Approxi- mation (GA). To understand why it was obtained here let us consider simply the second connected cumulant of the WD c 2 (w^)^ = — (u>2) . This cumulant however is not con- nected w.r.t. the u, and thus there is an exact relation: (w;2)2 Jxy (3.5) where here Gxy = u^Uy is the exact disorder averaged 2-point function. The first term is just Wick's theorem and would be the full result if the measure of the u were Gaussian. Analo- gous formulae exist for higher cumulants: the first term on the IGA 2"-i(n- 1 ./xi,..a GxiX2^X2X3 . . Gx (3.6) as a simple consequence of Wick's theorem. This is again easily resummed into (3.3) which would thus be exact if the measure of u were exactly Gaussian. Note that all results of [28, 29] for pure Gaussian theories can also easily be obtained by the present resummation method (temperature replacing disorder). Thus in the GA, the G{q) appearing in (3.3) is the exact 2-point function. It can be tested in a simulation by in- serting the measured 2-point function in (3.3). In the large but finite system size hmit it takes the scahng function form Gl((j) discussed above. When compiulng to the numerical results for the width dis- tribution, it turns out that the GA is a surprisingly good ap- proximation even down to c? = 1. This is discussed in details in [16]. However, we do not expect the GA to be exact. It is thus interesting to compute the deviation D = j'^^ ^^ w^m^ for the second cumulant of in (3.5). It is computed below and found to be of order while the GA contribution (first term in (3.5)) is of order (since G ^ e). Similarly the GA contri- bution to (3.5) is 0(e"), while the deviations are found to be 0(e^"). This can be summarized asu — ^/euo+eui where uq is a Gaussian random variable of 0(1) and ui a non-Gaussian oneof 0(1). Computing deviations from the GA is thus one motivation to compute higher point correlations. C. Laplace inversion Before doing so, let us discuss how the distribution P{w^) is obtained through an inverse Laplace-transform, as (3.7) Noting that in d = 1 (3.3) can also be written as W{X) = l[{l + 2XG{q))-\ (3.8) 9>0 this is equivalent to P{u p>0 Q>0,q¥=P G{p)J (3.9) This formula shows that for large w^, the distribution is dom- inated by the first term p = 1, and in practice summing the first few terms gives an excellent approximation. It is instructive to apply (3.9) to a random walk, where G{q) = Xjq^. Using that (n G N) ^ sin(7rx) n>0 n 1-^. one finds in terms of the width (3.10) (3.11) n>0 5 For d> 1 the situation is more complicated. Writing P{w^) = / ^ e-'^ n (1 + 2AG(g))-i . (3.12) we have e.g. at least multiphcity 2d for each factor in (3.3), as long as no component vanishes, but this multiphcity may even be higher, as can be seen from the following solutions of the diophantic equation (for 2 dimensions) 1^ + 7^ =5^ + 5^, 6^ + 7^ = 9^ + 2^, and many more. Let us define the class C{q) as peCiq) if p^^q^. (3.13) Let us index these classes by a, and introduce an order Ca<Ca' if 9 e Ca and q' € Ca' \q\ < \q'\ (3.14) The number of elements of each class is defined as \Ca I := number of elements in Ca (3.15) We further define qa ■= any element out of Ca ■ (3.16) Note that since forp e Ca and p ^ (by definition we exclude p = 0) \Ca\ is always even. (3.12) can then be rewritten as P{w') = / ^ e-^^ 11(1 + 2AG(g<,))-|C«l/2 . (3.17) a There are poles at A = —[2G{qa)]~^ ■ The sum over these poles can after partial integration be written as (|C„|/2-l)! \2G{qa) |Ca|/2 d_ |C„|/2-l H (l + 2AG(g„))- \Cc\/2 (3.18) A=-[2G(gc)]- IV. HIGHER POINT CORRELATIONS IN DEPINNING FIELD THEORY In this section we analyze how one can compute higher cor- relations in the depinning field theory and obtain simple dia- grammatic rules for doing so. These are illustrated on the 4- point function. Specific calculations and results will be given in the next Section. A. Preliminaries We want to compute at T = 0, using the dynamical action S in (2.3) , e.g. the 4-point expectation value, connected w.r.t. disorder (and u) as defined in the previous sections: / ultult' = / (ultult)^ . (4.1) Jxy ^ xy Similar formulas hold for higher correlation functions. This is identical to a connected expectation value w.r.t. the action S, denoted hereafter {■ ■ The first step is to show that correlations can all be ex- pressed as: i.'^xt'^yt) ~ I '^xt,x\t\T^xt,X2t2'^yt,xst2,'^yt,Xiti Jxi,ti<t 2:2*2, 2:3^3, 2:4*4) (4.2) Here TZ is the exact response function and T the exact effec- tive action (sum of IPI graphs) (with the choice e~'" for the probability and F^'') is symmetric). This is the standard rela- tion between connected correlation functions and the effective action (i.e. 1 -particle irreducible vertex functions IVF). The simplification here is that a priori the exact 2-point correlation function and vertices such as ri^L,-, could also contribute, but their contribution vanishes for T = at the depinning thresh- old. This is because {uxiUxm) time independent there, and then statistical tilt symmetry implies that all IVF's with at least one external u leg carrying frequency u> vanish when this fre- quency is set to zero (see section II A in [10]). The above formula (4.2) generaUzes straightforwardly to any cormected 2n-point correlation function of the field u in terms of rl^"]j. Next one can compute F^'^uu ir* perturbation, using the di- agrammatic rules for the non-analytic action arising from the expansion (2.6). Let us denote Ea the number of u external legs, riy the number of (unsplitted) vertices Aq, nj the num- ber of internal lines (response functions) and L the number of (momentum) loops. Then one has 2ni, — n/ — Eu and L = \ + ni — riy. Here one has Ey, = 4 and thus the lowest order contribution has n„ = 4, nj = 4, i.e. it is the 1-loop square (since two vertices n„ = 2 imphes n/ = and is disconnected, three vertices n„ = 3 imphes n/ = 2 and is 1 -particle reducible). The two loop corrections are diagrams with 5 vertices and so on. Similarly F?"! to lowest order is the 1-loop 2n polygon diagram. Since two u fields must come out (arrows) of each vertex, there are four possible diagrams corresponding to the 1-loop square, shown on figure 2. Each line entering a vertex cor- responds to one derivative of the A{u) function of the vertex. Thus from (2.6) we see that (a) is proportional to A'(0+)'', (b) and (c) to A"(0+)A'(0+)2 A(0) and (d) to A"{0+ fA{(]f. However, as we will show below using the so-called mount- ing property only (a) is non-zero. We found two ways to compute diagram (a) (as well as any other similar diagrams), a systematic but complicated way, and a simple way which uses a very important property and not yet fully elucidated of the field theory of depinning, the "quasistatic property" described below. To appreciate the ex- tent of the cancellations involved in this drastic simplification, we start by sketching the systematic method. To perform an actual calculation, since each A vertex in- volves fields at two times, one must switch to the splitted dia- grammatics, as described in [10]. Diagram (a) then becomes the sum of 16 splitted diagrams (two choices per vertex) repre- sented in Fig. 3 . Note that the last one is zero since it involves an acausal loop. That leaves 15 non-zero and non-trivial dia- grams. These diagrams correspond to the following. One first ex- pands <S^/4! using (2.6), which gives, schematically: A'(0+)'*. . — ^j^j — W1W2S12M12M3U4S34U34M5W6S56W56W7W8S78M78 • (4.3) In shorthand notations ui = iixi.ti, ^2 = u^-i.i-j' '^12 = ui — U2, S12 = sgn(fi — ^2) omitting all space and time integrals. On then carries the Wick contractions, yielding 4! 2U1U3U5U7S12S34S56S78 X (-R32 — R42){R54 — Re4){R76 — R86){Rl8 — R28) ■ (4.4) Rij = (uiUj) is the free response function. The factor of 2 comes from the two possible time orientations of the loop. Expanding the product of response functions yields the 16 di- agrams represented in Fig. 3, where space and time labels are ordered turning clockwise around the momentum loop. For il- lustration let us indicate the full expression of the first diagram in Fig. 3, in momentum space: ^ia.l)iPl2,ti;P23,t3;P3i,h;P41,t7) = (4.5) A'(0+)* / sgn(ti - t2)sgn(<3 - t4)sgn{t5 - h) J t2t4^tQtg XSgn(t7 — ts)Rpj^^t^-t2Rp2,t5-U^P3,t7-t6^PA,ti-t8 (a) (c) (b) (d) FIG. 2: The four 1-loop diagrams with unsplitted vertices which con- tribute to the 4-point irreducible vertex function F^*]^^ FIG. 3: The 16 1-loop diagrams with splitted vertices which cor- respond to diagram (a) in Fig. 2. The last one, which contains an acausal loop and thus vanishes, is added here for future convenience. with pij = p^ — pj the entering momenta and Rp^r = 9{T)e^P the free response function. Because of the sign functions the evaluation of these integrals, and of the other 14 non-vanishing diagrams is very tedious and was handled us- ing Mathematica. Adding all diagrams, massive cancellations occur The final result is very simple and given below. Let us now explain the simple method and the properties of the theory which lead to it. B. Theorems and other properties The simple way to compute the 4-point correlations (and higher ones) at depinning is based on the very important fol- lowing conjectured property: (4.6) Quasistatic property 1 : All correlation functions computed us- ing the diagrammatic rules of the quasistatic field theory of depinning at zero temperature and exactly at threshold, are independent of all time arguments. Using relations such as (4.2) for arbitrary times shows that an equivalent way to state this property is the following: Quasistatic property 2 : All T^?"'\{xiti, . . . ,X2nt2n) are independent of tl, ■ ■ ■ , t2n- This property, which appears as a physical requirement for the correct field theory of depinning, implies non-trivial prop- erties of the diagrammatics. Although we will not attempt to prove it here in full generality, we have checked it on many 7 examples, and believe that it works. We encourage the reader to contribute a valid proof. We will however state and prove some easier and useful properties below. Once the properties 1 and 2 are accepted, evaluation of the diagrams drastically simplifies, thanks to the following trick: Since the result does not depend on external times, one can take these times mutually infinitely separated, with some fixed (and arbitrary) ordering. Then one can integrate easily over all internal times since the order at each vertex is then specified and each sign function has a fixed value. One recalls that in the splitted diagrannmatics all non-vanishing T = diagrams are trees (see section 11 A in [10]). This can be seen on the fifteen non-vanishing diagrams of Fig. 3. Thus integrating in- dependently along each tree starting from the leaves yields one correlation function per link, since J^, Rq,t-t' = l/?^- Performing this calculation on all fifteen diagrams of Fig. 3 shows that they cancel pairwise since they differ only by a global sign, with the exception of the acausal diagram which is zero. Thus the final result is the same as if one had kept only (—1) times the acausal graph! Before giving the final result below, let us now state the easier properties. Theorem 1 (mounting trick): A diagram which contributes to an u-independent vertex function is if it contains a vertex, into which no response- function enters. Examples are diagrams (b), (c) and (d) on figure 2. This theorem thus ensures that only (a) is non- vanishing. Proof; The following figure demonstrates the principle. Note that it may be part of a larger diagram. Especially, there may be more response- functions entering into the upper disorder. The statement is that (4.7) Since no response-function enters into the lower disorder A{u — u'), due to the assumptions this gives A(0), with no dependence on time. Thus one can freely integrate over the response-function starting at time t. This integral for both di- agrams is / dtR{k,t) = The difference in sign comes from deriving the two different ends of the upper disorder ver- tex. Thus both contributions exactly cancel. Thus T^{^lf^,^-^{xiti, X2t2, xsts, X4ti) is given only by graph (a). We have not found a complete proof that it is independent of external times, but we can prove the weaker Lemma 1: ^uiuui^iti, X2t2, X2.h, Xiti) (4.8) (see graph (a) on figure 2) is independent of the most advanced time. Proof; First suppose that a response-function enters at the most- advanced time t' . Then there is the following cancellation = (4.9) The mechanism is the same as in the proof of theorem 1 ; since by assumption t' is the most advanced time, the argument of the right-most disorder can never change sign, and can be in- tegrated over, even though it is odd, i.e. A'(0+). Thus the remaining diagrams have the structure (4.10) This diagram is independent of t', as long as t' is the largest (external) time. Note that for a loop made out of two disorders, there is only one diagram remaining, namely It is manifestly time independent. Lemma 2; r^fifiC^^i^i v^2i2, x:,h,XiU) (4.12) (see graph (a) on figure 2) is independent of the differences in time, if those are very large. Proof: By inspection, one finds that by taking the external times, i.e. one time at each disorder, (infinitely) far apart, the remain- ing integrals become unambiguous. Thus integrating over the response-functions does not leave any time-dependence. 8 <II1> FIG. 4: The two contributions to {u^) ^ at leading order. As mentioned above, calculation of the diagram (a) be- comes then possible and one finds the Property (missing acausal loop): The diagram (a) is given by (—1) times the acausal loop, if there one replaces each response-function by a correlation- function. Intuitively this means, that if the acausal loop would give a contribution, then all diagrams would cancel. This seems to be a general mounting theorem. Check: One can calculate diagram A on Fig. 4 explicitly using none of the above theorems or conjectures. The result is a formidable expression, which has to be integrated over mo- menta. By evaluating it for given values of the momenta (not even necessarily conserving momentum), one can compare with the prediction of property 5. We found both expressions to be equal for any randomly chosen values of the momenta. The properties described here suggest the following Property (any loop): All graphs can be computed to any number of loops, using generalizations of the above rules. This is not attempted here but prehminary investigation suggests that the same mechanism holds for two loops with some simple end result related to the signs of possible "fermion loops". V. FINAL RESULT FOR THE FOURTH CUMULANT AND UNIVERSAL RATIO In this Section we compute the fourth cumulant D J XV as well as the ratio (kurtosis): R D 2 f Jxy (5.1) (5.2) which, according to the discussion in Sections II and III is expected to be universal and characterizes the deviations from the Gaussian approximation (for which R = Q). Below, we compute R at depinning both for short range (re- spectively long range) elasticity to lowest non-trivial order in e = 4 — d (respectively e = 2 — cf), i.e within a 1-loop cal- culation. However, since it turns out that the momentum inte- grals involved in the calculation depend very strongly on the dimension, we found it useful, and sometimes more accurate. to carry a 1-loop approximation directly in fixed dimension d. Also, since there is one exact result for a massive propagator, we also give the result in that case. We denote: 9{D = ^ (5.3) with the obvious change for long range elasticity g{q) = , and (see below) the massive propagator. The final result in the continuum is given by the sum of the two diagrams in Fig. 4: Jxy d'^q d'^k d'^p MQfgikfgipfgip + <i)9{p + k) +g{qfg{kfg{p)g(p + q)g{p + k)g{p + k + q) (5.4) d'^q_ (5.5) The combinatorics can be done as follows. There is a factor l/(4!2^). There are 4! ways to associate each one of the four external u to an unsplitted vertex. Say 1, 2 are now linked to and 3, 4 to M^. At each vertex a u field comes out. There are 2* ways to assign them to each splitted vertex. Then there is a unique set of four splitted points, one at each vertex, en- tering the acausal graph (which - in effect - is the only one arising, with the minus sign). But there are still 3 ways to join these four points in a loop: Two give the first integral, one the second, and finally for each case the orientation can be chosen in two ways. Let us go to the discrete model with periodic BC [34]. We recall that where here and in the following the term with n = is always excluded. In the limit of large L/a one finds: (5.7) -2A'(0+)^i-^''( E + \16 2tt' 1 1 (n2)2(TO2)2(p)(; + n)2(; + m)2(/ + m + n)^ 2n y- (5.8) 9 One can see that the ratio R will be universal since the 1- loop FRG fixed point equation taken at u = yields: (e - 2C)A(0) = (e/)A'(0+)2 (eJ) = LOlI ^ (47r)'^/2r(3- 87r2 (5.9) (5.10) The last identity is valid for d = 4. In fact, since the 1-loop FRG equation is universal, it holds as well for d = 4 as for d < 4. For d < 4, we use f d^p 1 1 U ^ 1 LdrJ = el = e y- (5.11) (5.12) This is all we need to compute the universal ratio. A. Calculation to lowest order in e = 4 - d The ratio i? = D/(2E,,G2^) is: E 1 1 (27r)^ J2neZ'' n» 1 1 (n2)2(TO2)2(Z2)(; + nf^i + rnf{i + ^ + nf _ (5.13) One finds, using ^e-*"' =6(3,0, e- (5.14) that in d = 4 1 1 ^ ^ = g y i'(e(3, 0, e-*)'^ - 1) = 10.2454 (5.15) Noting /(p) = Y.a^7d in2yiL^„yi One has: E E 1 E (^/W'-1850. (5.16) (n2)2(m2)2(Z2)(/ + nf{l + mf{l + m + nf w 980. (5.17) The final result is: R = -l.n^e^ « -0.13e2 (5.18) This results shows that R is quite small near d = A, but in- creases quite fast as the dimension is lowered. However the sums over the momenta depend very strongly on d (see below) and one should expect significant higher order corrections in e. Thus the result (5.18) is hkely to drastically overestimate the (absolute value of the) result in lower dimensions, which is why we now turn to an estimate in fixed dimension. B. 1-loop estimate in general dimension One can perform an estimate in general dimension, based on an arbitrary truncation on (i) 1-loop graphs (ii) neglect of the finite size scaUng function. In general dimension one has the dimensionless ratio: R = -{e-2Cf E 1 2 (n2)2(m2)2(i2)2(; _^ „)2(; + 1 (n2)2(TO2)2(;2)(; + nf{i + to)2(2 + ^ + (5.19) Let us give a table of values: ~ 45 2.16465 (d = 1) = 6.0268120 = 16.5 y- E 4725 2.00815 4.28143066080578 6.9458079272 (d = 2) (d = 3) (d=l) (d = 2) (d = 3) (5.20) (5.21) (n2)2(m2)2(Z2)2(; + nY^i + „)2 = 0.37342751117 (d = = 26.567 (d = = 240 (d = 1) 2) 3) (5.22) E n,m.,lGl (n2)2(m2)2(/2)(/ _l_ „)2(; _|_ „i)2(-/ -\-rn + ny = 0.069672062794960 = 14.138 = 143 In d = 1 one finds: R = -0.0867759691287 (^1-2^ « -0.00964177 (C = e/3) « -0.00347 (C = 1.2) (d=l) (d = 2) (d = 3) (5.23) (5.24) (5.25) (5.26) 10 In = 2 one finds: R = -0.4326 I 1 - 2^ -0.0481 (C = e/3) (5.27) In d = 3 one finds: R = -0.3297 1 -2 X -0.0366 (C = e/3) (5.28) where we have inserted various choices for C, including the 1-loop result C = e/3. One finds that already in d = 3 the 1-loop approximation is significantly lower than the extrapo- lation from (5.18) as discussed above. C. Long range elasticity in general dimension For LR elasticity the upper critical dimension is d^c = 2. The general expression for R is: R = -{e-2Cf E 2 {n^){m^){F)\l + n\\l + m\ 1 {n'^){m'^)\l\\l + n\\l + m\\l + m + n\_ (5.29) It is interesting to compute Rmd= 1. Using that: (5.30) (5.31) 3.847 (5.32) V ^ = — = 3.28987 ^ n2 3 Y — = ^ = 2.16465 E 1 E n,m.,l^Z one finds: {n^){m'^){P)\l + n\\l + m\ 1 {n'^){m'^)\l\\l + n\\l + m\\l + m + n\ 1.934 (5.33) R = -0.4109 1-2 w -0.04566 (C = e/3) D. Long range epsiion expansion (5.34) Similarly, one can perform an expansion in e = 2 — d. In d = 2 one has: e/= l/(27r) (5.35) 45 El tt' E : 6.0268120 {n'^){m'^){P)\l + n\\l + m\ E n,m,leZ^ (n2)(m2)|/||Z + n\\l + m\\l + m + n (5.36) 550 ± 20 (5.37) = 370 ± 10 (5.38) This yields i? = -e^ (1 - 2Ci)' (27r)2 ^ E 1 + {n^){m,^){l2)\l + n\\l + m\ 1 (n^)(m^ n||l + m||l + m + n| The result is R -6.17-e2 9 -0.686e^ (5.39) (5.40) again, a probable overestimation of the result if naively ex- trapolated tod= 1. E. Harmonic well, SR elasticity It is interesting to compare with the calculation in a massive scheme, i.e. an interface in a harmonic well. Setting g{q) := 1/(1-1- g'^), we have for d = 1 g{qf = - = 1.5708 g{qf = ^ = 1.62596 /// J q J h J p 16 9{<lf9{kf9{pfg{p + q)g{p + k) 163l7r3 = 1.62588 (5.41) (5.42) (5.43) /// J q J k J p _ 2457r3 ~ 5184 This gives the ratio 2366 31104 g{<lf'9{kf9{p)9{p + Q)g{p + k)g{p + k + q) 1.46538 (5.44) R = — (1 - 2C/e)2 = = -0.216369455 1215' " ' 10935 (5.45) It is interesting to compare the present result to one case (to our knowledge the only one apart from mean field models) where the full distribution ofu is known in a non- trivial disor- dered problem [30]. This is the static random field model in (i = in a harmonic well (i.e. the massive case), the so-called toy model. The exact result there is: = -0.080865 . . . (5.46) 11 The present 1-loop approximation for the problem of de- pinning, continued to d = would give the larger result R = —1/3. It is unclear at present whether the difference be- tween the two results indicates that the 1-loop approximation is unsatisfactory so far from d = 4, or if statics and depinning have radically different values of R. VI. BEHAVIOR AT THE CRITICAL DIMENSION In this Section we reexamine isotropic depinning, and stat- ics, exactly in d = 4. We solve the RG equations in d = 4 and obtain the behavior of the correlation function. Contrarily to periodic systems at the upper critical dimension [31], non- periodic objects such as interfaces submitted to either random bond or random field type disorder exhibit a roughness, which is a power of a logarithm. The FRG flow equation (/3-function) for the (renormalized) force correlator has a good limit for d = 4. If one defines Aiiu) = 87r2;2Ci-iA,(u/-Ci) (6.1) with I = ln(A/m) (A some UV cutoff) then the function A;(u) satisfies diA(u) = (1 - 2Ci)A(m) + CiuA'iu) 1 ~2 + 1-^I32{A) . (6.2) (A(«)-A(0))^ where, for depinning, /32(A) = [(A(,i) - A(0))A'(w)2]" + A'(0+)2A"(,i) . and Ci = C/c — 1/3 is the 1-loop value, see e.g. [10]. It is then easy to see that the function A; {u) converges towards the 1-loop fixed point with the foUowing asymptotic corrections: Ai{u) = A*{u)+J2l~^"Uu) + j(}[[AT'p2iA*) n (6.3) where (J3[)^^ is the inverse of the linearized 1-loop (3- function and the a>„ are the 1-loop eigenvalues. Using (6.1) with = 1/3 yields the result for the correla- tion function at g' = 0: {ugU_g)\g^rn = m-^A(O) (6.4) = cm"^ ln(A/m)-3 (1 + 0(1 / ln(A/m))) as m ^ 0, with c = 87r^A*(0), both for statics and depin- ning; the difference hes in the subdominant piece. Within the present approach using the renormaUzation scheme atq = 0, the 2-point correlation function at non-zero q can be computed from the renormalized (uniform) effective action by resum- ming an infinite set of diagrams. Using the standard finite size scaling ansatz allows to obtain the other Umit of the scaling function, where A » g » m. To lowest order (one loop) in the renormalized disorder one has [11] (g2 + m^fiugu.g) = (A(0) - A'{0+f{I{q) - 7(0)) + . . .) ^^^^ (p'+m2)((p + (7)2+m2) • ^^-^^ This gives: {ugU-g} = ciq"^ + m^y^ ln( — ] |1 3ln(^) + ... Assuming scaling, i.e. that the function (1 — + . . .) {1 + x)~i one finds that: (6.6) and thus: {Ugu.g)^q-' [in(^) (u,-uo)2~(lna;)2/3 1/3 (6.7) (6.8) VII. STABILITY OF THE 1 -LOOP FIXED POINT Here we analyze the stabiUty of a functional fixed point. Two cases have to be distinguished: (a) There is the freedom to rescale the field u while at the same time rescaling the disorder correlator. This includes the random bond and random field interface models. (b) There is no such freedom, since the period is fixed by the microscopic disorder This is the case for a charge density wave (random periodic problem), but also for the random field bulk problem, in its treatment via a non-Unear sigma model. We first analyze the simpler case (b). A. Stability of a functional fixed point; periodic case Be the flow-equation given by diR{u) = p[R]{u) = eR{u) + f[R, R]{u) , (7.1) where / is some biUnear form of R, which contains at least one derivative for each R. A similar equation of course exists for A(m) = —R"{u) and the corresponding /3[A](u). Suppose R* (u) is the non-trivial fixed point of order e, i.e. P[R*] = 0. Two eigenfunctions and eigenvalues above R* {u) can be identified. (i) The constant mode dR{u) = 1 with eigen-value e. (As long as it is permissible physically.) (ii) The first subleading eigenfunction 5R{u) = R*{u) with eigenvalue — e. Proof: For case (i), we have for k <^ 1 di{R*{u) + k) = ^[R* + k]{u) = eK , (7.2) since / does not couple to the constant by assumption. This proves (i). 12 For case (ii), we have {k <C 1): di {R*{u) + kR*{u)) = (3[R*{1 + k)]{u) = eR*{u){l + k) + (1 + Kff[R*, R*]{u) (7.3) Subtracting /3[i?*](u) = on the r.h.s. and expanding for small K yields di {R*{u) + kR*{u)) = eKR*{u) + 2Kf[R*,R*]{u) = -eKR*{u) , (7.4) where in the last equation we have again used the fixed-point condition l3[R*]{u) = 0, i.e. eR*{u) + f[R*,R*]{u) = 0, to eliminate f[R*, R*]{u). This proves (ii). Using the same line of arguments, it is easy to see that when starting from R{u) ^ R* (u) with some arbitrary ampUtude, the flow is always remaining on the critical manifold spanned by R*{u). Of course, there are in general more eigen-functions and eigen-values. See [10] for an exphcit example. B. Perturbations of the fixed point in presence of the freedom to rescale We state the following Theorem: The differential equation of the form -mdmAiu) = /3[A] /3[A] = (e - 20 A(«) + CuA'iu) + /[A, A] (7.5) where the symmetric functional /[A, A] transforms under A(u) — > K^^A(K?i) in the same way as A, has the two eigen- functions and eigenvalues of perturbations around the fixed point /3[A*] = zoiu) = wA'(w) - 2A(u) (7.6) Ao = . (7.7) zi{u) = CuA'{u) + {e- 2C)A(m) (7.8) Ai = -e . (7.9) (We note A instead of A* for the fixed point for simplicity of notations.) Note that the assumptions are satisfied by the 1-loop flow-equation (RF-case). Proof: Be (3[A]{u) = 0. Due to the assumptions, for all k P[k-^A]{ku) = . (7.10) Deriving w.r.t. k gives with (7.5) at k = 1 (e - 20 (uA'(u) - 2A(u)) + (u (uA'(u) - 2A{u))' +2f {A{u), uA'{u) - 2A{u)) = (7.11) This is nothing but the eigenvalue equation for the perturba- tion zo{u) about the fixed point /3[A] = and proves the existence of the solution zo{u) with eigenvalue Aq = 0. Note that this redundant operator with eigenvalue persists to all orders in perturbation theory. Let us turn to the next solution. Multiplying (7.11) with C, and adding 2e times /3[A] =0 gives (e - 20 (CuA'(?i) + (2e - 20A(w)) + Cu (CuA'(w) + (2e - 20 A(u))' + 2/ {A{u)XuA'{u) + (e - 20A{u)) = , (7.12) where we used the bilinearity of A{u). Rearranging yields (e-20 {CuA'{u) + (e-20A{u)) + Cu{C,uA'{u) + {€-20A{u))' + 2/ (A(m), CuA'iu) + (e - 2QA{u)) = -e(CuA'(u) + (e - 20A(w)) . (7.13) This equation is nothing but the eigenvalue-equation for the perturbation zi (m) = (uA'{u) + (e—2C,)A{u) about the fixed point P[A] = with eigen-value Ai = — e. q.e.d. Remark: Consider the case of short-range disorder, i.e. that A{u) falls off rapidly, and is monotonic. Usually, the leading fixed point solution has no knot (no u such that A{u) = 0). Then zq{u) has no knot and zi (u) has one knot. Eigenvalues should be ordered (like in quantum mechanics) due to their number of knots. Thus we should have found the two leading solutions for short-range disorder. This is confirmed by the numerical analysis given in subsection Vll C. Corollary: (Random bond case) The differential equation of the form -mdmR{u) = [3[R]{u) (7.14) (3[R]{u) = (e - 40i?(u) + C,uR'{u) + f[R,R]{u) where the synnmetric functional f[R,R] transforms under R{u) k~'^R{ku) in the same way as R, has the two eigen- functions and eigenvalues of perturbations around the fixed point f3[R*] = zoiu) = uR'{u) ~ AR{u) (7.15) Ao = . (7.16) zi{u) = CuR'{u) + {€- 40i?(u) (7.17) Ai = -e . (7.18) Note that the assumptions are satisfied by the 1-loop flow- equation (random bond case). Proof: This can either be proven along the same lines as for (7.5) or by deriving zq (u) and zi (u) twice w.r.t. u and then using the theorem for the random field case (7.5). C. Numerical analysis of the RF-fixed point We start from the 1-loop flow equation -mdmA{u) = (e - 20A(u) -|- CuA'{u) -l[{A{u)-AiO))^]" . (7.19) 13 1, 0.8- 0. 6- 0.4: 0.2 -0.2 FIG. 5: The solution zi(u) (in black, with 2:(0) — 1) and asymptotic behavior as given in (7.24) (blue/bright). It has the following solution [5, 10] yi{u) - Inyi(u) = 1 + (7.20) Perturbations around this solution satisfy the differential equa- tion -mdm [A(u) (l-3A)z(u) = ~ \ez{u) uz'{u) - [{yi{u) (7.21) 2/i(0))(z(u)-z(0))]" (7.22) In order to have a criterion for the numerical integration, one has to determine the behavior at infinity. Using that for u oo yiiu) (7.23) and assuming exponentially fast decay for z{u), one finds that the asymptotic behavior is z{u) w z(0) (^2 - i)e-i-"V2 2 + 3A (7.24) Quite surprisingly, the asymptotic behavior is fixed, including its amplitude[38]. In any case, slower power law decay is ruled out on physical ground since we are considering short range disorder. The solutions for A = and A = — e are given in Eqs. (7.6) ff. The solution zi [u) is dominant, and gives the correction to scaling exponent w = — e. Note that this exponent is mi- nus the engineering dimension of the coupling, as is the case in standard field-theory[32], and also for the random periodic class [10]. Ref. [26] cites the value oj = — e/3. We find the corresponding numerical solution to decay as u^^, incompat- ible with (7.24) and physically unacceptable [34]. The question arises, of whether there are more solutions with fast decay. Intuitively, one would expect this: Making A more negative, the solution overshoots and one might think of fine-tuning A such that it approaches the axis for large u from above. However this is incompatible with the asymptotic behavior in Eq. (7.24), which predicts, that all solutions for A < —2/3 converge from below. In fact, we have not been able to find any further solution, and we conjecture that there is none. It would be interesting to prove this rigorously. This behavior is in contrast to the random periodic case, solved in [10], which has infinite many subleading contributions. VIII. CONCLUSION In this paper we have explored further properties of the field theory of depinning. We have defined and computed universal observables, such as the distribution of the interface width and the ratio R of the connected 4-point cumulant to the square of the 2-point one (kurtosis). This ratio measures the deviations from a Gaussian approximation which we have also used to obtain the universal distribution. Higher order connected cu- mulants can be computed in a similar way to one loop using polygon diagrams, and one should be able to reconstruct the full distribution from them. Other properties of the theory such as the behavior at the upper critical dimension and the finite size scaling behavior have been clarified. All calcula- tions in the present paper are of interest for comparison with numerical simulations, existing ones [16] or in the near future. In the process of computing the 4-point function we dis- covered massive cancellations between diagrams. We traced this back to the physically expected property that correlations exactly at depinning should be time independent. A diagram- matic proof of this property is still incomplete, but we have provided some convincing elements in that directions. As a result the correlation functions can be computed in a much simpler way. Thus there seems to be an underlying theory, with "quasi-static" diagrams (i.e. not containing time explic- itly), with some additional rules. We have understood these rules to lowest ( 1 -loop) order and it would be of high interest to understand -and prove- them to all orders. It is even pos- sible that there exists a simpler formulation of the theory at depinning in terms of, e.g. effective fermions. The fermionic character is suggested by the cancellation of all diagrams ex- cept for the "acausal loops" with a minus sign. We thus en- courage further examination of this fascinating question and full elucidation of the field theory which describes depinning. APPENDIX A: DIFFERENT BOUNDARY CONDITIONS a. Periodic boundary conditions A periodic function u with period L satisfies u{x + Lei) — u{x) , (A.l) where the ii are orthonormal. It can be written as Uk e ' (A.l) where summation runs over all k, such that ki = niiii/L, rii g Z. 14 b. Open boundary conditions To simplify the notation, we give all formulas for one dimen- sion; generalizations are straightforward. Suppose the function f{x) is defined on [0, L]. Then a func- tion g{x) can be defined by the following prescription / \ _ / /(a;) for 0<x<L 9{x) - I ^^2L -x) for L<x<2L g{x) can be prolonged to a periodic function with period 2L, i.e. g{x+2L) = g{x), and satisfies by construction in addition g{2L -x)= g{x) . (A.4) In the basis needed to construct functions with period 2L, we have to restrain ourselves to since the sin do not satisfy (A.4). The such constructed set of functions f{x), x G [0, L] has Neumann-boundary condi- tions at X = and x = L. From (A.5), we infer that the number of modes is reduced by a factor of 2 (compared to the case of closed boundary conditions), but the construction does not change any observable constructed from w"^ or any energy, all based on the symmetry relation (A.4). Importantly, the modes have all mean 0, which is not the case for other basis, e.g. when using anti-periodic functions. Also note, that this ansatz reproduces the formula in [33]. As an interesting consequence, we observe that the following systems lead to the same distribution • An elastic object with TV degrees of freedom, and closed boundary conditions. • An elastic object with 2N degrees of freedom, and open boundary conditions. The simplest example is a 1 -dimensional random walk with closed boundary conditions, and a 2-dimensional random walk with open boundary conditions, which can be checked numerically [39]. Acknowledgments We are grateful to W. Krauth and A. Rosso for an ongoing collaboration and numerous stimulating discussions, and we thank E. Brezin, O. Narayan and J. M. Schwarz for very useful remarks. [1] A.P.Young, Spin glasses and random fields. World Scientific, Singapore, 1997. [2] M. Kardar, Nonequilibrium dynamics of interfaces and lines, Phys. Rep. 301 (1998) 85-112. [3] D.S. Fisher, Collective transport in random media: from super- conductors to earthquakes, Phys. Rep. 301 (1998) 113-150. [4] G. Blatter, M.V. Feigel'man, V.B. Geshkenbein, A.I. Larkin and V.M. Vinokur, Vortices in high-temperature superconductors. Rev. Mod. Phys. 66 (1994) 1125. [5] D.S. Fisher, Interface fluctuations in disordered systems: 5 — e expansion, Phys. Rev. Lett. 56 (1986) 1964-97. [6] T. Nattermann, S. Stepanow, L.-H. Tang and H. Leschhom, Dy- namics of interface depinning in a disordered medium, J. Phys. II (France) 2 (1992) 1483-8. [7] O. Narayan and D.S. Fisher, Critical behavior of sliding charge-density waves in 4- epsilon dimensions, Phys. Rev. B 46(1992) 11520-49. [8] O. Narayan and D.S. Fisher, Threshold critical dynamics of driven interfaces in random media, Phys. Rev. B 48 (1993) 7030-42. [9] P. Chauve, P. Le Doussal and K.J. Wiese, Renormalization of pinned elastic systems: How does it work beyond one loop?, Phys. Rev. Lett. 86 (2001) 1785-1788, cond-mat/0006056. [10] P. Le Doussal, K.J. Wiese and P. Chauve, 2-loop functional renormalization group analysis of the depinning transition, Phys. Rev. B 66 (2002) 174201, cond-mat/0205108. [1 1] R Le Doussal, K. Wiese and R Chauve, Two loop FRG study of pinned manifolds, in preparation. [12] A. Rosso and W. Krauth, Monte carlo dynamics of driven flux lines in disordered media, cond-mat/0107527 (2001). [13] A. Rosso and W. Krauth, Origin of the roughness exponent in elastic strings at the depinning threshold, Phys. Rev. Lett. 87 (2001) 187002. [14] A. Rosso and W. Krauth, Roughness at the depinning thresh- old for a long-range elastic string, Phys. Rev. E 65 (2002) 025101/1^. [15] A. Rosso, Depiegeage des varietes elastiques en milieu aleatoire, These de Doctoral de I'Universite Paris VI, Paris 30 Sept. 2002. [16] A. Rosso, W. Krauth, P. Le Doussal, J. Vannimenus and K.J. Wiese, Universal interface width distributions at the depinning threshold, cond-mat/0301464 (2003). [17] D.S. Fisher, Sliding charge-density waves as a dynamical criti- cal phenomena, Phys. Rev. B 31 (1985) 1396-1427. [18] DS. Fisher, Random fields, random anisotropics, nonlinear sigma models and dimensional reduction, Phys. Rev. B 31 (1985) 7233-51. [19] J. Vannimenus and B. Derrida, A solvable model of interface depinning in random media, J. Stat. Phys. 105 (2001) 1-23. [20] P. Le Doussal and K.J. Wiese, Functional renormalization group at large N for random manifolds, Phys. Rev. Lett. 89 (2002) , cond-mat/0109204vl. [21] L.-H. Tang, M. Kardar and D. Dhar, Driven depinning in anisotropic media, Phys. Rev. Lett. 74 (1995) 920-3. [22] R. Albert, A.-L. Barabasi, N. Carle and A. Dougherty, Driven interfaces in disordered media: determination of universality classes from experimental data, Phys. Rev. Lett. 81 (1998) 2926-9. [23] P. Le Doussal and K.J. Wiese, Functional renormalization group for anisotropic depinning and relation to branching pro- cesses, Phys. Rev E ??? (2003) ???, cond-mat/0208204. [24] A. Rosso, A.K. Hartmann and W. Krauth, Depinning of elastic manifolds, cond-mat/0207288 (2002). [25] L. Roters, S. Lubeck and K. D. Usadel, Functional renormal- 15 ization group for anisotropic depinning and relation to branch- ing processes, cond-mat/0207137 (2002). [26] S. Ramanathan and D.S. Fisher, Onset of propagation of planar cracks in heterogeneous media, Phys. Rev. B 58 (1998) 6026- 46. [27] J. M. Schwarz and Daniel S. Fisher, Depinning with dynamic stress overshoots: A hybrid of critical and pseudohysteretic be- havior, cond-mat/0204623 (2002). [28] M. Plischke, Z. Racz and R.K.R Zia, Width distribution of curvature-driven interfaces: a study of universality, Phys. Rev. E 50 (1994) 3589-93. [29] G. Foltin, K. Oerding, Z. Racz, R.I. Workman and R.K.P. Zia, Width distribution for random-walk interfaces, Phys. Rev. E 50 (1994) R639^2. [30] P. Le Doussal and Cecile Monthus, Exact solutions for the statistics of extrema of some random ID landscapes. Appli- cation to the equilibrium and the dynamics of the toy model, cond-mat/0204168. [31] R. Chitra, T. Giamarchi and Le Doussal, Disordered periodic systems at the upper critical dimension, Phys. Rev. B 59 (1999) 4058-65. [32] J. Zinn- Justin, Quantum Field Theory and Critical Phenomena, Oxford University Press, Oxford, 1989. [33] T. Antal, M. Droz, G. Gyorgyi and Z. Racz, Roughness distri- butions for 1 / f" signals, Phys. Rev. E 65 (2002) 046140/1-12. [34] To be precise, the calculations of this paper assume that the zero mode u is exactly set to zero. In that case momentum sums in internal lines (loops) can exclude g = in any diagram. (They are excluded by construction in external lines.) This restriction does not affect the discussion of the OA but is important for the calculations of R in section IV. The present calculation would apply directly to a numerical simulation where J2x '"^ = is enforced in each disorder configuration (pinned zero mode). We will discuss the precise connection between the boundary con- ditions chosen in simulations such as [16] and field theoretical calculations in a forthcomming publication. [35] This has been tested for the data of [27]. Even though there are no conclusive results, ui = — e is favored over ui = — e/3. We thank J. M. Schwarz for this private communication. [36] A subtle point in that construction is that if one defines Ao per- turbatively from A to a given order then Ao is not the original bare action (which is analytic) thus there is no contradiction in Ao being non-analytic. In a sense introducing Ao is just a trick to express a closed equation for the flow of A to the same order. [37] To extend the 2-loop ERG (2.5) to that periodic b.c. scheme, one would have to reevaluate 2-loop integrals with discrete sums, or carry other changes in the scheme. This however goes beyond the present 1-loop study. [38] We thank O. Narayan for this observation. [39] W. Krauth and A. Rosso, private communication.