Gravity with Explicit Diffeomorphism Breaking

Modified theories of gravity that explicitly break diffeomorphism invariance have been used for over a decade to explore open issues related to quantum gravity, dark energy, and dark matter. At the same time, the Standard-Model Extension (SME) has been widely used as a phenomenological framework in investigations of spacetime symmetry breaking. Until recently, it was thought that the SME was suitable only for theories with spontaneous spacetime symmetry breaking due to consistency conditions stemming from the Bianchi identities. However, it has recently been shown that, particularly with matter couplings included, the consistency conditions can also be satisfied in theories with explicit breaking. An overview of how this is achieved is presented, and two examples are examined. The first is massive gravity, which includes a nondynamical background tensor. The second is a model based on a low-energy limit of Ho\v rava gravity, where spacetime has a physically preferred foliation. In both cases, bounds on matter--gravity interactions that explicitly break diffeomorphisms are obtained using the SME.


I. INTRODUCTION
Diffeomorphism invariance is a fundamental symmetry in Einstein's General Relativity (GR). As a result of this symmetry and the Bianchi identities, when matter and gravity interact in GR, the energy-momentum of the matter fields is guaranteed to be covariantly conserved when the matter fields obey their equations of motion. Mathematically, this linkage between matter dynamics, geometrical constraints, and energy-momentum conservation is a result of Noether identities derived from the local diffeomorphism invariance of the Einstein-Hilbert and matter action [1,2].
As elegant and experimentally accurate as GR is as a classical theory, it is widely believed that it must ultimately merge with quantum physics at a very high-energy scale, such as the Planck scale. It is hoped that in the context of a new quantum theory of gravity, which would supersede GR, unsolved problems having to do with quantization might find solutions and unknown entities such as dark energy and dark matter might find physical explanations. Amongst the many avenues of research looking for clues of a new quantum theory of gravity, a number of them involve breaking local diffeomorphism invariance. Examples include mechanisms in string field theory [3][4][5], massive gravity [6,7], Chern-Simons gravity [8], Bumblebee [9][10][11][12][13][14][15][16][17] and Einstein-Aether models [18], and Hořava gravity [19][20][21][22][23]. The Standard-Model Extension (SME) provides a phenomenological theoretical framework used in experimental tests of gravity and matter interactions that break spacetime symmetries [24][25][26][27][28][29][30][31], and in recent decades, many high-precision bounds on possible symmetry-breaking terms have been obtained [32].
Diffeomorphism breaking can occur either spontaneously or explicitly. In the case of spontaneous breaking, a dynamical tensor field acquires a nonzero vacuum expectation value, which breaks the symmetry, or more precisely hides the symmetry. In this case, since all of the fields are dynamical, and excitations in the form of Nambu-Goldstone (NG) or Higgs-like excitations occur, the usual linkage between matter dynamics, the Bianchi identities, and covariant energy-momentum conservation still holds. On the other hand, with explicit breaking, nondynamical backgrounds can appear or the action can contain terms that are not scalars under the full diffeomorphism group. With explicit breaking, the usual relations between dynamics, geometric identities, and covariant energy-momentum conservation get disrupted, and viable theories must ensure that there are no inconsistencies that would rule them out. To avoid this possibility, the gravity sector of the SME was originally constructed assuming that the symmetry breaking was spontaneous [27].
However, more recently, the question of whether inconsistencies must occur or can be evaded when the symmetry breaking is explicit has been examined more closely [33][34][35][36][37][38][39][40][41]. It turns out that in many cases, the SME can still be applied to gravity theories with explicit diffeomorphism breaking, in particular when matter-gravity couplings are included. At the same time, in more restrictive limits, such as pure gravity in a first-order perturbative treatment, evading inconsistency puts very restrictive constraints on the geometry or the background fields, which makes using the SME not so useful in these limits. One goal of this paper is to give an overview of how this plays out. First, some general features are described. Then, specific examples of gravity models with explicit diffeomorphism breaking are examined, including massive gravity and a model based on a low-energy or infrared (IR) limit of Hořava gravity. Matter-gravity couplings involving scalar and tensor fields as well as fermions are investigated, and by matching them to terms in the SME, experimental bounds are obtained.
It should be pointed out that recently, the gravity sector of the SME has been generalized to include a much wider class of background fields at the level of effective field theory, including backgrounds that explicitly break spacetime symmetry [41]. Many new phenomenological tests are likely to be performed using this broader set of interaction terms. However, in this paper, the question under consideration is whether previous tests of spacetime symmetry breaking that placed bounds on terms in the original formulation of the SME (based on spontaneous breaking) can be applied to the case of explicit breaking. The answer is that, in some cases, they can be when matter-gravity couplings are included, as can be shown using a Stückelberg approach.
The outline of this paper is as follows. To begin, some background on explicit diffeomorphism breaking is given in Section II. This is followed in Section III by an examination of the consistency conditions that must hold with explicit breaking. Section IV looks at the question of whether the SME based on spontaneous breaking can be used to investigate gravity theories with explicit breaking. Examples of how it can be used when matter-gravity couplings are included are given in Section V. A summary and concluding remarks are then presented in Section VI.

