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 $\alpha$ and a logarithmic energy scale $L$ --- a reordering of terms as $G(\alpha,L) = 1 \pm \sum_{j \geq 0} \alpha^j H_j(\alpha L)$ is the corresponding log-expansion. In a first part of this paper, we derive the leading-log order $H_0$ and the next-to$^{(j)}$-leading log orders $H_j$ from the Callan-Symanzik equation. In particular, $H_j$ only depends on the $(j+1)$-loop $\beta$-function and anomalous dimensions. 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), our formulas reproduce the known expressions for the next-to-next-to-leading log approximation in the literature. 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 if 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 experiment above. 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 anti-particleψ (i.e. the positron). It contains one kinetic term for the photon and one for the electron and positron together; as well as 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 Electron and Positron

Photon
Interaction vertex 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, whose 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 via Feynman rules to a divergent integral. 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, which 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], where 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.
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 QFT's 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 etc. In this case, there exists a redefinition of the coupling constant as α = g l , such that renormalized Feynman rules map each kloop graph to a term proportional to α k . We show this explicitly in Section 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 QFT's with more than one coupling constant is presented in Section 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 graphs 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 Θ and Θ 0 is hidden in the coefficients. With this definition of renormalized Feynman rules, we will shortly explain the notion of Green's functions, for which we 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 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. Therefore, 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 zeroloop graph, in contrast to the 'propagator circles'. Indeed, the argument of φ R on the lhs above 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 corresponding 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. Therefore, 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: 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) LL) 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.
In this work, we present the following results: • First, we derive the N (j) LL approximation of the Green's functions H r in Eq. (1.10) with φ R (r) = 1 10 from the Callan-Symanzik equation. The latter is a simple first-order partial differential equation for G r that describes its dependence on the energy scale L, see Eq. (2.7). In particular, H r j only depends on the (j + 1)-loop β-function and anomalous dimension, which are defined in Section 2.1 and only depend on the Feynman amplitudes of at most (j +1)-loop graphs. In principle, this calculation must be long known in particle physics. However, this is barely exhibited in sufficient clarity.
• As an application, we compute the NNLL approximation for the propagator Green's functions in a class of QFT toy models, in which Feynman graphs with vertex sub-divergences are neglected. An example is the photon propagator Green's function H in QED, thanks to the Ward-Takahashi identities (see e.g. [3][4][5]). For these toy models, NNLL approximations have already been derived using chord diagrams [6] or the shuffle Hopf algebra of words [7]. Our results (see Eq. (2.21)) agree with [6], but differ from [7]. The reason for the mismatch is a mistake in [7], where the shuffle products of commutators were computed using a wrong algorithm.
• In the second part of this paper, we review the connection between the Callan-Symanzik equation and Dyson-Schwinger equations (DSE's) [8,9], i.e. fixedpoint relations for the infinite sums of Feynman graphs in Eqs. (1.5) and (1.6). This can be seen by extending H to the Hopf algebra of Feynman graphs [10,11] and writing the DSE's using so-called insertion operators [12], see Eq. (3.17) 11 .
• Finally, we comment on possible extensions of our results. For the Green's functions H r in Eq. (1.10) with φ R (r) = 0 (such as H in QED), DSE's are barely discussed in details in the literature. When these are found, our results can be generalized in a simple way to find the corresponding log-expansions.
On the other hand, there are DSE's for Green's functions in QFT's with more than one interaction term [11,[14][15][16] -we generalize our results to such cases in Section 4.
The paper is organized as follows: In the next Section, we derive the logexpansions of the Green's functions H r with φ R (r) = 1 for QFT's 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 Section 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 Section 4, we generalize our results to QFT's with more than one interaction term in the Lagrangian and conclude in Section 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 Section 2.2. Therefore, we first introduce the anomalous dimension and the β-function in Section 2.1, as they appear in the final results. Finally, Section 2.3 gives a simple example of the logexpansion for propagator Green's functions in QFT toy models, where all Feynman graphs with vertex sub-divergences are neglected.

The Anomalous Dimensions and β-Function
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 righthandsides 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. Theg 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 12 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 (in our case, there is only one) is the L-linear part of the invariant charge, Indeed, it follows from Eq. (2.2) that it is just 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 quite 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 vertex type 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. 12 Here, v is the vertex residue that corresponds to the interaction term in the Lagrangian.

Log-Expansions
We start with the following formula for the Green's functions: A detailed derivation from combinatorial Dyson-Schwinger equations is collected in Section 3 (for a short proof, see Section 3.4). Taking a linear combination of the above relations 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 The derivation of the log-expansions for the Green's functions requires another differential equation for the invariant charge Q. Therefore, one integrates Eq. (2.8) using separation of variables. With the above definition of the running coupling parameter, this 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) 13 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.
14 The usual definition of the β-function in the literature slightly differs from ours: E.g., 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.
With these considerations, we can now give the main equations that determine the log-expansions for the Green's functions. One therefore defines a log-expansion for the invariant charge 15 , and 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 16 . 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 (z) with N > n. Hence, collecting all terms including the function R n (z), 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, 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 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. The computation is quite technical but can always be done with the help of a computer algebra program. Here, we contentedly give the NNLL approximation, i.e. (2.19)

