Log expansions from combinatorial Dyson–Schwinger equations

We give a precise connection between combinatorial Dyson–Schwinger equations and log expansions for Green’s functions in quantum field theory. The latter are triangular power series in the coupling constant α\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\alpha $$\end{document} and a logarithmic energy scale L—a reordering of terms as G(α,L)=1±∑j≥0αjHj(αL)\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$G(\alpha ,L) = 1 \pm \sum _{j \ge 0} \alpha ^j H_j(\alpha L)$$\end{document} is the corresponding log expansion. In a first part of this paper, we derive the leading log order H0\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$H_0$$\end{document} and the next-to(j)\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$^{(j)}$$\end{document}-leading log orders Hj\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$H_j$$\end{document} from the Callan–Symanzik equation. In particular, Hj\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$H_j$$\end{document} only depends on the (j+1)\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$(j+1)$$\end{document}-loop β\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\beta $$\end{document}-function and anomalous dimensions. In two specific examples, our formulas reproduce the known expressions for the next-to-next-to-leading log approximation in the literature: for the photon propagator Green’s function in quantum electrodynamics and in a toy model, where all Feynman graphs with vertex sub-divergences are neglected. In a second part of this work, we review the connection between the Callan–Symanzik equation and Dyson–Schwinger equations, i.e., fixed-point relations for the Green’s functions. Combining the arguments, our work provides a derivation of the log expansions for Green’s functions from the corresponding Dyson–Schwinger equations.


