Emergent Gravity from Vanishing Energy-Momentum Tensor

A constraint of vanishing energy-momentum tensor is motivated by a variety of perspectives on quantum gravity. We demonstrate in a concrete example how this constraint leads to a metric-independent theory in which quantum gravity emerges as a nonperturbative artifact of regularization-scale physics. We analyze a scalar theory similar to the Dirac-Born-Infeld (DBI) theory with vanishing gauge fields, with the DBI Lagrangian modulated by a scalar potential. In the limit of a large number of scalars, we explicitly demonstrate the existence of a composite massless spin-2 graviton in the spectrum that couples to matter as in Einstein gravity. We comment on the cosmological constant problem and the generalization to theories with fermions and gauge fields.


Introduction
The idea that gauge interactions and gravitation may arise as emergent phenomena has a long history. The earliest compelling model of emergent electromagnetism was due to Bjorken [1], who proposed a scenario in which a four-fermion interaction of the form L int = G (ψγ µ ψ)(ψγ µ ψ) gives rise to a massless spin-1 composite state that interacts like the photon in electrodynamics. Bjorken argued that the dynamics of electromagnetism emerges in this scenario if the electromagnetic current, J µ = ψγ µ ψ, develops a nonvanishing condensate in the vacuum. It was further argued by Eguchi that Nambu-Jona-Lasinio-type models with emergent gauge interactions, such as Bjorken's model, may be renormalizable despite naive dimensional power counting as a consequence of the interactions between fundamental fermions and the collective gauge field excitations [2].
During the development of the theory of the strong interactions, when substantial evidence for color confinement as a theoretical consequence of quantum chromodynamics was still lacking, the possibility was briefly considered that color confinement should be imposed via a constraint of vanishing color current on field configurations in an otherwise free theory of fermions [3,4]. The full dynamics of quantum chromodynamics appeared to emerge as a consequence of confinement so implemented, rather than the other way around.
Following the successful description of the electroweak and strong interactions by the Standard Model, a renormalizable gauge theory, the existence of emergent gauge interactions was evidently unnecessary in order to explain existing experiments and observations. Still, it remains a possibility that some or all of the Standard Model gauge interactions are absent in a more fundamental description apparent only at short distances, yet emerge