II. EXPLICIT DIFFEOMORPHISM BREAKING
Infinitesimal diffeorphisms in GR are mappings of a spacetime manifold back to itself, where points on the manifold are locally translated along a small vector ξ µ , such that x µ → x µ + ξ µ . Changes in vector and tensor fields at a specified point x µ are given by their Lie derivatives, and collectively the transformations act like a gauge symmetry in GR. For example, the metric tensor changes as where L ξ g µν = D µ ξ ν + D ν ξ µ is the Lie derivative defined using covariant derivatives D µ . The convention used in this paper is that the metric has signs (−, +, +, +) on its diagonal, and it is assumed that there is no torsion. Hence, the covariant derivative of a spacetime tensor is defined with a connection given by the Christoffel symbol. A generic dynamical matter field, denoted here as f µ , transforms as where for a contravariant vector, An action in GR with a Lagrangian L(g µν , f ν ) depending on the metric and matter fields is then given generically as where g is the determinant of the metric. When L is a scalar, such an action in GR is invariant under local diffeomorphism gauge transformations, obeying δS = 0. In GR, infinitesimal general coordinate transformations, x µ → x µ ′ (x) = x µ − ξ µ , are often called diffeomorphisms as well, though technically they are the passive versions of the gauge-like transformations (with the direction of ξ µ reversed). Nonetheless, changes of vector and tensor fields under these transformations at a given coordinate x µ are also given by their Lie derivatives. Since a covariant theory will have an action that is unchanged under general coordinate transformations, δS = 0 holds as well under these transformations.
At the level of effective field theory, diffeomorphisms are broken when a fixed background tensor, denoted generically here ask χ , appears in the Lagrangian. The backgroundk χ is a pre-specified field that in general has spacetime dependence, but it does not respond dynamically when it interacts with other fields in the theory. As a fixed background,k χ does not transform under infinitesimal diffeomorphisms. Instead, it obeysk χ →k χ . Hence, an action containing the metric, dynamical matter fields, and a fixed background may not be invariant, depending on how the symmetry is broken. If the symmetry breaking is spontaneous, then the backgroundk χ is a vacuum expectation value of a fully dynamical field k χ . In this case,k χ is accompanied in the action by NG and Higgs-like excitations, and as a result the action remains invariant under diffeomorphisms. With spontaneous diffeomorphism breaking,k χ is still a dynamical field, since it is a solution of the vacuum equations of motion. In contrast, explicit diffeomorphism breaking occurs when the fixed background is nondynamical. In this case,k χ is not a solution of the equations of motion. Instead, it is a pre-determined background that is inserted into the action, sometimes referred to as an object with prior geometry. With such a background tensor included, the action takes the form and δS = 0 under diffeomorphisms. To include fermion fields in a gravity theory, a vierbein formalism is used. Effectively, the vierbein, e a µ , gives tensor fields their components, at each spacetime point, with respect to a local Lorentz basis, where Latin letters are used to denote the local indices. For the metric, can be viewed as a defining relation for the vierbein, where in this case its local components are those of the Minkowski metric η ab . A generic matter field f µ can be rewritten as f µ = e µ a f a , where f a are its components with respect to the local Lorentz frame, and e µ a is the inverse vierbein obeying e µ a e b µ = δ a b and e µ a e a ν = δ µ ν . The Lagrangian for a fermion field ψ, in the absence of diffeomorphism breaking, can then be written as eL = eie µ aψ γ a D µ ψ, where e is the determinant of the vierbein and γ a is a Dirac matrix. The covariant derivative D µ acting on a fermion field is D µ ψ = ∂ µ ψ + 1 4 iω ab µ σ ab ψ, where σ ab = γ a γ b − γ b γ a and ω ab µ is the spin connection. In addition to diffeomorphism invariance, local Lorentz invariance becomes a second fundamental spacetime symmetry in gravity theories when a vierbein formalism is introduced, where the fermion fields and local tensor fields transform as representations of the local Lorentz group.
When a fixed nondynamical background is included in a gravity theory in a vierbein description, its components in both the spacetime frame and with respect to the local Lorentz basis remain fixed under diffeomorphisms and local Lorentz transformations. For example,k χ →k χ andk a →k a under both types of spacetime transformations. With both of these components being nondynamical and remaining fixed, they cannot be connected using the physical vierbein, e.g.,k χ = e a χk a . Instead, they are related by a nondynamical vierbein denoted here asē a χ , which does not transform under either diffeomorphisms or local Lorentz transformations. Depending on whetherk χ ,k a orē a χ appears in a gravitational action, either or both of local diffeomorphism invariance and local Lorentz invariance can be explicitly broken.

