Location via proxy:   [ UP ]  
[Report a bug]   [Manage cookies]                

Cosmological constant as an integration constant

Justin C. Feng    Pisin Chen Leung Center for Cosmology and Particle Astrophysics, National Taiwan University, No.1 Sec.4, Roosevelt Rd., Taipei 10617, Taiwan, R.O.C.
Abstract

The discrepancy between the observed value of the cosmological constant (CC) and its expected value from quantum field theoretical considerations motivates the search for a theory in which the CC is decoupled from the vacuum energy. In this article, we consider theories in which the CC is regarded as an integration constant. These theories include trace-free Einstein gravity, theories constructed from the Codazzi equation (which includes Cotton gravity and a gauge-gravity inspired theory), and conformal Killing gravity. We remark on a recent debate regarding the viability of Cotton gravity and find that while the Codazzi equation is indeed underdetermined, the solutions of the Codazzi equation trivialize to λgab𝜆subscript𝑔𝑎𝑏\lambda g_{ab}italic_λ italic_g start_POSTSUBSCRIPT italic_a italic_b end_POSTSUBSCRIPT on generic backgrounds, and that in principle, one can close the system with the divergence-free condition and an appropriate choice of initial data. We also propose a full variational principle (full in the sense that variations in all variables are considered) for each of the aforementioned theories that can incorporate the matter sector; in this manner, we can obtain the trace-free Einstein equations without a unimodular constraint.

Introduction — Since Einstein introduced the cosmological constant (CC) to his equations in 1917 [1], the debate about its nature has never stopped. Before the discovery of the accelerating expansion of the universe in 1999 [2, 3], the main consideration was why the CC is not exactly zero (See, for example, the famous review by Weinberg[4]). In the aftermath of the discovery of the accelerating expansion of the universe, the problem became more acute. The notion of dark energy introduced to account for this effect behaves like the CC. Indeed, in most analyses in cosmology, such as that of the cosmic microwave background, the CC is commonly invoked in the modeling. If true, then the observed value for the CC is tiny compared with the microscopic interpretation of the CC a la quantum field theory. Why is the CC not only nonzero but also so tiny? Such a seeming dilemma further motivates the reconstruction of Einstein gravity.

Such a reformulation should somehow decouple the CC from microscopic physics. One strategy is to seek a formulation of Einstein gravity in which the CC arises as a constant of integration. At the level of the field equations, there are several proposals in the literature. A proposal originating from Einstein [5] is that of the trace-free Einstein equations, which are obtained in unimodular gravity [6, *vanderBij:1981ym, *Buchmuller:1988wx, *Buchmuller:1988yn, *Henneaux:1989zc, 4] (see [11, *Alvarez:2023utn] for recent overviews of the topic), but can be considered on their own without the constraint on the metric determinant [13, 14]. Another proposal by Cook and Chen [15, 16], inspired by a set of third-order equations developed by Kilmister, Newman [17, *Kilmister:1961], and later by Yang [19] is to construct a gravity theory from the Bianchi identities and derivatives of the energy-momentum tensor. We focus primarily on the aforementioned approaches, though we mention here other proposals in which one may obtain the unconstrained Einstein field equations with the CC arising as an integration constant, such as the thermodynamical approach of Padmanabhan [20, 21] (see also [22]), a modified metric-affine theory [23, *Tapia:1998qt, *Tapia:1998as], and contributions from a three-form action 111However, we note that in this case the vacuum energy still contributes to the cosmological term. [26, 27]. More recently, Harada proposed two theories of gravity, termed Cotton gravity [28] and the so-called conformal Killing gravity [29, 30]. Cotton gravity is constructed from the Cotton tensor, and it has been shown that one can rewrite the field equations of Cotton gravity as a Codazzi equation [aCb]c=0\nabla_{[a}C_{b]c}=0∇ start_POSTSUBSCRIPT [ italic_a end_POSTSUBSCRIPT italic_C start_POSTSUBSCRIPT italic_b ] italic_c end_POSTSUBSCRIPT = 0 [31, 32]; in this regard, Cotton gravity is in the same class as the earlier theory of Cook and Chen, as we will show that the field equations for the latter can also be reformulated as a Codazzi equation. Conformal Killing gravity, on the other hand, has a different structure altogether, as the field equations have the form of the defining equation for a gradient conformal Killing tensor.

Cotton gravity was recently the focus of a debate [33, 34, 35, 36] over its viability as a physical theory, as it was discovered that the symmetry-reduced equations of Cotton gravity are underdetermined. It is been argued in [34, 36] that while the equations of Cotton gravity are indeed underdetermined in cases of high symmetry (for instance, cosmological or spherically symmetric spacetimes), this is not necessarily the case in generic spacetimes. Similar concerns apply to the theory of Cook and Chen, as they share the same Codazzi structure. Another concern that applies to Codazzi type theories and conformal Killing gravity, is the fact that they are only defined at the level of the field equations; these theories suffer from the lack of a complete variational principle, as the variational principles proposed so far require holding some of the variables (such as the spacetime metric gabsubscript𝑔𝑎𝑏g_{ab}italic_g start_POSTSUBSCRIPT italic_a italic_b end_POSTSUBSCRIPT) fixed. One difficulty is that the energy-momentum tensor is obtained from variations of the metric tensor in the matter sector, but derivatives of the energy-momentum tensor appear in the field equations. The primary aim of this article is to shed light on the debate over the Codazzi equation and to propose a variational principle (including matter) for the trace-free Einstein equations, the Codazzi equation, and the conformal Killing equation without introducing additional constraints.

Trace-free Einstein equations — The trace-free Einstein equations may be obtained by assuming that the determinant of the metric is held fixed under the variation of the action; in this form, the theory is called unimodular gravity. However, the equations themselves can be written without assuming such a constraint [13, 14]:

Rμν14Rgμν=κ(Tμν14Tgμν),μTμν=0,formulae-sequencesubscript𝑅𝜇𝜈14𝑅subscript𝑔𝜇𝜈𝜅subscript𝑇𝜇𝜈14𝑇subscript𝑔𝜇𝜈superscript𝜇subscript𝑇𝜇𝜈0R_{\mu\nu}-\frac{1}{4}Rg_{\mu\nu}=\kappa(T_{\mu\nu}-\frac{1}{4}Tg_{\mu\nu}),% \qquad\nabla^{\mu}T_{\mu\nu}=0,italic_R start_POSTSUBSCRIPT italic_μ italic_ν end_POSTSUBSCRIPT - divide start_ARG 1 end_ARG start_ARG 4 end_ARG italic_R italic_g start_POSTSUBSCRIPT italic_μ italic_ν end_POSTSUBSCRIPT = italic_κ ( italic_T start_POSTSUBSCRIPT italic_μ italic_ν end_POSTSUBSCRIPT - divide start_ARG 1 end_ARG start_ARG 4 end_ARG italic_T italic_g start_POSTSUBSCRIPT italic_μ italic_ν end_POSTSUBSCRIPT ) , ∇ start_POSTSUPERSCRIPT italic_μ end_POSTSUPERSCRIPT italic_T start_POSTSUBSCRIPT italic_μ italic_ν end_POSTSUBSCRIPT = 0 , (1)

where R¯μνsubscript¯𝑅𝜇𝜈\underline{R}_{\mu\nu}under¯ start_ARG italic_R end_ARG start_POSTSUBSCRIPT italic_μ italic_ν end_POSTSUBSCRIPT is the Ricci tensor and μTμν=0superscript𝜇subscript𝑇𝜇𝜈0\nabla^{\mu}T_{\mu\nu}=0∇ start_POSTSUPERSCRIPT italic_μ end_POSTSUPERSCRIPT italic_T start_POSTSUBSCRIPT italic_μ italic_ν end_POSTSUBSCRIPT = 0 is imposed as an independent condition on the energy-momentum tensor. The contracted Bianchi identity μ(Rμν12Rgμν)=0superscript𝜇subscript𝑅𝜇𝜈12𝑅subscript𝑔𝜇𝜈0\nabla^{\mu}(R_{\mu\nu}-\tfrac{1}{2}Rg_{\mu\nu})=0∇ start_POSTSUPERSCRIPT italic_μ end_POSTSUPERSCRIPT ( italic_R start_POSTSUBSCRIPT italic_μ italic_ν end_POSTSUBSCRIPT - divide start_ARG 1 end_ARG start_ARG 2 end_ARG italic_R italic_g start_POSTSUBSCRIPT italic_μ italic_ν end_POSTSUBSCRIPT ) = 0 combined with Eq.(1) implies μ(R+κT)=0subscript𝜇𝑅𝜅𝑇0\nabla_{\mu}(R+\kappa T)=0∇ start_POSTSUBSCRIPT italic_μ end_POSTSUBSCRIPT ( italic_R + italic_κ italic_T ) = 0, which has the solution R+κT=4Λ¯𝑅𝜅𝑇4¯ΛR+\kappa T=4\bar{\Lambda}italic_R + italic_κ italic_T = 4 over¯ start_ARG roman_Λ end_ARG, where Λ¯¯Λ\bar{\Lambda}over¯ start_ARG roman_Λ end_ARG is an integration constant. Solving for T𝑇Titalic_T and plugging this back into Eq. (1) yields:

Rμν12Rgμν+Λ¯gμν=κTμν.subscript𝑅𝜇𝜈12𝑅subscript𝑔𝜇𝜈¯Λsubscript𝑔𝜇𝜈𝜅subscript𝑇𝜇𝜈R_{\mu\nu}-\frac{1}{2}Rg_{\mu\nu}+\bar{\Lambda}g_{\mu\nu}=\kappa T_{\mu\nu}.italic_R start_POSTSUBSCRIPT italic_μ italic_ν end_POSTSUBSCRIPT - divide start_ARG 1 end_ARG start_ARG 2 end_ARG italic_R italic_g start_POSTSUBSCRIPT italic_μ italic_ν end_POSTSUBSCRIPT + over¯ start_ARG roman_Λ end_ARG italic_g start_POSTSUBSCRIPT italic_μ italic_ν end_POSTSUBSCRIPT = italic_κ italic_T start_POSTSUBSCRIPT italic_μ italic_ν end_POSTSUBSCRIPT . (2)

Later, we consider a variational principle that does not impose a unimodular constraint on the metric.

Bianchi gravity — The works [15, 16] consider a field theory in which the connection and curvature tensors are treated as analogues of the respective gauge field and field strength tensor of electrodynamics, yielding:

ρRμνστ+τRμνρσ+σRμντρ=0subscript𝜌subscript𝑅𝜇𝜈𝜎𝜏subscript𝜏subscript𝑅𝜇𝜈𝜌𝜎subscript𝜎subscript𝑅𝜇𝜈𝜏𝜌0\nabla_{\rho}R_{\mu\nu\sigma\tau}+\nabla_{\tau}R_{\mu\nu\rho\sigma}+\nabla_{% \sigma}R_{\mu\nu\tau\rho}=0∇ start_POSTSUBSCRIPT italic_ρ end_POSTSUBSCRIPT italic_R start_POSTSUBSCRIPT italic_μ italic_ν italic_σ italic_τ end_POSTSUBSCRIPT + ∇ start_POSTSUBSCRIPT italic_τ end_POSTSUBSCRIPT italic_R start_POSTSUBSCRIPT italic_μ italic_ν italic_ρ italic_σ end_POSTSUBSCRIPT + ∇ start_POSTSUBSCRIPT italic_σ end_POSTSUBSCRIPT italic_R start_POSTSUBSCRIPT italic_μ italic_ν italic_τ italic_ρ end_POSTSUBSCRIPT = 0 (3)
σRσρ=μνκ2Jμνρ\nabla_{\sigma}R^{\sigma\rho}{}_{\mu\nu}=\tfrac{\kappa}{2}J_{\mu\nu}{}^{\rho}∇ start_POSTSUBSCRIPT italic_σ end_POSTSUBSCRIPT italic_R start_POSTSUPERSCRIPT italic_σ italic_ρ end_POSTSUPERSCRIPT start_FLOATSUBSCRIPT italic_μ italic_ν end_FLOATSUBSCRIPT = divide start_ARG italic_κ end_ARG start_ARG 2 end_ARG italic_J start_POSTSUBSCRIPT italic_μ italic_ν end_POSTSUBSCRIPT start_FLOATSUPERSCRIPT italic_ρ end_FLOATSUPERSCRIPT (4)