JHEP03(2017)134
in an effective description applicable to the low-energy environment in which experiments have so far been performed [5].
In the context of gravity, for which there still does not exist a fundamental description known to be consistent with the Standard Model, the paradigm that general relativity emerges as the effective description of a massless composite spin-2 state remains compelling and has motivated a large number of investigations. A review of some of the approaches to emergent gravity, including some of the ideas discussed here, was presented in ref. [6]. Much of the activity in this area has been inspired by Sakharov's observation that the assumption of general covariance in a quantum field theory, including its regulator, is generally sufficient to produce a curved-space effective action which contains the Einstein-Hilbert term [7]. The Einstein-Hilbert action is induced by radiative corrections even if it is absent at tree level, and the gravitational coupling is related to the regularization scale. Numerous explicit calculations in theories for which the spacetime metric is treated as an auxiliary field have demonstrated the validity of Sakharov's claim (see, for example, refs. [8][9][10][11]). However, in Sakharov's interpretation, Einstein's equations arise semiclassically, with the vacuum expectation value of the energy-momentum tensor acting as the source for an otherwise classical metric; the implications of Sakharov's observation for quantum gravity remain far from clear [12].
The observation that gauge interactions can emerge from a constraint of vanishing current suggests by analogy that gravitational interactions might emerge via a constraint of vanishing energy-momentum tensor. This is the possibility that we put forward in this paper, though our implementation of the constraint does not precisely parallel the related proposal for quantum chromodynamics proposed in ref. [3]. Vanishing of the energymomentum tensor is also motivated by a treatment of the spacetime metric as an auxiliary field, as in the Polyakov form of bosonic string theory [13]. Furthermore, a constraint of vanishing energy and momentum is in the spirit of the Wheeler-deWitt equation [14], which promotes the vanishing of energy and momentum densities on states, constraints which follow from time-and space-reparametrization invariance, to the status of an axiom underlying quantum gravity. However, vanishing of the full energy-momentum tensor is a stronger set of constraints than vanishing of the energy and momentum densities.
The main purpose of this paper is to analyze a diffeomorphism-invariant theory of N + D scalar fields which realizes the constraint of vanishing energy-momentum tensor and gives rise to an emergent gravitational interaction in D spacetime dimensions. At the linearized level, the composite graviton couples to the energy-momentum tensor of N of these scalars, while the remaining D scalars are gauged away. A similar model was discussed in ref. [15]. Our work differs from ref. [15], not only in the initial assumptions that lead to the general form of the action, but also in that we directly demonstrate the existence of a graviton pole in a scattering amplitude at the nonperturbative level. Avoiding the auxiliary field approach of ref. [15] allows us to address the issue of gauge fixing in a more transparent way and to directly identify the physical degrees of freedom that couple to the composite graviton in terms of the fields of the original theory.
Before proceeding, we should consider the feasibility of emergent gravitation in the first place. On the one hand, Weinberg demonstrated that a massless spin-2 state in a relativistic JHEP03(2017)134 field theory must couple to matter as in Einstein gravity [16], so our task is to demonstrate the existence of such a state in the spectrum and check that it gives rise to gravitational interactions consistent with Einstein gravity. However, Weinberg and Witten demonstrated the impossibility of massless spin-2 states in Lorentz-invariant field theories with an Smatrix and a nonvanishing conserved Lorentz-covariant energy-momentum tensor [17,18]. Any theory which claims to contain such massless spin-2 states must somehow violate the assumptions of the Weinberg-Witten theorem. The Weinberg-Witten theorem is proven by considering matrix elements of the conserved energy-momentum tensor between oneparticle states of definite helicity. By identifying the effect of a spatial rotation of the energy-momentum tensor on the one hand, and rotation of the states on the other, one finds that the two rotations can be equivalent only if the massless states have helicity ≤ 1. One possibility for evading the theorem is that the theory is diffeomorphism invariant from the outset, so that the energy-momentum tensor does not satisfy ∂ µ T µν = 0. This is how general relativity escapes the conclusions of the Weinberg-Witten theorem. The theory presented here is diffeomorphism invariant, but evades the Weinberg-Witten theorem even more directly: the energy-momentum tensor vanishes by construction.
To demonstrate the emergent gravitational interaction, we analyze the model at large N and find a nonperturbative massless spin-2 pole in the Fourier transform of the 2 → 2 scattering amplitude. While the classical theory is metric-independent and does not contain a nontrivial conserved energy-momentum tensor, after eliminating from functional integrals the integration over redundant degrees of freedom, we find that the composite graviton couples at the linearized level to the usual energy-momentum tensor of the surviving physical degrees of freedom. Up to regularization-scale-suppressed corrections, the effective lowenergy theory contains scalar fields with a potential, coupled to Einstein gravity.
We present the theory in section 2. We demonstrate the existence of a nonperturbative composite graviton for large N in section 3. We discuss nonlinear matter-gravity couplings in section 4. We comment on the cosmological constant problem (with more details in an appendix) and the possible generalization to a realistic theory of fermions and gauge interactions in section 5.

The theory
We begin by considering the action for a collection of N + D scalar fields in a curved D-dimensional spacetime described by a metric g µν , where µ, ν ∈ {0, 1, . . . , D − 1}: where g = det(g µν ), g µν is the inverse of the metric g µν , and η IJ are constants which we take to have the values of the Minkowski metric in D dimensions. (We use the mostly-minus convention for the signature of η IJ and the metric g µν .) The fields X 1 , X 2 , . . . , X D−1 have "wrong sign" kinetic terms, but due to general coordinate invariance those fields may be gauge fixed and are not independently dynamical, as we will discuss below. We assume that

