Abstract

Perturbative canonical quantum gravity is considered, when coupled to a renormalizable model for matter fields. It is proposed that the functional integral over the dilaton field should be disentangled from the other integrations over the metric fields. This should generate a conformally invariant theory as an intermediate result, where the conformal anomalies must be constrained to cancel out. When the residual metric is treated as a background, and if this background is taken to be flat, this leads to a novel constraint: in combination with the dilaton contributions, the matter lagrangian should have a vanishing beta function. The zeros of this beta function are isolated points in the landscape of quantum field theories, and so we arrive at a denumerable, or perhaps even finite, set of quantum theories for matter, where not only the coupling constants, but also the masses and the cosmological constant are all fixed, and computable, in terms of the Planck units.

## 1 Introduction

In a previous paper[1], it was argued that the functional integral in canonical quantum gravity,

(1.1) |

where is the Einstein-Hilbert action and is the action of the matter fields, here abbreviated as , should be considered to be taken in two steps:

(1.2) |

and the integral over the dilaton field should be taken together with the integrations over the matter fields .

Rewriting the Einstein-Hilbert action (including a possible cosmological constant) in terms of and , one finds, in four dimensions,

(1.3) |

where , and is the Ricci scalar associated to . It is convenient to split the lagrangian of the matter fields into conformally invariant kinetic parts, mass terms, and interaction terms. In a simplified notation (later we will be more precise), one has

(1.4) | |||||

(1.5) | |||||

(1.6) |

Here, ; the kinetic term for the Dirac fields is shorthand for the corresponding expression using a vierbein field for the metric with its associated connection field. stands short for the scalar masses and for the Dirac masses. The term could be removed by a field redefinition, but is kept here for convenience, making the scalar lagrangian conformally invariant.

Interaction between the matter fields must now be written in the form

(1.7) |

where the Yukawa couplings could be matrices in the indices labeling the fermion species, and and could be 4- and 3-index tensors in the scalar field indices.

If now we rescale the field:

(1.8) |

we notice two things. First, the action for the field is now nearly identical to the kinetic term for the field, both having the same conformal dimension:

(1.9) |

and secondly, the mass terms turn into conformally invariant quartic coupling terms between matter fields and the field:

(1.10) |

where has the dimension of an inverse mass. Also the scalar 3-field coupling, which originally was not conformally invariant, now turns into a conformally invariant 4 field coupling (see Eq. (1.7)):

(1.11) |

and a new quartic interaction term for the field is generated by the cosmological constant:

(1.12) |

It is important to observe that the lagrangian (1.9) for the dilaton field has an overall sign opposite to that of ordinary scalar fields . For any other field theory, this would be disastrous because it would violate causality. Here however, the unconventional sign is a necessary consequence of the canonical structure of the theory. Since it is an overall sign, it has no net effect on the Feynman rules; this we will exploit by rotating the field in the complex plane:

(1.13) |

so that the new field will be indistinguishable from other scalar fields, with one important exception: the conformal interaction terms from the original mass terms, Eq. (1.10), as well as the interaction from the original 3-field interaction, Eq. (1.11), get unconventional factors -1 or .

Another important thing to observe is that there is no term at all in the lagrangian that could serve as a kinetic term
for the field^{1}^{1}1Of course, this term reappears if one substitutes , where
is a small oscillating field, by expanding .. Any such kinetic terms should arise solely from
higher order effects due to the interactions with the matter fields. Naturally, since the entire lagrangian now is
conformally invariant, we should expect the effective lagrangian for to be conformally invariant as
well. It is emphasized in Ref. [1] however that, as is well known [2], this conformal invariance
is mutilated by conformal anomalies.

