Renormalization Group Equation for Weakly Power Counting Renormalizable Theories

We study the renormalization group flow in weak power counting (WPC) renormalizable theories. The latter are theories which, after being formulated in terms of certain variables, display only a finite number of independent divergent amplitudes order by order in the loop expansion. Using as a toolbox the well-known SU(2) non linear sigma model, we prove that for such theories a renormalization group equation holds that does not violate the WPC condition: that is, the sliding of the scale $\mu$ for physical amplitudes can be reabsorbed by a suitable set of finite counterterms arising at the loop order prescribed by the WPC itself. We explore in some detail the consequences of this result; in particular, we prove that it holds in the framework of a recently introduced beyond the Standard Model scenario in which one considers non-linear St\"uckelberg-like symmetry breaking contributions to the fermion and gauge boson mass generation mechanism.


I. INTRODUCTION
Arguably, the most general way to look at the (perturbative) renormalization of quantum field theories is the one introduced two decades ago by Weinberg and Gomis [1]. In this modern look at the subject, the boundary between what constitutes a renormalizable or a non-renormalizable theory gets blurred. Indeed, in [1] it was proven that, if one includes in the tree-level vertex functional all possible Lorentz-invariant monomials in the fields, the external sources and their derivatives, while respecting, at the same time, the symmetries of the theory (encoded in the Batalin-Vilkovisky master equation [2]), then it is possible to subtract all ultraviolet (UV) divergences by a redefinition of the tree-level parameters.
What distinguishes between the renormalizable/non-renormalizable cases is the stability of the classical action under radiative corrections. In fact, in the renormalizable case the finite number of terms already present in the tree-level action is sufficient to absorb all UV divergences irrespectively of the loop order. This ceases to be true for non-renormalizable theories, where no matter how many terms are added to the tree-level action, there will always exist a sufficiently high order operator in the loop expansion which will give rise to a UV divergence that cannot be absorbed into the classical action. Such an action would therefore be unstable against radiative corrections, and one has to allow for infinitely many (symmetry preserving) terms in order to absorb all divergences.
A prototype non-renormalizable theory (or, said in the positive, renormalizable in the modern sense), is the d-dimensional non-linear sigma model (NLσM for short) in which the massless pseudoscalar "pion" fields φ a form, together with the scalar field φ 0 (with φ 0 = m d > 0), a chiral multiplet subjected to the (nonlinear) constraint φ 2 a + φ 2 0 = m 2 d ; in this way the global symmetry SU(2) L × SU(2) R is non-linearly realized 1 [see Eq. (2. 7)].
Already at the one-loop level this theory shows an infinite number of divergent one-particle irreducible (1-PI) amplitudes, which in turn make the consistent subtraction of UV divergences problematic. However, by embedding the global NLσM into a locally symmetric formulation in terms of a classical source corresponding to a certain (flat) connection J a µ , it was shown in [3] that one can acquire full control over the UV divergences of the model. Specifically, it can be shown that in such a theory a Weak Power Counting (WPC) ex-ists [4,5] which implies that the number of independent divergent amplitudes stays finite at each loop order n, even though it increases with n. These "ancestor" amplitudes, are not the ones associated with the quantum pion fields φ a , but are rather written in terms of the connection J a µ and the composite operator K 0 implementing the aforementioned non-linear constraint. Then, the presence of the so-called Local Functional Equation (LFE) together with a suitable change of field variables called "bleaching", allows one to generate from the ancestor amplitudes all the (infinite) descendant (off-shell) amplitudes involving the pion fields, uniquely fixing en route their divergent part at any order in the loop expansion [5].
Within the aforementioned Weinberg-Gomis approach to renormalization, one can interprete the WPC as a condition dictating at which loop order the coefficient of a particular monomial in the ancestor variables and their derivatives becomes non-vanishing. While it has been proven in [6] that the WPC allows for the definition of a symmetric (i.e., compatible with the symmetries of the theory) subtraction scheme, the question has remained open of whether a Renormalization Group (RG) equation compatible with the WPC exists. The problem can be stated in the following terms. WPC renormalizable theories possess a RG flow associated to the RG equation for ancestor amplitudes. Imagine then that a change in the scale µ of the radiative corrections is inequivalent to a rescaling of the coefficients of the ancestor amplitudes counterterms, at the given loop order prescribed by the WPC: this would imply that the RG flow mixes up the hierarchy of UV divergences encoded in the WPC, thus making it impossible to slide the scale µ between different energies.
As we will prove in this paper, fortunately this is not the case, as indeed the RG equation of ancestor amplitudes turns out to be compatible with the WPC. This is an important result which acquires particular relevance in the context of the non-linearly realized electroweak theory introduced in [7,8], in which the WPC has been used as a model building principle.
In particular, in this model the classical source J a µ is promoted to a local dynamical field responsible for generating part of the mass of the W ± and Z gauge bosons through the Stückelberg mechanism; this leads in turn to many unique features which cannot be found in theories describing physics beyond the Standard Model (e.g., it is impossible to add a scalar singlet without breaking the WPC, so that the minimal number of physical scalar resonances in the model is 4). As we will see, the result proven here for the compatibility between the RG flow and the WPC applies also in this case. This entails the possibility of evolving the scale µ in a mathematically consistent way, thus allowing to obtain predictions for the relevant observables applicable in different energy regimes and thus paving the way for a systematic study of their deviations from the expected SM results.
The paper is organized as follows. By using the aforementioned example of the ddimensional SU(2) NLσM in Sect. II we review the embedding of the model in a local formulation in terms of a flat connection, leading to the LFE and the WPC condition. The concepts of ancestor and descendant amplitudes as well as of bleached variables are introduced together with the corresponding symmetric scheme for subtracting the UV divergences. In Sect. III we derive the equation governing the RG flow in the local NLσM . Next, after specializing this equation to ancestor amplitudes, we discuss under which conditions the RG flow preserves the WPC; we then prove that at the one-loop level these conditions are indeed satisfied. Sect. IV contains the central result of the paper: there we prove the general theorem stating that the RG-flow of ancestor amplitudes preserves the WPC. In Sect. V we introduce the notion of weak stability, while in Sect. VI we derive the consequences of the RG flow theorem for the aforementioned non-linearly realized SU(2)×U(1) electroweak theory. Our conclusions and outlook are finally presented in Sect. VII.

