Quantization of Hořava gravity in 2+1 dimensions

We study quantum corrections to projectable Hořava gravity with z = 2 scaling in 2 + 1 dimensions. Using the background field method, we utilize a non-singular gauge to compute the anomalous dimension of the cosmological constant at one loop, in a normalization adapted to the spatial curvature term.


Introduction
Considerable effort has been devoted to the study of Hořava gravity since it was introduced in [1,2]. Being renormalizable by naïve power counting, Hořava's theory constitutes a candidate for an ultraviolet-complete theory of quantum gravity. In spite of some work [3][4][5][6], nonetheless, as yet there have been no fully satisfactory quantum computations; in fact, perturbative renormalizability of one version -the "projectable" model -was established only recently [7]. The purpose of the present paper is to take a step forward in understanding quantum corrections to Hořava gravity by making a careful computation of a one-loop quantity working in non-singular gauges. (What we mean by this is explained in sections 2.2 and 3.1.) More specifically, the model we consider is z = 2 projectable Hořava gravity in 2 + 1 dimensions, and the quantity we compute is the anomalous dimension of the cosmological constant. 1 1 In anisotropic models, the effective coefficients of the temporal and spatial kinetic terms can scale differently -i.e., the dispersion relation runs with scale. This running can be captured by fixing the form in which either the energy or the spatial momentum appears in the dispersion relation. We compute the anomalous dimension of Λ with respect to a normalization condition that fixes the form of the spatial momentum contribution.

JHEP06(2017)004
Hořava gravity is constructed so that, at high energies, the classical action has anisotropic scale invariance with the dynamical critical exponent z: x → b x. (1.1) As our interest is in z = 2, we take the engineering dimensions of the time and space coordinates to be [t] = −1, In this convention, energy is of dimension one. The z = 2 theory is renormalizable in 2 + 1 dimensions [7]. The spacetime manifold is equipped with a foliation by leaves of codimension one, corresponding to the surfaces of constant time. Its geometry is naturally parametrized using the ADM variables -a spatial scalar N (the lapse), a spatial vector N i (the shift), and a spatial metric g ij . The classical scaling dimensions of the fields are The gauge symmetries are the diffeomorphisms that preserve the foliation. We parametrize the infinitesimal transformations by (Z, X i ), that act on the fields by A proper understanding of Hořava gravity requires a careful treatment of its gauge fixing. To this end, it is useful to begin with the simplest model possible. It is tempting to begin with the conformal case in 2 + 1 dimensions, because it has no local propagating degrees of freedom. Unfortunately, not only does it require the "non-projectable" version of the theory, which has second class constraints and their attendant difficulties, but also it raises the thorny issue of gauge anomalies for the Weyl symmetry.
A more modest starting point is "projectable" Hořava gravity in 2 + 1 dimensions. Projectability is the condition that N = N (t) be a function of time but not of space, so that it is constant on each spatial slice. We assume this condition for the remainder of the paper. The 2 + 1 dimensional projectable case is more than just a toy model for understanding the qualitative behavior of the more realistic 3+1 dimensional non-projectable theory. Mapping out the renormalization group (RG) structure of the projectable theory is important to further understand the phases of gravity, both in the context of Hořava gravity and the Causal Dynamical Triangulation approach to quantum gravity [8,9].
The action is written in terms of quantities invariant under those diffeomorphisms preserving the foliation of spacetime, namely scalars built from the intrinsic and extrinsic JHEP06(2017)004 curvatures of the leaves of the foliation and their covariant derivatives. The intrinsic curvature of a two-dimensional leaf is completely determined by its spatial Ricci scalar R. The extrinsic curvature is captured by the tensor (1.6) where ∇ i is the covariant derivative with respect to g ij . The most general z = 2 action invariant under (1.5) is where K = g ij K ij . Since d 2 x √ gR is a topological invariant in two dimensions, ρ does not appear in the local equations of motion, but only in the global Hamiltonian constraint arising from time reparametrization symmetry. As a result ρ cannot contribute to the perturbative beta function, and so we drop this term in what follows. 2 In general dimension, projectable Hořava gravity has a transverse traceless tensor mode and a scalar mode. Requiring the tensor polarizations to have a good dispersion relation around flat space then implies that γ > 0. Requiring the dispersion of the scalar also to be healthy imposes the constraint λ < 1 2 or λ > 1. (1.8) In 2+1 dimensions, however, there are no tensor modes. We then have the option of setting γ to be negative when 1 2 < λ < 1. The propagating spectrum of the theory is then healthy, at least classically. We do not worry about this explicitly in what follows, although our final result makes sense in this parameter region.
In this paper, we will compute contributions to the effective action using the background field method. In this method, fields are split into a sum of two terms: a classical background value, and quantum fluctuations of typical size 1/2 . For the action (1.7), the role of is played by κ 2 . This leads us to expand 9) where N , N i and g ij are background fields and n, n i and h ij are fluctuations around the given background. Gauge transformations can also be expanded in powers of κ,

