arxiv:7.0597v [cond-mat.stat-mech] 6 Dec 07 Diagram Reduction in Problem of Critical Dynamics of Ferromagnets: 4-Loop Approximation L. Ts. Adzhemyan, E. V. Ivanova, M. V. Kompaniets and S. Ye. Vorobyeva Saint-Petersburg State University, 7-9 Universitetskaya nab. Saint-Petersburg, Russian Federation E-mail: l.adzhemyan@spbu.ru, ellekspb@gmail.com, m.kompaniets@spbu.ru svetlana.e.vorobyeva@gmail.com, Abstract. Within the framework of the renormalization group approach to the models of critical dynamics, we propose a method for a considerable reduction of the number of integrals needed to calculate the critical exponents. With this method we perform a calculation of the critical exponent z of model A at 4-loop level, where our method allows to reduce number of integrals from 66 to 7. The way of constructing the integrand in Feynman representation of such diagrams is discussed. Integrals were estimated numerically with Sector Decomposition technique. Keywords: renormalization group, ε-expansion, multi-loop diagrams, critical exponents
Diagram Reduction in Critical Dynamics. Introduction Model A of critical dynamics describes critical slowing down effect for systems with non-conserved order parameter [,, 3]. Usually this model is used as a theoretical model for critical behavior of ferromagnets [4]. Recently the new classes of materials were investigated and it was found that this model describes phase transitions in multiferroics [5] and in the systems with ordering phase transitions [6] as well. An additional motivation to study model A is also the fact that it is the simplest model of critical dynamics and new technical methods can be tested on it. Despite the fact that renormalization group is one of the well acknowledged theoretical methods for investigation of continuous phase transitions, the application of this method to the problems of critical dynamics faces much greater difficulties in comparison with the problems of critical statics. The analytic results obtained here are limited in the best case to the third order of perturbation theory [7], whereas in the static ϕ 4 theory the six-loop result [8, 9, 0] is currently reached, and for the anomalous dimension of the field seven-loop []. In problems of critical dynamics a noticeable lag occurs also in numerical calculations, in which the Sector Decomposition technique of the calculation of Feynman diagrams [] proved to be very effective in critical statics problems (5 loops and partially 6 loops in the theory ϕ 4 [0, 3]), while in the dynamic problems this method has so far been used only in the two-loop approximation [4]. The calculation of multi-loop diagrams in critical dynamics, taking into account their complexity, requires considerable time, thus the problem of reducing the number of calculated diagrams arises. In this paper, a numerical calculation of the renormalization group functions of the model A [, ] is performed in the fourth order of perturbation theory. We present a method that allows to reduce significantly the number of Feynman diagrams to be calculated ( reduction of diagrams) by appropriate grouping of the original diagrams of the theory. The rules for constructing these diagrams and their integrand in the Feynman representation directly from the graph are formulated. The subsequent numerical four-loop calculation was carried out using the Sector Decomposition method. The paper is organized as follows: in section we recall renormalization procedure for the model A with the use of dimensional regularization (d = 4 ε) and minimal subtraction scheme (MS). In the next section we present the diagrammatic representation for model A. In section 4 we describe the method of diagram reduction. In subsequent section we present the four loop results for the dynamical critical exponent z, which are followed by conclusion. In Appendix A we present one more example of the diagram reduction and in Appendix B we discuss the Feynman representation for diagrams in models of critical dynamics.
Diagram Reduction in Critical Dynamics 3. Renormalization of the model Non-renormalized action of model A of critical dynamics in the space with the dimension d = 4 ε is determined by the set φ 0 of two non-renormalized fields φ 0 {ψ 0, ψ 0} and has the form [3]: S 0 (φ 0 ) = λ 0 ψ 0ψ 0 + ψ 0[ t ψ 0 + λ 0 δs st 0 /δψ 0 ] = = λ 0 ψ 0ψ 0 + ψ 0[ t ψ 0 + λ 0 ( ψ 0 τ 0 ψ 0 3! g 0ψ 3 0)] () with non-renormalized static action S st 0 (φ 0 ) = ( ψ 0 ) / τ 0 ψ 0/ 4! g 0ψ 4 0. () Renormalized action S R obtained by multiplicative renormalization of parameters and fields can be represented as the sum S R = S B + S of basic action S B and counterterms S [3]: S B = λψ ψ + ψ [ t ψ + λ( ψ τψ 3! µε gψ 3 )], (3) S R = Z λψ ψ + ψ [ Z t ψ + λ(z 3 ψ Z 4 τψ 3! Z 5µ ε gψ 3 )], (4) where λ 0 = λz λ, τ 0 = τz τ, g 0 = gµ ε Z g, ψ 0 = ψz ψ, ψ 0 = ψ Z ψ, (5) Z = Z λ Z ψ, Z = Z ψ Z ψ, Z 3 = Z ψ Z λ Z ψ, (6) Z 4 = Z ψ Z λ Z τ Z ψ, Z 5 = Z ψ Z λ Z g Z 3 ψ. It follows from the multiplicative renormalizability of the models (), (), that the renormalization constants Z ψ, Z τ, Z g in this model coincide with the static ones (i.e. of the model ()) Z ψ = (Z ψ ) st, Z τ = (Z τ ) st,, Z g = (Z g ) st (7) and the relation Z ψ Z λ = Z ψ fulfilled [3]. This means that the renormalization constants Z 3, Z 4, Z 5 are purely static and Z = Z. (8) The only new renormalization constant is Z λ = Z Z ψ = Z Z ψ. (9) For our purposes, it is convenient to calculate it through the renormalization constant Z, which is determined from the diagrams of the -irreducible function Γ ψ ψ = ψ ψ irr /(λ) on the zero external frequency ω and the momentum p. For the expansion of this function into a perturbation theory series, we will use the coupling constant u = S d g, where S (π) d d = πd/ is the area of the d-dimensional unit sphere. This Γ(d/) expansion has the form Γ R ψ ψ ω=0,p=0 = Z ( + u Zg (µ /τ) ε Zτ ε A () + u 3 Zg 3 (µ /τ) 3ε/ Zτ 3ε/ A (3) +...). (0)
Diagram Reduction in Critical Dynamics 4 For renormalization of the model (4) we will use the minimal subtraction (MS) scheme where counterterms contain only poles in ε. The renormalization constants Z g and Z τ in (0) are known from the statics, while Z at the 4-loop approximation in the MS scheme has the following form Z = + z ε u + ( z3 ε + z 3 ε ) ( u 3 z43 + ε + z 4 3 ε + z ) 4 u 4 + O(u 5 ), () ε coefficients z nk can be found from the condition of the absence of poles in ε in the function Γ R ψ ψ ω=0,p=0, thus the main technical problem is to calculate the coefficients A (), A (3), A (4) in (0). 3. Diagrammatic representation after integration over internal time variables Propagators of the model (3) in the time-momentum (t, k) representation have the form: t t = ψ(t )ψ(t ) = E k exp λe k t t, t t = ψ(t )ψ (t ) = θ(t t ) exp λe k(t t ), t t = ψ (t )ψ (t ) = 0, () where E k k + τ. The simple exponential dependence of these propagators on time makes it easy to integrate diagrams in (t, k) representation over internal time variables and reduce the problem to integration over momenta (momentum representation). The result of integration over time can be expressed in a diagram language using the technique of time versions (see, for example, [3]). Let us remind this technique with the following diagram, considered at zero external frequency ω: t0 t t ω=0. (3) Taking into account the θ function in the propagator ψψ (), the domain of integration in (3) can be represented as 3 contributions (time versions) (t 0 > t > t ) + (t 0 > t > t ) + (t > t 0 > t ). (4) Explicit integration over internal times for each time version (4) can be represented with new diagrammatic technique: ω=0 = 0 + 0 + 0. (5)
Diagram Reduction in Critical Dynamics 5 In this diagrammatic technique we associate factor /E ki with each solid line ( ψψ ), factor with each solid line with dash ( ψψ ) and factor / i E k i with dotted line, the last sum is going over all energies () of the diagram lines, which are crossed by the dotted line, where k i are the momenta of the corresponding lines. Than integrand in the momentum representation for the first diagram on the righthand side (5) has the form: 3 4 5 E E E 4 E 5 (E + E + E 3 ) (E + E 4 + E 5 ). (6) Here and in the following, we denote E i E ki. The integrands of the remaining time versions, shown in the figure (5), are constructed in a similar way. 4. Reduction of diagrams A complicating circumstance in the problems of critical dynamics, in comparison with the static case, is a significantly larger number of momentum integrals arising as a result of integration over time (and a more complicated form of them). In a number of papers [5, 6, 7, 8] the fact that, turning to certain sums of diagrams, one can appreciably simplify the integrands was used. We propose a systematic procedure for such a reduction of diagrams, that makes it possible to automate the calculations, which is necessary for calculations in the higher orders of perturbation theory. The possibility of such a reduction is actually seen from the relation (7). The equality (7) of static and dynamic counterterms is a consequence of a more general statement about the coincidence of -irreducible static functions ψψ irr st, ψψψψ irr st and dynamic functions ψ ψ irr, ψ ψψψ irr at zero frequency: ψ ψ irr ω=0 = ψψ irr st, ψ ψψψ irr ω=0 = ψψψψ irr st. (7) In the diagram language the equalities (7) mean that for these functions the sum of the dynamic diagrams is reduced to a simpler object to the sum of the static diagrams. We will consider examples of the technical implementation of such a procedure, and then apply similar techniques to simplify the function of interest to us ψ ψ irr ω=0. Let us prove the equality: ω=0 =, (8) 6 where and are symmetry coefficients of diagrams. Performing integration over time, 6 we can write the left-hand side as 3 ω=0 = 3 st 0 E E 3 (E + E + E 3 ) (9)
Diagram Reduction in Critical Dynamics 6 Here the digits denote the integration momenta, and the integrand stands on the righthand side of the formula. Symmetrizing this expression with respect to the momenta of integration, we obtain E E 3 (E + E + E 3 ) (0) ( + + ) 6 E E 3 E E E E 3 (E + E + E 3 ) =, 6 E E E 3 which coincides with the integrand of the right hand side of (8). In the diagram language, the symmetrization procedure can be written in the form: = + +. () 6 An analogous symmetrization for arbitrary diagrams can be written in the form of a symbolic equation: () Obviously, the following two equations are also valid (3) (4) Using these equalities, we consider a more complicated example of the sum of three diagrams J = t0 t t ω=0 + t0 t t ω=0 + t0 t t ω=0 (5) Calculating this sum using time versions and performing symmetrization, we obtain: J = ( 0 0 0 ) + 0 0 (6) The first diagram in (5) has one time version and is divided into half the sum of the first two diagrams in (6), the second diagram in (5) has two time versions corresponding to diagrams 3 and 4 in (6), and the last diagram has one time version corresponding
Diagram Reduction in Critical Dynamics 7 to the fifth contribution in (6). Using the relations ()-(4) in (6), we find J = 0 + 4 0 = 4 (7) To simplify the first brackets in (6) we used the equality (), for the second one (4) and for the last transition (). As a result, the sum of the dynamic diagrams (5) has been reduced to a single static one. This example shows how this technique reduces the number of diagrams. As for diagrams of the one-irreducible function Γ ψ ψ, a complete reduction to static diagrams is not possible. However, the use of the relations ()-(4) allows one even in this case to reduce significantly the number of contributions and to simplify their form. Let us consider the sum of the last two diagrams in (5), which we rewrite as: 0 J = 0 + (8) Taking into account that the values of the diagrams do not depend on the numbering of the vertices, making symmetrization and using the relation (4), we obtain: J = 0 + 0 + 0 + 0 = So one can see that even in this more complicated case reduction is possible. Now let us formulate a general recipe for the reduction of diagrams, illustrating it with the example of the sum of the following two diagrams, containing in aggregate ten time versions: The result of the reduction is the sum of the diagrams constructed according to the following recipe. (9) (30) Step I. Draw the diagrams of the static theory so that vertices with external legs are extreme left and right, while other (internal) vertices ordered in all possible ways so that nearest vertices are connected to each other. Possible order of vertices: Forbidden one: (3) (3)
Diagram Reduction in Critical Dynamics 8 Step II. On basis of the diagrams from step I, draw a set of diagrams with dashed lines (from one to (number of vertices -) sections) starting from the left vertex: Step III. In all the lines coming from the left to the vertex located between the two sections, we arrange the strokes in all possible ways. If there are similar lines, on which it is possible to arrange strokes, then we put only one, the remaining variant is taken into account by the symmetry coefficient. Thus, the original sum of ten time versions has been reduced to the sum of three effective diagrams (34). (See another example in the Appendix A) (33) (34) 5. Results The result of the reduction of the diagrams of Γ ψ ψ up to the four-loop approximation is depicted on Figs. -. The diagrams were calculated in the Feynman representation (a) (b) Figure : Diagrams of Γ ψ ψ after reduction: (a) two loops and (b) three loops using the Sector Decomposition method []. The required number of terms of their ε- expansion is given in the table, which also shows the corresponding symmetric factors S and additional weight factors f(n), which allow one to turn from the results for the one-component field with n = to the results for an n-component O(n)-symmetric model. The coefficients A (m) in the expansion (0) are determined from the data of the table by the relation A (m) = j S (m) j f (m) j (n)d (m) j (ε). (35) Here k = n + 8 9, k = n + 6n + 0 7, k 3 = n + 3, k 4 = 5n + 7. (36)
Diagram Reduction in Critical Dynamics 9 A A A3 B B B3 C C C3 C4 D D D3 D4 D5 D6 D7 Figure : Four loop diagrams of Γ ψ ψ after reduction grouped by graph topology Renormalization constant Z g in (0) is known from the statics and with the required accuracy is equal to Z g = + u 8 + n ( ) (8 + n) + u 4 + 3n + O(u 3 ). (37) 6ε 36ε 4ε Regarding the value of Z τ in (0), we need to make the following remark. The use of the static renormalization constants Z g and Z τ in (0) implies that in calculating of the function Γ ψ ψ one takes into account all diagrams, including ones that contain tadpole subgraphs. It is known that such diagrams can be ignored (which we did), if we do not take into account the tadpoles in counterterms as well. So, if tadpoles are not taken into account while calculating Γ ψ ψ, than to be consistent we should remove corresponding contributions from the renormalization constant Z τ. The resulting renormalization constant Z τ with the necessary accuracy is given by the expression ( ) + n Z τ = u 5( + n) + + O(u 3 ). (38) ε 44ε Note, that the problem of tadpoles is absent if the calculations of the renormalization constants are carried out in the massless theory with τ = 0, in which the tadpoles are
Diagram Reduction in Critical Dynamics 0 S (m) D (m) f (m) (n) m = ( loop) ε ε 0 ε /6 0.576556(0) 0.385354975(5) 0.56984660(35) k 3 m = 3 (3 loop) ε ε ε 0 /4 0.4384039(8) 0.33063368(30) 0.6007974(9) / 0.079054(4) 0.963588(6) 0.4485450(4) k k 3 m = 4 (4 loop) ε 3 ε ε A /8 0.0788085(9) 0.303746(4) 0.6479956(6) A /4 0.05394036(5) 0.749376(0) 0.45885(8) k k 3 A3 / 0.0697086() 0.0990954(0) 0.689090(4) B / 0.00489(4) 0.43086(8) B / 0.0376508(8) 0.0893734(4) k3 B3 /4 0.036365(4) 0.03958(8) C /4 0.05394039(4) 0.466740(0) 0.359704(0) C / 0.044906(8) C3 / 0.063848(3) 0.053546(7) k 3 k 4 C4 0.03485088() 0.04864(5) 0.779878() D /4 0.053940377(5) 0.46679(3) 0.886340(8) D / 0.05394044(5) 0.749404(0) 0.407878(8) D3 /4 0.0394088(8) D4 / 0.03485098() 0.048646(5) 0.4655896() k 3 k 4 D5 / 0.03485095(4) 0.048648(8) 0.39644(4) D6 / 0.007758853(0) 0.0305986(5) D7 /4 0.007758844(6) 0.0730478(0) Table : Values of ε-expansion of diagrams from Figs. - defined by zeros. In this theory, instead of (0), the value Γ ψ ψ ω=0,τ=0 is calculated for which the factor (µ /(τz τ )) ε/ on the right side of (0) is replaced by (µ/p) ε. However, with this approach, the integrands in the Feynman representation are slightly more complicated. Substituting the expressions (), (37), (38) in (0), and calculating values of A (m) with (35), we can find the coefficients z nm from the requirement of the cancellation of the pole contributions in ε. According to renormalization theory, the coefficients z nm at the highest poles in ε (m > ) are expressed in a certain way in terms of the coefficients z n of the first poles, which guarantees a cancellation of pole contributions in (0). This fact can be used as an additional self-consistency check for the multi-loop renormalization group calculations.
Diagram Reduction in Critical Dynamics The renormalization constant Z is associated with the RG-function γ γ (u) = β(u) u log Z. (39) The expression for the β-function is currently known with six-loop accuracy [9, 0]. We do not need its explicit form, since the connection mentioned above between the coefficients at the higher poles with the coefficient at the first pole makes it possible to represent γ (u) in a simpler form γ (u) = u u (z u + z 3 u 3 + z 4 u 4 +...). (40) The dynamic critical exponent z is expressed in terms of the value γ γ (u ) of the function γ (u) at the fixed point u and the Fisher exponent η by the relation [3] z = + γ η. (4) Substituting the values of z n into (40) and normalizing the result to the value of the two-loop contribution, we obtain: where γ = k 3 h u [ ] + b k u + (b k + b 3 k 3 + b 4 k 4 )u + O(u 5 4 ), (4) h = 6 ln(4/3).7609433, (43) b = 0.4939306(5), b = 0.5043(9), (44) b 3 = 0.69990(9), b 4 =.806593(30). (45) The value u, determined by the condition β(u ) = 0, with the required accuracy is given by the expression u = 6 8(3n + 4) ε + ε 3 ( + + 33n 3 + 0n + 760n + n + 8 (n + 8) 3 4(n + 8) ) 5 + 4544 96(n + 8)(5n + )ζ(3) ε 3 + O(ε 4 ). (46) The results of the dynamic exponent are usually presented in the form z = + Rη. (47) The value of η can be written in a form similar to (40), with the required accuracy u [ ] η = k 3 + a k u + (a k + a 3 k 3 + a 4 k 4 )u + O(u 5 4 ), (48) where a = 3 8, a = 5 64, a 3 = 5 3, a 4 = 45 3. (49) From (47) and (4) taking into account (40), (48), (46) we get [ ( R = (6 ln(4/3) ) + c ε + c + (c ) ] 3 + c 4 n) ε + O(ε 4 ), (50) (n + 8)
Diagram Reduction in Critical Dynamics where the coefficients c i are determined by the relations c = h 3 h (b a ), ( ) 4h c = 3(h ) 3 a (a b ) + (b a ), (5) 4h c 3 = 3(h ) ((b a ) 44(b a ) + 8(b 3 a 3 ) + (b 4 a 4 )), h c 4 = 3(h ) (9(b a ) 0(b a ) + 8(b 3 a 3 ) + 0(b 4 a 4 )). The first two terms of the ε-expansion (50) do not depend on the number of components of the field n. The first of them was calculated in the work [9], the second in the work [7], where the expression for b was obtained b = π /8 F (/4) ln(4/3) 3 4 + 3 8 F (x) = ln 4 8 which, according to (49), (5) and (5) corresponds to x ln 3 0.4939303, (5) ln t dt, (53) t c 0.8848346. (54) In the work [9] the value R was calculated in the leading order of the /n expansion for an arbitrary dimension d: ( ) R = 4 dγ (d/ )/Γ(d ) 4 d 8 / dx[x( x)]. (55) 0 d/ The first terms in the expansion of this quantity with respect to ε = 4 d have the form R = (6 ln(4/3) ) ( 0.8848347 ε 0.0999596 ε + O(ε 3 ) ). (56) Taking into account that the first two terms of the ε-expansion (50) do not depend on n, they coincide with the corresponding contributions to (56), which is confirmed by the results of [9], [7]. The expansion (56) also determines the coefficient c in the quadratic by ε contribution in (50): c = 0.0999596. (57) The values of the coefficients c i in (50) obtained in this paper are c = 0.884840(7), c = 0.09995(6), c 3 =.54(34), c 4 = 4.7847(8). (58) Results obtained are in full agreement with three loop calculations [7], as well as with /n-expansion [9], as for the four loop contribution, the coefficients c i were first calculated in [0] by a different method with much less accuracy: c = 0.884(9), c = 0.00(4), c 3 =.5(4), c 4 = 4.78(6), (59)
Diagram Reduction in Critical Dynamics 3 as one can see our results are in agreement with this calculations as well. Here we also give the resulting ε-expansion directly for the dynamic index z for n = : z = + 0.0344656ε + 0.003673(0)ε 3 0.005579(5)ε 4 + O(ε 5 ). (60) 6. Conclusion In this paper we performed four loop calculation of the critical exponent z in the framework of ε-expansion and renormalization group. To perform this calculation we developed a method of reduction of the diagrams in the models of critical dynamics which allows to significantly reduce a number of diagrams to be calculated. This method combined with Feynman representation and Sector Decomposition technique [] allows us to obtain high precision numbers for four-loop contribution to dynamic exponent. The necessity of high loop calculations for model A was pointed out in [] where Borel resummation of the results of the work [0] was performed. It was shown that results of resummation are very sensitive to particular realizations of the summation, which must be a consequence of the insufficient number of terms of the ε-expansion. As it was noted, the model A is the simplest model of critical dynamics, for more complicated models renormalization group calculations are limited at maximum by two loop order. The lack of perturbative information in this models does not allow to make solid theoretical predictions, moreover in some models (e.g. model E) it is not possible to confidently distinguish concurrent asymptotic regimes. The method discussed in this paper allows to significantly reduce the total calculation time for such problems and opens the possibility to extend this calculations to higher loops and more complicated models. For example, in model A at 5 loop level number of diagrams is reduced from 05 to 0 and with more simple integrands, which gives us a possibility to reach high accuracy of numerical calculations. While for more complicated theories like model E of critical dynamics [, 4, 6, ] (where RG calculations are limited only by loop order) our preliminary estimations show that this factor may be even greater than 5 and this gives us a hope that 3 and 4 loop calculations in this models can be feasible. Acknowledgements Research was carried out using computational resources provided by Resource Center Computer Center of SPb (http://www.cc.spbu.ru/en). Figures of Feynman graphs in this article were created with JaxoDraw [3]
Diagram Reduction in Critical Dynamics 4 Appendix A. Example of diagram reduction As a second example of the diagram reduction, let us consider the sum of 8 time versions of the diagrams of the Γ ψ ψ of the following type: (A.) Step I. Step II. (A.) Step III. (A.3) As a result, 8 versions of the diagrams (A.) were reduced to 8 diagrams (A.4). (A.4) Appendix B. Feynman presentation The dependence on the integration momenta in the diagrams after the integration over time has a structure that makes it possible to turn to the Feynman representation. This allows us to use the Sector Decomposition method [], as in problems of critical statics [3]. As an example we will consider the second diagram from (34): 7 4 6 (B.) 3 5 The numbers on the lines denote the integration momenta flowing from left to right. The integral that corresponds to diagram (B.) looks as follows: δ(k + k + k 7 ) δ(k 5 + k 6 + k 7 ) δ(k + k k 3 k 4 ) J = dk... dk 7 E k E k3 E k4 E k5 E k6 E k7 (E k + E k + E k7 ) (E k3 + E k4 + E k7 ). (B.)
Diagram Reduction in Critical Dynamics 5 Associating each of the 8 factors in the denominator of (B.) to the Feynman parameter v i and using the Feynman formula, we obtain: where F ({v}) = dk... J = 0 ( dvi δ vi ) F ({v}), dk 7 δ(k + k + k 7 ) δ(k 5 + k 6 + k 7 ) δ(k + k k 3 k 4 ) Q α Q = v E k + v 3 E k3 + v 4 E k4 + v 5 E k5 + v 6 E k6 + v 7 E k7 + (B.3), (B.4) + v 8 (E k + E k + E k7 ) + v 9 (E k3 + E k4 + E k7 ), (B.5) α = 8 the number of factors in the denominator of (B.). Writing Q in the form where Q = u E k + u E k + u 3 E k3 + u 4 E k4 + u 5 E k5 + u 6 E k6 + u 7 E k7, u = v + v 8, u = v 8, u 3 = v 3 + v 9, u 4 = v 4 + v 9, u 5 = v 5, u 6 = v 6, u 7 = v 7 + v 8 + v 9, (B.6) (B.7) from (B.4), (B.5) we obtain δ(k + k + k 7 ) δ(k 5 + k 6 + k 7 ) δ(k + k k 3 k 4 ) F ({v}) = dk...dk 7 ( 7 j= u. (B.8) je kj ) α Choosing in (B.8) as independent variables a certain set {k i, k i, k i3, k i4 }, and performing the integration with the help of δ-functions, we arrive at an expression of the form 7 F ({v}) = dk i dk i dk i3 dk i4 (C + V ij,i l k ij k il ), C τ u α j. (B.9) Calculating the integral of the power of the quadratic form, we obtain F ({v}) = π dl/ C 3d/ α Γ(α dl/) (det V ) d/, Γ(α) j= (B.0) where L is the number of loops in the diagram, in the case under consideration L = 4. The value of the determinant det V in (B.0) does not depend on the choice of the variables of integration {k i, k i, k i3, k i4 } and can be determined directly from the diagram view. By construction, det V is the sum of products of four factors u i. For any set of independent variables of integration {k i, k i, k i3, k i4 } the diagonal elements of the matrix V are equal to u i, u i, u i3, u i4, their product contributes to det V with coefficient one. The nondiagonal elements of matrix V do not contain the parameters u i, u i, u i3, u i4, consequently, det V does not contain the highest powers of u i. Obviously, we can not choose as independent variables some sets of {k i, k i, k i3, k i4 } which form conservation laws. As a result such a products of u i will not appear in det V.
Diagram Reduction in Critical Dynamics 6 For our particular diagram, in addition to the sets {k, k, k 7 }, {k 5, k 6, k 7 } and {k, k, k 3, k 4 }, which are defined by the δ-function arguments in (B.), conservation laws are also formed by sets {k 3, k 3, k 7 }, {k 3, k 4, k 5, k 6 } and {k, k, k 5, k 6 }. Thus, out of 35 possible quadruples of products of the parameters u i 5 products do not contribute to det V, and det V for the diagram (B.) is given by the expression det V = u u u 3 u 5 + u u u 3 u 6 + u u u 4 u 5 + u u u 4 u 6 + u u 3 u 4 u 5 + + u u 3 u 4 u 6 + u u 3 u 5 u 6 + u u 3 u 5 u 7 + u u 3 u 6 u 7 + u u 4 u 5 u 6 + + u u 4 u 5 u 7 + u u 4 u 6 u 7 + u u 3 u 4 u 5 + u u 3 u 4 u 6 + u u 3 u 5 u 6 + (B.) + u u 3 u 5 u 7 + u u 3 u 6 u 7 + u u 4 u 5 u 6 + u u 4 u 5 u 7 + u u 4 u 6 u 7 in which u i must be expressed in terms of v i according to (B.7). For the first and third diagrams of the formula (34), the expression (B.) is preserved, only the connections (B.7) of the variables u i and v i will change accordingly. This can be easily found by the form of the diagram. References [] Hohenberg P.C.,Halperin B.I. 977, Theory of dynamic critical phenomena Rev. Mod. Phys. 49 435 [] Folk R., Moser G. 006, Critical dynamics: a field-theoretical approach, J. Phys. A: Math. Gen. 39 R07 R33 [3] Vasil ev A.N. 004, The Field Theoretic Renormalization Group in Critical Behavior Theory and Stochastic Dynamics, Chapman and Hall, Chapman& Hall/CRC Boca Raton [4] Marinelli M., Mercuri F., Belanger D.P. 995, Specific heat, thermal diffusivity, and thermal conductivity of F ef at the Néel temperature, Phys. Rev. B 5 8897 [5] Niermann D., Grams C.P., Becker,Bohaty L., Schenck H., Hemberger J. 05, Critical Slowing Down near the Multiferroic Phase Transition in MnW O 4, Phys. Rev. Lett. 4 03704 [6] Livet F., Fevre M., Beutier G., Sutton M. 05, Ordering fluctuation dynamics in AuAgZn, Phys. Rev. B 9 0940 [7] Antonov N.V., Vasil ev A.N. 984, Critical dynamics as a field theory, Theor. Math. Phys. 60 67 679 [8] Batkovich D.V., Kompaniets M.V., Chetyrkin K.G. 06, Six loop analytical calculation of the field anomalous dimension and the critical exponent η in O(n)-symmetric ϕ 4 model, Nucl. Phys. B 906 47 67 (arxiv:60.0960 [hep-th]) [9] Kompaniets M.V., Panzer E. 06, Proc. of Loops and Legs in Quantum Field Theory (4-9 April 06 Leipzig Germany), PoS(LL06)038 (arxiv:606.090 [hep-th]) [0] Kompaniets M.V., Panzer E. 07, Minimally subtracted six loop renormalization of O(n)- symmetric ϕ 4 theory and critical exponents, Phys. Rev. D 96 03606 (arxiv:705.06483 [hepth]) [] Schnetz O. 06, Numbers and functions in quantum field theory, arxiv:606.08598 [hep-th] [] Binoth T., Heinrich G. 000, An automatized algorithm to compute infrared divergent multi-loop integrals, Nucl. Phys. B 585 74 759 [3] Adzhemyan L. Ts., Kompaniets M.V. 04, Five-loop numerical evaluation of critical exponents of the ϕ 4 theory, J. Phys.: Conf. Ser. 53 0049 [4] Adzhemyan L.Ts., Danco M., Hnatic M., Ivanova E.V., Kompaniets M.V. 06, Multi-Loop Calculations of Anomalous Exponents in the Models of Critical Dynamics, EPJ Web of Conferences 08 0004
Diagram Reduction in Critical Dynamics 7 [5] Siggia E.D., Halperin B.I., Hohenberg P.C. 976, Renormalization-group treatment of the critical dynamics of the binary-fluid and gas-liquid transitions, Phys. Rev. B 3 0 [6] De Dominicis C., Peliti L. 978, Field-theory renormalization and critical dynamics above T c : Helium, antiferromagnets, and liquid-gas systems, Phys. Rev. B 8 353 [7] Adzhemyan L. Ts., Vasiliev A. N., Kabrits Yu.S., Kompaniets M.V. 999, H-model of critical dynamics: Two-loop calculations of RG functions and critical indices, Theor. and Math. Phys. 9: 454 470 [8] Folk R., Moser G. 00, Dynamic Critical Behavior Near the Superfluid Transition in 3 He 4 He Mixtures in Two Loop Order, Phys. Rev. Lett. 89 530 [9] Halperin B.I., Hohenberg P.C., Ma S. 97, Calculation of dynamic critical properties using Wilson s expansion methods, Phys. Rev. Lett. 9 548 55 [0] Adzhemyan L. Ts., Novikov S. V. and Sladkoff L. 008, Calculation of dynamical exponent in model A of critical denamics to order ε 4, Vestnik SPbU Ser. 4 4 09 [in Russian] [] Nalimov M.Yu., Sergeev V. A., Sladkoff L. 009, Borel resummation of the ε-expansion of the dynamical exponent z in model a of the ϕ 4 (O(n)) theory, Theor. Math. Phys. 59: 499 508 [] Dohm V., 979, Density correlation function and dynamic transient exponents for liquid helium at and above T λ, Z. Physik B 33, 79-95 [3] D. Binosi and L. Theußl,004, JaxoDraw: A graphical user interface for drawing Feynman diagrams, Comput. Phys. Commun. 6, 76