In this paper, we now propose to postpone any attempts to describe the functional integrals over . Not
only do we have anomalies there, but also there are difficulties with unitarity and Landau ghosts.^{2}^{2}2In
Ref. [3],[4] claims are made that unitarity can be restored. This however requires the integration
contours to be rotated in the complex plane such that becomes complex. Such approaches are
interesting and may well serve as good starting points for generating promising theories, but they will not be pursued
here. As was explained in Ref.[1], the effective interactions with matter and dilaton fields generate
an action for that largely coincides with the familiar conformally invariant action obtained from the
square of the Weyl curvature, but with an infinite numerical coefficient, which would have to be renormalized.
The difficulties associated to that are sufficient reason for us now to postpone this sector of the theory entirely.

At first sight, one may well find logical objections to such a procedure: why not also first integrate over before examining whether the amplitudes obtained obey conformal constraints? We argue however that the integration is very different from the rest; determines the location of the local light cones, so that it determines the causal relationships between points in space-time. It may well be that quantum interference of states with light cones at different places will require treatments that differ in essential ways from the standard functional integral.

In any case, it is worth investigating what happens if we follow this procedure. The implications are quite remarkable, as we will show.

## 2 Renormalization

Let us assume that the matter fields consist of Yang-Mills fields , Dirac fields
and scalar fields , the latter three sets being in some (reducible or irreducible,
chiral or non chiral) representation of the local Yang-Mills gauge group. For brevity, we will write complex
scalar fields as pairs of real fields, and if Weyl or Majorana fermions occur, the Dirac fields can be
replaced by pairs of these.^{3}^{3}3A single Weyl or Majorana fermion then counts as half a Dirac field. Let us
rewrite the lagrangian for matter interacting with gravity more precisely than in the previous, introductory
section:

(2.1) | |||||

where is the (non Abelian) Yang-Mills curvature, and and are covariant derivatives containing the Yang-Mills fields; and also contain the vierbein fields and connection fields associated to ; the Yukawa couplings , and fermion mass terms are matrices in terms of the fermion indices. The scalar self interactions, and must be a third and fourth degree polynomials in the fields :

(2.2) | |||||

(2.3) |

In Eq. (2.1), like , also are mass matrices, in general. Furthermore, stands for . Of course, all terms in (2.1) must be fully invariant under the Yang-Mills gauge rotations. They must also be free of Adler Bell Jackiw anomalies.[6]

Now that the dilaton field has been included, the entire lagrangian has been made conformally invariant. It is so by construction, and no violations of conformal invariance should be expected. This invites us to consider the beta functions of the theory. Can we conclude at this point that the beta functions should all vanish?

Let us not be too hasty. In the standard canonical theory, matter fields and their ineractions are renormalized. Let us consider dimensional renormalization, and the associated anomalous behavior under scaling. In dimensions, where is infinitesimal, the scalar field dimensions are those of a mass raised to the power , so that the couplings have dimension . This means that in most of the terms in the lagrangian (2.1) the integral powers of will receive extra factors of the form or , which will then restore exact conformal invariance at all values for . If we follow standard procedures, we accept that is close to , so that singularities at or are not considered to be of any significance. Indeed, the limit may be seen to be the small-distance limit. This is the limit where gravity goes wrong anyway, so why bother?

However, now one could consider an extra condition on the theory. Let us assume that the causal structure, that is, the location of the light cones, is determined by , and that there exist dynamical laws for . This was seen to be a very useful starting point for a better understanding of black hole complementarity[5]. The laws determining the scale should be considered to be dynamical laws, and the canonical theory of gravity itself would support this: formally the functional integral over the fields is exactly the same as that for other scalar fields.

In view of the above, we do think it is worthwhile to pursue the idea that the field must be handled just as any other scalar component of the matter fields; but then, after renormalization, fractional powers, in the limit would lead to terms, and these must clearly be excluded. Renormalization must be done in such a way that no traces of logarithms are left behind. Certainly then, a scale transformation, which should be identical to a transformation where the fields are scaled, should not be associated with anomalies. Implicitly, this also means that the region is now assumed to be regular. This is the small distance region, so that, indeed, our theory says something non-trivial about small distances. This is why our theory leads to new predictions that eventually should be testable. Predictions follow from the demand that all beta functions of the conformal “theory” (2.1) must vanish.

