BREAKING NEWS
Massive gravity

## Summary

In theoretical physics, massive gravity is a theory of gravity that modifies general relativity by endowing the graviton with a nonzero mass. In the classical theory, this means that gravitational waves obey a massive wave equation and hence travel at speeds below the speed of light.

Massive gravity has a long and winding history, dating back to the 1930s when Wolfgang Pauli and Markus Fierz first developed a theory of a massive spin-2 field propagating on a flat spacetime background. It was later realized in the 1970s that theories of a massive graviton suffered from dangerous pathologies, including a ghost mode and a discontinuity with general relativity in the limit where the graviton mass goes to zero. While solutions to these problems had existed for some time in three spacetime dimensions,[1][2] they were not solved in four dimensions and higher until the work of Claudia de Rham, Gregory Gabadadze, and Andrew Tolley (dRGT model) in 2010.

One of the very early massive gravity theories was constructed in 1965 by Ogievetsky and Polubarinov (OP).[3] Despite the fact that the OP model coincides with the ghost-free massive gravity models rediscovered in dRGT, the OP model has been almost unknown among contemporary physicists who work on massive gravity, perhaps because the strategy followed in that model was quite different from what is generally adopted at present.[4] Massive dual gravity to the OP model[5] can be obtained by coupling the dual graviton field to the curl of its own energy-momentum tensor.[6][7] Since the mixed symmetric field strength of dual gravity is comparable to the totally symmetric extrinsic curvature tensor of the Galileons theory, the effective Lagrangian of the dual model in 4-D can be obtained from the Faddeev–LeVerrier recursion, which is similar to that of Galileon theory up to the terms containing polynomials of the trace of the field strength.[8][9] This is also manifested in the dual formulation of Galileon theory.[10][11]

The fact that general relativity is modified at large distances in massive gravity provides a possible explanation for the accelerated expansion of the Universe that does not require any dark energy. Massive gravity and its extensions, such as bimetric gravity,[12] can yield cosmological solutions which do in fact display late-time acceleration in agreement with observations.[13][14][15]

Observations of gravitational waves have constrained the Compton wavelength of the graviton to be λg > 1.6×1016 m, which can be interpreted as a bound on the graviton mass mg < 7.7×10−23 eV/c2.[16] Competitive bounds on the mass of the graviton have also been obtained from solar system measurements by space missions such as Cassini and MESSENGER, which instead give the constraint λg > 1.83×1016 m or mg < 6.76×10−23 eV/c2.[17]

## Linearized massive gravity

At the linear level, one can construct a theory of a massive spin-2 field ${\displaystyle h_{\mu \nu }}$ propagating on Minkowski space. This can be seen as an extension of linearized gravity in the following way. Linearized gravity is obtained by linearizing general relativity around flat space, ${\displaystyle g_{\mu \nu }=\eta _{\mu \nu }+M_{\mathrm {Pl} }^{-1}h_{\mu \nu }}$, where ${\displaystyle M_{\mathrm {Pl} }=(8\pi G)^{-1/2}}$ is the Planck mass with ${\displaystyle G}$ the gravitational constant. This leads to a kinetic term in the Lagrangian for ${\displaystyle h_{\mu \nu }}$ which is consistent with diffeomorphism invariance, as well as a coupling to matter of the form

${\displaystyle h^{\mu \nu }T_{\mu \nu },}$

where ${\displaystyle T_{\mu \nu }}$ is the stress–energy tensor. This kinetic term and matter coupling combined are nothing other than the Einstein–Hilbert action linearized about flat space.

Massive gravity is obtained by adding nonderivative interaction terms for ${\displaystyle h_{\mu \nu }}$. At the linear level (i.e., second order in ${\displaystyle h_{\mu \nu }}$), there are only two possible mass terms:

${\displaystyle {\mathcal {L}}_{\mathrm {int} }=ah^{\mu \nu }h_{\mu \nu }+b\left(\eta ^{\mu \nu }h_{\mu \nu }\right)^{2}.}$

Fierz and Pauli[18] showed in 1939 that this only propagates the expected five polarizations of a massive graviton (as compared to two for the massless case) if the coefficients are chosen so that ${\displaystyle a=-b}$. Any other choice will unlock a sixth, ghostly degree of freedom. A ghost is a mode with a negative kinetic energy. Its Hamiltonian is unbounded from below and it is therefore unstable to decay into particles of arbitrarily large positive and negative energies. The Fierz–Pauli mass term,

${\displaystyle {\mathcal {L}}_{\mathrm {FP} }=m^{2}\left(h^{\mu \nu }h_{\mu \nu }-\left(\eta ^{\mu \nu }h_{\mu \nu }\right)^{2}\right)}$

is therefore the unique consistent linear theory of a massive spin-2 field.

## The vDVZ discontinuity

In the 1970s Hendrik van Dam and Martinus J. G. Veltman[19] and, independently, Valentin I. Zakharov[20] discovered a peculiar property of Fierz–Pauli massive gravity: its predictions do not uniformly reduce to those of general relativity in the limit ${\displaystyle m\to 0}$. In particular, while at small scales (shorter than the Compton wavelength of the graviton mass), Newton's gravitational law is recovered, the bending of light is only three quarters of the result Albert Einstein obtained in general relativity. This is known as the vDVZ discontinuity.

We may understand the smaller light bending as follows. The Fierz–Pauli massive graviton, due to the broken diffeomorphism invariance, propagates three extra degrees of freedom compared to the massless graviton of linearized general relativity. These three degrees of freedom package themselves into a vector field, which is irrelevant for our purposes, and a scalar field. This scalar mode exerts an extra attraction in the massive case compared to the massless case. Hence, if one wants measurements of the force exerted between nonrelativistic masses to agree, the coupling constant of the massive theory should be smaller than that of the massless theory. But light bending is blind to the scalar sector, because the stress-energy tensor of light is traceless. Hence, provided the two theories agree on the force between nonrelativistic probes, the massive theory would predict a smaller light bending than the massless one.

## Vainshtein screening

It was argued by Vainshtein[21] two years later that the vDVZ discontinuity is an artifact of the linear theory, and that the predictions of general relativity are in fact recovered at small scales when one takes into account nonlinear effects, i.e., higher than quadratic terms in ${\displaystyle h_{\mu \nu }}$. Heuristically speaking, within a region known as the Vainshtein radius, fluctuations of the scalar mode become nonlinear, and its higher-order derivative terms become larger than the canonical kinetic term. Canonically normalizing the scalar around this background therefore leads to a heavily suppressed kinetic term, which damps fluctuations of the scalar within the Vainshtein radius. Because the extra force mediated by the scalar is proportional to (minus) its gradient, this leads to a much smaller extra force than we would have calculated just using the linear Fierz–Pauli theory.

This phenomenon, known as Vainshtein screening, is at play not just in massive gravity, but also in related theories of modified gravity such as DGP and certain scalar-tensor theories, where it is crucial for hiding the effects of modified gravity in the solar system. This allows these theories to match terrestrial and solar-system tests of gravity as well as general relativity does, while maintaining large deviations at larger distances. In this way these theories can lead to cosmic acceleration and have observable imprints on the large-scale structure of the Universe without running afoul of other, much more stringent constraints from observations closer to home.

## The Boulware–Deser ghost