Prologue: Green's functions in quantum field theory
Quantum field theory (QFT) predicts probabilities for certain particle processes. For example, initialize an experiment, where a virtual photon decays into an electronpositron pair. One cannot predict whether and when the photon decays, and if it does, where the electron and positron will go in the end. However, one can put two detectors D 1 and D 2 somewhere and predict the probability that an electron enters D 1 and a positron enters D 2 . This situation is totally different from deterministic classical mechanics, which predicts the exact time evolution for given initial conditions.
Such probabilities can be found in the following way: Each QFT is defined by a Lagrangian, which is a functional of the described particle fields. It consists of a kinetic part (terms that are quadratic in the fields) and interaction terms (of cubic or higher order in the fields). From all these, one can read off Feynman rules in a simple way. Now, a given particle process corresponds to a diagram (called Feynman graph) that translates via these rules to the so-called Feynman amplitude. Finally, the latter translates via the LSZ formula to a probability amplitude that squares to the actual probability. 1 Consider, for example, the above experiment. The QFT that describes this process is quantum electrodynamics (QED). The corresponding Lagrangian is a functional of the photon field A μ , the electron ψ and its antiparticleψ (i.e., the positron). It contains one kinetic term for the photon and one for the electron and positron together. It also includes one interaction term, i.e., gψ A μ γ μ ψ. 2 Here, the coupling constant g 1 is the electric charge of the electron and γ μ is a Dirac matrix. Each process involving these particles translates to a graph that consists of edges and vertices: A photon corresponds to a wiggly line and a straight arrow line indicates an electron or positron (depending on the direction of the arrow). An interaction between these three particles is represented by a vertex. Table 1 shows all these basic ingredients of QED Feynman graphs. The virtual photon decay described above corresponds to the vertex graph in  Table 1, and Feynman rules state that the respective Feynman amplitude is proportional to the electric charge g. It also depends on other factors, e.g., on the relative positions of the detectors.
Quantum mechanics now tells us that the obtained result is not exact. Indeed, one has to consider all possible ways in which the final state (an electron in D 1 and a positron in D 2 ) is achieved. For example, both particles could actually interchange another photon before entering the detectors. The respective Feynman diagram is the 1-loop graph , (1.1) and one has to add the resulting Feynman amplitude to the previous one before computing the actual probability. Note that the latter amplitude is of order g 3 , because by Feynman rules, each vertex constitutes one factor of g. Thus, the contribution to the overall probability is small compared to the initial one. All in all, one has to consider infinitely many Feynman graphs with more and more loops and vertices. However, their contributions to the final amplitude (the quantum corrections) become smaller and smaller. At some point, one may truncate the resulting Feynman amplitude when it is accurate enough.
There are two problems in this calculation: First, each loop diagram translates to a divergent integral via Feynman rules. In order to extract the correct quantum corrections, one has to regularize the integrals, i.e., to keep track of the different kinds of divergences. Then, one applies a renormalization scheme, i.e., one introduces 'counter terms' into the Lagrangian. These lead to additional Feynman graphs, whose amplitudes exactly cancel the divergences. In this way, the resulting Feynman amplitudes become finite. 3 A second problem are infrared divergences that occur when the probability of a particle process is computed from the Feynman amplitude. We do not discuss these any further, because our work is focused on the computation of Feynman amplitudes.
After regularization and renormalization in a momentum subtraction scheme, each Feynman graph contributes a quantum correction φ R ( ) to the initial Feynman amplitude. This correction depends on the coupling constants g k and the scalar products p i · p j between the momenta p μ i of the particles that enter the process. Here, it is very convenient to factor out an energy scale S and to define dimensionless scattering angles as = {p i · p j /S}. Now, renormalization requires some sort of boundary condition. In our case, we assume that φ R ( ) is known for certain values of S = S 0 and = 0 . One usually calls S 0 renormalization scale and {S 0 , 0 } renormalization point. 4 Then, it turns out that φ R ( ) is a polynomial in the coupling constants g k and the logarithm L = log(S/S 0 ) of the energy scale [2]. Hence, renormalized Feynman rules can be written as a linear map φ R from the set of Feynman graphs to the polynomial algebra A[g k , L]. The dependence of φ R ( ) on the scattering angles and 0 is hidden in the coefficients of the polynomial. In the following, we will make this more explicit.
First, φ R is an algebra homomorphism-the domain of definition is the algebra H of all one-particle irreducible (1PI) Feynman loop graphs 5 (of a given QFT) and their disjoint unions. The disjoint union is an associative product, and the empty graph I is the corresponding unit element. Furthermore, Feynman rules are normalized such that a 1PI Feynman graph contributes a change of φ R ( ) relative to the corresponding zero-loop Feynman amplitude. 6 This implies that 1PI graphs without loops would map to 1 (the unit element in A). Note that these are not in H, but are identified with the empty graph I. For example, the QED graphs in Table 1 would map to 1 by renormalized Feynman rules 7 and the 1-loop graph in Eq. (1.1) maps to a term that is proportional to g 2 (and not ∝ g 3 ). With this normalization, a one-particle reducible graph (1PR) corresponds to the disjoint union of its 1PI parts, e.g., (1.2) Secondly, the algebra H acquires a grading. In general, it is graded by the number of vertices of each kind in a graph. However, for QFTs with only one interaction term, it is more convenient to grade the algebra H by the number of loops, For example, H 0 = {I}, the graph in Eq. (1.1) belongs to H 1 and the disjoint union in Eq. (1.2) is an element in H 3 . In this case, there exists a redefinition of the coupling constant as α = g l , such that renormalized Feynman rules map each k-loop graph to a term proportional to α k . We show this explicitly in Sect. 3.1. For a vertex with three outgoing edges (e.g., QED), one finds l = 2, whereas for a vertex with four outgoing edges (e.g., in φ 4 theory), l = 1. For simplicity, we state our results using this redefinition. The generalization to QFTs with more than one coupling constant is presented in Sect. 4. All in all, renormalized Feynman rules are an algebra homomorphism and the quantum corrections of a certain particle process are found by applying φ R to an infinite sum of elements in H (1PI graphs and disjoint unions thereof). This results in a series expansion in the coupling constant α and the logarithm L, where the dependence on the scattering angles is hidden in the coefficients. Moreover, the normalization of Feynman rules implies that φ R evaluates Feynman graphs for scattering angles at the renormalization point, → 0 (see [3]). We will shortly explain the notion of Green's functions, using the above definition of renormalized Feynman rules. However, we first need some more notations: The residue r of a Feynman graph is its external leg structure. For example, the graph in Eq. (1.1) has residue r = and those in Eq. (1.2) have residue r = . Furthermore, its degree |r | denotes the number of external edges-Feynman diagrams with |r | = 2 are called 'propagator graphs' and those with |r | > 2 'vertex graphs'. Finally, we define sgn(r ), which is +1 for vertex-and −1 for propagator graphs. Now, in order to compute the exact Feynman amplitude of a certain particle process, one has to find the Green's functions of the respective QFT. These encode all necessary information about the full quantum corrections of all particle processes under consideration. For example, consider another QED experiment: a source S emits an electron, which is absorbed at a later time by a detector D. The simplest way for the electron would be to go straight from S to D. Then, the quantum corrections to the corresponding Feynman amplitude are obtained by applying renormalized Feynman rules to the sum of all Feynman diagrams with residue r = . We denote this sum by an ellipse. It is given by where the circle denotes the sum of all 1PI loop graphs. Thus, thanks to the geometric series, This works for all propagator graphs. Hence, let X r = I − 1PI loop graphs with residue r , |r | = 2. (1.5) Then, the corresponding propagator Green's function G r is defined as G r (α, L) = φ R (X r ). It is the inverse of the quantum corrections to the respective 1-particle process.
(The particle goes from S to D.) For interaction processes, the situation is similar, but there is no geometric series to compute. One simply defines X r = I + 1PI loop graphs with residue r , |r | > 2, (1. 6) and calls G r with G r (α, L) = φ R (X r ) the vertex Green's function of residue r . Knowing these functions allows to compute the correct Feynman amplitude of any interaction process that is described by the theory. Let us illustrate this again at the example of QED. First, consider the photon decay described at the beginning. The full quantum corrections are given by the Green's functions as follows: (1.7) Here, the circle on the vertex denotes the sum of all 1PI graphs including the zero-loop graph, in contrast to the 'propagator circles'. Indeed, the argument of φ R on the above lhs denotes the sum of all Feynman graphs with residue r = . Secondly, consider a scattering between two electrons. The corresponding Feynman graphs have four external legs, and in this case, the zero-loop graph is not 1PI. But again, the full quantum corrections to the corresponding Feynman amplitude are given by the Green's functions: (1.8) Note that the argument of φ R is the sum of all Feynman graphs with residue . As before, the circle on the vertex with four outgoing edges denotes the sum of all corre-sponding 1PI graphs. However, it does not contain the respective zero-loop Feynman graph, because it does not exist. 8 Let us summarize: One can compute the Feynman amplitude of any particle process when the Green's functions of the corresponding QFT are known. These are a power series in the coupling constant α and the logarithm of the energy scale L = log S/S 0 . Hence, G r ∈ A[α, L]. But renormalized Feynman rules tell us slightly more, namely that a k-loop graph maps to a polynomial of degree smaller or equal than k in L. We conclude that G r can be written as a triangular power series, where G r k is a degree k polynomial that encodes the contributions of all k-loop graphs to G r . 9 In particular, these polynomials vanish at the renormalization scale (for L = log S/S 0 = 0), such that The Green's functions are infinite power series, and the best one can do is to approximate them perturbatively. One possibility is to truncate the sum in Eq. (1.9) at some loop order. Since α 1, this is a reasonable approach, at least for energies near the renormalization scale S ∼ S 0 . However, the perturbative computation breaks down at energies far away from S 0 , where L 1. Then, it is better to reorder the terms in Eq. (1.9) as follows: (1. 10) This is called the log expansion for the Green's function G r . H r 0 is the leading log order (LLO), H r 1 is the next-to-leading log order (NLLO) and in general, H r j is the next-to ( j) -leading log order (N ( j) LLO). A truncation of the above sum, such that terms of order O(α n+1 ) are neglected, is called the next-to (n) -leading log (N (n) L L) approximation of H r . Perturbatively, this gives accurate results, as long as the functions H r j are regular at z = αL, which may well be the case for energies even far away from the renormalization scale.