Note that one set of terms is absent in Eq. (2.1): the terms linear in and hence cubic in . This, of course, follows from the fact that, usually, no terms linear in the scalar fields are needed in the standard matter lagrangians; such terms can easily be removed by translations of the fields: for some constants . Thus, the classical lagrangian is stationary when the fields vanish: is a classical solution. In our present notation, this observation is equivalent to the observation that fields may be freely transformed into one another without modifying the physics. One such transformation is a rotation of one of the scalar fields, say , with the fields:

(2.4) |

where stands for the original shift of the field . The transformation is taken to be a hyperbolic rotation because the “kinetic term” in Eq. (2.1) has to be invariant.

In most cases, these transformations need not be considered since terms linear in will in general not be gauge invariant.

Notice also that the Yang-Mills fields are not directly coupled to the field. In the “classical limit”, , we can see why this is so. Since not , but is real, the invariant quantity is . Rotating fields with fields would therefore be a non-compact transformation, and Yang Mills theories with non-compact Lie groups usually do not work. There is food for thought here, but as yet we will not pursue that.

When constructing a complete theory, the next step should be to consider just any configuration of , formulate the renormalized theory in this background, and finally consider functional integrals over . Unfortunately, this is still too difficult. In a non-trivial background, there will be anomalies depending on the derivatives of ; there are divergences[1] proportional to the Weyl curvature squared, which can be seen to correspond to the field combination , and subtracting those leads to new conformal anomalies[2]. It was suggested in Ref. [1] to keep the conformal infinity, which would turn gravitons into “classical” particles, or more precisely, particles that cannot interfere quantum mechanically, but whether this can be held up as a theory remains to be seen. To avoid further complications, we now decided to look at the case when , or, space-time is basically flat, though we just keep the field .

## 3 The functions

Thus, we return to a theory to which all known quantum field theory procedures can be applied, the only new thing being the presence of an extra, gauge neutral, spinless field , and the perfect local scale invariance of the theory.

We arrived at the lagrangian (2.1), and we wish to impose on it the condition that all its beta functions vanish, since conformal invariance has to be kept. As the theory is renormalizable, the number of beta functions is always exactly equal to the number of freely adjustable parameters. In other words: we have exactly as many equations as there are freely adjustable unknown variables, so that all coupling constants, all mass terms and also the cosmological constant, should be completely fixed by the equations . They are at the stationary points. Masses come in the combination and the cosmological constant in the combination , so all dimensionful parameters of the theory will be fixed in terms of Planck units.

In principle, there is no reason to expect any of these fixed points to be very close yet not on any of the axes, so neither masses nor the cosmological constant can be expected to be unnaturally small, at this stage of the theory. In other words, as yet no resolution of the hierarchy problem is in sight: why are many of the physical mass terms 40 orders of magnitude smaller than the Planck mass, and the cosmological constant more than 120 orders of magnitude? We have no answer to that in this paper, but we will show that the equations are quite complex, and exotic solutions cannot be excluded.

The existence of infinitely many solutions cannot be excluded. This is because one can still adjust the
composition and the rank(s) of the Yang-Mills gauge group, as well as one’s choice of the scalar and (chiral)
spinor representations.^{4}^{4}4which of course must be free of Adler Bell Jackiw
anomalies[6][7]. These form infinite, discrete sets. However, many choices turn out not to
have any non trivial, physically acceptable fixed point at all: the interaction potential terms
must be real and properly bounded, for instance. Searches for fixed points then automatically lead to
vanishing values of some or more of the coupling parameters, which would mean that the symmetries and
representations have not been chosen correctly.