III. CONSISTENCY CONDITIONS WITH EXPLICIT BREAKING
Consider a gravity theory that contains a nondynamical backgroundk χ that explicitly breaks diffeomorphism invariance, where χ generically denotes an arbitrary number of spacetime indices that can be covariant, contravariant, or mixed. The action is where an Einstein-Hilbert term is used for the pure-gravity sector.
Sincek χ is not dynamical, it does not have equations of motion, and its Euler-Lagrange expression does not need to vanish, In contrast, the metric g µν and matter fields f ν are dynamical, and their equations of motion hold on shell. For the metric, its equations of motion yield the Einstein equations, where G µν is the Einstein tensor, and T µν is the energy-momentum tensor that depends on the matter fields and the backgroundk χ . Since the contracted Bianchi identities require that D µ G µν = 0 holds off shell, consistency with the Einstein equations requires that D µ T µν = 0 must hold when the metric and matter fields are on shell. However, the fact that the backgroundk χ appears in T µν and does not itself obey equations of motion, calls into question whether consistency with the condition of energy-momentum conservation, D µ T µν = 0, can hold. General conditions that must hold for consistency in a theory with explicit diffeomorphism breaking can be obtained as mathematical Noether identities associated with general coordinate invariance. Any physical theory should be independent of the choice of coordinates, which requires that the action must be unchanged under a general coordinate transformation. Although the backgrounds remain fixed under diffeomorphisms, they do transform as tensors under observer-general coordinate transformations. As noted above, when an infinitesimal general coordinate transformation x µ → x µ ′ (x) = x µ − ξ µ is performed, vector and tensor fields transform with changes given by their Lie derivatives. This includes any fixed background tensors that are included in the action. Thus,k χ →k χ + L ξkχ under infinitesimal general coordinate transformations. With the action remaining unchanged under these transformations, obeying δS = 0, an off-shell identity can be found, which relates the Euler-Lagrange expressions for the metric, matter field, and the background. Putting the matter field on shell imposes that δS δf ν = 0.
Combining this with the contracted Bianchi identity results in four equations that must hold as a result of general coordinate invariance: Here, for simplicity, the identity is written out for the case of a background tensork χ that is a covariant vector. For other types of tensors, the Lie derivatives will differ and a different number of terms with different indices can appear, with the corresponding Euler-Lagrange expression appearing in each term. However, for the purpose of pointing out the important features, a vector background suffices.
Consistency with covariant energy-momentum conservation evidently requires that the right-hand side of (10) must vanish. However, the fact that the background does not satisfy Euler-Lagrange equations, as indicated in Equation (7), makes this appear problematic.
Nonetheless, the right-hand side of (10) can vanish quite generally. This is because, with explicit diffeomorphism breaking, four local symmetries are lost, which results in the theory having four additional degrees of freedom compared to a theory with unbroken symmetry. These extra degrees of freedom become additional metric components, which can take values that make D µ T µν = 0 hold in (10). Note that the metric will in general appear on the right-hand side of (10), because the covariant derivatives depend on it. However, if the theory is limited in some way that does not permit a sufficient number of degrees of freedom to appear, then imposing D µ T µν = 0 will put constraints on the background or on the geometry, resulting in a theory that is either inconsistent or uninteresting. Examples of when this occurs include theories where an ansatz form of the metric is used that does not have a sufficient number of degrees of freedom to satisfy the four equations in (10). Alternatively, if a linear perturbative treatment is used that suppresses the appearance of the metric in nonlinear terms in (10), the result again can be that the theory is inconsistent in such a limit.
A Stückelberg approach is often used in theories with explicit diffeomorphism breaking to reveal the extra degrees of freedom that appear. This approach consists of introducing four dynamical scalar fields φ A with A = 0, 1, 2, 3 and rewriting the background tensor in terms of these [42]. For example, the backgroundk µν used in massive gravity is replaced byk wherever it appears in the action. This introduces four extra degrees of freedom in the form of the scalars φ A , but it also restores the four diffeomorphism symmetries. The net result is that there is no change in the number of degrees of freedom compared to the original theory with explicit breaking. For infinitesimal excitations, the scalar can be expanded as , where x µ are the coordinates and π µ are the Stückelberg excitations. As can be seen when π µ → 0, the Stückelberg expression fork µν in (11) reduces to the original nondynamical tensor, which shows that the explicit breaking theory is equivalent to a gauge-fixed form of the Stückelberg modification. However, if the excitations are nonzero, the expression for the background to first order in the excitations becomes This expression can be recognized as the Lie derivative acting onk µν , with π µ playing the role of the transformation vector. In a theory with spontaneous diffeomorphism breaking, the excitations π µ generate the broken symmetry transformations and are therefore the NG excitations. What the Stückelberg approach does in a theory with explicit breaking is that it adds in the same NG modes as in the corresponding theory with spontaneous breaking. It also restores diffeomorphism invariance. However, note that there are no Higgs-like excitations generated in the Stückelberg approach, which distinguishes it from the case of spontaneous breaking. It is also the case that with explicit breaking, k µν is still not a dynamical solution, whereas with spontaneous breaking,k µν is a vacuum solution for a dynamical tensor field.
As a further example, the Stückelberg approach can be used in the action in (6) involving a vector backgroundk χ . In this case, the replacement isk χ → (∂ χ φ A )k A (φ). Since the scalars φ A are dynamical in the Stückelberg approach, the action in (6) can be varied with respect to these fields to obtain their equations of motion. If the diffeomorphism gauge freedom is then fixed by setting φ A = δ A µ x µ , which sets the NG modes to zero, the resulting equation of motion is Notice that this is the same expression as on the right-hand side of the identity in (10). Thus, when the equations of motion for the Stückelberg scalars are satisfied, the identity in (10) sets D µ T µν = 0. However, there is no guarantee that the Stückelberg equations of motion in (13) must hold. Since the Stückelberg scalars themselves have been gauge fixed so that they equal the coordinates, it is the metric field, including four extra components that have been left unfixed by the gauge choice, that must take values that satisfy (13). In general, this can occur, but in models with an ansatz choice of metric or where a perturbative expansion is used, it may not be possible for the metric to satisfy these equations. Therefore, the consistency issues can still persist when a Stückelberg approach is used.
The same links as in GR between the Bianchi identities, dynamcs, and energy-momentum conservation hold in the SME when the background coefficients are assumed to be dynamical and arising as a result of spontaneous spacetime symmetry breaking. For this reason, the SME originally assumed that the symmetry breaking was spontaneous and that the theory also contains NG and Higgs-like excitations in addition to the background fields. These small fluctuations play an important role in developing a useful post-Newtonian limit of the SME and in incorporating matter-gravity couplings in a systematic way.
Field redefinitions also must be taken into account in gravity tests of spacetime symmetry breaking, because the SME coefficients are not all physical or independent [27,44,45,71,72]. For example, changes of coordinates or field redefinitions can be used to move some of the sensitivity to the symmetry breaking from one particle sector to another. With gravity, ten components of the SME coefficients are unphysical, which corresponds to having four coordinates and six local Lorentz bases to choose. As a result of these field redefinitions, experiments with matter and gravity that are designed to test spacetime symmetry breaking need to have sensitivity to more than one particle sector.
For the purposes of this paper, it is sufficient to consider a fermion of mass m and a photon field in gravity, where the SME coefficients are limited to those in the minimal sector with s µν , c µν , and (k F ) κλµν . With the Lorentz-invariant terms included, the action is The coefficients s µν , c µν , and (k F ) κλµν consist of background values and excitations, where, in this notation, bars indicate the background values and tildes indicate the excitations. With spontaneous symmetry breaking, the backgrounds are vacuum expectation values of dynamical tensors, and the excitations include NG and Higgs-like modes. In perturbative treatments, the metric can also be expanded about a Minkowski background, g µν = η µν + h µν , which permits the construction of a post-Newtonian framework.
To consider the question of whether these terms in the SME can be used to investigate gravity models with explicit diffeomorphism breaking, a number of different features must be examined, including how to maintain consistency between the Bianchi identities, equations of motion, and covariant energy-momentum conservation. When the symmetry breaking is explicit, the backgrounds are fixed and have no excitations, and hence the SME coefficients in (14) reduce to The backgroundss µν ,c µν , and (k F ) κλµν are not solutions for dynamical tensors in the theory. Instead, they are pre-determined fixed tensors put directly into the action. While the backgroundss µν ,c µν , and (k F ) κλµν have spacetime dependence, it is common in a first-order perturbative treatment about a flat Minkowski background to approximate them as constant.
Although the backgroundss µν ,c µν , and (k F ) κλµν have no excitations, using a Stückelberg approach effectively introduces excitations. Moreover, the excitations take the same form as the NG modes in the corresponding and original version of the SME based on spontaneous breaking. These can be denoted ass µν ,c µν , and (k F ) κλµν , but unlike the SME based on spontaneous breaking, they do not also include Higgs-like excitations. However, insofar as maintaining the consistency of the model is concerned, it suffices that the same NG modes appear as in the original SME with spontaneous breaking. With the Stückelberg NG modes included, D µ T µν = 0 can hold while maintaining consistency with the Bianchi identities. Nonetheless, the same caveats apply as described previously. If enough degrees of freedom in the metric are suppressed, either due to using a perturbative treatment or an ansatz solution, then the geometry of the theory is severely restricted or the theory can be outright inconsistent.
An example of when the SME with explicit diffeomorphism breaking does not provide a useful framework is provided by the pure-gravity sector in a post-Newtonian limit, wheres µν is the relevant SME coefficient. When the Stückelberg approach is applied, this introduces the NG modes π µ . By gauge-fixing the Stückelberg scalars, the NG modes can be moved into the metric. The additional excitations in h µν in a linearized approximation take the form ∂ µ π ν + ∂ ν π µ . However, when these additional excitations are substituted into the linearized Ricci tensor R µν in (14), the NG modes drop out. This is because the linearized curvature is gauge-invariant, and the excitations ∂ µ π ν + ∂ ν π µ have a puregauge form. Without the NG modes, it is not possible to satisfy the consistency conditions without putting severe restrictions on the geometry. For example, withs µν = 0, consistency of the theory requires that ∂ λ R µν = 0 must hold in a post-Newtonian limit, and such a condition placed on the curvature does not give a workable model.
However, if matter-gravity couplings are added, the situation changes. In this case, all three of the coefficientss µν , c µν , and (k F ) κλµν are included in the action in (14). Moreover, the NG modes do not drop out in this case, since they couple to more than just the linearized curvature. With matter-gravity couplings included, the SME can be used when there is explicit diffeomorphism breaking. This follows because the same methodology used in the case of spontaneous breaking is applicable with explicit breaking using a Stückelberg approach. The Stückelberg method introduces the NG modes, but not the Higgs-like modes. However, in the investigation in [45], it was assumed that the Higgs-like modes were either frozen out or have negligible excitations. Thus, the same perturbative treatment can be applied to the SME with explicit breaking as was worked out for the case of spontaneous breaking. In this method, the symmetry properties allow the NG modes to be eliminated in terms of the gravitational excitations and background SME coefficients. The dominant signals of spacetime symmetry breaking can then be extracted regardless of the details of the underlying theory.
With all of the SME coefficientss µν ,c µν , and (k F ) κλµν included in the action, ten of them can be removed using field redefinitions. Because of this, it is important that experiments testing for spacetime symmetry breaking in gravity have sensitivity to at least two of these sectors. Otherwise, a test in only one sector could involve a set of coefficients that are unphysical. With sensitivity to more than one set of SME coefficients, experiments can obtain bounds assuming that one of the sets has vanishing coefficients while making measurements on an independent set. Alternatively, bounds can be placed on combinations of coefficients from two different particle sectors. An example of the latter type is atom interferometry tests, which are sensitive to both the gravity and photons. The relevant coefficients in this case ares µν and the symmetric trace (k F ) αµ ν α . For comparisons across different experiments, a Sun-centered celestial equatorial frame is used. The components of the SME coefficients are obtained with respect to this frame, where indices are labeled using letters JK · · · for the spatial directions. The combination of gravity and electromagnetic parameters that are bounded in these experiments are denoted as σ JK , and measured bounds of order 10 −9 have been obtained for them.
Ultimately, the full consistency of a gravity theory with explicit diffeomorphism breaking requires that the additional metric modes satisfy the identities stemming from general coordinate invariance for both the matter and gravity sectors at the same time. For the gravity sector, this might involve nonlinear terms in a perturbative approach. However, to first order the results shown here for the SME with explicit breaking indicate that the leading-order sensitivity to spacetime symmetry breaking will be in the matter-gravity sector, since the pure-gravity sector does not have a useful first-order post-Newtonian limit.