JHEP03(2017)134
the fields X I do not appear in the potential V (φ a ) so that the action maintains a global shift symmetry X I (x) → X I (x) + ∆ I . The action (2.1) is also reparametrization invariant with the fields X I , φ a transforming as scalars and the background field g µν transforming as a metric tensor. We want to emphasize at this point that classically the background metric g µν has no dynamics. We proceed to define the theory by functional integral quantization over the fields φ a (x), X I (x) and g µν (x), subject to the constraint that the energy-momentum tensor vanishes: T µν (x) = 0. The partition function for the theory is (2. 2) The energy-momentum tensor is given by where the Lagrangian L is defined by the action in eq.
The constraint T µν = 0 determines the spacetime metric in terms of the scalar field configuration. Each field configuration contributes to the functional integral as if it propagates in a unique spacetime unrelated to Einstein's equations. An explicit solution to T µν (x) = 0 for the metric is For generic field configurations and potentials, g µν in eq. (2.5) is nonsingular. In the perturbative expansion about a Minkowski-space background that we employ later, difficulties associated with singularities in the metric do not appear. 1 With the spacetime metric determined by eq. (2.5), the action for the theory resembles the Dirac-Born-Infeld action with vanishing gauge field, modulated by the scalar-field potential function V (φ a ): (2.6) 1 A similar issue regarding the singular nature of the induced metric appears in string theory. Though most of the time the D-brane action is taked to be non-singular, Gibbons and Ishibashi [19] have shown that there exist classical D-brane configurations with vanishing gauge fields (specifically D3 branes embedded in a 5-dimensional Lorentzian spacetime) whose induced metric exhibits signature change. The authors further noted that despite the curvature singularities exhibited by the induced metric, from the point of the bulk the D-brane geometry is smooth everywhere. Here we will ignore such singular solutions.

JHEP03(2017)134
The action is now independent of the spacetime metric, by construction. As a consequence, the energy-momentum tensor defined as in eq. (2.3) vanishes identically. Equivalently, one can confirm that the Noether currents associated with space-and time-translation invariance vanish. Note that eq. (2.6) is diffeomorphism invariant. If one were to add other diffeomorphism-invariant terms constructed out of the composite metric operator g µν by hand, doing so would modify the relationship between the Planck scale and the parameters of the theory, discussed later, but generically would not affect physics at long distances otherwise. The functional integrals that determine the partition function and correlation functions in this theory include the redundant integration over field configurations related by coordinate transformations. In order to compute correlation functions in the theory, we gauge fix the spacetime parametrization by identifying the fields X I with the corresponding spacetime coordinates, in analogy with the static gauge condition in string theory, up to an overall constant factor: Here, V 0 is a dimensionful parameter that will be chosen later. Note that the gauge condition eq. (2.7) would not be possible if there were fewer than D scalar fields. With precisely D scalar fields, in the static gauge the composite metric is gauge-fixed to the Minkowski metric, and the theory does not have dynamical content. 2 In the functional integral we can include the identity in the form, ). Note that the Fadeev-Popov determinant is trivial here, and there are no ghosts associated with the static gauge condition: where the last identity in eq. (2.9) is a consequence of the gauge-fixing condition. The integrals over coordinate-transformation functions α µ (x) then factor out of correlation functions of coordinate-invariant observables. Locally we can always transform to spacetime coordinates x µ which satisfy the static gauge condition (2.7). However, globally this choice imposes a topological constraint on the space of allowed field configurations, a constraint 2 In fact, any Lorentz-invariant action can be written in a generally covariant form by replacing the Minkowski metric with the composite metric of eq. (2.5) including only the D fields X I . The consequences of general covariance, including the emergence of a massless spin-two state, depends on the existence of additional degrees of freedom in the composite metric. By analogy, a U(1) global symmetry can be given a gauge-invariant description by minimal coupling to the derivative of a scalar field ∂µϕ(x) rather than to a vector field A µ . Without the additional degrees of freedom in A µ , ∂µϕ(x) may be gauge-fixed to zero, and the dynamical content of the theory is not altered.

JHEP03(2017)134
which we will impose in our perturbative analysis, but which should not be necessary for the emergence of gravitation in the theory. We also note that in static gauge the field X 0 inherits the role of an internal clock [20], and the fields X i , i = 1, . . . D − 1 rulers. The clock and ruler fields provide physical meaning to the spacetime backdrop in which dynamics takes place. It is sensible to consider correlation functions of operators built out of the local fields φ a (x µ ) after gauge fixing, even though those correlation functions take a different functional form in other gauge choices. 3 The classical equations of motion for the fields φ a and X I following from the action eq. (2.6) may be written, where the composite operator g µν that plays the role of the spacetime metric is precisely the solution to T µν = 0 in eq. (2.5). As a result of the assumed shift symmetry which acts on the fields X I , the static gauge choice eq. (2.7) satisfies the equations of motion, eqs. (2.10) and (2.11), with fields φ a uniform and fixed to the minimum of the potential V (φ a ), and with g µν a constant tensor proportional to η µν . This suggests that an expansion in the fields φ a may remain perturbative in this gauge choice. We will find this to be the case, up to the expected instability of the Minkowski-space solution if a regularization-scale-dependent cosmological constant is not tuned to zero.