where Eq. (3) is the Bianchi identity and Rμνστsubscript𝑅𝜇𝜈𝜎𝜏R_{\mu\nu\sigma\tau}italic_R start_POSTSUBSCRIPT italic_μ italic_ν italic_σ italic_τ end_POSTSUBSCRIPT is the curvature tensor. One postulates the source JμνρJ_{\mu\nu}{}^{\rho}italic_J start_POSTSUBSCRIPT italic_μ italic_ν end_POSTSUBSCRIPT start_FLOATSUPERSCRIPT italic_ρ end_FLOATSUPERSCRIPT to be integrable, so that for some symmetric rank-2 tensor T¯μνsubscript¯𝑇𝜇𝜈\bar{T}_{\mu\nu}over¯ start_ARG italic_T end_ARG start_POSTSUBSCRIPT italic_μ italic_ν end_POSTSUBSCRIPT, it may be written as:

Jμν=ρq[μT¯ν].ρJ_{\mu\nu}{}^{\rho}=-q\nabla_{[\mu}\bar{T}_{\nu]}{{}^{\rho}}.italic_J start_POSTSUBSCRIPT italic_μ italic_ν end_POSTSUBSCRIPT start_FLOATSUPERSCRIPT italic_ρ end_FLOATSUPERSCRIPT = - italic_q ∇ start_POSTSUBSCRIPT [ italic_μ end_POSTSUBSCRIPT over¯ start_ARG italic_T end_ARG start_POSTSUBSCRIPT italic_ν ] end_POSTSUBSCRIPT start_FLOATSUPERSCRIPT italic_ρ end_FLOATSUPERSCRIPT . (5)

The field equation (4) may be rewritten:

[μRν]ρκ[μT¯ν]=ρ0,{\nabla}_{[\mu}R_{\nu]}{}{{}^{\rho}}-{\kappa}\nabla_{[\mu}\bar{T}_{\nu]}{{}^{% \rho}}=0,∇ start_POSTSUBSCRIPT [ italic_μ end_POSTSUBSCRIPT italic_R start_POSTSUBSCRIPT italic_ν ] end_POSTSUBSCRIPT start_FLOATSUPERSCRIPT italic_ρ end_FLOATSUPERSCRIPT - italic_κ ∇ start_POSTSUBSCRIPT [ italic_μ end_POSTSUBSCRIPT over¯ start_ARG italic_T end_ARG start_POSTSUBSCRIPT italic_ν ] end_POSTSUBSCRIPT start_FLOATSUPERSCRIPT italic_ρ end_FLOATSUPERSCRIPT = 0 , (6)

where T¯μν=Tμν12gμνTsubscript¯𝑇𝜇𝜈subscript𝑇𝜇𝜈12subscript𝑔𝜇𝜈𝑇\bar{T}_{\mu\nu}={T}_{\mu\nu}-\tfrac{1}{2}g_{\mu\nu}Tover¯ start_ARG italic_T end_ARG start_POSTSUBSCRIPT italic_μ italic_ν end_POSTSUBSCRIPT = italic_T start_POSTSUBSCRIPT italic_μ italic_ν end_POSTSUBSCRIPT - divide start_ARG 1 end_ARG start_ARG 2 end_ARG italic_g start_POSTSUBSCRIPT italic_μ italic_ν end_POSTSUBSCRIPT italic_T is identified with the trace-reversal of the energy-momentum tensor. One class of solutions for Eq. (6) corresponds to the trace-reversed Einstein equation, up to some integration tensor Xμνsubscript𝑋𝜇𝜈X_{\mu\nu}italic_X start_POSTSUBSCRIPT italic_μ italic_ν end_POSTSUBSCRIPT:

Rμν=κ(Tμν12gμνT+Xμν),subscript𝑅𝜇𝜈𝜅subscript𝑇𝜇𝜈12subscript𝑔𝜇𝜈𝑇subscript𝑋𝜇𝜈R_{\mu\nu}=\kappa({T}_{\mu\nu}-\tfrac{1}{2}g_{\mu\nu}T+{X}_{\mu\nu}),italic_R start_POSTSUBSCRIPT italic_μ italic_ν end_POSTSUBSCRIPT = italic_κ ( italic_T start_POSTSUBSCRIPT italic_μ italic_ν end_POSTSUBSCRIPT - divide start_ARG 1 end_ARG start_ARG 2 end_ARG italic_g start_POSTSUBSCRIPT italic_μ italic_ν end_POSTSUBSCRIPT italic_T + italic_X start_POSTSUBSCRIPT italic_μ italic_ν end_POSTSUBSCRIPT ) , (7)

where Xμνsubscript𝑋𝜇𝜈X_{\mu\nu}italic_X start_POSTSUBSCRIPT italic_μ italic_ν end_POSTSUBSCRIPT satisfies the following:

Xμν=Xνμ,μXμν=0,[μXν]ρ=0.X_{\mu\nu}=X_{\nu\mu},\qquad\nabla^{\mu}X_{\mu\nu}=0,\qquad\nabla_{[\mu}X_{\nu% ]\rho}=0.italic_X start_POSTSUBSCRIPT italic_μ italic_ν end_POSTSUBSCRIPT = italic_X start_POSTSUBSCRIPT italic_ν italic_μ end_POSTSUBSCRIPT , ∇ start_POSTSUPERSCRIPT italic_μ end_POSTSUPERSCRIPT italic_X start_POSTSUBSCRIPT italic_μ italic_ν end_POSTSUBSCRIPT = 0 , ∇ start_POSTSUBSCRIPT [ italic_μ end_POSTSUBSCRIPT italic_X start_POSTSUBSCRIPT italic_ν ] italic_ρ end_POSTSUBSCRIPT = 0 . (8)

The tensor Xμνsubscript𝑋𝜇𝜈X_{\mu\nu}italic_X start_POSTSUBSCRIPT italic_μ italic_ν end_POSTSUBSCRIPT is assumed to be determined by boundary and initial data; if the initial/boundary data yields Xμν=λgμνsubscript𝑋𝜇𝜈𝜆subscript𝑔𝜇𝜈X_{\mu\nu}=\lambda g_{\mu\nu}italic_X start_POSTSUBSCRIPT italic_μ italic_ν end_POSTSUBSCRIPT = italic_λ italic_g start_POSTSUBSCRIPT italic_μ italic_ν end_POSTSUBSCRIPT, then the CC λ𝜆\lambdaitalic_λ may be interpreted as an integration constant (we discuss later how such a solution arises from initial data).

One might imagine rewriting Eq. (6) as:

[μKν]=ρ0,{\nabla}_{[\mu}K_{\nu]}{}{{}^{\rho}}=0,∇ start_POSTSUBSCRIPT [ italic_μ end_POSTSUBSCRIPT italic_K start_POSTSUBSCRIPT italic_ν ] end_POSTSUBSCRIPT start_FLOATSUPERSCRIPT italic_ρ end_FLOATSUPERSCRIPT = 0 , (9)

where:

Kμν=Rμνκ(Tμν12gμνT+Xμν).subscript𝐾𝜇𝜈subscript𝑅𝜇𝜈𝜅subscript𝑇𝜇𝜈12subscript𝑔𝜇𝜈𝑇subscript𝑋𝜇𝜈K_{\mu\nu}=R_{\mu\nu}-\kappa({T}_{\mu\nu}-\tfrac{1}{2}g_{\mu\nu}T+{X}_{\mu\nu}).italic_K start_POSTSUBSCRIPT italic_μ italic_ν end_POSTSUBSCRIPT = italic_R start_POSTSUBSCRIPT italic_μ italic_ν end_POSTSUBSCRIPT - italic_κ ( italic_T start_POSTSUBSCRIPT italic_μ italic_ν end_POSTSUBSCRIPT - divide start_ARG 1 end_ARG start_ARG 2 end_ARG italic_g start_POSTSUBSCRIPT italic_μ italic_ν end_POSTSUBSCRIPT italic_T + italic_X start_POSTSUBSCRIPT italic_μ italic_ν end_POSTSUBSCRIPT ) . (10)

Equation (9) has the form of a Codazzi equation [31, 32], which is the subject of recent interest in the Cotton gravity program [28].

Cotton gravity — Cotton gravity is based on the Cotton tensor:

Cabc:=bRaccRab16(gacbRgabcR)assignsubscript𝐶𝑎𝑏𝑐subscript𝑏subscript𝑅𝑎𝑐subscript𝑐subscript𝑅𝑎𝑏16subscript𝑔𝑎𝑐subscript𝑏𝑅subscript𝑔𝑎𝑏subscript𝑐𝑅C_{abc}:=\nabla_{b}R_{ac}-\nabla_{c}R_{ab}-\tfrac{1}{6}\left(g_{ac}\nabla_{b}R% -g_{ab}\nabla_{c}R\right)italic_C start_POSTSUBSCRIPT italic_a italic_b italic_c end_POSTSUBSCRIPT := ∇ start_POSTSUBSCRIPT italic_b end_POSTSUBSCRIPT italic_R start_POSTSUBSCRIPT italic_a italic_c end_POSTSUBSCRIPT - ∇ start_POSTSUBSCRIPT italic_c end_POSTSUBSCRIPT italic_R start_POSTSUBSCRIPT italic_a italic_b end_POSTSUBSCRIPT - divide start_ARG 1 end_ARG start_ARG 6 end_ARG ( italic_g start_POSTSUBSCRIPT italic_a italic_c end_POSTSUBSCRIPT ∇ start_POSTSUBSCRIPT italic_b end_POSTSUBSCRIPT italic_R - italic_g start_POSTSUBSCRIPT italic_a italic_b end_POSTSUBSCRIPT ∇ start_POSTSUBSCRIPT italic_c end_POSTSUBSCRIPT italic_R ) (11)

Defining a new tensor 𝒯abcsubscript𝒯𝑎𝑏𝑐\mathscr{T}_{abc}script_T start_POSTSUBSCRIPT italic_a italic_b italic_c end_POSTSUBSCRIPT