As a response to Freund–Maheshwari–Schonberg finite-range gravity model,[22] and around the same time as the vDVZ discontinuity and Vainshtein mechanism were discovered, David Boulware and Stanley Deser found in 1972 that generic nonlinear extensions of the Fierz–Pauli theory reintroduced the dangerous ghost mode;[23] the tuning ${\displaystyle a=-b}$ which ensured this mode's absence at quadratic order was, they found, generally broken at cubic and higher orders, reintroducing the ghost at those orders. As a result, this Boulware–Deser ghost would be present around, for example, highly inhomogeneous backgrounds.

This is problematic because a linearized theory of gravity, like Fierz–Pauli, is well-defined on its own but cannot interact with matter, as the coupling ${\displaystyle h^{\mu \nu }T_{\mu \nu }}$ breaks diffeomorphism invariance. This must be remedied by adding new terms at higher and higher orders, ad infinitum. For a massless graviton, this process converges and the end result is well-known: one simply arrives at general relativity. This is the meaning of the statement that general relativity is the unique theory (up to conditions on dimensionality, locality, etc.) of a massless spin-2 field.

In order for massive gravity to actually describe gravity, i.e., a massive spin-2 field coupling to matter and thereby mediating the gravitational force, a nonlinear completion must similarly be obtained. The Boulware–Deser ghost presents a serious obstacle to such an endeavor. The vast majority of theories of massive and interacting spin-2 fields will suffer from this ghost and therefore not be viable. In fact, until 2010 it was widely believed that all Lorentz-invariant massive gravity theories possessed the Boulware–Deser ghost.[24]

## Ghost-free massive gravity

In 2010 a breakthrough was achieved when de Rham, Gabadadze, and Tolley constructed, order by order, a theory of massive gravity with coefficients tuned to avoid the Boulware–Deser ghost by packaging all ghostly (i.e., higher-derivative) operators into total derivatives which do not contribute to the equations of motion.[25][26] The complete absence of the Boulware–Deser ghost, to all orders and beyond the decoupling limit, was subsequently proven by Fawad Hassan and Rachel Rosen.[27][28]

The action for the ghost-free de Rham–Gabadadze–Tolley (dRGT) massive gravity is given by[29]

${\displaystyle S=\int d^{4}x{\sqrt {-g}}\left(-{\frac {M_{\mathrm {Pl} }^{2}}{2}}R+m^{2}M_{\mathrm {Pl} }^{2}\sum _{n=0}^{4}\alpha _{n}e_{n}(\mathbb {K} )+{\mathcal {L}}_{\mathrm {m} }(g,\Phi _{i})\right),}$

or, equivalently,

${\displaystyle S=\int d^{4}x{\sqrt {-g}}\left(-{\frac {M_{\mathrm {Pl} }^{2}}{2}}R+m^{2}M_{\mathrm {Pl} }^{2}\sum _{n=0}^{4}\beta _{n}e_{n}(\mathbb {X} )+{\mathcal {L}}_{\mathrm {m} }(g,\Phi _{i})\right).}$

The ingredients require some explanation. As in standard general relativity, there is an Einstein–Hilbert kinetic term proportional to the Ricci scalar ${\displaystyle R}$ and a minimal coupling to the matter Lagrangian ${\displaystyle {\mathcal {L}}_{\mathrm {m} }}$, with ${\displaystyle \Phi _{i}}$ representing all of the matter fields, such as those of the Standard Model. The new piece is a mass term, or interaction potential, constructed carefully to avoid the Boulware–Deser ghost, with an interaction strength ${\displaystyle m}$ which is (if the nonzero ${\displaystyle \beta _{i}}$ are ${\displaystyle {\mathcal {O}}(1)}$) closely related to the mass of the graviton.

Th principle of gauge-invariance renders redundant expressions in any field theory provided with its corresponding gauge(s). For example, in the massive spin-1 Proca action, the massive part in the Lagrangian ${\displaystyle {\frac {1}{2}}mA_{\mu }A^{\mu }}$ breaks the ${\displaystyle U(1)}$ gauge-invariance. However, the invariance is restored by introducing the transformations:

${\displaystyle A_{\mu }\to A_{\mu }+\partial _{\mu }\pi }$. The same can be done for massive gravity by following Arkani-Hamed, Georgi and Schwartz effective field theory for massive gravity.[30] The absence of vDVZ discontinuity in this approach motivated the development of dRGT resummation of massive gravity theory as follows.[26]

The interaction potential is built out of the elementary symmetric polynomials ${\displaystyle e_{n}}$ of the eigenvalues of the matrices ${\textstyle \mathbb {K} =\mathbb {I} -{\sqrt {g^{-1}f}}}$ or ${\textstyle \mathbb {X} ={\sqrt {g^{-1}f}}}$, parametrized by dimensionless coupling constants ${\displaystyle \alpha _{i}}$ or ${\displaystyle \beta _{i}}$, respectively. Here ${\textstyle {\sqrt {g^{-1}f}}}$ is the matrix square root of the matrix ${\displaystyle g^{-1}f}$. Written in index notation, ${\displaystyle \mathbb {X} }$ is defined by the relation ${\displaystyle X^{\mu }{}_{\alpha }X^{\alpha }{}_{\nu }=g^{\mu \alpha }f_{\nu \alpha }.}$ We have introduced a reference metric ${\displaystyle f_{\mu \nu }}$ in order to construct the interaction term. There is a simple reason for this: it is impossible to construct a nontrivial interaction (i.e., nonderivative) term from ${\displaystyle g_{\mu \nu }}$ alone. The only possibilities are ${\displaystyle g^{\mu \alpha }g_{\alpha \nu }=\delta _{\nu }^{\mu }}$ and ${\displaystyle \operatorname {det} g}$, both of which lead to a cosmological constant term rather than a bona fide interaction. Physically, ${\displaystyle f_{\mu \nu }}$ corresponds to the background metric around which fluctuations take the Fierz–Pauli form. This means that, for instance, nonlinearly completing the Fierz–Pauli theory around Minkowski space given above will lead to dRGT massive gravity with ${\displaystyle f_{\mu \nu }=\eta _{\mu \nu }}$, although the proof of absence of the Boulware–Deser ghost holds for general ${\displaystyle f_{\mu \nu }}$.[31]

