Cosmological perturbations and non-Gaussianities in Hoˇrava-Lifshitz gravity

We investigate primordial perturbations and non-gaussianities in the Ho\v{r}ava-Lifshitz theory of gravitation. In the UV limit, the scalar perturbation in the Ho\v{r}ava theory is naturally scale-invariant, ignoring the details of the expansion of the Universe. One may thus relax the exponential inflation and the slow-roll conditions for the inflaton field. As a result, it is possible that the primordial non-gaussianities, which are ``slow-roll suppressed'' in the standard scenarios, become large. We calculate the non-gaussianities from the bispectrum of the perturbation and find that the equilateral-type non-gaussianity is of the order of unity, while the local-type non-gaussianity remains small, as in the usual single-field slow-roll inflation model in general relativity. Our result is a new constraint on the the Ho\v{r}ava-Lifshitz gravity.


I. INTRODUCTION
A renormalizable theory of gravity was proposed by Hořava [1][2][3]. This theory reduces to Einstein's general relativity (GR) for large scales, and may be a candidate for the UV completion of general relativity. This theory is renormalizable in the sense that the effective coupling constant is dimensionless in UV. The essential point of this theory is the anisotropic scaling of temporal and spatial coordinates with dynamical critical exponent z, In 3 + 1 spacetime dimension, the Hořava theory has an ultraviolet fixed point with z = 3. Since the Hořava theory is analogue to the scalar-field model studied by Lifshitz, in which the full Lorentz symmetry emerges only at the IR fixed point, the Hořava theory is also called the Hořava-Lifshitz theory. Because of this anisotropic scaling, time plays a privileged role in the Hořava theory. In other words, spacetime has a codimension-one foliation structure, which leaves the foliation hypersurfaces of constant time. Thus, contrary to GR, full diffeomorphism invariance is abandoned, and only a subset (of the form of local Galilean invariance) is kept. More precisely, the theory is invariant under the foliation-preserving diffeomorphism defined by t →t(t) , In the infrared (IR), due to a deformation by lower dimensional operators, the theory flows to z = 1, corresponding to the standard relativistic scale invariance under which the full deffeomorphism, and thus general relativity, is recovered. Non-relativistic scaling allows for many non-trivial scaling theories in dimensions D > 2. The Hořava-Lifshitz theory allows a theory of gravitation that is scale-invariant in UV, while the standard GR with full diffeomorphism emerges at the IR fixed point. The Hořava-Lifshitz gravity theory has been intensively investigated [4][5][6][7][8][9][10][11][12][13][14][15][16][17][18] (see also [19,20] for reviews and more references therein). In particular, cosmology in the Hořva theory has been studied in [7][8][9][10][11][12][13]. Homogeneous vacuum solutions in this theory were obtained in [12], and scalar and tensor perturbations were studied in [8,10]. In [7,9], the cosmological evolution in Hořava gravity with scalar-field was extensively studied, and the matter bounce scenario in the Hořava theory was investigated by Brandenberger [11].
As was pointed out in [7], the Hořava theory has at least two important properties. The first is its UV renormalizability, while the second is more interesting for cosmology. The fact that the speed of light is diverging in UV implies that exponential inflation is not necessary for solving the horizon problem. Moreover, the short distance structure of perturbations in the Hořava-Lifshitz theory is different from the standard inflation in GR. In particular, in the UV limit, the scalar field perturbation is essentially scale-invariant and is insensitive to the expansion rate of the Universe, as has been addressed in [7][8][9][10]. The key point is that the UV renormalizability indicates that the Lagrangian for the non-relativistic scalar field should contain up to six spatial derivatives. Thus, in the UV limit, the dispersion relation is ω 2 ∼ k 6 , which is contrary to ω 2 ∼ k 2 in standard GR. This phenomenon causes different k-dependence of the two-point correlation function and thus the scalar perturbation in Hořava gravity is naturally scale-invariant in the UV limit.
In this paper, we extend the previous works on cosmological perturbation theory in Hořava gravity, including nongaussianities. The 6-parameter ΛCDM model provides an accurate description of the Universe [21]. In particular, the primordial perturbations are assumed to be Gaussian. Deviation from the Gaussian distribution, i.e., primordial non-gaussianity, has not been observed. This puts a strict constraint on any model of the early Universe. The non-Gaussian features of Hořava gravity have not been studied in detail, except in [22] and in [23], which investigated the non-gaussianities of the scalar field and of the gravitational waves in Hořava gravity, respectively. Actually, one of the essential differences of Hořava gravity from Einstein's general relativity is that it contains quadratic curvature terms in the theory. Moreover, the foliation-preserving diffeomorphism does not allow, unfortunately, to choose a spatially-flat gauge as in GR. Thus, in general, the perturbation theory in Hořava gravity in quite involved. On the other hand, the number of dynamical degrees of freedom in the spatial metric is 3 in Hořava gravity (contrary to 2 in GR), with 2 tensor degrees of freedom, as usual, and an additional scalar dynamical degree of freedom (see also [24][25][26] for the general framework of spatially covariant theories of gravity).
In this work, we focus on the perturbation of the scalar field. Thus, for simplicity, we neglect the spatial metric perturbation. We introduce a scalar field, following the strategy in [7,9]. We pay special attention to the nongaussianities in this scalar field model in Hořava gravity. The basic idea is that, as has been addressed before, the divergence of the speed of light and the scale-invariance of the scalar perturbation in Hořava gravity indicate that there is no need to assume an exponential expansion of the Universe. Moreover, the traditional slow-roll conditions are not necessary. While it is well-known that in slow-roll inflationary models non-gaussianity is suppressed by slow-roll parameters [27], and thus too small to be detected (see e.g. [28] for a review of non-gaussianities in cosmological perturbations. Various models have been investigated to generate large non-gaussianities by introducing more complicated kinetic terms [29][30][31][32] or more fields [33][34][35][36][37][38]). However, in Hořava gravity, there are no slow-roll conditions, and thus the "slow-roll suppressed" non-gaussianities can become large. In this work, we focus on the non-gaussianity from the bispectrum, which is defined from the three-point correlation function of the perturbation. We find that the equilateral-type non-gaussianity is roughly of the order of unity, while the local-type non-gaussianity remains very small, as in the usual single-field slow-roll inflation in GR.
The paper is organized as follows. In Section 2, we briefly review the Hořava gravity and set out our conventions. In Section 3, we couple the scalar field to Hořava gravity, and describe the cosmological evolution of the Hořva gravity/scalar matter system. In Section 4, we calculate the scalar field perturbation, including the gravity perturbations. We get the full second-order perturbation action, which reduces to the standard result in the IR limit. In the UV limit, the scalar perturbation is essentially scale-invariant. In Section 5, we calculate the non-gaussianities. Finally, we make a conclusion and discuss several related issues.