A. Global NLσM
The d-dimensional action of a NLσM is conventionally written as where Ω = Ω(x) (with x the space-time coordinates) represents a matrix belonging to a symmetry group G. m d = m d/2−1 is the mass scale of the theory (m has dimension 1).
We will consider in particular the case where G is SU(N) (and then specialize to the case : In terms of a suitable basis of fields φ a (x), parametrizing the G matrix Ω, i.e., Ω(x) = Ω(φ a (x)), the action (2.1) reads The fields φ a will be generically referred to as "pion" fields. Geometrically they represent the coordinates of the group manifold G.
In what follows we will mostly deal with the case G ≡ SU(2), where one can set where a = 1, 2, 3, and τ a are the usual Pauli matrices. The action (2.3) will then contain non-polynomial derivative interactions involving the pion fields, reading It is exactly the presence of two derivatives in the interaction term that in d > 2 causes severe UV divergences, which are ultimately responsible for the non-renormalizability of the corresponding quantized theory.
When writing the action as in (2.1), it is immediate to show that the theory is invariant under a non-linearly realized SU(2) L × SU(2) R global symmetry: In terms of the pion fields, the left infinitesimal transformation of constant parameters ω a reads δφ 0 (x) = − 1 2 In matrix form the infinitesimal left transformation is For any unitary matrix Ω it is possible to define a flat connection 2 so that the action (2.1) can be cast in the form (2.10) Specialize now to the SU(2) case, and consider a local SU(2) L transformation on Ω; this will induce a gauge transformation on the flat connection F µ , namely Clearly, the action (2.10) is not invariant under these local transformations.
However, let us introduce an additional classical source J µ transforming as a gauge connection under the local SU(2) L group. At this point the difference F µ − J µ will transform in the adjoint representation so that the action is invariant under a local SU(2) L symmetry. In coordinates one has The local infinitesimal transformations are The global NLσM is embedded in the local formulation we have just provided. Specifically, the terms JF and J 2 are separately invariant under a global SU(2) L transformation (that is when the ω gauge parameters are kept constants); therefore we can set J directly to zero to obtain where Γ (0) is given by

19)
K 0 being an SU(2) L invariant source associated to the non-linear constraint (2.4). Thus, K 0 is associated with the auxiliary external field required to define the composite operator entering in the non-linear symmetry transformation. It plays the same role as the antifields [2] in gauge theories and the tree-level dependence of the vertex functional on K 0 is fixed by the form of the non-linear transformation (2.16): since the only composite operator entering in δφ a is φ 0 , one only needs an external source K 0 [5].
As the symmetry is non-anomalous, the LFE (2.18) is satisfied by the full vertex functional Γ: (2.20) In addition, notice that since the rhs of Eq. (2.20) is linear in the quantized fields and thus remains classical, this term will be present only at tree-level.
The LFE encodes at the quantum level the classical local SU(2) transformation, whose One of its main consequences is the separation of the 1-PI amplitudes into two classes. On the one hand, there are the amplitudes involving only the insertion of the SU(2) connection J a µ and of the source of the non-linear constraint K 0 : these are called ancestor amplitudes [3,4]. On the other hand, we have the so-called descendant amplitudes, i.e., those involving at least one external φ-leg. These amplitudes are not independent, as they are uniquely determined by the LFE once the ancestor amplitudes are known [3,4,9].

