Renormalizable and unitary Lorentz invariant model of quantum gravity

We analyze the R + R2 model of quantum gravity where terms quadratic in the curvature tensor are added to the General Relativity action. This model was recently proved to be a self-consistent quantum theory of gravitation, being both renormalizable and unitary. The model can be made practically indistinguishable from General Relativity at astrophysical and cosmological scales by the proper choice of parameters.


Introduction
Creation of quantum gravity still remains a prominent task of modern physics.
The problem is due to well known perturbative non-renormalizability of Einstein General Relativity. In the work [1] it was shown by direct calculations that at the one loop level General Relativity is renormalizable without matter fields but becomes unrenormalizable after inclusion of matter fields. Then also by explicit calculations it was demonstrated [2] that General Relativity is non-renormalizable at the two-loop level even without matter fields.
We mean here that the theory is perturbatively non-renormalizable. There are of course Quantum Field Theory examples of perturbatively non-renormalizable theories that lead to clear calculable predictions, such as for example the nonlinear sigma model above two dimensions. We will work within perturbation theory in the present paper and will not further consider non-perturbative aspects.
In [3] renormalizability of the R+R 2 -theory was proved. The proof used a specific covariant gauge for simplicity. For general gauges an assumption was made that ultraviolet divergences have the so called cohomological structure. This hypothesis was proved for a class of background gauges in the work [4]. Hence we consider renormalizability of R + R 2 gravity with four derivatives of the metric as well established.
We will also call this model briefly as quadratic gravity. But in the works [3,5] it was also stated that quadratic gravity is not physical because it violates unitarity or causality. So this model was commonly considered as unphysical.
Quite recenyly quadratic quantum gravity was proved to be in fact unitary [6]. Thus the R + R 2 model is a candidate for the quantum theory of gravitation.
In the present paper we discuss in detail the exact form of the Lagrangian of quadratic gravity, the questions of unitarity, stability of the vacuum state and the behavior of the model at astrophysical and cosmological scales.

