Instruction: Write an abstract for the following acadamic paper.

Input:
Paper: we have recorded the flight of a fly during take off and landing using digital high speed photography . it is shown that the dynamics of flexible wings are different for these procedures . during this observation fly flew freely in a big box and it was not tethered .
Abstract: in this fluid dynamics video , we demonstrated take off and landing of a fly . the deformation of wings is in focus in this video .

Paper: quenched disorder induces an extrinsic inhomogeneity in otherwise translationally invariant pure systems . as a result , the equilibrium and out - of - equilibrium physics of disordered systems may be influenced by rare collective events , such as `` avalanches '' , statistically unlikely regions , such as `` droplets '' or `` griffiths regions '' , and by the proliferation of `` metastable states '' . our objective is to develop a theory describing the long - distance physics of such systems through a nonperturbative functional renormalization group ( np - frg ) method and our first focus is the equilibrium behavior of the random - field model . we showed in previous papers , which we refer to as i@xcite and ii@xcite , that the effect of avalanches and droplets can be captured by an approach based on the rg flow of the cumulants of the renormalized disorder , provided that the full functional dependence of the latter is accounted for.@xcite a functional rg is therefore required to let the singular behavior due to rare events or regions emerge along the flow , as also shown in the case of the equilibrium and forced behavior of manifolds in a random environment.@xcite    the issue of metastable states , which , at zero temperature where the concept is well defined , refers to the presence of many minima of the microscopic hamiltonian ( bare action ) not simply related by symmetry transformations , is a recurring conundrum in theories of disordered systems . in the case of the random - field ising model ( rfim ) under study , it has a clear manifestation . the critical behavior of the model being controlled by a zero - temperature fixed point,@xcite the long - distance physics can be described through the properties of the ground state which , due to the presence of the random field , is obtained as the solution of a stochastic field equation . by standard field - theoretic manipulations , this leads to a theory expressed in terms of superfields . parisi and sourlas@xcite showed that a supersymmetry of the theory , more specifically the invariance under rotations of the underlying superspace , implies the property of dimensional reduction , according to which the critical behavior of the rfim in @xmath4 dimensions is identical to that of the pure ising model in @xmath5 dimensions . the property however has been proven to be wrong in low enough dimension.@xcite at the same time , it has been understood that the superfield construction breaks down because of the presence of many solutions of the stochastic field equation.@xcite however , in the parisi - sourlas formalism , it is not possible to disentangle breaking of superrotational invariance and collapse of the formalism due to the appearance of multiple solutions .    in the companion paper , henceforth referred to as paper iii,@xcite we have shown how to resolve the above conundrum and to combine an extended superfield approach with the np - frg formalism . in particular , ground - state selection can be achieved by adding a weighting factor involving an auxiliary temperature and letting at the end of the manipulations the auxiliary temperature go to zero in the exact np - frg equations for the cumulants of the renormalized disorder . the resulting property of `` grassmannian ultralocality '' then allows one to specifically investigate supersymmetry ( superrotational invariance ) and its spontaneous breaking along the rg flow . such an investigation is the purpose of the present paper .    the outline of the article is as follows . in sec . [ sec : model ] , we briefly summarize the main steps of the np - frg in a superfield formalism that we have developed in paper iii . this allows us to recall definitions , notations , and rg equations that will be used in the present article . in particular , we stress two important , and distinct , formal properties of the superfield theory :  grassmannian ultralocality " and  superrotational invariance " . we also consider the rg flow of the ward - takahashi identities associated with the latter and the consequences for the choice of the infrared regulator . in sec . [ sec : surot invariance ] , we show through our np - frg formalism that the superrotational invariance nonperturbatively leads to dimensional reduction . we conclude the section by building a scenario for a spontaneous breaking of the superrotational invariance along the rg flow that is based on our previous results on the appearance of a  cusp " in the functional dependence of the cumulants of the renormalized random field,@xcite and we propose a continuation of the np - frg flow equations when spontaneous breaking has taken place . [ sec : approximation scheme ] is devoted to the development of a supersymmetry - compatible nonperturbative approximation scheme for the exact np - frg equations . it relies on combined truncations of the cumulant expansion and the expansion in spatial derivatives of the field . we also provide details on the numerical resolution . in sec . [ sec : results ] we present the results . we show in particular that breakdown of the dimensional reduction predictions for the critical exponents of the rfim occurs below a critical dimension @xmath0 . we compute the critical exponents as a function of dimension down to @xmath2 and we find good agreement with the best available estimates in @xmath2 and @xmath3 . finally , we show that scaling is described by three independent exponents , contrary to a proposed conjecture.@xcite    a short account of this work has appeared in ref . we start by briefly recalling the main features of the formalism presented in the preceding article , with the associated definitions and notations .    in paper iii,@xcite we have developed a np - frg theory for describing the equilibrium long - distance physics of the rfim that is based on an extension of the parisi - sourlas@xcite supersymmetric formalism . the latter relies on the fact that the critical behavior of the model is dominated by disorder - induced fluctuations , thermal fluctuations being subdominant , and can therefore be studied by looking directly at zero temperature . the equilibrium properties are then described by the ground state which is solution of the following stochastic field equation @xmath6}{\delta \varphi(x ) } = j(x),\ ] ] where we have added an external source ( a magnetic field ) @xmath7 conjugate to the @xmath8 field and the action @xmath9 $ ] is given by @xmath10 where @xmath11 , @xmath12 , with @xmath13 a random ( magnetic ) field sampled from a gaussian distribution of zero mean and variance @xmath14 . two key ingredients of our extension of the parisi - sourlas formalism are :    \(1 ) the need to consider multiple copies ( or replicas ) of the original system with the same disorder , each copy being coupled to a different applied source , in order to generate cumulants of the renormalized disorder with their full functional dependence , thereby allowing for the emergence of a nonanalytic behavior in the field arguments ,    \(2 ) the introduction of a weighting factor involving an auxiliary temperature @xmath15 to the solutions of the stochastic field equation , so that when @xmath15 appraches 0 only the ground state contributes to the generating functional .    by using standard field - theoretical techniques,@xcite one then ends up with a superfield theory for multiple copies in a curved superspace , with @xmath16=\exp ( \mathcal w^{(\beta)}[\{\mathcal j_a\}])=\overline{\exp(\sum_{a=1}^n\mathcal w_h^{(\beta ) } [ \mathcal j_a ] ) } \\&= \int \prod_{a=1}^{n}\mathcal q\phi_a   \exp \bigg(-s^{(\beta)}[\{\phi_a\ } ] +   \sum_{a=1}^{n } \int_{\underline x }   \mathcal j_a(\underline x )   \phi_a(\underline x)\bigg ) , \end{aligned}\ ] ] where the multicopy action is given by @xmath17 = \sum_{a=1}^{n } \int_{\underline{x } } \left [ \frac 12 ( \partial_{\mu}\phi_a(\underline{x}))^2+u_{b}(\phi_a(\underline{x}))\right ] \\&- \frac{\delta_b}{2}\sum_{a_1=1}^{n}\sum_{a_2=1}^{n } \int_{x}\int_{\underline{\theta}_1\underline{\theta}_2 } \phi_{a_1}(x,\underline{\theta}_1)\phi_{a_2}(x,\underline{\theta}_2 ) . \end{aligned}\ ] ] in the above equations , we have introduced a superspace ( coordinates @xmath18 ) comprising the @xmath4-dimensional euclidean space ( coordinates @xmath19 ) and a @xmath20-dimensional grassmannian space ( anticommuting coordinates @xmath21 ) ; the metric is flat in the euclidean sector and curved in the grassmannian one with the curvature proportional to @xmath22 . ( for instance , the integral over superspace is defined as @xmath23 . ) we have also defined superfields @xmath24 , with one auxiliary bosonic (  response " ) field @xmath25 and two auxiliary fermionic (  ghost " ) fields @xmath26 and @xmath27 , and associated supersources @xmath28 . the action in eq . ( [ eq_superaction_multicopy ] ) is invariant under a large group of ( bosonic and fermionic ) symmetries . when expanded in increasing number of sums over copies , the generating functional of the connected green s functions , @xmath29 $ ] , gives access to the cumulants of the random generating functional @xmath30 $ ] . we have next applied the np - frg formalism to this superfield theory . this proceeds by first adding an infrared ( ir ) regulator that enforces a progressive account of the fluctuations to the bare action , @xmath31 where the two cutoff functions @xmath32 and @xmath33 are related through @xmath34 with @xmath35 the euclidean momentum , @xmath36 the strength of the renormalized random field , and @xmath37 the field renormalization constant . this ensures that all symmetries of the theory are satisfied . this includes the superrotational invariance found when the theory is restricted to a single copy and to @xmath38 : it then corresponds to @xmath39 , a property that is valid at the microscopic ( uv ) scale @xmath40 ( see also below ) . one next introduces the effective average action@xcite @xmath41=-\mathcal w_k^{(\beta)}[\{\mathcal j_a\}]+\sum_a\int_{\underline x } \phi_a(\underline x ) \mathcal j_a(\underline x ) -\delta s_k^{(\beta)}[\{\phi_a\}]$ ] , whose dependence on the ir cutoff @xmath42 is governed by an exact renormalization - group equation ( erge ) . the effective average action can also be expanded in increasing number of sums over copies , each term of the expansion being then related through the legendre transform to the cumulants of @xmath30 $ ] . we refer to the @xmath43th order term of the expansion of @xmath41 $ ] as the @xmath43th `` cumulant of the renormalized disorder '' . these expansions in increasing number of sums over copies lead to systematic algebraic manipulations that allow one to derive from the erge for @xmath41 $ ] a hierarchy of coupled erge s for the cumulants of the renormalized disorder . we have unveiled an important property of the random generating functional @xmath30 $ ] : the latter satisfies a specific dependence on the grassmann coordinates , which we have called  grassmannian ultralocality " , if and only if a unique solution of the stochastic field equation is included in its computation . this translates into an  ultralocal " property of the cumulants of the renormalized disorder .  grassmannian ultralocality `` becomes a property of the superfield theory when @xmath44 ; it is also asymptotically found for finite @xmath22 when @xmath45 ( after going to dimensionless quantities ) . it then reflects the desired ground - state dominance . when this property is satisfied ( _ e.g. _ by setting @xmath46 ) , the erge s for the cumulants simplify and only involve  ultralocal '' quantities that can be evaluated for physical fields @xmath47 , _ i.e. _ for superfields that are uniform in the grassmann subspace .    for illustration , we give below the erge s for the first two cumulants under the property of  grassmannian ultralocality " , which will be needed in the following : @xmath48= \\&- \dfrac{1}{2 } \tilde{\partial}_t \int_{x_2x_3}\widehat{p}_{k;x_2x_3}[\phi_1 ] \big(\gamma_{k2;x_2,x_3}^{(11)}\left[\phi_1,\phi_1\right ] - \widetilde{r}_{k;x_2x_3}\big ) \end{split}\ ] ] and @xmath49=\\ & \dfrac{1}{2 } \tilde{\partial}_t \int_{x_3x_4}\big \{- \widehat{p}_{k;x_3x_4}\left[\phi_1\right ] \gamma_{k3;x_3,.,x_4}^{(101)}\left[\phi_1,\phi_2,\phi_1\right ] + \\ & \widetilde{p}_{k;x_3x_4}\left[\phi_1,\phi_1\right ] \gamma_{k2;x_3x_4,.}^{(20)}\left[\phi_1,\phi_2\right ] + \frac{1}{2}\widetilde{p}_{k;x_3x_4}\left[\phi_1,\phi_2\right ] \\ & \times \left ( \gamma_{k2;x_3,x_4}^{(11)}\left[\phi_1,\phi_2\right ] - \widetilde{r}_{k;x_3x_4}\right ) + perm(12)\big \ } , \end{split}\ ] ] where @xmath50 denotes the expression obtained by permuting @xmath51 and @xmath52 , @xmath53 , @xmath54 is a short - hand notation to indicate a derivative acting only on the cutoff functions ( _ i.e. _ , @xmath55 ) , and the propagators @xmath56 and @xmath57 ( denoted by @xmath58}$ ] and @xmath59}$ ] in paper iii@xcite ) are defined as @xmath60=\left(\gamma _ { k,1}^{(2 ) } [ \phi ] + \widehat r_k\right ) ^{-1}\ ] ] and @xmath61= \widehat { p}_{k } [ \phi_1 ] ( \gamma _ { k,2}^{(11)}[\phi_1 , \phi_2 ] -\widetilde r_k ) \widehat { p}_{k } [ \phi_2 ] , \ ] ] and superscripts indicate functional differentiation with respect to the field arguments . generically , the flow of @xmath62 $ ] involves three types of quantities : the propagators @xmath56 and @xmath57 , second functional derivatives of @xmath63 in which all the arguments are different , and second functional derivatives of @xmath64 with two of their arguments equal to each other . a graphical representation of the hierarchy of erge s is provided is appendix c of the companion paper iii .    in the present paper we will focus on these flow equations , which correspond to the @xmath44 limit and describe the equilibrium properties associated with the ground state ( recall that the `` bath '' temperature is also at @xmath65 ) . we will briefly comment on the rg flow of corrections to  grassmannian ultralocality " when @xmath22 is large but finite . finally , it is worth stressing that to obtain the flow equation for @xmath66 $ ] with its full functional dependence on the @xmath43 field arguments , one needs to consider at least @xmath43 copies . if one considers less than @xmath43 copies , the equation necessarily involves @xmath63 and its derivatives in which several of the arguments are equal . as we will show later , breaking of supersymmetry and breakdown of dimensional reduction are precisely related to the analyticity properties of the functionals when taking the limit of equal arguments . formally , the whole hierarchy of flow equations for the cumulants can thus be obtained by considering an arbitrary large number of copies . when  grassmannian ultralocality " is satisfied , the curvature @xmath22 disappears from the flow equations [ see eqs . ( [ eq_flow_gamma1_ulapp],[eq_flow_gamma2_ulapp ] ) ] . formally , the latter are then the same as those obtained from considering a flat grassmann subspace . we have shown in paper iii that in this case , and provided one restricts the supersources such that one effectively recovers a one - copy theory , the theory is invariant under  superrotations " that mix the euclidean and grassmannian sectors of the superspace . the associated generators are @xmath67 and @xmath68 . as any linearly realized continuous symmetry , superrotational invariance leads to a set of ward - takahashi ( wt ) identities for the one - particle irreducible ( 1pi ) generating functional , _ i.e. _ at scale @xmath42 the effective average action . for a flat superspace and a restriction to one copy , the wt identity for the effective average action reads @xmath69=0\ ] ] and similarly with @xmath70 . one can now check that the above wt identities are _ a priori _ stable under the rg flow , _ i.e. _ @xmath71=0,\ ] ] where @xmath72 generically indicates a component @xmath73 and , since there is no curvature , @xmath74 . to prove eq .  ( [ eq_flow_ward ] ) , we rewrite the erge for the effective average action when the superfield theory is restricted to one copy and is considered in a flat superspace ( @xmath38 ) : @xmath75=\frac 12 \int_{\underline{x}_1\,\underline{x}_2 }    \left ( \partial_t \mathcal r_{k;\underline x_1,\underline x_2}\right )   \mathcal p_{k;\underline x_1,\underline x_2}[\left\lbrace \phi_a \right\rbrace ] , \ ] ] where the regulator function @xmath76 with @xmath77 and the modified ( full ) propagator @xmath78 is the inverse in the sense of operators of @xmath79 . ( inserting the `` multilocal '' expansion in the grassmann coordinates and using the assumption of  grassmannian ultralocality " then leads to erge s for the cumulants which are just eqs . ( [ eq_flow_gamma1_ulapp],[eq_flow_gamma2_ulapp ] ) and their extension to higher orders with all copy superfields taken as equal in both sides of the equations . )    the flow of the left - hand side of eq . ( [ eq_flow_ward ] ) can then be expressed as @xmath80\\ & = \frac{1}{2}\tilde{\partial}_t \int_{\underline{x}_1\,\underline{x}_2\ ,    \underline{x}_3}\phi(\underline{x}_1)\mathcal q_{\underline{x}_1 } \left ( \mathcal p_{k;\underline{x}_2\underline{x}_3}[\phi]\gamma_{k;\underline{x}_3\underline{x}_2\underline{x}_1}^{(3)}[\phi]\right ) \\ & = \frac{1}{2}\tilde{\partial}_t \int_{\underline{x}_2\,\underline{x}_3}\mathcal p_{k;\underline{x}_2\underline{x}_3}[\phi ] \int_{\underline{x}_1}\phi(\underline{x}_1)\mathcal q_{\underline{x}_1 } \gamma_{k;\underline{x}_3\underline{x}_2\underline{x}_1}^{(3)}[\phi ] , \end{aligned}\ ] ] where we have used that @xmath81 does not act on @xmath82 . assuming that the wt identity , eq . ( [ eq_ward_k ] ) , is satisfied down to the scale @xmath42 ( we study its further evolution ) , one has @xmath83= - \left ( \mathcal q_{\underline{x}_2 } + \mathcal q_{\underline{x}_3}\right ) \gamma_{k;\underline{x}_3\underline{x}_2}^{(2)}[\phi ] , \end{aligned}\ ] ] which is obtained by differentiating twice eq . ( [ eq_ward_k ] ) . after an integration by parts and a relabel of the dummy variables @xmath84 , one then finds @xmath85\\&= -\frac{1}{2}\tilde{\partial}_t \int_{\underline{x}_2\,\underline{x}_3}\gamma_{k;\underline{x}_3\underline{x}_2}^{(2)}[\phi]\ ; \left(\mathcal q_{\underline{x}_2 } + \mathcal q_{\underline{x}_3}\right ) \mathcal p_{k;\underline{x}_2\underline{x}_3}[\phi ] . \end{aligned}\ ] ] provided that one chooses the infrared cutoff function such that @xmath86 one can replace @xmath87 $ ] by @xmath87 + \mathcal r_{k;\underline{x}_3\underline{x}_2}$ ] in eq . ( [ eq_dem_susy_dr ] ) . by using the fact that the latter is equal to the inverse modified propagator @xmath88)_{\underline{x}_3\underline{x}_2}$ ] , one can rewrite the right - hand side of eq . ( [ eq_dem_susy_dr ] ) ( up to the trivial factor of @xmath89 ) as @xmath90)_{\underline{x}_3\underline{x}_2}\ ; \mathcal q_{\underline{x}_2 } \mathcal p_{k;\underline{x}_2\underline{x}_3}[\phi]\\&=\tilde{\partial}_t \int_{\underline{x}_2\,\underline{x}_4}\delta_{\underline{x}_2\underline{x}_4}\int_{\underline{x}_3}(\mathcal p_{k}^{-1}[\phi])_{\underline{x}_3\underline{x}_4}\ ; \mathcal q_{\underline{x}_2 } \mathcal p_{k;\underline{x}_2\underline{x}_3}[\phi]\\&=\tilde{\partial}_t \int_{\underline{x}_2\,\underline{x}_4}\delta_{\underline{x}_2\underline{x}_4 } \mathcal q_{\underline{x}_2}\ ; \int_{\underline{x}_3}(\mathcal p_{k}^{-1}[\phi])_{\underline{x}_3\underline{x}_4 } \mathcal p_{k;\underline{x}_2\underline{x}_3}[\phi]\\&=\tilde{\partial}_t \int_{\underline{x}_2\,\underline{x}_4}\delta_{\underline{x}_2\underline{x}_4 } \mathcal q_{\underline{x}_2}\ ; \delta_{\underline{x}_2\underline{x}_4}. \end{aligned}\ ] ] this last expression is easily shown to be identically zero , so that one finally obtains that eq . ( [ eq_flow_ward ] ) is satisfied at scale @xmath42 . the same reasoning can be repeated with the identity associated with @xmath91 . therefore , if one starts with an initial condition that is superrotationally invariant , which can be enforced by suppressing fluctuations ( see section ii - a and paper iii ) , the above derivation proves that the wt identities associated with the superrotational invariance are _ a priori _ preserved along the rg flow if the regulator satisfies eq . ( [ eq_condition_cutoff ] ) . the latter condition can be reexpressed as @xmath92=0 , \end{aligned}\ ] ] where we have used the fact that @xmath93 and @xmath94 are translationally invariant functions of the euclidean coordinates , with therefore @xmath95 and similarly for @xmath94 . going to fourier space , this implies that @xmath93 and @xmath94 are related through @xmath96 which , as already stated , corresponds to eq . ( [ eq_cutoffs_relation ] ) with @xmath39 , a condition that we choose to enforce at the beginning of the flow when @xmath97 . such a choice of cutoff then avoids an _ explicit _ breaking of superrotational invariance . an important property of the theory is that the superrotational invariance ( which for simplicity will often be denoted by the acronym susy , for supersymmetry , in the following ) leads to dimensional reduction . we consider here the case where  grassmannian ultralocality " is satisfied so that the curvature @xmath22 drops out of the flow equations as in eqs . ( [ eq_flow_gamma1_ulapp],[eq_flow_gamma2_ulapp ] ) , and we assume that susy is indeed obeyed when the theory is restricted to one copy ( see above ) . we showed in paper iii@xcite that susy implies nontrivial wt identities relating cumulants of different orders . more specifically , the following wt identity at the scale @xmath42 will be needed below : @xmath98 & - \frac{\delta_b}{2 } ( x_1^\mu - x_2^\mu)\gamma_{k1;x_1x_2}^{(2)}[\phi ]   = \\&- \int_{x_3}\phi(x_3 ) \partial_{3\mu}\gamma_{k2;x_1x_3,x_2}^{(21)}[\phi,\phi ] . \end{aligned}\ ] ]    to prove the dimensional reduction property , it is then useful to single out a @xmath20-dimensional subspace of the @xmath4-dimensional euclidean space . we define @xmath99 with @xmath100 and @xmath101 and we consider superfields that are uniform in the 2-dimensional grassmannian and the 2-dimensional euclidean subspaces , _ i.e. _ @xmath102 . for such fields , the wt identity in eq . ( [ eq_wt_susy ] ) implies that @xmath103 = \delta_b\ , \partial_{p^2}\gamma_{k1;y_1y_2}^{(2)}[p^2;\phi],\ ] ] where @xmath43 is the @xmath20-dimensional momentum obtained from a spatial fourier transform over @xmath104 [ by using translational invariance we have as usual defined @xmath105 and similarly for @xmath106 . the same type of relation holds between @xmath107 and @xmath108 .    with the help of the above relations , the erge for the first cumulant @xmath109 $ ] , eq . ( [ eq_flow_gamma1_ulapp ] ) , can be reexpressed as @xmath110&=- \big(\frac{\delta_b}{2 } \big ) \tilde \partial_t   \int d^{d-2}y_1 \int d^{d-2}y_2 \int \frac{d^{2}p}{(2\pi)^2}\\&\times \widehat{p}_{k;y_1y_2}[p^2 , \phi]\ ; \partial_{p^2}\left ( \gamma_{k1;y_1y_2}^{(2)}[p^2,\phi ] +   \widehat{r}_k(p^2)\right ) \\&= - \big(\frac{\delta_b}{2 } \big ) \tilde \partial_t   \int d^{d-2}y_1 \int \frac{d(p^2)}{4\pi}\\&\times \partial_{p^2}\left [ \log \left ( \gamma_{k1}^{(2)}[p^2,\phi ] +   \widehat{r}_k(p^2)\right)\right ] _ { y_1y_1 } , \end{aligned}\ ] ] where , in the last expression , @xmath111 and @xmath32 are functions of @xmath112 but operators in the @xmath113-dimensional euclidean space spanned by the coordinate @xmath114 ( since @xmath115 is _ not _ uniform ) . the integral over @xmath112 is easily performed . choosing an infrared cutoff function that becomes independent of @xmath42 when its argument @xmath112 is at the uv scale , it only remains : @xmath116=&\big(\frac{\delta_b}{4\pi } \big ) \frac 12 \tilde{\partial}_t \int d^{d-2}y \\&\times \big [ \log\left ( \gamma_{k1}^{(2)}[p^2=0,\phi ] +   \widehat{r}_k(p^2=0)\right)\big ] _ { yy}. \end{aligned}\ ] ] up to the trivial factor @xmath117 and with the identifications @xmath118\equiv \gamma_{k}[\phi ] , \gamma_{k1}^{(2)}[p^2=0,\phi]\equiv \gamma_{k}^{(2)}[\phi ] ,   \widehat{r}_k(p^2=0)\equiv r_k,\ ] ] eq .  ( [ eq_dim_red ] ) coincides with the erge for the effective average action of a standard scalar @xmath119 theory@xcite in dimension @xmath5 . this provides another nonperturbative demonstration@xcite of the property of dimensional reduction for supersymmetric scalar field theories , first put forward by parisi and sourlas.@xcite      we have shown within a nonperturbative implementation of the rg that the supersymmetry of the theory , and more specifically the superrotational invariance for one copy ( susy ) , imply the property of dimensional reduction . as one knows that dimensional reduction does not hold in low enough dimension , what then goes wrong in the formalism ? the answer is that some ( super)symmetries or identities must be spontaneously broken along the rg flow . we make the assumption , which will be supported by actual calculations , that only the superrotational invariance ( for the theory restricted to a single copy ) may be broken along the rg flow , all other ( super)symmetries and properties , most significantly the `` grassmannian ultralocality '' encoding ground - state dominance ( when @xmath44 ) , remaining unaltered . from the above proof , it follows that failure of dimensional reduction implies that the wt identity between @xmath120 and @xmath111 breaks down : a singularity must occur along the flow , which invalidates the wt identity ; the equality in eq . ( [ eq_flow_ward ] ) and some step in its derivation must go wrong at some scale @xmath42 . from our previous work,@xcite we anticipate that this arises due to the appearance of a strong enough nonanalytic dependence of the second cumulant of the renormalized random field @xmath121 $ ] in its field arguments when the latter two , @xmath51 and @xmath52 , become equal . ( one should keep in mind that the cumulants are invariant under permutations of their arguments ; consequently , @xmath120 is an even function(al ) in @xmath122 . ) as a result , some higher - order derivative ( _ i.e. _ , a 1pi vertex deriving from @xmath120 ) blows up and the whole hierarchy of coupled erge s and wt identities for one copy ceases to be valid .    in order to envisage the possible scenario for spontaneous susy breaking and breakdown of dimensional reduction , it is instructive to study the structure of the exact hierarchy of coupled rg flow equations , whose lowest - order examples are given in eqs . ( [ eq_flow_gamma1_ulapp],[eq_flow_gamma2_ulapp ] ) . the following exposition is not meant to be rigorous , but only to provide some heuristic arguments . consider first the hierarchy of erge s for the proper vertices evaluated for uniform fields ( in euclidean space ) that is obtained by repeated functional differentiation of the erge for the effective average action . for this part of the reasoning , the case of a simple scalar field theory is sufficiently illustrative . the erge for the effective average action then reads @xmath123=\frac{1}{2}\int_q \partial_t r_k(q^2 ) \left(\gamma_k^{(2)}[\phi]+r_k \right)_{q ,- q}^{-1}.\ ] ] after differentiating this equation @xmath124 times , it is easily realized that for @xmath125 the erge for any generic 1pi vertex @xmath126 depends on all lower - order proper vertices ( of order @xmath127 ) as well as on some @xmath128 s and @xmath129 s . two cases then arise : when @xmath130 , the flow equation is _ linear _ in @xmath126 itself , whereas in the other case when @xmath131 , the flow equation may be _ nonlinear _ in @xmath126 itself . the same structure applies to the present theory . a consequence is that , if no other 1pi vertices diverge first , @xmath126 can only diverge at the end of the flow ( @xmath45 ) when @xmath130 ( due to the linearity of the corresponding erge ) whereas a divergence may occur at a finite scale when @xmath132 .    to get more insight , we have to look more thoroughly into the structure of the rg equations and consider the 1pi vertices associated with the cumulants separately [ see eqs . ( [ eq_flow_gamma1_ulapp],[eq_flow_gamma2_ulapp ] ) for illustration ] . we further make the following plausible assumptions ( in order to have a renormalizable theory ) :    \(i ) the @xmath133-copy 1pi vertices @xmath134 always stay finite , as do the cumulants of the renormalized random field @xmath135 and their first - order derivatives . the latter may however be discontinuous ; the condition on the first derivatives amounts to assuming that no `` supercusp''@xcite stronger than linear appears in the cumulants when two arguments become equal . \(ii ) since we expect a nonanalytic behavior occuring as the copy fields become equal in the reduction to an effective @xmath133-copy system , we also anticipate that derivatives with respect to symmetric combinations of the field arguments are always bounded ( _ e.g. _ , @xmath136 \vert_{\phi } \vert < + \infty$ ] , where the derivative is evaluated for all replica fields equal to @xmath137 ) . a `` cusp''@xcite in @xmath135 with @xmath138 implies that one of its second derivatives , which means a proper vertex of order at least @xmath139 , blows up . similarly , for having a `` subcusp''@xcite in @xmath120 , at least one of its third derivatives , hence a proper vertex of order at least @xmath139 , should blow up . as we have seen above that the associated flow equations are linear , we conclude that the divergence of such proper vertices ( with @xmath130 ) can not be the triggering event for a singularity at a finite scale . on the other hand , due to the nonlinearity of the erge s for the second derivatives of @xmath120 , a `` cusp '' may appear in the latter at a finite scale .    to summarize this discussion : one expects that breakdown of dimensional reduction requires the presence of a _ cusp _ in the field dependence of @xmath120 , cusp that should appear _ at a finite scale _ , which by analogy with what occurs for a manifold in a random environment@xcite we call the `` larkin '' scale , during the rg flow . by contrast , weaker nonanalyticities in @xmath120 as well as nonanalyticities in higher - order cumulants ( if no cusp has appeared in @xmath120 ) can only appear at the fixed point , in the limit @xmath140 , thereby preserving the dimensional - reduction property . indeed , the derivation leading to eq . ( [ eq_dim_red ] ) then remains valid . another nontrivial question is then raised : if superrotational invariance ( susy ) is spontaneously broken , _ how can one continue the rg flow for the effective average action ? _ the answer , which again should be verified in calculations , is as follows :    \(i ) keeping in mind that superrotational invariance is explicitly broken in the presence of multiple copies ( see section  [ sec : model ] ) and that it is only recovered in the process of reducing the whole problem to a single - copy system , spontaneous breaking may then be described by rejecting any implicit assumption of analyticity of the field dependences in the latter process . one should therefore restrict the hierarchy of erge s to those equations for cumulants that are considered for generic ( nonequal ) field arguments , so that a putative nonanalytic dependence in these arguments can freely emerge along the rg flow.@xcite    \(ii ) one assumes that except for the superrotational invariance , all of the properties and symmetries of the effective average action @xmath141 remain valid . in particular , this applies to the property of  grassmannian ultralocality " of the random generating functional and to its consequences for the cumulants of the renormalized disorder . as discussed in the preceding paper , this property is obeyed when taking the limit @xmath46 in the erge s for the cumulants , which leads to eqs . ( [ eq_flow_gamma1_ulapp],[eq_flow_gamma2_ulapp ] ) and their extensions to higher orders . it keeps track of the fact that only the ground state is considered for each copy . the validity of the property of  grassmannian ultralocality " is therefore distinct from that of superrotational invariance , a distinction that can not be made in the original parisi - sourlas formalism . \(iii ) the cutoff functions satisfy the relation in eq . ( [ eq_cutoffs_relation ] ) which contains as a special case the explicit superrotational invariance when @xmath39 . in practice , the typical strength of the renormalized random field @xmath36 is defined from @xmath142 evaluated for some specific uniform field configuration and the field renormalization constant @xmath37 is obtained from @xmath143 also evaluated for some uniform field configuration . finally , the cusp in @xmath144 , which is associated with spontaneous susy breaking , must be stable upon further evolution with the infrared scale @xmath42 and must not generate `` supercusps '' . from the structure of the erge for @xmath144 , it is expected that a `` linear cusp''@xcite satisfies these requirements and provides a mechanism for dimensional reduction _ if _ it persists at the fixed point ( see appendix [ appendixa ] ) . up to this point , many of the previous considerations remain plausible conjectures . we now provide an implementation of the np - frg that allows us to check their validity . the hierachy of erge s of course can not be solved exactly and the next step is to provide a susy - compatible nonperturbative approximation scheme . from our previous work,@xcite we know that an efficient scheme relies on a joint truncation in the _ derivative expansion _ , which approximates the long - distance behavior of the 1pi vertices,@xcite and in the _ expansion in cumulants of the renormalized disorder_. the truncations however must be combined in a way that does not explicitly break the supersymmetry , more precisely the superrotational invariance ( susy ) . to implement this requirement , we use the wt identities .    as is clear from the discussion in the preceding section , a minimal truncation must at least include the second renormalized cumulant . when susy is not broken , the wt identity [ see eq . ( [ eq_wt_susy_vertex ] ) ] imposes that for a uniform field @xmath145 is given by the derivative of @xmath146 with respect to @xmath147 . ( a similar relation is satisfied by the infrared cutoff functions . ) consider now the derivative expansion . in the lowest order ( lpa for `` local potential approximation '' ) , the derivative of @xmath111 is a field - independent constant , which would force @xmath120 to be also field independent : this is clearly too crude to describe the physics of the rfim . the first order of the derivative expansion leads to @xmath146 of the form @xmath148 , with @xmath149 the `` effective average potential '' ( `` gibbs free energy '' for a magnetic system ) and @xmath150 a field renormalization function . when the wt identity is satisfied , this in turn requires that @xmath151 which corresponds to a local approximation for the second cumulant . it is easy to generalize this reasoning by taking into account the whole set of wt identities for the proper vertices obtained from the cumulants : a susy - compatible approximation at order @xmath152 consists of taking @xmath153 at the order @xmath152 of the derivative expansion , @xmath154 at the order @xmath155 of the derivative expansion , ... , @xmath156 in the local approximation , and all higher - order cumulants equal to zero.@xcite this provides a scheme of successive approximations that in principle can be used to check the robustness of the results and if necessary improve them . the minimal truncation that can already describe the long - distance physics of the rfim and does not explicitly break susy is then the following : @xmath157=\int_{x}\left [ u_k(\phi(x))+\frac{1}{2}z_k(\phi(x))(\partial_{\mu}\phi(x))^2 \right ] , \end{aligned}\ ] ] @xmath158=\int_{x}v_k(\phi_1(x),\phi_2(x)),\ ] ] with the higher - order cumulants set to zero.@xcite    inserted in the erge s for the cumulants , eqs . ( [ eq_flow_gamma1_ulapp ] ) and ( [ eq_flow_gamma2_ulapp ] ) , the above ansatz provides @xmath159 coupled flow equations for the @xmath133-copy effective average potential @xmath149 ( or its derivative ) that describes the thermodynamics of the system , the field renormalization function @xmath150 , and the @xmath20-copy effective average potential @xmath160 from which one obtains the second cumulant of the renormalized random field at zero momentum , @xmath161 . susy is obeyed when @xmath162 , which is easily satisfied at the uv scale @xmath97 .      in order to search for the fixed point that controls the critical behavior , the flow equations must be recast in a scaled form . this can be done by introducing appropriate scaling dimensions ( see refs . [ ] and paper iii@xcite ) . near a zero - temperature fixed point,@xcite one has the following scaling dimensions : @xmath163 and a renormalized temperature is introduced as @xmath164 , with @xmath165 and @xmath166 related through @xmath167 . more precisely , we define running anomalous dimensions @xmath168 and @xmath169 as @xmath170 @xmath171 where @xmath172 and @xmath36 have been introduced in sec . one also has @xmath173 the dimensionless counterparts of @xmath174 will be denoted by lower - case letters , @xmath175 . the resulting equations are @xmath176 \big [ z_{k}(\varphi )    -\delta_{k,0}(\varphi)\big ]   + l_1^{(d)}(\varphi)\big[-\frac{d}{2 }    z_k'(\varphi ) + \delta_{k,0}'(\varphi)\big]\\&+\frac{1}{2}(\eta_k-\bar\eta_k)\big[z_k'(\varphi)j_2^{(1,d+2)}(\varphi)+u'''_k(\varphi)j_2^{(1,d)}(\varphi)\big ] \big\ } , \end{split}\ ] ]    @xmath177 + 4 l_{3}^{(d)}(\varphi ) z_{k}'(\varphi ) u'''_k(\varphi)[z_{k}(\varphi ) \\&- \delta_{k,0}(\varphi ) ] - l_{2}^{(d)}(\varphi ) \big[\frac{2d+1}{2 } z_{k}'(\varphi)^2 + z_{k}''(\varphi)(z_{k}(\varphi ) - \delta_{k,0}(\varphi ) ) \\ & - 2 z_{k}'(\varphi ) \delta_{k,0}'(\varphi)\big ] + m_{4}^{(d-2)}(\varphi)u'''_k(\varphi)^2 + \frac{2}{d}m_{4}^{(d)}(\varphi ) u'''_k(\varphi)\times \\&[(d+2)z_{k}'(\varphi)- 2 \delta_{k,0}'(\varphi ) ] + \frac{1}{d}m_{4}^{(d+2)}(\varphi)z'_k(\varphi)[(d+4)z_{k}'(\varphi)- \\ & 4 \delta_{k,0}'(\varphi ) ] - \frac{8}{d}[m_{5}^{(d+4)}(\varphi ) z'_k(\varphi)^2 + 2 m_{5}^{(d+2)}(\varphi ) z'_k(\varphi)u'''_k(\varphi ) + \\ & m_{5}^{(d)}(\varphi ) u'''_k(\varphi)^2][z_{k}(\varphi ) - \delta_{k,0}(\varphi ) ] -\frac{1}{2}(\eta_k-\bar\eta_k)\big[z_k''(\varphi)\times\\&j_2^{(1,d)}(\varphi)-(4+\frac 2d)z_k'(\varphi)^2j_3^{(1,d+2)}(\varphi)-4z_k'(\varphi)u'''_k(\varphi)\times\\ & j_3^{(1,d)}(\varphi)+\frac 2dz_k'(\varphi)^2h_4^{(d+2)}(\varphi)+\frac 4dz_k'(\varphi)u'''_k(\varphi)h_4^{(d)}(\varphi)+\\ & \frac 2d u'''_k(\varphi)^2h_4^{(d-2)}(\varphi)\big ] \big\ } , \end{split}\ ] ]    @xmath178 ^ 2 + l_{2}^{(d)}(\varphi_1)v_k^{(20)}(\varphi_1,\varphi_2)\\ & \times [ \delta_{k,0}(\varphi_1)-z_k(\varphi_1 ) ] +   \frac{d-2}{2 } l_{1}^{(d-2)}(\varphi_1)v_k^{(20)}(\varphi_1,\varphi_2 ) + \\ & \frac{1}{2 } m_{1,1}^{(d-2)}(\varphi_1,\varphi_2 ) - n_{1,1}^{(d-2)}(\varphi_1,\varphi_2)[\delta_k(\varphi_1,\varphi_2)-z_k(\varphi_1 ) ] \\ & -\frac{1}{2}(\eta_k-\bar\eta_k ) \big[\delta_k(\varphi_1,\varphi_2)j_{1,1}^{(1,d)}(\varphi_1,\varphi_2)+j_{1,1}^{(2,d)}(\varphi_1,\varphi_2)\\&+v_k^{(20)}(\varphi_1,\varphi_2)j_{2}^{(1,d)}(\varphi_1)\big ] + perm(12 ) \big\ } , \end{split}\ ] ]    where @xmath179 , a prime denotes a derivative with respect to the field ( when only one argument is present ) , @xmath180 , @xmath181 , and @xmath182 denotes the terms obtained by permuting @xmath183 and @xmath184 .    the functions @xmath185 , @xmath186 , @xmath187 , @xmath188 , @xmath189 , @xmath190 , @xmath191 and @xmath192 appearing in the above flow equations are `` dimensionless threshold functions '' which are defined from the infrared cutoff functions that satisfy@xcite @xmath193 and @xmath194 with the prime again denoting a derivative . the dependence of the dimensionless threshold functions on the dimensionless fields comes through @xmath195 and @xmath196 . their definitions are given in appendix [ appendixb ] . the properties of these threshold functions have been extensively discussed.@xcite they decay rapidly when their arguments @xmath197 become large , which , since @xmath198 is the square of a renormalized mass , ensures that only modes with a mass smaller than @xmath42 contribute to the flow in eqs . ( [ eq_u_ising]-[eq_v_ising ] ) . these threshold functions essentially encode the nonperturbative effects beyond the standard one - loop approximation.@xcite    note that @xmath199 and @xmath196 are even functions of @xmath8 and that , due to the @xmath200 symmetry and the permutational symmetry of the cumulants in their field arguments , @xmath201 satisfies @xmath202 , and similarly for @xmath203 . finally , the initial conditions for eqs . ( [ eq_eta_z],[eq_etabar_delta],[eq_u_ising],[eq_z_ising],[eq_v_ising ] ) at the uv scale @xmath40 ( @xmath204 ) are obtained from the bare action in eq . ( [ eq_superaction_multicopy ] ) for @xmath38 and @xmath205 , _ i.e. _ , @xmath206 the last equality being valid in the ( relevant ) region where the bare potential has a nontrivial minimum;@xcite after a trivial rescaling of the fields so that the explicit dependence on @xmath207 disappears , one also has @xmath208 which implies @xmath209 . this initial condition trivially satisfies @xmath210 which ensures that invariance under superrotations is obeyed at the beginning of the flow . an important property of the present theory occurs in the limit @xmath211 when @xmath203 is regular enough . for convenience , we introduce the variables @xmath212 and @xmath213 . the property can be stated as follows : if the behavior of @xmath214 is regular enough when @xmath215 , _ i.e. _ , is such that @xmath216 with @xmath217 when @xmath215 , then the flow of @xmath218 coincides with that of @xmath196 . this is precisely the wt relation in eq . ( [ eq_wt_susy_vertex ] ) when evaluated at zero momentum ; indeed , in dimensionless form , this equation is just @xmath219 . within the present ansatz for the effective average action , dimensional reduction then follows .    to prove the equality of the two flow equations , we first derive eq . ( [ eq_v_ising ] ) with respect to @xmath183 and @xmath184 to obtain the rg equation for @xmath214 , change variables to @xmath220 and take the limit @xmath221 by assuming the above small @xmath114 behavior [ eq . ( [ eq_deltak_small_y ] ) ] . the derivation makes use of the properties of the dimensionless threhold functions that are summarized in appendix [ appendixb ] . this leads to @xmath222 \\ & - \frac{2d-3}{2 } l_{2}^{(d)}(\varphi ) \delta_{k,0}'(\varphi)^2   + m_{4}^{(d-2)}(\varphi)u'''_k(\varphi)^2 + 2 m_{4}^{(d)}(\varphi ) \times \\&u'''_k(\varphi)\delta_{k,0}'(\varphi ) + m_{4}^{(d+2)}(\varphi)\delta_{k,0}'(\varphi)^2 +   l_k[\varphi;\delta_{k,0}-z_k]\big\ } , \end{split}\ ] ] where @xmath223 $ ] is a polynomial of degree two in @xmath224 $ ] and its derivatives , with coefficients that are function of @xmath8 and such that @xmath225\equiv 0 $ ] . subtracting from the above flow equation the one for @xmath226 in eq . ( [ eq_z_ising ] ) leads to an equation for the difference @xmath227 in the form @xmath228 = m_k[\varphi ; \delta_{k,0}-z_k ] - ( \bar\eta_k-\eta_k ) z_k(\varphi),\ ] ] where @xmath229 $ ] is another polynomial of degree two in @xmath224 $ ] and its derivatives such that @xmath230\equiv 0 $ ] . [ note that due to their definitions , @xmath231 and @xmath168 are equal as soon as @xmath219 . ] since the initial condition satisfies @xmath232 and @xmath233 , the only solution to eq . ( [ eq_diff_ising ] ) is @xmath234 at all scales , which provides the desired demonstration . by inserting this result in eq . ( [ eq_u_ising ] ) , one obtains @xmath235 which coincides with the flow equation for @xmath236 in the @xmath237 theory without random field in dimension @xmath5 ( use the relation @xmath238 and compare with eq . ( 4.33 ) of ref .  [ ] ) . the same exercise can be repeated for the flow of @xmath196 , thereby completing the proof of dimensional reduction in this nonperturbative approximation .    on the other hand , a spontaneous breaking of susy and of the associated wt identity occurs whenever @xmath239 diverges and @xmath214 has a cusp - like singularity in the form @xmath240 as @xmath241 . we expect that such a cusp appears at a finite `` larkin '' scale @xmath242 through the formation of a boundary layer as @xmath243 and @xmath215 in the form @xmath244^{2\varpi } + y^2},\ ] ] with @xmath245 an ( _ a priori _ unknown ) exponent , as can be checked by inserting in the flow equations . the boundary layer provides the appropriate initial conditions to further continue the flow below @xmath242 with a cuspy function . finally , it is worth comparing the above np - frg equations with those we have previously derived without the superfield formalism@xcite in a minimal approximation based on the same ansatz as in eqs . ( [ eq_ansatz1],[eq_ansatz2 ] ) . the two main differences are that :    \(i ) the present equations are directly considered at zero temperature ( and zero auxiliary temperature ) and subdominant terms proportional to @xmath246 do not appear in the beta functionals . \(ii ) more importantly , since from the previous formalism we had no insight into the requirements for avoiding an _ explicit _ breaking of the underlying superrotational invariance at the single - copy level , we had for simplicity set the infrared cutoff function @xmath94 to zero . as a result , the wt identity associated with superrotational invariance was violated at _ all _ scales . the dimensional - reduction property could therefore not be exactly recovered when it was expected to be valid.@xcite      we are interested in the long - distance physics of the rfim near its critical point and in the associated ( once unstable ) fixed point . such information can be obtained from the solution of the above coupled frg flow equations with appropriate boundary conditions . we have carried out the three following types of resolution . \(1 ) we have investigated the location of the breakdown of susy associated with the appearance of a cusp in @xmath203 . to this end , we have solved the set of coupled flow equations which is obtained when @xmath214 has a sufficiently regular ( no cusp ) behavior as @xmath247 , _ i.e. _ @xmath215 : see eq . ( [ eq_deltak_small_y ] ) . as already shown , the flow of @xmath218 then coincides with that of @xmath196 , and the np - frg equations for @xmath248 and @xmath196 form a closed set . in addition , we have solved the equation for the second derivative of @xmath249 with respect to @xmath250 , @xmath239 , which depends on @xmath248 and @xmath196 . a divergence in @xmath239 signals the appearance of a cusp in @xmath203 . the numerical task of solving this system of three coupled partial differential equation is reasonable . we discretize the problem both in the time and field directions . we use finite differences to evaluate the field derivatives that appear in the beta functions and use an explicit scheme to find the ( rg ) `` time '' evolution of the functions . \(2 ) when a cusp is present , solving in their full glory the set of coupled partial differential equations for functions of up to two field arguments is a hard task . we have discretized the field variables , using a 2-dimensional grid to solve the equation for @xmath203 . for convenience , we have switched variables from @xmath251 to @xmath252 and @xmath213 and restricted the grid to a trapezoidal region with @xmath253 and @xmath254 . observe that the flow of @xmath203 depends on @xmath255 and @xmath256 and their derivatives . in this sense , the flow equation for @xmath214 is nonlocal since its evolution depends not only on its property in a neighborhood of @xmath257 but also on its values at points faraway in the @xmath257 plane . it is therefore important to choose the trapezoidal grid such that , for all point in the grid , the information necessary for computing the evolution of the function at this point is inside the trapezoid . in the larger simulations we had 60 points in the @xmath258 direction and 50 in the @xmath114 direction . we have solved the set of coupled equations with two kinds of _ modus operandi _ :    ( 2-i ) first , for selected values of the dimension @xmath4 we have searched the `` cuspy '' fixed point by dichotomy , starting from `` cuspy '' initial conditions ( @xmath259 ) . depending on the dimension and the regulator parameters , there appears from time to time numerical instabilities . such instabilities showed up for instance in the ( physically not very interesting ) large - field region , where the threshold functions are very small . to keep the latter under control , we have therefore used slightly modified flow equations : this is discussed in appendix  [ appendix c ] . ( 2-ii ) secondly , once a proper `` cuspy '' fixed point has been identified in some given @xmath4 , we have followed this fixed point by continuously decreasing or increasing @xmath4 . this was done by looking directly for the roots of the beta functions by means of the hybrid algorithm , using as an initial guess the fixed point found at the previous iteration . we have studied the cuspy fixed point down to @xmath1 ; it is indeed difficult to go down to lower @xmath4 because the anomalous dimensions @xmath260 and @xmath166 become large in low dimensions : @xmath261 ( see also ref . \(3 ) when approaching the critical dimension @xmath262 for both susy and dimensional - reduction breaking from below , the procedures ( 2-i ) and ( 2-ii ) above become inefficient . we have used instead an expansion of the function @xmath263 in @xmath264 with , in practice , @xmath265 note that the expansion can not be applied in the close vicinity of @xmath262 . indeed , one then expects the approach @xmath266 to proceed nonuniformly in @xmath114 through a boundary layer in @xmath267 . it is easy to check that a similar boundary - layer phenomenon is found in the more easily handled 1-loop equation for the rf@xmath268 m near @xmath3,@xcite around the critical value @xmath269 , _ i.e. _ when @xmath270 .    to compute the anomalous dimensions , we have defined the amplitudes @xmath37 and @xmath36 appearing in eqs . ( [ eq_eta_z],[eq_etabar_delta ] ) by choosing the specific configuration of the fields as equal to zero . this amounts to imposing @xmath271 @xmath272 which fixes @xmath168 and @xmath169 . in addition to the anomalous dimensions @xmath260 and @xmath166 ( obtained as the fixed - point values @xmath273 and @xmath274 ) , we have computed the critical exponent @xmath275 that controls the escape of the rg flow away from the fixed point along the unstable ( relevant ) direction . this was done by diagonalizing the flow equations around the fixed point and extracting the negative eigenvalue which identifies with @xmath276 . finally , we have chosen a dimensionless cutoff function @xmath277 [ see eqs . ( [ eq_scaled_hatr ] ) and ( [ eq_scaled_wider ] ) ] of the form @xmath278 where @xmath279 [ not to be confused with @xmath280 . this function satisfies the requirements @xmath281 and @xmath282 . it is a modification of the function used in refs . the parameters @xmath283 , @xmath284 and @xmath285 can be optimized via stability considerations ( see refs . [ , , ] ) . they are varied to find a region of values for which the computed critical exponents are stable under changes of the parameters . for more details , see appendix  [ appendix c ] . our main result is a numerical confirmation of the scenario put forward above : the np - frg equations allowing us to continuously follow the critical behavior and the associated fixed point as a function of dimension @xmath4 , we find that there is a critical dimension @xmath262 separating a region @xmath286 where susy is valid all along the rg flow ( except possibly right at the fixed point ) and dimensional reduction applies from a region @xmath287 where susy is broken at a finite ir scale along the rg flow and a breakdown of dimensional reduction takes place . by using the procedure ( 1 ) detailed in the above subsection , we have numerically located this critical dimension as @xmath288 ( the precise value has a residual dependence on the chosen cutoff function ) . note that the value @xmath0 obtained here is consistent with the value found in our previous , and somewhat cruder , np - frg approach of the rf@xmath268 m when extrapolating the transition line @xmath289 down to @xmath290 ( see figure 4 of paper ii@xcite ) .    for initial conditions of the rg flow at or near the critical point , the second derivative @xmath291 blows up at a finite rg `` time '' @xmath292 for @xmath293 , whereas it stays finite up to the fixed point for @xmath286.@xcite we illustrate the difference of behavior between these two cases in fig . [ fig_larkin ] for a field configuration @xmath294 . np - frg flow of @xmath295 in the regime where susy is valid . the initial conditions at @xmath97 ( _ i.e. _ , @xmath204 ) for @xmath296 and @xmath297 , with @xmath298 , are taken at the fixed - point solution , @xmath299 and @xmath300 [ @xmath301 , and those for @xmath302 and @xmath303 are chosen as @xmath304 and @xmath305 . the upper ( color online blue ) curve corresponds to @xmath306 and one observes that @xmath295 tends to a finite fixed - point value . the lower ( color online red ) curve corresponds to @xmath307 shows a divergence at a finite rg `` time '' @xmath308 . ]    the divergence of the full function @xmath239 when @xmath287 is shown in fig . [ fig_larkin_3d ] . ( due to the @xmath200 symmetry , it is more convenient to represent the functions in terms of @xmath298 . ) np - frg flow of @xmath309 for @xmath307 . the initial conditions at @xmath97 ( _ i.e. _ , @xmath204 ) for @xmath296 and @xmath310 are taken at the fixed - point solution , @xmath299 and @xmath300 , and those for @xmath302 and @xmath303 are chosen as @xmath304 and @xmath305 . one observes that the divergence first takes place for small values of @xmath311 . ] we find , as can be anticipated from an analysis of the np - frg equations for @xmath312 , @xmath218 and @xmath239 , that the latter stays finite at @xmath313 and that its fixed - point value for @xmath314 behaves as @xmath315 , as seen in fig . [ fig_d_dr_dpf ] . fixed - point solution @xmath316 for dimensions @xmath317 . note the square - root behavior @xmath318 as a function of dimension . ]    when considering dimensions smaller than @xmath262 , one must study the full dependence of the function @xmath319 . we show in fig . [ fig_larkin_y ] the evolution of this function for @xmath320 . starting from a constant function , one can clearly see that a linear cusp in @xmath114 appears in a finite rg time @xmath321 , close to 0.7 for the case shown . np - frg flow of @xmath322 for @xmath323 . the initial conditions at @xmath97 ( _ i.e. _ , @xmath204 ) for @xmath296 and @xmath310 are taken at the fixed - point solution , @xmath299 and @xmath300 , and that for @xmath249 is chosen as @xmath324 . one observes that a linear cusp appears at a finite rg time . by construction @xmath325 along the flow . ] the appearance of a cusp along the flow leads to a breakdown of the superrotational invariance and of the associated wt identities . this is illustrated for the @xmath2 fixed point in fig  [ fig_d0z3d ] : there , @xmath326 , which implies a breaking of the wt identity in eq . ( [ eq_wt_approx ] ) . fixed - point solution in @xmath2 for @xmath300 ( solid line ) and @xmath327 ( dashed line ) . the two functions differ for a large enough field ( by construction , they coincide at @xmath320 ) . ] the different asymptotic behaviors at large @xmath311 are easily deduced from eqs . ( [ eq_z_ising],[eq_v_ising ] ) , from which we show that @xmath328 and @xmath329 . in @xmath2 , @xmath330 is very small , which implies that the function @xmath327 decreases slowly to @xmath331 . fixed - point solution @xmath332 for dimensions ranging from 3 ( steepest ) to 4.8 ( smoother curve ) . the variable @xmath311 has been rescaled by a factor @xmath333 so that all curves can be represented on the same plot . ] fixed - point solution @xmath300 for dimensions ranging from 3 ( steepest ) to 4.8 ( smoother curve ) . the variable @xmath311 has been rescaled by a factor @xmath333 so that all curves can be represented on the same plot . ] we also display in figs . [ fig_uppf ] and  [ fig_zpf ] the fixed - point solutions @xmath332 and @xmath300 for different dimensions . one can see that the lower the dimension the steeper the curves . this means that for a numerical study of the critical properties in low dimensions , we need to discretize the field dependence in the np - frg equations with a small mesh , which entails a large number of points . the numerical integration is therefore more difficult and even becomes intractable in practice . finally , we show in fig . [ fig_deltapf ] the full dependence of the fixed - point solution @xmath334 . fixed - point solution @xmath334 in @xmath2 . note that the point @xmath335 corresponds to the lower right corner . the function is even in @xmath114 ( with a cusp around @xmath221 ) and the part corresponding to @xmath336 is not shown . ] we now turn to the results concerning the critical exponents . we begin with the anomalous dimensions @xmath260 and @xmath337 which are determined at the ( critical ) fixed point . ( we recall that @xmath260 characterizes the spatial decay of the `` connected '' pair correlation function and @xmath166 that of the `` disconnected '' pair correlation function at criticality . ) their dependence on the spatial dimension @xmath4 is shown in fig . [ fig_eta ] . anomalous dimensions @xmath260 and @xmath338 as functions of the spatial dimension @xmath4 . the dashed lines are lower bounds for the anomalous dimensions [ @xmath339 and @xmath340 for @xmath260 and @xmath166 respectively ] and the crosses correspond to predictions from computer simulations and ground - state determinations . ]    above @xmath0 , we find that @xmath260 and @xmath166 rigorously coincide , @xmath341 , which is the signature of dimensional reduction . ( their value is moreover equal to that found in the same order of the derivative expansion for the pure ising model in dimension @xmath5 . ) below @xmath262 , the two anomalous dimensions bifurcate and @xmath342 . note that the predicted values satisfy the required bounds , _ i.e. _ , @xmath343 , @xmath344 , @xmath345 . we give in table  [ tab_results_eta_etab ] our estimates for @xmath260 and @xmath166 in @xmath2 and @xmath3 , where a comparison is possible with numerical determinations from computer studies ( monte carlo calculations and @xmath65 ground - state determinations ) . .[tab_results_eta_etab ] comparison between the anomalous dimensions obtained in the present work and in computer studies ( a : [ ] ; b : [ ] ; c : [ ] ; d : [ ] ; e : [ ] ; f : [ ] ; g : [ ] ; h : [ ] ; i : [ ] ; j : [ ] ; k : [ ] ; l : [ ] ; m : [ ] ; n : [ ] ) . the error bars on the results from simulations and ground - state computations can be found in the original papers.the values of the temperature exponent @xmath346 that we obtain are @xmath347 in @xmath2 and @xmath348 in @xmath3 , in excellent agreement with the other studies . [ cols="<,<,<,<",options="header " , ]     one can see that the agreement is good or very good . ( real - space rg studies in @xmath2@xcite provide results in the same range , with an anomalous dimension @xmath260 between @xmath349 to @xmath350 and an exponent @xmath165 between @xmath351 and @xmath352 . ) as is well known from the study of simpler models such as the pure @xmath268 model , the way to further improve the accuracy of the exponents would be to consider higher orders of the approximation scheme presented in section iv - a . we have not been able to perform our np - frg calculation down to the lower critical dimension @xmath353 , where the values @xmath354 and @xmath355 are exactly known . the numerical resolution of the flow equations become extremely arduous in low dimension where the anomalous dimensions become large and approach their lower bound , @xmath339 and @xmath340 respectively . in addition , one also encounters numerical difficulties as one approaches the critical dimension @xmath262 from below . indeed , as mentioned in section iv - d , it is anticipated that the limit @xmath356 , @xmath215 is nonuniform with a boundary layer in @xmath267 . more work will be needed to solve this boundary - layer problem numerically . [ fig_eta ] also provides evidence that the claim according to which the two exponents @xmath260 and @xmath166 are related by a fixed ratio , @xmath357,@xcite can not be right in general . it is true that the relation @xmath357 is exact in @xmath353 ( and according to a phenomenological rg,@xcite also at first order in @xmath358 ) and is very closely obeyed by the numerical estimates in @xmath2 and @xmath3 ( see table  [ tab_results_eta_etab ] ) . however , the latter type of `` numerical evidence '' can always be challenged , and it is actually impossible to reach a definite conclusion on the sole basis of numerical results in selected dimensions ( @xmath359 ) . on the other hand , the overall and continuous dependence on spatial dimension that we provide through the np - frg leads to a different and firmer type of answer : since one goes from @xmath341 for @xmath360 to @xmath361 above , the relation @xmath357 can not be always valid . finally , we also display our results in @xmath2 and @xmath3 for the correlation length exponent @xmath275 , which is associated with the relevant direction around the critical fixed point , in table i. the agreement with the available data from computer studies is excellent in @xmath3 and fair in @xmath2 . ( note that in the real - space rg studies in @xmath2@xcite , the exponent @xmath275 is in general too large , with values ranging from @xmath362 to @xmath363 . ) all the other critical exponents are obtained in the np - frg through the expected relations : @xmath364 , @xmath365 , @xmath366 , etc.@xcite    the validity of the theoretical description is clearly confirmed by both the overall consistency and the quantitative accuracy of the predictions in @xmath2 and @xmath3 .      to assess the effect of a nonzero auxiliary temperature and of a deviation from `` grassmannian ultralocality '' , we have considered the simplest `` non - ultralocal '' contribution to the cumulants introduced in section vii - b of paper iii . combined with the above truncation of the `` ultralocal '' contributions of effective average action [ eqs . ( [ eq_ansatz1],[eq_ansatz2 ] ) ] , this amounts to expressing the first cumulant as ( see eq . ( 136 ) of paper iii ) @xmath367 =   \int_{\underline{\theta}_1 } \bigg ( \gamma_{k1 } [ \phi_1(\underline{\theta}_1)]+\frac{1}{2\beta}(1 + \beta \bar \theta_1 \theta_1)\times \\ & \int_x y_{k}(\phi_1(x,\underline{\theta}_1))\partial _ { \theta_1 } \phi_1(x,\underline{\theta}_1)\partial_{\bar \theta_1 } \phi_1(x,\underline{\theta}_1 ) \bigg ) , \end{aligned}\ ] ] with @xmath109 $ ] given in eq . ( [ eq_ansatz1 ] ) , whereas the second - order cumulant is taken as purely `` ultralocal '' and given in eq . ( [ eq_ansatz2 ] ) . the frg flow equation for the function @xmath368 is directly obtained from eq . ( 139 ) of paper iii . in a dimensionless form , it reads @xmath369 l^{(d)}_2(\varphi)+(\bar \eta_k-\eta_k)j^{(1,d)}_2(\varphi)\big\}+ o(t_k ) ,    \end{split}\ ] ] where the threshold functions are defined in appendix  [ appendixb ] . in deriving the above expression from eq . ( 139 ) of the preceding paper , we have used the dimensionless quantity @xmath370 defined through @xmath371 and the fact that as @xmath372 , @xmath373 the above equation is valid provided that @xmath374 which asymptotically goes as @xmath375 is subdominant with respect to @xmath376 . as discuss just below , this is indeed true , although @xmath370 goes to zero at the fixed point . the flow of the `` ultralocal '' quantities in the limit @xmath46 has already been solved and we can introduce the corresponding solution in eq . ( [ eq_flow_ydimensionless ] ) . we have solved the latter equation and found that @xmath377 goes to zero as one approaches the fixed point , in such a way that the dimensionful quantity @xmath378 goes to a constant @xmath379 : this is illustrated in fig . 10 for @xmath2 .     of the dimensionless field with @xmath380 . after a transient regime , the function tends to a constant . ]    in addition , we can investigate the effect of the `` non - ultralocal '' contribution to the flow of the `` ultralocal '' quantities when the auxiliary temperature @xmath15 is different from zero . this can be done by considering the erge s for the `` ultralocal '' components of the first and the second cumulants given in eqs . ( 134 ) and ( 135 ) of paper iii . from the above results ( in particular , the fact that @xmath376 goes to a finite constant when @xmath372 ) , we know that @xmath381 is asymptotically equal to @xmath382 : see eq . ( [ eq_hatqp_approx ] ) and discussion below . the flow of the dimensionless first cumulant therefore follows eq . ( 151 ) of paper iii and the contribution of the `` non - ultralocal '' piece goes as @xmath246 times a well - behaved function(al ) , with @xmath383 and @xmath384 . the flow of the second cumulant is also of the form of eq . ( 151 ) in the preceding paper , but the dimensionless beta function due to the `` non - ultralocal '' piece is now potentially singular when a cusp appears in the limit of zero ( auxiliary ) temperature of the second cumulant . more specifically , one finds that @xmath385+perm(12 ) \big\ } , \end{split}\ ] ] where @xmath386 is the zero ( auxiliary ) temperature beta functional obtained by taking derivatives with respect to @xmath183 and @xmath184 of the right - hand side of eq . ( [ eq_v_ising ] ) and we have omitted subdominant terms coming from the difference between @xmath387 and @xmath388 ( see above ) .    we can now follow the derivation done for the effect of the ( bath ) temperature in the formalism without superfields ( see paper ii@xcite ) . after changing variables from @xmath389 to @xmath390 and @xmath391 , one obtains the solution for @xmath392 in the close vicinity of the fixed point , when both @xmath246 and @xmath114 go to zero , in the form of a boundary layer @xmath393 + o(t_k^2 ) , \end{split}\ ] ] where @xmath394 , @xmath395 , @xmath396 are well - behaved functions of @xmath19 and @xmath397 is equal to the absolute value of the coefficient of the cusp in the zero - temperature fixed - point function @xmath398 .    the connection between this boundary - layer phenomenon and the presence of rare low - energy excitations above the ground - state known as `` droplets''@xcite is the same as that already discussed in our previous work@xcite ( see also refs . [ , ] ) and the discussion is not repeated here . in this paper and the preceding one @xcite we have extended our np - frg approach of disordered systems , which was initiated in the previous articles of this series.@xcite the objective was to discuss the property of dimensional reduction and its breakdown in the rfim from the parisi - sourlas@xcite perspective of an underlying supersymmetry and its breaking . we have reformulated the np - frg in a superfield formalism and , to do so , we have proposed a solution for properly selecting the ground state among the many metastable states that are present at zero temperature . through the introduction of an appropriate regulator and a supersymmetry - compatible nonperturbative approximation , we have been able to follow the supersymmetry ( more precisely , the superrotational invariance ) and its spontaneous breaking along the rg flow .    despite the fact that the effective hamiltonian ( or microscopic action ) in the presence of a random field has numerous minima in the region of interest near the critical point , dimensional reduction need not be systematically broken . by implementing a np - frg flow for the cumulants of the renormalized disorder , our work shows that there is a finite range of dimension below the upper critical dimension , @xmath399 , for which this multiplicity of minima has no effect on the long - distance properties of the model . more precisely , the scaling behavior around the critical point conforms to the dimensional - reduction predictions . the associated fixed point is characterized by a nonanalytic dependence of the effective action in the dimensionless fields , at odds with a naive description based on perturbation theory , but the nonanalyticity is too weak to alter the spectrum of critical exponents . it is only by going below a critical dimension @xmath262 , which we numerically find around @xmath400 , that spontaneous breaking of the supersymmetry ( superrotational invariance ) takes place . the robustness of our theoretical description , which explains the pending puzzles concerning the critical behavior of the rfim , is supported by the good agreement obtained between the np - frg predictions for the critical exponents and the available values from computer studies . finally , one may wonder whether the formalism developed here can be useful in a different context . an obvious extension is a study of the long - distance physics of self - avoiding branched polymers and the associated property of reduction to the yang - lee edge singularity problem in two fewer dimensions.@xcite dimensional reduction is known to hold in this case,@xcite and this provides a benchmark model to test the ability of our nonperturbative approximation scheme to reproduce this property . another extension concerns the hysteresis behavior and out - of - equilibrium phase transitions of the rfim at zero temperature when driven by a slow change of the external magnetic field.@xcite work in those directions is under way . more challenging , and far more speculative too , is the question known as the `` gribov ambiguity '' in the nonperturbative quantization of nonabelian gauge field theories.@xcite the standard faddeev - popov gauge - fixing procedure aims at restricting the functional integral over the gauge field to nonequivalent gauge configurations ( _ i.e _ , configurations that are not obtainable one from another by a gauge transformation ) . unfortunately , the gauge - fixing conditions that respect lorentz invariance and internal ( color ) symmetry , such as the landau gauge condition , are inconsistent because they admit solutions that are equivalent up to a gauge transformation ( `` gribov copies '' ) . there has been a strong theoretical effort to overcome this problem . this `` gribov ambiguity '' does not affect the perturbative regime , so that calculations at high energy can be performed in the standard faddeev - popov approach . however , very little is known on its influence in the nonperturbative regime . the gribov - zwanziger model  @xcite enables one to reduce the number of copies taken into account in the functional , but not to a single one , and , up to now , there is no unambiguous gauge - fixing procedure . the strong connections between the faddeev - popov and the parisi - sourlas formalisms , and their common failure when the relevant field equation has a multiplicity of solutions , have been underlined.@xcite the relevance of the tools and concepts developed in the present and the preceding papers to the problem of the gribov ambiguity has yet to be investigated . to illustrate that a linear cusp in @xmath144 does not generate a supercusp and may break dimensional reduction , we consider the erge for the fourier transform of the 2-point proper vertex @xmath401 ( for uniform fields ) in the limit @xmath402 . we switch to the variables @xmath403 and @xmath404 and we study a situation in which a cusp is present : @xmath405 where @xmath406 and @xmath407 ( recall that @xmath120 is even in @xmath137 and @xmath408 separately ) . with the assumption concerning the behavior of the proper vertices detailed in sec . iii - b , we find that the flow of @xmath120 , where for simplicity we take @xmath409 , can be written in the @xmath406 limit : @xmath410 , \end{aligned}\ ] ] where @xmath411 denotes the terms that could be obtained by directly considering the flow of @xmath412 and by assuming a regular behavior of its beta - functional when @xmath406 . the contribution due to the derivatives of the third cumulant @xmath413 is an even function of @xmath408 ; the @xmath414 term contributes to the regular component of the beta - functional and the next term is a @xmath415 ) with @xmath416 . the above equation indicates that a `` supercusp '' with @xmath417 leads to an ill - defined flow for @xmath412 . for @xmath418 ( cusp ) one should also study the erge for the third cumulant to determine the leading nonanalyticity ( _ i.e. _ , @xmath419 ) and therefore study the stability of the cusp under further rg flow . in any case , eq . ( [ eq_flow_cusp_stab ] ) shows that a cusp weaker than linear , with @xmath420 , has no effect on the flow of @xmath412 and , since @xmath421 is obtained from the latter , is not expected to modify dimensional reduction . a linear cusp ( @xmath422 ) on the other hand provides a possible mechanism for the failure of dimensional reduction while not _ a priori _ generating stronger nonanalyticities . in order to actually lead to breakdown of dimensional reduction , the linear cusp must of course remain in the renormalized second cumulant up to the appropriate fixed point . the various dimensionless threshold functions correspond to the various @xmath133-loop integrals involving the infrared cutoff functions and the propagators ( after account of the scaling dimensions ) . more specifically , they are defined as follows:@xcite @xmath423 @xmath424 @xmath425 with @xmath426 , @xmath279 , and @xmath427 acts only on the dimensionless function @xmath428 that is contained in the @xmath429s@xcite(by definition , @xmath430 $ ] ) . note that , in addition to the dependence on @xmath431 and @xmath432 , the threshold functions explicitly depend on the scale @xmath42 via the running anomalous dimension @xmath168 . ( note also that the functions @xmath433 and @xmath434 are not symmetric in the exchange of the indices @xmath133 and @xmath20 . ) the generic properties of these dimensionless threshold functions are discussed in detail in ref . [ ] .    in the present study , it is also necessary to introduce additional threshold functions , which accounts for the fact that the cutoff function @xmath435 has an anomalous scaling compared to @xmath436 when dimensional reduction is broken . these threshold functions then always appear multiplied by @xmath437 in the flow equations . we define : @xmath438 and @xmath439    we summarize here a number of relations that are useful in the developments of sec . one has @xmath440 , @xmath441 , @xmath442 , @xmath443 , with @xmath444 . in addition , the functions @xmath445 s and @xmath446 s are related by @xmath447    from the above definitions , one also straightforwardly finds that    @xmath448    @xmath449    and    @xmath450 @xmath451    @xmath452    similar relations hold for the functions @xmath453 and @xmath454 . from these expressions one can obtain the derivatives with respect to the field arguments @xmath455 by using that @xmath431 and @xmath456 . dependence of the anomalous dimensions @xmath260 ( upper panel ) , @xmath166 ( middle panel ) and of the critical exponent @xmath275 ( lower panel ) on the parameter @xmath283 of the ir cutoff function in eq . ( [ eq_chosen_r ] ) . the `` principle of minimal sensitivity '' leads to a determination of the anomalous exponents @xmath457 , @xmath458 , @xmath459 ( at the minima ) . the variations of the anomalous dimensions are very small over a wide range of @xmath283 and enables us to evaluate the precision to @xmath460 and @xmath461 . @xmath275 is more sensitive to the regulator parameters and we estimate @xmath462,title="fig : " ]   dependence of the anomalous dimensions @xmath260 ( upper panel ) , @xmath166 ( middle panel ) and of the critical exponent @xmath275 ( lower panel ) on the parameter @xmath283 of the ir cutoff function in eq . ( [ eq_chosen_r ] ) . the `` principle of minimal sensitivity '' leads to a determination of the anomalous exponents @xmath457 , @xmath458 , @xmath459 ( at the minima ) . the variations of the anomalous dimensions are very small over a wide range of @xmath283 and enables us to evaluate the precision to @xmath460 and @xmath461 . @xmath275 is more sensitive to the regulator parameters and we estimate @xmath462,title="fig : " ]   dependence of the anomalous dimensions @xmath260 ( upper panel ) , @xmath166 ( middle panel ) and of the critical exponent @xmath275 ( lower panel ) on the parameter @xmath283 of the ir cutoff function in eq . ( [ eq_chosen_r ] ) . the `` principle of minimal sensitivity '' leads to a determination of the anomalous exponents @xmath457 , @xmath458 , @xmath459 ( at the minima ) . the variations of the anomalous dimensions are very small over a wide range of @xmath283 and enables us to evaluate the precision to @xmath460 and @xmath461 . @xmath275 is more sensitive to the regulator parameters and we estimate @xmath462,title="fig : " ]    in the exact formulation of the np - frg , the results of the flow equations are independent of the ir regulator , which is only an intermediate means that does not affect the long - distance physics . once approximations are introduced , a residual dependence on the choice of regulator however remains . ( a similar situation occurs in perturbation theory where a residual dependence on the parameters of the borel resummation procedure is observed . ) one can in some sense optimize the choice of ir cutoff function by demanding that the output , say the critical exponents , satisfies a property of `` minimal sensitivity '' such that by varying the characteristics of the cutoff function around the optimum , a minimal variation of the exponents results . this guarantees the stability of the results . we apply this procedure to the type of function given in eq . ( [ eq_chosen_r ] ) where the parameters @xmath283 , @xmath284 , and @xmath285 can be varied . we illustrate in fig . [ fig_etapms ] the dependence of the anomalous dimensions in @xmath2 on the parameter @xmath283 when @xmath284 and @xmath285 are kept fixed at the values @xmath463 and @xmath464 respectively ( these latter values were determined by preliminary variational studies to find minimal sensitivity ) . one can see that there is a large domain of values in which @xmath260 and @xmath166 vary little ( say by less than @xmath465 ) . we use the `` principle of minimal sensitivity '' to determine the best estimate of the critical exponents [ note that the minima of @xmath466 and @xmath467 happen at very close values of @xmath283 ] . we also estimate the error bars from the observation that when @xmath283 varies in the range [ 1.5,3 ] , the maximum variation of @xmath260 is of 0.03 and that of @xmath338 is of 0.04 . these are the values that we have reported in sec . [ sec : results ] .      during the numerical integration of the flow , we encountered numerical instabilities in physically unimportant regions ( typically , at large fields ) . the origin of these instabilities can be understood as follows . in the region of large fields , the threshold functions rapidly decrease to zero . in this limit , the flow of a function is given by the dimensional part , _ i.e. _ the flow obtained by setting @xmath469 [ _ e.g. _ @xmath470 . this flow no longer depends on the second derivative of the associated function with respect to the field . however , this diffusive - like term [ @xmath471 in our previous example ] is very important to stabilize the numerical integration . in practice , when necessary , we have therefore added by hand a small diffusive contribution to the flow equations and checked that the numerical results are stable under varying the strength of these extra terms . the physics associated with  griffiths phases " is however unlikely to be captured by such an approach . in this case the rare regions that are involved are exponentially suppressed as their size increases and produce very weak ( essential ) singularities in the thermodynamics . such a behavior does not show up easily once the average over disorder has been performed , as done in the cumulants . droplets on the other hand are `` power - law rare '' regions in their size and have a clear signature in the functional dependence of the cumulants of the renormalized disorder ( see _ e.g. _ the detailed work by balents and ledoussal on the random manifold model.@xcite ) the same is true for  avalanches " , which are collective discontinuous events occurring in _ typical _ samples at zero temperature for exceptional values of the applied source . we define `` supercusp '' , `` cusp '' , and `` subcusp '' in @xmath472 according to the order of the nonanalyticity in @xmath473 when @xmath402 : @xmath474 corresponds to a supercusp , @xmath475 to a `` cusp '' , and @xmath419 a noninteger @xmath476 to a subcusp . a `` linear cusp '' is associated with @xmath477 . this can be generalized to higher - order cumulants of the renormalized random field . note that there does not seem to be any alternative route , for instance by introducing new sources conjugate to operators that break susy explicitly.we have not found any useful additional wt identities when the invariance under superrotations is broken . to check this , we have followed zinn - justin s procedure,@xcite which amounts to performing transformations on the action and adding sources to the new generated independent operators . however , introducing new operators to compensate for the terms breaking susy breaks additional symmetries , including rotations in euclidean subspace , and an increasing number of independent operators are generated at each new iteration of the procedure . if one also wishes to describe microscopic details associated with a specific definition of the bare action and keep a finite uv cutoff @xmath40 , one has to generalize the function @xmath277 to @xmath478 , with @xmath479 when @xmath480 and @xmath481 ; a possible choice of @xmath482 is @xmath483 where @xmath484 is the heaviside step function . when @xmath485 for @xmath147 fixed , the infrared cutoff function @xmath486 diverges as @xmath487 whereas @xmath488 goes to @xmath489 , which fulfills the requirements described in paper iii.@xcite    the constraint of no explicit susy breaking and the associated wt identities are not easily implemented in other approximation schemes . for instance one could envisage to replace the derivative expansion by a more powerful scheme to describe the full momentum dependence of the proper vertices , as proposed in ref . [ ] . however , it is not obvious then to combine this approximation with the expansion in cumulants . note that if the bare random field is not gaussian distributed , one can not bluntly neglect all higher - order cumulants . one should then check that , as expected , the non - gaussian behavior that is present at the bare level does not change the long - distance physics . to assess this point , it is necessary to include in the description at least the third cumulant and then consider the next order of the approximation scheme . when @xmath317 , one anticipates a subcusp in @xmath490 as @xmath491 , with no implication for the dimensional reduction . numerically however , it is very hard to find evidence for such subcusps as this would involve computing high orders in the derivatives with respect to the fields : for instance , one should be able to follow the fourth derivative of @xmath214 with respect to @xmath114 and check for its divergence at the fixed point above @xmath262 . one can also check that the `` magnetic - field eigenvalue '' @xmath492 corresponding to the relevant direction associated with odd eigenfunctions in the field arguments is analytically obtained in any dimension as @xmath493 , which as expected is the dimension of the fields ; the corresponding eigenfunctions @xmath494 , @xmath495 , and @xmath496 are expressed in terms of the fixed - point functions , namely : @xmath497 , @xmath498 , and @xmath499 .
Abstract:

Answer: we apply the nonperturbative functional renormalization group ( np - frg ) in the superfield formalism that we have developed in the preceding paper to study long - standing issues concerning the critical behavior of the random field ising model . through the introduction of an appropriate regulator and a supersymmetry - compatible nonperturbative approximation , we are able to follow the supersymmetry , more specifically the superrotational invariance first unveiled by parisi and sourlas [ phys . rev . lett . * 43 * , 744 ( 1979 ) ] , and its spontaneous breaking along the rg flow . breaking occurs below a critical dimension @xmath0 , and the supersymmetry - broken fixed point that controls the critical behavior then leads to a breakdown of the `` dimensional reduction '' property . we solve the np - frg flow equations numerically and determine the critical exponents as a function of dimension down to @xmath1 , with a good agreement in @xmath2 and @xmath3 with the existing numerical estimates .