D. Weak Power Counting
Despite the fact that they do not involve external legs of the quantized fields of the theory φ a , the truly fundamental Green's functions of the local NLσM are the ancestor amplitudes. Such Green's functions display an UV behaviour that is significantly better than the one of their descendants: namely, there exists a choice of the tree-level action, compatible with the symmetries of the theory, such that only a finite number of divergent ancestor amplitudes arises order by order in the loop expansion. This property is dubbed the Weak Power-Counting (WPC) condition [4].
Indeed, one can show that in d-dimensions an n-loop ancestor amplitude G with N J (N K 0 ) external J a µ (K 0 ) legs has a superficial degree of divergence given by [4] Thus, at every-loop order only a finite number of superficially divergent ancestor amplitudes exists, i.e., the ones for which D(G) ≥ 0; obviously, the local NLσM is still non renormalizable, as Eq. (2.21) shows that as n grows bigger the number of UV divergent amplitudes increases. For example, in the 4-dimensional case Eq. (2.21) tells us that at one-loop the UV divergent amplitudes involve up to four external J a µ legs and/or two K 0 legs. Accordingly, the one-loop 1PI functional for these ancestor amplitudes reads where the dots stand for ancestor amplitudes that are not UV-divergent at one loop.

E. Bleached variables
Ancestor amplitudes per-se are not a solution of the LFE (2.20) as they carry no information for amplitudes involving pion fields. To achieve this, it is necessary to introduce invariant combinations in one-to-one correspondence to the ancestor variables J a µ and K 0 .
These so-called bleached variables are found to be [4] where the action S appearing in the definition of K 0 is given in Eq. (2.13).
In terms of bleached variables, the one-loop version of the LFE can be cast in the form [9] and one can prove that the complete solution of the one-loop LFE is given by where A (1) is the 1PI functional of the ancestor amplitudes (2.22). Thus one finds that at this order, all the dependence on the pion fields is enclosed in the bleached variables.
This prescription allows one to write down all the descendant amplitudes depending on the pion fields in terms of the ancestor amplitudes; all one has to do is to expand the bleached variables up to the relevant order in the pion fields. A direct computation using the definition (2.23) shows that K 0 starts with two φ's, with For the variable j a µ one has instead the result (2.28) Thus, for example, one has for the two-and three-point pion sector (2.29) At higher orders in the loop expansion the bilinearity of the LFE (2.20) implies that there is an explicit dependence on the pion fields, governed by the equation [9] Notice that the above form of the LFE holds provided that Γ is written as a functional of the variables K 0 , j aµ .
The general solution of the LFE becomes then [9] Γ where G is the functional solving Eq. (2.30); as such it is uniquely fixed by the ancestor amplitudes, depends explicitly on φ a , and, finally, vanishes at φ a = 0. The existence of G can be proven by exploiting cohomological tools [9].

F. Renormalization
Summarizing, the combination of the LFE and the WPC, expressed in Eqs. (2.20) and (2.21) respectively, allows one to express the infinite number of divergent amplitudes involving the pion fields (descendant amplitudes) in terms of a finite number of ancestor amplitudes involving the connection J a µ and the source of the non-linear constraint (2.4) K 0 . It turns out that it is also possible to renormalize the theory in a symmetric fashion, that is, in a way that preserves the LFE [4,6].
Consider first the one-loop ancestor amplitudes. Taking into account Lorentz and global SU(2) R invariance, the list of UV divergent amplitudes reduces to the following eight (integrated) local monomials For example, for the one-loop two-point function one finds [4] Γ (1) where m ≡ m d=4 , and the dots indicate finite (µ-dependent) pieces. Then, we can dispose of this divergence by requiring that the monomials M 1,2 enters the counterterm action in the combination (2.34) Carrying out this procedure for all of the one-loop divergent ancestor amplitudes, one can fix the one-loop counterterm action S ct = i ρ B i M i thus rendering finite the NLσM at this level in the loop expansion.
The one-loop counterterms for the descendant amplitudes are then generated by expressing the monomials M i appearing in S ct in terms of the bleached variables, whence giving (2.35) in which all the covariant derivatives are defined with respect to the flat connection: It turns out that there is no one-loop counterterm associated with I 0 , as the theory is massless [4]. Since at φ a = 0 one has the normalization conditions given in the second line of Eq. (2.23), clearly we recover the counterterms introduced for the one-loop ancestor amplitudes; in addition, however, the above invariants generate the correct one-loop counterterms for all pion amplitudes, solving completely the hierarchy imposed by the LFE at this order.
At higher orders, say n > 1, the situation is slightly more complicated. The bilinear term in the LFE results in the term (2.30); however, this can only give rise to the mixing of lower order counterterms and therefore does not lead to new ones. As a consequence, this term will not appear in the evaluation of the n th -order counterterms for the ancestor amplitudes. Thus the symmetric subtraction procedure at order n is the following. One starts by computing the divergent part of the ancestor amplitudes that are superficially divergent according to the WPC condition (2.21). This will then fix the coefficients of the local monomials M i appearing at this order. Then one converts these monomials into the invariants I i by writing them in terms of the bleached variables. This will then give rise to the full set of counterterms required to make the theory finite at order n in the loop expansion.