JHEP06(2017)004
In this gauge, the action of ζ and ξ i (to linear order in κ) is Here, ∇ i denotes the Christoffel connection for g ij . We can use ζ to set n ≡ 0; since n is independent of space and so has only one degree of freedom per spatial slice, it does not contribute divergences. For our purposes, therefore, we can ignore the contribution from n to the path integral.
In the following, we will work only on backgrounds that are time-independent. We express the partition function in terms of functional determinants by integrating out the quantum fluctuations n i and h ij , and the gauge-fixing ghost modes. The one-loop effective action is then evaluated using heat kernel techniques. This will allow us to compute some (but not all) of the one-loop beta functions in the theory. To fully understand the RG properties of the theory at weak coupling (and in particular, determine whether the theory is asymptotically free), it is necessary to evaluate the heat kernel on background geometries with a time-dependent metric. We leave this to future work.
Previous work on the one-loop effective action in gravity with anisotropic scaling [10] overlooked crucial contributions from the gauge-fixing sector of the theory, a problem exacerbated by dropping from the partition function altogether singular determinants that did not cancel out in their analysis. We show that such confusion can be avoided by an appropriate choice of gauge. The gauge-fixing methods we developed have, in the meantime, appeared in a more general form in the work of [7], which applied them to show the renormalizablility of projectable Hořava gravity. We take advantage of their more general gauge in section 2.2 for reasons of clarity, although the bulk of our computation uses our more restrictive original gauge.
Section 2 develops the gauge-fixing method and field parametrizations we use in the remainder of the paper in the simpler context of linearized theories. Before embarking on the gravitational calculation, we begin in section 2.1 with a warm-up -free U(1) gauge theory in D + 1 dimensions with z = 2 scaling at short distances. One natural choice in this context, used in [11], is temporal gauge. Here, we utilize a gauge choice that manifestly respects the z = 2 scaling symmetry. Generalizing this gauge-fixing procedure to the gravitational case will lead us in section 2.2 to the same sort of gauge-fixing condition used by [7] in proving perturbative renormalizability of projectable Hořava gravity. Section 2.3 uses these results to compute the dependence of the one-loop effective action on the cosmological constant, which illustrates how the effective action can depend on gauge, and how to extract the correct gauge-invariant effective action.
Section 3 turns to computations in curved space using the background field method. There, we compute the partition function on static on-shell curved backgrounds (R = const, ∂ t g ij = 0) supported by non-vanishing Λ. Working with an on-shell background enables us to systematically disentangle the physical and unphysical modes and observe explicitly the cancellation of the unphysical modes among themselves. We give an explicit expression JHEP06(2017)004 for the physical dispersion relation, which generalizes the flat space result. We normalize the gravitational field such that γ/κ 4 is constant at all energy scales. With respect to this choice of normalization condition, we are able to determine the anomalous dimension of Λ. Extracting the beta functions for γ, λ and κ requires working on backgrounds that depend on time, which we leave to future research.

Gauge fixing in theories with anisotropic scaling
In gauge theories exhibiting an anisotropic scaling symmetry of the form (t, x) → (b t, b 1/z x), it is desirable to choose a gauge-fixing condition that respects this symmetry. This is especially true in models at their critical dimension, for which standard gauges -in particular, Lorenz gauge -may not be renormalizable.
In some simple cases (e.g., free Maxwell theory), there is no problem with singular gauges, such as the temporal or Coulomb gauges, which are in fact invariant under the scaling symmetry for any value of z. When the theory is coupled to gravity, however, such gauges can become problematic. For example, in temporal gauge the Faddeev-Popov determinant is det(∂ t ). While in the flat case this determinant can be dropped, in the gravitational case it couples non-trivially and should not be ignored. However, such operators have no dependence on large spatial momenta, leading to uncontrolled ultraviolet divergences. Moreover, this problem persists in both dimensional regularization and heat-kernel based methods. Such gauges therefore give rise to ambiguities, which need to be resolved in a manner consistent with BRST symmetry. From a more pedestrian perspective, our strategy ensures that the gauge-fixing Lagrangian, which is quadratic in the gauge-fixing condition, is of the same order in derivatives as the original Lagrangian. Thus, the two can be combined more seamlessly.
In this section, our goal is to introduce 3 such gauges in linearized z = 2 gravity. We first illustrate the process in free anisotropic U(1) gauge theory. This serves as a warm-up to the second case of z = 2 projectable Hořava gravity in 2+1 dimensions linearized around flat space. We apply these results to make a simple quantum computation. Section 3 will be concerned with the generalization to static backgrounds in the background field method.

U(1) gauge theory
We begin with free U(1) gauge theory in D + 1 dimensions exhibiting z = 2 scaling in the ultraviolet (UV) and z = 1 in the infrared (IR). The gauge field is a U(1) connection on Aristotelian spacetime [12]. The time and space components, A 0 and A i (i = 1, · · · , D), have gauge transformations, δA 0 =ζ, with ζ(t, x) an arbitrary scalar function. The invariant field strengths are

JHEP06(2017)004
At the ultraviolet z = 2 Gaussian fixed point, the engineering dimensions of the gauge fields are The basic most generic action with this scaling in the UV that is invariant under both the spacetime and gauge symmetries (including parity and time-reversal) is where v is the "speed of light" in the infrared. In components, the action becomes and O is the generalized d'Alebertian operator, (2.7) A natural z = 2 generalization of the Lorenz gauge is given by the gauge-fixing functional 4 f To quantize the theory with this gauge-fixing, we should further introduce a pair of fermionic ghosts (b, c), a bosonic auxiliary field Φ, and the fermionic BRST differential s acting on the fields as A generalized R ξ gauge-fixing action based on (2.8) can now be obtained from a gaugefixing fermion of the form (2.10) Note that, unlike standard R ξ gauge, if we wish to avoid introducing dimensionful parameters then D must be a differential operator of dimension one. The BRST-exact action is

