Embeddings of the “New Massive Gravity”

Here we apply different types of embeddings of the equations of motion of the linearized “New Massive Gravity” in order to generate alternative and even higher-order (in derivatives) massive gravity theories in D=2+1\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$D=2+1$$\end{document}. In the first part of the work we use the Weyl symmetry as a guiding principle for the embeddings. First we show that a Noether gauge embedding of the Weyl symmetry leads to a sixth-order model in derivatives with either a massive or a massless ghost, according to the chosen overall sign of the theory. On the other hand, if the Weyl symmetry is implemented by means of a Stueckelberg field we obtain a new scalar–tensor model for massive gravitons. It is ghost-free and Weyl invariant at the linearized level around Minkowski space. The model can be nonlinearly completed into a scalar field coupled to the NMG theory. The elimination of the scalar field leads to a nonlocal modification of the NMG. In the second part of the work we prove to all orders in derivatives that there is no local, ghost-free embedding of the linearized NMG equations of motion around Minkowski space when written in terms of one symmetric tensor. Regarding that point, NMG differs from the Fierz–Pauli theory, since in the latter case we can replace the Einstein–Hilbert action by specific f(R,□R)\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$f(R,\Box \, R)$$\end{document} generalizations and still keep the theory ghost-free at the linearized level.


Introduction
Massive spin-2 particles can be covariantly described by means of a symmetric rank-2 tensor. Although this is not the only possible tensor structure, it is very convenient. It is closely connected with a geometrical point of view (fluctuation about some metric) and it is a minimal description in the sense that we need just one auxiliary field, i.e., the trace of the tensor which vanishes on shell. If we further a e-mail: dalmazi@feg.unesp.br b e-mail: elias.fis@gmail.com require a second-order theory, in derivatives, we end up with a unique answer: the Fierz-Pauli (FP) theory [1]. Almost all developments in massive gravity, from earlier works [2][3][4][5] until recent developments, see [6,7] for review articles, are built up on top of the free FP theory. It is remarkable that the absence of ghosts [4] and of mass discontinuity [2,3] have been both achieved in recent theories with one [8,9] and two [10] dynamic metrics. A good question concerns the uniqueness of those massive gravities; see for instance [11][12][13][14]. Moreover, which features of those theories are model independent? Can we still dream of a renormalizable massive gravity? One way of addressing those questions is to abandon the FP paradigm as a starting point and allow for higher-derivative kinetic terms. In fact, the reader can find higher-derivative massive gravities and discussions regarding their physical consequences in [15][16][17][18][19].
At the linearized level, those higher-derivative models are specific generalizations of the FP theory which, however, have the same spectrum: massive spin-2 particles without ghosts. This is remarkable, since usually higher-derivatives introduce ghosts. The above result requires of course a systematic investigation of all possible higher-derivative generalizations of the FP model. In [20] we have addressed this problem.
One way of producing higher-derivative models dual to some "lower-derivative" theory is by means of the embedding of its Euler tensor (equations of motion). The method consists of adding to the action of some lower-derivative starting theory, quadratic terms in its equations of motion with coefficients which are functions of = ∂ μ ∂ μ . This guarantees that the equations of motion of the original theory also minimize the new action. We say that the equations of motion of the starting theory are embedded into the new theory. The particle content of the lower-order theory is inside the new theory. Then the coefficients are fixed such that the new theory contains exactly the same particle content of the lowerorder one without extra propagating modes. Such method produces alternative dual models of higher order which by themselves are not very useful as free models; however, in some cases those free models correspond to quadratic truncations of nonlinear models which are not equivalent in general to nonlinear versions of the starting lower-order theory. In other words, the embedding procedure and the addition of nonlinear terms do not commute. As a successful application of this method we can mention the "New Massive Gravity" (NMG) theory of [21] which is of fourth-order in derivatives and can be obtained from the usual FP theory (of second order) by the addition of quadratic terms in the FP equations of motion with specific constant coefficients as explained in [22]. Another example is the higher-derivative topologically massive gravity of [22,23] which is obtained from the topologically massive gravity (TMG) of [24]. The TMG itself, of third order, can be obtained at the linearized level from consecutive embeddings of the first-order spin-2 self-dual model of [25].
In [20] we have investigated, to all orders in derivatives, all possible embeddings of the equations of motion of the FP theory which are ghost-free at the linearized level. We have found a system of equations allowing for several solutions for the coefficients. Those theories are in general of higher order in derivatives but still ghost-free at the linearized approximation. They can all be nonlinearly completed with the help of a fiducial metric. Most of them are f (R, R) modifications of the Einstein-Hilbert kinetic term plus the FP mass term. Although a complete analysis has not yet been carried out, at least a subset of such theories can be modified with an appropriate non-derivative potential of the type suggested in [8,9] and become, apparently, ghost-free beyond the linearized approximation, as shown in [17].
It turns out that all higher-order models obtained in [20] in arbitrary D-dimensions correspond to modifications of the FP theory in the spin-0 sector of the propagator. Any modification in the spin-2 sector leads to ghosts. There is however, one exception in D = 2 + 1. Namely, the "New Massive Gravity" (NMG) theory 1 of [21]. Although of fourth order in derivatives, the NMG is ghost-free, though it is apparently still not renormalizable [27]. Here we take the NMG model as our "lower-order" starting point and try different types of embeddings in order to produce alternative massive gravity models in D = 2 + 1. In this sense the present work complements the work [20] by taking care of the special case of D = 2 + 1.
Sometimes the embedding of equations of motion leads to new theories with gauge symmetries not present in the lowerorder starting theory. In those cases, as in [22], one may call it a Noether gauge embedding (NGE). We have found useful to split the embeddings here into two categories. Namely, the ones based on requiring Weyl symmetry of the new model (Sect. 2) and the other ones (Sect. 3). In the first part of Sect. 2 we require Weyl symmetry without introducing any extra field. In the second part we introduce a scalar Stueckelberg field and obtain a new scalar-tensor theory for massive gravitons in D = 2 + 1. In Sect. 3 we look at rather general local embeddings where quadratic terms in the NMG equations of motion are added to the NMG theory with coefficients which are arbitrary functions of . We examine the propagator of the final higher-order model and require equivalence of the particle content, i.e., the new dual theory must describe massive spin-2 particles and nothing else. Differently from the FP case we show here (Sect. 3) to all orders in derivatives that there is no local ghost-free embedding of the NMG equations of motion. In particular, the NGE embedding of the Weyl symmetry of the first part of Sect. 2 also leads to a ghost. Only the Stueckelberg approach leads to a ghost-free model which, after elimination of the scalar field, becomes nonlocal. The work in Sect. 3 is based on the analytic structure of the propagator. In Sect. 4 we present some final comments and our conclusions explaining that the negative results of Sect. 3 are related with the linearized reparametrization invariance of the NMG theory which makes it hard to be embedded.