II. BRIEF REVIEW OF HOŘAVA-LIFSHITZ GRAVITY
In this section we briefly review the Hořava-Lifshitz theory [1]. The dynamical variables in Hořava-Lifshitz gravity are the spatial scalar N , spatial vector N i and spatial metric g ij . This is similar to the ADM formalism of the metric in standard general relativity, while in Hořava gravity, N , N i and g ij are related to the space-time metric as where c is the speed of light. The action of Hořava-Lifshitz gravity contains a "kinetic" part and a "potential" part, with where is the extrinsic curvature and K = g ij K ij . The potential terms are given in the "detailed-balance" form where C ij is the Cotton tentsor defined by Note that in (5) λ is a dimensionless coupling of the theory, and therefore runs. As mentioned in the introduction, the essential point of the Hořava theory is the anisotropic scaling of temporal and spatial coordinates: t → ℓ z t and x i → ℓx i . The classical scaling dimensions of various quantities in the Hořava theory are summarized in Tab.I.

III. COSMOLOGY WITH SCALAR FIELD MATTER
A. Coupling the scalar field to Hořava-Lifshitz gravity In this work, we couple the scalar-field matter with Hořava-Lifshitz theory following the strategy in [7,9]. The general structure of the action of the scalar-field matter and Hořava-Lifshitz gravity contains two parts: a quadratic kinetic term with the foliation-preserving diffeomorphisms and a potential term: The "potential term" is where ξ i and their properties in UV/IR are summarized in Tab.II. In Tab.II,  and g 1 , g 11 etc. can be constant, or in general they can be functions of φ. We assume g 3 > 0 in order to guarantee the stability of the perturbation in UV.