Graviton pole
In order to analyze the theory perturbatively, we write V (φ) = V 0 + ∆V (φ a ) and expand the action eq. (2.6) in powers of 1/V 0 . We also assume that N , the number of fields φ a in the theory, is large, and keep only leading terms in a 1/N expansion. For the moment we identify the constant V 0 in the potential with the V 0 that appeared in the gauge-fixing condition eq. (2.7). We modify this condition later. We write the gauge-fixed action as follows: and

JHEP03(2017)134
in the static gauge eq. (2.7). Expanding the determinant using det M = exp(Tr ln M ), we have the expansion where index contractions are done with the Minkowski metric η µν . Expanding the factor involving ∆V /V 0 , and using eq. (3.2), we find that the action can be written as For now, we assume a free theory with O(N )-symmetric potential in which case the first line of eq. (3.5) is the free part of the action. In terms of the energy-momentum tensor for free fields φ a in Minkowski space, 4 the interacting terms to O(1/V 0 ) in the Lagrangian take the simple form, We will now demonstrate that this theory includes a massless spin-two graviton state that mediates a gravitational interaction between φ a particles. We consider two-into-two scattering of φ bosons in the large N limit to determine first whether the scattering amplitude has a massless pole; we then study the tensor structure of the associated propagator to show that it has the form appropriate for Einstein gravity.
The scattering amplitude is represented by  The reader may worry at this point that we have omitted many other possible diagrams that contribute at the same order in 1/N ; for example, one could take any diagram shown in figure 1a and simply append a single loop to any one of the internal lines. We will show that all other diagrams that might contribute to the scattering amplitude at leading order in 1/N are cancelled as a consequence of the way in which we impose the condition of vanishing cosmological constant. This constraint is necessary for the stability of the theory given that we have chosen a Minkowski-space background. We discuss the case where the cosmological constant is non-vanishing in an appendix.
Let us parameterize the two degrees of freedom in our theory that are relevant to this discussion. First, the background value of the fields X I need not have been chosen as in eq. (2.7); we now allow for a rescaling of X I , where c 1 is a constant. This change modifies the form ofh µν in eq. (3.2) so that the operator ∂ µ φ a ∂ ν φ a is replaced by ∂ µ φ a ∂ ν φ a − c 1 η µν . In addition, while still expanding the JHEP03(2017)134 action in powers of 1/V 0 , we add an additional constant to the scalar potential, so that the original ∆V is replaced by ∆V = m 2 φ a φ a /2 − c 2 . (3.10) Since the action of our theory is exclusively a function ofh µν and ∆V , these modifications assure that every occurrence of the operator ∂ µ φ a ∂ ν φ a is via ∂ µ φ a ∂ ν φ a − c 1 η µν and every occurrence of m 2 φ a φ a /2 is via m 2 φ a φ a /2 − c 2 . We may think of c 1 and c 2 as counterterms: any diagram involving a loop that is created by connecting a φ a line within the operator φ a φ a or within ∂ µ φ a ∂ ν φ a will be paired with another diagram that has a counterterm vertex in place of that loop. We impose the renormalization condition that the sum of these diagrams vanishes. This is merely a convenient way to organize the perturbative calculation, since a different choice of the counterterms is equivalent to a change in gauge and a shift in V 0 , which we have not yet fixed. In our scattering problem, this eliminates all other diagrams that are leading order in 1/N that are not shown in figure 1a. Our renormalization condition is equivalent to replacing the operators φ a φ a and ∂ µ φ a ∂ ν φ a in eqs. (3.7) and (3.8) with normal-ordered operators : φ a φ a : and : ∂ µ φ a ∂ ν φ a :, respectively.
What is convenient about the problem as we have set it up is that the diagrams in figure 1a can be represented by the recursive formula shown graphically in figure 1b. 5 Each diagram in the infinite sum involves common functions of the external momenta that we will factor out in doing the resummation. We define for inwardly (or outwardly) directed external momenta p 1 and p 2 . At leading order in 1/N , the Feynman rules for the external lines follow solely from the action of T µν on the external states; hence the form of eq. (3.11) is determined by eq. (3.7), summing over the ways in which the fields can annihilate (or create) incoming (or outgoing) scalar bosons. We then write the scattering amplitude in the form where The recursive relation represented by figure 1b is given by (3.14) The first term on the right-hand-side follows from the tree-level amplitude and is given by while the kernel K λκ ρσ (q) can be found by doing a one-loop calculation corresponding to the portion of the second diagram in figure 1b that connects to the shaded blob. We find .
(3.16) 5 A similar analysis in a theory with emergent electromagnetic interactions was performed in ref. [5].