Weyl embedments of the "New Massive Gravity"
By a systematic Lagrangian procedure, called Noether gauge embedment (NGE), one can deduce a gauge invariant massive theory out of a non gauge invariant one [28]. The gauge symmetry of part of the initial Lagrangian is extended to the whole final theory. However, there is no guarantee that the particle spectrum is preserved. In [22] we have shown that the linearized NMG theory can be obtained via NGE from the usual Fierz-Pauli theory via embedding of linearized reparametrizations. We have also shown that a linearized higher-derivative topologically massive gravity is obtained from the usual linearized topologically massive gravity of [24] via NGE of the Weyl symmetry. In all those cases the particle content is preserved. One could wonder what would be the gauge invariant action obtained from the linearized New Massive Gravity theory. Since part of the action of NMG, the curvature square term, is invariant under Weyl transformation δ W h μν = φη μν , one might try to embed this symmetry into a new theory. This is what we next do. The linearized NMG theory can be written, up to an overall constant, as where G μν (h) is the usual linearized self-adjoint Einstein tensor and S μν (h) is the linearized Schouten tensor in D = 3 defined as R μν (h) − η μν R(h)/4 hh which has the useful property of "commutativity" with the Einstein tensor in the sense that inside integrals In the NGE procedure an important ingredient is the Euler tensor, With the help of an auxiliary field a μν such that δ W a μν = −δ W h μν we implement a first iteration of the form The Weyl variation of (3) can be written as Therefore we end up with the Weyl invariant theory, Noticing that the Euler tensor (2) can be written in terms of the Einstein tensor, i.e., After the shift a μν →ã μν − H μν in (6), theã μν auxiliary field decouples. We can safely discard the last term,ã μν G μν (ã), which is a linearized Einstein-Hilbert term without particle content. Thus we have a sixth-order Weyl invariant action which turns out to have a nonlinear completion, Since the tensor structure of (7) is the same as the curvature square term of the NMG theory, it is clear that S W is invariant under Weyl transformations. The particle content of S W will be examined in the next section. The theory S W contains a ghost.
Another way to embed the Weyl symmetry in the "New Massive Gravity" is to introduce a scalar Stueckelberg field in (1) by substituting h μν → h μν + η μν φ. Since the fourthorder term of S NMG is Weyl invariant we end up with By construction, the linear theory It is clear that the equations of motion δ S N L φ = 0 contain the trivial solution g μν = η μν , φ = 0. Expanding about such vacuum up until quadratic terms in the fluctuations we recover S L φ . Therefore, the particle content of S N L φ consists, at tree level, of one massive spin-2 particle just like the NMG of [21]. However, as in the K-model (massless limit of NMG) studied in detail in [29,30], we might have problems at nonlinear level since the linearized Weyl symmetry is probably broken at nonlinear level and consequently the scalar field stops being pure gauge in the full model (51). In particular, the phenomenon of bifurcation of constraints found in [30] might also be present here. A detailed study of the constraint structure should be carried out.