JHEP06(2017)004
The resulting BRST-invariant gauge-fixed action is giving the quantum partition function As in the case of the standard R ξ -gauge procedure, the partition function is independent of D. We demonstrate this explicitly in appendix A. Setting D = −∂ 2 + v 2 is a particularly nice choice, as it eliminates the cross-terms between A 0 and A i , after integrating out Φ. Redefining the A 0 field via A 0 → √ D A 0 results in a Jacobian J A 0 = (det D) 1/2 , which cancels the factor of (det D) −1/2 produced by the integral over Φ. The action then becomes 14) The overall sign in front of the piece quadratic in A 0 in (2.14) is negative, so we must Wick rotate A 0 when we rotate t. The partition function evaluates to This represents D − 1 physical propagating modes with dispersion relation Before we move on to Hořava gravity, we make the following comment. As in the Lorentz-invariant theory, one can diagonalize the kinetic operator of (2.5) explicitly in field space, without gauge-fixing. There is one pure gauge mode on which the operator vanishes completely. There is also one unphysical mode which gets a wrong-sign dispersion relation. The rest of the modes should then reproduce the correct dispersion relation (2.16). We perform this exercise for illustrative purposes in appendix B. 5

Hořava gravity around flat space
We now turn to the linearization of z = 2 projectable Hořava gravity in (2 + 1) dimensions around flat space, with g ij in (1.9) set to δ ij . The flat background is on-shell if the cosmological constant Λ is set to zero. However, since we are also interested in the Λ dependence of the off-shell effective action, we will allow for a nonzero Λ.
The action is that of (1.7), with ρ = 0. The quadratic part 6 of the action is 5 K.T.G. would like to thank Laure Berthier for this point. 6 Since we are interested in the effective action, we drop the linear part, which is non-vanishing when Λ = 0.

JHEP06(2017)004
where h ≡ h ii , and the matrix S (2) is the second functional derivative of the action. Explicitly, Intuitively, one can think of this theory as "adding a spatial index" to the U(1) gauge theory of the previous section: n i is analogous to A 0 , and h ij to A i . Likewise, the gaugefixing functional f i , ghost fields b i and c i , and bosonic auxiliary field Φ i all carry a spatial index. The BRST differential s acts as In analogy with the U(1) theory, we choose the gauge-fixing fermion where D ij is some spatial differential operator of dimension one, and f i is a gauge-fixing functional. As pointed out in [7], the most general such operator is where u 1 and u 2 are constants. The analog of the gauge-fixing condition (2.8) reads f i =ṅ i − D ijk h jk , where D ijk is a spatial differential operator of energy dimension 3 2 (e.g., containing three spatial derivatives). Forcing the cross-terms between n i and h ij to vanish upon integrating out Φ i uniquely determines D ijk to be where u = u 1 +u 2 . As before, the final result is independent of the particular choice of D ijk . The analog of the BRST-exact action (2.11) is We can immediately read off the ghost partition function, The rest of the action, after integrating out Φ i , called the "effective" part, reads Note that various field components need to be Wick-rotated as well as the time when performing the path integral. The contribution of n i to the partition function is Next, we compute the contribution to the partition function from h ij . We first decompose h ij as where H ii = 0. We decompose H ij further as with constraints

JHEP06(2017)004
In two spatial dimensions, the transverse traceless component H ⊥ ij has no local degrees of freedom, and in flat space is forced by boundary conditions to vanish. Furthermore, in two dimensions, one can parametrize η i as (2.34) Thus, H ij is finally parametrized as We require the Jacobian for this change of variables. The Jacobian for (2.31) is a constant. The Jacobian for (2.35) is computed in appendix D, We can eliminate this Jacobian altogether by changing variables from η and σ to The action for η is just Therefore, the contribution of η to the partition function is . (2.40) Meanwhile, h and σ remain coupled via the action The matrix S h σ is diagonal only for the choice For general gauge parameters, S h σ can't be diagonalized locally. When Λ = 0, however, the determinant itself factorizes neatly, The operator O phys is independent of the gauge-fixing parameters u 1 and u 2 .
In summary, the modes corresponding to the various dispersion relations are h, σ, n i and ghost : η, n i and ghost : For these to have the right sign dispersion relation requires Note that the "nice gauge" of [7] is when all three of the dispersions, including the unphysical ones, are actually identical. This condition is satisfied if and only if Finally, the total on-shell partition function is the product of (2.26), (2.30), (2.40), (2.44) and the extra factor of (det D ij ) −1/2 from integration over Φ i . The result simplifies greatly, This represents one physical degree of freedom, with dispersion relation This dispersion relation is healthy when λ > 1 or λ < 1 2 . Note that this degree of freedom is a linear combination of h and σ. Therefore, it will not be captured entirely if one neglects everything except the trace component of h ij , as was done in [10]. When λ > 1 (and γ > 0), the overall sign in front of (2.49) is negative and we must Wick rotate the field 7 corresponding to O phys .
Once again, this dispersion relation can be derived without regard to a specific gaugefixing procedure, as in the case of the U(1) gauge theory. Details are given in appendix B.

One-loop effective action with a nonzero cosmological constant
Let us calculate the contribution to the determinants of first order in Λ. Our object of study is the effective action Γ(ϕ), where ϕ denotes the expectation values of all fields Φ of the gravitational theory. Expanding in , (2.51) 7 In general, this field is some combination of h and σ. In the "nice" gauge (2.48), this field is just h.