Summary of the results
In this section, we outline our work and summarize the findings. For the explicit derivations, see the next sections.
In the first part of the paper, we derive the N ( j) LL approximation of the Green's functions H r in Eq. (1.10) with φ R (r ) = 1. These are the propagator as well as the vertex Green's function that belongs to the interaction term in the Lagrangian. The general result is given in Eq. (2.18) in its most compact form. It is obtained by solving the Callan-Symanzik equation, which is a simple first-order partial differential equation for G r describing its dependence on the energy scale L; see Eq. (2.7). It turns out that H r j only depends on the ( j + 1)-loop β-function and anomalous dimension. The latter are defined in Sect. 2.1 and only depend on the Feynman amplitudes of at most ( j + 1)-loop graphs.
In particle physics, it is long known how to solve the Callan-Symanzik equation in order to find the corresponding Green's function. See, for example, a standard textbook in quantum field theory (Peskin and Schroeder [4] in Section 12.3). There, the equations (12.76) and (12.79) show the two-and four-point Green's functions in a scalar φ 4 theory. These formulas are specific examples of our Eq. (2.18), but written in a way that is more adapted to particle physics.
We find the explicit log expansions of the Green's functions by integrating and expanding the exponent in Eq. (2.18). This was not done before in the literature 10 for the following reason: In QFT toy models and in many examples such as QED, socalled Ward-Takahashi identities (see, e.g., [5][6][7]) relate the Green's functions to the running coupling parameter. The latter can be found by solving the renormalization group equation (RGE), which is a linear first-order differential equation similar to the Callan-Symanzik equation. 11 Now, the standard procedure in particle physics is to seek the log expansion of the running coupling parameter and to deduce the one for the Green's functions from there. In contrast, our approach is more general and can be directly used, even without Ward-Takahashi identities.
In order to check our formulas, we apply them to a class of QFT toy models, in which Feynman graphs with vertex sub-divergences are neglected; see Sect. 2.3. An example is the photon propagator Green's function H in QED, thanks to the Ward-Takahashi identities. Here, we explicitly compute the NNLL approximations for the propagator Green's functions in Eq. (2.22) and give the precise connection to the log expansion of the running coupling parameter, see Eq. (2.21).
In the mathematical literature, there are two different direct approaches to obtain the N ( j) LL expansions in these toy models. This gives us the opportunity to compare our result to those derivations. In [8], the infinite sum X r in Eq. (1.5) is mapped to an element in the shuffle algebra of words. This paper describes an algorithm to filtrate a word into shuffles of letters and (anti-)symmetrized combinations of letters. The quasi-shuffle map θ operates as a symmetrizer, whereas the commutator [,] takes care of the anti-symmetrization. In the end, each term in the filtration belongs to a certain N ( j) LL order of the corresponding Green's function G r . Our results do not agree with the general findings of that work. The reason is a serious mistake in the filtration algorithm in [8]. More explicitly, the (anti-)symmetrization operators θ and [,] have been implemented self-consistently. However, combinations thereof lead to the wrong filtration. The error can be traced back to the computation of a shuffle product of a letter a with a commutator of two letters b and c. In particular, the resulting expression depends on the order in which the steps are performed: replacing the commutator by [b, c] = bc − cb before computing the shuffle product yields another result than doing it the other way around. 12 Another approach to derive the log expansions in such toy models is given in [10] using chord diagrams. Our expression in Eq. (2.22) agrees with [10], and we give an explicit translation of the notation in Eq. (2.23).
In the second part of this paper, we review the connection between the Callan-Symanzik equation and Dyson-Schwinger equations (DSEs) [11,12]. The latter are fixed-point relations for the infinite sums of Feynman graphs in Eqs. (1.5) and (1.6). The connection is established by extending H to the Hopf algebra of Feynman graphs [13,14] and writing the DSEs using so-called insertion operators [15]; see Eq. The paper is organized as follows: In the next section, we derive the log expansions of the Green's functions H r with φ R (r ) = 1 for QFTs with only one interaction term in the Lagrangian. We also give the example of the propagator Green's functions in toy models that neglect Feynman graphs with vertex sub-divergences. The results are obtained from the Callan-Symanzik equation (2.7). In Sect. 3, we rederive that relation from Dyson-Schwinger equations for the infinite sums of Feynman graphs in Eqs. (1.5) and (1.6). There, we also review some known aspects about the combinatorics of Feynman graphs (we extend H to a Hopf algebra and rederive the essential properties of renormalized Feynman rules φ R ). In Sect. 4, we generalize our results to QFTs with more than one interaction term in the Lagrangian and conclude in Sect. 5.