Generalized Euler tensor embedment of "New Massive Gravity"
Our starting point is the linearized NMG theory with the addition of quadratic terms in its equations of motion: The first term in (10) is the linearized NMG theory. The NMG Euler tensor K μν is given in (2). The coefficients d( ) and f ( ) are so far arbitrary functions of = ∂ μ ∂ μ such that the Lagrangian L G remains local. Due to the conservation law ∂ μ K μν = 0 which holds identically due to the linearized reparametrization invariance of the NMG (δh μν = ∂ μ ξ ν + ∂ ν ξ μ ), other terms which might show up in (10) like (∂ μ K μν ) 2 and ∂ μ K μν ∂ ν K do not contribute. 2 In terms of the original field h μν we can rewrite L G as The coefficients c i ( ) are given by 3 The Lagrangian L G can be further rewritten in terms of a four index differential operator L G = h μν G μναβ h αβ whose inverse G −1 does not exist due to the linearized reparametrization symmetry. We choose the de Donder gauge fixing term: which amounts to the shift c 1 → c 1 + λ , c 2 → c 2 − λ in Eqs. (12) and (13) but not in (14) and (15). In (14) we make c 3 → c 3 − λ /4. Before we display G −1 we take a closer look at the coefficients c i ( ) in (12)- (16). A local Lagrangian density in terms of h μν requires whered( ) andf ( ) are analytic functions of while a is an arbitrary real constant. In terms of the spin-s projection operators P (s) I J given in the appendix A and suppressing the four indices we have 2 We have restricted ourselves here to parity invariant theories, thus avoiding terms like μνα K μγ ∂ ν K γ α . 3 Henceforth, unless otherwise stated, we replace d( ), f ( )) by d, f respectively, though they are still arbitrary functions of .
Regarding the local symmetries of L G there is one special case, since we have a zero in the denominator of the spin-0 sector which indicates a spin-0 symmetry. Indeed, under a Weyl transformation δ W h μν = η μν we have from (10) after integrations by parts (22) In the special case (21) we need to add another (Weyl) gauge fixing term, we may choose L W G F = ξ h 2 . This implies the shift c 3 → c 3 +ξ . Consequently, in the Weyl symmetric case we have Next we analyze the particle content of L G from the analytic structure of G −1 and G −1 W . In momentum space we can calculate the gauge invariant two-point amplitude A(k) by saturating G −1 or G −1 W with external sources. For instance, Here G −1 (k) = G −1 (∂ μ → i k μ ). Due to the linearized reparametrization symmetry, the source must be transverse k μ T μν = 0, consequently, where d =d(−k 2 ) and f =f (−k 2 ) are analytic functions of k 2 = k μ k μ and The quantity T = η μν T μν = −T 00 + T ii is the trace of the external source in momentum space. If the two conditions for Weyl invariance (21) hold we must have T = 0.
A key role is played by the imaginary part of the residue of A(k) at each pole. For instance, at k 2 = −m 2 we have If and only if I m > 0 we have a physical particle. If I m = 0 we have a non-propagating mode while I m < 0 or no definite sign for I m signalizes the presence of ghost. In order to verify the sign of I m we fix a convenient coordinate frame splitting the cases of massless and massive poles. In the massless case we fix a frame such that k μ = (k 0 , , k 0 ), thus k μ k μ = 2 . We will take → 0 at the end. This caution is necessary for the analysis of double poles. From the three conditions k μ T μν = 0 we have in this frame Consequently, |T | 2 = |T 11 | 2 − 2 k 0 (T 12 T * 11 + T * 12 T 11 ) In the case of massive poles we choose the frame k μ = (m, , 0) such that k 2 + m 2 = 2 . From k μ T μν = 0 we have Thus, Since d(k 2 ) and f (k 2 ) are arbitrary analytic functions, there might be double poles in the denominator of A(k). We first examine those poles. From (26), (27), (35), and (36) we see that is impossible to take linear combinations of S (0) and S (2) in order to end up only with terms of order 2 . Therefore a massive double pole 1/(k 2 + m 2 ) 2 = 1/ 4 cannot be reduced to a simple pole by any fine tuning of the functions d and f . So henceforth we assume that all massive poles must be simple poles. The conclusion remains the same for the Weyl symmetric case. In the latter case T 00 = T 11 + T 22 and (34) imply T 22 = −T 11 + ( 2 /m 2 )T 11 , which does not help canceling the term |T 12 | 2 in (35). The massless case is a bit different. From (26), (27), (32), and (33) we see that we do have one special combination of order 2 which may turn double poles into simple ones, namely, However, the term (T 11 T * 22 + T * 11 T 22 ) has no definite sign. So we end up with a ghost. In the special case of Weyl symmetry, using (31) in T 00 = T 11 + T 22 we have T 11 = 2( /k 0 )T 12 + O( 2 ). Thus, the dangerous term of (37) becomes of order 3 and will not contribute to the residue. So we may hope to turn a double massless pole 1/k 4 into a physical pole only in the Weyl invariant case after a specific fine tuning of d and f . In particular, this is the mechanism behind the fourthorder K-term which describes a physical massless particle as explained in [29] via decomposition of h μν in orthogonal modes and in [31] via the analytic structure of the propagator.
In summary, multiple poles lead us to ghosts in general except for the double massless pole in the Weyl symmetric case which will be examined later on.
Henceforth we split our analysis in four cases: In the cases III and IV we have Weyl symmetry. The case IV corresponds to the S φ model of (8) after elimination of φ.
3.1 Case I: a = 1 As a warm up we start reproducing the results of [26] for the NMG theory. We take a = 0 = d = f . We have one massless and one massive pole, Taking → 0 in (32) and (33) we have a vanishing residue at the massless pole and consequently a non-propagating mode: The residue at the massless pole vanishes for the very same reason as it does in the Maxwell-Chern-Simons theory of [24], namely, the lowest-order term (in derivatives) of the theory (linearized Einstein-Hilbert) has no particle content. Taking → 0 in (35) and (36) we have a positive residue at the massive pole, a physical massive spin-2 particle, Now we go back to the general case a = 1. Except for the NMG case a = 0 = d = f , which will not be treated here anymore, we have in general extra massive poles stemming from the polynomial Q(k 2 ). Requiring that no tachyons show up we can write where N Q ≥ 1 is the number of extra massive poles coming from Q(k 2 ). Since we are specially interested in the massive poles, we can write from (26), (27) and (35), (36) at → 0, The quantity A(k 2 ) is an analytic real function of k 2 whose specific form is not important, it is defined by comparing (41) with (25). Defining the polynomial of degree N Q + 1: where m 2 0 ≡ m 2 is the mass squared already present in the NMG theory, it is clear that the sign of the residue I m j at some pole k 2 = −m 2 j depends essentially upon the sign of the ratio S A /P calculated at k 2 = −m 2 j , where P = dP/dk 2 . Since the derivative of a polynomial has alternating signs at its consecutive simple zeros, the only hope of having positive residues at the different massive poles is to require that S A also has alternating signs at such points. However, it is easy to prove that S A (−m 2 j ) either has no definite sign or is definite positive. The point is that if A(−m 2 j ) ≥ 0 we can guarantee that the last three terms of (42) add up to a nonnegative number, so in those cases S A (−m 2 j ) ≥ 0. On the other hand, since the three complex numbers T 12 , T 11 , T 22 are totally unconstrained, even if we take A(−m 2 j ) < 0, depending on the relative strength of those three complex numbers, the sign of (42) may change. Thus, we cannot guarantee that S A (−m 2 j ) < 0. In conclusion, whenever we have more than one massive pole we have ghosts and only the NMG case is safe at a = 1.