JHEP03(2017)134
In the appendices, we present the calculation of the kernel. The result may be written as eq. (B.5), where we have omitted terms that vanish as a consequence of the identity 6 q µ E µν (p 1 , p 2 ) = 0 . Defining the quantity λ by we may substitute our solution for the kernel into eq. (3.14), from which we obtain We wish to determine whether the existence of a massless pole can be inferred from eq. (3.20). We make the fine-tuned choice which corresponds to λ = 1. 7 In this case, the left-hand-side of eq. (3.20) vanishes, and it immediately follows that the amplitude contains a massless pole: where we have used the tree-level amplitude given by eq. (3.15). Notice that eq. (3.22) has the tensor structure that one expects for the graviton propagator of Einstein gravity in de Donder gauge. Our choices of c 1 , c 2 and V 0 are tuned to give a vanishing cosmological constant. Notice that if we were to hold c 1 and c 2 fixed while replacing V 0 → V 0 + δV 0 , we would find that δV 0 appears in two ways in the action: (1) in a constant that multiplies every occurrence of the Minkowski metric, which can be removed by a rescaling of the fields and of m 2 and (2) in a shift of the scalar potential. In appendix C, we show that the latter indicates the presence of a cosmological constant and we discuss the consequences for the theory.
We now return to the physical interpretation of our renormalization conditions on the counterterms c 1 and c 2 . Consider first how gravity couples to matter in our theory. Using the solution for the metric g µν , we may identify the composite operator that represents the fluctuation about Minkowski space; it can be expressed in terms of T µν as

JHEP03(2017)134
where we define the tensor structure To determine the couplings to φ that are linear in h αβ in an effective theory in which h αβ is treated as a fundamental field, we note that there are two possible identifications of the composite operator h µν in eq. (3.8) that are relevant. Including both, we invert eq. (3.23) and substitute, yielding which matches our expectation for the graviton coupling implied by linearized general relativity. At the order in 1/V 0 that we are working, a cosmological constant would be present if there were a non-vanishing vacuum expectation value proportional to η αβ for the operator T αβ . However, we have imposed renormalization conditions on c 1 and c 2 so that diagrams formed by closing a φ a line in either the operator φ a φ a or ∂ α φ a ∂ β φ a vanish identically. The diagrams that would contribute to T αβ are of this form; hence our renormalization choice fixes the cosmological constant to zero, to the order in 1/V 0 that we are working. In appendix C we consider the consequences of not imposing this tuning and we compare to the expectations of general relativity when one tries to expand about a Minkowski-space background in the presence of a cosmological constant. Finally, we may use our results to compute the D-dimensional Planck mass M P . The corresponding scattering amplitude computed via linearized general relativity is given by 8 . (3.27) It is important to keep in mind that we use dimensional regularization here as a placeholder for a generally covariant, physical regulator, associated with some high scale Λ. We do so by writing D = 4 − , with small, but finite and positive. In this case, the fourdimensional Planck mass is proportional to N/ m, which can be made large even when m is comparable to the masses of known particles. We comment on other physical regulators in the Discussion and Conclusions of this paper. 8 In general relativity the scattering amplitude follows from the linearized gauge-fixed action where we parametrized the metric fluctuation around flat space by δgµν = 2M 1−D/2 P hµν and introduced the source coupling via the standard definition T µν = (2/ |g|)δL/δgµν . This action is obtained by linearizing the Einstein-Hilbert term d D x( 1 2 ∂ρhµν ∂ ρ h µν − ∂µhνρ∂ ν h µρ + ∂µh µν ∂ν h − 1 2 ∂µh∂ µ h) and adding to it the de Donder gauge-fixing term d D x(∂ ν hµν − 1 2 ∂µh) 2 , and the coupling to matter. The propagator in the de Donder gauge is i/(D/2 − 1) P µν|αβ /q 2 , where P µν|αβ was previously defined in eq. (3.24).