Results
In this section, we present the derivation of the log expansion H r in Eq. (1.10) with φ R (r ) = 1 from the Callan-Symanzik equation (2.7); see Sect. 2.2. The result contains the anomalous dimension and the β-function, which we introduce in Sect. 2.1. Finally, Sect. 2.3 gives a simple example of the log expansion for propagator Green's functions in QFT toy models, where all Feynman graphs with vertex sub-divergences are neglected.

The anomalous dimensions andˇ-function
For each vertex residue r = v, define Here, e ∼ v means that the edge e is incident to the vertex v. The square root must be expanded in a Taylor series, similar to the geometric series of propagator Feynman graphs. For example, in QED: The square roots are graphically represented by half ellipses. Since renormalized Feynman rules φ R are an algebra homomorphism, one can easily evaluate these infinite sums of Feynman graphs: In physics literature, G v amp is often referred to as the amputated Green's function for the vertex v. In more mathematical literature, these functions are called the invariant charges.
Let us shortly motivate the above definition by showing how the amputated Green's functions appear in physics computations. For example, consider any particle interaction process that is described by the QFT. Let v be the residue of the respective Feynman graphs. Then, one can easily show that the corresponding Feynman amplitude (including all quantum corrections) is given by e∼v (G e ) −1/2 times a polynomial of the amputated Green's functions. For example, the right-hand sides of Eqs. (1.7) and (1.8) read Secondly, it is common in particle physics to redefine the coupling constants of the theory asg k = g k G v k amp . This makes them dependent on the energy scale and has the advantage that a perturbative expansion for the Green's functions in terms ofg k and the logarithmic energy scale L may behave much better. This is an alternative way to approximate the Green's functions accurately without using the log expansion. Thẽ g k are called the running coupling parameters. In our case, there is only one coupling constant and we redefined it such that α = g l counts the number of loops in a Feynman graph. Hence, we must also redefine the respective invariant charge as 13 Then, the running coupling parameter isα = α Q. At the renormalization scale (L = 0), one hasα = α, because Q(α, 0) = 1. Now, the anomalous dimension γ r is defined as the negative of the L-linear part of the Green's function G r , Similarly, the β-function 14 is the L-linear part of the invariant charge, The latter identity follows from Eq. (2.2): the β-function is a linear combination of the various anomalous dimensions. For example, l = 2 in QED: Two remarks are in order: First, the notion of an anomalous dimension γ r for vertex-type residues is unusual in particle physics. Indeed, it is redundant, since γ v can be expressed by the β-function and the various propagator anomalous dimensions γ e via Eq. (2.5). In our case, however, the definition in Eq. (2.4) including vertextype residues is very convenient because the formulas for the log expansions of the Green's functions H r generalize to all kinds of residues. Secondly, the minus sign in Eq. (2.4) is chosen such that it cancels the sign in Eq. (1.5) for propagator-type residues r = e. In this case, the anomalous dimension γ e is given by the L-linear part of the Feynman amplitude that corresponds to the sum of all respective 1PI propagator graphs.
The functions γ r and β can be computed perturbatively. Neglecting terms of order O(α n+1 ) in Eqs. (2.4) and (2.5) gives the so-called n-loop anomalous dimensions and n-loop β-function. These approximations can be obtained by only evaluating Feynman graphs with at most n loops. It turns out that the N (n) LL approximation of the Green's functions H r only depend on the (n + 1)-loop anomalous dimensions γ r and β-function. We show this explicitly in the next section.