Case II: a = 1 and d + 2 f = 0
In this case we have in principle a double massless pole and massive poles: We can choose d(0) + f (0) = 1 2 m 2 and turn the double massless pole into a simple one. However, since we have no Weyl symmetry, as explained in the paragraph of Eq. (37), we are doomed to have a massless ghost. Therefore we go to the next case. where The above Lagrangian can be nonlinearly completed in terms of the square of curvatures in the form of a K-term of the NMG theory, i.e., In the casef = 0 we recover the Weyl embedding of the previous section; see (7). Back to the two-point amplitude A(k): we have again a double massless pole and massive poles in general. Assuming that the analytic function f (k 2 ) =f ( → −k 2 ) is such that we have no tachyons, we can write where P(k 2 ) is defined in (43). From (27), (35), and the fact that T = 0, due to the Weyl symmetry, we have at each massive pole .
Due to the alternating signs of P at its consecutive single zeros, it is impossible to have I m j > 0 for all j = 0, . . . , N Q . We are forced to assume N Q = 0, i.e., f = 0. In this subcase P(k 2 ) = k 2 + m 2 , so P = 1 and the massive pole is a physical one I m = I m 0 = 2(|T 12 | 2 + |T 11 | 2 ) > 0. Regarding the massless double pole, we have already seen that due to the Weyl symmetry we have T 11 = 2( /k 0 )T 12 + O( 2 ), substituting back in (32) we obtain [31] T μν T * μν = −2 Consequently, although the apparent double pole has become a simple pole, we still have a ghost due to the negative sign of the residue: Therefore, the Weyl invariant theory (7) of the previous section will unavoidably contain a ghost. In this last case we have a nonlocal theory corresponding to the Weyl invariant action S φ given in (8) after the elimination of φ: After adding a gauge fixing term like (17) plus another one for the Weyl symmetry L W G F = ζ h 2 the action acquires the form (11) with the coefficients The propagator, suppressing indices, is given by After saturating the propagator with transverse and traceless sources as in (24) we are left with two simple poles, one massive and one massless which come both from the pure spin-2 sector: As expected, the dependence on the gauge parameters λ and ζ disappear which guarantees gauge invariance of the twopoint amplitude. When we look closer at the massless pole using k μ = (k 0 , , k 0 ), it is easy to see that its residue vanishes with power 2 . From (32), (33), and T = η μν T μν = 0 we have The massive pole is a physical one (positive residue). From (35) we have Therefore, the linearized version of the model (51) is unitary and contains only massive gravitons in the spectrum.

Conclusion
Recent work on massive gravity [8][9][10] has shown how to overcome longstanding problems like the appearance of ghosts at the nonlinear level [4] and the vDVZ [2,3] mass discontinuity. Those models can be described in terms of one [8,9] or two [10] dynamic metrics. In the linearized limit those theories reduce, respectively, to the old paradigmatic Fierz-Pauli (FP) theory [1] and to the addition of the linearized Einstein-Hilbert plus the FP theory. In fact, those new massive gravity models are built up on top of the second-order FP theory which describes free massive spin-2 particles. It is expected that if we change the starting point (underlying free theory) we might end up with alternative massive gravities. Those alternative massive gravities might help us to understand which physical features of the new theories are really model independent and eventually we might still hope of finding a renormalizable massive gravity model. In [20] we have shown that there are several higherderivative modifications of the Fierz-Pauli theory which describe free massive spin-2 particles and are still ghostfree at linear level. So they could be used as new starting points for alternative massive gravities. Indeed, they can be identified with quadratic truncations of f (R, R) modifications of the Einstein-Hilbert theory plus the Fierz-Pauli (FP) mass term; see [20]. Some of those models have been further changed, see e.g. [19], by the addition of a convenient non-derivative nonlinear potential of the type found in [8,9] in order to account for the absence of ghosts at nonlinear level. They define massive gravity theories with interesting cosmological properties.
All modifications of the FP model found in [20] occur in the spin-0 sector of the theory. Modifications in the spin-2 sector lead in general to a ghost, with the exception of the so-called "New Massive Gravity" of [21] which only exists in D = 2 + 1. In section III we have used the embedding procedure used in [20] in order to search for arbitrary higherorder modifications of the NMG model similarly to what has been done in [20] for the case of the FP theory. Although the embedding procedure does not guarantee full equivalence of the initial and final models, it ensures that the particle content of the initial theory is contained in the new theory. At the end one always has to check both particle contents and fix the embedding coefficients, see (10), in order to achieve full equivalence.
Contrary to the FP case, we conclude here that there is no local embedding of the NMG theory which remains ghost-free. We have carried out a thorough calculation of the residues at all possible massive and massless poles in the propagator. There is always a pole with negative residue (ghost). A key technical point is the fact that the Euler tensor K μν , which defines the NMG equations of motion K μν = 0, satisfies ∂ μ K μν = 0 identically due to the linearized reparametrization invariance δh μν = ∂ μ ξ ν + ∂ ν ξ μ . Consequently, the quadratic terms (∂ μ K μν ) 2 , (∂ μ ∂ ν K μν ) 2 , and ∂ μ K μν ∂ ν K identically vanish too and strongly restrict the embedding Ansatz (10), contrary to the embedding of the FP theory where we have five independent quadratic terms in the equations of motion which can be added to the FP theory. The more local symmetries we have, the more difficult is to embed it into a higher-order model.
In Sect. 2 we have used the Weyl symmetry as a guiding principle for the embedding of the NMG. First we have looked at the Noether gauge embedment of the Weyl symmetry, which leads to the sixth-order model (7). In this case we have one massive and one massless pole. Unfortunately, their residues have opposite signs. If we reverse the overall sign of (7), we have a physical massless graviton and a massive ghost. In this case we might use the Weyl invariant model (7) with reversed sign as a phenomenological toy model along the lines of [32]. Namely, we can have a consistent unitary theory if the ghost mass stays above the energy cut-off of the theory. Although the model (7) is of sixth order in derivatives, the would-be double massless pole is reduced to a simple pole and the analytic structure of the propagator is similar to some curvature square modifications of general relativity in D = 4.
At the end of Sect. 2 we have obtained a promising candidate for a consistent massive gravity different from the NMG model. From the introduction of a scalar Stueckelberg field in the linearized version of the NMG theory we have derived the linearized Weyl invariant model given in (8). As in the case of (7) we have one massless and one massive pole but differently from (7) both poles are simple poles. The Weyl symmetry now kills the residue at the massless pole such that the "would-be" massless ghost does not propagate at all as shown at the end of Sect. 3 (Case IV). The residue at the massive pole is positive and we are left with physical massive gravitons, at least in the linearized approximation. The model can be nonlinearly completed leading to the scalartensor theory (51) which might be an alternative to the usual NMG model. However, as in the case of the pure K-term (massless limit of NMG) analyzed in [30], the Weyl symmetry is probably broken beyond the linearized level which might lead to a ghost in the full theory. This is probably true also for the higher-derivative topologically massive gravity of [22,23]. Those examples require detailed investigations of the constraints structure which are beyond the scope of the present work.
Finally, we mention Ref. [33] where stability and unitarity of fourth-order (quadratic in curvatures) three dimensional gravities, including the parity breaking gravitational Chern-Simons term of [24], have been investigated around maximally symmetric spaces. Whenever our results overlap (fourth-order theories around Minkowski space without parity breaking terms) we have agreement.
Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecomm ons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. Funded by SCOAP 3 .

Appendix
Here we display the operators P (s) I J , the coefficients A i j ( ). They make use, as building blocks, of the spin-0 and spin-1 projection operators acting on vector fields, respectively, we define the spin-s operators P They satisfy the symmetric closure relation