Every advantage has its disadvantage. Since all parameters of the theory will be fixed, we cannot apply perturbation theory. However, we can make judicious choices of the scalar and spinor representations in such a way that the existence of a fixed point for the gauge coupling to these fields can be made virtually certain. The function for gauge theories with fermions and complex scalars in the elementary representation is

(3.1) | |||||

(3.2) | |||||

(3.3) |

Choosing one scalar extra, or one missing, we can have as small as , while a quick inspection in the literature[8][9] reveals that, in that case, may still have either sign:

(3.4) |

depending on further details, such as the ratio of fermions and scalars, the presence of other representations, and the values of the Yukawa couplings. Choosing the sign of opposite to that of , one then expects that a fixed point can be found at

(3.5) |

This, we presume, is close enough to zero that the the following procedure may be assumed to be reliable. Let there be physical constants, the one being the gauge coupling , which is determined by the above equation (3.5). If we take all other coupling parameters to be of order or , then the beta function equations are reliably given by the one-loop expressions only, which we will give below. Now these are equations for the remaining coupling parameters, and they are now inhomogeneous equations, since the one coupling, , is already fixed. All we have to do now, is find physically acceptable solutions. We already saw that non-Abelian Yang-Mills fields are mandatory; we will quickly discover that, besides the fields, both fermions and other scalar matter fields are indispensable to find any non-trivial solutions.

One trivial, yet interesting solution must be mentioned: super Yang-Mills. We take its lagrangian, and add to that the field while postulating that this field does not couple to the matter fields at all. Then indeed all functions vanish[10]. However, since the field is not allowed to couple, the physical masses are all strictly zero, which disqualifies the theory physically. Note, however, that also the cosmological constant is rigorously zero. Perhaps the procedure described above can be applied by modifying slightly the representations in this theory, so that a solution with masses close to zero, and in particular a cosmological constant close to, but not exactly zero, emerges.

The one loop functions are generated by an algebra[11], in which one simply has to plug the Casimir operators of the Yang Mills Lie group, the types of the representations, the quartic scalar couplings and the fermionic couplings. If we take the scalar fields and together as , the generic lagrangian can be written as

(3.6) |

where and are in general in reducible representations of the gauge group, is the gauge covariant derivative, is a gauge-invariant quartic scalar potential, and and are matrices in terms of the fermion flavor indices. Everything must be gauge invariant and the theory must be anomaly free[6][7].

The covariant derivatives contain the hermitean representation matrices , and :

(3.7) | |||||

(3.8) |

The gauge coupling constant(s) are assumed to be included in these matrices and . The operators and are projection operators for the left- and right handed chiral fermions.

The group structure constants are also assumed to include a factor , and they are defined by

(3.9) |

Casimir operators and will be defined as

(3.10) |

All beta functions are given by writing down how the entire lagrangian (3.6) runs as a function of the scale [11]:

(3.11) | |||||

(3.12) |

Here,

(3.13) | |||||

where .

It is convenient to define the complex matrices as

(3.14) |

Then,

(3.15) | |||||

If now we write the collection of scalars as , taking due notice of the factors in all terms odd in , we can apply this algebra to compute all functions of the lagrangian (2.1).

The values of the various functions depend strongly on the choice of the gauge group, the representations, the scalar potential function and the algebra for the Yukawa terms, and there are very many possible choices to make. However, the signs of most terms are fixed by the algebra (3.12)—(3.15). By observing these signs, we can determine which are the most essential algebraic constraints they impose on possible solutions. As we will see, they are severely restrictive.

## 4 Adding the dilaton field to the algebra for the functions.

Consider the dilaton field added to the lagrangian, as in Eq. (2.1). This requires extending the indices in the lagrangian (3.6) to include a value referring to the field. The unusual thing is now that the terms odd in are purely imaginary, while all terms in Eq. (2.1) are of dimension 4. Let us split off this special component. In Eq. (3.6), we then write

(4.1) | |||||

(4.2) |

