-
The principle of renormalization [1-6] in quantum field theory [7-10] is to introduce divergent counterterms in the interaction (perturbation) Hamiltonian (or Lagrangian). These terms naturally yield new terms in the S matrix [11]. They are the counterterms to eliminate the ultraviolet divergence of Feynman integrals for all Feynman diagrams to any given order [4].
Thus there are two concepts: 1. What are the counterterms in the Feynman integrals for a given Feynman diagram [12], and are they enough to give a convergent result [4]? 2. What are the counterterms in the perturbation Hamiltonian? Can they precisely give the counterterms for the Feynman diagrams to a given order? The first question was answered by the BPHZ renormalization scheme [1,4]. The second question can easily be answered for Feynman diagrams without a symmetry factor. However, for symmetric Feynman diagrams, the situation becomes more involved. To the best of our knowledge, the consistency of these two concepts would still need an explicit proof, since the non-trivial symmetry factor is an important issue in perturbation field theory [13-15].
In this paper, we study the symmetry group
GΓ of a Feynman diagramΓ . We find that there is a subgroupG(1)Γ associated with a reduced Feynman diagram˜Γ which keeps the union of the reduced part ofΓ:⋃τγτ invariant. Furthermore, the symmetry group of the reduced Feynman diagrams˜Γ is the quotient groupG(1)Γ/∏τGγτ . We then explicitly give the counterterms in the perturbation Hamiltonian. We further prove that they give the counterterms of Feynman integrals for Feynman diagrams with symmetry factors. We remark that the idea of the current paper arose when the authors were editing a textbook on quantum field theory [16], and this paper is a further discussion of the material regarding the BPHZ scheme in that textbook.The paper is organized as follows. In the next section, we introduce the symmetry group of a Feynman diagram and derive the symmetry factor appearing in the S matrix via Wick's theorem. In Section III, we review the theory of the BPHZ renormalization scheme, then give the symmetry group of
˜Γ in Section IV. In Section V, we prove that we can consistently introduce new interaction terms (details are given in Appendix A) corresponding to each renormalization part of the Feynman diagram of a field theory, which will naturally give the counterterms required in the BPHZ scheme. -
The S matrix of a perturbation field theory is
Sfi=∑n(−i)nn!⟨f|T{∫d4x1⋯d4xnHI(x1)⋯HI(xn)}|i⟩≡∑nS(n)fi .
(1) Here
HI(x) is the perturbation Hamiltonian density, and it is a polynomial of derivatives of the fields. By Wick's theorem, the quantity⟨f|T{∫d4x1⋯d4xnHI(x1)⋯HI(xn)}|i⟩
(2) is expressed as a sum of the product of field (and derivative) contractions for all possible combinations. Each combination of a pairing of the fields corresponds to a Feynman diagram, where
HI(x) is related to some vertices. For each Feynman diagram, one can calculate a Feynman integral in either coordinate space or momentum space [17].Since in a vertex of
HI(x) several fields may be equivalent, the corresponding vertex lines are equivalent. For definiteness, we may assign numbers to the fields, i.e. give a number to each vertex line. Changing the assignment will, in general, give a new term in the Wick's expansion. Its value is equal to the original one. In this way, for a vertex with a group of l equivalent lines inϕl theory, we get a factoriall! . In Eq. (2), changing the assignment ofx1⋯xn for a given contraction in Wick's expansion will give a new term in this expansion. After integrating overx1⋯xn , these terms give the same result. In this way, we get a factorialn! . This means we haven! equal terms, each of them corresponding to an assignment of verticesxi for a Feynman diagram.A given perturbation Hamiltonian
HI(x) may have several vertices. Each vertex is composed of a vertex point x and several vertex lines. Each of them represents a field operator. Some vertex lines of the same vertex may be equivalent (Fig. 1). Once we assign a coordinate number to a vertex point:x=xi , and assign numbers1,2,⋯ to each set of equivalent vertex lines (of the chosen vertex), we actually fix all of the field operators in a vertex inHI(x) in Eqs. (1) and (2). Thus when all these vertex lines are connected pair-wise (including those in "in" and "out" states), we get a certain term in Wick's expansion. In this way, we get a Feynman diagram [12] with all vertices given coordinate numbers and all equivalent vertex lines of each vertex given line numbers.Figure 1. Vertex with vertex point and vertex lines.
(a) Three vertex lines are equivalent.(b) The two vertex linesα andβ are equivalent, whileγ andδ are equivalent.Alternatively, for a blank Feynman diagram, we first assign a coordinate number (
x1,⋯,xm ) to each vertex point (p1,⋯,pm ). Say,p1=x1,p2=x2⋯,pm=xm . Next, for each vertexpa (andpb ), we assign numbers1,2,⋯ to each end (attached to each vertex) of those equivalent vertex linesLabσ . Then we get a "marked" Feynman diagram, which corresponds to one term in Wick's expansion. Notice that only for topologically different "marked" Feynman diagram, we get different terms in Wick's expansion. While for different assignment of a blank Feynman diagram but topologically equivalent, they correspond to the same term in Wick's expansion (Fig. 2).Figure 2. Marking a Feynman diagram
Γ .(a) A blank Feynman diagramΓ before marking.(b) ,(b′) ,(c) , and(c′) are marked Feynman diagrams of(a) .(b) and(b′) are topologically equivalent (lines1 ,2 , and3 ofx1 connect to lines1 ,2 , and3 ofx2 , respectively).(c) and(c′) are topologically equivalent (lines1 ,2 , and3 ofx1 connect to lines2 ,1 , and3 ofx2 , respectively), while(b) and(c) are not equivalent.The operation of marking a blank Feynman diagram can be related to a group G. Changing the assignment (marking) of a blank (unmarked) Feynman diagram can be expressed as a group element
g=gn×∏igi.
(3) Here
gn is an element of the permutation groupSn , which corresponds to the reassignment of the coordinates of vertices (the permutation group is a symmetry of the Feynman diagram only after the integration over all coordinates is performed; this is also true for vertex lines).gi is an element of the symmetry group of each vertex which corresponds to the reassignment of the equivalent vertex lines of vertex pointpi . The operations of reassignment which do not change the topological structure of a marked Feynman diagram form a setGΓ , and it satisfies the properties of a group.GΓ is the symmetry group of a blank Feynman diagramΓ .We have
GΓ⊂G . For elementgΓ∈GΓ ,g∈G , and the element of G:g′=gΓg gives a topologically equivalent marked Feynman diagram to that given by g. Thus the number of cosetm=|G|/|GΓ| is the actual number of different terms in Wick's expansion which give the same blank Feynman integral form (2).We use the following (Figs. 3 and 4) to further illustrate the action of g,
gΓg , andggΓ .Figure 3. The action of
g (Eq. (3)). The first factor of g acts on vertices, and the rest of factors act on attached vertex lines.(a) A blank Feynman diagram.(b) The vertex points (p1 ,p2 ,p3 , andp4) and vertex lines (p1:α1β1γ1 ,p2:α2β2γ2 ,p3:α3β3γ3δ3 ,p4:α4β4γ4δ4 ), where all vertex lines are equivalent for each vertex.(c) An assignment of(a) , which gives a marked Feynman diagram(d) .(e) Another assignment of(a ), which gives a marked Feynman diagram(f) . The group element g is an operation changing the assignment, and it acts on(d) giving(f) . The group element g acts on(d) , while the indices1 ,2 ,⋯ in the above expression of g refer to the indices of the original blank Feynman diagram(a) (and(b) ). For example, the first factor of g swapsP2 withP3 in(a) , which means exchangingx2 andx4 in(c) (also(d) ). The aim of this way of labeling the Feynman diagram is to better classify different kinds of vertices and vertex lines, and to avoid confusion. This concept is in accordance with the first paragraph in Section 4.Figure 4. The actions of
gΓg,g, andggΓ. This is similar to the case in Fig. 3. Their first factors act on vertices, and the rest act on vertex lines. Similarly, the indices1 ,2 ,⋯ in the above expression for g andgΓ refer to the indices (P1 ,P2 , etc.) of the original blank Feynman diagram(a) .(a) Blank Feynman diagram(b) A marked Feynman diagramΓ (c) The action of g onψ : it swapsP2 withP3 in(a) , which swapsx2 withx3 in(b) . Still1 ,2 , and3 in the first factor of g correspond toP1 ,P2 , andP3 , respectively.(d) The action ofgΓ onψ : the first factor exchangesP1 withP2 (x2 andx1 in(d) ), and the last factor exchanges vertex lineα3 with lineβ3 (2 and1 ofx3 in(d) ).(e) The action ofgΓg onψ . Based on(c) , further action ofgΓ ongψ , which means swappingx1 withx3 in(c) (P1 andP2 in(a) ), and swapping1 with2 ofx3 in(c) (vertex lineα3 and lineβ3 in(a) ).(f) The action ofggΓ onψ . Based on(d) , further action of g ongΓψ means to swapx1 andx3 in(d) (P2 andP3 in(a) ). Here(c) and(e) are topologically equivalent, while(c) and(f) are topologically unequivalent. ThusgΓgψ andgψ correspond to the same term in Wick's expansion, but(ggΓ)ψ andgψ do not.We denote
L(Γ) as the set of internal lines ofΓ andV(Γ) as the set of vertices ofΓ . By Wick's expansion of Eq. (2) we get∑Γ′(n)∫d4x1⋯d4xn∏Labσ∈L(Γ)˜ΔabσF(xa−xb)∏Va∈V(Γ)Pa.
(4) Here
Pa is a numerical function depending on vertexVa , and˜ΔabσF(xa−xb) is a polynomial of derivatives of the Feynman propagator with respect to 4-coordinatesxa andxb . The summation is over all topologically different marked n point Feynman diagramsΓ′(n) . After integrating over allxi′s , one obtains an integral expression in momentum space:S=∑n1n!∑Γ′(n)∫∏Labσ∈L(Γ)d4labσ(2π)4ΔabσF(labσ)∏Va∈V(Γ)×{Pa({lab′σ′})×(2π)4δ4(∑b′σ′lab′σ′−qa)} ,
(5) where
Pa is a polynomial oflab′σ′ (the momentum of lineLab′σ′ ),qa is the total incoming momentum atVa , whileΔF is the Feynman propagator in momentum space. In Eq. (4) the(−i)n factor is absorbed into the∏Pa term. After integration of the momenta to eliminateδ functions, each term of connectedΓ′(n) leaves an overall(2π)4δ4(Σqa) function and an integral over independent momenta, since for the same blank Feynman diagram, integrals corresponding to different marked Feynman diagrams are equal. The contribution of the S matrix for a connected blank Feynman diagramΓ(n) with n vertices is:SΓ(n)=1n!|G||GΓ|∫dk1⋯dkl(2π)lI0Γ(k)(2π)4δ4(∑aqa)=1n!|Gn|×|∏iGi||GΓ|∫dk1⋯dkl(2π)l×I0Γ(k)(2π)4δ4(∑aqa) .
(6) Thus one has:
SΓ(n)=1|GΓ|×∫dk1⋯dkl(2π)l|∏iGi|×I0Γ(k)(2π)4δ4(∑aqa)≡1|GΓ|×∫dk1⋯dkl(2π)lIΓ(k)(2π)4δ4(∑aqa) .
(7) The factor
(−i) is absorbed into each vertex ofI0Γ(k) , while the term|∏iGi| is absorbed into each vertex ofIΓ(k) .The total S matrix is then:
S=∑n∑Γ(n)SΓ(n) .
(8) The second summation is over all blank Feynman diagrams with n vertices.
In Eq. (6), the momenta
k1⋯kl are essentially those remaining after eliminating theδ functions. The set{k1⋯kl} is denoted by k. -
In the following, we deal mainly with diagrams which are connected. If we cut any internal line of a connected diagram and the diagram is still connected, we call such a diagram a proper diagram. A proper diagram
γ with superficial dimension [7-9]d(γ)⩾0 is called a renormalization part [4].We have for a blank proper Feynman diagram
Γ (7):SΓ(n)=1|GΓ|∫dk1⋯dkl(2π)lIΓ(k)(2π)4δ4(∑qa) .
(9) A Feynman integral associated with a proper Feynman diagram is defined as:
JΓ=∫∏Labσ∈L(Γ)d4labσ(2π)4∏Va∈V(Γ){(2π)4δ4(∑labσ−qa)}IΓ(q,l)=∫dk1⋯dkl(2π)lIΓ(k,q)×(2π)4δ4(∑aqa) ,
(10) where
V(Γ) is the set of vertices inΓ ,qa is the total incoming momentum atVa ,L(Γ) is the set of internal lines ofΓ , and the lineLabσ is one of the lines from a vertexVa to a vertexVb with 4-momentalabσ . The integrandIΓ isIΓ=∏Labσ∈L(Γ)ΔabσF∏Va∈V(Γ)Pa ,
(11) where
ΔF is the Feynman propagator, andPa is a polynomial of{labσ} determined by the vertex ofHI(x) .For further derivation, we need to specify the integral parameters
ki′s . Due to theδ functions at each vertex, we have nonhomogeneous linear equations for the incoming 4-momentaqa at the vertexVa ,∑bσlabσ=qa ,a=1,⋯,n,
(12) where
labσ=−lbaσ
(13) is the 4-momentum on one of the lines
Labσ from the vertexVa to the vertexVb . The solution is not unique for∑Va∈Γqa=0 .
(14) When
Γ is connected, we define a solution for which the quantity [4]M=∑Labσ∈L(Γ)l2abσ=∑(qabσ)2
(15) is minimal. We call it a canonical distribution of incoming momenta
{qa} forΓ . One can prove thatqabσ is a linear combination ofqa′s . Then generallylabσ consists of two parts. One part comes fromqa , and the other part comes from integral variables likekabσ (they are momenta that form inner loop flows in the connected diagram)labσ=qabσ+kabσ.
So here
kabσ satisfies the homogeneous linear equations (compared with Eqs. (12) and (13)):∑bσkabσ=0 ,kabσ=−kbaσ,a=1,⋯,n.
(16) It is important to point out that for a subdiagram
γ⊂Γ with a lineLabσ∈γ (alsoLabσ∈Γ ),qγabσ andqΓabσ (we use the superscript to indicate the belonging) are different. From the linearity ofqabσ with respect toqa , we have the following statement.Proposition 1. The difference of
qγabσ andqΓabσ is the canonical distributionΔqa=−∑bσkΓabσ,Labσ∈L(Γ),Labσ∉L(γ),
(17) and thus they are linear functions of
kΓabσ of lines inΓ .We can choose (not arbitrarily) some lines of
Γ to fix the parameterskΓabσ and to further fix the solution. WhenΓ has m loops, we can fix m lines. Then the integralsk1⋯kl can be chosen as(kΓabσ)μ,μ=0,1,2,3 (forϕ4 theory) of these lines. Forγ1,⋯γc⊂Γ and whenγ1,⋯γc are disjoint diagrams, we can choose them in such a way that when the chosen lines are inγτ , they also form independent lines in eachγτ . In this case, each subdiagramγτ can be reduced to a vertex. We denote the lines in the reduced diagram˜Γ as:L˜Γabσ∈Γ/γ,Lγτa′b′σ′∈γτ.
Then due to proposition 1, we can properly choose independent lines such that the integral over independent momenta is:
∫ind∏d4l˜Γabσc∏τ=1ind∏d4lγτa′b′σ′=∫ind∏d4k˜Γabσc∏τ=1ind∏d4kγτa′b′σ′ .
(18) The Feynman integral (10) is generally divergent at large k for
Γ with loops. The BPHZ renormalization scheme [7,8] gives a finite partRΓ=RΓ(k,qa) such that the renormalized integral converges:FΓ=∫dk1⋯dkl(2π)lRΓ(k,q)(2π)4δ4(∑aqa) .
(19) Zimmermann proved the convergence of Eq. (19) [4] by an application of Weinberg's theorem [7]. The integrand
RΓ(k,q) of Eq. (19) is defined by:RΓ(k,q)=IΓ(k,q)+∑γ1⋯γcIΓ/γ1⋯γc(kq)c∏τ=1Oγτ(kγτ,qγτ) .
(20) The sum is over all sets of renormalization parts
γτ ,τ=1,⋯,c , which are mutually disjoint. The parametersk={k1⋯kl} are essentially the remaining momenta after integration eliminates theδ functions in Eq. (4). In Eq. (20) we have definedIΓ/γ1⋯γτ=IΓ/∏cτ=1Iγτ , which is only determined by vertices and lines contained inΓ but not in anyγτ . The functionsOγ are recursively defined asOγ(kγ,qγ)=−td(γ)qγ{Iγ(kγ,qγ)+′∑γ′1⋯γ′c′Iγ/γ′1⋯γ′c′(kγ,qγ)×c∏τ=1Oγτ(kγτ,qγτ)} .
(21) The sum extends over all sets of mutually disjoint renormalization parts
{γ′1,⋯,γ′c′} (assuming they are subdiagrams ofγ ), but does not include{γ} . The functionOγ is a Taylor series with respect to the incoming independent momentumqγ≡{qa} ,tdqf=d∑l=01l!∑q1⋯qlq1⋯ql(∂∂q1⋯∂∂qlf(g))q=0 .
(22) The sum of
qj extends over all components of independent incoming 4-momenta of{qa} .Since
Oγ is recursively defined, one must calculate it step by step, from smaller interior renormalization parts to bigger outer ones. Due to the fact that forγ′⊂γ , we have proposition 1, we get:qγ′a′b′σ′=qγa′b′σ′+linear combinations of {kγabσ}kγ′a′b′σ′=kγa′b′σ′+linear combinations of {kγabσ} .
(23) One can conclude that
O_{\gamma} is a polynomial of independent components of\{q^{\gamma}_{a}\} . The coefficients are functions of\{k^{\gamma}_{ab\sigma}\} . Equation (19) can be rewritten as\begin{aligned}[b] F_{\Gamma} =& \Bigg\{\int\frac{{\rm d}k_1\cdots {\rm d}k_l}{(2\pi)^l}I_{\Gamma}+\sum\limits_{\gamma_1\cdots\gamma_c}\int\frac{{\rm d}k_1\cdots {\rm d}k_l}{(2\pi)^l}I_{\Gamma/\gamma_1\cdots\gamma_c}(k,q)\\&\times\prod\limits^c_{\tau = 1}O_{\gamma_{\tau}}(k^{\gamma_\tau},q^{\gamma_{\tau}}) \Bigg\}(2\pi)^4\delta^4\left(\sum_a q_a\right) \equiv J_{\Gamma}+\sum\limits_{\gamma_1\cdots\gamma_c}\widetilde{J}_{\Gamma/\gamma_1\cdots\gamma_c}\ , \end{aligned}
(24) where
\widetilde{J}_{\Gamma/\gamma_1\cdots\gamma_c} are counterterms ofJ_{\Gamma} . From Eq. (18) we have:\begin{aligned}[b] \widetilde{J}_{\Gamma/\gamma_1\cdots\gamma_c} = & \int\frac{{\rm d}k_1\cdots {\rm d}k_l}{(2\pi)^l}I_{\Gamma/\gamma_1\cdots\gamma_c}(k,q)\prod\limits^c_{\tau = 1}O_{\gamma_\tau}(k^{\gamma_\tau},q^{\gamma_\tau})(2\pi)^4\delta^4 \left(\sum\limits_a q_a\right) = \int\prod\limits^{ind}_{L_{ab\sigma}\in{\Gamma/\gamma_1\cdots\gamma_c}}\frac{{\rm d}^4k^{\Gamma}_{ab\sigma}}{(2\pi)^4}I_{\Gamma/\gamma_1\cdots\gamma_c}(k^{\Gamma},q) \prod\limits^{c}_{\tau = 1}\int\prod\limits^{ind}_{L_{a'b'\sigma'}\in{\gamma_\tau}}\\ & \times\frac{{\rm d}^4k^{\gamma_\tau}_{a'b'\sigma'}}{(2\pi)^4}O_{\gamma_\tau}(k^{\gamma_\tau},q^{\gamma_\tau})(2\pi)^4\delta^4\left(\sum\limits_a q_a\right) \equiv \int\prod\limits^{ind}_{L_{ab\sigma}\in{\Gamma/\gamma_1\cdots\gamma_c}}\frac{{\rm d}^4k^{\Gamma}_{ab\sigma}}{(2\pi)^4}I_{\Gamma/\gamma_1\cdots\gamma_c} (k^{\Gamma},q)\times \prod\limits^{c}_{\tau = 1}Q_{\gamma_\tau}(2\pi)^4\delta^4\left(\sum\limits_a q_a\right)\ , \end{aligned}
(25) where
\prod\limits^{ind} extends to all independent lines. The integralQ_{\gamma_\tau} is divergent, but it is convergent after regularization. Thus we may regardQ_{\gamma_\tau} as a new vertex which corresponds to a renormalization part\gamma_\tau . We call it a vertex of reduction.Equation (25) is a Feynman integral for Feynman diagram
\widetilde{\Gamma} . It is a reduced diagram from\Gamma , in which all lines and vertices of\gamma_{\tau} are shrunk to a vertex\widetilde{Q}_{\gamma_\tau} .Can we introduce new interactions (with divergent coefficients) into interaction Lagrangian
\Delta{\cal{L}}_I which gives precisely these counterterms to a given order? This is a consistency problem for the BPHZ renormalization scheme. The answer is affirmative. We would like to prove this in the coming sections. It can be proved that the functionQ_{\gamma_{\tau}} has symmetry of\gamma_\tau (Appendix A). We next study the symmetry group of\widetilde{\Gamma} denoted byG_{\widetilde{\Gamma}} . -
In the derivation in Section II, one group element "
g^i " ing\in G for changing the number of equivalent vertex lines in a vertex is attached to a certain pointP_a of a blank Feynman diagram. Generally, the vertices are of different type at different points. For the symmetry group of a Feynman diagram\Gamma , we can alternatively attach such a group element to a running coordinatex_j . This is because ing_{\Gamma}\in G_{\Gamma} , we always have the same type of vertices for the samex_j before and after reassigning. In this way, we may regard the element ofg_{\Gamma}\in G_{\Gamma} as a continuous mapping of a Feynman diagram into itself. This maps a vertex (including the vertex point together with all its vertex lines) into another vertex of the same type. In the following, we always referg_{\Gamma} to the second meaning.We define a set A as the union of all
{\gamma_{\tau}}'s A = \bigcup\limits^{c}_{\tau = 1}\gamma_\tau\ ,
(26) which includes all internal lines and vertex points of all
{\gamma_{\tau}}'s . We define setsB_{\tau} as all internal lines and vertex points of\gamma_{\tau} ; actually it is\gamma_{\tau} itself.In the group
G_{\Gamma} , those elements which do not change the set A (i.e. which map A to A) form a subgroup ofG_{\Gamma} . We denote it asG^{(1)}_{\Gamma}(\gamma_1\cdots\gamma_c) . This is because the set of such elements is closed under multiplication and inversion. Similarly, we have a subgroupG^{(2)}_{\Gamma}(\gamma_1\cdots\gamma_c) = \prod^{c}_{\tau = 1} G^{(2)}_{\Gamma}(\gamma_{\tau}) , whereG^{(2)}_{\Gamma}(\gamma_{\tau})\subset G_{\Gamma} maps eachB_{\tau} = \gamma_{\tau} to itself, andG^{(2)}_{\Gamma}(\gamma_{\tau}) keeps all lines and vertex points outside\gamma_{\tau} invariant (identity). Since the mapping keeps the connection relation between the vertex points and the vertex lines of any vertex,G^{(2)}_{\Gamma}(\gamma_{\tau}) also keeps the "boundary points" which connect exterior lines of\gamma_{\tau} invariant. We can prove thatG^{(2)}_{\Gamma} is a normal subgroup ofG^{(1)}_{\Gamma} since an element ofG^{(2)}_{\Gamma} does not change the exterior of A. Then we can further prove that the symmetry group of the reduced Feynman diagram\widetilde{\Gamma} is the quotient groupG_{\widetilde{\Gamma}} = G^{(1)}_{\Gamma}/G^{(2)}_{\Gamma}\ .
(27) This is because it maps the reduced vertices to themselves with their symmetry. We use (Fig. 5) to illustrate the relation between a Feynman diagram
\Gamma and its reduced Feynman diagrams. The configuration of\Gamma with\gamma_1\cdots\gamma_c can be mapped to\Gamma with\gamma_1',\cdots,\gamma_c' (in the same blank configuration) by each element ofG_{\Gamma} . Those elements which keep the set A invariant form subgroupG^{(1)}_{\Gamma} . Letg_{\Gamma}\in G_{\Gamma} , andh\in G^{(1)}_{\Gamma} . Then for a certaing_{\Gamma} and anyh\in G^{(1)}_{\Gamma} , the elementg' = g_{\Gamma}h maps the set A to a setA' = \bigcup_{\tau}\gamma_{\tau}' . Thus the cosetG_{\Gamma}/G^{(1)}_{\Gamma} is characterized by different "A' "s which are also sets of disjoint renormalization parts of\Gamma . Therefore the renormalized integrand (20) can be written as:Figure 5. Diagram of
\Gamma with\widetilde{\Gamma}'s . (a) Diagram of\Gamma whereA,B,C,D are renormalization parts.(b) Choose A and C as{\gamma_{\tau}}'s to get\check{\Gamma}_1 .(c) \widetilde{\Gamma}_1 from\check{\Gamma}_1 (shrink all lines and points of A (and C) into a point).(d) Another configuration equivalent to(b) , which gives\check{\Gamma}_1' .(e) \widetilde{\Gamma}_1' from\check{\Gamma}_1' , equivalent to\widetilde{\Gamma}_1 .(f) \{\gamma_{\tau}\} = \{A,B,C,D\} gives\check{\Gamma}_2 .(g) \widetilde{\Gamma}_2 . R_{\Gamma}(k,q) = I_{\Gamma}(k,q)+\sum\limits_{\check{\Gamma}}\!_1\sum\limits_{\rm coset}\!_2I_{\Gamma/\gamma_1\cdots\gamma_c} (k,q)\prod\limits^{c}_{\tau = 1}O_{\gamma_\tau}(k^{\gamma_\tau},q^{\gamma_\tau})\ .
(28) The first sum extends over all configurations of
\check{\Gamma}:\Gamma with different\gamma_1\cdots\gamma_c , which are topologically different from each other. The second sum extends over configurations produced by coset elementsG_{\Gamma}/G^{(1)}_{\Gamma} (acting on\Gamma ). After integration one has the renormalized integral from Eq. (24):\begin{aligned}[b] F_{\Gamma} =& J_{\Gamma}+\sum\limits_{\check{\Gamma}}\!_1m_{\check{\Gamma}}\int\prod\limits^{ind}_{L_{ab\sigma}\in\Gamma/\{\gamma_1\cdots\gamma_c\}} \frac{{\rm d}^4k^{\Gamma}_{ab\sigma}}{(2\pi)^4}I_{\Gamma/\gamma_1\cdots\gamma_c}(k^{\Gamma},q)\\&\times\prod\limits^{c}_{\tau = 1}O_{\gamma_\tau} (2\pi)^4\delta^4\left(\sum_a q_a\right), \end{aligned}
(29) where
m_{\check{\Gamma}} = \dfrac{|G_\Gamma|}{|G^{(1)}_{\Gamma}|} . In the derivation, we need the fact thatQ_{\gamma_{\tau}} is symmetric under the mapping ofG_{\Gamma} . -
When we introduce new reduced vertices
\hat{Q}_{\gamma} to\Delta{\cal{L}}_I = -\Delta{\cal{H}}_I , the S matrix gets many new terms due to these new vertices. In a reduced vertex, the order of the original perturbation constant is more than one since it contains more than one original vertex. Thus if we collect terms according to the orders of original perturbation constants in the S matrix, say, n-th order, the term with respect to a reduced diagram\widetilde{\Gamma} may be with a "perturbation order" m,m<n . We denote such a termS_{\widetilde{\Gamma}(n,m)} . Assume we have all reduced vertices associated with all\gamma_{\tau} (renormalization part) forn_{\gamma_\tau}\leqslant n , and assume these vertices have the symmetry of\gamma_\tau . We denote the exact symmetry group for exterior lines of such vertices asg_{O_{\gamma_\tau}} (see Appendix A for details). Due to the perturbation theory, there will be a term associated with\widetilde{\Gamma} in the S matrix. From Eq. (6) we have:\begin{aligned}[b] S_{\widetilde{\Gamma}(n,m)} =& \frac{1}{m!}\frac{ |G_m|\times|\prod^1 G_{O_{\gamma_\tau}}|\times|\prod^2G^i|}{|G_{\widetilde{\Gamma}}|}\\&\times\int\frac{{\rm d}k_1\cdots {\rm d}k_s}{(2\pi)^s}I^0_{\widetilde{\Gamma}}(2\pi)^4\delta^4 \left(\sum\limits_a q_a \right), \end{aligned}
(30) where
G_m is the permutation groupS_m ,\prod^{1} extends over all reduced vertices, and\prod^{2} extends over the remaining original vertices in\widetilde{\Gamma} . Since|G_m| = m! , we obtain:\begin{aligned}[b] S_{\widetilde{\Gamma}(n,m)} =& \frac{ |\prod^1G_{O_{\gamma_\tau}}|\times|\prod^2G^i|}{|G_{\widetilde{\Gamma}}|} \\&\times\int\frac{{\rm d}k_1\cdots {\rm d}k_s}{(2\pi)^s}I^0_{\widetilde{\Gamma}}(2\pi)^4\delta^4\left(\sum\limits_a q_a\right)\ . \end{aligned}
(31) Together with
S_{\Gamma(n)} , this gives:\begin{aligned}[b] S_{\Gamma(n)}+S_{\widetilde{\Gamma}(n,m)} =& \frac{ |\prod\limits^2G^i|\times|\prod\limits^3G^i|}{|G_{\Gamma}|} \int\frac{{\rm d}k_1\cdots {\rm d}k_l}{(2\pi)^l}I^0_{\Gamma}(2\pi)^4\delta^4\left(\sum_a q_a\right)+ \frac{ |\prod\limits^2G^i|\times|\prod\limits^3G^i|}{|G_{\Gamma}|}\times \frac{ |G_{\Gamma}|\times|\prod\limits^1G_{O_{\gamma_\tau}}|}{ |G_{\widetilde{\Gamma}}|\times|\prod\limits^3G^i|}\\&\times\int\frac{{\rm d}k_1\cdots {\rm d}k_s}{(2\pi)^l}I^0_{\widetilde{\Gamma}}(2\pi)^4\delta^4\left(\sum\limits_a q_a\right)\ , \end{aligned}
(32) where
\prod^{3} extends over all vertices in\bigcup_{\tau}\gamma_{\tau}\subset\Gamma . Thus one has|\prod^{2}G^{i}|\times|\prod^{3}G^{i}| = \prod G^{i} , which extends over all vertices of\Gamma . The termS^{(n,m)}_{\widetilde{\Gamma}} should match the counterterms in the BPHZ renormalization scheme. We denote the vertex of\gamma_\tau as\int\prod_{ind}k^{\gamma_\tau}_{ab\sigma}O^0_{\gamma_\tau} . Based on Eq. (25), Eq. (32) can be written as:\begin{aligned}[b] S_{\Gamma(n)}+S_{\widetilde{\Gamma}(n,m)} =& \frac{|\prod G^i|}{|G_{\Gamma}|}\int\frac{{\rm d}k_1\cdots {\rm d}k_l}{(2\pi)^l}\left\{I^0_{\Gamma}+\frac{|G_{\Gamma}|\prod^1|G_{O_{\gamma_\tau}}|}{|G_{\widetilde{\Gamma}}|\prod^3|G^i|} I^0_{\Gamma/\{\gamma_1\cdots\gamma_c\}}\prod\limits_{\tau}O^0_{\gamma_\tau}\right\}(2\pi)^4\delta^4\left(\sum_a q_a\right)\\ = & \frac{|\prod G^i|}{|G_{\Gamma}|}\int\frac{{\rm d}k_1\cdots {\rm d}k_l}{(2\pi)^l}\left\{I^0_{\Gamma}+\frac{|G_{\Gamma}||G^{(2)}_{\Gamma}|\prod_{\tau}|G_{O_{\gamma_\tau}}|}{|G^{(1)}_{\Gamma}|\prod^3|G^i|}I^0_{\Gamma/\{\gamma_1\cdots\gamma_c\}} \prod O^0_{\gamma_\tau}\right\}(2\pi)^4\delta^4\left(\sum_a q_a\right)\\ =& \frac{|\prod G^i|}{|G_{\Gamma}|}\int\frac{{\rm d}k_1\cdots {\rm d}k_l}{(2\pi)^l}\left\{I^0_{\Gamma}+\frac{|G_{\Gamma}|\prod_\tau|G_{\gamma_\tau}|\prod_{\tau}|G_{O_{\gamma_\tau}}|}{|G^{(1)}_{\Gamma}|\prod^3|G^i|}I^0_{\Gamma/\{\gamma_1\cdots\gamma_c\}} \prod O^0_{\gamma_\tau}\right\}(2\pi)^4\delta^4\left(\sum_a q_a\right), \end{aligned}
(33) where
|\dfrac{G^{(2)}_{\Gamma}}{G^{(1)}_{\Gamma}}| = \dfrac{1}{|G_{\widetilde{\Gamma}}|} , and|\dfrac{G_{\Gamma}}{G^{(1)}_{\Gamma}}| is the number of cosets, which indicates the number of different\widetilde{\Gamma}'s whose configurations are the same as\widetilde{\Gamma} .Taking into account all reduced diagrams of
\Gamma , we have:\begin{aligned}[b] S_{\Gamma(n)}+\sum\limits_{l}S_{\widetilde{\Gamma}_l(n,m)} = \frac{\prod|G^i|}{|G_{\Gamma}|}\int\frac{{\rm d}k_1\cdots {\rm d}k_l}{(2\pi)^l}\left\{I^0_{\Gamma}+\sum\limits_l\frac{\prod_{\tau}|G_{O_{\gamma_\tau}}||G_{\gamma_\tau}|}{\prod^3|G^{i}|}\right\}\times I^0_{\Gamma/\{\gamma_1\cdots\gamma_c\}}\prod\limits_{\tau}O^0_{\gamma_\tau}(2\pi)^4\delta^4\left(\sum\limits_a q_a\right)\ , \end{aligned}
(34) where l is specified by the set
\{\gamma_1\cdots\gamma_c\} of mutually disjoint renormalization parts in\Gamma .Next we absorb all factors
|G^{i}| of symmetry groups of the original vertices into these vertex constants, and require\frac{1}{|G_{\gamma_\tau}||G_{O_{\gamma_\tau}}|}O_{\gamma_\tau} = \frac{1}{|G_{\hat{\gamma}_\tau}|}O_{\gamma_\tau} = O^0_{\gamma_\tau}\ ,
(35) which corresponds to
O^0_{\gamma_\tau} in Eq. (34). We have:\begin{aligned}[b]& S_{\Gamma(n)}+\sum\limits_{l}S_{\widetilde{\Gamma}_l(n,m)}\\ =& \frac{1}{|G_{\Gamma}|}\int {\rm d}k_1\cdots {\rm d}k_l\left\{I_{\Gamma}+\sum\limits_{\{\gamma_1\cdots\gamma_c\}}I_{\Gamma/\{\gamma_1\cdots\gamma_c\}}\prod\limits_{\tau}O_{\gamma_{\tau}}\right\}\\&\times(2\pi)^4\delta^4\left(\sum\limits_a q_a\right) = \frac{1}{|G_\Gamma|}F_{\Gamma}\ . \end{aligned}
(36) This is just the BPHZ formula. The corresponding operator of
\hat{O}^0_{\gamma_\tau} in Eq. (35) is given in Appendix A, and it only depends on the renormalization part\gamma_{\tau} . -
In conclusion, we have considered the procedure for introducing a new vertex
\hat{Q}^0_{\gamma_\tau} into\Delta{\cal{L}}_I for all renormalization parts\gamma_\tau (proper diagram withd(\gamma_\tau)\geqslant0 ) with number of verticesm\leqslant n . When we collect all terms of original parameter orderm'\leqslant n (the order for original perturbation parameters) in the S matrix, then it will automatically give the counterterms of Feynman integrals as the BPHZ scheme requires for any Feynman diagram\Gamma .Thus the BPHZ scheme is consistent with adding counterterms in
\Delta{\cal{L}}_{I} . -
All Feynman diagrams in this paper were drawn by the program JaxoDraw [18]
-
Similar to
G_\Gamma , we define a mapping groupG_{\hat{\gamma}} for the Feynman diagram\gamma . In this mapping, we allow the external lines to be mapped to equivalent lines of another vertex of the same type. This is different fromG_\Gamma , where the external lines always keep fixed. We find that the groupG_\gamma is a normal subgroup ofG_{\hat{\gamma}} . We denote\hat{\gamma} with external lines of\gamma belonging to{\cal{L}}(\hat{\gamma}) . We have:\tag{A1} O_{\gamma} = -t^{\gamma}\left(I_{\gamma}+\sum\limits_{\{\gamma_1\cdots\gamma_c\}}\!'I_{\gamma/\{\gamma_1\cdots\gamma_c\}}\prod_{\tau}O_{\gamma_\tau}\right),
and
\tag{A2} Q_{\gamma} = \int\prod\limits^{ind}\frac{{\rm d}^4k^{\gamma}_{ab\sigma}}{(2\pi)^4}O_{\gamma},
due to Eq. (25). Assume for all proper renormalization parts
\gamma_{\tau}\subset\gamma ,O_{\gamma_{\tau}} are symmetric functions with the symmetry of\gamma_{\tau} . That means, underG_{\hat{\gamma}_{\tau}} , the value ofQ_{\gamma_{\tau}} is invariant. We can write Eq. (A2) as:\tag{A3} Q_{\gamma} = -t^{\gamma}V_{\gamma}\ .
Here
V_{\gamma} is defined as:\tag{A4} \begin{aligned}[b] V_{\gamma}\delta^4\left(\sum q^{\gamma}_a\right) =& \int\prod\limits_{L_{ab\sigma}\in{\cal{L}}(r)}\frac{{\rm d}^4k^{\gamma}_{ab\sigma}}{(2\pi)^4}\check{I}^0_{\gamma}\\ & + \sum\!'\int\prod\limits_{L_{ab\sigma}\in{\cal{L}}(\gamma/\{\gamma_1\cdots\gamma_c\})}\frac{{\rm d}^4k^{\gamma}_{ab\sigma}}{(2\pi)^4} \check{I}^0_{\gamma/\{\gamma_1\cdots\gamma_c\}}\\&\times\prod\limits_{\tau}\left\{Q_{\gamma_\tau}(2\pi)^4\delta^4\left(\sum_{V_a\in\gamma_{\tau}}q^{\gamma_{\tau}}_{a}\right)\right\}, \end{aligned}
in which:
\tag{A5} \check{I}^0_{\gamma} = \prod\limits_{L_{ab\sigma}\in{\cal{L}}(\gamma)}\Delta^{ab\sigma}_{F}\prod\limits_{V_a\in{\cal{V}}(\gamma)}\left\{P_a(2\pi)^4\delta^4\left(\sum\limits_{L_{ab'\sigma'}\in{\cal{L}}(\hat{\gamma})}l_{ab'\sigma'}\right)\right\}\ .
The situation is similar for
\check{I}_{\gamma/\{\gamma_1\cdots\gamma_c\}} , and Eq. (A4) is invariant under the mapping ofG_{\hat{\gamma}} . Thus we can always chooseV_{\gamma}(s) which is invariant underG_{\hat{\gamma}} . On the other hand, we can also chooseV_{\gamma}(ind) which is only a function of independentq^{\gamma}_a . Due to\delta^4(\sum^a q_a) in Eq. (A4), the number of independent{q_{a}^{\gamma}} is less than the number of all{q_{a}^{\gamma}} . We have:\tag{A6} V_{\gamma}(ind)\delta^4\left(\sum q^{\gamma}_a\right) = V_{\gamma}(s)\delta^4\left(\sum q^{\gamma}_a\right).
One can show by truncating the Taylor series in Eq. (22) that:
\tag{A7} -t^{\gamma}{(ind)}V_{\gamma}(ind)\delta^4\left(\sum q^{\gamma}_a\right) = -t^{\gamma}V_{\gamma}(s)\delta^4\left(\sum q^{\gamma}_a\right),
where
-t^{\gamma}(ind)V_{\gamma}(ind) is theQ_\gamma in Eq. (A2). We then prove thatQ_\gamma can also be chosen as a symmetric (invariant) function under the mapping ofG_{\hat{\gamma}} . It is fixed underG_\gamma . Thus as a reduced vertex, we have the symmetry group,\tag{A8} G_{O_{\gamma}} = G_{\hat{\gamma}}/G_{\gamma}\ .
Here we make a remark about invariants. If the momenta of external lines change following the mapping, the function
Q_{\gamma} is invariant. For example, under the mapping depicted in Fig. A1, if the functionf(x_1, x_2,x_3) = f(x_1', x_2',x_3') , we say function f is invariant under the mapping g.We denote
Q_{\gamma} as\tag{A9} \delta^4\left(\sum\limits_{a}q^{\gamma}_{a}\right)Q_{\gamma} = -t^{\gamma}V_{\gamma}(\{q^{\gamma}_{a}\})\delta^4\left(\sum\limits_{a}q_a^{\gamma}\right)\ .
In the integral equations (1) and (2), the derivative
\left(-{\rm i}\dfrac{\partial}{\partial x_a}\right) operator produces a momentumq_a factor. Thus the functionQ_\gamma can be realized in\Delta{\cal{L}}_I by the operator proportional to\tag{A10} \begin{aligned}[b] \hat{Q}_{\gamma} =& {\rm i}\sum\limits^{d(\gamma)}_{l = 0}\frac{1}{l!}\sum\limits_{i_i\cdots1_l}\left(-{\rm i}\frac{\partial}{\partial x_{i_1}}\right)\cdots\left(-{\rm i}\frac{\partial}{\partial x_{i_l}}\right)\\&\times\left(\prod\limits_{i}\varphi_{a_i}(x_a)\prod\limits_{j}\varphi_{b_j}(x_b)\cdots\right)\Big|_{x_a = x_b = \cdots = x}\\& \times \frac{\partial}{\partial q_{i_1}}\cdots\frac{\partial}{\partial q_{i_l}}V_{\gamma}(q^{\gamma}_{a})|_{q = 0}, \end{aligned}
where
i_1,\cdots,i_l in summation run over alla = 1,\cdots,n_{\gamma} ,(n_\gamma = \text{number of vertices in }\gamma) , and\mu = 0,1,2,3. In Eq. (A10), fields\{\varphi_{a_i},i = 1,\cdots,l_a\} produce external lines at the vertexV_a,\cdots .From Eq. (35) the operator of the reduced vertex introduced in
\Delta{\cal{L}}_I is then\tag{A11} \hat{Q}^0_{\gamma_\tau} = \frac{1}{|G_{O_{\gamma_{\tau}}}||G_{\gamma_\tau}|}\hat{Q}_{\gamma_{\tau}} = \frac{1}{|G_{\hat{\gamma}_{\tau}}|}\hat{Q}_{\gamma_{\tau}}\ .
G_{{\gamma}_\tau} is a normal subgroup ofG_{\hat{\gamma}_\tau} , see Fig. A2 below.Figure A2. A renormalization part
\gamma_\tau and associatedG_{\hat{\gamma}_\tau} andG_{{\gamma}_\tau} .(a) The lines belonging to\hat{\gamma}_\tau .(b) The lines (only internal lines!) belonging to\gamma_{\tau} .(c) A mappingg_{\hat{\gamma}_\tau} , including the change of points and external lines of\gamma_\tau . This mapping does not belong toG_{\gamma_\tau} .(d) A mappingg_{\gamma_{\tau}}\in G_{\gamma_{\tau}}\subset G_{\hat{\gamma}_{\tau}} .