Log expansions
We start with the following formula for the Green's functions: A detailed derivation from combinatorial Dyson-Schwinger equations is collected in Sect. 3 (for a short proof, see Sect. 3.4). Taking a linear combination of the above rela-tions for different residues r implies an ordinary differential equation for the invariant charge, i.e., ∂ log Q ∂ L = β(α Q). (2.8) Together with the initial conditions, Equations (2.7) and (2.8) determine the full L-dependence of the Green's functions G r 15 . In particle physics, the above relations are usually written in terms of the (Ldependent) running coupling parameterα = α Q. Then, Eq. (2.8) corresponds to the renormalization group equation (RGE) 16 and Eq. (2.7) represents the Callan-Symanzik equation for the Green's functionG r withG r (α, L) = G r (α, L): The derivation of the log expansions for the Green's functions requires another differential equation for the invariant charge Q. The latter can be found as follows: Integrating Eq. (2.8) via separation of variables and using the above definition of the running coupling parameter leads to Here, the integration constant is well chosen such thatα = α for L = 0. Since we cannot solve this for Q, it is more convenient to take the derivative with respect to α. This results in (2.10) With these considerations, we can now give the main equations that determine the log expansions for the Green's functions. Here, it is convenient to define a log expansion for the invariant charge as well, 17 (2.11) 15 Note that this only requires to know the anomalous dimensions γ r and β-function of the QFT, i.e., to say the L-linear parts of the Green's functions. 16 The usual definition of the β-function in the literature slightly differs from ours: For example, the QED 1-loop β-function is β lit (α) = β 0α 2 . The notations in this paper are consistent, such that our 1-loop β-function is linear in the (loop counting) coupling parameter. 17 Accordingly, the log expansion for the running coupling is given byα(α, z) = α R(α, z).
A change of variables in Eqs. (2.7-2.10) finally results in (2.14) Now, one uses Eqs. (2.13) and (2.14) in order to find the N (n) LL approximation of the invariant charge. 18 The obtained result can then be used in Eq. (2.12) to deduce the N (n) LL approximation of the Green's functions H r . The first step is simply achieved by a Taylor series expansion of Eq. (2.14) in α. The 0th order reads which gives the LLO of the invariant charge. Now, consider all terms of order α n in Eq. (2.14) with n ≥ 1 and note that they do not contain any of the functions R N with N > n. Hence, collecting all terms that involve R n , one finds where the rhs only consists of the functions R 0 , . . . , R n−1 . This is a nice recursive formula for the various log orders of the invariant charge, except for n = 1. Two remarks are in order: First, the α-linear terms in Eq. (2.14) vanish completely, because an α-derivative was taken in order to obtain Eq. (2.10). However, R 1 can be found from the α-linear part of Eq. (2.13), which reads The solution of this ordinary first-order differential equation with the respective initial condition R 1 (0) = 0 is given by Secondly, an explicit form of the recursive relation for R n is not required. This is because Eq. (2.14) with cleared fractions can be given directly to a computer algebra program in order to obtain the various log orders. For example, the α 2 order in Eq. (2.14) reads 2188 O. Krüger In order to obtain the N (n) LL approximation of H r , one needs to expand the exponent up to order α n+1 . As promised earlier, this requires knowledge of the (n + 1)-loop anomalous dimension γ r and β-function. For example, the LL approximation requires to compute the nominator in the exponent up to first order in α. Using the one-loop anomalous dimension (see Eq. 2.4) and the explicit form of R 0 in Eq. (2.15), we find This expression agrees with the result already obtained in [9] using the shuffle algebra of words. We close this section by contentedly giving the NNLL approximation. Here, the exponent in Eq. (2.18) must be computed up to third order in α. This requires knowledge of the three-loop anomalous dimension and the invariant charge up to NNLL order (R 0 , R 1 and R 2 ). The integration is technical, but can always be done with the help of a computer algebra program. We used FORM [21] to obtain (2.20) 19 In [20], the author discussed the behavior of the Green's functions at fixed points. In this case, the invariant charge evaluates to unity, R(α, z) = 1 and Eq. (2.18) reproduces the known results.

Examples: truncated propagator Green's functions and QED
Let us give a short example of the above formula (2.20). In QED, the Ward-Takahashi identities state that the Green's functions for the fermion propagator G and the 3point vertex G coincide [22]. Therefore, the invariant charge in Eq. (2.3) is given by the inverse of the photon propagator Green's function. Furthermore, the β-function equals the respective anomalous dimension, see Eq. (2.6): The above relation is a special case of a broader class of QFT toy models that describe only one particle (edge-type e in Feynman graphs). In such models, the invariant charge and β-function are given by

Mathematical background
This section collects the essential properties of Feynman graphs and renormalized Feynman rules that lead to the Callan-Symanzik equation (2.7). The derivation of this relation is given in Sect. 3.4 and requires two main ingredients: The first one is the well-known exponential formula for φ R ; see Eq. (3.8). We motivate this property of renormalized Feynman rules and give a short introduction to the Hopf algebra of Feynman graphs in Sect. 3.2. The second ingredient is an (also well-known) expression for the coproduct of the infinite sums X r in Eqs. (1.5) and (1.6); see Eq. (3.10). We review in Sect. 3.3, how this relation is obtained from fixed-point equations for X r , i.e., so-called combinatorial Dyson-Schwinger equations (DSEs). In the literature, the latter are written down in full detail only for those sums X r with φ R (r ) = 1, 21 which is the reason why we restrict our results to that case. For presentation purposes, we consider QFTs with only one interaction term in the Lagrangian. (We explain the redefinition of the coupling constant in Sect. 3.1.) However, DSEs also exist in theories with more than one coupling parameter-the generalization to that case is given in Sect. 4.

The grading of H
For QFTs with only one interaction term in the Lagrangian, one redefines the coupling constant as α = g l , which implies that a K -loop graph maps to a term proportional to α K under renormalized Feynman rules. Let us show this now. First, Euler's formula states that for a graph with K loops, V vertices and E internal edges, the following identity holds: Secondly, let |v| be the degree of the vertex v 22 that corresponds to the interaction term in the Lagrangian. For example, |v| = 3 in QED and |v| = 4 in φ 4 theory. Then, a Feynman graph with residue r has |r | external edges; hence, Together with Euler's formula and the definition of l := 2/(|v| − 2), one finally finds V = l K + 0, for |r | = 2 1, for |r | = |v| . (3.1) One concludes that with the correct normalization of renormalized Feynman rules, a graph with K loops and V vertices maps to a term that is proportional to g l K = α K . 21 Hence, |r | = 2 or r = v is the vertex that corresponds to the interaction term in the Lagrangian. 22 That is, the number of edges which are incident to v.