III. A RENORMALIZATION GROUP EQUATION FOR WPC RENORMALIZ-ABLE THEORIES
From the discussion of the previous section we know that the LFE holds true for the effective action Γ of the theory, which comprises the tree-level Feynman rules plus counterterms. In addition, in the case of zero pion fields, the full bare effective action can be decomposed on a basis of integrated local monomials involving only the variables K 0 and J a µ and their derivatives In the above equation the sum spans all possible (infinite) local monomials, compatible with Lorentz invariance. By expressing the bare parameters ρ B i in terms of the renormalized ones ρ i and of the scale µ of the radiative corrections, one gets the effective action Γ, which yields a finite theory. Thus, from the point of view of Weinberg and Gomis renormalizability, the WPC selects which coefficients ρ B i must be zero in the tree-level approximation, and prescribes the loop order at which the counterterms of a given local monomial starts to appear (or, in other words, at which order a particular coefficient ρ B i = 0). Knowledge of the rhs of Eq. (3.1) completely fixes (through the LFE) the dependence on the pion fields of the complete bare action Γ B . Then, one can reabsorb the dependence on the scale µ of the radiative corrections into the renormalized parameters ρ i by expressing the bare parameters in terms of the renormalized ones: This is always possible, since due to the linearity of Γ B 0 on the bare parameters, one can reabsorb the divergences associated with the monomial M i by redefining the bare parameter ρ B i . In addition, Eq. (3.2) entails that the same result holds for the full vertex functional Γ, namely Next, by differentiating Eq. (3.3) with respect to the scale µ, we get the following RG which holds in full generality for all 1-PI Green's functions (including those involving an explicit dependence on the pion legs) 3 .
The very important question we are addressing in this paper is whether the WPC con- in the bare parameters ρ B i , and therefore a change in the scale µ, affecting the n th order action, can be accommodated by a change of the finite part of the n th order (ancestor amplitudes) counterterms, so that it appears like that the RG equation does not mix up the WPC hierarchy 5 . However this is not sufficient to prove compatibility with the WPC, as 3 In the Chiral Lagrangian approach (momentum expansion) a RG equation has been derived in [10]. 4 This is definitely not the case for the complete effective action Γ, as in this case the functional G appearing in the complete LFE solution (2.31) contains a product of lower order terms and therefore it has a complicated dependence on the ρ i . 5 At the level of amplitudes with the explicit dependence on the pion legs things are in general much more complicated and one has to resort to the LFE in order to fix them in a way compatible with the symmetry of the theory. the redefinition of the renormalized parameters ρ i still spans in principle at a given order n all (infinite) integrated local monomials corresponding to divergent ancestor amplitudes and compatible with Lorentz and SU(2) R global symmetry (under whichJ aµ is in the adjoint and K 0 is a singlet) on the space of ancestor variables. What we need to prove is that i spans only those monomials that are required by the WPC at order n and nothing else.