V. EXAMPLES WITH MATTER COUPLINGS
Two examples of gravity theories with explicit diffeomorphism breaking are examined in this section. They are massive gravity and a model based on the IR limit of Hořava gravity. In massive gravity, a fixed background tensor explicitly breaks diffeomorphism invariance, while in Hořava gravity, diffeomorphism invariance is broken by intro-ducing a physical distinction between time and space in the form of a preferred foliation. In both of these models, subject to some assumptions that must hold, matter-gravity interactions can appear that can be matched to SME terms in (14), and experimental bounds can be placed on the relevant couplings.
It is important to keep in mind that both massive gravity and Hořava gravity have experimental constraints in the pure-gravity sector and that other consistency conditions involving stability, strong coupling in the gravity sector, black holes, and cosmology must be addressed as well. In fact, different formulations or modifications of the original versions of these theories have been proposed to address these other types of issues. (See, for example, the reviews in [6,7,[20][21][22][23].) However, the purpose of this paper is to examine how the SME can be applied, and in this context, the matter-gravity couplings that are introduced in the following subsections are independent of the original couplings in the pure-gravity sector, and therefore any bounds that are obtained for these new couplings will be independent of previously obtained bounds in the pure-gravity sector.

A. Massive Gravity
A ghost-free massive gravity theory has been found by de Rham, Gabadadze, and Tolley (dRGT) [6,7]. To form mass terms for the metric, a symmetric background tensor, denoted here ask µν , is coupled to the metric in a nonlinear potential, where the form of this potential eliminates the ghost mode. As a result of introducingk µν , the dRGT Lagrangian explicitly breaks diffeomorphism invariance. The dRGT action can be divided into a gravitational sector and a matter sector, and either a metric or vierbein formalism can be used. In a metric description, the action in the gravity sector has the form where µ is the graviton mass and U(X) is the potential that provides a mass term for the graviton. The potential is formed in terms of square roots X µ ν , which are defined as and which obey X µ α X α ν = g µαk αν . In a vierbein description, a background vierbein,v a µ , is introduced, which is defined byk µν =v a µv b ν η ab The potential U can then be written as a sum of elementary symmetric polynomials formed from products and traces of matrices If the vierbein and background vierbein obey a symmetry condition, then the pure-gravity metric and vierbein descriptions are equivalent, and X µ ν = γ µ ν [73,74]. Couplings of matter to both the metric g µν and the backgroundk µν become likely when quantum corrections are included. To avoid reintroducing the ghost, the matter fields couple with an effective metric g (eff) µν that is formed out of both the metric and the background field [75][76][77], while the gravity sector remains unchanged. The form of the effective metric is where γ µν = g µσ γ σ ν and α and β are constants. Note that with lowered indices, the square root matrix is symmetric, obeying γ µν = γ νµ . With β = 0, matter couplings to g (eff) µν break local Lorentz symmetry and diffeomorphisms. However, since Lorentz breaking is known to be small, to good approximation α ≃ 1 while β ≪ 1, and the effective metric therefore has the form It follows from this that the effective vierbein is µ a ≃ e a µ + βe a α γ α µ .
The matter terms in the action have the form (29) or if fermions are included, it is given as µ a , f ν , ψ).
Here, f ν represents a tensor matter field, while ψ denotes a fermion field. These terms can be matched to the SME using field redefinitions that change the effective metric back to the physical metric. These fields' redefinitions also result in the introduction of SME coefficients.
As an example, consider interactions in massive gravity involving photons with a fermion field. The action in this case is given as where m is the fermion mass. Here, the gravity sector depends on the physical metric g µν or vierbein e a µ , while the matter terms and D (eff) µ depend on the effective metric or vierbein. Field redefinitions can be used on g (eff)µν and e (eff)µ a , which at leading order reduce them to the physical metric and vierbein, while introducing SME coefficients defined in terms of the matrices γ µν . In this way, the massive gravity action S dRGT can be expressed entirely in terms of the physical metric or vierbein, but with additional interactions involving SME coefficients. The result is where the relevant SME coefficients that arise are given as Notice that (k F ) κλµν has been partially traced and that the fermion field ψ has been rescaled and relabeled as χ so that the action remains in a standard Dirac form. Notice also that the Lorentz-violating couplings to the effective metric g (eff) µν have been replaced by Lorentz-violating interactions with SME coefficients (k F ) α µαν and c µν . Because the field redefinitions were only on the matter sector, the coefficients u or s µν in (32) do not appear. Thus, with β = 0, it is not possible to make a further field redefinition to remove (k F ) α µαν and c µν without generating u and s µν . The potential term U is important in cosmology and gravitational radiation; however, it has a negligible effect on matter-gravity tests, because the graviton mass is so small with bounds on the order of µ ∼ < 10 −29 eV [78]. Thus, the effects of µ can be ignored at leading order in matter-gravity tests. Furthermore, to leading order, g µν andf µν can be approximated as Minkowski backgrounds, in which case the contributions of γ µν are also of order one. With these approximations, the spacetime symmetry breaking is determined by the parameter β in (27) and (28), and therefore β serves as a measure of the symmetry-breaking matter couplings in massive gravity.
To obtain physical bounds, experiments that are sensitive to two sets of SME coefficients are needed. In particular, matter-interferometry experiments, which have sensitivity in both the gravity and electromagnetic sectors of the SME [48,49], are suitable for placing bounds at the level of 10 −9 on combinations of s µν and the symmetric trace (k F ) αµ ν α . As a result, a bound of β ∼ < 10 −9 is obtained on the matter-gravity couplings in dRGT massive gravity.