Here, is the fermionic mass matrix, are matrices representing the scalar Yukawa couplings, and are the pseudoscalar Yukawa coupling matrices. is the scalar mass matrix, which is allowed to have negative eigenvalues, so we allow the Higgs mechanism to take place. We henceforth choose modified Planck units by setting .

Note, that in the Standard Model, there is no gauge-invariant Dirac mass matrix and no gauge-invariant cubic scalar interaction, so there and are zero, but we will need to be more general.

We write , and now . The algebra (3.12) — (3.15) is then found to become

where , etc, and we apply summation convention: double indices are summed over starting with 1, except the index in the scalar mass matrix , which is only summed over if it occurs twice elsewhere as well, or if this is explicitly indicated. is the expression that we already had, in Eq. (3.13).

The beta coefficients for the Yukawa couplings follow from adding the index 0 to Eq. (3.15):

(4.4) | |||||

(4.5) | |||||

where stands for the standard function for the corresponding dimension 4 interaction terms.

Before demanding that the functions all vanish, we observe that we can still allow for an infinitesimal field transformation of the form (2.4) in the original lagrangian. This adds to the counter terms:

(4.6) |

and a similar rotation in the Yukawa couplings. This can be used to eliminate the term (4 e) by adjusting ; it corresponds to a field shift in the non-gravitational case. In most cases, however, such as in the Standard Model, the terms in (4 e) are forced to vanish anyhow due to gauge invariance. In a similar way, an infinitesimal chiral rotation among the fermions can be used to eliminate the last term in Eq. (4.5).

Thus, after term (e) has been made to vanish by hand, the demand that all functions vanish, for all values of , applies in particular to the terms in Eq. (4 a — d) and (f), and to Eqs. (4.4) and (4.5).

We have already assumed that the non-Abelian Yang Mills field coupling(s) have a small but non-vanishing fixed point. Through the effects of the group matrices and , the coupling(s) determine the values of the other parameters, by as many coupled non-linear equations as there are unknowns. It follows that, in this theory, we must have non-Abelian Yang-Mills fields. In contrast, Abelian components are not allowed since those do not have fixed points close to the origin.

Next, let us consider the requirement that the term (f) in Eq. (4) vanishes:

(4.7) |

The r.h.s. of this equation resembles a supertrace. Since its sign must be positive, we read off right away that there must be fermions. If furthermore we like to have a very small or vanishing cosmological constant , we clearly need that the sum of the fourth power of the Dirac fields (approximately) equals the sum of the fourth powers of the masses of the real scalar particles divided by 4.

Can we do without the scalar fields ? Eq. (4.7) would have a solution, although the cosmological constant would come out fairly large. However, now there is only one more equation to consider: Eq. (4.5), with all Yukawa couplings and replaced by 0. That gives:

(4.8) |

Whenever has a real, non vanishing eigenvalue, this would imply that the trace of is negative, an impossible demand. Therefore, our theory also must have scalars , besides the dilaton field .

It appears that in today’s particle models not only the cosmological constant but also the mass terms are quite small, in the units chosen, which are our modified Planck units. Also, if there is a triple scalar coupling, , it appears to be small as well. This is the hierarchy problem, for which we cannot offer any solution other than suggesting that we may have to choose a very complex group structure — as in the landscape scenarios often proposed in superstring theories. Perhaps, the small numbers in our present theory are all related.

If the masses are indeed all small, then the only large terms in our equations are the ones that say how the coupling constants and masses run with scale. Our theory suggests that they might stop running at some scale; in any case, a light Higgs particle indeed follows from the demand that the Higgs self coupling is near an UV fixed point.

The author did not (yet) succeed in finding a physically interesting prototype model with a non-trivial fixed point; this is a very complex, but interesting technical problem. Let us briefly set out a strategy.