Main part
We consider the relativistic R + R 2 action including all terms quadratic in the Riemann tensor R µνλρ and its simplifications here the R-term is the Einstein-Hilbert Lagrangian. The Λ-term is not essential in perturbation theory which we consider. M 2 P l = 1/(16πG) is the Planck mass squared, R µνρσ is the Riemann tensor, R µν is the Ricci tensor and R is the Ricci scalar. α, β and δ are coupling constants, D = 4 − 2ǫ is the dimension of the space-time within dimensional regularization [7]. ǫ is the regularizarion parameter and µ is the parameter with the dimension of a mass in dimensional regularization.
The Riemann tensor reads here are the Christoffell symbols Let us underline that dimensional regularization [7] is presently the only known continuous (not discrete like lattice) regularization of ultraviolet divergences appropriate for perturbative calculations and preserving gauge invariance of gravity.
The term containing the coupling δ in the Lagrangian (1) is usually omitted in the literature, see e.g. [3,4,8]. This is because of the Gauss-Bonnet identity The identity is valid only in four-dimensional space. But the dimension of the space-time in dimensional regularization is 4−2ǫ. Thus it seems that the term with the coupling δ must be preserved in the action to have renormalizability. From the other side it is in principle possible that one will invent fourdimensional continuous regularizatiom which preserves gauge invariance of gravitational Lagrangian. Then the term with the coupling δ should be omitted. The number of coupling constants in the Lagrangian most probably should not depend on the choice of regularization. In this case the term with the coupling δ should be omitted in dimensional regularization also. The point can be checked with direct calculations of counterterms of the Lagrangian. To establish the full picture it is most probably necessary to perform two-loop calculations, as it was with the establishing perturbative non-renormalizability of pure gravity mentioned in the introduction. Corresponding calculations are rather involved even at the one-loop level. This is a subject for a seperate publication.
It should be also mentioned that there is the known Regge-Wheeler lattice regularization which preserves a form of lattice diffeomorphism invariance.
We will work within perturbation theory. Thus a linearized theory is considered around the flat space metric here the convention in four dimensions is η µν = diag(+1, −1, −1, −1). Within dimensional regularization η µν η µν = D. Indexes are raised and lowered by means of the tensor η µν . Gauge transformations of gravity are generated by diffeomorphisms x µ → x µ + ζ µ (x) and have the form whith arbitrary functions ζ µ (x). Following standard Faddeev-Popov quantization [9], see also [10], one adds to the Lagrangian a gauge fixing term which can be chosen e.g. in the form here F µ = ∂ ν h νµ , ξ is the gauge parameter. Physical results, of course, do not depend on the allowed choice of the form of the gauge fixing term. One should also add the ghost term where C and C are ghost fields. Then one gets the generating functional of graviton Green functions here N is the normalization factor of the functional integral in the usual notation, J µν is as usual the source of gravitons. We work within perturbation theory, hence one makes the shift of the fields Perturbative expansion is in inverse powers of the Plank mass or in other words in powers of the Newton coupling constant G ∝ 1/M 2 P l . Let us obtain the graviton propagator. One takes the part of the Lagrangian quadratic in h µν and makes the Fourier transform µνρσ being projectors to the spin-2, spin-1 and spin-0 components of the field h µν : Here Θ µν = η µν −k µ k ν /k 2 and ω µν = k µ k ν /k 2 are transverse and longitudinal projectors correspondingly. We note that the expression (11) differs from the similar expression presented in [8] by the absence of ǫ-dependent terms.
To get the graviton propagator D µνρσ one inverts the matrix in square brackets of the expression (11): Then the propagator has the form Then one performs partial fractioning. The propagator takes the forn In four dimensions one obtains the following graviton propagator We will now consider classical quadratic gravity. In this case for a point particle having the energy-momentum tensor T µν = δ 0 µ δ 0 ν Mδ 3 (x) one gets the gravitational field [5] m 2 2 = M 2 P l /(−α − 4δ) and m 2 0 = M 2 P l /(2α + 6β + 2δ) are squared masses correspondingly of massive spin-2 and spin-0 gravitons. Cupling constants α, β and δ can be chosen to obtain positive masses. In [3,5] it was noted that masses can be chosen large enough to have an agreement with experiments.
Our propagator (19) reproduces the expression (20). One can see it by means of the calculation of the tree level Feynman diagram corresponding to an exchange of two point-like particles by a graviton.
The graviton propagator in the work [3] dos not produce the expression (20). It contains some technical errors. To see this one puts in the R + R 2 Lagrangian coupling constants equal to zero except the Newton coupling. The Lagrangian is reduced then to General Relativity. Hence the graviton propagator should also be reduced to one of General Relativity: where the gauge condition with ξ = 0 is taken for simplicity. Our propagator (19) reproduces the propagator (21) in this limit. The propagator of the work [3] has the factor 1 instead of 1/2 in the third term of the numerator of (21) in the corresponding limit.
The second term in the expression (19) for the graviton propagator has the non-standard minus sign. Hence one considers it as the massive spin-2 ghost. For renormalizability of quadratic gravity one must shift poles of all propagators in Feynman diagrams in the same way k 2 → k 2 + i0. Hence the spin-2 ghost must be considered as a state with negative metric [3]. That is why violation of either unitarity or causality within the R + R 2 model was claimed in [3,5].
But this massive spin-2 ghost is unstable. It unavoidably decays in two or more physical massless gravitons. The width of the decay is small. However independently of the value of this decay width spin-2 ghost particles do not appear as asymptotic states of the S-matrix elements. Only physical gravitons appear as external particles of the S-matrix amplitudes. Thus one concludes that unitarity is preserved in the R + R 2 model.
There is a statement about instability of theories with ghosts , i.e. their Hamiltonians are unbounded from below and they do not have stable vacuum states. This question was raised in [11] within Quantum Mechanics, see also [12] for the brief review. But this statement is proved only for Quantum Mechanical systems. Quantum Field Theory is quite a different story and renomalizability plus unitarity is enough to have a consistent theory.
To see this let us consider the graviton propagator in the opeator formalism: One transforms it to the momentum space and inserts the sum over the complete set of momentum eigenstates between two graviton fields. The states with negative norms in the sum have the extra factor −1. It gives the negative residue for the massive spin-2 ghost pole. There is another way to produce the negative residue for the spin-2 ghost. One can prescribe negative energy to this ghost. The expansion of the graviton fields into the creation and annihilation operators produces normalization factors 1/ √ −2k 0 . This is the reason for the negative residue for the spin-2 ghost. In this case of negative energy the Hamiltonian would be indeed unbounded and the vacuum state would be unstable.
But as it was mentioned above, one should choose the variant with negarive metric in order to have renormalizability in the theory [3]. Thus one has the consistent theory with the stable ground state. There are no reasons for a Hamiltonian to be unbounded from below if there are no states with negative energies.
It should be mentioned that the S-matrix by construction automatically satisfies the unitarity relation in theories having Hermitian Lagrangians [13]. To see it one considers the S-matrix in the operator formalism One introduces a function g(x) having the values in the interval (0, 1). This function describes intensity of interactions. Interactions are switched off if g(x) = 0. If g(x) = 1 then interactions are switched on. Interactions are switched on partly if 0 < g(x) < 1. One substitutes the product L(x)g(x) for the Lagrangian L(x). The S-matrix becomes the functional One splits the interaction region characterised by the function g(x) into an infinitely large number of infinitely thin segments ∆ i using the space-like surfaces t = const. Then one gets S(g) is defined as the limit The r.h.s. of (27) is a product taken in the chronological order of the segments ∆ j . Each factor in this product is unitary up to small terms of higher orders for sufficiently small ∆ j . These higher orders can be neglected in the considered limit. Hence the whole product is unitary. Unitarity of S(g) and of the matrix is proved. Sometimes one understands the following thing under unitarity. One derives from (23) the famous optical theorem stating that imagenary part of an amplitude of some forward scattering coincides up to a factor with the corresponding total annihilation crosssection where |i > is the scattering state, T is the scattering matrix: S = 1 + iT , and one assumes that all physical states |n > form a complete set in the theory n |n >< n| = 1.
From the other side one can calculate Im < i|T |i > directly from Feynman diagrams using Cutcosky cuts. Then one assumes that the result should coincide with (29). But if it does not happen it does not mean violation of unitarity. It only means that physical states in the theory do not form a complete set (30) and the complete set is formed by physical plus unphysical states. Unitarity of theories with negative metric states was previously considered in [14,15], see also references therein. Question of causality were also considered there.
We would like to note that the tree level graviton propagator (19) is modified by the summation of the chain of one-loop insertions. As it was already mentioned above the second term of the propagator (19) has the minus sign. Therefore the summation of the one-loop insertions with the massless graviton in the loop will shift the pole of the spin-2 ghost from the value Here Γ is the width of the spin-2 ghost decay into the pair of massless physical gravitons. This complex pole is located on the unphysical Riemann sheet. It is analogous to the known virtual level of the neutron-proton system with opposite spins of nucleons [16].
It should be noted that one loop corrections in quadratic gravity were studied in [17].
We would like to underline that we consider not pure R 2 theory but the R + R 2 theory where the R 2 terms are added to the Einstein-Hilbert Lagrangian. Gravitational constants α, β and δ of these terms in the Lagrangian can be chosen suffisiently small to ensure that quadratic gravity will be practically indistinguishable from General Relativity at astrophysical and cosmological scales. This is independently of the discussed above in detail question wether the coupling δ is exactly zero or not. The R 2 terms are introduced only to have renormalizability of quadratic gravity which is valid in particular for arbitrary small couplings α, β and δ.
We have analyzed only purely gravitational R + R 2 action. The inclusion of the matter fields in the Lagrangian is straightforward and does not change conclusions.

Conclusions
We have proved unitarity of quantum gravity with the R + R 2 action. This model was previously shown to be renormalizable in the work [3]. The parameters of quadratic gravity can be adjusted to ensure that the theory will be practically indistinguishable from General Relativity at astrophysical and cosmological scales.
One can conclude that the R + R 2 model is an appropriate candidate for the fundamental quantum theory of gravity.

Acknowledgments
The author is grateful to the collaborators of the Theory Division of the Institute for Nuclear Research for valuable discussions.