On the Lagrangian structure of metric nonlinear gravity theories: the redundant Weyl tensor dependence

In the framework of metric gravity theories we investigate whether the Lagrangian may explicitly depend on the Weyl tensor, through an additive term being an arbitrary function of its quadratic invariant. We assume that for a physically viable Lagrangian one should be able to identify the degrees of freedom with a multiplet of fields with definite masses and spins. Using a well-established method to obtain a (fully equivalent and covariant) second-order formulation of the field equations by introducing a multiplet of conjugate fields, we study whether these can be interpreted as a system of interacting fields which are physically independent (in the sense that in flat spacetime each field could be in an excited state independently of the other fields). Lagrangians depending only on the curvature scalar R are known to fulfill this criterion, being equivalent to a model with a scalar field interacting with the metric. For a quadratic Lagrangian depending on the full Ricci tensor, instead, the criterion is fulfilled only if the coefficients are such that the scalar momentum is excluded and gravitation is described by the metric and a massive spin-two field. The presence of a Weyl tensor contribution does not add physical degrees of freedom in the form of classical fields with definite mass and spin in flat spacetime; this outcome is independent of whether a Ricci term contribution to the Lagrangian is present. In conclusion, though a generic (purely metric) Lagrangian has eight degrees of freedom, in the second-order equivalent picture the physical fundamental fields can either be a metric and a scalar field, or a metric and a spin-two field, hence the maximum number of attainable field-theoretical degrees of freedom is seven. In the perspective of the equivalence with a second-order field theory, any Weyl tensor contribution to the Lagrangian is unphysical and redundant.