JHEP06(2017)004
by standard methods the one-loop quantum effective action takes the form 8 is the second functional derivative of S. Since this is a gauge theory, we must also include ghost contributions after gauge-fixing, leading to the standard expression Note that the only dimensionful parameter present is Λ itself, with dimension [Λ] = 2. As a result, the only contribution Λ can have to the logarithmic divergence (and therefore to the one-loop beta functions) is proportional to Λ. To evaluate it, it therefore suffices to evaluate the first derivative of the partition function Z with respect to Λ. Separating out the Λ dependence of S (2) , Evaluating the relevant traces, we obtain the integral form (2.61) and constants Here, I = 1 corresponds to the η contribution, while I = 2, 3 arise from the (h, σ) sector. Later on in this paper we will use heat kernel methods, which preserve diffeomorphism invariance. It is difficult to use the heat kernel here, however, because we have not diagonalized M h, σ . (Note, however, that in the diagonal "nice" gauge this is not a problem.) Although it breaks gauge symmetry and modifies the infrared behavior of the theory, to extract the coefficient of the logarithmic divergence it suffices to use a cutoff regularization. We integrate over all ω and introduce a cutoff k * in k. In addition, the denominators have an implicit +iε, specifying the appropriate Wick rotatation ω = iω E . The integrals evaluate to giving the final result The contribution of Λ to the effective action (2.54) is therefore This is obviously gauge-dependent. As we will discuss in section 3.6, this gauge dependence should be eliminated by a field strength redefinition for the background metric g ij , We will utilize this field redefinition in the subsequent section in order to extract the key gauge-independent information.

Time-independent curved background
In time-independent backgrounds (ġ ij = 0), the background values of the extrinsic curvature and the Riemann tensor are K ij = 0 and R i jk = R i jk (x), respectively. In two spatial dimensions, the Riemann tensor is determined by the scalar curvature,

JHEP06(2017)004
By dimensional analysis, the only divergence sensitive to ∇ i R that can appear in the effective action is R, which is a total derivative. Therefore, it suffices to take R to be constant. Consider the action (1.7) with the coupling constant ρ set to zero, We now expand each term in this action to quadratic order in κ. With 5) and The action can be put on-shell by imposing the equation of motion for the background field g ij . This essentially sets the cosmological constant to be Plugging (3.7) back into the action eliminates the tadpole terms linear in κ. The on-shell action is Considering such an on-shell action enables us to observe explicit cancellations between ghosts and non-physical modes, reducing the computation of the effective action to that of a single scalar functional determinant.

JHEP06(2017)004
We offer one caveat: since we do not impose the constraint equation generated by the gauge choice n = 0, by "on shell" we actually mean the background satisfies the g ij equations of motion. To render the background (3.7) fully on-shell requires imposing the further condition ρ = 2γR. As noted in the introduction, however, neither the lapse nor the value of ρ affects the local divergences, and therefore we can ignore both in the computation at hand.

BRST quantization
We now turn to the problem of gauge-fixing. We will apply the BRST formalism. Instead of classifying the most general gauge-fixing conditions, let us take a more minimalistic approach and construct a gauge-fixing condition such that the cross terms between n i and h ij cancel in the BRST action. For this purpose, it is sufficient to set the gauge-fixing Here D 1 and D 2 are local operators of dimension one, which we will take to be linear combinations of the diffeomorphism-invariant objects R and ≡ ∇ i g ij ∇ j . As we reviewed in section 2.2, equation (3.9) is not the most general gauge choice consistent with background diffeomorphism invariance: for example, one may also include in f i terms of the form (3.10) In the zero curvature limit, this extra term is the same as the u 2 term in (2.23). In the flat case, if one requires that nonphysical modes have a right sign dispersion relation, the conditions derived in (2.47) must be satisfied. For λ > 1, a nonzero u 2 is indispensable for these conditions to hold. On an on-shell background, the one-loop contributions from nonphysical modes cancel exactly in the partition function, and it is not necessary to include (3.10) in the gauge-fixing condition. When on-shell, we can focus on λ < 1 2 and adopt the simpler gauge-fixing condition in (3.9). Evaluating the partition function will result in a gauge invariant expression that is analytically continuable to λ > 1.
BRST quantization proceeds by introducing a ghost field c i associated to the generator of infinitesimal diffeomorphisms. The BRST differential s acts on the physical fields in the same way as the linearized diffeomorphisms in (1.12): We also require a cohomologically trivial BRST pair (b i , Φ i ), with fermi and bose statistics respectively. The ghost sector has BRST variations Gauge-fixing actions are given by the BRST differential of a gauge-fixing fermion. We take the gauge-fixing fermion

JHEP06(2017)004
which gives the BRST-exact action This action is associated to a gauge-fixing condition of the form (3.9), except that we have replaced the δ-function type by a gauge of generalized R ξ type. We have introduced auxiliary fields Φ i of dimension 1 2 and a local operator D of dimension 1. We choose the following expression for the operator D: where and The BRST partition function is Next, let us integrate out the auxiliary field Φ i in S BRST . Since Φ i only appears in S g.f. , we can focus on the following piece in the partition function, Here, Φ i is not dynamical and can be eliminated by imposing its equation of motion, The resulting action after eliminating Φ i in S g.f. is From now on, we will take S g.f. to denote the expression (3.22), even though it is different from the original expression in (3.17).