B. Hořava-Based Model with a Preferred Spacetime Foliation
Hořava gravity explicitly breaks diffeomorphism invariance by introducing a physical preferred foliation of spacetime [19]. Using coordinates x µ = (t, x i ), with i = 1, 2, 3, where constant values of t distinguish the preferred foliations and x i labels spatial points, Hořava gravity introduces an anisotropic scaling in t and x i . This permits terms to be added to the action that maintain two time derivatives, but which give a theory with power-counting renormalizability. The preferred foliation explicitly breaks the full diffeomorphism group to a subgroup consisting of time reparametrizations and three-dimensional spatial diffeomorphisms, These are called foliation-preserving diffeomorphisms (FPdiffs). Notice that under FPdiffs, with x µ → x µ + ξ µ , the vector ξ µ obeys ∂ i ξ 0 = 0.
The pure-gravity sector of Hořava gravity has been explored extensively, and a number of issues have been found that go beyond maintaining consistency with the Bianchi identities, and as a result, a number of modified versions of Hořava gravity have been proposed [20][21][22][23]. In particular, the role of the extra-metric modes in Hořava gravity have been shown to play an important role in the pure-gravity sector [79][80][81]. However, for the purposes of this paper, it is the matter-gravity couplings in a low-energy limit that can most readily be matched to the SME, making the precise details of how the pure-gravity sector merges with GR less critical to consider.
If the low-energy or IR limit of Hořava gravity is to merge with GR and the SM, there should be a covariant formulation of Hořava gravity in this limit. There are different approaches that can be followed in constructing a covariant theory in the IR limit. The approach examined here (see [82][83][84] as an example of an alternative) involves projecting covariant vectors and tensors in four dimensions into the preferred foliation. Such projections use a timelike normal vector n µ , obeying n µ n µ = −1, which is orthogonal to the spatial foliations at each point. These vectors are used in GR in a Hamiltonian formulation, where the metric g µν is replaced by ADM variables (N, N i , g ij ), where N is the lapse, N i is the shift, and g ij is the three-dimensional spatial metric. In GR, the normal vectors n µ can be defined with respect to any foliation. However, in an IR limit of Hořava gravity, the vectors n µ become fixed with respect to the preferred foliation. As a result, n µ can transform under FPdiffs, which maintain it as a normal vector, but not under diffeomorphisms with ∂ i ξ 0 = 0, since these would point it in a nonorthogonal direction. A covariant theory constructed using n µ therefore explicitly breaks diffeomorphism invariance to FPdiffs, and it can be used to model Hořava gravity in the IR limit. If a mapping of such a model to terms in the SME can be found, then existing experimental bounds on the SME terms can be used to put limits on the corresponding couplings in the Hořava-based model.
The action in Hořava gravity in the IR limit can be divided into two sectors, consisting of a gravity sector S grav and a matter sector S matter , where in principle both sectors can break diffeomorphism invariance to FPdiffs, though with different coupling parameters. The pure-gravity sector S grav has been well studied, including how it can be matched up with the gravity sector of the SME [40]. It is defined using ADM variables for the metric, the extrinsic curvature K ij , its trace K = g ij K ij , and the curvature tensor in the three-dimensional spatial slice. Its kinetic terms involve a combination (K ij K ij − λ g K 2 ), where λ g is a running coupling constant. Here, each of the products K ij K ij and K 2 is separately a scalar under FPdiffs. However, unless λ g = 1, their difference breaks diffeomorphism invariance in four dimensions. There are also potential terms in S grav consisting of contractions of higher-dimensional spatial components that make Hořava gravity power-counting renormalizable in the high-energy limit. However, in the IR limit, the leading-order potential terms reduce to the three-dimensional curvature scalar R (3) and a cosmological constant term. For Hořava gravity to match with GR in the IR limit, the running coupling λ g must approach 1, in which case the kinetic and potential terms at leading order reproduce the usual Einstein-Hilbert and Λ terms in GR in ADM formalism. However, with small residual values of (1 − λ g ) in the IR limit, there can still be an additional symmetrybreaking term. With these assumptions, an action modeling the gravity sector of Hořava gravity in the IR limit can be written as where g (3) is the determinant of the three-dimensional metric. In this way, (1 − λ g ) ≪ 1 becomes a parameter that measures the spacetime symmetry breaking in the gravity sector, since it is the second term in S grav that gives rise to the spacetime symmetry breaking. It can be put in covariant form using the timelike unit vector n µ for the preferred foliation to redefine the four-dimensional metric g µν as where g ij is the three-dimensional inverse spatial metric. The Kronecker delta functions appear as a result of using coordinates (t, x i ), where different t values label the spacelike preferred foliations. The extrinsic curvature can be written in terms of a covariant derivative of n µ as K ij = δ µ i δ ν j D ν n µ . It is reasonable to assume that in the matter sector, S matter , diffeomorphism invariance is also explicitly broken to FPdiffs. This allows terms similar to those in the gravity sector, where individually they are scalars under FPdiffs, but where they break four-dimensional diffeomorphisms. In a low-energy limit that agrees with GR, some subset of these terms would need to combine in such a way that reproduces standard covariant matter-sector terms of dimension four or less. Any residual diffeomorphism-breaking terms would have to have small couplings to be consistent with phenomenology.

