^{1}

^{*}

^{2}

^{*}

The first mathematically consistent exact equations of quantum gravity in the Heisenberg representation and Hamilton gauge are obtained. It is shown that the path integral over the canonical variables in the Hamilton gauge is mathematically equivalent to the operator equations of quantum theory of gravity with canonical rules of quantization of the gravitational and ghost fields. In its operator formulation, the theory can be used to calculate the graviton S-matrix as well as to describe the quantum evolution of macroscopic system of gravitons in the non-stationary Universe or in the vicinity of relativistic objects. In the S-matrix case, the standard results are obtained. For problems of the second type, the original Heisenberg equations of quantum gravity are converted to a self-consistent system of equations for the metric of the macroscopic space time and Heisenberg operators of quantum fields. It is shown that conditions of the compatibility and internal consistency of this system of equations are performed without restrictions on the amplitude and wavelength of gravitons and ghosts. The status of ghost fields in the various formulations of quantum theory of gravity is discussed.

In the works [1,2] it is shown that the observed accelerated expansion of the Universe (the Dark Energy effect) [3,4] may be a consequence of macroscopic effect of the quantum theory of gravity which is the condensation of gravitons on the horizon scale of the non-stationary Universe. The theory of this macroscopic quantum effect is based only on the most general properties of quantum gravitational field in the non-stationary Universe, which are Bose-Einstein statistics, zero rest mass of gravitons, their conformal non-invariance and one-loop finiteness of quantum gravity off the graviton mass shell. (Oneloop finiteness on the graviton mass shell was proved by T’Hooft and Veltman in 1974 [

This work aims to show that the equations of one-loop quantum gravity obtained in [

For curvatures that are less than Planck curvature, quantum cosmology is represented as a theory of gravitons in the macroscopic space time with the self-consistent geometry. The quantum state of gravitons is determined by their interaction with a macroscopic field, and the macroscopic (background) geometry, in turn, depends on the state of gravitons. The background metric and graviton operator appearing in the self-consistent theory are extracted from the unified gravitational field, which initially satisfies exact equations of quantum gravity. The classical component of the unified field is a function of coordinates and time by definition. The quantum component of the same unified field is described by a tensor operator function, which also depends on coordinates and time. Under such a formulation of the problem, the original exact equations should be the operator’s equations of quantum theory of gravity in the Heisenberg representation. The present work is devoted to obtaining such equations and their identity transformation to the self-consistent system of classical and quantum equations with no restrictions on the wavelength and amplitude of the graviton field.

The Heisenberg representation of quantum theory can be introduced only by its Hamilton formulation. For the theory with constraints, we deal with a generalized Hamilton formalism. The generalized Hamilton formalism of the theory of gravity was constructed by Dirac [7,8] and Arnowitt, Deser and Misner [

At the present time, an operator version of equations proposed in [7-9] is usually considered as the equations of quantum gravity in the Heisenberg representation. One should, however, bear in mind that local gauge conditions used in [7-9] do not remove the degeneracy completely, i.e. quantum fields of inertia are present in the equations of theory. Therefore, perturbative S-matrixes calculated in these theories do not satisfy the unitarity condition. As is known, the existence of this problem drew the attention of Feynman [

Despite this fact, the operator Einstein equations in the Heisenberg representation with the local gauges (without ghosts) that do not remove the degeneracy completely and do not satisfy the test for unitarity S-matrix, were used in discussing of cosmological applications of quantum gravity, i.e. in the quantum theory of gravitons in the non-stationary Universe (see, e.g., [15,16]).

In connection with the situation that exists in the scientific literature, we emphasize the following circumstance. Except the test on a unitarity of S-matrix, there is no other criterion of correctness of the equations of quantum gravity, and, as a consequence, the correctness of the results obtained from these equations. Therefore, the Hamilton equations in the Heisenberg representation with the canonical quantization rules must be strictly deduced from the exact equations of quantum gravity satisfying the above test. In this paper, we show that these equations do exist under condition that the Faddeev-Popov ghosts are among the elements of the Hamilton formalism of quantum gravity.

We obtain the Hamilton equations in the Heisenberg representation by identity transformation of the path integral over the canonical variables. Such an integral was proposed by Faddeev [

The procedure for obtaining of these equations is as follows. In the first step, the gauge of Faddeev’s path integral [

The Heisenberg representation in the quantum theory of gravity (outside of perturbation theory) exists in the synchronous gauge in the normal coordinates. This gauge is an analogue of Hamilton gauge in the Yang-Mills theory. Basically, the operator Hamilton formalism with the canonical quantization of gravitational and ghost fields is to solve problems in quantum cosmology but can be also used for the calculation of S-matrix where it leads to standard results.

Upon obtaining of the exact equations of the theory of gravity in the Hamilton form, the gravitational field can be regarded as the deviation of the metric from the metric of Minkowski space. Normal coordinates of gravitational field of are given by exponential parameterization of the density of the contravariant metric [

The density of the gravitational Lagrangian as a function of normal coordinates reads

where dots denote a full derivative which does not contribute to the equations of motion. A variation of the action of the gravity theory over the normal coordinates leads to the Einstein equation

The same equations can be rewritten in the form

In accordance with the general properties of Einstein equations, are equations of motion , and are equations of constraints. Note also the following form in which the Bianci identity can be presented:

The Hamilton gauge is given by additional conditions

In normal coordinates, the conditions (6) are reduced to. The specific of the Hamilton gauge is that the Bianchi identity (5) takes the differential form of conservation laws over the solutions of the equations of motions. It reads:

After the gauge is applied, the Lagrangian and the equations of motions take the form:

(9)

The equations (9) are obtained by variation of the gauged action over the variables. After the gauge is applied, Lagrangian (8) does not contain variables whose variations would produce equations of constraints. However, these constraints are restored in solutions of the equations of motion (9) by imposition of conditions on the numerical values of integrals of motion. As follows from (7), first integrals of the equations of motion (9) are of the form

where C_{i} are formally numerical constants whose values are fixed by initial and boundary conditions. The general covariant nature of the equations of constraints imposes the obvious limitations over values of C_{i}: only those initial and boundary conditions are allowable for the metric in the Hamilton gauge in which constants C_{i}, appearing in (10), are equal to zero.

Thus, the formalism of the theory, based on a gauged Lagrangian (8), differs from the formalism of non-degenerate field theories only by limitations that are related to the initial and boundary conditions but not to the equations. This fact allows to construct the Hamilton formalism easily. Generalized momentums are of a standard definition over functional derivatives:

where. The gauged action is reduced to the form

where

is the density of canonical Hamiltonian by which the full Hamiltonian is calculated. The variation of action (12) over the generalized coordinates and momentums (as over independent variables) leads to the Hamilton system:

As it should be, the Hamilton formalism is mathematically equivalent to the Lagrange formalism, and this can be seen from the fact that the system of Equations (14) is easily reduced to the Lagrange Equations (9). The Hamilton constraint (10) with i = 0, C_{0} = 0 becomes the formula for calculating of the energy of the gravitational field through surface integral:

The general formulae for infinitesimal transformations reads

For the Hamilton gauge (6), the expression (16) leads to the equations for the parameters of residual transformations:

The solution of Equations (17) is presented in the form

where

is a function of spatial coordinates defined by arbitrarily;

is the solution of the d’Alambert wave equation

Function (18) defines the fields of inertia:

According to (19), fields of inertia

are trivial in the sense that and do not interact with gravity. These fields describe the given motions of the reference frame. The motion corresponding to a field is a given stretch of coordinate grid, and the motion corresponding to a field is a given rotation. The lack of connection with the metric implies that these fields of inertia do not affect the solution of Einstein’s equations; for this reason they can be excluded from consideration.

The function (satisfying the wave Equation (21)) defines a nontrivial field of inertia

. The nature of this field is following.

The condition synchronizes the clocks, and component sets the rate of a clock. According to (6), the Hamilton gauge is synchronous, and the rate of a clock is rigidly connected with the dynamics of 3 volumes:, where is the determinant of spatial metric. Effects of continuous clock adjustment to a new physical situation are reflected in the measured values of 3-dimensional component of the metric in the form of wave field of inertia. It is propagating with the same speed with which a clock located at the point of measurement is synchronized with the clocks located at the neighboring points.

Note that in the theory of gravity in the Hamilton gauge the non-trivial field of inertia is of one only independent internal degree of freedom.

The procedure of hamiltonization of classical theory of gravity in the Hamilton gauge is carried out in Section 2.2. Its content and sequence of procedures exactly match to the hamiltonization procedure of Yang-Mills theory [

Differences between the Yang-Mills theory and the theory of gravity appear during the phase of quantization. According to [

where

is the operator of the system of Equations (17) for the parameters of residual transformations. The expression (23) is obtained from Faddeev’s path integral (see formula (5.6) in [

The computation of the determinant of operator (24) leads to

The structure (25) completely corresponds to the classification of the fields of inertia in the Hamilton gauge. The localization leads to the Lagrangian containing only one non-trivial ghost field, which does interact with the gravity. The direct localization leads to the path integral over the configuration space of ghost fields:

where

is the density of the ghost Lagrangian. A path integral over the phase space is a mathematical equivalent to (26). It reads:

where

is the density of the ghost Hamiltonian. Substitution of (28) to (23) gives a path integral over the extended phase space of gravitational and ghost variables:

Expression (30) is the standard definition of the matrix element of the operator of evolution that defines the dynamics of state vectors in the Schrödinger representation or, equivalently, the dynamics of operators in the Heisenberg representation. Therefore (30) allows immediately writing down the terms of the canonical quantization and equations of quantum gravity in the Heisenberg representation.

The graviton sector is:

The ghost sector is:

In (32), (34)

is the full canonical Hamiltonian of gravitational and ghost fields. The canonical quantization (31) and (33) and solution of operator Equations (32) and (34) are procedures that are of mathematically equivalence to the computation of path integral (30).

It is easy to check that differential conservation laws are satisfied for the solutions of Equations (32) and (34)

where is the density of Einstein tensor,

is the density of the energy-momentum tensor of the ghosts. Integral conservation laws are presented as conditions imposed on the Heisenberg state vector, i.e. on the vector of initial quantum state:

The Hamilton constraint acquires a meaning of the equation for eigenvalues and eigenvectors of the total Hamiltonian of graviton and ghost fields:

A comparison of quantum Einstein’s equations (32) in the Hamilton gauge in the Heisenberg representation with the classical Equations (14) shows that the transition from classical to quantum theory is not limited to the replacement of C-numeric functions with the operator functions satisfying the canonical commutation relations. In the operator equations of quantum gravity, the operator functions of Faddeev-Popov ghosts automatically appear. They satisfy the canonical anti-commutation relations. The role of ghosts in the Heisenberg equations of motion (32) is exactly the same as in the path integral (30): the ghosts compensate for the effects of interaction of a true gravitational field with a field of inertia. However, in the complete system of Equations (31)-(34) the ghosts formally are of the status of additional dynamic variables, i.e. their status is equal to the status of dynamical variable of gravitational field.

At the operator level, the structure of equations of the theory in the Heisenberg representation with canonical quantization of gravitons and ghosts ensures the existence and carrying out routine operations on the transition from (31)-(34) to the representation of interaction, and the construction of S-matrix in the framework of the theory of perturbations. In this theory, gravitons are described by the deviation of the metric from the metric of Minkowski space. Asymptotic states for gravitons are constructed over the vectors of the Fock space for the quanta of 3D-transverse tensor modes, and asymptotic states for the ghosts are set as vacuum states. Ghosts are present only in the Green function that compensates for the contribution of the wave field of inertia in the Green function of gravitons. The above S-matrix is identically coincides with S-matrix constructed from the path integral over canonical variables (30) or from the path integral over the Lagrange variables (40). Match of S-matrixes constructed in various ways is a direct illustration of the existence of the Heisenberg representation for the quantum theory of gravity in the Hamilton gauge. Application of this theory to the description of a macroscopic number of gravitons in the Riemannian space with a self-consistent geometry is discussed in Section 5.

The quantum theory of gravity in the Heisenberg representation formulated in (31)-(34) is mathematically equivalent to the calculation of gauged path integral over the Lagrange variables of gravitational and ghost fields:

where are unit time-like vectors defining a Hamilton gauge. Of course, (40) is a specific case of the Faddeev-Popov-DeWitt gauged path integral [14,20,21]. Structurally, (40) is obtained from the gauge invariant expression

through the expansion of unit in the Hamilton gauge and factorization of volume of the gauge group. In the context of existence of the Hamilton formalism, we note that (40) can be regarded as the result of identity transformations of the expression (30). We do not describe these transformations due to their rather obvious.

The expression (40) clearly reflects the unique property of the Hamilton gauge: the effective Lagrangian of gauged theory is formally defined by the covariant expression

Because of this property, the variation of action over normal coordinates and Lagrange variables, , leads to a system of generally covariant equations for the gravitational and Grassmann scalar fields:

The left-hand-side of Einstein’s Equations (43) satisfies the Bianchi identities, and the energy-momentum tensor of Grassmann fields appearing in the right-handside is conservative over the equations of motion (44). The existence of these differential identities allows separating Einstein’s equations to the equations of motion of gravitational field and equations of constraints. It is also clear that after an explicit account of the Hamilton gauge, the gauged equations of motion are reduced to the form (32) and (34). Thus, the Hamilton gauge ensures the existence of Heisenberg representation for the equations of quantum theory of gravity both in Hamilton and in the Lagrange form.

The exact equations in the Heisenberg representation in the Hamilton gauge (32), (34) together with postulates of canonical quantization of (31), (33) claim the status of the theory formulated only on the basis of the first principles of quantum gravity. We turn to the approximate theory of macroscopic system of gravitons in the macroscopic space time with the self-consistent geometry. In this theory, the interaction of gravitons with the classical gravitational field is taken into account exactly, and the perturbation theory over the amplitude of quantum fluctuations is only used when describing the graviton-graviton interaction.

The averaging of normal coordinates over the Heisenberg state vector in the general case provides non-zero C-numeric function:. A quantum fluctuation is defined as the difference between the operator and its average value:. Substituting in the formula (1) for the operator of the density of covariant metric gives

A matrix exponential built over C-numeric functions, will be called the density of contravariant metric of the macroscopic space time. In this space, we normally enter the covariant metric, the connection, the curvature and define covariant derivatives. The quantum fluctuations are endowed by the properties of a tensor in the macroscopic space time. The Hamilton gauge is partitioned into gauge of the macroscopic metrics and gauge of the quantum field:

Transformations of equations can be made in the Hamilton formalism (32), (34) as well as in the Lagrange formalism (43), (44). To elucidate the general properties of the theory, the Lagrange formalism is more convenient. Substituting (45) in (43) gives:

where

is the energy-momentum tensor of gravitons;

is the energy-momentum tensor of ghosts. The Lagrange equations for the ghosts read

Semicolons in (47)-(50) denote covariant derivatives in the macroscopic spacetime.

The transformation of shift is of physical meaning if the true degree of freedom of the quantized gravitational field are contained in the operator function, and the macroscopic gravitational field is self-consistent and serves as a way to describe collective interactions in the system consisting of a large (macroscopic) number of particles. The ability to carry out specific calculations in the approximation of a self-consistent field appears when symmetry properties of the macroscopic space time allow the introducing of a system of basis state vectors, and on this base to construct a state vector of gravitons and ghosts of a general form. The Einstein equations for the macroscopic geometry are obtained from (47) by averaging over the state vector:

The equations for gravitons in the Heisenberg representation are obtained by subtraction (51) from the exact Equations (47):

In terms of the path integral formalism, the existence of Equations (50)-(52), is provided by factorization of the measure of integration and consistent integration, first over the quantum fields of gravitons and ghosts, and then over the classical field. It can be seen from following. After an additive transformation is used, the measure of integration in normal coordinates and Hamilton gauge is represented as:

Formally, the exact integration of quantum fields, , corresponds to the solution of quantum Equations (50) and (52), and the approximate integration of the classical field corresponds to the solution of the classical Einstein Equation (51).

The system of classical and quantum equations satisfies three differential identities. The left-hand-side of (51) by definition satisfies the contracted Bianchi identities in the macroscopic space time. In the right-hand-side, the condition for conservation of the expectation value of energy-momentum tensor of ghosts satisfies to the equations of motion for the ghosts (50). The expectation value of energy-momentum tensor of gravitons is of the similar properties. The calculation of a covariant divergence yields:

As seen from (54), the condition of conservation is performed over the quantum equations of motion (52).

In the self-consistent theory of gravitons, the nonperturbative effect of the formation of the geometry of macroscopic spacetime is described in semi-classical/semiquantum level: in the Einstein Equation (51) the source of the classical field is the expectation value of energy-momentum tensor of quantum fields. In the quantum sector, the order of a perturbation theory is given by maximum degree “n” in which the operators of quantum fields of gravitons and ghosts are in the quantum equations of motion (52) and (50). In the same order of the perturbation theory, terms of the order of “n + 1” are accounted in the energy-momentum tensor. As seen from (54), the condition of conservativeness of the graviton energy-momentum tensor takes place in each order of perturbation theory separately. This fact is the basis for the application of the perturbation theory.

A specific feature of quantum gravity in the Heisenberg representation is the inevitable presence of ghost fields in the canonical formalism of the theory. Formally (in mathematical sense), the ghost sector is a consequence of the lack of gauges fully removing the degeneracy with respect to the residual transformations of the group of diffeomorphisms. Physically, this means that in the theory of gravity the fields of inertia are described together with the true gravitational field. In the classical theory of gravity, this does not lead to any problems with interpretation because the fields of inertia and gravitational fields are locally separable. Therefore, the deterministic evolution of each of these fields always allows finding invariant characteristics of local effects. The specificity of the quantum gravity lies in the fact that because of the lack of gauges completely removing the degeneracy, the existing formalism leads to the quantization of fields of both types (i.e., the uncertainty principle is turned out to be applied not only to the gravitational field but also to the field of inertia).

Of course, in the quantum gravity the question about the status of ghost fields arises as well as the question about the rules of computing the contribution of ghosts to the observables. Answers to these questions are contained in the formalism of the theory. The general answer is that in the quantum gravity the ghosts and fields of inertia together form a physical effect. Details of this procedure depend, however, from the specific of a physical problem.

Let’s start with the situation, which is well known and understood. The Faddeev-Popov path integral [10,14] is initially defined as a method of calculating of the transition amplitude between asymptotic states in the framework of perturbation theory. Here should be noted that at the very formulation of the problem, the fields of inertia are unobservable in the asymptotic states, which is formally provided by appropriate boundary conditions. The impossibility of observing of the inertia fields is automatically accompanied by the lack of asymptotic states with the nonzero ghost occupation numbers. Formally, the lack of ghosts is ensured by the stability of the ghost vacuum and by the possibility of conducting of the various “expansions of unit” during the transformations of the gauged path integral. However, the ghosts as well as the fields of inertia are present in the virtual states in the region of interaction, and their quantitative role in the formation of the observable amplitude of the graviton scattering is given by the relevant Feynman diagrams. In various gauges, different inertial and ghost fields are appeared but the total contribution of these fields in the graviton sector of S-matrix is gauge invariant.

In quantum gravity in the Heisenberg representation there is a possibility of the task, which physically is fully equivalent to the problem of the S-matrix calculation. Equations (31)-(34) can be used to describe the gravitational field whose potentials are only weakly perturb the metric of Minkowski space. In this situation, one can also enter the asymptotic state, and then, specifying the flux of gravitons in the initial asymptotic state to calculate the flux of gravitons in the final asymptotic state. The ghosts together with the inertia fields influence the formation of finite flux of gravitons, and this effect is taken into account in the process of solution of the operator equations of motion.

There are two specifics in the solution of the problem described above. First, the Heisenberg representation in the quantum gravity as well as in the Yang-Mills theory exists only in special gauges, so that the quantum fields of inertia and ghost sector are not arbitrary. Second, in the Equations (31)-(34), there actually are no mathematical indications of any special status of ghost fields in contrast to the original path integral where they are. Formally, in these equations, the quantum wave field of ghosts acts as a second dynamic subsystem. Therefore, the exclusion of the ghosts of the asymptotic states is an additional condition. The mathematical consistency of the selection rule of the ghost-free asymptotic states is provided by the stability of the ghost vacuum in the theory of perturbations on the background of the Minkowski space.

We now turn to the most non-trivial situation, which is the role of the ghosts in the quantum theory of the graviton macroscopic system. Let us consider the Equations (50)-(52), in the one-loop approximation. In the macroscopic Einstein equations, the energy-momentum tensors of the gravitons and ghosts are quadratic in quantum fields:

The quantum fields of gravitons and ghosts are described by the following linear equations:

As is seen from (55)-(57), the one-loop quantum gravity is not confined to the theory of free fields with the spin J = 2 in the curved space time. The objects of the theory are necessarily two fields which are tensor field of spin J = 2 and the scalar Grassmann field with spin J = 0. Moreover, as is seen from the solutions of the equations (55)-(57), the gravitons and ghosts must have an equal status of dynamical subsystems of macroscopic quantum gravity system [1,2]. This somewhat unusual structure of the theory (in comparison with the theory of physical fields with spin) needs for additional comments.

The first thing that should be stressed is the fact that the inclusion of ghosts in the list of “physical” fields is the result of regular mathematical transformations. These transformations are as follows: the Faddeev path integral over the canonical variables (23) ® equations of quantum gravity in the Heisenberg representation (31)-(34) ® the allocation from the normal coordinates of the gravitational field of their mean values and the identity transformation of the exact equations in the Heisenberg representation to the self-consistent system of classical and quantum Equations (50)-(52). There is no question about the existence of such transformations because they obviously are in a formal mathematical sense. The question is whether or not these transformations lead to a new physical content of quantum theory of gravity in the last stage of computations (they are initially absent in the gauged path integral)?

The answer to this question is “Yes”. The physical content of the theory is really changed but the reason for this change is absolutely transparent. The source path integral is defined as the mathematical object on the background of Minkowski space with the asymptotic boundary conditions and under the assumption of stability of graviton (and ghost) vacuum. From a mathematical point of view, the transition from this integral to the theory of the macroscopic quantum gravity system (50)- (52), represents an extrapolation of the theory to a new physical area with the properties that differ from the area in which the theory was formulated initially.

First of all, it should be noted that the formation of the metric of the classical curved spacetime with the selfconsistent geometry is a significantly non-perturbative effect of collective interactions in the macroscopic system. Cosmological applications of quantum gravity make the abandonment of the Minkowski space (appearing in the theory of S-matrix) inevitable in favor of the real existing curved spacetime. Direct consequences of a realistic formulation of the problem of describing of a macroscopic system of massless and conformal noninvariant quantum fields are disappearance of asymptotic states and instability of the graviton-ghost vacuum. Any observer located inside such a system (i.e., in fact, inside the real Universe) is in the area of interaction. It is rather obvious that for such an observer the non-perturbative effects of vacuum instability are priority subjects of research.

Next, you need to bear in mind that in General Relativity the observer is represented by the fields of inertia which are inevitably quantized (in quantum gravity) in the absence of gauges completely removing the degeneracy. Therefore, in terms of standard quantum gravity, among of locally observable physical quantities, there necessarily are the observables formed by quantum fields of inertia and ghosts jointly. The instability of the vacuum in the spectral region where the wavelength of quantum fields is comparable to the radius of 4-curvature does not allow to fix the quantum state of these fields by zero occupation numbers. Exactly for this reason, the ghost fields are beginning to perform the role of the second quantum subsystem which is dynamically equal to the graviton subsystem.

Thus, on the one hand, the ghosts perform their standard function in the theory of macroscopic quantum gravity systems, i.e. the ghost contributions together with the contributions of the inertia fields form the observables. But, on the other hand, this function occurs in the non-standard conditions of the absence of asymptotic states and the instability of the vacuum.

This specific of the ghost sector is fully reflected in the equations of one-loop quantum gravity (55)-(57). In these equations, the interaction of gravitons with each other is taken into account through a self-consistent field. According to the most general concepts of the quantum theory of gravity (existing on today), the one-loop effects of the interaction can not be properly described without ghosts. In the Equations (55)-(57), the inertia fields are accumulated in the macroscopic self-consistent field. Therefore, the ghost contribution to the formation of observables is described by the obvious and only a way possible through the influence of ghosts on the self-consistent field. This contribution is taken into account through the presence of an averaged energy-momentum tensor of the ghosts in the macroscopic Einstein Equations (55).

The exact solutions of the Equations (55)-(57) describing macroscopic effects of quantum gravity in the homogeneous and isotropic non-stationary Universe are obtained in [1,2]. A possible role of graviton, ghost and instanton condensates in the formation of observable Dark Energy is also discussed in [1,2].

In this paper we show that the quantum theory of gravity can be formulated as operator equations in the Hamilton gauge in the Heisenberg representation. The postulate of quantization in this theory is given by canonical commutation and anti-commutation relations for the generalized coordinates and momentums of gravitational and ghost fields. The equations of theory (31)-(34) (or the equations obtained from these Equations (50)-(52)) are used to describe a macroscopic system of gravitons forming a macroscopic spacetime with the self-consistent geometry.

There is the reason served as one of the main motives in writing this paper. During the last 30 years (since 1977 to 2008) many papers were published on the quantum theory of gravitons taking into account their backreaction on the cosmological background (quantum backreaction theory). To the best of our knowledge, there are no mathematically consistent papers correctly taking into account the structure of the existing quantum theory of gravity. In particular, the papers [15,16,22-24] published in Physical Review D are erroneous. The typical errors are following:

1) Use of linear parameterization of the metric fluctuations leading to the non-self-consistent system of classical and quantum equations. This is a problem facing the authors of all works on the theory of gravitons, but they tried not to direct attention to the discussion of this matter. Existence of the problem is exhaustively documented in the work [

2) Incorrect gauging (a gauge is imposed not on the full metric but on its fluctuations only). Such a method is a serious and obvious error if it is used in applications in which the macroscopic background has no a global 4- symmetry, permitting the separation of inertia fields. Formally, the mathematical nature of the error is as follows. The group of diffeomorphisms acting in the original Riemannian space with the metric (1), cannot be factorized into the group of diffeomorphisms of the background space and the gauge group of gravitational waves after the separation in this metric of background and wave modes according to (45). However, this error does not affect the erroneous results of [15,16,22-24] because in the quantum theory the main consequence of gauging is the emergence of the ghost sector, which in these works is still ignored.

3) It is an unpardonable error of [15,16,22-24] to completely disregard the most important property of pure quantum gravity (without matter fields) which is its one-loop finiteness. Meanwhile, the one-loop finiteness is a rigorous and well-known T’Hooft and Veltman result [

4) Incorrect use (or ignore) the ghost sector. The lack of a one-loop finiteness in [15,16,22-24] versions of the theory, is mathematically uniquely associated with either ignoring ghost sector [15,16,22-24], or with a mathematically incorrect work with this sector [

All these works completely lose any meaning after renormalization of gravitational Lagrangian by quadratic counterterms off the mass shell. This procedure modifies the original definition of the graviton and makes the theory mathematically inconsistent. This is an unavoidable and unacceptable internal contradiction of one-loop quantum gravity (see Section 12 of [

We would like to express our deep appreciation to Ludwig D. Faddeev of the Steklov Mathematical Institute for graciously agreeing to read our manuscript and verifying the correctness of our approach before submission to the journal. We are deeply grateful to Mikhail Shifman and Arkady Vainshtein of the University of Minnesota for discussions of the structure and content of the theory. Also, we would like to express our deep appreciation to our friend and colleague Walter Sadowski for invaluable advice and help in the preparation of the manuscript.