We search for a solution where all mass terms, and of course also the cosmological constant, are small. Start with a theory that has nearly, but not quite, a set of functions that vanish at one-loop. Assume that it has a fixed point near the origin. It is known that Eqs. (3.12)—(3.15) allow this. For simplicity, let us assume that there are no triple scalar couplings, , and no gauge invariant scalars, so that the terms (4 c) and (e) are forbidden by gauge invariance. If we deviate slightly from the fixed point, term (b) dictates that we must be at a point where the function for the scalar self couplings is negative. This gives us a first guess for , the absolute values of the scalar masses, but not the signs of .

Knowing the approximate values of the Yukawa couplings and allows us to fix the Dirac mass matrix using Eq. (4.5). Since we want these masses to be small, we must assume that the pseudoscalar Yukawa couplings largely cancel against the scalar ones in this equation. Then, Eq. (4.4) can be obeyed by moving slightly away from the original fixed point in that direction.

The only remaining equation is then the vanishing of term (d), the running of the scalar mass-squared terms. Knowing and , and assuming that is very small, then gives us an equation for . Notice that its sign was still free, so that there is some freedom here. Various further attempts to refine this procedure may well lead to interesting models with fixed points. We do note that, apparently, relatively large pseudoscalar Yukawa couplings are wanted. Also, we are talking about primary, gauge invariant Dirac masses, which have to be there due to the term (4 f), while they do not occur in the presently known Standard Model.

## 5 Discussion.

Our theory derives constraints from the fact that matter fields interact with gravity. The basic assumption could be called a new version of relativity: the scalar matter fields should not be fundamentally different from the dilaton field . Since there are no singular interactions when a scalar field tends to zero, there is no reason to expect any singularity when tends to zero at some point in space-time. Standard gravity theory does have singularities there: this domain refers to the short distance behavior of gravity, which is usually considered to be “not understood”. What if the short distance behavior of gravity and matter fields is determined by simply demanding the absence of a singularity? Matter and dilaton then join smoothly together in a perfectly conformally invariant theory. This, however, only works if all functions of this theory vanish: its coupling parameters must be at a fixed point. There are only discrete sets of such fixed points. Many theories have no fixed point at all in the domain where physical constants are real and positive — that is, stable. Searching for non trivial fixed points will be an interesting and important exercise.

Indeed, all physical parameters, including the cosmological constant, will be fixed and calculable in terms of the Planck units. This may be a blessing and a curse at the same time. It is a blessing because this removes all dimensionless freely adjustable real numbers from our theory; everything is calculable, using techniques known today; there is a strictly discrete set of models, where the only freedom we have is the choice of gauge groups and representations. It is difficult to tell how many solutions there are; the number is probably infinite.

This result is also a curse, because the values these numbers have in the real world is strange mix indeed: the range of the absolute values cover some 122 orders of magnitude:

(5.1) |

The question where these various hierarchies of very large, or small, numbers come from is a great mystery called the “hierarchy problem”. In our theory these hierarchies will be difficult to explain, but we do emphasize that the equations are highly complex, and possibly theories with large gauge groups and representations have the potential to generate such numbers.

Our theory is a “top-down” theory, meaning that it explains masses and couplings at or near the Planck domain. It will be difficult to formulate any firm predictions about physics at energies as low as the TeV domain. Perhaps we should expect large regions on a logarithmic scale with an apparently unnatural scaling behavior. There is in principle no supersymmetry, although the mathematics of supersymmetry will be very helpful for constructing the first non-tivial models.

What is missing furthermore is an acceptable description of the dynamics of the remaining parts of the metric field. In Ref. [1], it was suggested that this dynamics may be non quantum mechanical, although this does raise the question how can back react on the presence of quantum matter. Standard quantum mechanics possibly does not apply to because the notion of energy is absent in a conformal theory, and consequently the use of a hamiltonian may become problematic. A hamiltonian can only be defined after coordinates and conformal factor have been chosen, while this is something one might prefer not to do. The author believes that quantum mechanics itself will have to be carefully reformulated before we can really address this problem.