Scalar and Vector Fields
As an example, consider a scalar field of mass m in a gravitational theory. In terms of ADM variables for the metric, the standard kinetic term for a scalar can be rewritten as whereφ = ∂ 0 φ. Together, the two terms on the right are diffeomorphism-invariant, while individually, each term is separately invariant under FPdiffs. Thus, in a theory with a preferred foliation where FPdiffs is the fundamental symmetry, a different linear combination of these terms can be used in the Lagrangian. In this case, the action for a scalar field in a Hořava-based model in the IR limit is where c 1 and c 2 are introduced as weighting parameters. Making a projection with the normal vector, the second term can be rewritten as Then, by combining (40) and (42), the action S scalar becomes A rescaling of the field φ and the mass m can then be made, which effectively sets c 2 = 1 and renames c 1 as λ φ , giving the final form of the scalar action as Thus, the parameter (1 − λ φ ) can be used to measure the extent of the diffeomorphism breaking in the matter sector, similar to how (1 − λ g ) is used in the gravity sector. Since λ φ can run as the energy scale changes, just as λ g does in the gravity sector, agreement with GR and the SM in the IR limit requires (1 − λ φ ) ≪ 1.
The Higgs sector of the SME has a coefficient (k φφ ) µν that can be matched with this coupling. The correspondence is However, with three unbroken spatial diffeomorphisms, the spatial components of n µ can be gauged to zero in any coordinate frame, including the Sun-centered celestial equatorial frame. With this choice, only the purely timelike component of the SME coefficient (k φφ ) µν remains nonzero and is given as While bounds have been obtained for (k φφ ) 00 in the Higgs sector, the experiments done to date have all assumed a Minkowski spacetime and ignore gravitational effects. For this reason, there are currently no experimental bounds that can be placed on (1 − λ φ ) for the Higgs sector, since sensitivity to both gravity and the Higgs is required.
A similar construction can be made for photon fields in the matter sector. First, a decomposition of the usual four-dimensional Lagrangian term for a massless photon can be made using the ADM formalism, and for simplicity, the unbroken spatial diffeomorphisms can be used to set N i = 0. This reduces the usual photon Lagrangian term to Next, since the decomposition splits the Lagrangian into sectors that are individually scalars under FPdiffs, a linear combination of the terms in (47) can be used in a gravity theory with a preferred foliation. The resulting action in the IR limit is Using projections defined with respect to the normal to the preferred foliation, the following expressions are obtained, Combining these with (47) allows (48) where the SME coefficient (k F ) κλµν in this context is given as (k F ) κλµν = 1 2 (1 − λ γ )[g κµ n λ n ν − g κν n λ n µ − g λµ n κ n ν + g λν n κ n µ ].
Here, n ν acts effectively as a partially fixed background that does not transform under diffeomorphisms with ∂ i ξ 0 = 0. With the gauge fixing setting n µ = ( 1 N , 0), the only nonzero SME components are Thus, with |g ij | ≃ 1 and N ≃ 1, gravity experiments with sensitivity to (k F ) 0i0j can be used to place bounds on the parrameter (1 − λ γ ). In this example, there are experiments with sensitivity to both gravity and electromagnetism. These are the atom interferometry tests, which place bounds at a level of 10 −9 on the quantities σ JK in the Suncentered celestial equatorial frame [48,49]. With these experimental results, the bounds on photon-gravity couplings in the Hořava-based model are estimated as |1 − λ γ | ∼ < 10 −9 .