The reference metric transforms like a metric tensor under diffeomorphism ${\displaystyle f_{\mu \nu }\to f'_{\mu \nu }(X(x))\equiv f_{\alpha \beta }\partial _{\mu }X^{\alpha }\partial _{\nu }X^{\beta }.}$ Therefore ${\displaystyle [\mathbb {X} ^{2}]=X^{\mu }{}_{\alpha }X^{\alpha }{}_{\mu }=g^{\mu \alpha }f_{\mu \alpha }}$, and similar terms with higher powers, transforms as a scalar under the same diffeomorphism. For a change in the coordinates ${\displaystyle x_{\mu }\to x_{\mu }+\xi _{\mu }}$, we expand ${\displaystyle X^{\mu }=x^{\mu }-\phi ^{\mu }}$ with ${\displaystyle f_{\mu \nu }=\eta _{\mu \nu }}$ such that the perturbed metric becomes ${\displaystyle h_{\mu \nu }\to h'_{\mu \nu }=h_{\mu \nu }+\partial _{\mu }\phi _{\nu }+\partial _{\nu }\phi _{\mu }-\partial _{\mu }\phi ^{\alpha }\partial _{\nu }\phi _{\alpha }}$, while the potential-like vector transforms according to Stueckelberg trick as ${\displaystyle \phi _{\mu }\to \phi _{\mu }+\xi _{\mu }}$ such that the Stueckelberg field is defined as ${\displaystyle \phi ^{\mu }=\eta ^{\mu \nu }(A_{\nu }-\partial _{\nu }\pi )}$.[32] From the diffeomorphism, one can define another Stueckelberg matrix ${\displaystyle \mathbb {Y} _{~b}^{a}\equiv f_{bc}g^{\mu \nu }\partial _{\mu }X^{a}\partial _{\nu }X^{c}}$, where ${\displaystyle \mathbb {Y} _{~b}^{a}}$ and ${\displaystyle \mathbb {X} _{~b}^{a}}$ have the same eigenvalues.[33] Now, one considers the following symmetries:

• ${\displaystyle \delta h_{\mu \nu }=\partial _{\mu }\xi _{\nu }+\partial _{\nu }\xi _{\mu }+{\frac {2}{M_{\text{Pl}}}}{\mathcal {L}}_{\xi }h_{\mu \nu }}$,
• ${\displaystyle \delta A_{\mu }=\partial _{\mu }\pi -m\xi _{\mu }+{\frac {2}{M_{\text{Pl}}}}\xi ^{\alpha }\partial _{\alpha }A_{\mu }}$,
• ${\displaystyle \delta \pi =-m\pi }$,