𝒯abc:=bTaccTab16(gacbTgabcT),assignsubscript𝒯𝑎𝑏𝑐subscript𝑏subscript𝑇𝑎𝑐subscript𝑐subscript𝑇𝑎𝑏16subscript𝑔𝑎𝑐subscript𝑏𝑇subscript𝑔𝑎𝑏subscript𝑐𝑇\mathscr{T}_{abc}:=\nabla_{b}T_{ac}-\nabla_{c}T_{ab}-\tfrac{1}{6}\left(g_{ac}% \nabla_{b}T-g_{ab}\nabla_{c}T\right),script_T start_POSTSUBSCRIPT italic_a italic_b italic_c end_POSTSUBSCRIPT := ∇ start_POSTSUBSCRIPT italic_b end_POSTSUBSCRIPT italic_T start_POSTSUBSCRIPT italic_a italic_c end_POSTSUBSCRIPT - ∇ start_POSTSUBSCRIPT italic_c end_POSTSUBSCRIPT italic_T start_POSTSUBSCRIPT italic_a italic_b end_POSTSUBSCRIPT - divide start_ARG 1 end_ARG start_ARG 6 end_ARG ( italic_g start_POSTSUBSCRIPT italic_a italic_c end_POSTSUBSCRIPT ∇ start_POSTSUBSCRIPT italic_b end_POSTSUBSCRIPT italic_T - italic_g start_POSTSUBSCRIPT italic_a italic_b end_POSTSUBSCRIPT ∇ start_POSTSUBSCRIPT italic_c end_POSTSUBSCRIPT italic_T ) , (12)

one can write the governing equation for Cotton gravity as:

Cabc=2κ𝒯abc.subscript𝐶𝑎𝑏𝑐2𝜅subscript𝒯𝑎𝑏𝑐C_{abc}=2\kappa\mathscr{T}_{abc}.italic_C start_POSTSUBSCRIPT italic_a italic_b italic_c end_POSTSUBSCRIPT = 2 italic_κ script_T start_POSTSUBSCRIPT italic_a italic_b italic_c end_POSTSUBSCRIPT . (13)

This equation can be rewritten in the form:

[a𝒞b]=c0,{\nabla}_{[a}\mathcal{C}_{b]}{}{{}^{c}}=0,∇ start_POSTSUBSCRIPT [ italic_a end_POSTSUBSCRIPT caligraphic_C start_POSTSUBSCRIPT italic_b ] end_POSTSUBSCRIPT start_FLOATSUPERSCRIPT italic_c end_FLOATSUPERSCRIPT = 0 , (14)

where:

𝒞ab=GabκTab13(GκT)gab+Xab.subscript𝒞𝑎𝑏subscript𝐺𝑎𝑏𝜅subscript𝑇𝑎𝑏13𝐺𝜅𝑇subscript𝑔𝑎𝑏subscript𝑋𝑎𝑏\mathcal{C}_{ab}=G_{ab}-\kappa T_{ab}-\frac{1}{3}(G-\kappa T)g_{ab}+X_{ab}.caligraphic_C start_POSTSUBSCRIPT italic_a italic_b end_POSTSUBSCRIPT = italic_G start_POSTSUBSCRIPT italic_a italic_b end_POSTSUBSCRIPT - italic_κ italic_T start_POSTSUBSCRIPT italic_a italic_b end_POSTSUBSCRIPT - divide start_ARG 1 end_ARG start_ARG 3 end_ARG ( italic_G - italic_κ italic_T ) italic_g start_POSTSUBSCRIPT italic_a italic_b end_POSTSUBSCRIPT + italic_X start_POSTSUBSCRIPT italic_a italic_b end_POSTSUBSCRIPT . (15)

“Conformal” Killing gravity — The equations for the so-called conformal Killing gravity theory proposed in [29] may be written in the form [30]:

6(Kbc;a+Kca;b+Kab;c)=gbcK,a+gcaK,b+gabK,c,\displaystyle 6(K_{bc;a}+K_{ca;b}+K_{ab;c})=g_{bc}K_{,a}+g_{ca}K_{,b}+g_{ab}K_% {,c},6 ( italic_K start_POSTSUBSCRIPT italic_b italic_c ; italic_a end_POSTSUBSCRIPT + italic_K start_POSTSUBSCRIPT italic_c italic_a ; italic_b end_POSTSUBSCRIPT + italic_K start_POSTSUBSCRIPT italic_a italic_b ; italic_c end_POSTSUBSCRIPT ) = italic_g start_POSTSUBSCRIPT italic_b italic_c end_POSTSUBSCRIPT italic_K start_POSTSUBSCRIPT , italic_a end_POSTSUBSCRIPT + italic_g start_POSTSUBSCRIPT italic_c italic_a end_POSTSUBSCRIPT italic_K start_POSTSUBSCRIPT , italic_b end_POSTSUBSCRIPT + italic_g start_POSTSUBSCRIPT italic_a italic_b end_POSTSUBSCRIPT italic_K start_POSTSUBSCRIPT , italic_c end_POSTSUBSCRIPT , (16)
Kab:=κTabGab.assignsubscript𝐾𝑎𝑏𝜅subscript𝑇𝑎𝑏subscript𝐺𝑎𝑏\displaystyle K_{ab}:=\kappa T_{ab}-G_{ab}.italic_K start_POSTSUBSCRIPT italic_a italic_b end_POSTSUBSCRIPT := italic_κ italic_T start_POSTSUBSCRIPT italic_a italic_b end_POSTSUBSCRIPT - italic_G start_POSTSUBSCRIPT italic_a italic_b end_POSTSUBSCRIPT .

The equation above has the form of the equation for a conformal Killing tensor:

Kbc;a+Kca;b+Kab;c=uagbc+ubgca+ucgab,subscript𝐾𝑏𝑐𝑎subscript𝐾𝑐𝑎𝑏subscript𝐾𝑎𝑏𝑐subscript𝑢𝑎subscript𝑔𝑏𝑐subscript𝑢𝑏subscript𝑔𝑐𝑎subscript𝑢𝑐subscript𝑔𝑎𝑏K_{bc;a}+K_{ca;b}+K_{ab;c}=u_{a}g_{bc}+u_{b}g_{ca}+u_{c}g_{ab},italic_K start_POSTSUBSCRIPT italic_b italic_c ; italic_a end_POSTSUBSCRIPT + italic_K start_POSTSUBSCRIPT italic_c italic_a ; italic_b end_POSTSUBSCRIPT + italic_K start_POSTSUBSCRIPT italic_a italic_b ; italic_c end_POSTSUBSCRIPT = italic_u start_POSTSUBSCRIPT italic_a end_POSTSUBSCRIPT italic_g start_POSTSUBSCRIPT italic_b italic_c end_POSTSUBSCRIPT + italic_u start_POSTSUBSCRIPT italic_b end_POSTSUBSCRIPT italic_g start_POSTSUBSCRIPT italic_c italic_a end_POSTSUBSCRIPT + italic_u start_POSTSUBSCRIPT italic_c end_POSTSUBSCRIPT italic_g start_POSTSUBSCRIPT italic_a italic_b end_POSTSUBSCRIPT , (17)

where uasubscript𝑢𝑎u_{a}italic_u start_POSTSUBSCRIPT italic_a end_POSTSUBSCRIPT is a one-form; Kabsubscript𝐾𝑎𝑏K_{ab}italic_K start_POSTSUBSCRIPT italic_a italic_b end_POSTSUBSCRIPT is a symmetric Killing tensor when ua=0subscript𝑢𝑎0u_{a}=0italic_u start_POSTSUBSCRIPT italic_a end_POSTSUBSCRIPT = 0. However, as pointed out in [37] a simple redefinition K¯ab:=6KabgabKassignsubscript¯𝐾𝑎𝑏6subscript𝐾𝑎𝑏subscript𝑔𝑎𝑏𝐾\bar{K}_{ab}:=6K_{ab}-g_{ab}Kover¯ start_ARG italic_K end_ARG start_POSTSUBSCRIPT italic_a italic_b end_POSTSUBSCRIPT := 6 italic_K start_POSTSUBSCRIPT italic_a italic_b end_POSTSUBSCRIPT - italic_g start_POSTSUBSCRIPT italic_a italic_b end_POSTSUBSCRIPT italic_K reveals that any nontrivial solution of Eq. (16) requires the existence of a Killing tensor given by K¯absubscript¯𝐾𝑎𝑏\bar{K}_{ab}over¯ start_ARG italic_K end_ARG start_POSTSUBSCRIPT italic_a italic_b end_POSTSUBSCRIPT. It follows that the solutions of Eq. (16) either trivialize to Kab=λgabsubscript𝐾𝑎𝑏𝜆subscript𝑔𝑎𝑏K_{ab}=-\lambda g_{ab}italic_K start_POSTSUBSCRIPT italic_a italic_b end_POSTSUBSCRIPT = - italic_λ italic_g start_POSTSUBSCRIPT italic_a italic_b end_POSTSUBSCRIPT, or admit a Killing tensor. For this reason, it is perhaps more appropriate to refer to this theory as “Killing tensor gravity”, but we retain here the common terminology in the literature. Moreover, we note that perturbations of Killing tensor solutions satisfying Kabλgabsubscript𝐾𝑎𝑏𝜆subscript𝑔𝑎𝑏K_{ab}\neq-\lambda g_{ab}italic_K start_POSTSUBSCRIPT italic_a italic_b end_POSTSUBSCRIPT ≠ - italic_λ italic_g start_POSTSUBSCRIPT italic_a italic_b end_POSTSUBSCRIPT must also admit a Killing tensor; since our universe (being inhomogeneous and anisotropic on small scales) does not have the required symmetries to admit a global Killing tensor, we regard such solutions as unphysical.

The Codazzi equation — Equations (9) and (14) have the form of a Codazzi equation:

[aCb]=c0,{\nabla}_{[a}C_{b]}{}{{}^{c}}=0,∇ start_POSTSUBSCRIPT [ italic_a end_POSTSUBSCRIPT italic_C start_POSTSUBSCRIPT italic_b ] end_POSTSUBSCRIPT start_FLOATSUPERSCRIPT italic_c end_FLOATSUPERSCRIPT = 0 , (18)

where Cab=Cbasubscript𝐶𝑎𝑏subscript𝐶𝑏𝑎C_{ab}=C_{ba}italic_C start_POSTSUBSCRIPT italic_a italic_b end_POSTSUBSCRIPT = italic_C start_POSTSUBSCRIPT italic_b italic_a end_POSTSUBSCRIPT is termed a Codazzi tensor, following the terminology of [31, 32]. One can perhaps imagine a class of theories defined by an equation of the form in Eq. (18). Now there is some debate over whether equations of this type are predictive [33, 34, 35, 36]; in particular, it has been claimed in [35] that the Codazzi equations are underdetermined, and in a responding comment it was claimed that this only occurs for highly symmetric spacetimes [36]. In the following, we attempt to shed some light on this dispute.

Contracting Eq. (18), one obtains:

cCacaC=0.{\nabla}_{c}C_{a}{}{{}^{c}}-{\nabla}_{a}C=0.∇ start_POSTSUBSCRIPT italic_c end_POSTSUBSCRIPT italic_C start_POSTSUBSCRIPT italic_a end_POSTSUBSCRIPT start_FLOATSUPERSCRIPT italic_c end_FLOATSUPERSCRIPT - ∇ start_POSTSUBSCRIPT italic_a end_POSTSUBSCRIPT italic_C = 0 . (19)

where C:=CaaC:=C_{a}{}{{}^{a}}italic_C := italic_C start_POSTSUBSCRIPT italic_a end_POSTSUBSCRIPT start_FLOATSUPERSCRIPT italic_a end_FLOATSUPERSCRIPT. Taking the divergence of Eq. (18), one obtains a wave-like equation:

CabacCcbRacCcbRacdCcdb=0.\Box C_{a}{}{{}^{b}}-\nabla_{a}{\nabla}_{c}C{}^{cb}-R_{ac}C^{cb}-R_{acd}{{}^{b% }}C^{cd}=0.□ italic_C start_POSTSUBSCRIPT italic_a end_POSTSUBSCRIPT start_FLOATSUPERSCRIPT italic_b end_FLOATSUPERSCRIPT - ∇ start_POSTSUBSCRIPT italic_a end_POSTSUBSCRIPT ∇ start_POSTSUBSCRIPT italic_c end_POSTSUBSCRIPT italic_C start_FLOATSUPERSCRIPT italic_c italic_b end_FLOATSUPERSCRIPT - italic_R start_POSTSUBSCRIPT italic_a italic_c end_POSTSUBSCRIPT italic_C start_POSTSUPERSCRIPT italic_c italic_b end_POSTSUPERSCRIPT - italic_R start_POSTSUBSCRIPT italic_a italic_c italic_d end_POSTSUBSCRIPT start_FLOATSUPERSCRIPT italic_b end_FLOATSUPERSCRIPT italic_C start_POSTSUPERSCRIPT italic_c italic_d end_POSTSUPERSCRIPT = 0 . (20)

and the divergence of Eq. (19) is:

CabCab=0.𝐶subscript𝑎subscript𝑏superscript𝐶𝑎𝑏0\Box C-{\nabla}_{a}{\nabla}_{b}C^{ab}=0.□ italic_C - ∇ start_POSTSUBSCRIPT italic_a end_POSTSUBSCRIPT ∇ start_POSTSUBSCRIPT italic_b end_POSTSUBSCRIPT italic_C start_POSTSUPERSCRIPT italic_a italic_b end_POSTSUPERSCRIPT = 0 . (21)

Of course, Eq. (20) is underdetermined due to the term acCcb\nabla_{a}{\nabla}_{c}C{}^{cb}∇ start_POSTSUBSCRIPT italic_a end_POSTSUBSCRIPT ∇ start_POSTSUBSCRIPT italic_c end_POSTSUBSCRIPT italic_C start_FLOATSUPERSCRIPT italic_c italic_b end_FLOATSUPERSCRIPT. In an orthonormal basis, the derivative t2Cbt\partial^{2}_{t}C{{}_{t}}^{b}∂ start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT italic_C start_FLOATSUBSCRIPT italic_t end_FLOATSUBSCRIPT start_POSTSUPERSCRIPT italic_b end_POSTSUPERSCRIPT disappears from Eq. (20), so that the time evolution of the components CbtC{{}_{t}}^{b}italic_C start_FLOATSUBSCRIPT italic_t end_FLOATSUBSCRIPT start_POSTSUPERSCRIPT italic_b end_POSTSUPERSCRIPT is ill-defined—the equations are indeed underdetermined. However, Eq. (20) reveals a straightforward resolution; one should supply the Codazzi equation with a constraint equation, and a natural one is the divergence-free condition aCab=0subscript𝑎superscript𝐶𝑎𝑏0\nabla_{a}C^{ab}=0∇ start_POSTSUBSCRIPT italic_a end_POSTSUBSCRIPT italic_C start_POSTSUPERSCRIPT italic_a italic_b end_POSTSUPERSCRIPT = 0. Of course, an attempt to directly solve the Codazzi equations supplemented with the aCab=0subscript𝑎superscript𝐶𝑎𝑏0\nabla_{a}C^{ab}=0∇ start_POSTSUBSCRIPT italic_a end_POSTSUBSCRIPT italic_C start_POSTSUPERSCRIPT italic_a italic_b end_POSTSUPERSCRIPT = 0 may still lead to ambiguities, but such ambiguities correspond to freedom in choosing initial data for Cabsubscript𝐶𝑎𝑏C_{ab}italic_C start_POSTSUBSCRIPT italic_a italic_b end_POSTSUBSCRIPT and tCabsubscript𝑡subscript𝐶𝑎𝑏\partial_{t}C_{ab}∂ start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT italic_C start_POSTSUBSCRIPT italic_a italic_b end_POSTSUBSCRIPT for the equivalent system:

CabRacCcbRacdCcdb=0,aCab=0.\displaystyle\Box C_{a}{}{{}^{b}}-R_{ac}C^{cb}-R_{acd}{{}^{b}}C^{cd}=0,\qquad% \nabla_{a}C^{ab}=0.□ italic_C start_POSTSUBSCRIPT italic_a end_POSTSUBSCRIPT start_FLOATSUPERSCRIPT italic_b end_FLOATSUPERSCRIPT - italic_R start_POSTSUBSCRIPT italic_a italic_c end_POSTSUBSCRIPT italic_C start_POSTSUPERSCRIPT italic_c italic_b end_POSTSUPERSCRIPT - italic_R start_POSTSUBSCRIPT italic_a italic_c italic_d end_POSTSUBSCRIPT start_FLOATSUPERSCRIPT italic_b end_FLOATSUPERSCRIPT italic_C start_POSTSUPERSCRIPT italic_c italic_d end_POSTSUPERSCRIPT = 0 , ∇ start_POSTSUBSCRIPT italic_a end_POSTSUBSCRIPT italic_C start_POSTSUPERSCRIPT italic_a italic_b end_POSTSUPERSCRIPT = 0 . (22)

Now Eq. (20) follows from Eq. (18). It is perhaps appropriate to consider the application of the D’Alembertian to the following quantity:

Qab:=c2[aCb],cQ_{ab}{{}^{c}}:=2{\nabla}_{[a}C_{b]}{}{{}^{c}},italic_Q start_POSTSUBSCRIPT italic_a italic_b end_POSTSUBSCRIPT start_FLOATSUPERSCRIPT italic_c end_FLOATSUPERSCRIPT := 2 ∇ start_POSTSUBSCRIPT [ italic_a end_POSTSUBSCRIPT italic_C start_POSTSUBSCRIPT italic_b ] end_POSTSUBSCRIPT start_FLOATSUPERSCRIPT italic_c end_FLOATSUPERSCRIPT , (23)

which after using Eq. (20) has the form:

Qcde=subscript𝑄𝑐𝑑𝑒absent\displaystyle\Box Q_{cde}=□ italic_Q start_POSTSUBSCRIPT italic_c italic_d italic_e end_POSTSUBSCRIPT = RcdebaCab+2(bCa)[cRd]bea+4[eRa][cCd]a\displaystyle R_{cdeb}\nabla_{a}C^{ab}+2(\nabla^{b}C^{a}{}_{[c})R_{d]bea}+4% \nabla_{[e}R_{a][c}C_{d]}{}^{a}italic_R start_POSTSUBSCRIPT italic_c italic_d italic_e italic_b end_POSTSUBSCRIPT ∇ start_POSTSUBSCRIPT italic_a end_POSTSUBSCRIPT italic_C start_POSTSUPERSCRIPT italic_a italic_b end_POSTSUPERSCRIPT + 2 ( ∇ start_POSTSUPERSCRIPT italic_b end_POSTSUPERSCRIPT italic_C start_POSTSUPERSCRIPT italic_a end_POSTSUPERSCRIPT start_FLOATSUBSCRIPT [ italic_c end_FLOATSUBSCRIPT ) italic_R start_POSTSUBSCRIPT italic_d ] italic_b italic_e italic_a end_POSTSUBSCRIPT + 4 ∇ start_POSTSUBSCRIPT [ italic_e end_POSTSUBSCRIPT italic_R start_POSTSUBSCRIPT italic_a ] [ italic_c end_POSTSUBSCRIPT italic_C start_POSTSUBSCRIPT italic_d ] end_POSTSUBSCRIPT start_FLOATSUPERSCRIPT italic_a end_FLOATSUPERSCRIPT (24)
+2Cab[dRc]aebQdRcaaeQaceRda\displaystyle+2C^{ab}\nabla_{[d}R_{c]aeb}-Q_{d}{}^{a}{}_{e}R_{ca}-Q_{ace}R_{d}% {}^{a}+ 2 italic_C start_POSTSUPERSCRIPT italic_a italic_b end_POSTSUPERSCRIPT ∇ start_POSTSUBSCRIPT [ italic_d end_POSTSUBSCRIPT italic_R start_POSTSUBSCRIPT italic_c ] italic_a italic_e italic_b end_POSTSUBSCRIPT - italic_Q start_POSTSUBSCRIPT italic_d end_POSTSUBSCRIPT start_FLOATSUPERSCRIPT italic_a end_FLOATSUPERSCRIPT start_FLOATSUBSCRIPT italic_e end_FLOATSUBSCRIPT italic_R start_POSTSUBSCRIPT italic_c italic_a end_POSTSUBSCRIPT - italic_Q start_POSTSUBSCRIPT italic_a italic_c italic_e end_POSTSUBSCRIPT italic_R start_POSTSUBSCRIPT italic_d end_POSTSUBSCRIPT start_FLOATSUPERSCRIPT italic_a end_FLOATSUPERSCRIPT
2QabRcadbe+QdRcaebabQcRdaebab.2superscript𝑄𝑎𝑏subscriptsubscript𝑅𝑐𝑎𝑑𝑏𝑒subscript𝑄𝑑superscriptsubscript𝑅𝑐𝑎𝑒𝑏𝑎𝑏subscript𝑄𝑐superscriptsubscript𝑅𝑑𝑎𝑒𝑏𝑎𝑏\displaystyle-2Q^{ab}{}_{e}R_{cadb}+Q_{d}{}^{ab}R_{caeb}-Q_{c}{}^{ab}R_{daeb}.- 2 italic_Q start_POSTSUPERSCRIPT italic_a italic_b end_POSTSUPERSCRIPT start_FLOATSUBSCRIPT italic_e end_FLOATSUBSCRIPT italic_R start_POSTSUBSCRIPT italic_c italic_a italic_d italic_b end_POSTSUBSCRIPT + italic_Q start_POSTSUBSCRIPT italic_d end_POSTSUBSCRIPT start_FLOATSUPERSCRIPT italic_a italic_b end_FLOATSUPERSCRIPT italic_R start_POSTSUBSCRIPT italic_c italic_a italic_e italic_b end_POSTSUBSCRIPT - italic_Q start_POSTSUBSCRIPT italic_c end_POSTSUBSCRIPT start_FLOATSUPERSCRIPT italic_a italic_b end_FLOATSUPERSCRIPT italic_R start_POSTSUBSCRIPT italic_d italic_a italic_e italic_b end_POSTSUBSCRIPT .