JHEP03(2017)134 4 Graviton couplings beyond linear order
In this section, we begin addressing the question of how the non-linear graviton couplings in our model match one's expectations from general relativity. We focus primarily on the couplings that originate from the matter action; a treatment of the graviton selfinteractions, i.e., those originating from |g| R in general relativity, will be discussed below briefly, but studied in a longer publication [21].
Prior to the imposition of the condition that the energy-momentum tensor vanishes, our initial action is that of a real scalar field in curved space. Expanding about the Minkowskispace background, g µν = η µν + h µν , we find the usual h µν couplings to φ a that one would expect in general relativity. In our case, however, the subsequent constraint we impose on the energy-momentum tensor identifies h µν with a function of the scalar fields φ a . Using a diagrammatic argument below, we argue that it is justified to identify this composite operator as the operator that creates and annihilates gravitons; then the theory includes all the multi-graviton couplings that we expect to find in |g| L matter . The theory contains more than just gravitational interactions; however, we expect the higher-dimensional local effective interactions generated by integrating over scalar loops to be suppressed by the regularization scale. An analysis of these operators will be presented in ref. [21].
Consider any term in the expansion of |g| L matter written in terms of the φ a and h µν , prior to imposing our condition that the energy-momentum tensor vanishes. We wish to focus on one factor of h µν appearing in such a term at a time, so we will represent the generic interaction as where the factor of −1/2 is included as a matter of convenience and V µν ≡ V µν (φ, h αβ ), with the chosen h µν factored out. After we constrain the energy-momentum tensor of the original theory, we may express h µν in terms of the flat-space energy-momentum tensor, as given by eq. (3.23). Hence, we can write eq. (4.1) in the form which we should compare to the form of the quartic interactions in eq. (3.8), Now we consider the effect of attaching a chain of scalar loops to the T αβ that we obtained by substituting for the designated h µν . The relevant diagrams are shown figure 2. The circled cross represents the Feynman rule following from V µν (φ, h αβ ), which we will represent by the function E V (p 3 , . . . , p n ) µν , assuming there are n external lines in the diagram (with n − 2 not displayed in figure 2). Since each diagram will involve a factor of E ρσ on the left and E V µν on the right, we notice that the resulting amplitude has the form where A µν|ρσ is the same amplitude defined in our previous two-into-two scattering calculation. Note that the factor of E µν on the left is sufficient to eliminate the same terms dropped from the kernel in eq. (3.17). In comparing figure 1 and figure 2, the difference between eqs. (4.2) and (4.3) simply changes how each diagram in these figures terminates on the right. The part of A µν|ρσ that has a pole, eq. (3.22), was shown in the previous section to correspond to the graviton propagator. Hence, if one is interested only in the graviton coupling, the result in eq. (4.4) matches what one would find by treating the designated h µν in eq. (4.1) as a fundamental field. Note that the factor of E ρσ associated with the left-side of each diagram is what we expect given the form of eq. (3.25). Since the h µν appear in |g| L matter exactly as in general relativity, the effective interactions for the composite graviton match with those that we would obtain in the conventional formulation where h µν is identified as a fundamental field. Note that radiative corrections to the matter-graviton couplings are suppressed by powers of V 0 ∝ N ; in the large-N limit these are just the usual corrections expected in quantum gravity.
The preceding argument was made possible by the fact that we started with a form of the theory in which the composite graviton operator h µν appears in the tree-level action in the same way as in general relativity. We cannot apply the same reasoning in determining, for example, the three-graviton coupling, since we have no tree-level terms involving h µν that we could match to the cubic part of |g| R. The graviton self-interaction vertices necessarily arise via connecting scalar loops in the present model: for example, one would extract the three-graviton vertex from the diagram shown in figure 3. Since the form of the three-graviton vertex is dictated by general covariance, and should not be affected by our gauge fixing choice, at lowest order in 1/V 0 we expect to recover the same vertex as in general relativity. A similar argument has been applied to other theories where kinetic terms for an auxiliary-field metric are induced by radiative corrections [9]. This argument seems plausible since in a heat kernel approach the full nonlinear effective action may be obtained directly in terms of generally covariant operators [12], but that approach uses a UV heatkernel regulator rather than dimensional regularization. It should be possible to support this argument in the present context via an explicit evaluation of the amplitude given in figure (3). This calculation is more tedious than the one discussed in section 3 and in the appendix, and will be included with a more general study of loop corrections in ref. [21].