Our theory indeed is complex. We found that the presence of non-Abelian Yang-Mills fields, scalar fields and spinor fields is required, while gauge fields are forbidden (at least at weak coupling, since the function for the charges here is known to be positive). Because of this, one “prediction” stands out: there will be magnetic monopoles, although presumably their masses will be of the order of the Planck mass.

Finally, there is one other firm prediction: the constants of nature will indeed be truly constant. Attempts to experimentally observe variations in constants such as the finestructure constant or the proton electron mass ratio, with time, or position in distant galaxies, are predicted to yield negative results.

## Acknowledgements

The author benefited from discussions with C. Kounnas and C. Bachas at the ENS in Paris.

## References

- [1] G. ’t Hooft, Probing the small distance structure of canonical quantum gravity using the conformal group, arXiv:1009.0669v2 [gr-qc].
- [2] D.M. Capper and M.J. Duff, Conformal Anomalies and the Renormalizability Problem in Quantum Gravity, Phys. Lett. 53A 361 (1975); M.J. Duff, Twenty Years of the Weyl Anomaly, Talk givenat the Salamfest, ICTP, Trieste, March 1993, arXiv:hep-th/9308075.
- [3] P. D. Mannheim and D. Kazanas, Astrophys. J. 342, 635 (1989); D. Kazanas and P. D. Mannheim, Astrophys. J. Suppl. 76, 431 (1991); P. D. Mannheim, Prog. Part. Nucl. Phys. 56, 340 (2006), astro-ph/0505266; P. D. Mannheim, Intrinsically Quantum-Mechanical Gravity and the Cosmological Constant Problem, arXiv:1005.5108 [hep-th]; G.U. Varieschi, A Kinematical Approach to Conformal Cosmology, Gen. Rel. Grav. 42 929 (2010), arXiv:0809.4729.
- [4] C.M. Bender and P.D. Mannheim, No-ghost theorem for the fourth-order derivative Pais-Uhlenbeck oscillator model, Physical Review Letters 100, 110402 (2008), arXiv:0706.0207 [hep-th]; id., Exactly solvable PT-symmetric Hamiltonian having no Hermitian counterpart, Phys. Rev. D 78, 025022 (2008), arXiv:0804.4190 [hep-th].
- [5] G. ’t Hooft, Quantum Gravity without Space-time Singularities or Horizons, Erice School of Subnuclear Physics 2009, to be publ.; arXiv:0909.3426.
- [6] S.L. Adler, Phys. Rev. 177 (1969) 2426; J.S. Bell and R. Jackiw, Nuovo Cim. 60A (1969) 47.
- [7] S.L. Adler and W.A. Bardeen, Phys. Rev. 182 (1969) 1517; W.A. Bardeen, Phys. Rev. 184 (1969) 1848.
- [8] D.R.T. Jones, Nucl. Phys. B75 (1974) 531; W.E. Caswell, Phys. Rev. Lett. 33 (1974) 244; R. van Damme, Phys. Lett. 110B (1982) 239; 2-loop renormalization of the gauge coupling and the scalar potential for an arbitrary renormalizable field theory, Nucl. Phys. B227 (1983) 317.
- [9] A.G.M. Pickering, J.A. Gracey and D.R.T. Jones, Three loop gauge beta-function for the most general single gauge-coupling theory, arXiv:hep-ph/0104247.
- [10] L. Brink, J.H. Schwarz and J. Scherk, Supersymmetric Yang-Mills Theories, ERDA Research and development report, CALT-68-574 (1976); M. Sohnius Introducing supersymmetry, Physics Reports 128 (2-3) (1985) 39-204.
- [11] G. ’t Hooft, The birth of asymptotic freedom, Nucl. Phys. B254 (1985) 11; id., The Conceptual Basis of Quantum Field Theory, in Handbook of the Philosophy of Science. Philosophy of Physics, Vol. Eds. J. Butterfield and J. Earman, 2007 Elsevier B.V., p. 661.