Note that the right-hand side contains at most first derivatives of Cabsubscript𝐶𝑎𝑏C_{ab}italic_C start_POSTSUBSCRIPT italic_a italic_b end_POSTSUBSCRIPT (and also that the divergence-free condition removes the first term). One can then regard Eq. (24) as a wave equation for Qcdesubscript𝑄𝑐𝑑𝑒Q_{cde}italic_Q start_POSTSUBSCRIPT italic_c italic_d italic_e end_POSTSUBSCRIPT, with Cabsubscript𝐶𝑎𝑏C_{ab}italic_C start_POSTSUBSCRIPT italic_a italic_b end_POSTSUBSCRIPT being an independent variable satisfying Eq. (20) and the constraint (24). Now consider initial data consistent with Eq. (18), so that Qcde=0subscript𝑄𝑐𝑑𝑒0Q_{cde}=0italic_Q start_POSTSUBSCRIPT italic_c italic_d italic_e end_POSTSUBSCRIPT = 0 (but not its derivatives). For nonvanishing curvature, the right-hand side will generically produce a source unless Cab=λgabsubscript𝐶𝑎𝑏𝜆subscript𝑔𝑎𝑏C_{ab}=\lambda g_{ab}italic_C start_POSTSUBSCRIPT italic_a italic_b end_POSTSUBSCRIPT = italic_λ italic_g start_POSTSUBSCRIPT italic_a italic_b end_POSTSUBSCRIPT on a Cauchy hypersurface and cCab=0subscript𝑐subscript𝐶𝑎𝑏0\nabla_{c}C_{ab}=0∇ start_POSTSUBSCRIPT italic_c end_POSTSUBSCRIPT italic_C start_POSTSUBSCRIPT italic_a italic_b end_POSTSUBSCRIPT = 0. This indicates that if Cabλgabsubscript𝐶𝑎𝑏𝜆subscript𝑔𝑎𝑏C_{ab}\neq\lambda g_{ab}italic_C start_POSTSUBSCRIPT italic_a italic_b end_POSTSUBSCRIPT ≠ italic_λ italic_g start_POSTSUBSCRIPT italic_a italic_b end_POSTSUBSCRIPT and cCab0subscript𝑐subscript𝐶𝑎𝑏0\nabla_{c}C_{ab}\neq 0∇ start_POSTSUBSCRIPT italic_c end_POSTSUBSCRIPT italic_C start_POSTSUBSCRIPT italic_a italic_b end_POSTSUBSCRIPT ≠ 0 satisfies Eq. (24) on a Cauchy surface, the time evolution of Cabsubscript𝐶𝑎𝑏C_{ab}italic_C start_POSTSUBSCRIPT italic_a italic_b end_POSTSUBSCRIPT under Eq. (20) will in general generate violations of Eq. (24). There may also be special cases (highly symmetric spacetimes, for instance) in which Eq. (18) may admit nontrivial solutions Cabλgabsubscript𝐶𝑎𝑏𝜆subscript𝑔𝑎𝑏C_{ab}\neq\lambda g_{ab}italic_C start_POSTSUBSCRIPT italic_a italic_b end_POSTSUBSCRIPT ≠ italic_λ italic_g start_POSTSUBSCRIPT italic_a italic_b end_POSTSUBSCRIPT, but one cannot expect the existence of solutions to Eq. (18) for generic initial data. It should be noted that Eq. (24) and the conclusions outlined in the previous paragraph are independent of the constraint aCab=0subscript𝑎superscript𝐶𝑎𝑏0\nabla_{a}C^{ab}=0∇ start_POSTSUBSCRIPT italic_a end_POSTSUBSCRIPT italic_C start_POSTSUPERSCRIPT italic_a italic_b end_POSTSUPERSCRIPT = 0.

In short, Eq. (20) indicates that the quantity QabcQ_{ab}{{}^{c}}italic_Q start_POSTSUBSCRIPT italic_a italic_b end_POSTSUBSCRIPT start_FLOATSUPERSCRIPT italic_c end_FLOATSUPERSCRIPT evolves according to Eq. (24), and that generically, the condition Qab=c0Q_{ab}{{}^{c}}=0italic_Q start_POSTSUBSCRIPT italic_a italic_b end_POSTSUBSCRIPT start_FLOATSUPERSCRIPT italic_c end_FLOATSUPERSCRIPT = 0 can only be maintained for initial data also satisfying Cab=λgabsubscript𝐶𝑎𝑏𝜆subscript𝑔𝑎𝑏C_{ab}=\lambda g_{ab}italic_C start_POSTSUBSCRIPT italic_a italic_b end_POSTSUBSCRIPT = italic_λ italic_g start_POSTSUBSCRIPT italic_a italic_b end_POSTSUBSCRIPT and cCab=0subscript𝑐subscript𝐶𝑎𝑏0\nabla_{c}C_{ab}=0∇ start_POSTSUBSCRIPT italic_c end_POSTSUBSCRIPT italic_C start_POSTSUBSCRIPT italic_a italic_b end_POSTSUBSCRIPT = 0, where λ𝜆\lambdaitalic_λ is an arbitrary constant; other solutions will likely require background geometries with special properties, as argued in [34, 36]. We find that both sets of authors, [33, 35] and [34, 36], present valid points: theories based on the Codazzi tensors are indeed underdetermined, but the space of solutions is nonetheless highly restricted in generic settings.

Variational principle — The variational principles presented in the literature for the trace-free Einstein equations [4] and theories admitting a Codazzi equation formulation (Bianchi gravity [15, 16] and Cotton gravity [28]) typically require holding certain quantities fixed; for the trace-free Einstein equations, the metric determinant is held fixed (corresponding to an additional constraint), and for the Codazzi equation, the metric is held fixed (as variations in the metric will generally introduce additional unwanted constraints). We therefore seek a general variational principle that yields the desired field equations (including contributions from the matter sector) without unwanted constraints under a full variation—that is, a variation with respect to all variables appearing in the action.

Consider a simple action of the form for the case of the Codazzi equation:

SC:=assignsubscript𝑆𝐶absent\displaystyle S_{C}:=italic_S start_POSTSUBSCRIPT italic_C end_POSTSUBSCRIPT := Ω|g|d4x[Yabc[aCb]c+ΨabcdefYabcYdef].\displaystyle\int_{\Omega}\sqrt{|g|}d^{4}x[Y^{abc}\nabla_{[a}C_{b]c}+\Psi_{% abcdef}Y^{abc}Y^{def}].∫ start_POSTSUBSCRIPT roman_Ω end_POSTSUBSCRIPT square-root start_ARG | italic_g | end_ARG italic_d start_POSTSUPERSCRIPT 4 end_POSTSUPERSCRIPT italic_x [ italic_Y start_POSTSUPERSCRIPT italic_a italic_b italic_c end_POSTSUPERSCRIPT ∇ start_POSTSUBSCRIPT [ italic_a end_POSTSUBSCRIPT italic_C start_POSTSUBSCRIPT italic_b ] italic_c end_POSTSUBSCRIPT + roman_Ψ start_POSTSUBSCRIPT italic_a italic_b italic_c italic_d italic_e italic_f end_POSTSUBSCRIPT italic_Y start_POSTSUPERSCRIPT italic_a italic_b italic_c end_POSTSUPERSCRIPT italic_Y start_POSTSUPERSCRIPT italic_d italic_e italic_f end_POSTSUPERSCRIPT ] . (25)

where YcabY{{}^{ab}}_{c}italic_Y start_FLOATSUPERSCRIPT italic_a italic_b end_FLOATSUPERSCRIPT start_POSTSUBSCRIPT italic_c end_POSTSUBSCRIPT is an auxiliary field, and Ψabcdef=ΨdefabcsubscriptΨ𝑎𝑏𝑐𝑑𝑒𝑓subscriptΨ𝑑𝑒𝑓𝑎𝑏𝑐\Psi_{abcdef}=\Psi_{defabc}roman_Ψ start_POSTSUBSCRIPT italic_a italic_b italic_c italic_d italic_e italic_f end_POSTSUBSCRIPT = roman_Ψ start_POSTSUBSCRIPT italic_d italic_e italic_f italic_a italic_b italic_c end_POSTSUBSCRIPT is a Lagrange multiplier222ΨabcdefsubscriptΨ𝑎𝑏𝑐𝑑𝑒𝑓\Psi_{abcdef}roman_Ψ start_POSTSUBSCRIPT italic_a italic_b italic_c italic_d italic_e italic_f end_POSTSUBSCRIPT may be regarded as a metric on the field space for Yabcsuperscript𝑌𝑎𝑏𝑐Y^{abc}italic_Y start_POSTSUPERSCRIPT italic_a italic_b italic_c end_POSTSUPERSCRIPT., and Cab:=GabκZabassignsubscript𝐶𝑎𝑏subscript𝐺𝑎𝑏𝜅subscript𝑍𝑎𝑏C_{ab}:=G_{ab}-\kappa Z_{ab}italic_C start_POSTSUBSCRIPT italic_a italic_b end_POSTSUBSCRIPT := italic_G start_POSTSUBSCRIPT italic_a italic_b end_POSTSUBSCRIPT - italic_κ italic_Z start_POSTSUBSCRIPT italic_a italic_b end_POSTSUBSCRIPT, where Zabsubscript𝑍𝑎𝑏Z_{ab}italic_Z start_POSTSUBSCRIPT italic_a italic_b end_POSTSUBSCRIPT is a field that we require to be equal on shell to the energy-momentum tensor obtained from some matter action SMsubscript𝑆𝑀S_{M}italic_S start_POSTSUBSCRIPT italic_M end_POSTSUBSCRIPT. The variation of ΨabcdefsubscriptΨ𝑎𝑏𝑐𝑑𝑒𝑓\Psi_{abcdef}roman_Ψ start_POSTSUBSCRIPT italic_a italic_b italic_c italic_d italic_e italic_f end_POSTSUBSCRIPT yields Yabc=0superscript𝑌𝑎𝑏𝑐0Y^{abc}=0italic_Y start_POSTSUPERSCRIPT italic_a italic_b italic_c end_POSTSUPERSCRIPT = 0, and the variation with respect to Yabcsuperscript𝑌𝑎𝑏𝑐Y^{abc}italic_Y start_POSTSUPERSCRIPT italic_a italic_b italic_c end_POSTSUPERSCRIPT yields:

[aCb]c+2ΨabcdefYdef=0,\nabla_{[a}C_{b]c}+2\Psi_{abcdef}Y^{def}=0,∇ start_POSTSUBSCRIPT [ italic_a end_POSTSUBSCRIPT italic_C start_POSTSUBSCRIPT italic_b ] italic_c end_POSTSUBSCRIPT + 2 roman_Ψ start_POSTSUBSCRIPT italic_a italic_b italic_c italic_d italic_e italic_f end_POSTSUBSCRIPT italic_Y start_POSTSUPERSCRIPT italic_d italic_e italic_f end_POSTSUPERSCRIPT = 0 , (26)

which reduces to the Codazzi equation when Yabc=0superscript𝑌𝑎𝑏𝑐0Y^{abc}=0italic_Y start_POSTSUPERSCRIPT italic_a italic_b italic_c end_POSTSUPERSCRIPT = 0. The variation of the action SCsubscript𝑆𝐶S_{C}italic_S start_POSTSUBSCRIPT italic_C end_POSTSUBSCRIPT with respect to gabsuperscript𝑔𝑎𝑏g^{ab}italic_g start_POSTSUPERSCRIPT italic_a italic_b end_POSTSUPERSCRIPT and Zabsubscript𝑍𝑎𝑏Z_{ab}italic_Z start_POSTSUBSCRIPT italic_a italic_b end_POSTSUBSCRIPT is rather nontrivial, but ultimately consists of terms proportional to products of Yabcsuperscript𝑌𝑎𝑏𝑐Y^{abc}italic_Y start_POSTSUPERSCRIPT italic_a italic_b italic_c end_POSTSUPERSCRIPT and its derivatives; variations with respect to gabsuperscript𝑔𝑎𝑏g^{ab}italic_g start_POSTSUPERSCRIPT italic_a italic_b end_POSTSUPERSCRIPT and Zabsubscript𝑍𝑎𝑏Z_{ab}italic_Z start_POSTSUBSCRIPT italic_a italic_b end_POSTSUBSCRIPT vanish when Yabc=0superscript𝑌𝑎𝑏𝑐0Y^{abc}=0italic_Y start_POSTSUPERSCRIPT italic_a italic_b italic_c end_POSTSUPERSCRIPT = 0. To independently enforce the divergence-free constraint, one may add to the action the following:

SD:=assignsubscript𝑆𝐷absent\displaystyle S_{D}:=italic_S start_POSTSUBSCRIPT italic_D end_POSTSUBSCRIPT := Ω|g|d4x[UbaCab+ψabUaUb],subscriptΩ𝑔superscript𝑑4𝑥delimited-[]superscript𝑈𝑏superscript𝑎subscript𝐶𝑎𝑏subscript𝜓𝑎𝑏superscript𝑈𝑎superscript𝑈𝑏\displaystyle\int_{\Omega}\sqrt{|g|}d^{4}x[U^{b}\nabla^{a}C_{ab}+\psi_{ab}U^{a% }U^{b}],∫ start_POSTSUBSCRIPT roman_Ω end_POSTSUBSCRIPT square-root start_ARG | italic_g | end_ARG italic_d start_POSTSUPERSCRIPT 4 end_POSTSUPERSCRIPT italic_x [ italic_U start_POSTSUPERSCRIPT italic_b end_POSTSUPERSCRIPT ∇ start_POSTSUPERSCRIPT italic_a end_POSTSUPERSCRIPT italic_C start_POSTSUBSCRIPT italic_a italic_b end_POSTSUBSCRIPT + italic_ψ start_POSTSUBSCRIPT italic_a italic_b end_POSTSUBSCRIPT italic_U start_POSTSUPERSCRIPT italic_a end_POSTSUPERSCRIPT italic_U start_POSTSUPERSCRIPT italic_b end_POSTSUPERSCRIPT ] , (27)

where Uasuperscript𝑈𝑎U^{a}italic_U start_POSTSUPERSCRIPT italic_a end_POSTSUPERSCRIPT is an auxiliary field and ψab=ψbasubscript𝜓𝑎𝑏subscript𝜓𝑏𝑎\psi_{ab}=\psi_{ba}italic_ψ start_POSTSUBSCRIPT italic_a italic_b end_POSTSUBSCRIPT = italic_ψ start_POSTSUBSCRIPT italic_b italic_a end_POSTSUBSCRIPT is a Lagrange multiplier. The variation of SDsubscript𝑆𝐷S_{D}italic_S start_POSTSUBSCRIPT italic_D end_POSTSUBSCRIPT with respect to ψabsubscript𝜓𝑎𝑏\psi_{ab}italic_ψ start_POSTSUBSCRIPT italic_a italic_b end_POSTSUBSCRIPT implies Ua=0superscript𝑈𝑎0U^{a}=0italic_U start_POSTSUPERSCRIPT italic_a end_POSTSUPERSCRIPT = 0, and the variation with respect to Uasuperscript𝑈𝑎U^{a}italic_U start_POSTSUPERSCRIPT italic_a end_POSTSUPERSCRIPT yields:

aCab+2ψabUa=0.superscript𝑎subscript𝐶𝑎𝑏2subscript𝜓𝑎𝑏superscript𝑈𝑎0\nabla^{a}C_{ab}+2\psi_{ab}U^{a}=0.∇ start_POSTSUPERSCRIPT italic_a end_POSTSUPERSCRIPT italic_C start_POSTSUBSCRIPT italic_a italic_b end_POSTSUBSCRIPT + 2 italic_ψ start_POSTSUBSCRIPT italic_a italic_b end_POSTSUBSCRIPT italic_U start_POSTSUPERSCRIPT italic_a end_POSTSUPERSCRIPT = 0 . (28)

When Ua=0superscript𝑈𝑎0U^{a}=0italic_U start_POSTSUPERSCRIPT italic_a end_POSTSUPERSCRIPT = 0, the above reduces to the divergence-free condition, and the variation of SDsubscript𝑆𝐷S_{D}italic_S start_POSTSUBSCRIPT italic_D end_POSTSUBSCRIPT with respect to the remaining variables vanishes.

It is straightforward to apply this strategy to the other theories considered in this article. For conformal Killing gravity, we construct the action:

SK:=assignsubscript𝑆𝐾absent\displaystyle S_{K}:=italic_S start_POSTSUBSCRIPT italic_K end_POSTSUBSCRIPT := Ω|g|d4x[Yabc(aK¯bc)+ΨabcdefYabcYdef].\displaystyle\int_{\Omega}\sqrt{|g|}d^{4}x[Y^{abc}\nabla_{(a}\bar{K}_{bc)}+% \Psi_{abcdef}Y^{abc}Y^{def}].∫ start_POSTSUBSCRIPT roman_Ω end_POSTSUBSCRIPT square-root start_ARG | italic_g | end_ARG italic_d start_POSTSUPERSCRIPT 4 end_POSTSUPERSCRIPT italic_x [ italic_Y start_POSTSUPERSCRIPT italic_a italic_b italic_c end_POSTSUPERSCRIPT ∇ start_POSTSUBSCRIPT ( italic_a end_POSTSUBSCRIPT over¯ start_ARG italic_K end_ARG start_POSTSUBSCRIPT italic_b italic_c ) end_POSTSUBSCRIPT + roman_Ψ start_POSTSUBSCRIPT italic_a italic_b italic_c italic_d italic_e italic_f end_POSTSUBSCRIPT italic_Y start_POSTSUPERSCRIPT italic_a italic_b italic_c end_POSTSUPERSCRIPT italic_Y start_POSTSUPERSCRIPT italic_d italic_e italic_f end_POSTSUPERSCRIPT ] . (29)

where K¯ab:=6KabgabKassignsubscript¯𝐾𝑎𝑏6subscript𝐾𝑎𝑏subscript𝑔𝑎𝑏𝐾\bar{K}_{ab}:=6K_{ab}-g_{ab}Kover¯ start_ARG italic_K end_ARG start_POSTSUBSCRIPT italic_a italic_b end_POSTSUBSCRIPT := 6 italic_K start_POSTSUBSCRIPT italic_a italic_b end_POSTSUBSCRIPT - italic_g start_POSTSUBSCRIPT italic_a italic_b end_POSTSUBSCRIPT italic_K. For the trace-free Einstein equations, we have:

STF:=assignsubscript𝑆𝑇𝐹absent\displaystyle S_{TF}:=italic_S start_POSTSUBSCRIPT italic_T italic_F end_POSTSUBSCRIPT := Ω|g|d4x[Yab(R¯abκZ¯ab)+ΨabcdYabYcd],subscriptΩ𝑔superscript𝑑4𝑥delimited-[]superscript𝑌𝑎𝑏subscript¯𝑅𝑎𝑏𝜅subscript¯𝑍𝑎𝑏subscriptΨ𝑎𝑏𝑐𝑑superscript𝑌𝑎𝑏superscript𝑌𝑐𝑑\displaystyle\int_{\Omega}\sqrt{|g|}d^{4}x[Y^{ab}\left(\underline{R}_{ab}-% \kappa\underline{Z}_{ab}\right)+\Psi_{abcd}Y^{ab}Y^{cd}],∫ start_POSTSUBSCRIPT roman_Ω end_POSTSUBSCRIPT square-root start_ARG | italic_g | end_ARG italic_d start_POSTSUPERSCRIPT 4 end_POSTSUPERSCRIPT italic_x [ italic_Y start_POSTSUPERSCRIPT italic_a italic_b end_POSTSUPERSCRIPT ( under¯ start_ARG italic_R end_ARG start_POSTSUBSCRIPT italic_a italic_b end_POSTSUBSCRIPT - italic_κ under¯ start_ARG italic_Z end_ARG start_POSTSUBSCRIPT italic_a italic_b end_POSTSUBSCRIPT ) + roman_Ψ start_POSTSUBSCRIPT italic_a italic_b italic_c italic_d end_POSTSUBSCRIPT italic_Y start_POSTSUPERSCRIPT italic_a italic_b end_POSTSUPERSCRIPT italic_Y start_POSTSUPERSCRIPT italic_c italic_d end_POSTSUPERSCRIPT ] , (30)

where R¯ab:=Rab14Rgabassignsubscript¯𝑅𝑎𝑏subscript𝑅𝑎𝑏14𝑅subscript𝑔𝑎𝑏\underline{R}_{ab}:=R_{ab}-\tfrac{1}{4}Rg_{ab}under¯ start_ARG italic_R end_ARG start_POSTSUBSCRIPT italic_a italic_b end_POSTSUBSCRIPT := italic_R start_POSTSUBSCRIPT italic_a italic_b end_POSTSUBSCRIPT - divide start_ARG 1 end_ARG start_ARG 4 end_ARG italic_R italic_g start_POSTSUBSCRIPT italic_a italic_b end_POSTSUBSCRIPT and Z¯ab:=Rab14Zgabassignsubscript¯𝑍𝑎𝑏subscript𝑅𝑎𝑏14𝑍subscript𝑔𝑎𝑏\underline{Z}_{ab}:=R_{ab}-\tfrac{1}{4}Zg_{ab}under¯ start_ARG italic_Z end_ARG start_POSTSUBSCRIPT italic_a italic_b end_POSTSUBSCRIPT := italic_R start_POSTSUBSCRIPT italic_a italic_b end_POSTSUBSCRIPT - divide start_ARG 1 end_ARG start_ARG 4 end_ARG italic_Z italic_g start_POSTSUBSCRIPT italic_a italic_b end_POSTSUBSCRIPT. We emphasize that no constraint on the metric determinant is required.

The matter sector is a somewhat delicate matter. Naively, one can add a term to the matter action:

SM=SM+Ω|g|d4x(Z~abgab),subscriptsuperscript𝑆𝑀subscript𝑆𝑀subscriptΩ𝑔superscript𝑑4𝑥subscript~𝑍𝑎𝑏superscript𝑔𝑎𝑏S^{\prime}_{M}=S_{M}+\int_{\Omega}\sqrt{|g|}d^{4}x(\tilde{Z}_{ab}g^{ab}),italic_S start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_M end_POSTSUBSCRIPT = italic_S start_POSTSUBSCRIPT italic_M end_POSTSUBSCRIPT + ∫ start_POSTSUBSCRIPT roman_Ω end_POSTSUBSCRIPT square-root start_ARG | italic_g | end_ARG italic_d start_POSTSUPERSCRIPT 4 end_POSTSUPERSCRIPT italic_x ( over~ start_ARG italic_Z end_ARG start_POSTSUBSCRIPT italic_a italic_b end_POSTSUBSCRIPT italic_g start_POSTSUPERSCRIPT italic_a italic_b end_POSTSUPERSCRIPT ) , (31)