B. Equations of motion
The full equations of motion for N , N i and g ij have been derived in [7,9,12]. For our purpose we focus on the equations of motion for N and N i , which we write here for later convenience :

C. Cosmological Evolution
We now consider the cosmological background evolution in Hořava-Lifshitz gravity. We assume the background to be homogeneous and isotropic, and use the residual invariance under time re-parametrization to set N = 1. We focus on the flat 3-dimensional case. The background values are In this background, the action in the gravity sector is significantly simplified. In particular, R ij = C ij = 0 and the spatial covariant derivatives mostly vanish. The equation of motion for N to the 0-th order gives where we have denoted and H ≡ȧ/a is the familiar Hubble parameter. The equation of motion for g ij gives The equation of motion for the scalar field isφ

IV. COSMOLOGICAL PERTURBATION
We now consider cosmological perturbation in Hořava gravity coupled to scalar-field matter. As has been addressed before, the action of Hořava gravity is complicated due to the quadratic terms in spatial curvature. We recall that in GR, one can choose various gauges to simplify the calculations. However, the case is different in Hořava gravity (see Appendix A for a discussion of gauge transformation and gauge choice in Hořava theory). Since the "foliation-preserving" diffeomorphism is only a subset of the full diffeomorphism in GR, one may expect in general that there are less gauge modes and more physical modes 1 . Thus, the perturbation theory in Hořava-Lifshitz theory is very involved, but also interesting.
We consider the scalar-field perturbation by neglecting the spatial metric perturbation. We assume that the background scalar field is homogeneous φ 0 = φ 0 (t).
As we focus on the scalar perturbation, we write (in the background N = 1 gauge) 2 Here, the subscript "1" denotes the first-order in δφ ≡ δφ. The constraints in Eq. (11) become Solving equations (18) up to the first-order of δφ, we get 3 and θ 1i = 0, as usual in GR. Note that when λ = 1, α 1 also depends onQ, which is different from GR, where α 1 ∼ δφ.
It is useful to note that in the UV limit α 1 and β 1 become where we have used the background equations of motion in the UV limit.

A. Linear perturbation
As we are neglecting the spatial metric perturbation, the gravity sector is greatly simplified After a rather tedious but straightforward calculation, we get the quadratic part of the action for the scalar perturbation Q: where In (22), g i ≡ g i (φ 0 ). In deriving Eq. (22), we have used the background equations of motion. Moreover, no approximations were made in deriving Eq. (22) and thus it is exact. We can use Eq. (22) to analyze the behavior of perturbations both in the IR or UV limits and in the interpolation era.

IR limit
Taking the IR limit of the full second-order perturbation action Eq. (22) and choosing λ = 1, we get where we have defined the dimensionless parameters ǫ =φ 2 2H 2 , If we further set α = 1 2 and choose g 1 = − 1 2 , Eq. (24) reduces to the familiar result in the perturbation theory in GR. Especially, the perturbation Q is scale-invariant when the expansion of the Universe is exponential, and thus with an approximately constant Hubble parameter H ≈ const.