In the framework of metric gravity theories we investigate whether the Lagrangian may explicitly depend on the Weyl tensor, through an additive term being an arbitrary function of the quadratic invariant C α βµν C α βµν . We assume that for a physically viable Lagrangian one should be able to identify the degrees of freedom with a multiplet of fields with definite masses and spins. Using a well-established method to obtain a (fully equivalent and covariant) second-order formulation of the field equations by introducing a multiplet of "conjugate momenta", we study whether these can be interpreted as a system of interacting fields which are physically independent (in the sense that in flat spacetime each field could be in an excited state independently of the other fields). Lagrangians depending only on the curvature scalar R are known to fulfill this criterion, being equivalent to a model with a scalar field interacting with the metric. For a quadratic Lagrangian L = f (g µν , R, R αβ ) depending on the full Ricci tensor, instead, the criterion is fulfilled only if the coefficients are such that the scalar momentum is excluded and gravitation has seven degrees of freedom (carried by

Introduction
Alternative theories of gravity have appeared soon after the advent of general relativity, the first of them was Weyl's theory (1919) and since then they have been copiously created.
This phenomenon on one hand resembles the situation in other branches of physics, for example particle physics, but on the other is very different. In particle physics, after the discovery of the Higgs boson, there is no single direct experimental evidence that the Standard Model (SM) should be enlarged or modified. However, we know that it should: first, because there are indirect observational hints from cosmology like matter-antimatter asymmetry in the Universe and second, theoretical extrapolation of SM to high energies shows pathologies in the pure SM like instabilities or Landau poles. Therefore there are many proposals of extensions of SM, either within particle physics, mostly based on supersymmetry (or conformal symmetry), or more radical, like string theory, that were designed to cure these problems. Within particle physics, under the assumption of renormalizability, our choice of possibilities is limited to different gauge groups or different field content; within string theory the fundamental theory is essentially unique but through compactification we can have an almost continuous manifold of effective low-energy 4-dimensional theories. In either case we know that SM should be extended.
In gravitational physics the situation is different. From the experimental point of view we also don't need any competing theory to Einstein's GR because of lack of any experimental evidence pointing to a necessity of supplementing the Einstein-Hilbert action with higher order terms. However, reasoning behind any modification of the simplest GR is not as convincing and straightforward as in particle physics for two reasons.
First, higher order terms in Einstein's theory can be relevant only at gravitational curvatures many orders of magnitude bigger than what we could imagine as even indirectly observable in any foreseeable future (inside black holes or in the very early Universe). Second, although there are hints that such terms are necessary in quantum theory of fields coupled to gravity we don't have quantum theory of gravity itself to make such a picture fully consistent.
The theoretical arguments pointing to a necessity of adding higher order terms in gravity come from two sources. The first one is string theory or (not yet fully constructed) more general UV completion of the theory of particles, called M theory. The theory of closed strings has a graviton as its massless excitation. We can calculate N-point amplitudes involving gravitons as external legs and these amplitudes can be translated as coming from the effective theory of gravity with higher order terms in the Riemann curvature expanded around the flat geometry. There is however another approach leading directly to these corrections without any expansion: the calculation of the (2-dimensional) conformal anomaly of string propagating on a curved background. Then the requirement of vanishing of this anomaly leads directly to the Einstein equations with higher order corrections [1].
The second argument that points to higher order corrections in GR comes from the theory of particles coupled to gravity without invoking any UV completion. Then the loop calculations with massive particles on internal lines and gravitons as external legs lead at low curvatures to local higher order corrections in the Riemann tensor. These calculations produce infinities (in contrast to string theory where amplitudes are finite due to modular invariance) and require renormalization. Unfortunately, gravity is nonrenormalizable and at each level we have to introduce a new renormalization prescription what makes such a theory useless from the physical point of view. Another issue is a (D-dimensional) conformal anomaly that also contains higher powers of the Riemann tensor but even at low curvatures the anomaly is non-local and therefore has to vanish because of potentially disastrous consequences for observable theory of gravity [2].
All these arguments point to the fact that any quantum theory of particles coupled to gravity leads to an effective theory of gravity with higher order corrections. All effective theories have a limited application range -they are good approximations only in certain interval of energies. In this case (barring a possibility of conformal anomaly which is non-local even in the IR regime) the limiting energy is presumably close to the Planck scale where quantum gravity effects start to be important and/or more massive particles should be added. It seems worthwhile to consider consequences of adding higher order corrections below the Planck scale, i.e. at a purely classical level. For the purpose of this paper it is important to emphasize that the effective theories even with higher order terms or non-localities remain bona fide field theories (either classical or quantum). Therefore, assuming that there exists a quantum version of these theories, in some appropriate limit of coupling constants approaching 0 these theories should admit an interpretation in terms of masses and spins (except in 2 dimensions where there may exist some 'pathologies' associated with the fact that the little group has continuous representations). Even if the theory admits for non-vanishing coupling constant g nontrivial classical solutions, like solitons or instantons, that do not allow to take the limit g → 0, in the quantum interpretation they are not treated as separate particles but as an (infinite) collection of excitations out of a vacuum that are saddle points in the path integral (either Euclidean or Lorentzian).
In this paper we are dealing with one class of alternative theories, which differ from GR in only one axiom: in gravitational field equations. Instead of Einstein-Hilbert Lagrangian, L = R, one assumes that L is some smooth scalar function of the Riemann tensor. These theories have no commonly accepted name and they are dubbed either as 'metric nonlinear gravity' (NLG, Lagrangians nonlinear in the curvature) or 'higher derivative theory' (though the field equations need not be of fourth order); we shall use the first term. A revival of these theories has taken place after the discovery in 1999 of the apparent accelerated evolution of the universe. Most of the researchers attempting to account for this acceleration without resorting to the mysterious dark energy or the fine tuned cosmological constant, have applied various f (R) gravity theories (see [3] for a review). However, once one rejects the simplest Lagrangian L = R, one is not restricted to a function f (R) of the curvature scalar, and one may investigate the generic case of L = f (g µν , R α βµν ) and pass over any dependence on derivatives of the tensor. Hence, prior to any cosmological applications of any NLG theory one should firmly establish on physical grounds (by theoretical arguments since the observational data are insufficient) the correct Lagrangian. From a mathematical viewpoint there is nothing inappropriate in assuming that L depends on R, the Ricci tensor R µν and the Weyl tensor C α βµν apart from that the field equations become formidably complicated in comparison to Einstein's field equations. Yet from a physical viewpoint severe doubts arise. All metric nonlinear gravity theories are inherently ambiguous in their physical interpretation due to possibility of performing various redefinitions of their dynamical variables and only Einstein's GR is free of these ambiguities: this fact was gradually discovered in a series of works (see [4,5,6] for full references).
It has been shown by various authors (starting with Stelle [7]) that a Lagrangian depending on the Ricci tensor of a metric -in dimension four and without external matter -corresponds, generically, to a model with eight degrees of freedom (in the fieldtheoretical sense). For a quadratic Lagrangian, the inclusion of the Weyl tensor does not increase the number of DOF. In the case of a (generically nonlinear) Lagrangian f (R) depending only on the curvature scalar, the DOF reduce to three. According to common wisdom, a Lagrangian depending in a generic way on the Ricci tensor is equivalent to a second-order field theory including a massless spin-two field, a massive spin-two field and a massive scalar field. In the f (R) case the massive spin-two field is absent.
Apart from the f (R) case, nearly all theories with higher order terms in the Riemann tensor suffer from the presence of ghosts, i.e. fields with negative coefficients in front of the kinetic terms 1 . Such a situation is quite common in gauge theories where these fields are associated with gauge degrees of freedom. There are two approaches to this problem in quantum theory: either we gauge away these fields (what usually results in a much more complicated theory, since we explicitly break the symmetry to keep only positive kinetic terms), or we allow for quantization of these terms as well but later we decouple quantum states associated to these fields from the rest (Gupta-Bleuler quantization). The presence of ghosts is expected to lead to vacuum instability already at the classical level. However, it has been shown in [8] that for Lagrangian depending quadratically on the Ricci tensor the ground-state solution, where the massive spin-two field is set to zero, is linearly stable against small excitations of this field. The consequences of the ghost nature of the massive spin-two component of the multiplet are thus still an open question, but we shall not deal with this problem in this article; we shall only make a short comment on derivation of this property at the end of Sect. 4.
We shall instead focus on a different problem arising when one tries to identify NLG theories with second-order field theories involving a multiplet of fields, or -in a quantum perspective -to attach particle degrees of freedom to higher-derivative terms.
As we have already recalled, in standard field theory each fundamental field is expected to have a definite mass and spin: if definite values of these quantities cannot be assigned to the field under consideration, this simply means that the field is actually a unifying field (a kind of 'mixture') for a number of distinct fundamental fields. In the case of a NLG theory, this means that the metric appearing in its Lagrangian is just a unifying field and should be decomposed in a multiplet of gravitational fields of definite spin, that we have mentioned above. The physical metric of the spacetime will be only one element of the multiplet (and need not be identical to the original unifying metric). We have mentioned above this decomposition; for f (R) gravity theories the decomposition into a metric and a scalar field was performed in [9,10,11,12,13]. The case of the quadratic Lagrangian (without the Weyl tensor), L = κR + aR 2 + bR µν R µν , is also relatively well understood [8] (for a completely different approach to quadratic theory see [14]). From the physical viewpoint, the decomposition of the unifying metric into two or three component fields having definite masses and spins means establishing the particle content of the given NLG theory.
However, this decomposition does not ensure by itself that each field of the multiplet can be regarded as an independent physical entity. In ordinary field theories, in fact, it is (tacitly) assumed that physically distinct fields can, in principle, be excited independently from each other. We write "in principle", because if two fields are mutually interacting then an excitation of either of them indeed produces, in general, an excitation of the other: for instance, a charged matterfield cannot be excited without affecting the state of the electromagnetic field. Nevertheless, we can claim that the two fields are physically distinct -i.e. not merely two components of a single physical entitybecause, taking the appropriate limit of the model parameters (in the easiest situations, by setting some coupling constant equal to zero), we can produce a consistent version of the model in which both fields can live independently from each other.
When gravitation is involved the situation becomes more subtle, because any field entering the stress-energy tensor, when excited, may affect the spacetime metric and therefore alter the state of all other fields, even in the absence of direct coupling. Thus, as we discuss below in greater detail, in this case the physical content of the model should be investigated in a limit where both mutual interactions of fields and their interactions with the spacetime metric are suppressed (notice that we should not expect, in this limit, that all fields become free fields: there is no reason, in fact, to require that the self-interaction of each field be suppressed and in some cases it is not).
If it turns out that even in this limit it is impossible to get rid of constraints binding together the fundamental fields, then one should conclude that the theory cannot be consistently regarded as describing a multiplet of physically distinct fields.
It has been already observed in [8] that for a quadratic Lagrangian depending on the Ricci tensor the massive spin-two field and the scalar field cannot be excited independently. Here we address the question for a larger class of models, whereby a term nonlinearly depending on the quadratic invariant of the Weyl tensor is added to the Lagrangian. In section 2 we discuss the general setting of the problem.
NLG theories do not comprise an inherent method of decomposing the unifying metric into a gravitational multiplet of fields and in section 3 we discuss the best known decomposition procedure, based on a kind of generalized Legendre transformation. The procedure must be covariant and should not be confused with the (3+1) canonical formalism applied in general relativity; actually, we employ a far going generalization of known Legendre transformation to obtain second-order field equations (and in consequence our formalism cannot be applied to GR itself).
As a test of viability of the limit procedure, we apply it in section 4 to the quadratic Lagrangian L = κR + aR 2 + bR µν R µν . We stress that the procedure does not consist in expanding equations around the flat space solution: no truncation or linearization whatsoever is done. Instead, we push the three parameters entering the Lagrangian to a limit that decouples the two massive fields from the metric (this limit is singular for the Lagrangian, but not for the field equations). In this limit, the metric should be Ricciflat, as a consequence of the Einstein equation for the model. The (full) equations for the massive fields holding in this limit show that, although the two field together carry six DOF as expected, the scalar field cannot be excited independently of the massive spin-two field. Hence, one cannot say that this NLG model is equivalent to a theory where five DOF are carried by a massive spin-two field and one DOF is carried by a massive scalar field. We shall see that this procedure actually comprises performing successively two independent operations to the multiplet. The subsequent sections are devoted to Lagrangians explicitly depending on the Weyl tensor. For the fairly general form of the Weyl contribution that we are considering one can introduce, besides the tensorial momentum conjugate to the Weyl tensor itself, a scalar momentum which subsumes the nonlinear dependence on the quadratic invariant C α βµν C α βµν . The procedure also depends on whether the Lagrangian explicitly contains the Ricci tensor or not. Accordingly, sections 5 to 8 consider the corresponding four cases. The common outcome of these cases is that the presence of the conformal tensor makes the interpretation of the theory in terms of a multiplet of fields (particles) obscure or even impossible. The paper is focused on concepts, methods and results whereas all calculations are reduced to the essential minimum.
2 Physical fields, degrees of freedom and flat spacetime The first assumption that we make is that any tenable theory of gravitation is a specific field theory. This means that most concepts and notions developed in standard classical field theory in flat spacetime should adequately apply (except for energy!) to gravitation. This is the case of Einstein's GR which can be well approximated by a linear spin-two field in Minkowski space (there are also attempts to formulate GR as a nonlinear field in flat space). We therefore postulate that any metric gravity theory should obey the tenet of Lagrangian field theory (both classical and quantum): a fundamental field should have a definite (nonnegative) mass and spin. If a field appearing in the Lagrangian (and in this sense a fundamental one [15]) does not satisfy this condition, it should be decomposed into a multiplet of interacting fundamental fields. Yet if gravitation (NLG theory) does not obey this tenet, then very little can be asserted about its physical interpretation. Properties of fundamental fields (particularly masses and spins) must be determined in the coupling-free limit. By this expression we mean the following: in any field theory the fundamental fields can be interacting with each other, but if they are to be regarded as being physically distinguishable (and not just an artificial representation of different degrees of freedom of a single physical entity), then it should be possible to remove their mutual interaction, by taking a suitable limit of the coupling constants which appear in the Lagrangian.
In this limit, the fundamental fields should behave as independent test fields: in particular, it should be possible to excite any one of them while the others remain in their ground state. This should be true for the full model, and not only in the linear approximation. We remark that in the coupling-free limit each fundamental field can still have a nontrivial self-interaction: it is the mutual interaction between different fields that should disappear. As we shall see, in our case equations of motion for the fundamental fields become linear in this limit, however one should not assume it in advance.
If a system of interacting fields is considered in flat spacetime, studying the couplingfree limit is relatively easy. Suppose that fields φ 1 , φ 2 , . . . φ n interact. In many cases one gets φ 1 free just by assuming that all other fields are in their ground states, φ i = 0 or φ i are constant tensors for i = 2, . . . n. Then all equations of motion for the system (including possible constraints) are reduced to some equations for φ 1 and one determines the properties of the field in the standard way; next one repeats this procedure for all other fields of the system.
It may happen, however, that even after 'switching off' all interactions by canceling their corresponding coupling constants, some fields in the system, say φ 1 and φ 2 , are so strictly coupled that φ 1 cannot be excited if φ 2 is in its ground state. This may occur, for instance, due to a constraint relating φ 1 to φ 2 . We then say that φ 1 and φ 2 are not independent physical fields. We shall show that this case does occur in most NLG theories, particularly if their Lagrangian contains a Weyl tensor term. According to our viewpoint, in these cases the original fourth-order theory cannot be legitimately regarded as being equivalent to a physical second-order field theory: equivalence holds only provided the Lagrangian is restricted in such a way that the degrees of freedom are reduced to those carried by physically independent fields. We shall say that such restricted Lagrangian represents a physically tenable theory.
This procedure of identifying a dependent field may be equivalently expressed in terms of number of degrees of freedom (DOF). Counting of physical degrees of freedom in a theory with higher order terms in gravity is given for example in a seminal paper on conformal supergravity (i.e. with the Weyl squared term in the bosonic sector) [16]. The counting is quoted as an argument for the equal number of bosonic and fermionic degrees of freedom but it illustrates the main point advocated in the present paper for a general theory with higher order terms in the Riemann tensor. It should be emphasized that in supergravity the couplings cannot be arbitrarily put to 0 since the local supersymmetry provides strong constraints on such a procedure but still the counting works as described below.
The standard way of determining the number F of DOF carried by the system is by applying the Hamiltonian formalism and its power lies in that it not only counts the degrees of freedom, but that it first of all identifies them. However in the case of nonlinear field systems the formalism is intricate (it is more intricate for NLG theories, see [17], than the approach developed here) and for our purposes explicit identification of DOF is unnecessary; what we actually need is their number F . To this end one computes the number of independent initial data on a Cauchy surface, which is equal to the number of algebraically independent components of all the fields multiplied by the order of the Lagrange equations minus the number of all constraints imposed on these equations and following from them (as we recall in the next section, these include the conservation laws originating from coordinate invariance); F is one half of this number [18,19]. This approach is non-algorithmic (whereas the canonical Hamiltonian formalism is algorithmic, but in our case is too complex to be useful) because there is no general method of checking if all constraints have already been found; some constraints may be hard to recognize. For simple gravitational Lagrangians we study here this approach works effectively; moreover due to decomposition of gravitation into a multiplet of fields, some DOF may be precisely identified. Then one determines the number of DOF carried by each field of the system. To this end, as above, one first ensures that all the fields may be decoupled by canceling their interactions and this may require taking appropriate limits in all coupling constants. Next one takes one field of the system, assumes that all others are in their ground state and studies the resulting equation of motion for it. This equation together with possible constraints determines the number of DOF for the field. Repeating this procedure for all other fields one gets the total number of DOF of fundamental fields of the system. This analysis should be done in flat spacetime, for two reasons. Firstly, physical properties of ordinary fields -to which one should compare the properties of the gravitational multiplet -are currently formulated in the framework of special relativity, i.e. in Minkowski space. Secondly, gravitational interaction of a field (even if it is a test field) may increase the number of its DOF (usually it is believed that an interaction may decrease the number of DOF, yet curvature may do the opposite) [19]. This may occur for simple Lagrangians for vector fields since spacetime curvature may cancel some of the constraints satisfied by the field in flat spacetime. In any NLG theory the corresponding gravitational multiplet, whatever way it has been defined (we discuss this issue in the next section), is given together with its Lagrange equations in a curved spacetime and making a transition to the flat one is a subtle problem. In the framework of GR one recovers the standard physics as follows. Clearly one cannot merely replace g µν by Minkowski metric η µν in Einstein field equations (EFE) because they then imply T µν = 0. Physical fields have positive energy density and flat spacetime trivially implies that it is empty. Instead one must 'trivialize' EFE by breaking the coupling of any matter to geometry by taking the limit G → 0; this mathematically makes all matter a test one. Gravitation (curvature) is then free, but one cannot assume that the spacetime is any fixed solution of R µν = 0 or any other fixed metric. This is so because in a fixed but curved spacetime it is impossible to recover the standard physics, e.g. the integral energy and momentum conservation laws for continuous matter do not hold; moreover physics in a plane gravitational wave spacetime or in anti-de Sitter space is bizarre [20]. Gravitation may be also switched off in the framework of GR by two other methods in the full Lagrangian: i) by putting R = 0 and replacing g µν by η µν in the matter Lagrangian, ii) since the gravitational Lagrangian is c 3 R(16πG) −1 , one can take the limit 1/G → 0 and insert the flat metric. Clearly these two methods cannot be applied to NLG theories, where one is interested in the gravitational multiplet and introducing ordinary matter is a subtle problem (which in most papers is circumvented by an arbitrary assumption) and merely makes the picture more intricate.
We therefore investigate pure gravity with L = f (g µν , R α βµν ); this Lagrangian must depend on at least two dimensional constants (besides c). Take for instance an analytic function L = f (R) = n q n R n . Dimensionality of any Lagrangian is ML −2 T −1 (M-mass, L-length and T -time), hence dimensionality of the coefficients is [q n ] = ML 2n−2 T −1 = [c 3 /G]L 2n−2 and one needs one additional independent constant of dimensionality of length. Contributions of the Ricci and Weyl tensors to the Lagrangian may require further dimensional constants. All these coupling constants appear in the equations of motion for the gravitational multiplet of fields (these equations are equivalent to the fourth-order field Lagrange equations for the unifying metric g µν ). For both mathematical simplicity and physical interpretation the multiplet fields are so defined as to make their equations of motion second-order ones; in particular, the equations contain the Riemann tensor but are free of its derivatives and in this sense they are analogous to Einstein field equation (EFE) in GR. The gravitational multiplet consists of a spacetime metric (identical or not to the unifying g µν ) and a number of non-geometric fields 2 which act in these equations as sources for this metric. Crucial for the issue of tenability of the Lagrangian are properties of these non-geometric components of the multiplet and these features must be studied in flat spacetime, what means that one forgets their geometrical origin (as they arise from L depending solely on g µν ) and regards them as ordinary fields consistently defined in Minkowski space. As in the case of GR, one cannot simply insert flat metric into the equations of motion since this makes them trivial: the equations of motion "remember" the geometric origin of these fields and enforce them to remain in the ground state if the metric is flat. The transition to flat spacetime requires a procedure of breaking all couplings of these fields to the curvature. In GR the coupling is only in EFE and to break it, it is sufficient to put G → 0 whereas the form of matter propagation equations remains untouched, for example the electromagnetic field satisfies standard Maxwell equations. In a NLG theory, instead, the removal of the coupling affects the entire system of equations of motion (and some of them just disappear). This is accomplished by taking appropriate limits in all coupling constants. Furthermore one must assume that the constants tend to a limit at the same rate, so that their ratios remain finite in the limit.
To summarize, our method of determining the physical properties of a gravitational multiplet of fields 3 comprises two operations which are independent, consistent and unique in providing the required outcome if performed in the following order.
i) For the coupling constants of the multiplet fields (originally being essential dimensional parameters in the Lagrangian for the unifying metric) one takes the limits such that the multiplet fields cease to act as sources in all equations involving the Riemann tensor of the spacetime metric. In the model resulting from this limit, the multiplet fields (other than the metric itself) become devoid of geometrical interpretation and propagate as test fields in a fixed spacetime.
ii) The physical nature of these fields is investigated upon choosing, in particular, Minkowski spacetime (to avoid bizarre physical effects occurring e.g. in anti-de Sitter space). This amounts to merely identifying the fixed metric with the flat one and then studying equations of motion.
To avoid any misinterpretation we stress once more that the method is exact, there are no approximations.
Up to now these are qualitative principles: to make them precise one should introduce a concrete decomposition procedure and make some simplifying assumptions on the Lagrangian.

Decomposition, various generalized Legendre transformations and the Weyl tensor
The fourth-order equations arising from any NLG Lagrangian clearly show that one cannot assign a definite mass and spin to the metric appearing in it 4 , implying that it should be decomposed into a multiplet of gravitational fields. In the linear approximation it was first done by Stelle [7]. The method of decomposition in the exact theory is not inherent to it and should be selected in such a way as to provide the required outcome -a well defined collection of fields with definite masses and spins. The correctness of the method is verified by counting the degrees of freedom and this is closely related to showing the equivalence of the equations of motion for the multiplet to those for the unifying field. In other terms, first one should compute the number of DOF for the original NLG Lagrangian including the "unifying metric" alone (and its curvature components); we shall denote this number by F . After decomposing the model into field variables of definite spin and mass, one should find that the number of mathematical DOF carried by the whole multiplet (in curved space, counting the metric itself in the multiplet) coincides with F : this is a necessary condition for the correctness of the decomposition. Having verified this one should study individual fields of the multiplet, in flat space, to identify fields which are not physically independent (in the sense described in sect. 2). If such fields do exist, the underlying Lagrangian should be appropriately modified or rejected. By inspection of the Lagrange equations for a generic metric NLG theory with L = f (g µν , R αβµν ) one finds that the initial metric g µν carries F = 8 DOF (see [7,21,22,23] for quadratic Lagrangians, [24] for a perturbative method and [25] for a generic proof) and the above procedure is applied.
Here two comments on counting the number of DOF in this way are in order.
1. In general relativity, the metric carries two DOF; in vacuum, there are four constraint equations G 0 α = 0. If any matter is present these are replaced by G 0 α = T 0 α and in general the energy-momentum tensor does contain second time derivatives (of the metric and matter components) and these equations are no longer constraint ones. However if the matter is in the form of a scalar, vector or tensor (of second rank) field, it may be shown that the second time derivatives of the field may be removed from T 0 α with the aid of the Lagrange equations and the four constraints are restored [35]. A fully generic proof of that is lacking, nonetheless disregarding possible very exotic forms of matter being a counterexample, the metric has two DOF also in the presence of sources. The case of any NLG theory is more subtle. Due to the inherent ambiguity of these theories concerning their physical interpretation (i.e. which set, or 'frame', of field variables is physically measurable), it is unclear whether matter should be (minimally) coupled to the unifying metric or to the metric appearing in an appropriately redefined multiplet of fundamental fields (the latter choice is advocated in [12]). Depending on the choice, the matter source of gravity might be a priori different from T αβ (be a function of the tensor) and the required four constraints might disappear. Fortunately, it may be argued that this is not the case [12].
2. We emphasize that independently of the actual form of the Lagrangian, four equations among the Lagrange ones are always constraints. Let the fourth order field equations of any NLG theory be denoted E µν (g) = 0. The invariance of the corresponding action functional under any infinitesimal coordinate transformation gives rise to the strong Noether conservation law ('generalized Bianchi identity') ∇ β E β α ≡ 0 and this implies that the components E 0 α involve at most third time derivatives, then E 0 α = 0 are constraint equations.
To the best of our knowledge, the best (fully covariant) decomposition method derives from the idea, first introduced by J. Kijowski [26], of a "canonical formalism" where the metric is regarded as the "conjugate momentum" to the connection in a purely affine Lagrangian theory and the second-order form of the equations is obtained through some kind of generalized Legendre transformation. The idea to generalize the Legendre transformation in this direction is mathematically correct and was applied to physics in the book [27] and later appeared many times in the literature. The choice of this formalism is strongly supported by a theorem that the dynamics of the unifying metric with the generic Lagrangian is equivalent to equations of motion for the resulting multiplet of fields and these equations necessarily contain the Einstein field equations for (the original or a transformed) metric [4].
In the usual Legendre transformation, the first step consists in the definition of conjugate momenta as function of the "velocities" (i.e. the highest derivatives of the dynamical variables occurring in the Lagrangian). Then one should (at least locally) invert the mapping, so that all "velocities" in the Lagrangian can be replaced by suitable functions of the momenta.
However, if one starts from a covariant second-order Lagrangian and aims at a covariant decomposition, the role of "velocities"cannot be played by the second derivatives of the metric components. In fact, the corresponding "conjugate momenta" would not be tensors; moreover, it would be impossible to invert the relation between "velocities" and "momenta", because a covariant Lagrangian can depend on second derivatives of a metric only through the components of the Riemann curvature. Then, one should consider a more general setting and introduce the conjugate momenta to the curvature components.To this end we first briefly recall the procedure introduced in [28,29].
Consider a Lagrangian L = L(φ, ∂φ, ∂ 2 φ), where φ is a tensor field of some rank. L is assumed to depend on a set of k independent linear combinations ω (the Hessian matrix is of rank k), one defines momenta canonically conjugate to ω A , and these equations can be inverted (solved w.r.t. ω A ), ω A = Ω A (φ, ∂φ, π). Then one constructs an auxiliary function H playing here a role analogous to the Hamiltonian in the standard canonical formalism of a field theory (this function should not be confused with the physical Hamiltonian generating the time evolution of the fields), and a Helmholtz Lagrangian (originally introduced by Helmholtz in mechanics [30,31,32]) The dynamical variables are now φ and π A and L H is a specific Lagrangian depending on φ, ∂φ and π A , but not on ∂π A . The variation δπ A gives rise to second order equations whereas the fact that L H is linear in ω A (and in consequence in ∂ 2 φ) makes the equations δL H /δφ = 0 third order ones. The two sets of equations for L H are equivalent to fourth order Lagrange equations for the original L, showing that this generalization of a Legendre transformation makes sense. Nevertheless it raises the problem: in the standard Hamiltonian formalism the canonical momenta are derivatives of a Lagrangian with respect to time derivatives of dynamical variables and carry independent degrees of freedom. The conjugate momenta defined above comprise both first and second time and spatial derivatives, then can they be carriers of physical DOF, or are they merely mathematically convenient auxiliary fields?
We claim (and this becomes evident when the theory of gravitational multiplets is developed in the following sections) that the equations arising from the Helmholtz Lagrangian involve derivatives of these auxiliary fields that cannot be integrated to give algebraic equations. Therefore we have all the reasons to treat these conjugate momenta (except χ below) as bona fide physical degrees of freedom. The case of a generic NLG theory with L = f (g µν , R αβµν ) is more complicated than that above. The curvature tensor has altogether 14 independent scalar invariants and if one wishes to study NLG theories in full generality, then f is a generic (transcendental) smooth function of these 14 variables. It is clear that if a theory is to be of any practical use, its Lagrangian should be simple. Long ago Einstein expressed this idea as a principle stating that Lagrangians of verified theories are the simplest possible functions of the field variables. In the case of NLG theories this principle is broken by definition, nevertheless their Lagrangians should be possibly simple and our taskproving that the Weyl tensor should not be present in L -requires a proof for quite simple f , without entangling into a generic L. In the case of L = f (R) it was shown [12,25,33] using physical arguments that if f (R) is analytic at R = 0 (as it should be), then it should be of the form R + aR 2 + . . ., with a > 0. For this Lagrangian it is not easy to find out constraints following from the field equations and in order to determine the number of DOF one first applies a version of this canonical formalism to decompose the unifying metric into a doublet comprising the metric and a scalar; this system has three DOF. This system is dynamically consistent showing that no DOF have been lost due to the appropriate generalized Legendre transformation. In our opinion the fact that it is quite easy to count the DOF for the doublet shows the advantage of dealing with these two fundamental fields instead of the unifying metric.
If L explicitly depends on the Ricci tensor, one views it as an additive 'correction' to f (R): since any linear dependence on R µν is already incorporated in the term f (R), the simplest additional term containing the Ricci tensor is quadratic; the quadratic term R µν R µν , on the other hand, should necessarily be present to ensure that the Lagrangian be regular with respect to R µν around the ground state (flat) solution. The resulting Lagrangian L = κR + aR 2 + bR µν R µν was studied in [8]. It is clear that for this L the linear combinations ω A (∂ 2 g µν ) are R and the traceless S µν ≡ R µν − Rg µν /4.
Next, if one admits a further correction due to the Weyl tensor, this should be expected to depend on the simplest Weyl invariant, W ≡ 1 2 C αβµν C αβµν ; however, this invariant by itself does not provide a useful additional term, because (in dimension four) it can always be removed from the action integral by subtracting a full divergence (in accordance with the Gauss-Bonnet theorem). Even in higher dimensions, however, constraints arising from Bianchi identities are known to affect the DOF of a quadratic Lagrangian containing the Weyl tensor [29]). We shall therefore investigate L of the form where f is a nonlinear function (f ′′ (W ) = 0), and we shall argue that the term f (W ) is physically redundant. We claim that this simple Lagrangian is sufficient for our task and that it is unlikely that a more intricate dependence of L on R, R µν and C αβµν will make the Weyl tensor dependence physical. In general the term, f (W ) may be dealt with in two different ways: (i) either one defines a tensorial momentum conjugate to the linear combination of the second derivatives ω A = C α βµν (sections 5 and 7), (ii) or one defines a scalar momentum conjugate to W (sections 6 and 8). We shall show that none of these ways provides a satisfactory physical theory (independently of whether some solutions may fit the observational data).
The method referred as (ii) requires a preliminary explanation. The generalized Legendre transformation, described above in this section, relies on the assumption that the Lagrangian depends (in a regular way) on a linear combination of the highest (in our case, the second) derivatives of a dynamical variable. To our knowledge, the definition of a momentum conjugate to a nonlinear (e.g. quadratic) function of the highest derivatives has not been considered so far. Yet, it works: the procedure yields an equivalent Lagrangian, but there is a substantial difference w.r.to a genuine Legendre transformation: here, the highest derivatives do not disappear, but the dependence of the Lagrangian on these derivatives is "standardized". To clarify the procedure, let us show a simple (first-order) example from classical particle mechanics.
On a finite-dimensional configuration space Q, consider a holonomic system described by a Lagrangian of the form where K ≡ 1 2 g µνq µqν (the Lagrangian for an ordinary mechanical system corresponds to the case where f is linear). For simplicity, assume that the configuration space metric g µν is constant. The Euler-Lagrange equations are Let us now introduce the variable Provided f ′′ (K) = 0, this relation can be inverted, i.e. one can locally define a function S(p) = K such that f ′ (S(p)) ≡ p. Let us then introduce the Helmholtz Lagrangian where V (p) = f (S(p)) − p S(p). The Euler-Lagrange equations for L H are and are manifestly equivalent to (3), on account of f ′ (S(p)) = p.
We have not added a new degree of freedom, because L H is degenerate: it does not depend onṗ, and therefore the first of the E-L equations is actually a constraint. We have merely shown that the the original Lagrangian L, which depends in a nonlinear way on the kinetic energy K (and thus is not quadratic in the velocity componentsq λ ) is equivalent to a Lagrangian L H which is quadratic in the velocity components, but is written in an extended configuration space Q × R: the number of degrees of freedom is unchanged because all solutions in T (Q × R) entirely lie on the submanifold defined by the equations p = f ′ (K) andṗ = f ′′ (K)K, which has the same dimension as T Q.
The Helmholtz Lagrangian is now a function on T (T * Q). Apparently, the configuration space Q has been enlarged to T * Q, but the Lagrangian is now doubly degenerate: not only it does not depend onṗ µ , but it is also linear inq µ . Therefore, the Euler-Lagrange equations are first-order (instead of second-order) in both q λ and p λ . These equations define a vector field on T * Q, while the E-L equations for the original Lagrangian L (assumed to be hyperregular) define a vector field on T Q. In other words, after the Legendre transformation, the momenta p λ have completely replaced the velocity componentsq λ as coordinates in the space of physical states of the system. On the contrary, in the Helmholtz Lagrangian (5) the velocitiesq λ are still present, but the new Lagrangian depends linearly on K: the original nonlinear dependence f (K) in (2) has been "absorbed" by the auxiliary (non-dynamical) variable p; the coupling between p and K in (5) is "universal", and the only vestige of the original form of the Lagrangian is the "potential" V (p) which generates a constraint.
In sections 6 and 8 we exploit this method to show that the scalar momentum conjugate to W contributes to a mathematically correct dynamical description of the corresponding gravitational multiplet. It is a purely physical argument, rather than a mathematical one, that shows that the momentum conjugate to W should be rejected. In conclusion we state that the formalism developed in this work is applied to prove that if a metric nonlinear gravity theory such as that investigated here is not equivalent to GR plus a number of fields which represent non-geometric, physically distinct DOF, then the theory is not tenable from the field-theoretical viewpoint. In precisely this sense we shortly say that it is unphysical. And we emphasize that contrary to the common belief this equivalence does not hold for a generic gravitational Lagrangian.

Lagrangian quadratically depending on R and the Ricci tensor
This case was discussed in detail in [8] and here we present a modified version of those calculations focused on properties of the resulting gravitational triplet. There is no matter and the Lagrangian contains three coupling constants, Here Positivity of energy in Einstein frame for L = f (R) gravity requires a > 0 [12], whereas it is suggested in [8]  The fourth-order field equations were first (?) derived in [7] and then in many other works (see e.g. [8,34]); clearly they are equivalent to second-order equations derived below. One introduces one scalar and one tensor momentum, hence π µν g µν = 0. The fields χ and π µν are of the same dimensionality equal to [χ] = [π µν ] = [κ] = MT −1 . For convenience we set c = 1. In this way the original metric gets decomposed into a triplet forming the Helmholtz-Jordan frame (HJF) {g µν , π µν , χ}. The metric remains unchanged (the signature is −+++), only its dynamics will be described by equations of motion of different form. The function H (pseudo-Hamiltonian) is It generates the Helmholtz Lagrangian The Lagrange equations for the momenta χ and π µν are algebraic and recover the definitions of the fields, eq. (15) confirms that π µν is traceless. Equations (14) and (15) are together composed into quasi-Einsteinian field equations, and these give rise to four differential constraints (due to Bianchi identities) The metric variation of L H generates EFE G µν = 8πGT µν (g, χ, π) comprising an energymomentum tensor for the two momenta. One sees comparing (16) and EFE that for solutions these equations provide two expressions for T µν which may be used both as equations for the metric and for the momenta χ and π µν . Dividing (16) by 8πG one gets First the trace of (18) generates an equation of motion for χ. Applying the divergence to eq. (17) and inserting the resulting expression into this equation one gets For 3a + b > 0 the mass of χ is real, m 2 χ = [32πG(3a + b)] −1 . If 3a + b = 0 the scalar χ is constrained to vanish [8]. Now assuming 3a + b > 0 for a > 0 and b < 0 (negative b gives rise to real mass for π µν in the case 3a + b = 0 [8]) one can eliminate χ from (18). To this end one again takes the divergence of (17) which reads ;αβ (20) and applies it to eliminate ✷χ from (19) and finally arrives at This shows that χ is not an independent physical field, being completely determined by values of the double divergence of π µν . One would be lead to regard the scalar momentum as being spurious and remove it by the appropriate restriction of Lagrangian (8), 3a + b = 0. However, (21) was obtained using (16), so this relationship between χ and π µν might be an outcome of their interaction with the metric. As we have discussed in the previous section, to prove that this is not the case we need to envisage the behavior of these fields in the decoupling limit in flat spacetime. Before taking this limit we exhibit the complete dynamical structure of the model. To this end one replaces in eqs. (18) χ by (21), χ ;µν by the derivative of (17) and ✷χ by (20) and arrives at Notice that the symbol T µν is shown here merely to indicate the origin of this equation, but is otherwise irrelevant. After elimination of the scalar the system of equations of motion consists of: -9 quasi-linear eqs. (22) for π µν (their trace vanishes identically), -EFE eqs. (16) for g µν , -expression (21) for the auxiliary scalar, to be inserted into (16). These equations do not directly follow from a Lagrangian according to the standard formalism, instead one has to perform the above more involved procedure. The fact that the procedure generates two different expressions for T µν (g, χ, π) which are equal for solutions, plays essential role in it. These equations are independent. In fact, (16) successively generate (17) and (20), but one cannot derive (21) from them and from (22). Yet applying (21) to (16) and (22) one can invert the whole procedure. The fact that equation (19) can be integrated to yield (21) shows that the initial data for χ are determined by initial data for π µν and the scalar carries no its own degrees of freedom. It is likely that the possibility of integrating of (19) is due to the fact that the Lagrangian is quadratic in R and R µν ; in presence of higher order terms an analogous equation for χ is not expected to be integrable to an algebraic expression. Yet it is unclear at this point whether an algebraic relationship such as (21) does imply that the number of available DOF for the doublet {g µν , π µν } is decreased from 8 to 7. This issue must be determined in the decoupling limit in flat spacetime and we shall see that there are still 8 DOF (including the metric).
To get the physical interpretation we now consider the gravity-free system in flat spacetime. Since the system of equations (16), (22) and (21) is equivalent to (16), (17), (18) and (20), we take the decoupling limit in the first system and parallelly observe its outcomes in the latter system. 1. By analogy to taking the limit in GR coupled to ordinary matter, we first decouple all physical fields from the spacetime metric by putting G → 0 in EFE G µν = 8πGT µν . Here T µν is the variational energy-momentum tensor generated by L H ; in principle it might include the contribution from ordinary matter. In this limit the spacetime is a fixed solution to G µν = 0 and T µν describes test fields in it. 2. Vanishing of G µν implies that RHS of (16) vanishes for arbitrary π µν and χ(π). This is possible only in the limit 1/b → 0 ⇔ b → −∞. Formally it is sufficient and the constant a might remain finite, |a| < ∞, but it would break the condition 3a + b > 0. This shows that also a must go to infinity and at the same rate as b. Let a = −ξb for dimensionless ξ, then ξ > 1/3 and let ξ remain constant for b → −∞. In this limit eqs. (14) and (15) reduce to R(g) = 0 and S µν (g) = 0. 3. To take the limit of G, a and b in formula (21) (and respectively in (19)) one assumes that G and 1/b tend to 0 at the same rate, so that the product Gb is finite and negative in this limit. Let κ/b → −λ −2 < 0 with [λ] = L. One finds in the limit that eq. (19) reads ✷χ = χ 2(3ξ − 1)λ 2 (23) and formula (21) takes the form showing that also in the gravity-free system the scalar is not independent and is just a name for the double divergence of π µν . For further use we show the constraint (17) which reads now π µν 4. Finally one takes the limit of the three coupling constants in eq. (22). After some manipulations one arrives at The whole system gets reduced to a test linear massive field π µν subject to eqs. (26), whereas one preserves eqs. (23), (24) and (25) as auxiliary formulae to be applied later on.
Now one performs the second procedure of the two fundamental ones described in Sect.2, namely one studies the dynamics in flat spacetime. In this spacetime one may determine the physical properties of π µν by solving the above system of equations in terms of arbitrary initial data. In a chosen inertial reference frame and using Cartesian coordinates one gives initial data at t = 0: π µν (t = 0) = f µν (x) and ∂ 0 π µν (t = 0) = h µν (x). These are together 18 functions since π µν η µν = 0 = π µν,0 η µν imply f 00 = j f jj and h 00 = j h jj for i, j = 1, 2, 3. To show that these Cauchy data uniquely determine a solution to eqs. (26) one first solves eq. (23) for χ as if it were an independent variable, applying (25). The initial data define a scalar and a 3-vector at t = 0, these are known functions. The constraints (25) determine initial data for χ, One sees that K i must be a gradient and the functions f ij and h 0i are subject to K i,j = K j,i . If these constraints hold, the components K i uniquely determine (up to an additive constant) the function K(x) such that K i = ∂ i K. Then the initial value of χ is χ(t = 0) = (1 − 4ξ)K(x) and the integration constant is eliminated by imposing appropriate boundary conditions at infinity on h 0i , f ij and χ(0). The Cauchy data χ(t = 0) and χ ,0 uniquely determine the solution χ(x µ ) of (23). Now one returns to eqs. (26) and employing (24) and (25) one finds that in flat space they read Although χ fully depends on π µν , we have seen that the solution for χ can be fully determined from initial data, without knowing the solution for π µν : doing so, the term τ µν becomes an explicitly known source for unknown π µν . Physically this is a bizarre situation, yet mathematically there is nothing inconsistent here; (29) is a kind of tensorial Klein-Gordon equation and its solution is uniquely determined by the initial data.
Having this knowledge one counts the degrees of freedom for π µν . The constraints K i,j = K j,i reduce the number of independent functions K i to one, K(x), hence the number of independent data is diminished by 2. Then there are dynamical constraints since not all of E µν = 0 are hyperbolic propagation equations. From (25) one finds that hence these equations are constraints ones. However E 00 = 0 is not an independent equation, since η µν E µν ≡ 0 yields E 00 = j E jj and the independent dynamical constraints are E i0 = 0 and j E jj = 0. Together the number of independent initial data is 18 − 2 − 4 = 12 corresponding to six DOF for π µν . Yet a massive quantum spin-two particle has five DOF [36,37,38,39,40] indicating that the hypothetical particles of the quantized field π µν do not have a definite spin. This is unphysical. In the case 3a + b = 0 the scalar momentum is eliminated by its equation of motion and π µν is subject to the tensorial Klein-Gordon equations with the mass m π = 1/λ and constraints η µν π µν = 0 = ∂ ν π µν in flat space. This field carries 5 DOF [40] and the doublet {g µν , π µν } has 7 DOF (since for π µν = 0 the metric satisfies G µν = 0 and carries two DOF). This agrees with the fact that the unifying metric has seven DOF: the original fourth order equations, beyond the standard constraints giving rise to 8 DOF, satisfy in this case two additional constraints, R = 0 and ∂R/∂t = 0, which diminish the number to seven. For the Lagrangian (8) we previously assumed in [8], just for convenience, that 3a + b = 0, then we got dynamically χ = 0 and 7 DOF for the doublet. However, the choice 3a + b = 0 is not a matter of convenience, we emphasize that it is required by physics: though the scalar χ inevitably appears at an intermediate step in the Legendre transformation formalism, it is not only spurious as being determined by π µν , moreover it must vanish (due to equations of motion), otherwise one cannot assign a definite spin to π µν .
For the sake of completeness we make a comment on physical properties of π µν in this case. Eq. (16) shows that the field is not fully satisfactory physically: for solutions its energy-momentum tensor is linear, hence its energy density is indefinite. Its properties were investigated in [8]. Its ground state solution π µν = 0 is at least linearly stable (the mass is positive if b > 0); the problem of nonlinear stability (the exact nonlinear π µν interacts with the metric) is hard and open. Whether or not π µν is a ghost field cannot be established in the exact theory since its (Helmholtz) Lagrangian (13) does not contain at all derivative terms. The field turns out a ghost in a linear approximation [8] and this is proven in a rather intricate way. Thus π µν is not a typical ghost. Furthermore there are arguments [41] that being a ghost is not so disastrous as it was formerly believed.
Conclusion: the gravitational theory based on (8) is equivalent to GR interacting with a massive spin-2 field only in the case 3a + b = 0 and only in this case may be tenable from the field-theoretical viewpoint.

Tensorial momentum conjugate to the Weyl tensor
From the fact that both L = f (g µν , R αβµν ) and L = R + aR 2 + bR µν R µν (prior to excluding the scalar) have 8 DOF (in the mathematical sense) one would be tempted to conclude that the presence of the Weyl tensor is redundant. Yet, we have just shown that the latter Lagrangian cannot represent a physical field theory with 8 DOF, so one may wonder if adding to (8) (with 3a + b = 0) a term involving the conformal tensor would restore the expected 8 DOF.
We thus apply the approach developed in previous sections to a Lagrangian of the form (1), with 3a + b = 0: . We assume that f (W ) contains no dimensional constants and that f (W ) is analytic at W = 0, f (0) = 0 = f ′ (0) and f ′′ (W ) = 0.
In the sequel, to avoid making the computations pointlessly cumbersome we assume that f (W ) is a simple function, say with n some even positive integer: this allows us to invert f ′ (W ) explicitly and thus to spell out all terms in the subsequent computations. The theory for the momentum σ α βµν should be consistent for any such n.
The "traditional" strategy that we have adopted so far consisted in introducing conjugate momenta which carry definite spin. This suggests to introduce again the scalar and the tensor momenta conjugate to the Ricci tensor, χ and π µν respectively, as in (10) and (11) (the unwanted scalar momentum may be eliminated only at the level of equations of motion), and in addition a new tensorial conjugate momentum corresponding to the Weyl tensor [29]: All momenta have dimensionality as MT −1 . Following from their definitions, π µν is traceless and σ α βµν has all symmetries of the conformal tensor: thus, at first sight, σ α βµν would be expected to correspond to a spin-two field as well. Yet, we know that the model has at most 8 DOF, so it is evident that there is no place for two independent spin-two fields; however, we need to understand how π µν and σ α βµν are constrained to each other (in the coupling-free limit where they should behave as test fields, if they were physically independent) to assess whether non-independent components can be identified and eliminated, to obtain -if possible -a viable covariant field theory which may saturate the 8 mathematically available DOF.
To invert the relationship (33) one introduces the square of σ, Z ≡ 1 2 σ αβµν σ αβµν and finds that it equals We denote the inverse of the function (34) The function H is (36) and the Helmholtz Lagrangian reads L H does not contain derivatives of the momenta, hence the Lagrange eqs. for them are algebraic and recover their definitions, and (35). Eqs. (38) combine together in quasi-Einsteinian field equations, and Bianchi identities imply 4 differential constraints π µν ;ν + 3χ ;µ = 0.
While taking variations w.r.t. metric g µν one assumes, as always, that the components σ α βµν , as they are defined in (33), are independent of the metric, δ g σ α βµν ≡ 0. After a number of manipulations one arrives at the following system of equations of motion, again having the form of EFE. If we assume (32), these read All the terms in (41) explicitly depending on the three momenta should be interpreted as a collective energy-momentum tensor for them. In practice RHS of (39) is an effective total energy-momentum tensor 8πGT µν expressed in terms of solutions. The trace of equations (41) determines the scalar χ as a function of the scalar Z(σ), hence the scalar momentum is eliminated from the field equations (if f (W ) = 0, χ ≡ 0). We emphasize that this elimination is possible due to the presence of two expressions for the energy-momentum tensor T µν , which are equal for solutions. One can therefore employ (41) as propagation equations for the remaining two momenta, whereas eqs. (39) take the role of Einstein equations for the metric. (Whether or not this stress tensor has positive definite energy density may be determined if the space of solutions is known.) In the propagation equations one should replace G µν according to (39). However eqs. (35) form another set of equations for the metric where the momentum σ is a source for the Weyl tensor C αβµν (g). Hence the full Riemann tensor is expressed as an algebraic function of the momenta. This shows that the formulated here second-order dynamics is different from the standard GR dynamics comprising EFE and equations of motion for interacting matter fields. Equations (35) are consistent with (39) provided C αβµν satisfies the standard differential identities for the conformal tensor. The first order identity in four dimensions, upon employing (35), (38) and (32) reads Since the scalar χ is merely a function of the other fields, the degrees of freedom are associated with the gravitational triplet of fields {g µν , π µν , σ αµνβ } and are subject to the following equations of motion: -10 EFE (39), equivalent to (38), -10 eqs. (35) which are not viewed as propagation equations for the metric and instead are interpreted as generating equations for necessary constraints imposed on the momenta, -the constraints (43), -9 differential propagation eqs. (41) for 9 momenta π µν and 10 momenta σ αµνβ (one equation, the trace, has been used up to provide an algebraic expression eliminating the scalar χ), -expression (42) for χ.
Notice that for these equations there is no limit (such as n → 0) which would reduce this system to that studied in sect. 4. This is a determined system of equations of motion which are equivalent to the 10 fourth order equations for the metric g µν directly following from the Lagrangian (31). Mathematically this is OK, yet this is inconsistent with field theory in flat spacetime requiring that each physical field, when its couplings to other fields are switched off, should be determined by a propagation equation (possibly containing self-interaction terms). In fact, to elucidate the issue, one decouples the two momenta from the metric, then it turns out that the decoupling makes the spacetime flat. As we shall see below it follows that 9 + 10 fields π µν and σ αµνβ are underdetermined.
1. First, one decouples all the momenta from the curvature by putting G → 0 in eqs. (41) written in the form G µν − 8πGT µν = 0 and gets G µν = 0.
3. The momentum σ αµνβ is fully decoupled from the curvature if k → 0 in eq. (35). Then C αµνβ = 0 and the spacetime is flat, R αµνβ = 0. In general relativity, if matter does not gravitate, R µν = 0 and gravitational waves may exist. Here one sees that inclusion of the conformal tensor in the Lagrangian as an independent field in gravitational multiplet has very restrictive consequences.

5.
To proceed further one assumes that G, m 2 and k vanish at the same rate, so that both m 2 κ = m 2 /(16πG) ≡ λ −2 and m 2 /k ≡ µ remain constant and positive in this limit; λ has the dimension of length.
6. Under this assumption the constant factor in RHS of (43) vanishes in the limit, 2n−1 → 0 and these constraints are reduced to 7. Finally, eqs. (41) as the propagation ones, are reduced in these limits to The fields π and σ are subject to only nine equations (45) and to the constraints π µν ;ν = 0 and (44).
We now check if in this limit, where all gravitational couplings have been removed, the two fields can exist independently. If we set σ αµνβ = 0, constraints (44) vanish identically. The field π µν is subject to the standard tensorial Klein-Gordon equations and to the same constraints as in the case where the term f (W ) is absent in the Lagrangian: hence, it has the same properties. If we instead set π µν = 0, the tensorial momentum is subject to σ αµνβ ;αβ = 0 (46) and to the constraints (44). It is well known [42] that if a massless field H αµνβ has all algebraic symmetries of the conformal curvature tensor and satisfies linear field equations H αµνβ ;β = 0, then it has spin two. Therefore, eq.(46) would be compatible with σ αµνβ being a massless spin-two field, however equations (46) and (44) are inconsistent. In fact, take divergence ∇ α of (44) and symmetrize the resulting equation in βµ, then (46) may be applied. One arrives at a system of 9 nonlinear second order equations, These depend on the power n in f (W ) whereas solutions of (46) are n-independent.
Consistency of solutions to (47) with those to (46) requires the former be n-independent and this implies that the n-dependent term in (47) must vanish, then (47) is reduced to Z ;αβ σ αµνβ = 0.
Furthermore, multiplying (44) by Z ;α and applying (48) one gets a new constraint, Z ;α σ αβµν ;ν = 0. This constraint, together with the constraint (48) and equations (46) and (49), form a system which is consistent only for σ αµνβ = 0. Clearly this inconsistency is due to truncating eqs. (41) and (43) and invalidating eq. (35). The inconsistency shows that the momentum σ α βµν , unlike the momentum π µν , does not exist as an autonomous physical field in flat spacetime. It makes sense only as an auxiliary notion derived from the conformal curvature.
We conclude this section by stating that the description of the gravitational system (31) in terms of the momenta π µν and σ α βµν is defective. There are too few propagation equations for the momenta and the definition of σ gives rise to constraints which exclude nonzero σ in flat spacetime. In other terms NLG theory (31) is not equivalent to Einstein's gravity generated by fields π µν and σ α βµν carrying real masses and definite spins. Degrees of freedom associated with f (W ), if there are any, cannot be expressed by σ α βµν .
To investigate of whether π µν and σ may independently live in flat spacetime, one assumes, as in section 5, that G, m 2 and k tend to 0 at the same rate, so that in this limit their ratios m 2 κ ≡ λ −2 and m 2 /k ≡ µ are finite and positive. 1. Tensors t µν and τ µν get decoupled from the curvature for G = 0 in EFE. This implies G µν = 0. 2. Equations (39) are consistent with G µν = 0 if m 2 = 0. 3. The Weyl scalar W must be independent of the scalar momentum σ = 0 and for f (W ) = W n /n this is possible iff k = 0 in W = v(σ) given in (51). Then the spacetime satisfies R µν = 0 and C αβµν C αβµν = 0. This comprises Minkowski space, any planeparallel (p-p) gravitational wave and few other spacetimes. 4. For f (W ) as in (32), RHS of (59) is proportional to a positive power of k and vanishes for k = 0, hence χ = 0. In consequence, eqs. (62) are reduced to these are traceless. The constraints (40) are simplified to π µν ;ν = 0 (since G/m → 0 in (60)). Now the system comprises (63) and the constraints and is underdetermined since there are only 9 eqs. (63) for 10 functions π µν and σ. In flat spacetime (as this case is crucial for the physical interpretation) the scalar σ disappears from the field equations and it is clear that it has no physical existence. Yet π µν is the standard massive spin-2 field. Just for pure curiosity one may study the dynamics of σ in a plane gravitational wave for π µν = 0. It turns out that C αµνβ σ ;αβ = 0 is reduced in this spacetime to one physically bizarre propagation equation which does not uniquely determine σ for given initial data.
We conclude that the momentum σ cannot be interpreted as a classical counterpart of a quantum particle and should not be introduced into the theory. 7 The tensorial momentum in the case where the Ricci tensor is absent The discussion above has shown that the Ricci term contribution to the Lagrangian consumes all degrees of freedom available to nongeometric components of a gravitational multiplet leaving no space for the Weyl tensor contribution. We therefore investigate now the case where only R and the conformal tensor are present in a gravitational Lagrangian and assume The aR 2 term with a > 0 is necessary on physical grounds [12]. We shall also usually assume f (W ) = W n /n for integer n > 1. Again the momentum associated to f (W ) may be either tensorial or scalar one and in this section we study the tensorial case. One defines χ and σ α βµν according to (10) and (33) respectively, introduces the scalar Z as in section 5, then Z(W ) is given in (34) and the latter is inverted to W = v(Z) giving rise to (35) for C αβµν . The pseudo-Hamiltonian is and the Helmholtz Lagrangian reads The Lagrange equations read χ = 2aR and (35) whereas for the metric one finds where the potential term is defined as it is significant that k is here in a positive power. One sees a substantial difference to the case where the Ricci term is present: there is no analogue to quasi-Einsteinian equations (39) and one cannot eliminate G µν from (67). The trace of eqs. (67) provides a propagation equation for χ, χ is an independent variable and also cannot be eliminated. The triplet {g µν , χ, σ αβµν } is subject to χ = 2aR, ten eqs. (35), now interpreted as propagation equations for the metric with field σ α βµν as a source, nine eqs. (67) and one equation (69) for χ. There are no constraints.
To establish whether χ and σ are classical counterparts of quantum particles one considers the case of test fields by putting G = 0 = k and a −1 = 0. The first three steps in the sequence are the same as in section 5 and one gets R αβµν = 0 and V (Z) = 0. Then assuming that both G and 1/a vanish at the same rate, so that κ/a → λ −2 > 0, one finds that (69) is simplified to and χ is the standard massive scalar field. Yet the lack of equations analogous to (39) causes that in the flat space limit equations (67) do not form equations of motion for σ αµνβ . In a curved spacetime eqs. (67), written in the form of EFE, G µν = 8πGT µν (g, χ, σ), define the energy-momentum tensor subject to ∇ ν T µν = 0 and these are four equations equivalent to a subset of the full system of equations of motion for the two momenta. In flat spacetime the tensor reads (after applying (70)) and is subject to ∂ ν T µν = 0. It is easy to show that ∂ ν T µν is identically zero and hence gives rise to no equations of motion for σ αµνβ , which is completely arbitrary. The description of the system (64) in terms of the tensorial momentum conjugate to the Weyl tensor is defective.

The scalar momentum conjugate to the Weyl invariant if there is no Ricci tensor
Finally we envisage the scalar momentum conjugate to the Weyl tensor contribution to the gravitational system (64). Now one introduces two scalar momenta, χ as in (10) and σ according to (50); the inverse to (50) is W = v(σ) as in (51). In analogy to (52) and (53) one computes ) and The Lagrange equations comprise χ = 2aR, W (g) = v(σ) and and τ µν (g, σ) is given in (55) and the trace g µν τ µν = 0 is shown in (56) and (57). Tensor τ µν involves fourth order derivatives of the metric via C α(µν)β ;αβ which cannot be expressed in terms of derivatives of the momenta. This means that the alleged Einstein field equations (73) are inherently of fourth order.
Again the trace of (73) results in a propagation equation for χ, it is quasi-linear and contains a source. The full system of equations of motion consists of R = χ/(2a), W (g) = v(σ), one eq. (75) for χ and nine eqs. (73) for 11 components of g µν and σ; there are no constraints. Next one derives the flat spacetime limit for G → 0, a −1 → 0 and k → 0; this makes sense provided κ a → λ −2 > 0. In this way one arrives at R = 0 = G µν and 2W = C αβµν C αβµν = 0, as in section 6. The fact that t µν (η, χ) is linear indicates that both the metric and χ require a redefinition [12,25,33]. Furthermore, one gets the Klein-Gordon equation (70) for χ, whereas for σ one may seek for four equations arising from ∇ ν (t µν + τ µν ) = 0. In flat spacetime τ µν = 0 showing that σ carries no energy and has no equation of motion. Yet ∂ ν t µν ≡ 0 when (70) holds.

Summary and conclusions
It is natural to expect that gravitational interactions are akin to all known physical fields (fundamental fields of the Standard Particle Model) which on the quantum level are interpreted as particles with definite masses and spins. Even without an explicit quantization, the particle interpretation should hold for gravitation. The particle interpretation makes sense if each particle species (actually a classical field) may individually exist when its gravitational interactions, as well as the interaction with other particle species, are switched off (by taking a suitable limit of the coupling constants) in Minkowski space.
In our discussion, we focused on theories of gravity which are envisaged as possible variants to general relativity: we write "variants", rather than "alternatives", to mean that such models share with general relativity the fundamental postulates -in particular, the idea of describing gravitation in terms of metric curvature. In general relativity, one expects that the "particle spectrum" of gravitational interactions reduces to a single massless spin-two particle, the graviton, corresponding to excitations of the spacetime metric around the ground state. For metric nonlinear gravity theories, where one gets fourth-order equations for the metric, one is led to regard the metric itself (with its higher derivatives) as a unifying field, and gravitation should be described by a multiplet of fundamental fields, each with definite mass and spin. According to common wisdom earlier found in the lowest-order approximation to a nonlinear Lagrangian, these include, besides the graviton field (the physical spacetime metric), a massive scalar field and a massive spin-two field (the latter is a ghost, since its kinetic term has negative sign in its quadratic Lagrangian). Here we apply a different method, which is exact (no approximations) and again get the graviton and two massive fields, a scalar and a spintwo one (the latter may be shown to coincide dynamically, in a linear approximation, with the ghost field mentioned above, though the issue is subtle). Thus gravitation in these theories is not interpreted as pure curvature, since the gravitational multiplet comprises, besides the metric, other fields that should not be viewed as determined by the geometry. Most results in this direction concern Lagrangians which depend quadratically on the Ricci tensor. We explored the consequences of a possible Weyl tensor dependence, starting from the following assumptions: • each of the fundamental fields should carry a definite number of degrees of freedom and should correspond to a field which in flat spacetime, once gravitational interactions are switched off, can exist independently of the other ones (i.e. can be excited while other fields remain in their ground states); • to obtain (for the full nonlinear model) a proper decomposition of the dynamical variables into a multiplet of fields, one should not rely on ad hoc tricks: instead, one should rely on a method which can be formulated independently of any particular Lagrangian. We identify this method with a generalized Legendre transformation, following and extending Kijowski's original proposal [26] (although the version we exploit here is based on the variational formalism, not on the introduction of a symplectic structure). Still, for Lagrangians depending on curvature, the application of the method is not unique and one should consider different possibilities: the correctness of the choice is verified by counting degrees of freedom for both the unifying metric and the multiplet fields. (Once the multiplet fields have been generated from the unifying metric, then they may be subject to further redefinitions which, however, do not alter their masses and spins. The freedom of the redefinitions makes any experimental test of the theory ambiguous: it is unclear which set -or "frame" -of multiplet fields is actually measurable; thus, we stress that here we consider gravity theories from purely field-theoretical viewpoint and do not take into account observational tests); • we restrict our investigation to Lagrangians which are reasonably simple: in particular, we consider Lagrangians depending on linear and quadratic invariants of the metric curvature. It is known that if the Lagrangian depends linearly on both the square of the Ricci tensor and the square of the Weyl tensor, then the latter can be eliminated by subtracting a full divergence: therefore we allow the square of the Weyl tensor to enter the Lagrangian through a nonlinear term, added to the usual generic quadratic Lagrangian.
Under these assumptions we have shown that whatever generalized Legendre transformation is used to determine the fundamental field associated to the Weyl contribution, the resulting field has no independent existence in flat spacetime and carries no degrees of freedom. We conclude that within this framework the Weyl tensor does not contribute to the multiplet of gravitational fields. Only R and R µν may contribute and the number of degrees of freedom is either seven if the traceless part of Ricci is explicitly present, or three if it is not (and the R-dependence is nonlinear). The form of the Lagrangian cannot be fully general, since the requirement of forming a particle multiplet imposes strict restrictions on both R and R αβ dependence. The mathematically allowed maximal number of 8 DOF is physically unattainable for a second-order field-theoretical model satisfying the above assumptions. In other terms, we have shown that contrary to the common wisdom mentioned above, a generic nonlinear metric Lagrangian does not describe gravity which is equivalent to Einstein's general relativity comprising as a source of the metric a number of matter fields (in the sense of being non-geometric quantities) with real masses and definite spins. The equivalence only holds for specific Lagrangians and in particular it does not hold if there is any explicit Weyl tensor contribution.