Vierbein Description and Fermions
To include matter couplings to fermions, a vierbein formalism can be used. The vierbein e a µ provides a link between tensor components in a local Lorentz frame with those in the spacetime frame. It can be used on the normal vector, n µ , which is orthogonal to the preferred foliations in a Hořava-based model, to find the corresponding components n a in the local frames. Since there is freedom to choose and align the orientation of the local Lorentz bases, a simplifying choice is to pick one that gives n a = (1, 0, 0, 0) and n a = η ab n b = (−1, 0, 0, 0), which together obey n a n a = −1. Using these as projections, the vierbein in four dimensions can be decomposed into components with time indices and a three-dimensional spatial subblock. The local indices can be distinguished as a = 0, A with A = 1, 2, 3 labeling the local spatial components. With these, a three-dimensional spatial vierbein e A i is defined such that it equals e a µ when µ = i and a = A. The four-dimensional inverse vierbein can then be decomposed as e µ a = δ µ i δ A a e i A − n a n µ .
where e i A is the three-dimensional inverse of e A i . When this is substituted into g µν = e µ a e ν b η ab , it reproduces (39). Note also that e i A e j B η AB = g ij . In local Lorentz frames, n a remains orthogonal to the preferred spacelike foliation in the same way that n µ remains orthogonal to the preferred foliation in the spacetime frame. For this reason, n a must stay fixed under local Lorentz boosts, since otherwise its orthogonality to the spatial foliations cannot be maintained. In this way, local Lorentz transformations are explicitly broken when a vierbein description is used that introduces n a , similarly to how n µ explicitly breaks diffeomorphisms. However, spatial rotations in the local frame remain symmetrical, which is similar to how FPdiffs remain symmetrical in the spacetime description.
The usual relativistic kinetic Lagrangian term for a fermion ψ with mass m is given as To obtain terms that are suitable for describing a fermion in a gravity model based on the IR limit of Hořava gravity, this fully relativistic term can be subdivided into two parts, where each part, L 1 and L 2 , is invariant individually under FPdiffs and local spatial rotations. The contractions in the purely spatial subblock define L 2 as while the terms that come from summing out the remaining components in µ and a in (55) define L 1 , which for simplicity is not written out here. Using (54), L 2 becomes L 2 = i(e µ a + n a n µ )ψγ a D µ ψ = L ψ + in a n µψ γ a D µ ψ.
Since L 1 and L 2 are each separately scalars under FPdiffs and local spatial rotations, a linear combination of these terms can be used in a Hořava-based model for a fermion, This can be rewritten as L fermion = c 1 (L ψ − L 2 ) + c 2 L 2 = c 1 L ψ + (c 2 − c 1 )L 2 = c 1 L ψ + (c 2 − c 1 )(L ψ + in a n µψ γ a D µ ψ) (60) = c 2 L ψ + (c 2 − c 1 )in a n µψ γ a D µ ψ.
Matching this to the SME term identifies the SME parameter for the Hořava-based model as (1 − λ ψ )n a n µ = −e µ b c αβ e βb e α a , or c µν = −(1 − λ ψ )e a µ n a n ν .
Experiments with sensitivity to matter-gravity couplings involving protons, neutrons, and electrons have been analyzed in [45,85]. They obtain bounds on the SME parameter (c) T T for the time components in the Sun-centered celestial equatorial frame in the range of 10 −6 −10 −8 , depending on the type of experiment and assumptions concerning nuclear binding energies. These can be used to estimate a bound on the coupling (1 − λ ψ ), which limits it to the range |1 − λ ψ | ∼ < 10 −6 − 10 −8 .