A. One-loop analysis
To see where the problem resides and what need to be proven, let us consider the oneloop case n = 1. As already said, in this case there are five divergent ancestor amplitudes: K 0 K 0 , and Γ (1) K 0 J J . The one-loop topologies possibly contributing to these amplitudes are shown in Fig. 1 through Fig. 4.
The WPC-compatible tree-level couplings used to construct the ancestor amplitudes are the ones coming from the action (2.19). On the other hand, the effective action Γ (1) will contain also the eight monomials (2.32) with a coefficient 6 λ (1) i , i = 0 (as already noticed, the coefficient ρ (1) 0 is zero at one loop level; in addition, notice that M 0 is allowed by the WPC, so its coefficient is not of the λ-type), fixed by the divergent part of the corresponding ancestor amplitude. When trading the K 0 and J a µ variables for the bleached ones, the monomials M i will become the invariants I i of Eq. (2.35); expanding then the bleached variables in terms of the pion fields as in Eqs. (2.26) and (2.28), will generate new vertices with pions and external sourcesJ µ and K 0 . Contrary to the tree-level vertices however, these so-called λ-vertices (see also the definition given in Sect. IV) violate the WPC and are proportional to a parameter λ (1) i , which was zero at tree-level (from which the name). A Feynman graph constructed from this type of vertices might in principle contribute to the RG equation (3.4) due to the derivative term in ρ i , evaluated on the λ-type coefficients.
To understand how this can possibly happen, recall that for writing down the RG equation one writes down all possible tree-level couplings compatible with the LFE which are of two types: the coefficient m d , which is also compatible with the WPC, and all the λ-vertices. 6 We reserve the notation λ i (respectively λ B i ) for the renormalized (respectively bare) parameters that are bound to be zero at tree-level due to the WPC. This is to be contrasted with ρ i (respectively ρ B i ) which denote the renormalized (respectively bare) coefficients of all monomials (that is, including those that are non-zero according to the WPC).  According to our previous discussion, one of the vertices appearing there must come from the tree-level Feynman rules and it is given by Tr while the remaining vertex is a λ-vertex. Thus, we expand the first five invariants   Fig. 2 shows the possible topologies contributing to this amplitude. As before the topology (t 31 ) vanishes; for graphs of the type (t 32 ) one has instead that the vertex with two external sources J and two pion fields has to be a λ-vertex since this kind of vertices is not present in the tree-level Feynman rules. Therefore, we now expand Notice that the above vertices contain at least one derivative acting on a pion field. Therefore, the UV degree of divergence of the (t 32 ) graphs is at least 2 (the power counting is the same as that of graphs (t 22 ) in Fig. 1 apart for the fact that in this case the λ-vertex contains at least one derivative). Thus, D(t 32 ) > D(Γ J J J ) so that, as in the previous case, this kind of graphs cannot contribute to one-loop counterterms.
Finally, the UV degree of divergence of (t 33 ) graphs must be at least 2, since they all have one loop, three bosonic propagators, one derivative from every tree-level vertex and at least two derivatives from the λ-vertex. So, by the same token, also these diagrams do not appear in the one-loop RG equation for Γ   Fig. 3. Once again the massless seagull graphs (t 41 ) vanish. Next, in order to compute the UV degree of divergence of the topology (t 42 ), we need to expand the one-loop invariants in powers of φ keeping only terms of the form J J Jφφ. One finds This completes the analysis of the J-sector at the one-loop level and shows that for the two-, three-and four-point functions of the external source J a change in the scale µ is compensated by a change of the finite parts of genuinely one-loop invariants. Let us conclude, by observing that even though we have taken into account only λ-vertices stemming from the one-loop invariants (2.35), the argument is valid also, a fortiori, for monomials that appear as counterterms at higher loops, for the latter will contain either more derivatives or more bleached variables j (or both).

K 0 -sector
In the K 0 sector one needs to consider the invariants I 6 and I 7 of Eq. (2.35). However, recall that the tree-level dependence from K 0 is completely fixed by the coupling to the Finally, we show in Fig. 4 the topologies that contribute to the ancestor amplitude with one external scalar source K 0 and twoJ legs.
As discussed above, the λ-vertex in both topologies can never come from the K 0 source. Therefore for (t 51 ) graphs the λ-vertex contains at least one derivative; hence the UV degree of divergence of these graphs is at least 1. This is also the UV degree of divergence for the topologies (t 52 ). Therefore, we conclude that the one-loop invariant with one scalar source K 0 and two J, which does not contain derivatives, cannot receive contributions from these graphs.
Also in this sector the results obtained are valid for λ-vertices originating from higher order counterterms. This completes the analysis of the K 0 -sector and it allows us to conclude that at the one-loop level a change in the scale µ only requires a change in the finite parts of the WPC one-loop invariants to be compensated.

IV. A GENERAL THEOREM
We are now ready to tackle a general proof of the fact that the WPC is preserved by the RG-flow. The strategy followed for proving this will be different from the one adopted for illustrating the one-loop case, which obviously cannot be adapted to an all-order analysis.
Before dwelling on the detailed proof let us recall the precise definition of a λ-vertex and state in a precise form the theorem we would like to prove.

Definition.
A λ-vertex is an interaction vertex generated upon the expansion of a symmetric (i.e., fulfilling the LFE) local functional forbidden by the WPC in powers of the pion fields φ. It must contain at least one external J leg 7 .
One has then the following Theorem. There are no λ-vertex contributions to the RG-flow of a n th loop ancestor amplitude.
The proof is divided into several steps that we detail in the following five subsections.