Discussion and conclusions
We have demonstrated how a constraint of vanishing energy-momentum tensor can lead to a metric-independent theory in which quantum gravity emerges as a nonperturbative artifact of regularization-scale physics. Vanishing of the energy-momentum tensor is closely related to metric-independence of the action, which is desirable for a background-independent description of quantum gravity. As an example, we have constructed a scalar field theory with vanishing energy-momentum tensor that has a perturbative low-energy description including scalar fields with a potential. In a large-N limit, we have explicitly demonstrated that scattering of scalar particles includes a massless spin-2 pole, corresponding to exchange of a massless composite graviton that couples to matter as in Einstein gravity. The gravitational coupling is determined by the short-distance regularization of the theory, which for definiteness in our calculations we have defined by dimensional regularization. We conclude by commenting on a number of issues related to the interpretation of these results.
• Although our analysis of the theory was done in the large-N limit for simplicity, we expect that the existence of the massless spin-two state does not depend on this limit, due to the general covariance of the theory. Furthermore, although we have not yet computed the gravitational self couplings in the effective theory, we expect those couplings to agree with general relativity up to corrections suppressed by the regularization scale. These issues deserve further investigation.
• It is worth commenting on renormalizability in this approach. Unless Einstein gravity is asymptotically safe [22], quantization of the Einstein-Hilbert action is nonrenormalizable: arbitrarily high-dimension operators must be included as corrections to the action in order that correlation functions remain finite as the cutoff tends to infinity, and there is no unique choice of the coefficients of the higher-dimension operators. However, if gravity emerges as a nonperturbative regularization-scale artifact, then the Planck scale is linked to the regularization scale, and higher-dimension operators JHEP03(2017)134 may remain suppressed by the regularization scale. Hence, renormalizability is not required by emergent gravity models in the usual sense. However, whether or not functional integrals in the theory defined by eq. (2.6) are well defined nonperturbatively remains an open question. Similarly, the consequences of introducing additional field dependence in the metric-independent action requires further exploration.
• As a quantum theory of gravity, analysis of the theory at short distances requires that the regularization procedure be treated as physical. For dimensional regularization, this implies that functional integrals are to be analytically continued in the number of spacetime dimensions, which is then to be fixed at some D = 4 − . The theory determines correlation functions at arbitrarily short distances, but due to the analytic continuation in D those correlators are not expected to satisfy the usual axioms of 3+1-dimensional relativistic quantum field theory; the same could be said for regularization by Pauli-Villars fields with "wrong" spin-statistics properties.
If we allow for the possibility that short-distance physics violates one or another property usually held dear, such as unitarity of the scattering matrix, then we have an avenue by which to address questions related to the singularities that plague classical Einstein gravity. One possibility is that spacetime remains continuous at arbitrarily short distances, and correlation functions define the physical content of the theory, but the particle-like description breaks down at short distances and correlators lose their interpretation in terms of scattering amplitudes. In order to address issues in quantum gravity related to spacetime singularities and spacetimes with horizons, such as the black-hole firewall puzzle [23], we need to move beyond the perturbative analysis of this paper and consider background field configurations for which the composite metric is far from the Minkowski metric.
• The cosmological constant problem is not immediately resolved in this theory. The composite graviton state couples to the energy-momentum tensor in the same way as the sum of diagrams in figure 2. Closing the external lines corresponds to the coupling of the metric fluctuation to the vacuum energy, and leads to the same instability as in general relativity. Hence, even though we showed that the static gauge-fixing condition was consistent with a classical perturbative expansion about a flat composite spacetime, the requirement of vanishing cosmological constant is necessary for perturbative stability of the flat spacetime.
• It is necessary to understand the coupling of gravity to the Standard Model in this approach. Coupling the Standard Model to an auxiliary vielbein, the vanishing energymomentum-tensor constraint might determine the vielbein in terms of Standard-Model fields, perhaps similar to the description in ref. [11]. We would again expect the composite vielbein to give rise to emergent gravitational interactions. The clock and ruler fields X I might be included in a complete theory, as well as additional scalar fields φ a . We note that the φ a masses are arbitrary and may be taken large enough in such a scenario to evade potential phenomenological bounds.