B. Scale-invariant spectrum in Hořava-Lifshitz era
We now focus on the behavior of the perturbation theory in the UV-limit, whereQ 2 and Q∆ 3 Q terms dominate. The perturbation action Eq. (22) becomes rather simple in the UV limit, Note that γ is now constant. It is convenient to use a new variable u, defined as u ≡ a √ γQ. After changing into conformal time η, defined by dt = adη, and going into Fourier space, the second-order perturbation action reads The equation of motion for the perturbation reads Here, we assume for simplicity that g 3 is approximately constant. The mode function is The mode function is chosen such that it satisfies the Wronskian normalization condition: Moreover, the short-time behavior of the mode function Eq. (30) is analogue to that of a positive-frequency oscillator. The tree-level two-point correlation function of Q(k, η) is and thus the power spectrum of Q is given by with The so-called dimensionless power spectrum of Q is The power spectrum of the scalar perturbation is naturally scale-invariant in the UV limit, ignoring the details of the expansion of the Universe. This feature is contrary to that in GR, where a nearly constant Hubble expansion rate H is needed to guarantee the scale-invariance of the spectrum. Due to this fact, there is no need to take any "slow-roll"-type conditions for the scalar field. We would like to make some comments here. The crucial picture of the standard inflation is that quantum fluctuations are generated in the subhorizon region (k ≫ aH), and are then stretched to the cosmological size and become classical (k ≪ aH). The horizon-exiting point corresponds to k = aH. Thus, the "horizon-exiting" process exists only when aH is an increasing function of time. If we assume a power law inflation a ∝ t p , it requires p > 1. This is violated by the curvature, and is only satisfied when the equation of state w < −1/3. This is why in the standard inflation model we need a slow-rolling scalar field to mimic the cosmological constant and to drive an exponentially expanding background. However, the "horizon-exiting" process occurs generically in the Hořava-Lifshitz era for rather general cosmological backgrounds. From Eq. (29), it is obvious that the perturbation stops oscillating when k 6 a 6 ∼ H 2 , and thus requires that a 6 H 2 is an increasing function of time. Obviously, if we assume a power law expansion a ∝ t p , this demands p > 1/3. In terms of the comoving time a ∝ |η| p , we need p < 1/2. This condition can be satisfied by any matter component with the equation of state w < 1.
We get the scale-invariant power spectrum in the UV limit from the equation of motion Eq. (29). While one may consider the full equation of motion, the second-order action Eq. (22) can be used. In general, the equation of motion for the perturbation has the following functional form: where c 1 , λ 1 and λ 2 are dimensionless parameters, m is the effective mass parameter, and ℓ is the length scale of the whole theory. The functional form of the dispersion relation in (36) has been intensively studied in the investigation of trans-Planckian effects [39][40][41][42], and also in statistical anisotropy [43]. Although complicated, it is interesting and important to investigate Eq. (36) in order to understand the behavior of the perturbation not only in the UV limit, but also in the interpolation region between UV and IR.

V. NON-GAUSSIANITIES
In this section, we investigate the non-gaussianities, which characterize the interaction of the perturbations.

A. Bispectrum
We focus on the third-order perturbation action and the three-point correlation function of the perturbation Q. The third-order action in the gravity sector is and in the scalar field sector is where α 1 and β 1 are given in Eq. (19). We are interested in the UV behavior of the perturbation. In the UV limit, after a straightforward calculation, the third-order perturbation action reads, where which are dimensionless constants (recall that we assume g 3 to be approximately constant). In Eq. (39), the formal operator ∂i∂j ∂ 2 should be understood in momentum space. After changing into comoving time η and into Fourier space, we have where we denote k 123 ≡ k 1 + k 2 + k 3 . The three-point correlation function in cosmological context is evaluated in the so-called "in-in" formalism where η * is the time when perturbation modes exit the sound horizon, and H is the Hamiltonian which can be read from Eq. (41) by noting that in the third-order H (3) = −L (3) . Thus, for three-point interactions described by Eq. (41), we have where we have introduced the "shape factor" S(k 1 , k 2 , k 3 ) defined as For power law expansion a(η) ∝ |η| p (p = 0), the time-integral in Eq. (43) can be evaluated exactly, and the three-point correlation function reads, with which is the so-called bispectrum.