A. Loop expansion of the RG equation
Let us consider a λ i parameter that, according to the WPC, is bound to be zero up to the order n: (4.1) and suppose that there exists an order m < n such that a contribution to the RG equation arises at that order from the λ i -parameter.
Such a contribution to the second term in the lhs of the RG equation (3.4) for the ancestor amplitudes is given by Notice that one obtains a local contribution, as it should be since m is the lowest order where λ i is assumed to contribute. If such an integer m < n existed, one would clearly mix up 7 Again we remind that the dependence of the tree-level effective action on K 0 is fixed by the nonlinear SU(2) symmetry and hence no λ-vertices originating from invariants involving K 0 need to be considered. the WPC counting: a change in the scale µ would be reflected in lower order contributions, associated with counterterms that cannot appear at that order according to the WPC.
Let us now prove that this is indeed not the case.

B. Topologies
The amplitudes that we need to consider have no external pion legs and display an insertion of a single λ-vertex V λ , all other vertices being generated by the WPC tree-level action (2.19). We will denote with V λ K 0 (respectively, V λ J ) the number of K 0 (respectively, J a µ ) legs attached to the vertex V λ . Finally r will denote the number of pion legs attached to V λ .
The relevant n th -loop topologies can be classified according to the number of petals P composing a daisy diagram centered on V λ , see or n = n (1) + r − 1 for just one petal.

C. P-amplitudes degree of divergence
Let us indicate with P (i) φ (respectively, P (i) J ) the number of φ (respectively, K 0 , J a µ ) legs attached to the 1PI amplitude building up the i-th petal P (i) . Then, one has the following Lemma. The degree of divergence of a P (i) -amplitude satisfies the WPC bound (2.21), i.e., To prove this, let us analyze the two possible contributions to a P (i) -amplitude, that is the one coming from the bleached variables substitution and the one from the G term of the general solution of the LFE. Let be an ancestor amplitude that, upon the substitution of the N K 0 and N J legs contributes to the P (i) -amplitude under scrutiny Now, according to their definition (2.23) and the corresponding expansions (2.26) and (2.28), one observes that each substitution of a J-leg can give at most one derivative acting on the pion fields, while in the case of a K 0 -leg one gets at most two derivatives.
Then one has However this does not exhaust all the possible ancestor amplitudes eventually contributing to the given P-amplitude. In fact, observe that single derivatives of the pion fields coming from the replacement K 0 → K 0 /m d get accompanied by an extra J a µ ; the same is true for zero derivative pion fields terms coming from the replacement J a µ → −j a µ . Thus one has the following two cases: (i ) If we start from an ancestor amplitude and in the replacement of a J we get no derivative and a J, the number of external J of the descendant amplitude stays invariant; (ii ) If we start from an ancestor amplitude and in the replacement of a K 0 -leg we get a derivative and a J, then the additional contribution to the degree of divergence is through the term that is we get exactly the contribution of the additional J-leg generated.
Thus we conclude that for the part of the petal amplitude that is generated through the bleached variables the bound (4.4) holds.
Next consider the contribution coming from the bilinear term of the LFE. We know that the bound is satisfied when n = 1, so that we can proceed by induction. We then consider the differential operator O = δ ℓ+s+t δ J a 1 µ 1 (x 1 ) · · · δ J a ℓ µ ℓ (x ℓ )δK 0 (y 1 ) · · · δK 0 (y s )δφ b 1 (z 1 ) · · · δφ bt (z t ) (4.10) and we apply it to the bilinear equation (2.30). This operation will give rise to a term contributing to an amplitude involving ℓ J-legs, s K 0 -legs and t + 1 φ-legs. As there are many possible ways of acting on the rhs of Eq. (2.30), we then denote by n (I) J and n (II) J the number of J derivatives acting on δΓ (j) δK 0 (x) and δΓ (n−j) δφ a (x) respectively, and similarly for n (I) K 0 and n (II) K 0 .
Using the induction hypothesis, the UV degree of divergence of the two amplitudes obtained in this way are Summing everything up one gets again the result that is the bound holds at order n.
This concludes the proof of the lemma.

D. One-petal amplitudes
Consider now the n th -order amplitude Γ corresponding to a λ-vertex, with r φ-legs, V λ J J-legs and V λ K 0 K 0 -legs. Assume also that there is a single petal P (1) (see diagram (a) in Fig. 5), so that r ≡ P (1) φ and n = n (1) + P The UV degree of divergence of this amplitude is where δ ≥ 0 is the number of derivatives present in the λ-vertex. Since n (1) = n − r + 1, by means of the lemma previously proved we can write (4.14) A bound on δ can be then obtained by observing that the degree above cannot be greater than the one of a n th -order ancestor amplitude with V λ J + P J J-legs and V λ K 0 + P (1) thus we get the inequality One has then the following cases (i ) If V λ K 0 = 1, then V λ J = δ = 0 (i.e., the λ-vertex has no derivatives). In this case the λ-vertex is of the type K 0 φ 0 , already present in Γ (0) and allowed by the WPC.
(ii ) If V λ K 0 = 0, then either V λ J = 1, in which case δ = 1 so that these are the couplings of the type JF coming from S, or V λ J = 2, in which case δ = 0 so that these are the couplings of the type J 2 coming again from S.
Thus we find that the only possible λ-vertices are the ones allowed by the WPC, which proves the theorem at the level of a single petal amplitude.