JHEP03(2017)134
• Finally, we note that the Dirac-Born-Infeld action including the gauge field [24][25][26][27], which describes the dynamics of bosonic D p=D−1 branes (with the Ramond fluxes, the dilaton and the B-field turned off for simplicity), is also metric independent and diffeomorphism invariant. Quite manifestly it is similar to our starting point (2.6). The X's and φ's in (2.6) would be interpreted respectively as coordinates along the Dbrane and transverse to the D-brane. In this sense we would be looking at a D p=D−1 brane embedded in a D + N dimensional spacetime, with a flat bulk Minkowski metric. A difference between our model and the DBI action is that in the latter the transverse scalars are massless, since they are Goldstone bosons for the broken translational symmetry, whereas in our theory the masses of the φ's are nonvanishing.
It is tempting to speculate that taking to zero the scalar mass m will not change the conclusion that generic covariant regulators will generate emergent D-dimensional gravity. With the string scale playing the role of a regularization scale in the theory, we would therefore expect D-branes to support brane-localized gravity. Thus, string theory appears to provide another alternative to the realization of the brane-induced gravity scenario of Dvali, Gabadadze and Porrati [28], where dynamical bulk gravity gets localized on the brane. 9 This possibility was also noted in [30].

Acknowledgments
The work of CC and JE was supported by the NSF under Grant PHY-1519644. The work of DV was supported in part by the DOE grant DE-SC0007894. DV would also like to acknowledge the hospitality of William and Mary Physics Department during the completion of this paper.

A Integrals
In this section of the appendix we list the integrals used in the evaluation of the kernel defined in section 3. 9 In another related work, [29] computed the open string partition function in the presence of a potential that localized the open string endpoints on a D-brane, whose embedding was specified in terms of some transverse scalars Y i . The resulting effective action, as a function of the scalars Y i which were nondynamical, was expressed in terms of the induced metric gµν = ηµν + ∂µY i ∂ν Yi, and contained a volume term, and higher derivative terms which involved the extrinsic curvature of the brane. One of these terms was the Einstein-Hilbert action for the induced metric, leading to a picture of open-string induced gravity on the brane.

JHEP03(2017)134
where for the sake of clarity we mention that the parentheses on the left hand side denote the argument of the tensor. Switching back to the mostly minus Minkowski metric used in the paper, the integral we want to evaluate to order (q 2 /m 2 ) 2 is .

(B.2)
This is related to the kernel K as follows (3.16): To this end we employ the integrals computed in appendix A, and find The kernel becomes In the main text we drop the q µ,ν = (p 1 + p 2 ) µ,ν dependent terms from the kernel, since for the purpose of our calculation we need to contract the kernel with the external line factor E µν (p 1 , p 2 ) which satisfies (3.18) q µ E µν (p 1 , p 2 ) = 0.

C The cosmological constant
If the counterterm c 2 in eq. (3.10) is not tuned as described in section 3, so that c 2 is replaced with c 2 + ∆c 2 , then the effective energy-momentum tensor that couples to the composite graviton becomes : ∂ µ φ a ∂ ν φ a : −η µν 1 2 : ∂ α φ a ∂ α φ a : − 1 2 m 2 : φ a φ a : +∆c 2 . (C.1) We choose to leave the gauge choice eq. (3.9) unchanged; otherwise we would need to take into account the associated wavefunction renormalization of the fields φ a in this analysis, but would be led to the same conclusions. We assume that ∆c 2 V 0 . In the coupling L int of eq. (3.8), we have the following additional interaction, at the same order in 1/V 0 to which we have been working:
As illustrated in figure 5, ∆c 2 does not contribute to scattering amplitudes at leading order in 1/N . Hence, to this order the argument presented in section 3 for the massless graviton pole remains unchanged. However, the constant ∆c 2 in the scalar potential acts as a source for the composite metric h µν , as per the diagram in figure 6. This source corresponds to the tadpole instability of the Minkowski-space vacuum in the presence of the cosmological constant. Hence, the cosmological constant problem appears to remain in this emergent gravity scenario.
Open Access. This article is distributed under the terms of the Creative Commons Attribution License (CC-BY 4.0), which permits any use, distribution and reproduction in any medium, provided the original author(s) and source are credited.