Feynman rules from a Hopf algebraic point of view
This section represents a minimal introduction to the Hopf algebra of Feynman graphs [13,14]. As a motivation, consider a typical observation during the renormalization of Feynman amplitudes: 23 (3.2) Here, and only within this section, we write renormalized Feynman rules φ L i R with an upper index L i . It indicates the energy scale L i = log S i /S 0 at which Feynman graphs are evaluated. The goal of this section is to generalize the above observation (to the exponential formula). This can be achieved in a very elegant way, once we introduced the Hopf algebra of Feynman graphs.
Let us start with the algebra H, i.e., the vector space H equipped with an associative product m : H ⊗ H → H (the disjoint union) and a unit element I (the empty graph). In a more mathematical notation, where id denotes the identity map on H. In a first step, a coproduct and a co-unitÎ are defined, such that H extends to a bi-algebra. The important tool for our purposes is the coproduct. It is a linear map : H → H ⊗ H, which decomposes a graph into its subgraphs in a certain way. The explicit definition of its action on 1PI Feynman graphs is given by where the sum extends over (disjoint unions of) divergent subgraphs γ of . Furthermore, /γ corresponds to , but with the subgraphs γ shrunk to a point. For example, On a disjoint union of Feynman graphs, is defined such that it is compatible with the product (a requirement for bi-algebras), Here, Sweedler's notation for the coproduct is used, e.g., X = X ⊗ X . The co-unitÎ : H → K maps an element X ∈ H to the coefficient of I in X . Indeed, it is easy to show that the coproduct is co-associative and thatÎ is the co-unit element for , co-associativity: Hence, (H, m, I, ,Î) forms a bi-algebra. Secondly, one defines an antipode to further extend H to a Hopf algebra. It is the linear map S : H → H, which is recursively given by All in all, (H, m, I, ,Î, S) forms a Hopf algebra, called the Hopf algebra of Feynman graphs.
With the above-defined structures, the vector space Hom(H, A) of algebra homomorphisms extends to a group G A H , which is called the convolution or character group [23]. For example, renormalized Feynman rules are a character, φ R ∈ G A H . The conjunction in G A H is the convolution product : let φ, ψ ∈ G A H . Then, is also an algebra homomorphism; hence, φ ψ ∈ G A H . The convolution is associative because the coproduct is co-associative and the product m A in A[α, L] (multiplication) is associative. The unit element in G A H is the linear map e = 1 •Î, and the inverse of a character φ is given by φ • S.
We can now state a fundamental property of renormalized Feynman rules: This follows from a more general decomposition relation [2,3], which reduces to Eq. (3.7), because φ L R evaluates Feynman graphs for scattering angles at the renormalization point, → 0 . Note that Eq. (3.2) is a special case of Eq. (3.7).
The above relation implies that the vector space of maps φ L R , equipped with the convolution product, forms a group. This is the so-called renormalization group and it is isomorphic to the additive group. Hence, there exists a linear map σ : H → A, such that This is the exponential formula (see, e.g., [2,3]). Here, we dropped the upper index in φ R again, in order to adapt the notation to the rest of this paper. The map σ is an infinitesimal character, and we can give a meaning to it: it maps a (disjoint union of) Feynman graph(s) to the L-linear part of the corresponding Feynman amplitude, In particular, the infinite sums X r in Eqs. (1.5) and (1.6) are mapped to the negative of the anomalous dimension (see Eq. 2.4), From now on, we do not need the co-unit and antipode of the Hopf algebra anymore. These structures were only required for a proper definition of the convolution group G A H . The proof of Eq. (2.7) only uses the coproduct and the exponential formula in Eq. (3.8).