E. Daisy amplitudes
Let us now consider a full daisy graph composed by m P-amplitudes. For such an amplitude one has then where m is the number of petals attached to the λ-vertex V λ and δ ′ is the number of derivatives of the λ-vertex. The lemma tells us that for each one of the petal amplitudes In addition, the loop order of the P (m) amplitude is given by Then, after simple algebra, Eq. (4.17) yields At this point one obtains the same bound as before on δ ′ , and therefore the same conclusions hold.
This completes the proof of the RG-flow theorem.

V. STABILITY IN WPC RENORMALIZABLE THEORIES
The existence of a RG equation allows one to extend the notion of stability of the classical theory to the non-linearly realized models based on the WPC.
In fact, the WPC prescribes uniquely which coefficients ρ i are non-zero at tree-level (and therefore it defines the set of λ-invariants). In addition, the RG flow theorem just proven implies that the finite parts of the counterterms, needed to reabsorb a change in the µdependence, appear exactly at the order where the pole part of the corresponding ancestor amplitudes becomes non-zero according to the WPC.
Suppose now that one adds a free finite (µ-independent) part a (n) i at order n in the loop expansion, in a way to preserve the symmetries of the theory (and not violating the WPC).
As the loop order n of the finite coefficient a does not correspond to the topological loop order 8 , one can proceed as in the power-counting renormalizable theories, that is one rescales to obtain a finite physically equivalent theory (as = 1).
In the power-counting renormalizable case, the rescaling will give back a term already present at the tree-level: this is the well-known stability of the tree-level action against radiative corrections. For WPC renormalizable theories however, under the rescaling the addition of a finite free coefficient a Incidentally, notice that there is yet another reason why one cannot add the symmetric finite renormalization terms a (n) i just discussed. The WPC uniquely identifies the graphs in the expansion based on the topological loop number, thus selecting a particular Hopf algebra, as the latter is constructed as a dual of the enveloping algebra of the Lie algebra of the Feynman graphs associated to the theory under scrutiny [11,12]. On the other hand, it also guarantees that there exists a suitable exponential map on this Hopf algebra [13] which allows the removal of all the divergences. The addition of any a (n) i is equivalent to a change in the Hopf algebra of the model, as it would modify the set of 1-PI Feynman diagrams on which the Hopf algebra is constructed. This change destroys the compatibility between the WPC and the RG equation; therefore, the addition of such terms is not allowed.
The existence of a connection between the WPC preserving RG flow and the selected Hopf algebra clearly deserves further investigations.