so that the variation with respect to gabsuperscript𝑔𝑎𝑏g^{ab}italic_g start_POSTSUPERSCRIPT italic_a italic_b end_POSTSUPERSCRIPT yields Z~ab=Tab12Tgabsubscript~𝑍𝑎𝑏subscript𝑇𝑎𝑏12𝑇subscript𝑔𝑎𝑏\tilde{Z}_{ab}=T_{ab}-\tfrac{1}{2}Tg_{ab}over~ start_ARG italic_Z end_ARG start_POSTSUBSCRIPT italic_a italic_b end_POSTSUBSCRIPT = italic_T start_POSTSUBSCRIPT italic_a italic_b end_POSTSUBSCRIPT - divide start_ARG 1 end_ARG start_ARG 2 end_ARG italic_T italic_g start_POSTSUBSCRIPT italic_a italic_b end_POSTSUBSCRIPT. However, the difficulty with such a strategy is that the variation δZ~ab𝛿subscript~𝑍𝑎𝑏\delta\tilde{Z}_{ab}italic_δ over~ start_ARG italic_Z end_ARG start_POSTSUBSCRIPT italic_a italic_b end_POSTSUBSCRIPT will place a nontrivial constraint on gabsuperscript𝑔𝑎𝑏g^{ab}italic_g start_POSTSUPERSCRIPT italic_a italic_b end_POSTSUPERSCRIPT unless it becomes a surface term δZ~ab=cϕabc𝛿subscript~𝑍𝑎𝑏subscript𝑐subscriptsuperscriptitalic-ϕ𝑐𝑎𝑏\delta\tilde{Z}_{ab}=\nabla_{c}\phi^{c}_{ab}italic_δ over~ start_ARG italic_Z end_ARG start_POSTSUBSCRIPT italic_a italic_b end_POSTSUBSCRIPT = ∇ start_POSTSUBSCRIPT italic_c end_POSTSUBSCRIPT italic_ϕ start_POSTSUPERSCRIPT italic_c end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_a italic_b end_POSTSUBSCRIPT, or otherwise trivializes. However, the only known tensor with this property is a Ricci tensor formed from a metric-compatible connection, in which case we recover (by way of a Palatini variation) the Einstein equations with a specified CC, or the corresponding Einstein-Cartan equations if fermions are included.

Consider the action:

SM¯=Ω|g|d4x[(gabhab)Zab+(LhLg)+Ξ+JζJ],subscript𝑆¯𝑀subscriptΩ𝑔superscript𝑑4𝑥delimited-[]superscript𝑔𝑎𝑏superscript𝑎𝑏subscript𝑍𝑎𝑏subscript𝐿subscript𝐿𝑔Ξ𝐽𝜁𝐽S_{\bar{M}}=\int_{\Omega}\sqrt{|g|}d^{4}x\biggl{[}(g^{ab}-h^{ab})Z_{ab}+(L_{h}% -L_{g})\mathcal{E}+\Xi+J\cdot\zeta\cdot J\biggr{]},italic_S start_POSTSUBSCRIPT over¯ start_ARG italic_M end_ARG end_POSTSUBSCRIPT = ∫ start_POSTSUBSCRIPT roman_Ω end_POSTSUBSCRIPT square-root start_ARG | italic_g | end_ARG italic_d start_POSTSUPERSCRIPT 4 end_POSTSUPERSCRIPT italic_x [ ( italic_g start_POSTSUPERSCRIPT italic_a italic_b end_POSTSUPERSCRIPT - italic_h start_POSTSUPERSCRIPT italic_a italic_b end_POSTSUPERSCRIPT ) italic_Z start_POSTSUBSCRIPT italic_a italic_b end_POSTSUBSCRIPT + ( italic_L start_POSTSUBSCRIPT italic_h end_POSTSUBSCRIPT - italic_L start_POSTSUBSCRIPT italic_g end_POSTSUBSCRIPT ) caligraphic_E + roman_Ξ + italic_J ⋅ italic_ζ ⋅ italic_J ] , (32)

where ϕitalic-ϕ\phiitalic_ϕ is a matter field (assumed to be nonspinorial and minimally coupled), χ𝜒\chiitalic_χ is an extra field that is equal to ϕitalic-ϕ\phiitalic_ϕ on shell, habsuperscript𝑎𝑏h^{ab}italic_h start_POSTSUPERSCRIPT italic_a italic_b end_POSTSUPERSCRIPT and J𝐽Jitalic_J are auxiliary fields, ζ𝜁\zetaitalic_ζ is a Lagrange multiplier333In particular, one assumes ζ𝜁\zetaitalic_ζ to be defined so that JζJ𝐽𝜁𝐽J\cdot\zeta\cdot Jitalic_J ⋅ italic_ζ ⋅ italic_J is a scalar., and the following have been defined:

Lh:=L(χ,aχ,hab),Lg:=L(φ,aφ,gab),formulae-sequenceassignsubscript𝐿𝐿𝜒subscript𝑎𝜒superscript𝑎𝑏assignsubscript𝐿𝑔𝐿𝜑subscript𝑎𝜑superscript𝑔𝑎𝑏\displaystyle L_{h}:=L(\chi,\partial_{a}\chi,h^{ab}),\qquad L_{g}:=L(\varphi,% \partial_{a}\varphi,g^{ab}),italic_L start_POSTSUBSCRIPT italic_h end_POSTSUBSCRIPT := italic_L ( italic_χ , ∂ start_POSTSUBSCRIPT italic_a end_POSTSUBSCRIPT italic_χ , italic_h start_POSTSUPERSCRIPT italic_a italic_b end_POSTSUPERSCRIPT ) , italic_L start_POSTSUBSCRIPT italic_g end_POSTSUBSCRIPT := italic_L ( italic_φ , ∂ start_POSTSUBSCRIPT italic_a end_POSTSUBSCRIPT italic_φ , italic_g start_POSTSUPERSCRIPT italic_a italic_b end_POSTSUPERSCRIPT ) , (33)
Ξ:=(φχ)J,(gab,h¯ab):=exp[gabh¯ab2].formulae-sequenceassignΞ𝜑𝜒𝐽assignsuperscript𝑔𝑎𝑏subscript¯𝑎𝑏superscript𝑔𝑎𝑏subscript¯𝑎𝑏2\displaystyle\Xi:=(\varphi-\chi)\cdot J,\qquad\mathcal{E}(g^{ab},\underline{h}% _{ab}):=\exp\left[\sqrt{g^{ab}\underline{h}_{ab}}-2\right].roman_Ξ := ( italic_φ - italic_χ ) ⋅ italic_J , caligraphic_E ( italic_g start_POSTSUPERSCRIPT italic_a italic_b end_POSTSUPERSCRIPT , under¯ start_ARG italic_h end_ARG start_POSTSUBSCRIPT italic_a italic_b end_POSTSUBSCRIPT ) := roman_exp [ square-root start_ARG italic_g start_POSTSUPERSCRIPT italic_a italic_b end_POSTSUPERSCRIPT under¯ start_ARG italic_h end_ARG start_POSTSUBSCRIPT italic_a italic_b end_POSTSUBSCRIPT end_ARG - 2 ] .

The variation with respect to Zabsubscript𝑍𝑎𝑏Z_{ab}italic_Z start_POSTSUBSCRIPT italic_a italic_b end_POSTSUBSCRIPT yields hab=gabsuperscript𝑎𝑏superscript𝑔𝑎𝑏h^{ab}=g^{ab}italic_h start_POSTSUPERSCRIPT italic_a italic_b end_POSTSUPERSCRIPT = italic_g start_POSTSUPERSCRIPT italic_a italic_b end_POSTSUPERSCRIPT, the variation with respect to ζ𝜁\zetaitalic_ζ yields J=0𝐽0J=0italic_J = 0, and the variation with respect to J𝐽Jitalic_J yields φ=χ𝜑𝜒\varphi=\chiitalic_φ = italic_χ when evaluated on J=0𝐽0J=0italic_J = 0. The variation of \mathcal{E}caligraphic_E has the form:

δ=exp{gabh¯ab2}gabh¯ab[δgabh¯abgcdh¯ach¯bdδhab],𝛿superscript𝑔𝑎𝑏subscript¯𝑎𝑏2superscript𝑔𝑎𝑏subscript¯𝑎𝑏delimited-[]𝛿superscript𝑔𝑎𝑏subscript¯𝑎𝑏superscript𝑔𝑐𝑑subscript¯𝑎𝑐subscript¯𝑏𝑑𝛿superscript𝑎𝑏\delta\mathcal{E}=\frac{\exp\left\{\sqrt{g^{ab}\underline{h}_{ab}}-2\right\}}{% \sqrt{g^{ab}\underline{h}_{ab}}}\biggl{[}\delta g^{ab}\underline{h}_{ab}-g^{cd% }\underline{h}_{ac}\underline{h}_{bd}\delta h^{ab}\biggr{]},italic_δ caligraphic_E = divide start_ARG roman_exp { square-root start_ARG italic_g start_POSTSUPERSCRIPT italic_a italic_b end_POSTSUPERSCRIPT under¯ start_ARG italic_h end_ARG start_POSTSUBSCRIPT italic_a italic_b end_POSTSUBSCRIPT end_ARG - 2 } end_ARG start_ARG square-root start_ARG italic_g start_POSTSUPERSCRIPT italic_a italic_b end_POSTSUPERSCRIPT under¯ start_ARG italic_h end_ARG start_POSTSUBSCRIPT italic_a italic_b end_POSTSUBSCRIPT end_ARG end_ARG [ italic_δ italic_g start_POSTSUPERSCRIPT italic_a italic_b end_POSTSUPERSCRIPT under¯ start_ARG italic_h end_ARG start_POSTSUBSCRIPT italic_a italic_b end_POSTSUBSCRIPT - italic_g start_POSTSUPERSCRIPT italic_c italic_d end_POSTSUPERSCRIPT under¯ start_ARG italic_h end_ARG start_POSTSUBSCRIPT italic_a italic_c end_POSTSUBSCRIPT under¯ start_ARG italic_h end_ARG start_POSTSUBSCRIPT italic_b italic_d end_POSTSUBSCRIPT italic_δ italic_h start_POSTSUPERSCRIPT italic_a italic_b end_POSTSUPERSCRIPT ] , (34)

which on hab=gabsuperscript𝑎𝑏superscript𝑔𝑎𝑏h^{ab}=g^{ab}italic_h start_POSTSUPERSCRIPT italic_a italic_b end_POSTSUPERSCRIPT = italic_g start_POSTSUPERSCRIPT italic_a italic_b end_POSTSUPERSCRIPT simplifies to:

δ=12gab(δgabδhab),𝛿12subscript𝑔𝑎𝑏𝛿superscript𝑔𝑎𝑏𝛿superscript𝑎𝑏\delta\mathcal{E}=\frac{1}{2}g_{ab}\left(\delta g^{ab}-\delta h^{ab}\right),italic_δ caligraphic_E = divide start_ARG 1 end_ARG start_ARG 2 end_ARG italic_g start_POSTSUBSCRIPT italic_a italic_b end_POSTSUBSCRIPT ( italic_δ italic_g start_POSTSUPERSCRIPT italic_a italic_b end_POSTSUPERSCRIPT - italic_δ italic_h start_POSTSUPERSCRIPT italic_a italic_b end_POSTSUPERSCRIPT ) , (35)

Similarly, the Lagrangian variations satisfy:

δLh𝛿subscript𝐿\displaystyle\delta L_{h}italic_δ italic_L start_POSTSUBSCRIPT italic_h end_POSTSUBSCRIPT =Lhhabδhab+Lhχδχ+Lh(aχ)δaχ,absentsubscript𝐿superscript𝑎𝑏𝛿superscript𝑎𝑏subscript𝐿𝜒𝛿𝜒subscript𝐿subscript𝑎𝜒𝛿subscript𝑎𝜒\displaystyle=\frac{\partial L_{h}}{\partial h^{ab}}\delta h^{ab}+\frac{% \partial L_{h}}{\partial\chi}\delta\chi+\frac{\partial L_{h}}{\partial(% \partial_{a}\chi)}\delta\partial_{a}\chi,= divide start_ARG ∂ italic_L start_POSTSUBSCRIPT italic_h end_POSTSUBSCRIPT end_ARG start_ARG ∂ italic_h start_POSTSUPERSCRIPT italic_a italic_b end_POSTSUPERSCRIPT end_ARG italic_δ italic_h start_POSTSUPERSCRIPT italic_a italic_b end_POSTSUPERSCRIPT + divide start_ARG ∂ italic_L start_POSTSUBSCRIPT italic_h end_POSTSUBSCRIPT end_ARG start_ARG ∂ italic_χ end_ARG italic_δ italic_χ + divide start_ARG ∂ italic_L start_POSTSUBSCRIPT italic_h end_POSTSUBSCRIPT end_ARG start_ARG ∂ ( ∂ start_POSTSUBSCRIPT italic_a end_POSTSUBSCRIPT italic_χ ) end_ARG italic_δ ∂ start_POSTSUBSCRIPT italic_a end_POSTSUBSCRIPT italic_χ , (36)
δLg𝛿subscript𝐿𝑔\displaystyle\delta L_{g}italic_δ italic_L start_POSTSUBSCRIPT italic_g end_POSTSUBSCRIPT =Lggabδgab+Lgφδφ+Lg(aφ)δaφ.absentsubscript𝐿𝑔superscript𝑔𝑎𝑏𝛿superscript𝑔𝑎𝑏subscript𝐿𝑔𝜑𝛿𝜑subscript𝐿𝑔subscript𝑎𝜑𝛿subscript𝑎𝜑\displaystyle=\frac{\partial L_{g}}{\partial g^{ab}}\delta g^{ab}+\frac{% \partial L_{g}}{\partial\varphi}\delta\varphi+\frac{\partial L_{g}}{\partial(% \partial_{a}\varphi)}\delta\partial_{a}\varphi.= divide start_ARG ∂ italic_L start_POSTSUBSCRIPT italic_g end_POSTSUBSCRIPT end_ARG start_ARG ∂ italic_g start_POSTSUPERSCRIPT italic_a italic_b end_POSTSUPERSCRIPT end_ARG italic_δ italic_g start_POSTSUPERSCRIPT italic_a italic_b end_POSTSUPERSCRIPT + divide start_ARG ∂ italic_L start_POSTSUBSCRIPT italic_g end_POSTSUBSCRIPT end_ARG start_ARG ∂ italic_φ end_ARG italic_δ italic_φ + divide start_ARG ∂ italic_L start_POSTSUBSCRIPT italic_g end_POSTSUBSCRIPT end_ARG start_ARG ∂ ( ∂ start_POSTSUBSCRIPT italic_a end_POSTSUBSCRIPT italic_φ ) end_ARG italic_δ ∂ start_POSTSUBSCRIPT italic_a end_POSTSUBSCRIPT italic_φ .

One finds that the variations with respect to gabsuperscript𝑔𝑎𝑏g^{ab}italic_g start_POSTSUPERSCRIPT italic_a italic_b end_POSTSUPERSCRIPT and habsuperscript𝑎𝑏h^{ab}italic_h start_POSTSUPERSCRIPT italic_a italic_b end_POSTSUPERSCRIPT yield two constraints that coincide when evaluating on gab=habsuperscript𝑔𝑎𝑏superscript𝑎𝑏g^{ab}=h^{ab}italic_g start_POSTSUPERSCRIPT italic_a italic_b end_POSTSUPERSCRIPT = italic_h start_POSTSUPERSCRIPT italic_a italic_b end_POSTSUPERSCRIPT and φ=χ𝜑𝜒\varphi=\chiitalic_φ = italic_χ:

Zab=Lhhab12habLh=Lggab12gabLg,subscript𝑍𝑎𝑏subscript𝐿superscript𝑎𝑏12subscript𝑎𝑏subscript𝐿subscript𝐿𝑔superscript𝑔𝑎𝑏12subscript𝑔𝑎𝑏subscript𝐿𝑔\displaystyle Z_{ab}=\frac{\partial L_{h}}{\partial h^{ab}}-\frac{1}{2}h_{ab}L% _{h}=\frac{\partial L_{g}}{\partial g^{ab}}-\frac{1}{2}g_{ab}L_{g},italic_Z start_POSTSUBSCRIPT italic_a italic_b end_POSTSUBSCRIPT = divide start_ARG ∂ italic_L start_POSTSUBSCRIPT italic_h end_POSTSUBSCRIPT end_ARG start_ARG ∂ italic_h start_POSTSUPERSCRIPT italic_a italic_b end_POSTSUPERSCRIPT end_ARG - divide start_ARG 1 end_ARG start_ARG 2 end_ARG italic_h start_POSTSUBSCRIPT italic_a italic_b end_POSTSUBSCRIPT italic_L start_POSTSUBSCRIPT italic_h end_POSTSUBSCRIPT = divide start_ARG ∂ italic_L start_POSTSUBSCRIPT italic_g end_POSTSUBSCRIPT end_ARG start_ARG ∂ italic_g start_POSTSUPERSCRIPT italic_a italic_b end_POSTSUPERSCRIPT end_ARG - divide start_ARG 1 end_ARG start_ARG 2 end_ARG italic_g start_POSTSUBSCRIPT italic_a italic_b end_POSTSUBSCRIPT italic_L start_POSTSUBSCRIPT italic_g end_POSTSUBSCRIPT , (37)

which makes use of the fact that for gab=habsuperscript𝑔𝑎𝑏superscript𝑎𝑏g^{ab}=h^{ab}italic_g start_POSTSUPERSCRIPT italic_a italic_b end_POSTSUPERSCRIPT = italic_h start_POSTSUPERSCRIPT italic_a italic_b end_POSTSUPERSCRIPT and φ=χ𝜑𝜒\varphi=\chiitalic_φ = italic_χ, the integrand of Eq. (32) is constructed to vanish. The variations with respect to φ𝜑\varphiitalic_φ and χ𝜒\chiitalic_χ yield matter field equations that coincide on gab=habsuperscript𝑔𝑎𝑏superscript𝑎𝑏g^{ab}=h^{ab}italic_g start_POSTSUPERSCRIPT italic_a italic_b end_POSTSUPERSCRIPT = italic_h start_POSTSUPERSCRIPT italic_a italic_b end_POSTSUPERSCRIPT and φ=χ𝜑𝜒\varphi=\chiitalic_φ = italic_χ, and reduce to the usual field equations obtained from Lgsubscript𝐿𝑔L_{g}italic_L start_POSTSUBSCRIPT italic_g end_POSTSUBSCRIPT after applying the constraint J=0𝐽0J=0italic_J = 0.

Discussion — In this article, we considered some theories in which the CC arises as an integration constant, including trace-free Einstein gravity, the Codazzi equation, and conformal Killing gravity. We revisited a recent debate [33, 34, 35, 36] concerning the issue of nonuniqueness of solutions to equations of Cotton gravity in symmetry-reduced settiongs. By taking additional derivatives of the Codazzi equations, we find that solutions of the Codazzi equations must satisfy at least one of two conditions: either the solutions coincide with that of the Einstein equations (with arbitrary CC), or that certain contractions of the curvature tensor, the Codazzi tensor, and their derivatives must vanish (in particular, all terms on the right-hand side of Eq. (24) must vanish). The latter can be satisfied in situations of high symmetry, even if the solutions do not satisfy the Einstein equations, but if we regard such non-Einstein solutions as physical, the right-hand side of Eq. (24) must vanish exactly even under perturbations. It is unlikely that such conditions can be satisfied unless the perturbations are highly restricted, and one may in this manner argue that such non-Einstein solutions are unphysical. We find that the Codazzi equations can in principle be closed with the addition of the divergence-free constraint so that such non-Einstein solutions can be excluded at the level of initial data.

We formulated a rather general variational principle that can yield the trace-free Einstein equations, the Codazzi equation, and the equations of conformal Killing gravity. An action for trace-free gravity without a unimodular constraint [38] was recently proposed (another possible approach is described in [23, *Tapia:1998qt, *Tapia:1998as]), but the matter sector was not explicitly incorporated in their treatment; our action (30) provides an alternative that can readily incorporate the matter sector via the action (32). The appearance of Lagrange multipliers and auxiliary fields in the actions corresponds to introducing additional nondynamical degrees of freedom—one might prefer an action that depends exclusively on dynamical degrees of freedom. Moreover, there is some question of whether auxiliary fields in gravity theories are pathological [39, 40, 41]. Regarding the latter, there may be instances in which some of the pathologies (particularly those associated with sharp density gradients) can be mitigated [42, 43, 44], and one can perhaps regard auxiliary fields as effective descriptions for underlying dynamical degrees of freedom (with large masses). Regarding the former concern, it is well-known that one can express the Lagrangian for f(R)𝑓𝑅f(R)italic_f ( italic_R ) theories in the O’Hanlon [45] form ψRV(ψ)𝜓𝑅𝑉𝜓\psi R-V(\psi)italic_ψ italic_R - italic_V ( italic_ψ ) (see also [46, 47, 48, 49, 50] and references therein), where ψ𝜓\psiitalic_ψ is an auxiliary field, and one can easily recover the f(R)𝑓𝑅f(R)italic_f ( italic_R ) Lagrangian from the O’Hanlon Lagrangian by algebraic elimination. Of course, a naive algebraic elimination would trivialize the actions considered here, but the point is that auxiliary fields can provide a starting point for constructing an action exclusively from dynamical degrees of freedom. To clarify, we are not suggesting that the variational principle we have proposed is fundamental—indeed, one can construct a variational principle for virtually any set of consistent field equations in this manner.444Presumably one could attempt to use the same strategy to construct a variational principle for inconsistent field equations, but this will result in an action that does not admit extrema. Instead, we argue that it can serve as a first step towards the construction of a more fundamental action for these theories.

Ghosts are a concern; for instance, the Codazzi equation contains up to three derivatives of the metric. One might also recognize that the duplicate matter degrees of freedom χ𝜒\chiitalic_χ in Eq. (32) are necessarily ghosts. At the classical level, there is no issue in the matter sector, as the matter field equations ultimately reduce to those identical to that obtained in the standard variational principle—there is no runaway behavior (see also [51]). In the gravitational sector, one can place restrictions on the initial data such that the system becomes dynamically equivalent to the Einstein equations with a CC. Regarding the Hamiltonian, we note that on shell, both the matter and gravitational actions vanish in our approach, and it is not too difficult to convince oneself that the total matter Hamiltonian must also vanish on shell. Though there is no issue at the classical level, the appearance of ghosts may be of particular concern for the quantization of the theory. Quantization is beyond the scope of the present article (see the discussion in [16] regarding ghosts in higher curvature theories), but we leave the reader with a couple of remarks. First, as argued in the preceding paragraph, the action presented in this article may require further modification at high energies, where the structure of the resulting action may differ significantly. Second, the vanishing of the Hamiltonian indicates that these theories suffer from the problem of time, in which the Hamiltonian operator does not yield the time evolution of the quantum state of the field, a well-known problem in canonical approaches to quantum gravity [52, *Isham:1992ms].

Acknowledgments — We thank Marco de Cesare, Naresh Dadhich, and Victor Tapia for helpful comments and references. PC and JCF acknowledge support from the Taiwan National Science and Technology Council (NSTC) under project numbers 110-2112-M-002-031 and 112-2811-M-002-132, and also from the Leung Center for Cosmology and Particle Astrophysics (LeCosPA), National Taiwan University.

References