JHEP06(2017)004
Integrating out Φ i in the partition function (3.20) also contributes a functional determinant. To evaluate this determinant, we first make the change of variables with ε ij = √ g ij the covariant Levi-Civita symbol for g. The Jacobian is given by (D.10) in appendix D, In terms of φ and φ, the part of (3.17) that is quadratic in Φ i can be written as To derive that the cross term between φ and φ is zero, we used the form (3.15) of D and applied Identity 1 in appendix C. Therefore, where the functional determinant is evaluated to be (3.27) Therefore, the final expression for the Φ i contribution (3.20) is where S g.f is given by (3.22) and We applied Identity 1 for the second equality in (3.29). Finally, we determine the operators D 1 and D 2 in (3.13) by requiring that the cross terms between n i and h ij cancel in the sum S + S g.f. , with S g.f. set to the expression in (3.22). The kinetic contribution in the action S comes from The part contained in S K that is quadratic in terms of the fluctuations is The cross terms in S K are

JHEP06(2017)004
The contributions to the cross terms from S g.f. are These two contributions, (3.31) and (3.32), cancel if Since D has been defined in (3.15), this fixes both D 1 and D 2 .

The ghost sector
The integration over the ghosts in the partition function can be treated separately. From (3.18) we obtain We would like to evaluate the partition function  In terms of the fields b, b, c and c, the ghost action becomes where Therefore,

The non-ghost sector
Now, we would like to come back to examine the non-ghost part in the action S BRST , namely, the combined contribution from S + S g.f. . It is useful to take the following decomposition of the metric fluctuation h ij such that where H ij is a traceless 2-tensor, and Note that the quantum field H ⊥ ij is both traceless and divergenceless. In 2 + 1 dimensions, H ⊥ ij encodes only global information about the geometry of the spatial slice (the moduli of the Riemann surface), and carries no local degrees of freedom. Therefore, we can drop H ⊥ ij without affecting the β-functions. The constraint on η i can be solved by parametrizing η i as The Jacobian from the transformation (3.43) is computed in (D.15), Under the decomposition (3.43), we have The full one-loop BRST partition function can be written as For u 1 > 0 (and λ < 1 2 ), we must Wick rotate n i as well as the time when performing the path integral.
First, let us focus on Z n . We decompose n i into scalar components ν and ν as follows, We choose this particular decomposition in order to make the action (3.48a) local. The corresponding Jacobian is Under this parametrization, we obtain Collecting these results gives the partition function of n i , (3.54) Contributions from η, σ and h can be read off of the actions (3.48b)-(3.48c) (in the h, σ sector the differential operator is a 2 ×2 matrix, whose determinant we take directly) and give

Reduction to physical spectrum
Let us collect the results that we have derived above. The BRST partition function Z BRST is given by where,

JHEP06(2017)004
and (3.59b) The operators O g , O g and O phys take the form The full BRST partition function reduces to It is reassuring that the final result is gauge independent and all singular prefactors simply cancel. This partition function counts exactly one physical degree of freedom. On the other hand, on an off-shell background there is no reason to expect the result to reduce to a single functional determinant, and the analysis would be more difficult. While the preceding discussion is formally correct, some care must be taken with analytic continuation to ensure that the path integral converges properly. Requiring that O g give rise to a sensible dispersion relation gives u 1 > 0; for O g , this requires that λ < 1. However, both of these operators drop out in the final BRST partition function, and the singular behavior for O g (when λ > 1) can be fixed by modifying the gauge-fixing condition (3.9). Working on-shell gives us the luxury of ignoring this issue: both the operators O g and O g cancel out in the final BRST partition function.
All that remains is the determinant of O phys in (3.61), whose evaluation requires an appropriate choice of contour. The coefficient of ∂ 2 t in O phys has a healthy sign for λ < 1/2, in which case the standard contour will do. For λ > 1 on the other hand, when we perform Wick rotation we must also rotate the field; this is perhaps not surprising, since a similar rotation must be done for the scale factor in general relativity to get a well-defined Euclidean path integral.
In momentum space, we obtain the following dispersion relation for the physical degree of freedom:

JHEP06(2017)004
Note that there are values such that the right-hand side is negative, indicating instability. On the sphere (R > 0), at most one unstable mode can arise, namely the zero-momentum mode which is unstable for λ > 1. 9 More troubling is the case where λ > 1 and R < 0, since as λ → 1 + , the range of momenta with unstable dispersion will grow arbitrarily large. Nonetheless, provided that the UV scale is much larger than R/(1 − λ) this will not affect the divergences of the theory, and so for the purposes of computing the beta function we can ignore any instabilities in the low momentum modes.

Evaluation of the heat kernel
It remains to compute the determinant of (3.61), which we will do using zeta function regularization. The real time quantum effective action is . To evaluate divergences, we will use the standard heat kernel representation which gives us the following representation of the one-loop effective action,

JHEP06(2017)004
Our background is a product geometry R × M 2 , so we decompose |t, x = |t ⊗ |x . Expanding |t in Fourier modes allows us to write Here we have defined I O (τ ; x) = x|e −τ O |x for any spatial differential operator O and set . (3.72) Note that the ω-integral converges after Wick rotation ( t ≡ it, ω ≡ −iω). Performing the integral over the frequency, By rescaling τ → τ /A, we obtain The spatial term I V can be evaluated by using the results of [15], which computed the divergent contributions due to operators of the form In our case, Expanding I V in powers of τ defines the Seeley-Gilkey coefficients, The logarithmic divergence comes from the m = 2 term. The computation of the Seeley-Gilkey coefficient a 2 of [15] yields for T i = 0, The log divergence can be evaluated by introducing a cutoff µ −4 for the τ -integral, which gives d ds s=0

JHEP06(2017)004
Inserting this into the expression for Γ 1 gives the one-loop logarithmic divergence of the effective action on our background: (3.82)

Renormalization for γ and Λ
So far, we have evaluated the one-loop effective action over two different background geometries, both of which are described by a time-independent metric: • The Aristotelian spacetime with a nonzero cosmological constant Λ = 0. This background geometry is off-shell, i.e., the background metric does not satisfy the associated background equations of motion. The effective action was evaluated in (2.66). The covariant expression is contains gauge dependence. Although (3.84) was computed using a sharp cutoff, the coefficient of the logarithmic divergence is universal, so we can use this result in studying the logarithmic divergence that arose in zeta function regularization.
• A background geometry with a time-independent metric but a nonvanishing Riemann tensor. We study the on-shell action with Λ set to be The effective action is given in (3.82): This result is on-shell, and therefore guaranteed to be gauge-independent.
Since Y Λ is gauge-dependent, we cannot use Y Λ by itself to extract physically meaningful information. Our goal will be to eliminate this gauge dependence and identify a physical quantity that can be extracted from Y . We begin by examining the effective action evaluated on an off-shell time-independent background. We expand to one-loop order, keeping only the logarithmic divergence: Note that K ij = 0 for a time-independent background. From (3.83) and (3.86), we obtain The effective action Γ on a time-independent background can be written as As we noted, the naïve off-shell effective action (3.91) depends on our choice of gauge parameters. In fact, as a function on the space of background metrics, the effective action is gauge-independent, but the parametrization of field space can depend on gauge. Such dependence can therefore be removed by a field redefinition. (For example, see [14,16].) In general, these field redefinitions could include curvature terms. In our case, however, for dimensional reasons it suffices to rescale the metric. Under the rescaling, Under this rescaling, the effective action becomes To extract beta functions requires specifying a normalization condition that fixes the field rescaling. First, let us choose the normalization condition such that the coefficient of the R 2 term is set to one. This fixes the field rescaling C to be thereby turning the effective action into Indeed, the gauge-dependent contribution Y Λ drops out altogether from this last expression. The factor (1 − κ 2 Y ) can be absorbed into the renormalization of Ω. We are working in bare perturbation theory, so that the physical coupling Ω ph is related to the bare coupling Ω by Ω ph = (1 − κ 2 Y )Ω. Then, the anomalous dimension of Ω is (3.98)

JHEP06(2017)004
It is interesting to note that the running of Ω is independent of any field rescaling defined in (3.92). A simple analysis is helpful for understanding this observation. Throughout the paper, we have taken the scaling dimensions of time and spatial coordinates to be −1 and − 1 2 , respectively. In a more fundamental picture, however, we assign two independent dimensions, T to time, and L to length of space. In this latter convention, we have which suggests that Ω is independent of a rescaling of spatial coordinates. Further note that the rescaling of g ij can be absorbed completely into a rescaling of spatial coordinates.
Hence Ω should not change under the field redefinition of the spatial metric.
As we have seen in (3.98), an off-shell time-independent background provides us with only one piece of RG information. There are, however, three couplings, κ, γ and Λ, in the action evaluated on a time-independent background. Since we have the freedom of choosing a normalization condition to fix the field redefinition, not all these three couplings are independent. By an appropriate choice of the normalization condition, we can at least separate the flow of one coupling constant. Again, we would like to adapt a normalization condition to the spatial curvature term and extract the beta function for the cosmological constant.
Instead of using (3.95), let us first take C to be thereby turning the effective action (3.94) into (3.102) Note that dim(γ/κ 4 ) = 1 by (3.99), which motivates us to take a simple choice of the normalization condition by fixing γ/κ 4 to be constant at all scales. Then, and the gauge-independent effective action becomes In bare perturbation theory, we require that the physical couplings γ ph , κ ph and Λ ph satisfy Therefore, the beta function for γ/κ 4 vanishes, while the anomalous dimension for the cosmological constant is (3.106)

JHEP06(2017)004
For γ > 0 and λ > 1 or λ < 1 2 , δ Λ is real and positive. It is interesting to note that when λ = 1 2 , which is required for Weyl symmetry, δ Λ vanishes at one-loop order. When λ approaches 1, which is required for Lorentz symmetry to be realized, the one-loop expression for δ Λ blows up, reflecting the strong coupling problem of the λ → 1 limit [17]. Of course, we are still far from determining if the theory is asymptotically free. One will have to evaluate the heat kernel for time-dependent background geometries to map out the full RG structure.
As a final comment, we note that there is no logarithmically divergent contribution to the coupling in front of the term dt d 2 x N g R. (3.107) This can be seen as follows. Since the UV properties are controlled by the terms with the most derivatives, we can view Λ purely as a coupling constant and expand in a power series of Λ. Since ρ does not contribute to the differential operator O phys , Λ is the only dimensionful parameter that can arise in the one-loop divergence. The contribution of lowest dimension, linear in Λ, has dimension two, and so cannot appear in the coefficient for R. Hence, (3.107) cannot appear at all in the logarithmic divergence at one loop.

Discussion
This paper dealt with the computation of quantum corrections in the simplest version of critical Hořava gravity, the z = 2 projectable theory in 2 + 1 dimensions. Working in a gauge with two free parameters, we computed the quantum effective action in two different cases. The first was flat space with Λ = 0; this is an off-shell background, and we saw that the naïve result was gauge dependent. This gauge dependence is however ephemeral: the effective action in gauge theory can be gauge-dependent, provided the gauge dependence can be eliminated by a field redefinition. On the other hand, for an on-shell background field an infinitesimal field redefinition leaves the value of the action invariant (since the action is stationary under any variation), and therefore the result (if correct) must be gauge independent. Working on the timeindependent on-shell background R × S 2 or R × H 2 with γR 2 = 2Λ, we find a gaugeindependent effective action, as expected. Using this action, we are able to extract one of the one-loop beta functions.
The main result of our paper is therefore equation (3.106), which captures the flow of the cosmological constant Λ at one loop order in z = 2 Hořava gravity in 2 + 1 dimensions, as defined relative to a metric normalization such that γ/κ 4 is constant at all scales.
We focused on the flow of this variable for several reasons, which are all rooted in the fact that our computation is based on the effective action for on-shell, time-independent backgrounds. Working on-shell has several advantages, notably the automatic gauge invariance of the quantum effective action. We furthermore saw an explicit reduction of the partition function to only the physical degree of freedom in the one-loop partition function. This simplification can be traced to the on-shell condition. In this way, the computation of the on-shell effective action could be reduced to the functional determinant of a single scalar operator.

JHEP06(2017)004
Time independence had the further virtue of allowing us to reduce our computations to known properties of the heat kernels of higher order relativistic differential operators. And as a background field computation, of course, this can all be done using only the divergences in a single "vacuum bubble" diagram, without having to compute vertices explicitly.
Towards the full β function. One pays a price for working on time-independent backgrounds, however: divergences in the effective action proportional to K ij are invisible. This means that out of the four couplings 10 of the model -λ, κ, γ and Λ -that played a role here, we can only determine the flow of one. (Note that not all of these coefficients are physically meaningful. For example, in the text we rescaled g ij to make one coupling take a value of our choosing.) In order to compute the remaining beta functions, one must relax one of these restrictions. The full computation can in principle be done entirely on-shell, provided we allow time-dependent backgrounds. This approach runs into one of two possible difficulties. The first is that of finding explicit classical backgrounds on which to work. The simplest backgrounds are cosmological backgrounds of FLRW type, in which case K ij is pure trace. Imposing the trace constraint reduces the number of beta functions that can be computed by one; to obtain the complete flow of the theory would still require backgrounds on which K ij is not pure trace.
If we accept this limitation, we run into the second complication, that in pure Hořava gravity such backgrounds are de Sitter-like. As a result they suffer from large contributions to the effective action from temporal boundaries (the boundary area grows at about the same rate as the bulk volume), which makes it difficult to distinguish the boundary and bulk contributions to the effective action.
Even after overcoming these difficulties there remains a potentially troublesome point. Our methods expressed the determinant in the (h, σ) sector as a product, after which we cancel against O g coming from the ghost sector. This requires the product identity det(AB) = (det A)(det B), but this identity runs into difficulties in the infinitedimensional case. These can be surmounted straightforwardly when [A, B] = 0 (as was the case for us), but it is more problematic when [A, B] = 0, as occurs in the time-dependent case, and leads to ambiguities in the result. (For one discussion of this issue, see [18].) These problems point to a general need for more flexible methods to compute loop effects in Hořava gravity. In the end, it may turn out that the only viable method is to work on perturbative backgrounds, performing explicit expansions of the heat kernel of a matrix differential operator.
Generalization to non-projectable and conformal gravity. For many purposes, the most interesting class of Hořava gravities are the non-projectable theories, which relax the constraint ∇ i N = 0 and allow N = N (t, x) to depend on space. For example, in phenomenological applications the non-projectable variant requires much less fine tuning JHEP06(2017)004 to be consistent with observational constraints [19,20]. From a more conceptual point of view, the "conformal" variants -those invariant under anisotropic Weyl symmetry [1] are also of considerable interest. We here briefly summarize the extension of our methods to these models, and discuss some of the new challenges that arise.
The novelty arising in the non-projectable theory is that once N has local fluctuations, it gives rise to a new constraint. Because the number of additional constraints equals the number of additional fields (one in both cases) the number of propagating degrees of freedom remains unchanged, but the details of the spectrum and the gravitational interaction are modified.
In the computation of the one-loop effective action, the non-projectability leads to two new features that should be handled carefully. The first is that N cannot be set to 1 by a gauge transformation, and therefore needs to be incorporated appropriately into the gauge-fixing conditions. The second is that the second-class constraint is non-linear, and so its measure needs to be defined carefully. The question of whether the right approach is to solve directly for the Dirac bracket, or to use the ghost formalism of [21], or whether there exists a simple prescription giving the correct contributions to the path integral, we leave for future work. Now for the conformal case. For certain choices of parameters in the gravitational action, an additional local symmetry arises: anisotropic Weyl invariance. This is a symmetry under a Weyl scaling where Ω = Ω(t, x) is an arbitrary function. In this case, at the classical level the secondclass constraint of N is replaced by a first-class constraint, which eliminates the scalar degree of freedom entirely. The question of whether this symmetry can survive at the quantum level is of considerable interest, particularly in 2 + 1 dimensions, where conformal Hořava gravity has no propagating degrees of freedom and therefore provides a useful analog of three dimensional Einstein gravity, with its importance in addressing the conceptual issues of quantum gravity. In some ways, the conformal case bears similarities to the projectable theory, in that we can gauge fix N = 1 if we like. On the other hand, to answer questions about the preservation of conformal symmetry, it is important to choose a gauge-fixing condition that is invariant under background Weyl transformations. 11 In particular, if we want to study whether Weyl symmetry is anomalous, we should not gauge-fix N = 1, and instead work in a more general background gauge. This requires us to modify the gauge-fixing conditions.
One important difference in the conformal case is that to preserve background Weyl symmetry, the gauge fixing must respect z = 2 scaling. The type of gauge fixing used here and in [7] makes this possible. It is this consideration that initially led us to the gauge-fixing used in this paper. We note that background Weyl invariance requires some new features in the gauge-fixing condition, in particular in that N and n must be included to construct an appropriate Weyl-invariant object.

JHEP06(2017)004
Beyond its interest as a toy model, the study of the conformal theory is relevant to the problem of quantum membranes [1]. The path integral for relativistic quantum membranes is not renormalizable, putting a theory of fundamental relativistic quantum membranes out of reach. This is reflected in the Polyakov action formalism in the nonrenormalizability of three-dimensional gravity. With z = 2 scaling, on the other hand, the Polyakov action becomes power-counting renormalizable. In this picture, the critical membrane theory would become conformal Hořava gravity coupled to a z = 2 non-linear sigma model. The crucial question of whether such critical membrane theories exist, or whether a Weyl anomaly spoils criticality, we leave to future research.

JHEP06(2017)004
The contribution of the gauge fields to the partition function is therefore equal to where we recall that the (det D) −1/2 piece comes from integrating out the auxiliary field Φ. The operator S (2) is given in (2.6) and S (2) + S g.f. reads Here we see explicitly the virtue of the choice D = −∂ 2 + v 2 , or U = 1: and Combining this with (A.1) gives the total partition function To calculate det S (2) + S g.f. for general D, we write S (2) + S g.f. in ADM form, Comparing (A.12) with (A.7) immediately gives the "spatial metric" the inverse of which is given by The "shift" variables with lower and upper indices are

B Physical modes without gauge-fixing
The case of flat space is sufficiently simple that we may actually bypass the gauge-fixing procedure in either the U(1) gauge theory or the Hořava gravity theory and still determine the physical modes and dispersion relations.
In the case of U(1) gauge theory, we start with the action (2.5), which has not yet been gauge-fixed. This is written as At this point A 0 and A i have different dimensions. We must rescale one to remedy this. The solution was presented in section 2.1: redefine A 0 by a factor of √ D where D is some spatial differential operator of dimension one, namely some linear combination of ∂ 2 and v 2 . Thus, S (2) becomes While we will use the inspired choice D = −∂ 2 + v 2 , one could use any other linear combination, including just v 2 . The subsequent conclusions will not change, which is consistent with the D-independence shown in appendix A. Fourier transforming S (2) gives For an explicit example, take D = 2, in which case the above matrix is 3 × 3 and can be easily diagonalized. The unnormalized eigenvectors and eigenvalues are Eigenvector Eigenvalue The first is a zero mode and is unphysical. The second gives an unphysical dispersion relation and is thus also an unphysical mode. Only the third mode is physical. This mode propagates with the dispersion relation we expect from the gauge-fixing procedure and we can also see from the eigenvector that it is precisely the one transverse mode in the A i 's.

C Useful formulas
In this appendix, we prove a number of formulas which are useful in expanding out the action of Hořava gravity around curved space. The identities are understood to hold under integration and thus we set all total derivatives to zero. Finally, we take the background to have constant curvature R.
Recall that we decompose the spatial metric fluctuation as where h = g ij h ij and g ij H ij = 0. Furthermore, we decompose H ij as where g ij H ⊥ ij = 0, ∇ j H ⊥ ij = 0 and ∇ i η i = 0. In two dimensions, we set H ⊥ ij to zero and η i = ε ij ∇ j η for some scalar η.

JHEP06(2017)004 D Jacobians
Here we review the computation of a partition function under the change of variables Φ = F Ψ, for some linear differential operator F that is often, but not always, local. The path integral for Φ is defined with respect to a measure on the space of field configurations, which in the background field method we take to be covariant with respect to background diffeomorphisms. It is natural to define the measure in terms of a metric G Φ on this space. For example, if Φ i has a spatial index, the inner product of two infinitesimal variations δΦ (1) and δΦ (2) takes the form G Φ (δΦ (1) , δΦ (2) ) = δΦ (1) , δΦ (2) Φ = dt d 2 x g g ij δΦ (1) i δΦ (2) j . (D.1) The path integral can schematically be written 12 The metric is not covariant under a change of variables; instead, the measure transforms as The operator O F is computed by setting The Jacobian is then expressed as J F = √ det O F . As an example, let us consider the Jacobian for the transformation The natural metric for Φ i is the one given above, while that for φ and φ is δφ (1) , δφ (2) φ = dt d 2 x g δφ (1) δφ (2) , (D.7a) δ φ (1) , δ φ (2) φ = dt d 2 x g δ φ (1) δ φ (2) . (D.7b) In the text we are primarily interested in the case where g is time-independent and whose spatial slice is a symmetric space. Then δΦ, δΦ Φ = δφ, − δφ φ + δ φ, − δ φ φ . For another example, let us conside the Jacobian for the transformation defined in (C.2), (D.11) Note that H ij is traceless. A natural metric on the space of traceless tensors is δH (1) , δH (2) H = dt d 2 x g δH (1) ij g ik g j δH (2) k . (D.12) For the scalar modes η and σ we define δη (1) , δη (2) η = dt d 2 x g δη (1) δη (2) , (D.13a) δσ (1) , δσ (2) σ = dt d 2 x g δσ (1) δσ (2) . (D.13b) Then, applying Identity 8, we obtain δH, δH H = δη, 2 + R δη η + δσ, 1 2 + R δσ σ . (D.14) (In general there is an η-σ cross-term involving ∇ i R, but this vanishes on the backgrounds used in this paper.) Therefore, the associated Jacobian is 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.