Graph insertion and Dyson-Schwinger equations
In this section, we motivate the following formula for the coproduct 24 of the infinite sums X r in Eqs. (1.5) and (1.6): 25 Here, X Q ∈ H is the combinatorial invariant charge (φ R (X Q ) = Q), see Eq. (2.3). For example, in QED: Furthermore, | k denotes the projection onto H k ; hence, φ R (X | k ) ∝ α k for each X ∈ H. Let us start by introducing the notion of graph insertion. Therefore, consider, for example, the QED Feynman graphs and . (3.11) The first one is obtained, when the subgraph is inserted into . However, the second one does not contain any subgraph and cannot be constructed via graph insertion. Note that the coproduct of a 1PI Feynman graph without subgraphs is given by = ⊗ I + I ⊗ ; see Eq. (3.3). Hence, φ R ( ) ∝ L by the exponential formula (3.8) and all graphs with that property are called primitive graphs. On the other hand, each non-primitive 1PI Feynman graph can be obtained via graph insertions of primitive ones.
In the following, we make this more explicit. Therefore, let be a primitive graph. Each vertex of represents an insertion place for a 1PI vertex graph. Furthermore, each internal edge of type e in is an insertion place for 1PI propagator graphs with residue e. Note that more than one propagator graph can be inserted into an edge, whereas only one vertex graph can be inserted into a vertex of . Now, for each primitive , define an insertion operator as a linear map B + : H → H , such that a disjoint union of 1PI Feynman graphs maps to Here, the sum extends over all graphs˜ that can be obtained by insertion of 1 , . . . , n into the various insertion places of . The factor sym(˜ ) is assigned such that the insertion operator commutes with the coproduct in the following way: 26 The following examples should make Eq. (3.12) clear. For a more rigorous definition, see [15]: Note that by Eq. (3.3), these examples imply Eq. (3.13). Some remarks are in order: First, inserting zero-loop graphs into does not change it. This is consistent with our realization of zero-loop graphs as the unit element I ∈ H. For example, B + (I) = and B + ( 1 I) = B + ( 1 ). Secondly, inserting more than V vertex graphs into a primitive with V vertices does not give a well-defined object in H. For example, there is no element in H that corresponds to Finally, there are non-primitive graphs in H that are not in the image of some insertion operator. For example, contains a subgraph, but it is not in the image of B + .
However, the infinite sums X r in Eqs. (1.5) and (1.6) consist of such linear combi- 26 The B + are Hochschild 1-co-cycles.
nations of 1PI Feynman graphs that are in the image of some insertion operator B + . This will be essential in the following. The closed formula for the coproduct in Eq. (3.10) can be found from fixed-point relations (DSEs) for the infinite sums X r . Let us motivate this again at the QED example, where the sums of all 1PI Feynman graphs of a certain residue satisfy [32] = + . . . , = + · · · , = + + · · · (3.14) Here, we recall the notation of the introduction: A circle represents the sum of all respective 1PI graphs and an ellipse denotes the geometric series of propagator Feynman graphs, as in Eq. (1.4). 27 On the rhs, these infinite sums are inserted into the one-loop primitives , and .
Here, only one vertex of the propagator primitives is dressed with an infinite sum in order to avoid double counting of graphs. The dots in Eq. (3.14) represent all graph insertions into higher loop primitives. The QED example is somewhat special, because the propagator DSE can be formulated for the quenched limit. Then, the first two fixed-point relations in Eq. (3.14) are actually complete. Equation (3.14) translate to DSEs for the infinite sums X r in Eqs. (1.5) and (1.6). These can be written in a very elegant way using the insertion operators: Note that there is no double counting of Feynman graphs, although the structure differs from the one in Eq. (3.14). Indeed, the (negative) power of X (X e ) in the argument of an insertion operator B + corresponds to the number of vertices (type e edges) in . Hence, in contrast to Eq. (3.14), each insertion place is dressed by an infinite sum.
Here, the insertion operators take care that there is no double counting and the deeper reason for that is the commutation relation in Eq. (3.13) [15]. With these considerations, the above relations generalize in a simple way. Denote the sum of all insertion operators B + for K -loop primitives with residue r by Then, the general DSEs [11,12,29,32] for the infinite sums X r in Eqs. (1.5) and (1.6) with φ R (r ) = 1 read Before we proceed, we quickly show that the insertion places of each K -loop primitive graph are fully dressed by infinite sums in Eq. (3.17): 28 Consider a K -loop primitive graph . The argument of the corresponding insertion operator B + is X r X K Q , where by Eqs. (2.1) and (2.3), The power of X v in the argument of B + is l · K for |r | = 2 and l · K + 1, if r = v. This matches the number of vertices V in , given by Eq. (3.1). Secondly, the power of 1/X e in the argument of B + for some edge-type residue e corresponds to This equals the number of type e internal edges in , which completes the proof. Note that Eq. (3.15) coincide with the general form in Eq. (3.17), but only include the K = 1 terms.
To close this section, we show that the coproduct formula in Eq. (3.10) and the commutation relation in Eq. (3.13) are consistent with the DSE in Eq. (3.17). With slightly more effort, one can turn this into a proof of Eq. (3.10) by induction. For details, see, e.g., [24][25][26][27]. First, let M(X r ) = r (X r ) s r with integers s r be an arbitrary monomial in X r . Then, Eq. (3.10) implies that Now, we compute the coproduct of the rhs in Eq. (3.17). We commute with the insertion operators via Eq. (3.13) and use the above relation for the monomials: [rhs(3.17)] = X r ⊗ I + sgn(r ) Finally, we rearrange the sums as k≥1 k K =1 using an index shift, which results in Note that the terms on the rhs of the ⊗ sign correspond to X r | k (see Eq. 3.17). Hence, we are left with

Rederivation of the Callan-Symanzik equation
With the considerations of the previous sections, the proof of Eq. (2.7) is simple. It requires the exponential formula (3.8) and the coproduct relation (3.10). Together with the definitions of the convolution product in Eq. (3.6) and the σ function in Eq. (3.9), we find