VI. BEYOND THE STANDARD MODEL: THE WPC AS A MODEL-BUILDING PRINCIPLE
If one promotes the classical source J a µ to the status of a dynamical field, the NLσM action gives rise to the Stückelberg mass term. By formulating a non-linearly realized SU(2) Yang-Mills theory in the LFE framework [14][15][16] (with the pseudo-Goldstone fields taking over to the role of the pion fields) and imposing the WPC, one arrives to a somewhat surprising result.
Specifically, notwithstanding the fact that if the local gauge symmetry is realized nonlinearly the Yang-Mills action is not singled out on the basis of gauge invariance 9 , it turns out that if the WPC condition is satisfied, then the only admissible solution of the tree-level LFE is the Yang-Mills action plus the Stückelberg mass term: where G a µν is the field strength of the gauge field A a µ and F a µ the flat connection. An important consequence of this fact is that one can formulate a non-linearly realized theory based on the gauge group SU(2)×U(1) [17][18][19]. One combines the SU(2) gauge fields A µ = A a µ τa 2 and the U(1) Y gauge field B µ into the bleached variable g and g ′ are the SU(2) and U(1) coupling constants respectively. The bleached counterparts of the A, Z and W ± fields are given by where s W (c W ) is the sine (cosine) of the Weinberg angle the tangent of which is given by the ratio g ′ /g. The action of this non-linear version of the Eletroweak Theory is highly constrained if one requires the WPC to be satisfied: in this case the self-couplings of gauge bosons as well as the couplings between gauge bosons and fermions are the same as the conventional SM ones. However, the following combination of two independent mass invariants arise, and the Weinberg relation is broken This is a peculiar feature of nonlinearly realized electroweak theories [20].
The inclusion of physical scalar resonances in the non-linearly realized Electroweak Theory, while respecting the WPC, yields a definite beyond the Standard Model (bSM) scenario.
Indeed it turns out that it is impossible to add a scalar singlet without breaking the WPC condition [7]. The minimal solution requires a SU(2) scalar physical doublet, leading to a CP-even physical field (to be identified with the recently discovered scalar resonance at 125.6 GeV) and three additional heavier physical states, one neutral CP-odd and and two charged ones. Notice that this is a rather peculiar signature, since in Two Higgs-Doublet Models and the Minimal Supersymmetric Standard Model the number of physical scalar resonances is five (see, e.g., [21]).
While some preliminary phenomenology issues of this model have been addressed in [7], a full analysis and comparison with the experimental data can be carried out in a satisfactory and systematic way only in the presence of an RG-flow, as the possibility of running the scale µ in a mathematically consistent way would allow one to obtain physical predictions of the same observables applicable in different energy regimes.
The extension of the analysis carried out for the NLσM to the non-linearly realized electroweak theory with scalar resonances (NLSM for short) requires some care.
The mass generation mechanism for gauge bosons cannot be entirely of the Stückelberg type, for in this case the decays of the Higgs scalar h → W W * and h → ZZ * would be loopinduced and thus phenomenologically unacceptably small. Therefore the problem arises to assess whether a fraction of the mass is generated by the Stückelberg mechanism, in addition to the contribution associated with the linearly realized spontaneous symmetry breakingà la Higgs.
Moreover, since current LHC data are in very good agreement with the Standard Model [22,23], one can assume in a first approximation custodial symmetry and set κ = 0 in Eq. (6.5). The W and Z masses are therefore where G = g 2 + g ′ 2 and A is a parameter of mass dimension squared describing the fraction of the mass generated by the Stückelberg mechanism.
The Lagrange multiplier formulation of the non-linear constraint is particularly suited for studying the small A limit, which is the phenomenologically relevant regime since bSM effects are known to be small [22,23]. This is because in this formulation it is easy to derive the dominant contribution in the small A expansion of the 1-PI amplitudes (without the need of resummations) [8].
The Lagrange multiplier (BRST-invariant) implementation of the non-linear constraint is realized by introducing a pair of BRST variables B, c such that sB = c, sc = 0. B is the Lagrange multiplier field, c is the associated ghost (that is free). On the other hand, the BRST variation of the antighostc yields the (invariant) non-linear constraint: There is no WPC in the sector spanned byc and its BRST variation sc [8]. However, since they form a BRST doublet 11 [24], it turns out that they can only modify the BRST-exact sector of the theory, which, like the gauge-fixing, is not physical.
On the other hand, in the gauge-invariant (BRST-closed) sector of the theory, the WPC holds. In this sector the proof of the compatibility between the WPC and the RG flow given (like, for instance, the ghost fields). It turns out that the bleaching procedure for a variable that transforms in the fundamental or in the adjoint representation of the gauge group does not involve derivatives. To be sure, one has the expressions f = Ω † f (respectively, X = Ω † XΩ) for the bleached variable associated to a field that transforms in the fundamental (respectively, adjoint) representation of the gauge group. Therefore the substitution of these fields with their bleached counterparts does not modify the degree of divergence of the amplitude considered.
Moreover one should also note that in a massive theory the WPC provides only an upper bound on the degree of divergence of an ancestor amplitudes, with the bound being saturated only in the massless case. Yet the upper bound is sufficient in order to establish the validity of the RG equation, as can be seen from the analysis of Sect. IV.
One therefore sees that the RG equation is again compatible with the WPC in the physically relevant gauge-invariant sector of the theory.

VII. CONCLUSIONS AND OUTLOOK
The existence of a RG equation for the NLσM and the NLSM, compatible with the WPC, shows that the sliding of the scale µ on physical amplitudes can be reabsorbed by suitable finite counterterms, arising at the loop order prescribed by the WPC itself. As a result, the running with energy of physical observables becomes a consistent procedure also within non-linearly realized theories based on the LFE and for which the WPC holds.
Moreover one can formulate the notion of weak stability, in close analogy with the powercounting renormalizable case: the inclusion of free finite counterterms at higher order in the loop expansion alters the Hopf algebra of the model and moreover generates λ-vertices, thus mixing up the order of the counterterms needed to recover the effect of a change in the scale µ of the radiative corrections. This destroys the compatibility between the WPC and the RG equation.
It is a rather remarkable fact that in the NLSM the RG equation can be written for a theory with a finite number of parameters dictated by the WPC. Consequently, since the RG equation allows to run physical observables with energy, the parameterization of bSM physics embedded in the NLSM can be tested on a wide range of energy, from the GeV to the TeV scale.
This provides a consistent theoretical framework for the study of the non-linear Stückelberg-like symmetry breaking contribution (and their bSM implications) to the fermion and gauge bosons mass generation mechanism, that will be one of the main goals of the next LHC run.