B. Non-linear parameter f NL
In practice, it is convenient to introduce non-linear parameters to characterize the non-gaussianities. The dimensionless non-linear parameter f NL from the three-point correlation function is defined as where the power spectrum P (k) is given in Eq. (33). Note that although dimensionless, f NL is in general k-dependent.
A straightforward calculation gives In the equilateral limit (k 1 ≈ k 2 ≈ k 3 ), while in the squeezed limit ( Thus, in the UV limit, we find that the equilateral-type non-gaussianity is roughly ∼ O(1), while the local-type non-gaussianity is very small. This is not surprising, since in Hořava-Lifshitz gravity the scalar-field perturbation in the UV limit is naturally scale-invariant. Thus, no slow-roll condition is needed to guarantee the exponential expansion of the Universe. On the other hand, it is well known that in standard slow-roll inflationary models in GR non-gaussianities are suppressed by slow-roll parameters [27]. However, in Hořava gravity, no slow-roll parameters are needed. Thus, one may expect the slow-roll suppressed non-gaussianity to be of the order of unity. However, in our simplest scalar-field model with the action given by Eq. (8) and Eq. (9), the temporal kinetic term is canonical, and in general there is no enhancement of the non-gaussianity by non-canonical kinetic terms as in the K-inflation or DBI-inflation models [29][30][31][32]. Moreover, the scalar-field action Eq. (8) is mostly "derivative-coupled", thus the local-type non-gaussianity (which characterizes the local couplings of the perturbations in the real space) is small, as expected.

VI. CONCLUSION
In this work, we investigated the cosmological perturbation theory in Hořava-Lifshitz gravity and the nongaussianities from the bispectrum. The most interesting feature of Hořava gravity is that in the UV limit the scalar perturbation is essentially scale-invariant, ignoring the details of the expansion of the Universe. Moreover, together with the fact that the speed of light in the UV limit diverges, there is no need to assume exponential expansion of the early Universe, or the usual scalar-field driven slow-roll inflation. In particular, the slow-roll conditions are not necessary. Thus, one may expect that in the absence of slow-roll conditions, the non-gaussianities can become large. We calculated the three-point correlation function of the scalar perturbation and found that the equilateral-type non-gaussianities are of the order of unity due to the absence of the slow-roll-type conditions, while the local-type non-gaussianities remain small, as in the usual single field inflation in GR.
We focused on the scalar-field perturbation in the Hořava-Lifshitz theory, neglecting the spatial metric perturbations. However, the latter is obviously the most important and interesting part in the Hořava theory. Since in the Hořava theory the dynamical degrees of freedom in the spatial metric part are 3, especially, there is one additional scalar degree of freedom. It is important to investigate the property of this additional degree of freedom. Moreover, in this work, we only investigated the behavior of the perturbation in the UV limit at z = 3, while it is interesting to study the full equations of motion, especially in the interpolating region between UV and IR. We considered the scalar field with canonical temporal kinetic term. However, one could expect enhancement of non-gaussianities if more general kinetic terms are considered. The Hořava-Lifshitz theory, although originating from a renormalizable quantum gravity in 4 dimensions, may be a potential competitor to the standard inflation theory and deserves further investigation.
Due to the fact that ξ 0 is a function of η only, the gauge choice is fairly restricted. In particular, (or unfortunately), two familiar gauge-choices -the "longitudinal gauge" and "spatially-flat gauge" -are not allowed in the Hořava theory. We can use χ = χ(η, x i ) freely to set E = 0, while in general we cannot use ξ 0 to set B = 0 and get the longitudinal gauge, or to set ψ = 0 and get the spatially-flat gauge.
However, there is still a possible gauge choice in the Hořava theory. First, as we have mentioned, we can choose χ = a 2 E to getẼ = 0. This leaves the question: are we able to choose another gauge condition to set one of φ, ψ and B to 0? In Hořava's original formulation of the theory, N = N (t) is assumed to be a function of time only. In this case we can choose the proper value of ξ 0 to set the fluctuation of N to zero, i.e. φ = 0. Thus, after gauge transformation with ξ 0 = 1 a dη aφ(η) + c , we getφ = 0. This does not determine time-slicing unambiguously, but we are left with time reparametrization. If we relax the restriction that N has to be a function of time only, then there is no gauge condition we can choose. In this case, we are left with 3 non-vanishing variables φ, ψ and B.
The two Bardeen potentials are still gauge-invariant variables in Hořava gravity. This is also because the foliation-preserving diffeormorphism is a subset of the full symmetry in GR. Note that there is an infinite number of gauge-invariant variables, for example, combining Φ and Ψ gives another useful gauge-invariant variablẽ For the "space-time" scalar field φ (now "scalar" means invariant under "foliation-preserving" diffeomorphism), if we assume that the background value is homogeneous φ 0 = φ 0 (η), the gauge transformation for the scalar fluctuation is as usual The gauge-invariant variable for δφ is as in GR: