Conformal invariance versus Weyl invariance

The most general Lagrangian describing spin 2 particles in flat spacetime and containing operators up to (mass) dimension 6 is carefully analyzed, determining the precise conditions for it to be invariant under linearized (transverse) diffeomorphisms, linearized Weyl rescalings, and conformal transformations.


Introduction
Particle physics interactions, when considered at very high energy (probing then smaller and smaller distances), are expected to be independent of the individual masses of the particles themselves, which are negligible in comparia e-mail: enrique.alvarez@uam.es b e-mail: jesusanero@gmail.com c e-mail: raquel.santosg@uam.es son with the energy scale. Some sort of scale invariance is expected to be at work there. The same thing happens in second order phase transitions. The correlation length diverges, and again, scale invariance is at work. In fact in many known cases (in all unitary theories 1 ) this symmetry is upgraded to full conformal symmetry [2]. There is however a caveat. In quantum field theory the well-known need to renormalize the bare quantities implies that an arbitary mass must be introduced. This is the origin of the dependence of coupling constants with the energy scale, encoded in the correspondiong beta-functions.
It is however only recently that the precise relationship between scale invariance, conformal invariance and Weyl invariance has been clarified (cf. [3,4] and references therein). This includes the precise conditions for scale invariant theories to become conformal invariant and also the existence of the so-called a-theorem for renormalizable theories (cf. [5] and [6] for a recent review). Most of the work done so far has been in flat spacetime, where the gravitational field is absent, or at most, non-dynamic.
When such a gravitational field is present [7] there are two possible generalizations of scale invariance. The most direct of those is the algebra of conformal Killing vector fields (CKV), that is, those that obey L(ξ )g μν = φ(x)g μν (1) the fact that implies that the set of all CKV generate an algebra, which for Minkowski spacetime is the conformal algebra, SO (2, n). In fact, the maximal possible dimension of the conformal algebra is precisely d = (n + 1)(n + 2) 2 (3) which is attained by conformally flat spacetimes (the ones with vanishing Weyl tensor [8,9]). Unfortunately, however, this property is not generic; that is, an arbitrary metric does not support any CKV, and the corresponding algebra has to be studied for each particular spacetime by itself. The next most natural symmetry to study is Weyl invariance, the invariance of the action under local rescalings of the metric tensor.
this invariance, besides, can still be studied in the linear limit, when the gravitational field manifests itself as a perturbation of the Minkowski metric.
g μν (x) = η μν + κh μν (x) (5) Given that the previous expansion is exact, the linearized Weyl symmetry of the metric perturbation can be written as κ δh μν = 2ω(x)η μν + 2κ ω L ξ h μν (6) Let us remark the appearance of an order h piece in the variation of the perturbation that will be relevant in our analysis.
In order for a flat spacetime theory to be scale invariant, the (Rosenfeld) energy-momentum tensor must be a total derivative (on-shell) where V μ is the virial current [10][11][12]. This is enough to guarantee the existence of a conserved scale current In the particular case when the virial current is itself a divergence, that is, when V μ = ∂ ν σ μν (9) then the theory is conformally invariant under the group O (2, n), and the conserved current reads (10) which also implies that the energy-momentum tensor can be improved.
In the present paper we want to clarify the precise relationship between Weyl invariant theories (WIFT) and conformal invariant theories (CFT) for systems where the gravitational field is dynamic, but still approachable as a fluctuation of flat spacetime. Our analysis then concerns mostly spin 2 theories in flat spacetime as described by a rank two symmetric field in Minkowski space. Our plan is to do it systematically, determining the conditions for scale invariance (which is still meaningful in flat space), conformal invariance and Weyl invariance.
We analyze first the most general Lagrangian containing dimension 4 operators, and then we do the same analysis for dimension 6 operators, containing two and four derivatives respectively (operators appearing in the weak field expansion of gravitational theories linear and quadratic in the curvature). We then study dimension 5 and dimension 6 operators with two derivatives. The analysis is, in some sense, the continuation of the one in [13,14] and also in [15]. Recent works regarding conformal invariance and Weyl invariance include [16][17][18]. We are always (with the only exception of our discussion of the improvement of the energy-momentum tensor) refering to actions defined as spacetime integrals, so that we allow for integration by parts, in other words, when we claim that an expression vanishes, we mean this only up to total derivatives.

Symmetries of the low energy spin action
Let us revisit the possible symmetries of low dimension kinetic operators in spin 2 theories [13] in flat spacetime where the graviton is represented by a symmetric tensor h μν . Lorentz invariance will be assumed throughout the paper. If we want the field equations to be given by differential operators of (at least) second order, then the Lagrangian has got to incorporate at least two derivatives. Let us begin our study with the operators of lowest possible dimension.

Dimension 4 operators
Our building blocks are the gravitational field, h αβ (assumed to be of mass dimension 1) and the spacetime derivatives, ∂ μ . There are four different dimension four operators with two derivatives There is another operator which is equivalent to D 2 modulo total derivatives: Then, the most general action principle involving dimension 4 derivative operators reads First, we are going to consider invariance under linearized diffeomorphisms, LDiff gauge symmetry. This is the one implemented in the pioneering work by Fierz-Pauli The variation of the fragment of the action containing the D terms only (that is, with m 2 i = λ α = 0) and upon integration by parts, yields so that linearized diffeomorphisms (LDiff henceforth) imposes some relations among the coupling constants this is, α i = 1 ∀i. The second case we are analyzing is the invariance under transverse linearized diffeomorphisms (LTDiff henceforth), that is, diffeomorphisms such that their generating vector fields obey These conditions impose but allow for arbitrary values of α 3 and α 4 . In third place, under linearized Weyl transformations, LWeyl, the variation of the metric reads (20) and after integration by parts, 2 where we have multiplied by κ 2 for simplicity. The invariance under LWeyl puts further constraints on the coupling constants, namely In [13] we have dubbed WTDiff to the theory with TDiff invariance enhanced with linearized Weyl symmetry, LWeyl. This is the particular case of the above, corresponding to α 1 = α 2 = 1 and A consistent non-linear completion of the actions which fullfill these requirements are the ones explained in the Appendix (110) namely actions proportional to Finally, we could consider only traceless graviton fields, h αβ such that h ≡ η αβ h αβ = 0. Obviously, in this case, D 3 = D 4 = 0, and for consistency, we can only implement TDiff with the coupling constants fixed to

Scale and conformal invariance
The most general Lagrangian we are considering (without the mass terms) is obviously scale invariant under with the assigned scaling dimensions. In order to make a full analysis of the scale and conformal invariance of the theory we have to compute the energy momentum tensor of these theories. In this case, and neglecting total derivatives, the metric (or Rosenfeld's) energy momentum tensor has the form whose trace reads The equations of motion (eom) read In fact we know that, as explained in the introduction, scale invariance implies that the trace of the energy momentum tensor can be written, on-shell, as a total derivative, that is, as the divergence of the virial (7). What happens is that whereas two Lagrangians that differ by a total derivative still generate the same eom, the specific form of the virial does depend on the particular form of the Lagrangian. This, in turn, also determines whether the virial itself can be written as a total derivative. Examining the contributions of the whole set of total derivative operators, leads to the convenient action The trace of the energy momentum tensor can then be rewritten as so that on-shell, we can express the trace as T = ∂ μ V μ , where It is also well-known that a conformal current can be formally constructed in case the virial can be so expressed as V μ = ∂ ν σ μν . In our case, and upon using the eom, In that way we get such a tensor (on-shell), for any value of the coupling constants. 3 Once a particular form of an σ μν is found, there is a systematic way of improving the energy-momentum tensor [10]. The improvement consists on adding another piece to the initial energy momentum tensor, so that the trace of the new energy-momentum tensor is precisely cero, that is, we avoid the total derivative terms. The piece in [10] has the form where X λρμν is symmetric (μ, ν), and divergenceless. The precise for of the improvement reads 3 We insist that this result (because this tensor is after all a total derivative) depends on the boundary terms that are neglected in order to write down the original Lagrangian. The monomials have to be written down as indicated; to be specific, in order to split D 2 in the two pieces, there is an integration by parts so we have two total derivatives appearing in order to interchange the two derivatives. Although this does not contribute to the equations of motion, it does contribute to the virial, If we take this contributions into account one of the total derivatives cancels one of the pieces of the virial proportional to α 2 and the other one is summed with the other piece. We end up with and we could not write this as the derivative of a two-index tensor.
where σ μν + stands for the symmetric part of σ μν . The original analysis [11] was specific for n = 4 dimensions but it can be generalized to arbitrary dimension [19].

Dimension 6 operators
In this section, we want to study the dimension 6 operators, and among these, there are various possibilities. First we take operators with four derivatives and two h αβ . After integration by parts, those are There is a small caveat here. There are also many operators with two derivatives and four h αβ , that will be analyzed in the next section. It is the case that these operators do not appear as a limit of quadratic Diff or TDiff invariant theories; they can only appear as higher order contributions to Lagrangians linear in the curvature. Let us consider the general theory involving dimension 6 operators which can come from quadratic theories of gravity, namely Like in the previous section, we first study the LDiff symmetry, upon which the O transform as after having integrated by parts. Then the symmetry under L Di f f imposes the following relations between the coupling constants These still allow for arbitrary values of g 1 and g 2 , and In the second place we consider invariance under LTDiff, which imposes Finally for LWeyl symmetry, the variations read (multiplying by κ/2 again for simplicity) so that the action is invariant under such tranformations whenever Now it is interesting to combine LDiff and LWeyl. In the case of dimension 4 operators, actions which are invariant under both symmetries do not exist. For dimension 6 operators, we can have LWDiff invariant theories as long as the coupling constants are constrained to have the following relations These actions with LWDiff invariance are obtained as the weak field limit of the following quadratic theories Note that the term √ g is immaterial at the order we are working. The weak field expansion of the quadratic invariants is worked out in the Appendix A. For n = 4 spacetime dimension these theories can be rewritten as where E 4 is the four-dimensional Euler density It can be easily checked that the weak field expansion of E 4 around Minkowski spacetime vanishes. This is the origin of the arbitrary coefficient α in the above expression. The quantity W 4 is the square of the 4-dimensional Weyl tensor The actions which are precisely proportional to the Weyl squared tensor are the ones with 3g 2 = α. Nevertheless, at the linear level, the Euler density does not contribute, so that all the solutions will effectively correspond to actions proportional to W 4 . Also quite interesting are those actions that are LWTDiff invariant; that is LWeyl invariant, but LDiff invariant under transverse diffeomorphisms only. They are characterized by The most general quadratic WTDiff invariant Lagrangian is the one obtained by Weyl transforming the metric in the quadratic action withg μν = g −1/4 g μν (this transformation ensures TDiff and automatically introduces a Weyl invariance). The expansion around flat spacetime reads The weak field limit of these theories automatically satisfies the constraints needed for LWTDiff (51). The precise form of these theories after the Weyl transformation (120) is shown in the Appendix.

Scale and conformal invariance
Scale invariance is now lost with the assingment given to h μν (conformal weight one). If we have a theory incorporating dimension 6 operators only, it is possible to recover scale invariance, just by making the graviton inert (conformal weight 0). It is plain that this does not hold when we have both dimension 4 and dimension 6 operators in the theory.
On the other hand, the conformal invariance demands as usual, tracelessness of the (metric, or Rosenfeld) energymomentum tensor. In this case, the energy-momentum tensor takes the form and the trace reads The eom read δS and they imply that, on-shell, the Lagrangian indeed vanishes, L = 0 (up to total derivatives).
In order to study the virial in detail, let us start from the specific form of the Lagrangian Now it is a simple matter to show that on-shell, L = ∂ μ ∂ ν σ μν with This result is somewhat puzzling, because we have already indicated that this theory is not even scale invariant with the standard assignment of conformal weight for the graviton field (namely 1). The result is however logical if we remember that the low energy of the Weyl squared Lagrangian is of this form. The theory containing both dimension 4 as well as dimension 6 operators, should not be conformal however. This fact can be easily understood from a simpler example, namely a scalar Lagrangian, where all complications of indices can be avoided. Consider then the Lagrangian which is equivalent to (up to total derivatives) This is our starting point. The eom read and the energy-momentum tensor The trace of the above reads Even if we are working in n = 4, we leave n arbitrary to maintain the second piece and illustrate the point we want to make. Note that this is not proportional to the total Lagrangian, because the trace counts the number of derivatives. We can rewrite the trace as which fails to be a total derivative when both α and β are nonvanishing, because Note that this is true even if there are WTDiff (that is TDiff and Weyl invariant) theories linear as well as quadratic in the Riemann tensor. The weak field limit of those Weyl invariant theories fails to be conformal invariant.

Dimension 5 and dimension 6 operators (with 2 derivatives)
Next, we study dimension 5 and dimension 6 operators containing just two derivatives, so that they come from the weak field limit of theories linear in the curvature, when expanded to higher orders in the perturbation. In the previous sections, operators coming from the lowest (non-trivial) order of gravitational actions were analyzed. In that cases, the lowest order of (T)Diff and Weyl variations was enough to obtain the conditions for those actions to be invariant under such symmetries. In this case, however, different orders of the expansion are needed because of the two orders involved in the field variations This translates into dimension 4 operators mixing with dimension 5 ones, and dimension 5 operators with dimension 6 ones. A full list of the independent dimension 5 and dimension 6 operators (containing 2 derivatives) can be found in Appendix 1. The most general dimension 5 Lagrangian with such operators reads Again, the only diffeomorphism invariant combination corresponds to √ g R, in this case, to the order O(κ 3 ) expansion of it This piece then combines with the previous order of the expansion to attain diffeomorphism invariance We can also look for the most general Lorentz and Weyl invariant Lagrangian built with this kind of operators. Again, we need the dimension 4 operator part that will contribute with the O(h) piece of the Weyl variation, which already has two arbitrary constants appearing in it (22). Taking that piece into account, the most general Weyl invariant Lagrangian up to dimension 5 operators reads Let us insist on the fact that both pieces are needed so that Weyl invariance is attained, that is, up to certain order in the expansion, the previous order is needed for the computation of the invariance conditions. Thanks to the mixing of the different orders, more freedom is avalaible to attain invariance under the studied symmetries. In particular, there are 8 independent coupling constants a 1 , a 2 , a 3 , a 5 , a 9 , a 11 , c 1 and c 2 .
It is straightforward to see that W T Di f f (24) is a particular case of this general Weyl invariant Lagrangian, with the constants fixed to Accordingly, the most general dimension 6 Lagrangian with two derivatives can be written as In this case, the expansion of Finally, let us analyze the most general Weyl invariant Lagrangian up to dimension 6 operators. We have different pieces apearing in it. First of all, it contains the pieces up to dimension 5 that were computed in this section (69), together with the dimension 6 piece of two derivative operators, that combine with specific coefficients so that Weyl invariance is attained. Moreover, we have another Weyl invariant combination coming from dimension 6 operators containing four derivatives (45). Taking everything into account, the most general Weyl invariant Lagrangian up to dimension six operators is shown in Appendix 1.

Interaction terms
It would appear quite intuitive to think that there are no potential terms invariant under either Diff or Weyl invariance. This is based in our GR intuition, but let us get rid of those prejudices and carry on with our perturbative analysis. It is easy to systematize the perturbative expansion. Up to quartic interactions we have the monomials so that the most general potential up to dimension four will read We want to analyze the invariance under diffeomorphisms as if this h μν corresponds to the perturbation of the metric around flat spacetime (65), but we take another energy scale M instead of κ. The crucial point is that owing to the fact that the Diff variations contain an order zero piece and an order one piece in the perturbation, each order in the perturbative expansion of the variation of the potential contributes to both the lower and upper orders. Up to total derivatives and dimension four operators, it can be seen that the following interaction Lagrangian is diffeomorphism invariant In fact this is an iterative process, each term in the expansion determining the following. The final potential contains infinite terms depending on just one arbitrary constant with dimensions of mass, m. In fact this is exactly the weak field expansion of m M 4 √ g − 1 .
At this point our GR intuition strikes back and asks whether this is not precisely the expansion of the cosmological constant term. (In fact they do not quite fit).
Concentrating in the quadratic terms Not knowing anything on GR we would say that there is spontaneous symmetry breaking in the system and the ground state has shifted from h μν = 0 to h μν = −Mη μν , leaving behind a vacuum energy Fluctuations around the new vacuum state are damped (provided m M < 0) with a quadratic term as is not positive semidefinite except for traceless H T αβ when m M < 0. In order to reach a definite conclusion on positivity, higher order terms should be considered. To the extent that this is related to the weak field expansion of m M 2 √ g, we expect it to have a definite sign however. Similar reasoning as in the previous paragraph leads to a Lorentz and Weyl invariant potential In this case we have more freedom as more arbitrary constants appear with each order of the perturbative expansion. The quadratic piece can be written as Fluctuations around the minimum of the potential behave as which again is positive semisefinite only for traceless H T αβ or else for pure trace when m M < 0 as

Global Weyl invariance
There is another symmetry that can be studied in this context, which is global Weyl invariance, that is, when the Weyl scaling factor is just a constant δg μν = 2 g μν , ∂ μ = 0 ( 8 4 ) When we expand the metric around flat spacetime, g μν = η μν + κh μν , the linearized variation of the quantum fluctuation reads In the case of global (rigid) Weyl invariance where ω is constant, the variations of the operators quadratic in the fields have to be computed taking into account both terms in the above (that is, the linear order in the quantum field). If we just took the first piece, proportional to the Minkowski metric, all the variations computed in (21) and (44) would just be total derivatives, which have been neglected in this work.
In order to illustrate this point, let us take two simple actions. We know that the Einstein Hilbert action is not globally Weyl invariant in four dimensions, On the other hand, we can take the simplest quadratic action which is invariant in four dimensions These equalities have to be true order by order in the perturbation of the metric. In this case, the quadratic order in the variation together with the linear order in the Weyl variation, combines with the third order of the perturbation in the action and the lowest order in the Weyl variation. These terms are going to be of order O(κ 2 ) and have to match exactly the O(κ 2 ) part of the rhs of the equation. Namely, The quadratice xpansion of the Einstein Hilbert action reads It is straightforward to see that taking the variation δh μν = 2ωh μν we get In order to compute the other piece contributing to O(κ 2 ) we need the third order of the expansion of the Einstein Hilbert action which contains terms with three quantum fields h μν and two derivatives, which are shown in the Appendix 1.
Once we have this expansion, we perform the Weyl variation δh μν = 2ω η μν κ yielding Adding the two contributions which precisely yields the right hand side of (88) for n = 4.
We can see that in n = 2 the Einstein Hilbert action is globally Weyl invariant (as well as locally). In fact this is basically the reason all two-dimensional metrics are conformally flat.
In the case of the quadratic action we have As before, taking the Weyl variation proportional to the quantum field is straightforward and it yields For the other piece, we need the third order variation of the quadratic action which can be easily computed. After performing the Weyl transformation on the quantum field, δh μν = 2ω η μν κ , we get Summing both contributions,

Non-local extensions
There is a permanent temptation to avoid the Källen-Lehman spectral theorem (which states that the price to pay for having propagators that fall off at euclidean infinity faster than k −2 is to have negative norm states) by considering non-local theories. For example in [20,21] a non-local generalization of the dimension 4 operators has been considered, namely so that the general Lagrangian of this type will be (where c i (z) are analytic functions with dimensionless argument). The five functions c i (z), i = 1 . . . 5 (which are assumed to include the corresponding coupling constants) characterize the theory. The constants put in front are such that the LDiff Fierz-Pauli theory corresponds to The correspondence with the dimension 6 operators in (38) is as follows i.e. g 1 = M 2 c 5 (z), g 2 = −2M 2 c 3 (z), g 3 = 2M 2 c 2 (z),g 4 = −4M 2 c 1 (z) and g 5 = 4M 2 c 4 (z), in such a way that the conditions for LDiff invariance now translate into It is claimed in [20] that the theory is ghost-free provided that and the function c 1 (z) is chosen as an entire function, such as Note that both constraints, (102) and (103) are different and incompatible.
It is well-known, however, that non-local theories suffer from unitarity and causality problems, some of those can be sometimes hidden uunder the rug of experimental precision of the measurements [22]. However, in order to do that, the theory needs to be quasi-local, which means that the corresponding function has got to have bounded support, which seems to contradict other conditions. It is not clear at all that a consistent solution exists.
Outstanding problems in this respect according to [23] are first and foremost, the fact that the presence of the exponential damping factor in the propagator prevents analytic continuation from the riemannian theory to the lorentzian one, owing to the essential singularities in the complex energy plane. It must be stressed, however that such an analytic continuation is problematic in any theory involving the gravitational field. Another argument is that none of the theories proposed so far complies with reflexion-positivity, which is believed to be an essential requirement in order to get a consistent quantum field theory.

Conclusions
In this paper we have presented a complete analysis of operators up to (mass) dimension 6 describing spin 2 theories (e.g. weak field limit of theories linear and quadratic in the curvature), analyzing with some care the conditions for the theory to be (transverse) diffeomorphism invariant, scale invariant, conformal invariant and Weyl invariant. We have also identified a possible non-linear completion of those Lagrangians.
Conformality on shell is attained for any combination of the constants appearing in the dimension 4 and dimension 6 cases. The trace of the energy-momentum tensor is a total derivative, and besides the virial current for specific Lagrangians is also the derivative of a two-index tensor, leading to improved forms of the corresponding energymomentum tensors.
On the other hand, Weyl invariance instead does impose constraints on the coupling constants. Our main conclusion is to confirm [16][17][18]24] that Weyl invariance and confor-mal invariance are independent symmetries: not every Weyl invariant theory is conformal invariant in the weak field limit and conversely, not every conformal invariant theory is Weyl invariant in spite of the fact that it is always invariant under global such Weyl transformations. To illustrate the first part of this statement, let us take for example the following WTDiff invariant theory where the precise form of these terms after permoning the transformation of the metric can be found in (109) and (120). The weak field expansion of this theory will contain, at quadratic order in the perturbation, dimension 4 operators and dimension 6 operators coming from the linear and quadratic (in curvature) pieces respectively. Theories combining operators of different dimension are not scale invariant, as pointed out in the example in (58). The analysis of dimension 5 and dimension 6 operators does not bring anything new with respect to diffeomorphism invariant theories, as expected. However, we have given expressions for the most general Lorentz and Weyl invariant Lagrangians up to dimension 5 and dimension 6 operators, and we can clearly see that those theories contain an increasing number of arbitrary constants. We have also discussed global Weyl invariance and it is clear that this symmetry is less restrictive than the local one. An analysis of the interaction terms has been done. It can be seen that potentials with diffeomorphism and weyl invariance can be constructive iteratively, for every orther of the perturbative expansion.
To end up, let us stress that the conditions that are argued to be neccessary for a ghost free non-local theory [20] are not compatible with the ones stemming from diffeomorphism invariance.
We finally point out that our results prove that any Lorentz invariant Lagrangian for spin 2 particle up to quadratic order in the field is conformal invariant. under the Marie Sklodowska-Curie grants agreement No 674896 and No 690575. We also have been partially supported by FPA2016-78645-P(Spain). RSG is supported by the Spanish FPU Grant No FPU16/01595.

Data Availability Statement
This manuscript has associated data in a data repository. [Authors' comment: All data generated or analysed during this study are included in this published article.] Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, pro-vide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecomm ons.org/licenses/by/4.0/. Funded by SCOAP 3 .

A Weak-field limit of geometric scalars
We are interested in the expansion of the geometric invariants when we expand the metric around Minkowski spacetime g μν = η μν + κh μν (106) If we take the limit of linear theories of gravity up to quadratic order in the fluctuations we have When considering TDiff scalars one can also have terms of the type The existence of this operator gives one extra freedom. To build the action which is WTDiff invariant we perform a Weyl transformation in the usual Einstein Hilbert action takingg μν = g −1/4 g μν so that Expanding it up to quadratic order in the fluctuations and writting it in terms of the four dimensional operators (11) we get On the other hand, taking into account that we learn that For the Ricci tensor we have so that Finally the expansion of the Ricci scalar reads and it follows that A useful relationship is given by Using this it can easily be seen that the Euler density vanishes at this level of the expansion, whereas the Weyl squared tensor decomposes into If we again consider quadratic theories which are TDiff invariant, one would have terms of the type must add again, this yields in this case one extra freedom. We can make the same analysis for the quadratic invariants but when considering actions that are WTDiff invariant. This can be achieved by making a Weyl transformationg μν = 2 g μν on the usual quadratic action (117) and then taking 2 = g −1/n . For a general we have αR μνρσR μνρσ + βR μνR μν + γR 2 = −4 α R μνρσ R μνρσ + β R μν R μν + γ R 2 + −5 (−8α − 2(n − 2)β) R μν ∇ μ ∇ ν + −6 4α + (3n − 4)β + 4(n − 1) 2 γ ( ) 2 + −6 4(n − 2)α + (n − 2) 2 β ∇ μ ∇ ν ∇ μ ∇ ν + −6 (−4α − 2(n − 3)β −2(n − 1)(n − 4)γ ) R∇ μ ∇ μ + −6 (16α + 4(n − 2)β) R μν ∇ μ ∇ ν × −7 8(n − 3)α + 4(n 2 − 5n + 5)β Using that = g −1/2n (in order to have WTDiff) and keeping dimension six operators with four derivatives and two metric fluctuations, we get for n = 4 αR μνρσR μνρσ + βR μνR μν + γR 2 These are the most general theories that possess LWTDiff.

B Dimension 5 and dimension 6 operators with twoderivatives
This set of operators does not appear in the expansion of terms quadratic in the Riemann tensor, although they appear in the expansion of the Einstein-Hilbert Lagrangian. For dimension 5, there are 14 independendent operators (up to total derivatives), that form a basis to expand the most general Lorentz invariant Lagrangian containing such opeators