Multiple coupling constants
In this section, we rederive the Callan-Symanzik equation from Dyson-Schwinger equations in QFTs with multiple interaction terms in the Lagrangian. We assume that each such term comes with a separate coupling constant and gives rise to a certain vertex in Feynman graphs. We denote the different vertices by v 1 , . . . , v m and collect the coupling parameters to a vector g = (g 1 , . . . , g m ). The positive integer m is given by the number of interaction terms. In the following, we summarize the main changes to the one coupling case. First of all, it is inappropriate to grade the Hopf algebra H of Feynman graphs with respect to the loop number in this setting. Instead, H is multi-graded as [14,18,30,33] Here, k = (k 1 , . . . , k m ) denotes a multi-index (a vector of integers) and the subspaces H k contain the following Feynman graphs: Note that this grading respects the product and coproduct of H. Secondly, the Green's functions G r depend on all coupling constants g and the logarithmic energy scale L. In particular, G r (g, L) = φ R (X r ) as before, where the X r are given by the infinite sums in Eqs. (1.5) and (1.6). Hence, the corresponding anomalous dimensions γ r are power series in all coupling parameters, Furthermore, there is one amputated Green's function (invariant charge) for each vertex v i ; see Eq. (2.1). These define the running coupling parameters as before. In order to simplify the notation, we write Accordingly, there is one β-function for each vertex v i , with power series expansion Finally, the Dyson-Schwinger equation for the infinite sums in Eqs. (1.5) and (1.6) is given by is the generalization of Eq. (3.16); see, e.g., [18,34]. Note that in Eq. (4.1), each insertion place of a primitive is fully dressed with an infinite sum and there is no double counting of graphs because the insertion operators B k;r + are Hochschild 1-cocycles, i.e., they satisfy the commutation relation in Eq. (3.13). In full analogy to the end of Sect. 3.3, the DSE in Eq. (4.1) implies the coproduct formula (see, e.g., [30,31]) where k is the projection onto H k . This is the generalization of Eq. (3.10).
The Callan-Symanzik equation can now be derived in the same way as in Sect. 3.4. Therefore, the coproduct formula in Eq. (3.10) is replaced by Eq. (4.2) and the action of σ on the infinite sums X r in Eq. (3.9) must be generalized to σ X r k = −γ r k g k 1 1 · · · g k m m .
Then, using the exponential formula in Eq. (3.8) and the definition of the convolution product in Eq. (3.6) results in which can be written in the same form as in the one coupling case; see Eq. (2.7): ∂ ∂ L log G r (g, L) = −γ r (g),g = (g 1 Q 1 , . . . , g m Q m ). In particular, taking certain linear combinations of these relations yields a system of ordinary differential equations, which is the analogue to Eq. (2.8), ∂ ∂ L log Q i (g, L) = β i (g),g = (g 1 Q 1 , . . . , g m Q m ).
Secondly, the log expansions for the Green's functions can be derived as shown in Sect. 2.2 via a change of variables. Here, the main difference to the previous case is that we have several coupling constants, but only one energy scale. Hence, one has to choose one coupling g i and changes L → z = g i L. In particular, let G r (g, L) = H r (g, g i L). Then, the solution to Eq. (4.3) is given by dz . (4.5) This formula only requires the log expansion for the running coupling parametersg, which can be found log order by log order from the RGE. Note that the resulting log expansion will depend on the choice of g i . We hope to investigate such cases further and give examples in future work.

Conclusions
In this paper, we established a precise connection between Dyson-Schwinger equations for the Green's functions in a given QFT and their corresponding log expansions. We mainly discussed QFTs with only one interaction term in the Lagrangian (Sect. 2.2), but also expanded the results to the more general case (Sect. 4).
The formulas for the log expansions are given in Eqs. (2.18) and (4.5)-the nextto (n) -leading log approximation is obtained by an analytical integration that requires the solution of the renormalization group equation up to the next-to (n) -leading log order. In particular, it only depends on the (n + 1)-loop β-function and anomalous dimension. An explicit expression for the next-to-next-to-leading log approximation is given in Eq. (2.20), when only one coupling constant is involved. Restricting to special cases, our formulas reproduce the known expressions in the literature. These cases include the photon propagator Green's function in quantum electrodynamics and a toy model, where all Feynman graphs with vertex sub-divergences are neglected. Compare Eq. (2.19) with [9] and Eqs. (2.22, 2.23) with [10].
Our results may be used for an accurate perturbative calculation of the Green's functions, for energy scales far away from the renormalization scale S S 0 . However, the log expansions do only converge, if the coefficients do well behave perturbatively. The latter must not be true. So far, we have no control about the dependence of the Green's functions on the scattering angles 0 at the renormalization point. In order to understand this better, one needs to generalize the exponential formula in Eq. (3.8) to include also the angle dependence. Together with the coproduct formulas in Eqs. (3.10) and (4.2), this would lead to a generalization of the Callan-Symanzik equation, from which the correct perturbative approach may be obtained. However, this is a highly non-trivial problem, as it requires to understand the behavior of the Green's functions orthogonal to the renormalization group flow.
Our work gives the derivation of the log expansions for the Green's functions from the corresponding Dyson-Schwinger equations. However, if the latter are not known, we cannot comment on the corresponding log expansions. For example, in QED, a detailed Dyson-Schwinger equation for the 4-point function G in terms of the insertion operators is not explicitly written down in the literature. From there, one could generalize the coproduct formula in Eq. (3.10). Together with Eq. (2.7), this would lead to the respective Callan-Symanzik equation and then to the corresponding log expansion.
Finally, we would like to investigate the case with multiple coupling parameters. Here, we need to solve the system of renormalization group equations in Eq. (4.4) analytically, such that Eq. (4.5) can be integrated log order by log order.