Examples: Truncated Propagator Green's Functions and QED
Let us give a short example of the above formula (2.19). In QED, the Ward-Takahashi identities state that the Green's functions for the fermion propagator G and the 3point vertex G coincide. 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), where the invariant charge and β-function are given by (2.21)

Mathematical Background
This section collects the essential properties of Feynman graphs and renormalized Feynman rules (presented in a mathematical framework) that lead to the Callan-Symanzik equation (2.7). The derivation of this relation is given in Section 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 Section 3.2. The second ingredient is an (also well known) expression for the co-product of the infinite sums X r in Eqs. (1.5) and (1.6), see Eq. (3.10). We review in Section 3.3, how this relation is obtained from fixed-point equations for X r , i.e. so-called combinatorial Dyson-Schwinger equations (DSE's). In the literature, the latter are written down in full detail only for those sums X r with φ R (r) = 1 19 , which is the reason why we restrict our results to that case. For presentation purposes, we consider QFT's 17 Note that this is consistent with Eqs. (2.13, 2.18 and 2.20). 18 In [6], equivalent expressions have been derived using Chord diagrams. These formulas contain the symbols a i,j , which are related to the anomalous dimension γ e as follows: = −(a 3,0 + a 1,1 a 2,0 + a 2 1,1 a 1,0 + 3a 1,2 a 2 1,0 + 3a 2,1 a 1,0 ) . 19 Hence, |r| = 2 or r = v is the vertex that corresponds to the interaction term in the Lagrangian.
with only one interaction term in the Lagrangian (we explain the redefinition of the coupling constant in Section 3.1). However, DSE's also exist in theories with more than one coupling parameter -the generalization to that case is given in Section 4.

The Grading of H
For QFT's 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 20 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 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 lK = α K .

Feynman Rules from a Hopf Algebraic Point of View
This section represents a minimal introduction to the Hopf algebra of Feynman graphs [10,11]. As a motivation, consider a typical observation during the renormalization of Feynman amplitudes: 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 co-product ∆ and a co-unitÎ are defined, such that H extends to a bi-algebra. The important tool for our purposes is the co-product. It is a linear map ∆ : H → H ⊗ H, which decomposes a graph into its sub-graphs 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 sub-graphs γ of Γ. Furthermore, Γ/γ corresponds to Γ but with the sub-graphs γ shrinked 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 co-product 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 co-product 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 charactergroup [18]. 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 an algebra homomorphism, hence φ ⋆ ψ ∈ G A H . The convolution is associative because the co-product 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.
The observation in Eq. (3.2) can now be generalized in a beautiful way: Renormalized Feynman rules satisfy the property 21 It follows 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 σ ∈ G A H , such that This is the exponential formula. Here, we dropped the upper index in φ R again, in order to adapt the notation to the rest of this paper. We can give a meaning to the character σ: 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 co-product and the exponential formula in Eq. (3.8).

Graph Insertion and Dyson-Schwinger Equations
In this section, we motivate the following formula for the co-product 22 of the infinite sums X r in Eqs. (1.5) and (1.6) 23 : 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 sub-graph is inserted into . However, the second one does not contain any sub-graph and cannot be constructed via graph insertion. Note that the co-product of a 1PI Feynman graph Γ without sub-graphs 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 co-product in the following way 24 : The following examples should make Eq. (3.12) clear -for a more rigorous definition, see [12]: 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 sub-graph, 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 combinations 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 co-product in Eq. (3.10) can be found from fixed-point relations (DSE's) 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 [25] = + . . . , = + . . . , (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) 25 . On the rhs, these infinite sums are inserted into the one-loop primitives , and , where only one vertex of the propagator primitives is dressed with an infinite sum in order to avoid double counting of graphs. The dots in Eqs. (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 Eqs. (3.14) are actually complete. Eqs. (3.14) translate to DSE's 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 Eqs. (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 Eqs. (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) [12].
With these considerations, the above relations generalize in a simple way. Therefore, denote the sum of all insertion operators B Γ + for K-loop primitives Γ with residue r by Then, the general DSE's [8,9,24,25] for the infinite sums X r in Eqs. (1.5) and (1.6) with φ R (r) = 1 read Indeed, the insertion places of each K-loop primitive graph are fully dressed by infinite sums 26 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 To close this section, we show that the co-product 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. [19][20][21][22]. First, let M(X r ) = r (X r ) sr with integers s r be an arbitrary monomial in X r . Then, Eq. (3.10) implies that Now, we compute the co-product of the rhs in Eq. (3.17). Therefore, we commute ∆ with the insertion operators via Eq. (3.13) and use the above relation for the monomials: 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 co-product relation (3.10). Together with the definitions of the convolution product in Eq. (3.6) and the σ function in Eq. (3.9), we find B Γ + for some edge-type residue e is equal to This equals the number of type-e internal edges in Γ, which completes the proof.
In this section, we rederive the Callan-Symanzik equation from Dyson-Schwinger equations in QFT's with multiple interaction terms in the Lagrangian. Therefore, 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 [11,16,26,27] H = k H k .
Here, k = (k 1 , . . . , k m ) denotes a multi-index (a vector of integers) and the sub-spaces H k contain the following Feynman graphs: Note that this grading respects the product and co-product 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  27 . 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-co-cycles, i.e. they satisfy the commutation relation in Eq. (3.13). In full analogy to the end of Section 3.3, the DSE in Eq. (4.1) implies the co-product formula 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 Section 3.4. Therefore, the co-product 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 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 ) . ∂ ∂g i + γ r (g) G r (g, L) = 0 .
Secondly, the log-expansions for the Green's functions can be derived as in Section 2.2 via a change of variables L → z = g i L (for one of the coupling parameters g i ). In particular, let G r (g, L) = H r (g, g i L). Then, the solution to Eq. (4.3) is given by This formula only requires the log-expansion for the running coupling parametersg, which can be found log order by log order from the RGE. We hope to investigate this case 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 QFT's with only one interaction term in the Lagrangian (Section 2.2), but also expanded the results to the more general case (Section 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.19), when only one coupling constant is involved. In particular, 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) our formulas lead to the known expressions in the literature (Compare Eq. (2.21) with [6]).
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 Θ and Θ 0 . 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 co-product 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. We plan to do that in future work -with our analysis, this would lead to the respective Callan-Symanzik equation and then, to the corresponding log-expansion.
Finally, we would like to make the case with multiple coupling parameters more explicit. Therefore, 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.