such that the transformed perturbed metric becomes: ${\displaystyle h'_{\mu \nu }=h_{\mu \nu }+\partial _{\mu }A_{\nu }+\partial _{\nu }A_{\mu }-\partial _{\mu }A^{\alpha }\partial _{\nu }A_{\alpha }-\partial _{\mu }A^{\alpha }\partial _{\nu }\partial _{\alpha }\pi -\partial _{\mu }\partial _{\alpha }\pi \partial _{\nu }A^{\alpha }-2\partial _{\mu }\partial _{\nu }\pi -\partial _{\mu }\partial _{\alpha }\pi \ \partial _{\nu }\partial ^{\alpha }\pi .}$

The covariant form of these transformations are obtained as follows. If helicity-0 (or spin-0) mode ${\displaystyle \pi }$ is a pure gauge of unphysical Goldstone modes, with ${\displaystyle \Pi _{\mu \nu }=\nabla _{\mu }\nabla _{\nu }\pi }$,[34] the matrix ${\displaystyle {\displaystyle \mathbb {X} }}$ is a tensor function of the covariantization tensor ${\displaystyle {\displaystyle H_{\mu \nu }=\eta _{\mu \nu }+2{\Pi }_{\mu \nu }-\eta ^{\alpha \beta }{\Pi }_{\mu \alpha }{\Pi }_{\beta \nu }}=\eta _{\mu \nu }+h_{\mu \nu }-\eta _{ab}\nabla _{\mu }\phi ^{a}\nabla _{\nu }\phi ^{b}}$ of the metric perturbation ${\displaystyle {\displaystyle h_{\mu \nu }}}$ such that tensor ${\displaystyle H_{\mu \nu }}$ is Stueckelbergized by the field ${\displaystyle \phi ^{a}=A^{a}-\eta ^{a\mu }\nabla _{\mu }\pi }$.[35] Helicity-0 mode transforms under Galilean transformations ${\displaystyle \pi \to \pi +c+v_{\mu }x^{\mu }}$, hence the name "Galileons".[36] The matrix ${\displaystyle \mathbb {X} }$ is a tensor function of the covariantization tensor ${\displaystyle H_{\mu \nu }\equiv g_{\mu \nu }-f'_{\mu \nu }}$ of the metric perturbation ${\displaystyle h_{\mu \nu }}$ with components are given by ${\displaystyle {\displaystyle \mathbb {X} _{\mu \nu }=\eta _{\mu \nu }+2{\mathcal {K}}_{\mu \nu }-\eta ^{\alpha \beta }{\mathcal {K}}_{\mu \alpha }{\mathcal {K}}_{\beta \nu }}}$, where ${\textstyle {\mathcal {K}}_{\mu \nu }=\eta _{\mu \nu }-\left({\sqrt {\mathbb {X} }}\right)_{\mu \nu }=\eta _{\mu \nu }-{\sqrt {\eta _{\mu \nu }-H_{\mu \nu }}}}$ is the extrinsic curvature.[37]

Interestingly, the covariantization tensor was originally introduced by Maheshwari in a solo authored paper sequel to helicity-(${\displaystyle 2\oplus 0}$) Freund–Maheshwari–Schonberg finite-range gravitation model.[38] In Maheshwari's work, the metric perturbation obeys Hilbert-Lorentz condition ${\displaystyle m^{2}(\partial _{\nu }h_{\mu \nu }+q\partial _{\mu }h)=0}$ under the variation ${\displaystyle \delta ^{*}h_{\mu \nu }=\delta h_{\mu \nu }+\delta h_{\mu \nu }^{spin}=\partial _{\mu }\xi _{\nu }+\partial _{\nu }\xi _{\mu }+p\eta _{\mu \nu }h+\delta h_{\mu \nu }^{spin}}$ that is introduced in Ogievetsky–Polubarinov massive gravity, where ${\displaystyle p~{\text{and}}~q}$ are to be determined.[39] It is easy to notice the similarity between tensor ${\displaystyle \mathbb {X} }$ in dRGT and the tensor ${\displaystyle h_{(p)}^{\mu \nu }=(\eta ^{\mu \nu }-n\psi ^{\mu \nu })^{1/n}}$ in Maheshwari work once ${\displaystyle n=2}$ is chosen. Also Ogievetsky–Polubarinov model mandates ${\displaystyle n=-1/p}$, which means that in 4D ${\textstyle p=-1/n=-{\frac {1}{2}}}$, the variation ${\displaystyle \delta h_{\mu \nu }}$ is conformal.

The dRGT massive fields split into two helicity-2 ${\displaystyle h_{\mu \nu }}$, two helicity-1 ${\displaystyle A_{\mu }}$ and one helicity-0 ${\displaystyle \pi }$ degrees of freedom, just like those of Fierz-Pauli massive theory. However, the covariantization, together with the decoupling limit, guarantee that the symmetries of this massive theory are reduced to the symmetry of linearized general relativity plus that of ${\displaystyle U(1)}$ massive theory, while the scalar decouples. If ${\displaystyle v_{\mu }}$ is chosen to be divergenceless, i.e. ${\displaystyle \Box \pi =0}$, the decoupling limit of dRGT gives the known linearized gravity.[40] To see how that happens, expand the terms containing ${\displaystyle {\mathcal {K}}_{\mu \nu }}$ in the action in powers of ${\displaystyle H_{\mu \nu }}$, where ${\displaystyle H_{\mu \nu }}$ is expressed in terms of ${\displaystyle \phi ^{a}}$ fields like how ${\displaystyle h'_{\mu \nu }}$ is expressed in terms of ${\displaystyle A^{\mu }}$. The fields ${\displaystyle h_{\mu \nu }\ ,\ A_{\mu }\ ,\ \pi }$ are replaced by: ${\displaystyle {\tilde {h}}_{\mu \nu }=M_{\text{Pl}}h_{\mu \nu }\ ,\ {\tilde {A}}_{\mu }=M_{\text{Pl}}m\ A_{\mu }\ ,\ {\tilde {\pi }}=M_{\text{Pl}}m^{2}\pi \ ,\ {\hat {h}}_{\mu \nu }={\tilde {h}}_{\mu \nu }-\eta _{\mu \nu }{\tilde {\pi }}}$ . Then it follows that in the decoupling limit, i.e. when both ${\displaystyle M_{\text{Pl}}\to \infty \ ,m\to 0\ ,\ m^{2}M_{\text{Pl}}={\text{const.}}}$, the massive gravity Lagrangian is invariant under:

1. ${\displaystyle \delta h_{\mu \nu }=\partial _{\mu }\xi _{\nu }+\partial _{\nu }\xi _{\mu }}$ as in Linearized general theory of relativity,
2. ${\displaystyle \delta A_{\mu }=\partial _{\mu }\pi }$ as in Maxwell's electromagnetic theory,
3. ${\displaystyle \delta \pi =0}$.

In principle, the reference metric must be specified by hand, and therefore there is no single dRGT massive gravity theory, as the theory with a flat reference metric is different from one with a de Sitter reference metric, etc. Alternatively, one can think of ${\displaystyle f_{\mu \nu }}$ as a constant of the theory, much like ${\displaystyle m}$ or ${\displaystyle M_{\mathrm {Pl} }}$. Instead of specifying a reference metric from the start, one can allow it to have its own dynamics. If the kinetic term for ${\displaystyle f_{\mu \nu }}$ is also Einstein–Hilbert, then the theory remains ghost-free and we are left with a theory of massive bigravity,[12] (or bimetric relativity, BR) propagating the two degrees of freedom of a massless graviton in addition to the five of a massive one.

In practice it is unnecessary to compute the eigenvalues of ${\displaystyle \mathbb {X} }$ (or ${\displaystyle \mathbb {K} }$) in order to obtain the ${\displaystyle e_{n}}$. They can be written directly in terms of ${\displaystyle \mathbb {X} }$ as

{\displaystyle {\begin{aligned}e_{0}(\mathbb {X} )&=1,\\e_{1}(\mathbb {X} )&=[\mathbb {X} ],\\e_{2}(\mathbb {X} )&={\frac {1}{2}}\left([\mathbb {X} ]^{2}-\left[\mathbb {X} ^{2}\right]\right),\\e_{3}(\mathbb {X} )&={\frac {1}{6}}\left([\mathbb {X} ]^{3}-3[\mathbb {X} ]\left[\mathbb {X} ^{2}\right]+2\left[\mathbb {X} ^{3}\right]\right),\\e_{4}(\mathbb {X} )&=\det \mathbb {X} ,\end{aligned}}}

where brackets indicate a trace, ${\displaystyle [\mathbb {X} ]\equiv X^{\mu }{}_{\mu }}$. It is the particular antisymmetric combination of terms in each of the ${\displaystyle e_{n}}$ which is responsible for rendering the Boulware–Deser ghost nondynamical.

The choice to use ${\displaystyle \mathbb {X} }$ or ${\displaystyle \mathbb {K} =\mathbb {I} -\mathbb {X} }$, with ${\displaystyle \mathbb {I} }$ the identity matrix, is a convention, as in both cases the ghost-free mass term is a linear combination of the elementary symmetric polynomials of the chosen matrix. One can transform from one basis to the other, in which case the coefficients satisfy the relationship[29]

${\displaystyle \beta _{n}=(4-n)!\displaystyle \sum _{i=n}^{4}{\frac {(-1)^{i+n}}{(4-i)!(i-n)!}}\alpha _{i}.}$

The coefficients are of a characteristic polynomial that is in form of Fredholm determinant. They can also be obtained using Faddeev–LeVerrier algorithm.

## Massive gravity in the vierbein language

In 4D orthonormal tetrad frame, we have the bases:

{\displaystyle {\begin{aligned}e_{~\mu }^{0}&=(-1,0,0,0)\\e_{\mu }^{I}&=(0,e_{i}^{I})\end{aligned}}}

where the index ${\displaystyle i}$ is for the 3D spatial component of the ${\displaystyle \mu }$-non-orthonormal coordinates, and the index ${\displaystyle I}$ is for the 3D spatial components of the ${\displaystyle a}$-orthonormal ones. The parallel transport requires the spin connection ${\displaystyle e^{0\nu }\nabla _{e^{0\nu }}e_{~\mu }^{I}=0.}$ Therefore, the extrinsic curvature, that corresponds to ${\displaystyle {\mathcal {K}}_{\mu \nu }}$ in metric formalism, becomes

${\displaystyle K_{j}^{i}\equiv {\frac {1}{2}}\gamma ^{ik}\partial _{t}(\gamma _{kj})=e_{I}^{i}\partial _{t}(e_{j}^{I}),}$

where ${\displaystyle \gamma _{ij}}$ is the spatial metric as in the ADM formalism and initial value formulation.

If the tetrad conformally transforms as ${\displaystyle e_{i}^{I}\to {e'}_{i}^{I}\equiv a(t)~e_{i}^{I},}$ the extrinsic curvature becomes ${\displaystyle {K'}_{j}^{i}={\frac {a}{\dot {a}}}K_{j}^{i}=\delta _{j}^{i}-{\frac {a}{\dot {a}}}e_{I}^{i}\partial _{t}\!\!\left(e_{j}^{I}\right)}$, where from Friedmann equations ${\displaystyle {\frac {a}{\dot {a}}}={\frac {a}{\partial _{t}a}}\sim {\frac {1}{\sqrt {\Lambda }}}}$, and ${\displaystyle m\sim 1/{\sqrt {\Lambda }}}$ (despite it is controversial[41]), i.e. the extrinsic curvature transforms as ${\displaystyle K_{j}^{i}\to mK_{j}^{i}=\delta _{j}^{i}-m~e_{I}^{i}{\dot {e}}_{j}^{I}}$. This looks very similar to the matrix ${\displaystyle \mathbb {K} }$ or the tensor ${\displaystyle {\mathcal {K}}_{j}^{i}}$.

The dRGT was developed inspired by applying the previous technique to the 5D DGP model after considering the deconstruction of higher dimensional Kaluza-Klein gravity theories,[42] in which the extra-dimension(s) is/are replaced by series of N lattice sites such that the higher dimensional metric is replaced by a set of interacting metrics that depend only on the 4D components.[37]

The presence of a square-root matrix is somewhat awkward and points to an alternative, simpler formulation in terms of vierbeins. Splitting the metrics into vierbeins as

{\displaystyle {\begin{aligned}g_{\mu \nu }&=\eta _{ab}e^{a}{}_{\mu }e^{b}{}_{\nu }\\f_{\mu \nu }&=\eta _{ab}f^{a}{}_{\mu }f^{b}{}_{\nu }\end{aligned}}}

and then defining one-forms

{\displaystyle {\begin{aligned}\mathbf {e} ^{a}&=e^{a}{}_{\mu }dx^{\mu }\\\mathbf {f} ^{a}&=f^{a}{}_{\mu }dx^{\mu }\\\mathbf {i} ^{a}&=\delta ^{a}{}_{\mu }dx^{\mu }\end{aligned}}}

the ghost-free interaction terms in Hassan-Rosen bigravity theory can be written simply as (up to numerical factors)[43]

{\displaystyle {\begin{aligned}e_{0}(\mathbb {X} )\propto \epsilon _{abcd}\mathbf {e} ^{a}\wedge \mathbf {e} ^{b}\wedge \mathbf {e} ^{c}\wedge \mathbf {e} ^{d}\\e_{1}(\mathbb {X} )\propto \epsilon _{abcd}\mathbf {e} ^{a}\wedge \mathbf {e} ^{b}\wedge \mathbf {e} ^{c}\wedge \mathbf {f} ^{d}\\e_{2}(\mathbb {X} )\propto \epsilon _{abcd}\mathbf {e} ^{a}\wedge \mathbf {e} ^{b}\wedge \mathbf {f} ^{c}\wedge \mathbf {f} ^{d}\\e_{3}(\mathbb {X} )\propto \epsilon _{abcd}\mathbf {e} ^{a}\wedge \mathbf {f} ^{b}\wedge \mathbf {f} ^{c}\wedge \mathbf {f} ^{d}\\e_{4}(\mathbb {X} )\propto \epsilon _{abcd}\mathbf {f} ^{a}\wedge \mathbf {f} ^{b}\wedge \mathbf {f} ^{c}\wedge \mathbf {f} ^{d}\end{aligned}}}

In terms of vierbeins, rather than metrics, we can therefore see the physical significance of the ghost-free dRGT potential terms quite clearly: they are simply all the different possible combinations of wedge products of the vierbeins of the two metrics.

Note that massive gravity in the metric and vierbein formulations are only equivalent if the symmetry condition

${\displaystyle (e^{-1})_{a}{}^{\mu }f_{b\nu }=(e^{-1})_{b}{}^{\mu }f_{a\nu }}$

is satisfied. While this is true for most physical situations, there may be cases, such as when matter couples to both metrics or in multimetric theories with interaction cycles, in which it is not. In these cases the metric and vierbein formulations are distinct physical theories, although each propagates a healthy massive graviton.

The novelty in dRGT massive gravity is that it is a theory of gauge invariance under both local Lorentz transformations, from assuming the reference metric ${\displaystyle f_{\mu \nu }}$ equals the Minkowski metric ${\displaystyle \eta _{\mu \nu }}$, and diffeomorphism invariance, from the existence of the active curved spacetime ${\displaystyle g_{\mu \nu }}$. This is shown by rewriting the previously discussed Stueckelberg formalism in the vierbein language as follows.[44]

The 4D version of Einstein field equations in 5D is read

${\displaystyle G_{\mu \nu }n^{\mu }n^{\nu }={\frac {1}{2}}\left(R-K^{\mu \nu }K_{\mu \nu }+\left(K_{~\mu }^{\mu }\right)^{2}\right),}$

where ${\displaystyle n^{\mu }}$ is the vector normal to the 4D slice. Using the definition of massive extrinsic curvature ${\displaystyle mK_{j}^{i}=\delta _{j}^{i}-m~e_{I}^{i}{\dot {e}}_{j}^{I},}$ it is straightforward to see that terms containing extrinsic curvatures take the functional form ${\displaystyle (f^{a}-e^{a})\wedge (f^{b}-e^{b})\wedge e^{c}\wedge e^{d}}$ in the tetradic action.

Therefore, up to the numerical coefficients, the full dRGT action in its tensorial form is

${\displaystyle S={\frac {M_{\text{Pl}}^{2}}{2}}\int dx^{4}{\sqrt {g}}\left(R+2m^{2}[e_{2}({\mathcal {K}})+e_{3}({\mathcal {K}})+e_{4}({\mathcal {K}})]\right),}$

where the functions ${\displaystyle e_{i}({\mathcal {K}})}$ take forms similar to that of the ${\displaystyle e_{i}(\mathbb {X} )}$. Then, up to some numerical coefficients, the action takes the integral form

${\displaystyle S={\frac {M_{\text{Pl}}^{2}}{2}}\epsilon _{abcd}\int {\Big (}\mathbf {e} ^{a}\land \mathbf {e} ^{b}\land R^{cd}-m^{2}\left[\mathbf {e} ^{a}\land \mathbf {e} ^{b}\land \mathbf {e} ^{c}\land \mathbf {e} ^{d}+\mathbf {i} ^{a}\land \mathbf {e} ^{b}\land \mathbf {e} ^{c}\land \mathbf {e} ^{d}+\mathbf {i} ^{a}\land \mathbf {i} ^{b}\land \mathbf {e} ^{c}\land \mathbf {e} ^{d}+\mathbf {i} ^{a}\land \mathbf {i} ^{b}\land \mathbf {i} ^{c}\land \mathbf {e} ^{d}\right]{\Big )},}$

where the first term is the Einstein-Hilbert part of the tetradic Palatini action and ${\displaystyle \epsilon _{abcd}}$ is the Levi-Civita symbol.

As the decoupling limit guarantees that ${\displaystyle \Box \pi =0}$ and ${\displaystyle A_{\mu }\to x_{\mu }}$ by comparing ${\displaystyle X^{\mu }}$ to ${\displaystyle \phi ^{\mu }}$, it is legit to think of the tensor ${\displaystyle \partial _{\mu }\phi ^{a}=\delta _{~\mu }^{a}.}$ Comparing this with the definition of the 1-form ${\displaystyle \mathbf {i} ^{a},}$ one can define covariant components of frame field ${\displaystyle f_{~\mu }^{a}=\partial _{\mu }\phi ^{\nu }\delta _{~\nu }^{a},}$ i.e. ${\displaystyle e_{~\mu }^{a}={\frac {\partial x^{\nu }}{\partial \phi ^{\mu }}}\Lambda _{~b}^{a}e_{~\nu }^{b}}$, to replace the ${\displaystyle \mathbf {i} ^{a}}$ such that the last three interaction terms in the vierbein action becomes

${\displaystyle S=-{\frac {M_{\text{Pl}}^{2}}{2}}m^{2}\epsilon _{abcd}\int \left[\mathbf {f} ^{a}\land \left(\Lambda _{b'}^{b}\mathbf {e} ^{b'}\right)\land \left(\Lambda _{c'}^{c}\mathbf {e} ^{c'}\right)\land \left(\Lambda _{d'}^{d}\mathbf {e} ^{d'}\right)+\mathbf {f} ^{a}\land \mathbf {f} ^{b}\land \left(\Lambda _{c'}^{c}\mathbf {e} ^{c'}\right)\land \left(\Lambda _{d'}^{d}\mathbf {e} ^{d'}\right)+\mathbf {f} ^{a}\land \mathbf {f} ^{b}\land \mathbf {f} ^{c}\land \left(\Lambda _{d'}^{d}\mathbf {e} ^{d'}\right)\right].}$

This can be done because one is allowed to freely move the diffeomorphism transformations ${\displaystyle \partial _{\nu }\phi ^{\mu }}$ onto the reference vierbein through the Lorentz transformations ${\displaystyle \Lambda _{~b}^{a}}$. More importantly, the diffeomorphism transformations help manifesting the dynamics of the helicity-0 and helicity-1 modes, hence the easiness of gauging them away when the theory is compared with its version with the only ${\displaystyle U(1)}$ gauge transformations while the Stueckelberg fields are turned off.

One may wonder why the coefficients are dropped, and how to guarantee they are numerical with no explicit dependence of the fields. In fact this is allowed because the variation of the vierbein action with respect to the locally Lorentz transformed Stueckelberg fields yields this nice result.[44] Moreover, we can solve explicitly for the Lorentz invariant Stueckelberg fields, and on substituting back into the vierbein action we can show full equivalence with the tensorial form of dRGT massive gravity.[45]

## Cosmology

If the graviton mass ${\displaystyle m}$ is comparable to the Hubble rate ${\displaystyle H_{0}}$, then at cosmological distances the mass term can produce a repulsive gravitational effect that leads to cosmic acceleration. Because, roughly speaking, the enhanced diffeomorphism symmetry in the limit ${\displaystyle m=0}$ protects a small graviton mass from large quantum corrections, the choice ${\displaystyle m\sim H_{0}}$ is in fact technically natural.[46] Massive gravity thus may provide a solution to the cosmological constant problem: why do quantum corrections not cause the Universe to accelerate at extremely early times?

However, it turns out that flat and closed Friedmann–Lemaître–Robertson–Walker cosmological solutions do not exist in dRGT massive gravity with a flat reference metric.[13] Open solutions and solutions with general reference metrics suffer from instabilities.[47] Therefore, viable cosmologies can only be found in massive gravity if one abandons the cosmological principle that the Universe is uniform on large scales, or otherwise generalizes dRGT. For instance, cosmological solutions are better behaved in bigravity,[14] the theory which extends dRGT by giving ${\displaystyle f_{\mu \nu }}$ dynamics. While these tend to possess instabilities as well,[48][49] those instabilities might find a resolution in the nonlinear dynamics (through a Vainshtein-like mechanism) or by pushing the era of instability to the very early Universe.[15]

## 3D massive gravity

A special case exists in three dimensions, where a massless graviton does not propagate any degrees of freedom. Here several ghost-free theories of a massive graviton, propagating two degrees of freedom, can be defined. In the case of topologically massive gravity[1] one has the action

${\displaystyle S={\frac {M_{3}}{2}}\int d^{3}x{\sqrt {-g}}(R-2\Lambda )+{\frac {1}{4\mu }}\epsilon ^{\lambda \mu \nu }\Gamma _{\lambda \sigma }^{\rho }\left(\partial _{\mu }\Gamma _{\rho \nu }^{\sigma }+{\frac {2}{3}}\Gamma _{\mu \alpha }^{\sigma }\Gamma _{\nu \rho }^{\alpha }\right),}$

with ${\displaystyle M_{3}}$ the three-dimensional Planck mass. This is three-dimensional general relativity supplemented by a Chern-Simons-like term built out of the Christoffel symbols.

More recently, a theory referred to as new massive gravity has been developed,[2] which is described by the action

${\displaystyle S=M_{3}\int d^{3}x{\sqrt {-g}}\left[\pm R+{\frac {1}{m^{2}}}\left(R_{\mu \nu }R^{\mu \nu }-{\frac {3}{8}}R^{2}\right)\right].}$

## Relation to gravitational waves

The 2016 discovery of gravitational waves[50] and subsequent observations have yielded constraints on the maximum mass of gravitons, if they are massive at all. Following the GW170104 event, the graviton's Compton wavelength was found to be at least 1.6×1016 m, or about 1.6 light-years, corresponding to a graviton mass of no more than 7.7×10−23 eV/c2.[16] This relation between wavelength and energy is calculated with the same formula (the Planck–Einstein relation) that relates electromagnetic wavelength to photon energy. However, photons, which have only energy and no mass, are fundamentally different from massive gravitons in this respect, since the Compton wavelength of the graviton is not equal to the gravitational wavelength. Instead, the lower-bound graviton Compton wavelength is about 9×109 times greater than the gravitational wavelength for the GW170104 event, which was ~ 1,700 km. This is because the Compton wavelength is defined by the rest mass of the graviton and is an invariant scalar quantity.

Review articles
• de Rham, Claudia (2014), "Massive Gravity", Living Reviews in Relativity, 17 (1): 7, arXiv:1401.4173, Bibcode:2014LRR....17....7D, doi:10.12942/lrr-2014-7, PMC 5256007, PMID 28179850
• Hinterbichler, Kurt (2012), "Theoretical Aspects of Massive Gravity", Reviews of Modern Physics, 84 (2): 671–710, arXiv:1105.3735, Bibcode:2012RvMP...84..671H, doi:10.1103/RevModPhys.84.671, S2CID 119279950

## References

1. ^ a b Deser, Stanley; Jackiw, R.; Templeton, S. (1982). "Topologically Massive Gauge Theories". Annals of Physics. 140 (2): 372–411. Bibcode:1982AnPhy.140..372D. doi:10.1016/0003-4916(82)90164-6.
2. ^ a b Bergshoeff, Eric A.; Hohm, Olaf; Townsend, Paul K. (2009). "Massive Gravity in Three Dimensions". Phys. Rev. Lett. 102 (20): 201301. arXiv:0901.1766. Bibcode:2009PhRvL.102t1301B. doi:10.1103/PhysRevLett.102.201301. PMID 19519014. S2CID 7800235.
3. ^ Ogievetsky, V.I; Polubarinov, I.V (November 1965). "Interacting field of spin 2 and the einstein equations". Annals of Physics. 35 (2): 167–208. Bibcode:1965AnPhy..35..167O. doi:10.1016/0003-4916(65)90077-1.
4. ^ Mukohyama, Shinji; Volkov, Mikhail S. (2018-10-22). "The Ogievetsky-Polubarinov massive gravity and the benign Boulware–Deser mode". Journal of Cosmology and Astroparticle Physics. 2018 (10): 037. arXiv:1808.04292. Bibcode:2018JCAP...10..037M. doi:10.1088/1475-7516/2018/10/037. ISSN 1475-7516. S2CID 119329289.
5. ^ Ogievetsky, V. I; Polubarinov, I. V (1965-11-01). "Interacting field of spin 2 and the einstein equations". Annals of Physics. 35 (2): 167–208. Bibcode:1965AnPhy..35..167O. doi:10.1016/0003-4916(65)90077-1. ISSN 0003-4916.
6. ^ Curtright, T. L.; Alshal, H. (2019-10-01). "Massive dual spin 2 revisited". Nuclear Physics B. 948: 114777. arXiv:1907.11532. Bibcode:2019NuPhB.94814777C. doi:10.1016/j.nuclphysb.2019.114777. ISSN 0550-3213.
7. ^ Alshal, H.; Curtright, T. L. (2019-09-10). "Massive dual gravity in N spacetime dimensions". Journal of High Energy Physics. 2019 (9): 63. arXiv:1907.11537. Bibcode:2019JHEP...09..063A. doi:10.1007/JHEP09(2019)063. ISSN 1029-8479. S2CID 198953238.
8. ^ Nicolis, Alberto; Rattazzi, Riccardo; Trincherini, Enrico (2009-03-31). "Galileon as a local modification of gravity". Physical Review D. 79 (6): 064036. arXiv:0811.2197. Bibcode:2009PhRvD..79f4036N. doi:10.1103/PhysRevD.79.064036. S2CID 18168398.
9. ^ Deffayet, C.; Esposito-Farèse, G.; Vikman, A. (2009-04-03). "Covariant Galileon". Physical Review D. 79 (8): 084003. arXiv:0901.1314. Bibcode:2009PhRvD..79h4003D. doi:10.1103/PhysRevD.79.084003. S2CID 118855364.
10. ^ Curtright, Thomas L.; Fairlie, David B. (2012). "A Galileon Primer". arXiv:1212.6972 [hep-th].
11. ^ de Rham, Claudia; Keltner, Luke; Tolley, Andrew J. (2014-07-21). "Generalized Galileon duality". Physical Review D. 90 (2): 024050. arXiv:1403.3690. Bibcode:2014PhRvD..90b4050D. doi:10.1103/PhysRevD.90.024050. S2CID 118615285.
12. ^ a b Hassan, S.F.; Rosen, Rachel A. (2012). "Bimetric Gravity from Ghost-free Massive Gravity". JHEP. 1202 (2): 126. arXiv:1109.3515. Bibcode:2012JHEP...02..126H. doi:10.1007/JHEP02(2012)126. S2CID 118427524.
13. ^ a b D'Amico, G.; de Rham, C.; Dubovsky, S.; Gabadadze, G.; Pirtskhalava, D.; Tolley, A.J. (2011). "Massive Cosmologies". Phys. Rev. D84 (12): 124046. arXiv:1108.5231. Bibcode:2011PhRvD..84l4046D. doi:10.1103/PhysRevD.84.124046. S2CID 118571397.
14. ^ a b Akrami, Yashar; Koivisto, Tomi S.; Sandstad, Marit (2013). "Accelerated expansion from ghost-free bigravity: a statistical analysis with improved generality". JHEP. 1303 (3): 099. arXiv:1209.0457. Bibcode:2013JHEP...03..099A. doi:10.1007/JHEP03(2013)099. S2CID 54533200.
15. ^ a b Akrami, Yashar; Hassan, S.F.; Könnig, Frank; Schmidt-May, Angnis; Solomon, Adam R. (2015). "Bimetric gravity is cosmologically viable". Physics Letters B. 748: 37–44. arXiv:1503.07521. Bibcode:2015PhLB..748...37A. doi:10.1016/j.physletb.2015.06.062. S2CID 118371127.
16. ^ a b B. P. Abbott; et al. (LIGO Scientific Collaboration and Virgo Collaboration) (1 June 2017). "GW170104: Observation of a 50-Solar-Mass Binary Black Hole Coalescence at Redshift 0.2". Physical Review Letters. 118 (22): 221101. arXiv:1706.01812. Bibcode:2017PhRvL.118v1101A. doi:10.1103/PhysRevLett.118.221101. PMID 28621973. S2CID 206291714.
17. ^ L. Bernus; et al. (18 October 2019). "Constraining the Mass of the Graviton with the Planetary Ephemeris INPOP". Physical Review Letters. 123 (16): 161103. arXiv:1901.04307. Bibcode:2019PhRvL.123p1103B. doi:10.1103/PhysRevLett.123.161103. PMID 31702347. S2CID 119427663.
18. ^ Fierz, Markus; Pauli, Wolfgang (1939). "On relativistic wave equations for particles of arbitrary spin in an electromagnetic field". Proc. Roy. Soc. Lond. A. 173 (953): 211–232. Bibcode:1939RSPSA.173..211F. doi:10.1098/rspa.1939.0140.
19. ^ van Dam, Hendrik; Veltman, Martinus J. G. (1970). "Massive and massless Yang-Mills and gravitational fields". Nucl. Phys. B. 22 (2): 397–411. Bibcode:1970NuPhB..22..397V. doi:10.1016/0550-3213(70)90416-5. hdl:1874/4816.
20. ^ Zakharov, Valentin I. (1970). "Linearized gravitation theory and the graviton mass". JETP Lett. 12: 312. Bibcode:1970JETPL..12..312Z.
21. ^ Vainshtein, A.I. (1972). "To the problem of nonvanishing gravitation mass". Phys. Lett. B. 39 (3): 393–394. Bibcode:1972PhLB...39..393V. doi:10.1016/0370-2693(72)90147-5.
22. ^ Freund, Peter G. O.; Maheshwari, Amar; Schonberg, Edmond (August 1969). "Finite-Range Gravitation". The Astrophysical Journal. 157: 857. Bibcode:1969ApJ...157..857F. doi:10.1086/150118. ISSN 0004-637X.
23. ^ Boulware, David G.; Deser, Stanley (1972). "Can gravitation have a finite range?" (PDF). Phys. Rev. D. 6 (12): 3368–3382. Bibcode:1972PhRvD...6.3368B. doi:10.1103/PhysRevD.6.3368.
24. ^ Creminelli, Paolo; Nicolis, Alberto; Papucci, Michele; Trincherini, Enrico (2005). "Ghosts in massive gravity". JHEP. 0509 (9): 003. arXiv:hep-th/0505147. Bibcode:2005JHEP...09..003C. doi:10.1088/1126-6708/2005/09/003. S2CID 5702596.
25. ^ de Rham, Claudia; Gabadadze, Gregory (2010). "Generalization of the Fierz–Pauli Action". Phys. Rev. D. 82 (4): 044020. arXiv:1007.0443. Bibcode:2010PhRvD..82d4020D. doi:10.1103/PhysRevD.82.044020. S2CID 119289878.
26. ^ a b de Rham, Claudia; Gabadadze, Gregory; Tolley, Andrew J. (2011). "Resummation of Massive Gravity". Phys. Rev. Lett. 106 (23): 231101. arXiv:1011.1232. Bibcode:2011PhRvL.106w1101D. doi:10.1103/PhysRevLett.106.231101. PMID 21770493. S2CID 3564069.
27. ^ Hassan, S.F.; Rosen, Rachel A. (2012). "Resolving the Ghost Problem in non-Linear Massive Gravity". Phys. Rev. Lett. 108 (4): 041101. arXiv:1106.3344. Bibcode:2012PhRvL.108d1101H. doi:10.1103/PhysRevLett.108.041101. PMID 22400821. S2CID 17185069.
28. ^ Hassan, S.F.; Rosen, Rachel A. (2012). "Confirmation of the Secondary Constraint and Absence of Ghost in Massive Gravity and Bimetric Gravity". JHEP. 1204 (4): 123. arXiv:1111.2070. Bibcode:2012JHEP...04..123H. doi:10.1007/JHEP04(2012)123. S2CID 54517385.
29. ^ a b Hassan, S.F.; Rosen, Rachel A. (2011). "On Non-Linear Actions for Massive Gravity". JHEP. 1107 (7): 009. arXiv:1103.6055. Bibcode:2011JHEP...07..009H. doi:10.1007/JHEP07(2011)009. S2CID 119240485.
30. ^ Arkani-Hamed, Nima; Georgi, Howard; Schwartz, Matthew D. (June 2003). "Effective field theory for massive gravitons and gravity in theory space". Annals of Physics. 305 (2): 96–118. arXiv:hep-th/0210184. Bibcode:2003AnPhy.305...96A. doi:10.1016/S0003-4916(03)00068-X. S2CID 1367086.
31. ^ Hassan, S.F.; Rosen, Rachel A.; Schmidt-May, Angnis (2012). "Ghost-free Massive Gravity with a General Reference Metric". JHEP. 1202 (2): 026. arXiv:1109.3230. Bibcode:2012JHEP...02..026H. doi:10.1007/JHEP02(2012)026. S2CID 119254994.
32. ^ de Rham, Claudia; Gabadadze, Gregory; Tolley, Andrew J. (May 2012). "Ghost free massive gravity in the Stückelberg language". Physics Letters B. 711 (2): 190–195. arXiv:1107.3820. Bibcode:2012PhLB..711..190D. doi:10.1016/j.physletb.2012.03.081. S2CID 119088565.
33. ^ Alberte, Lasma; Khmelnitsky, Andrei (September 2013). "Reduced massive gravity with two Stückelberg fields". Physical Review D. 88 (6): 064053. arXiv:1303.4958. Bibcode:2013PhRvD..88f4053A. doi:10.1103/PhysRevD.88.064053. ISSN 1550-7998. S2CID 118668426.
34. ^ Hassan, S. F.; Rosen, Rachel A. (July 2011). "On non-linear actions for massive gravity". Journal of High Energy Physics. 2011 (7): 9. arXiv:1103.6055. Bibcode:2011JHEP...07..009H. doi:10.1007/JHEP07(2011)009. ISSN 1029-8479. S2CID 119240485.
35. ^ de Rham, Claudia; Gabadadze, Gregory; Tolley, Andrew J. (2011-06-10). "Resummation of Massive Gravity". Physical Review Letters. 106 (23): 231101. arXiv:1011.1232. Bibcode:2011PhRvL.106w1101D. doi:10.1103/PhysRevLett.106.231101. ISSN 0031-9007. PMID 21770493.
36. ^ Rham, Claudia de; Tolley, Andrew J (2010-05-14). "DBI and the Galileon reunited". Journal of Cosmology and Astroparticle Physics. 2010 (5): 015. arXiv:1003.5917. Bibcode:2010JCAP...05..015D. doi:10.1088/1475-7516/2010/05/015. ISSN 1475-7516. S2CID 118627727.
37. ^ a b de Rham, Claudia (December 2014). "Massive Gravity". Living Reviews in Relativity. 17 (1): 7. arXiv:1401.4173. Bibcode:2014LRR....17....7D. doi:10.12942/lrr-2014-7. ISSN 2367-3613. PMC 5256007. PMID 28179850.
38. ^ Maheshwari, A. (March 1972). "Spin-2 field theories and the tensor-field identity". Il Nuovo Cimento A. 8 (2): 319–330. Bibcode:1972NCimA...8..319M. doi:10.1007/BF02732654. ISSN 0369-3546. S2CID 123767732.
39. ^ Ogievetsky, V. I; Polubarinov, I. V (1965-11-01). "Interacting field of spin 2 and the einstein equations". Annals of Physics. 35 (2): 167–208. Bibcode:1965AnPhy..35..167O. doi:10.1016/0003-4916(65)90077-1. ISSN 0003-4916.
40. ^ Koyama, Kazuya; Niz, Gustavo; Tasinato, Gianmassimo (December 2011). "The self-accelerating universe with vectors in massive gravity". Journal of High Energy Physics. 2011 (12): 65. arXiv:1110.2618. Bibcode:2011JHEP...12..065K. doi:10.1007/JHEP12(2011)065. ISSN 1029-8479. S2CID 118329368.
41. ^ Pitts, J. Brian (August 2019). "Cosmological Constant $\Lambda$ vs. Massive Gravitons: A Case Study in General Relativity Exceptionalism vs. Particle Physics Egalitarianism". arXiv:1906.02115 [physics.hist-ph].
42. ^ Arkani-Hamed, Nima; Cohen, Andrew G.; Georgi, Howard (May 2001). "(De)Constructing Dimensions". Physical Review Letters. 86 (21): 4757–4761. arXiv:hep-th/0104005. Bibcode:2001PhRvL..86.4757A. doi:10.1103/PhysRevLett.86.4757. ISSN 0031-9007. PMID 11384341. S2CID 4540121.
43. ^ Hinterbichler, Kurt; Rosen, Rachel A. (2012). "Interacting Spin-2 Fields". JHEP. 1207 (7): 047. arXiv:1203.5783. Bibcode:2012JHEP...07..047H. doi:10.1007/JHEP07(2012)047. S2CID 119255545.
44. ^ a b Ondo, Nicholas A.; Tolley, Andrew J. (November 2013). "Complete decoupling limit of ghost-free massive gravity". Journal of High Energy Physics. 2013 (11): 59. arXiv:1307.4769. Bibcode:2013JHEP...11..059O. doi:10.1007/JHEP11(2013)059. ISSN 1029-8479. S2CID 119101943.
45. ^ Groot Nibbelink, S.; Peloso, M.; Sexton, M. (August 2007). "Nonlinear properties of vielbein massive gravity". The European Physical Journal C. 51 (3): 741–752. arXiv:hep-th/0610169. Bibcode:2007EPJC...51..741G. doi:10.1140/epjc/s10052-007-0311-x. ISSN 1434-6044. S2CID 14575306.
46. ^ de Rham, Claudia; Heisenberg, Lavinia; Ribeiro, Raquel H. (2013). "Quantum Corrections in Massive Gravity". Phys. Rev. D. 88 (8): 084058. arXiv:1307.7169. Bibcode:2013PhRvD..88h4058D. doi:10.1103/PhysRevD.88.084058. S2CID 118328264.
47. ^ de Felice, Antonio; Gümrükçüoğlu, A. Emir; Lin, Chunshan; Mukohyama, Shinji (2013). "On the cosmology of massive gravity". Class. Quantum Grav. 30 (18): 184004. arXiv:1304.0484. Bibcode:2013CQGra..30r4004D. doi:10.1088/0264-9381/30/18/184004. S2CID 118669165.
48. ^ Comelli, Denis; Crisostomi, Marco; Pilo, Luigi (2012). "Perturbations in Massive Gravity Cosmology". JHEP. 1206 (6): 085. arXiv:1202.1986. Bibcode:2012JHEP...06..085C. doi:10.1007/JHEP06(2012)085. S2CID 119205963.
49. ^ Könnig, Frank; Akrami, Yashar; Amendola, Luca; Motta, Mariele; Solomon, Adam R. (2014). "Stable and unstable cosmological models in bimetric massive gravity". Phys. Rev. D. 90 (12): 124014. arXiv:1407.4331. Bibcode:2014PhRvD..90l4014K. doi:10.1103/PhysRevD.90.124014. S2CID 86860987.
50. ^ B. P. Abbott; et al. (LIGO Scientific Collaboration and Virgo Collaboration) (2016). "Observation of Gravitational Waves from a Binary Black Hole Merger". Phys. Rev. Lett. 116 (6): 061102. arXiv:1602.03837. Bibcode:2016PhRvL.116f1102A. doi:10.1103/PhysRevLett.116.061102. PMID 26918975. S2CID 124959784.