VI. SUMMARY AND CONCLUSIONS
The question of whether the SME based on the process of spontaneous diffeomorphism breaking can also be used to analyze models with explicit spacetime symmetry breaking has been investigated. The answer hinges on whether consistency between the Bianchi identities, the equations of motion, and covariant energy-momentum conservation can be maintained. Identities relating these can be obtained for theories that maintain general coordinate invariance, and they reveal that, in general, consistency can be maintained as long as the additional metric modes that come from the diffeomorphism breaking are not suppressed in the consistency conditions. In perturbative treatments, such as in a first-order post-Newtonian limit, the additional metric excitations are suppressed, and therefore there is no useful post-Newtonian limit. However, when matter-gravity couplings are considered, the result is different. Using a Stückelberg approach, the additional excitations in a theory with explicit breaking due to a fixed background tensor are found to be the same as the NG excitations in a corresponding theory with spontaneous breaking. Therefore, the same methodology can be applied in theories with explicit diffeomorphism breaking, and the SME can be used as a phenomenological framework for making comparisons with experiments.
Two examples of theories with explicit diffeomorphism breaking have been examined, where possible matter-gravity couplings have been investigated using the SME. In massive gravity, experiments with sensitivity to gravitational and electromagnetic interactions place bounds on the order of 10 −9 on possible photon couplings. Similarly, a model based on the IR limit of Hořava gravity includes matter-gravity terms that can be compared to the SME. In this case, sensitivity to photon couplings on the order of 10 −9 and to fermion couplings in the range 10 −6 -10 −8 can be attained.