Fierz–Pauli theory reloaded: from a theory of a symmetric tensor field to linearized massive gravity

Modifying gravity at large distances by means of a massive graviton may explain the observed acceleration of the Universe without Dark Energy. The standard paradigm for Massive Gravity is the Fierz–Pauli theory, which, nonetheless, displays well known flaws in its massless limit. The most serious one is represented by the vDVZ discontinuity, which consists in a disagreement between the massless limit of the Fierz–Pauli theory and General Relativity. Our approach is based on a field-theoretical treatment of Massive Gravity: General Relativity, in the weak field approximation, is treated as a gauge theory of a symmetric rank-2 tensor field. This leads us to propose an alternative theory of linearized Massive Gravity, describing five degrees of freedom of the graviton, with a good massless limit, without vDVZ discontinuity, and depending on one mass parameter only, in agreement with the Fierz–Pauli theory.


Introduction Motivations
It is an observational fact that the Universe is expanding at an accelerated rate [1,2].Within the theory of General Relativity (GR), the cosmological constant Λ explains this phenomenon.In a picture of the Universe seen as a perfect fluid, the cosmological constant gives rise to an accelerated expansion by acting as a constant energy density (called Dark Energy) ρ ∼ ΛM 2 P providing a negative pressure.A fine tuning of the constant Λ enables us to match the expansion predicted by GR to the one which is observed.The estimate from observations is Λ M 2 P ∼ 10 −65 [3], where M P is the Planck mass, while the quantum field theoretical prediction on the cosmological constant seen as the vacuum energy density gives Λ M 2 P ∼ 6×10 54 .Unfortunately, these two estimates disagree by about 120 orders of magnitude.This huge tension is known as the cosmological constant problem [4], which motivates alternative descriptions for Dark Energy, with the request that any cosmological model should reproduce an Universe which, at our epoch, is almost perfectly flat and filled by matter and DE in the ratio of about 3/7, where the DE is effectively approximated by a constant [5,6,7,8,9,10].Modifying gravity at large distances by means of a massive graviton is a possible solution to the cosmological constant problem [11].In fact, the Yukawa potential for a massive field at large distances goes like ∝1 r e −αmr , where m is the mass of the field and α a dimensional constant.At scales comparable to 1  αm , the exponential factor suppresses the potential and the strength of interactions as well.This is the reason why long-range forces are associated with massless bosons and short-range forces with massive bosons.By means of a Yukawa-like potential, the gravitational effect of the vacuum energy density is exponentially suppressed at large scales, thus explaining the disagreement between the quantum field theory calculation and the observed cosmological value.Clearly, this exponential suppression of long-range gravitational interactions is constrained by experimental evidence and, consequently, the mass of the graviton is subject to restrictive upper limits.Another way of understanding this point is to remember that a mass acts as a momentum space cutoff of interactions, effectively damping the effect of low frequency sources.The vacuum energy density can be considered constant, therefore its contribution to the expansion of the Universe is greatly reduced by the introduction of a mass for the graviton, again confirming that in a picture where the graviton is massive the cosmological constant problem might be fixed by a suitable graviton mass.Theories in which the graviton is massive are referred to as Massive Gravity (MG) theories.We nonetheless remark that the introduction of a mass for the graviton as a potential solution of the cosmological constant problem was a hope more than a decade ago, but early experience with nonlinear MG (see, for example [12]) showed that this expectation is not realized in the nonlinear theory where a cosmological constant bends spacetime in the same way as in GR.While a large vacuum energy can be canceled by a MG contribution, this still requires fine tuning and does not amount to a degravitation mechanism. 1In some sense MG is a substantial modification of GR, even if the mass of the graviton is extremely close to zero.In fact, when adding a mass, the two degrees of freedom (DOF) carried by the massless boson become 2S + 1 massive DOF.For the graviton, which is the spin S = 2 representation of the Poincaré group, the two massless DOF become five massive ones.

Phenomenological Limits
As previously said, long range forces are usually carried by massless bosons because the Yukawa exponential factor limits the range of interactions mediated by massive particles.Nevertheless, choosing a sufficiently small value of the graviton mass, the interactions at scales of the observable Universe remain essentially identical to those of ordinary gravity, while only interactions at greater scales are affected.Indeed, the observed accelerated rate of the Universe gives an upper limit to the mass of the graviton, which in [13] is shown to be of order 10 −34 eV c 2 , roughly an order of magnitude smaller than the Hubble parameter at our epoch.With such a limit, for interactions at distances much smaller than the radius of the observable Universe (which is estimated to be about 46.6×10 9 ly [14]), the Yukawa exponential drops to one, and the classical Newton potential ∝ 1 r is recovered, as desired.Still, this cosmological upper limit becomes irrelevant when the mass is acquired through the condensation of some additional scalar field (see [13]).Another constraint on the mass of the graviton was recently obtained by the LIGO-Virgo collaboration through the analysis of binary black holes merger signals, measuring the phase shift between components of different frequencies of the gravitational waves detected by the interferometers [15].In this way, the upper limit for the mass of the graviton is estimated to be 4.7 × 10 −23 eV c 2 , which is several orders of magnitude higher than the limit derived in [13] by cosmological considerations.

State of the Art
In 1939 Fierz and Pauli (FP) proposed a relativistic theory for a massive particle with arbitrary spin f , described by a symmetric rank-f tensor field [16].They showed that "in the particular case of spin two, rest-mass zero, the equations agree in the force-free case with Einstein's equations for gravitational waves in GR in first approximation; the corresponding group of transformations arises from the infinitesimal coordinate transformations".This is what is known as the FP theory, which is considered as the standard paradigm for Linearized MG (LMG).Much later, it has been realized by van Dam, Veltman and Zakharov [17,18] that the massless limit of the FP theory exhibits a discontinuity with GR, known as vDVZ discontinuity.In particular, in [19] it is shown that the interaction of light with a massive body, like a star for instance, is 25% smaller in the massless limit of the FP theory than in GR.Consequently, the angle of deviation of light is also different by 25% in the two theories, a difference which can be observed by means of the time delay in the gravitational lensing phenomenon.The problems of the FP theory at small scales can be solved by the Vainshtein mechanism [20], which nonetheless requires a nonlinear modification of the FP theory.In principle, this measurable disagreement between GR and the FP theory of LMG would allow us to distinguish between a strictly vanishing mass and a very small one.Such a tension is crucial, since GR has been extensively tested and as a consequence we might rule out the FP theory as a theory of LMG.It was also pointed out in [21] that in order for Linearized Gravity (LG) to describe a pure spin two system, it is necessary to use a particular FP mass term, referred to as the FP tuning.Otherwise, there will be admixtures of lower spin, in general with negative-energy, called Boulware-Deser ghost.Moreover, Boulware and Deser argued that the ghost will reappear if nonlinear extension of LG are considered.Some possibilities to circumvent these problems have been proposed in [22,23,24].A detailed review on MG can be found in [25].

Strategy
The approach adopted in this paper is based on a field-theoretical treatment of LG.We consider LG as a gauge theory of a symmetric rank-2 tensor field, the gauge symmetry being the infinitesimal diffeomorphism invariance [26].Before adding a mass term to the action, which breaks the gauge symmetry, the massless theory should be well defined, in the sense that a well defined partition function must exist.In gauge field theory this is achieved by gauge fixing the action.This procedure (gauge fixing the massless action first, adding a breaking mass term after) leads to a gauge theory of LMG describing five DOF for the massive graviton, with a good massless limit and without the vDVZ discontinuity.In the FP theory, instead, the mass term is added directly to the invariant action and effectively plays the double and unnatural role of both mass term and gauge fixing.Basically, this is the reason why the FP theory displays an ill-defined massless limit, which shows itself in a divergent propagator and, physically, in the vDVZ discontinuity with GR.We stress that any theory of LMG should display a good massless limit since the phenomenological limits on the mass of the graviton are such that the mass term must be seen as a small perturbation of massless LG and, equivalently, as a small breaking of diffeomorphism invariance.

Summary
This paper is organized as follows: in Section 2 the gauge fixed massive theory is presented.We also review the FP theory and the problems related to its massless limit which motivate this work, and we remark that, curiously, the original 1939 paper by Fierz and Pauli does not exactly describe what is generally known as the FP theory.In Section 3 we show, in a gauge independent way, that the propagating DOF of the ten components of the rank-2 symmetric tensor describing the massive graviton are indeed five, and depend on one mass parameter only, in agreement with the FP theory.In Section 4 we compute the propagator of the theory, we show that its massless limit is regular and that the same observable which is used to unveil the vDVZ discontinuity presents an unique pole in the same mass parameter which characterizes the five DOF of the theory.In Section 5 we show the absence of the vDVZ discontinuity, for any gauge choice.Our results are summarized and discussed in the concluding Section 6.

The massive action
The weak field expansion of GR around the flat Minkowskian background η µν = diag(−1, 1, 1, 1) is given by the LG action where h µν (x) is a symmetric rank-2 tensor field representing the graviton, and h(x) ≡ η µν h µν (x) is its trace.The action S LG [h] (2.1) is the most general functional invariant under the infinitesimal diffeomorphism transformation where ξ µ (x) is a local vector parameter.The transformation (2.2) represents the gauge symmetry of the action S LG (2.1).

The Fierz-Pauli theory
The most general mass term which can be added to the invariant action S LG (2.1), respecting Lorentz invariance and power counting, is where m 2 1 and m 2 2 are massive parameters.The presence of a mass term breaks the diffeomorphism invariance (2.2), as usual in any gauge field theory.It can be shown (see for instance [25]) that the action S = S LG + S m .
(2.4) describes the propagation of five DOF only if otherwise a sixth ghost mode with negative energy appears [21], and the theory does not describe a massive graviton.The choice (2.5) is generally referred to as FP tuning, and the FP theory is defined by the action Following [11], we now show that the theory described by the action S F P (2.6) does indeed display five DOF.The Equations of Motion (EOM) obtained from (2.6) read which, saturated with ∂ µ , yield the constraint Plugging (2.8) into (2.7)we get Saturating (2.9) with η µν we find h = 0 , (2.10) which, together with (2.8), implies Therefore, the EOM (2.7) imply the following set of equations Eq. (2.12) is the Klein-Gordon equation for the field h µν (x), while (2.13) and (2.14) represent five constraints (transversality and tracelessness) which reduce the ten independent components of h µν to five.These five components carry the five massive DOF of the graviton.

Problems with the Fierz-Pauli theory
The propagator of the FP theory (2.6) is where P µν is the transverse massive projector defined as A crucial remark is that the propagator (2.15) exists only thanks to the presence of the mass term (2.3) on the FP point (2.5).As a consequence of this fact, it is apparent from (2.16) that the FP theory has a divergent massless limit.Moreover, the massless limit of the FP theory is flawed by the vDVZ discontinuity [17,18], which basically consists in the fact that the correlator involving two energy-momentum tensors, computed in the FP theory in the limit m 2 1 → 0, does not match the GR prediction.We sketch here the proof (for details see [11,19]).The gravitational interaction between two non relativistic energy-momentum tensors T (1) µν and T (2)  µν (which are conserved, i.e. p ν T (1) µν = p ν T (2) µν = 0) is described by the introduction in the action of an interaction term of the type where T µν denotes a generic energy-momentum tensor, coupled to h µν (x) through a constant λ, which we call λ LG and λ F P for LG and FP theory, respectively.Therefore, the interaction strength between T (1) µν and T (2 µν can be computed by contraction with the propagator of the graviton.In LG, the propagator G µν,αβ LG is obtained from the LG action (2.1) after a gauge fixing, as done in [11].The FP propagator G µν,αβ F P , on the other hand, is given by (2.15).In the non relativistic limit, only the 00-components of the energy-momentum tensors are non negligible.Two cases are considered: in the first, T µν and T (2) µν are associated with massive objects and therefore have non vanishing trace.In the second case, T µν still has a non vanishing trace whereas T (2) µν is traceless, representing, for instance, electromagnetic radiation (e.g.light).Concerning the first case, the interaction strength in LG is λ2 LG LG while, in the second, i.e. when T µν is traceless, we have LG LG On the other hand, the interaction strengths corresponding to (2.18) and (2.19) obtained using the FP propagator (2.15), in the massless limit m 1 → 0 are, respectively The vDVZ discontinuity arises from the observation that it is not possible to choose an unique FP coupling λ F P which matches both (2.LG matches the first pair but not the second.The interaction strength between a massive body and an electromagnetic wave in the massless limit of the FP theory turns out to be 3 4 of the LG prediction, hence the discontinuity.

Back to the origins
In [27] Fierz studied the general relativistic theory of force-free particles with any (integer or halfinteger) spin f , described by symmetric "world" tensors A ik...l (x) 2 .The starting point in [27].The starting point in [27] is the request that the fields A ik...l (x) satisfy the massive wave equation (which Fierz called of the "Schrödinger-Gordon type") where κ is a constant with the dimension of the inverse of a length (to which corresponds the mass m = hκ c ).In addition, the fields A ik...l (x) are asked to satisfy the additional constraints (called "secondary conditions") A ii...l = 0 (2.23) As explained in [27], the secondary conditions (2.23) and (2.24) are introduced in order to ensure that only particles with the spin f and not also those of smaller spins could be assigned to the tensor field, and that if the field A ik...l (x) describes a particle with spin f satisfying the massive "Schrödinger-Gordon" wave equation (2.22), then the number of linearly independent plane waves is 2f + 1, which differ by the orientation of the spin.Shortly after, in [16] Fierz and Pauli applied the general formalism described in [27] to the particular case of spin f = 2, simply taken as an example.Hence, they looked for the theory of a symmetric tensor field A ij (x) satisfying the massive wave equation together with the secondary conditions

26)
Imposing the condition of tracelessness (2.26) by hand, they derived the above equations (2.25) and (2.27) from the EOM of the following Lagrangian where an additional scalar field C(x) is introduced, and a i , i = 1, 2, 3 are dimensionless constant parameters.They showed that the desired results (2.25) and (2.27) are obtained if In (2.28) the introduction of the scalar field C(x) is an artifice which enables one to derive the transversality condition (2.27).In a more modern language, we would call this scalar field a Nakanishi-Lautrup Lagrange multiplier [28,29] introduced to implement the gauge condition (2.27).
It is therefore interesting (and surprising !) to remark that, in the original article [16], Fierz and Pauli wrote an action for the graviton which included both a mass term and a gauge fixing term, implemented by means of the Lagrange multiplier C(x).It appears that the problem of treating the mass term as a gauge fixing was absent in the original formulation.Moreover, the massless sector of the lagrangian L orig (2.28) does not coincide with the LG action S LG (2.1), since the trace A ii (x) (a.k.a.h(x)) is set to zero a priori.As a consequence, the only mass parameter appearing in (2.28) is the one coupled to A ij A ij (a.k.a.h µν h µν ), which therefore corresponds to m 2 1 in (2.3).It seems that, despite the fact that the original FP approach to LMG was realized in term of a gauge fixed action, later developments left the gauge fixing term behind, leaving room for the divergent massless limit and for the vDVZ discontinuity.
In this paper we adopt a conservative policy, recovering the original FP approach to LMG, with a few important differences.

Adding masses to the gauge fixed action
One of the reasons of the divergent massless limit of the FP theory is that the mass term S m (2.3) serves as gauge fixing as well, in the sense that its presence is necessary to define the propagator (2.15), as shown in Section 2.2.The standard way to proceed in gauge field theory, instead, is first to gauge fix the invariant action, in order to have a well defined partition function Z[J].This latter generates all the correlation functions of the theory, starting from the 2-points green function, a.k.a. the propagator.Only after having obtained that, the fields of the theory can be given masses through various procedures.In LMG this is achieved by adding to the gauge fixed action the mass term S m (2.3) [30].So, let us first proceed by gauge fixing the invariant action S LG [h] (2.1).The gauge field is represented by a rank-2 symmetric tensor h µν (x).Hence, the most general covariant gauge fixing is which, by means of the usual Faddeev-Popov (ΦΠ) exponentiation [31], yields the gauge fixing action term (2.32) Notice that the gauge fixing term (2.32) depends on two gauge parameters k and κ, which play different roles.In fact, k determines how the gauge fixing condition (2.31) is enforced.It can be seen as a kind of primary gauge fixing parameter, which corresponds to the standard gauge fixing parameter of Yang-Mills theory: k = 0 corresponds to the Landau gauge, for instance.On the other hand, the parameter κ fine-tunes the class of gauge fixing identified by k.It plays a secondary role.
As an example, the harmonic, or Lorenz, gauge is obtained with the choice κ = −1/2.Hence, it makes sense to talk about harmonic-Landau gauge, for instance, meaning by that the choice k = 0 and κ = −1/2.Once the action S LG [h] (2.1) has been gauge fixed by the gauge fixing term (2.32), we can add the mass term (2.3), so that our starting point for a theory of LMG is given by the action The action (2.33) is the starting, rather than the arrival, point, because the road ahead of us is still long.We have indeed to face the problems which affect the FP theory: in particular the massless limit and the absence of the vDVZ discontinuity.But, first, we have to deal with the main feature of LMG: the five DOF which must characterize a spin-2 massive particle.Our comparison is the FP theory, which reaches this goal with one mass parameter only, because of the FP tuning (2.5).The action (2.33), instead, depends on two masses, for now.This will be done in the next Section.

Degrees of freedom
A realistic theory of MG needs five propagating massive DOF.The easiest way to see this is to notice that the massive graviton is a spin S = 2 particle, which displays 2S + 1 independent components.Hence, given that the graviton is described by a symmetric rank-2 tensor h µν (x), only five out of its ten components correspond to physical DOF.Therefore, a necessary condition for a gauge theory of a massive rank-2 symmetric tensor to be promoted to a theory of LMG is to recover the five linear equations represented by the constraints of tracelessness (2.26) (or (2.14)) and of transversality (2.27) (or (2.13)), in order to lower the number of independent components of h µν (x) from ten to five.The realization of this necessary condition, in the framework of a well defined gauge field theory of a symmetric rank-2 tensor, is the main aim of this paper.To reach our goal, we shall now consider the EOM of the action (2.33) and we shall manage to restrict the mass parameters m 2 1 and m 2 2 to the cases in which we can find enough constraints to ensure the propagation of five massive DOF.Moreover, in order to deal with a physically consistent theory, we must also require that the propagating DOF do not depend on the gauge parameters k and κ, allowing instead, of course, a dependence on the mass parameters m 2 1 and m 2 2 .The action S LM G (2.33) in momentum space reads: where the kinetic operator Ω is In order to study the propagating DOF, we saturate the EOM (3.5) with η µν and e µν (3.4), to get From these two equations we deduce that, if m 2 1 = 0 (we will consider the case m 2 1 = 0 later), the only solution is (3.9) since (3.6) and (3.7) form a homogeneous system of two linear equations which has a non trivial solution only if the determinant of the coefficients matrix vanishes.For m 2 1 = 0 this determinant cannot vanish, independently of the choice of κ, k and m 2 2 .Substituting (3.8) and (3.9) into the EOM (3.5) we get which, saturated with p ν and using again (3.9), yields The above equation is satisfied if but the requirement that the physical masses should not depend on the gauge parameters, implies that the only allowable solution of (3.11) is (3.13).This, together with (3.8), gives the five constraints Substituting m 2 1 = 0 and h = 0 in the EOM (3.5) we notice that all the dependence on the mass parameters vanishes.The resulting EOM would therefore describe the propagation of a massless field or, alternatively, the propagation of a field with a mass dependent on the gauge parameters k and κ.Both cases do not represent acceptable descriptions of a LMG theory, therefore we conclude that it must be To summarize, we showed that the gauge fixed action S LM G (2.33), with the mass term S m (2.3), displays five DOF, provided that m 2 1 = 0, for any value of m 2 2 and for arbitrary gauge parameters k and κ.

Propagators
The momentum space propagator G µν,αβ (p) is defined by the condition where Ω µν,αβ (p) is the kinetic operator introduced in (3.2) and I µν,ρσ is the rank-4 tensor identity (3.3).In order to find G µν,αβ (p) it is convenient to introduce the rank-2 projectors e µν (3.4)  • idempotency : • orthogonality of A, B, C and D : • contractions with E : where, after a lenghty but straightforward calculation, the coefficients are given by Similarly, we can expand the propagator G µν,αβ (p) : where t, û, v, ẑ and ŵ are functions of the momentum p and depend on the gauge parameters k and κ appearing in S gf (2.32) and on the masses m 2 1 and m 2 2 of the mass term S m (2.3).Solving the defining equation (4.1), we find : where we used the shorthand notation for the denominator The propagator (4.23) displays poles that depend on the gauge parameters k and κ and might even be tachyonic.This does not come as a surprise, given that the theory describes the dynamics of a rank-2 symmetric tensor field, and the pole structure of its propagator is more complicated than the usual scalar, spinor or vector cases.Neither it should be seen as a problem, since in the previous Section we proved that only five of the ten components of h µν represent independent DOF, which satisfy the massive wave equation of the Klein-Gordon type (3.16) which depends on the mass m 2 1 only, in agreement with the FP theory.Hence, looking at the whole propagator is neither helpful nor correct in order to identify the physical pole in this case, since we may allow for non physical poles located in non physical sectors of the propagator.Rather, what should be done in order to select the physical pole, is to look to the observables related to the propagator (4.23).An important example is the one already considered in the analysis of the vDVZ discontinuity: the scattering amplitude of light and a massive object, mediated by the gravitational interaction, which is responsible for the observed time delay in gravitational lensing.Formally, this observable can be traced back to the more general interaction amplitude of two conserved energy-momentum tensors T µν and T µν , of which the one corresponding to light (which in the following we choose to be T (2) µν ) is traceless : Substituting the propagator (4.23) into (4.30),we get T (1) µν tA µν,αβ + ûB µν,αβ + vC µν,αβ + ẑD µν,αβ + ŵE µν,αβ T ( The above expression contains contractions of the tensor projectors ( µν .These are greatly simplified thanks to the fact that the energy-momentum tensors are conserved: In particular: T (1) µν B µν,αβ T (2) αβ = 0 (4.33) Eqs. (4.32)-(4.36)can be further simplified by using the fact that T (2) µν is traceless (η µν T (2) µν = 0), which means that the only non-vanishing contraction left is which confirms that the theory describes a graviton with mass m 2 1 , in agreement with the massive wave equation (3.16), which was obtained by studying the EOM deriving from the action S LM G (2.33), and also with the FP theory, (2.12).Notice that the same argument may be repeated to show that also the scattering amplitude between two radiation-like objects, both described by conserved and traceless energy-momentum tensors, isolates (4.39) as the unique physical pole of the propagator (4.23).

Absence of the vDVZ discontinuity
In Section 2.2 we reviewed the vDVZ discontinuity, which concerns the mismatch between the gravitational interaction of two energy-momentum tensors in the pure massless theory described by the action S LG (2.1) and in the massless limit of FP theory S F P (2.6).In [32] it has been shown that in a particular gauge (namely the k = − 1 2 , κ = − 1 2 harmonic one) the theory described by the action S LM G (2.33) is free of the vDVZ discontinuity.In this Section we generalize the result of [32] to all possible gauge choices, i.e. for every value of the gauge parameters k and κ.With the propagator G µν,αβ (4.23) derived from the gauge fixed action S LM G (2.33), hence for generic k and κ gauge parameters, we can compute the interaction amplitudes between two non-relativistic conserved energy-momentum tensors T µν and compare them to the corresponding LG amplitudes, (2.18) and (2.19).As explained in Section 2.2, the interaction of the graviton field with an external energy-momentum tensor is obtained by means of the interaction term S int (2.17).To distinguish the gauge fixed theory S LM G (2.33) from the LG and FP theories, we denote its coupling constant with λ LM G .Following the steps described in [19], if neither T (1) µν are traceless, which corresponds to the scattering between two massive bodies, the resulting amplitude in the massless limit (m 2 1 , m 2 2 ) → 0 of the theory is whereas, if T µν is traceless, corresponding to the scattering between light and a massive body, like in the gravitational lensing, we have It is straightforward to see that if we set the LG amplitudes (2.18) and (2.19) exactly match (5.1) and (5.2) respectively, which proves the absence of the vDVZ discontinuity in a gauge independent way.We stress that, although the massive propagator G µν,αβ (4.23) depends on the gauge parameters k and κ through the coefficients (4.24)-(4.28), the massless limit of the interaction amplitudes (5.1) and (5.2) is gauge independent and coincides with the LG prediction.This is related to the fact that the massless limit of the action S LM G (2.33) is the LG gauge fixed action where S LG and S gf are given by (2.1) and (2.32), respectively.Hence, in the massless limit, the propagators of the two theories coincide.This is the basic reason of the absence of the vDVZ discontinuity.In other words, the vDVZ discontinuity is a direct consequence of the structure of the FP action, where the mass term S m (2.3) is added directly to the invariant action S LG (2.1), therefore acting effectively as a gauge fixing term.Without the mass term, the invariant action does not have a propagator.Therefore, the mass parameters, which also play the role of gauge fixing parameters, cannot be physical.In the massless limit, the FP theory is not dynamical (i.e.does not have a propagator), and therefore it is not surprising that a discontinuity arises.Throughout this paper, we repeatedly claimed that the vDVZ discontinuity may be related to the nature of the FP mass term as a gauge fixing term.The massive gravity propagator in (2.15) is not well defined in the zero mass limit due to the lack of gauge fixing in the massless theory.However, the vDZV discontinuity appears in the analysis of the propagator contracted with energy-momentum tensors like in (4.30).This quantity is gauge invariant in linearized GR (for conserved sources).
Hence one might expect that the ill defined terms in the zero mass limit of massive FP propagator drop out of equations like (4.30) which contribute to the vDVZ discontinuity (this quantity indeed remains finite in the zero mass limit).The apparent contradiction between this observation and the explanation provided in this paper is explained by observing that when we say that the mass term S m (2.3) at the FP point (2.5) has the role of a gauge fixing in the FP theory, we are only referring to the fact that it enables us to invert the quadratic part of the action.In fact, not all quadratic terms which make it possible to find a propagator are also legitimate gauge fixing terms.The gauge fixing procedure is much more than this.Indeed, it consists in limiting the functional integral in order to choose one representative for each gauge orbit (modulo Gribov copies).This is realized through the highly non trivial ΦΠ trick.In the case of an abelian theory such as linearized gravity, where the ghost sector decouples, the ΦΠ gauge fixing is realized by adding to the action the square of a functional F µ [h αβ ], where F must carry the same number of Lorentz indices as the gauge parameter of the theory, defined in our case in (2.2).Since the mass term cannot be written as the square of any such functional, it follows that S m is not a true gauge fixing term and therefore no contradiction arises when comparing the linearized GR and FP propagators.

Summary of results and Discussion
In this paper we presented a gauge field theory of a massive symmetric tensor field describing a massive spin-2 particle.We summarize here the main features which allow to interpret this theory as an alternative to the FP theory of LMG. 4. continuity with LG Any candidate theory for LMG must display a finite massless limit, which implies a non divergent propagator, and, in this limit, it should also provide physical predictions in agreement with GR.Physically, LMG should represent only a small correction of GR, therefore the mass term (2.3) should be seen as a perturbation.In particular, the effects of a mass of the graviton should become relevant only at very large distances, while at smaller scales GR predictions must be restored.As discussed in Section 2.2, a serious flaw of the FP theory of LMG is the vDVZ discontinuity [17,18].This well known issue is fixed by the Vainshtein mechanism [11,20,25], which recovers the non linearities of GR in order to shield the effect of the extra scalar DOF introduced by the Stückelberg formalism [33] in the massless limit of the FP theory.Our gauge fixed theory of LMG restores the continuity with GR in a more natural way, without the introduction of additional fields, cutoffs or non linearities of the theory.In [32] it has been shown that in the particular (k = κ = − 1 2 ) harmonic gauge, the theory described by the action S LM G (2.33) does not display the vDVZ discontinuity.In Section 5 this result has been generalized to any gauge choice.The fact that the gauge fixed massive theory described by the action S LM G (2.33) is not affected by the vDVZ discontinuity with GR encourages to believe that the way of approaching LMG presented in this paper is correct.Moreover, although the absence of the vDVZ discontinuity is a specific test and it does not prove in general that every measurable quantity is continuous with GR, the way the theory itself was constructed might imply this result.In fact, the evaluation of any observable in quantum field theory needs a gauge fixing in order to define a propagator, which is guaranteed by the ΦΠ procedure, adopted in this paper.The key feature of our approach is that the massive action S LM G (2.33) in the zero-mass limit becomes exactly the ΦΠ gauge fixed LG action, i.e. the one we would have used to calculate the propagator and every other quantity in massless LG as well.Therefore, while it does not come as a surprise that the gravitational coupling between non relativistic matter and light turns out to be continuous with GR, we may also infer that every other observable is, indeed, continuous with LG, in the massless limit.Moreover, the striking difference between the FP theory of LMG, where the vDVZ discontinuity is present, and the massive ΦΠ gauge fixed theory described in this paper, where it is absent, gives evidence of the fact that, as anticipated, the FP theory is not a sub-case of our theory despite the fact that our approach contains the FP tuning (2.5) ;

role of the mass parameters
We stress as a remarkable fact that, although by means of a completely different approach, we recover here the main result of the FP theory, i.e. that only the mass parameter m 2 1 associated to h µν h µν in the action S m (2.3) appears in the propagation of the physical DOF through the massive wave equation (3.16), which is the starting request (2.25) of the original FP theory [16].Therefore, the mass of the graviton should be identified by m 2 1 .

2 ,
propagation of five DOF.Conditions (3.14) and (3.15) inserted into the EOM (3.5) give the massive propagation of h µν (p 2 + m 2 1 ) hµν = 0 , (3.16) which, interestingly, does not depend on the mass parameter m 2 .If, on the other hand, m 2 1 = 0, the two equations (3.6) and (3.7) have non trivial solutions only if κ = −1 ; k = − 3 In terms of the projectors e µν and d µν the operators X µν,αβ read µα e νβ + d µβ e να + e µα d νβ + e µβ d να ) • decomposition of the identity I µν,αβ (3.3) : A µν,αβ + B µν,αβ + C µν,αβ + D µν,αβ = I µν,αβ ; (4.11) µν,αβ = tA µν,αβ + uB µν,αβ + vC µν,αβ + zD µν,αβ + wE µν,αβ , (4.17) 1. finite massless limitThe theory described by the action S LM G(2.33)yields the propagator G µν,αβ (4.23) which is regular in the limits (m 2 1 , m 2 2 ) → 0 ; 2. five massive DOF In Section 3 we showed that only five of the ten components of the massive symmetric tensor h µν (x) actually represent DOF, as required for a massive graviton .The remarkable fact is that the mass associated with the propagation of the five DOF is m 2 1 , while the value of m 2 2 is irrelevant, under this respect.Although obtained in a quite different way, this result is in agreement with the FP theory, which is characterized by one mass parameter only ; 3. gauge independence In gauge field theory, any claim concerning physical properties should not depend on a particular gauge choice.In our case, this requirement is met in the determination of the physical DOF of the theory: none of the steps leading to the five constraints represented by (3.14) and (3.15) rely on particular values of the gauge parameters k and κ ;