{'fraction_non_alphanumeric': 0.030063929044519957, 'fraction_numerical': 0.011288371825145424, 'mean_word_length': 4.409190031152648, 'pattern_counts': {'":': 0, '<': 8, '<?xml version=': 0, '>': 1, 'https://': 0, 'lorem ipsum': 0, 'www.': 0, 'xml': 0}, 'pii_count': 0, 'substrings_counts': 0, 'word_list_counts': {'cursed_substrings.json': 5, 'profanity_word_list.json': 0, 'sexual_word_list.json': 0, 'zh_pornsignals.json': 0}}
An effective field theory of holographic dark energy

A general covariant local field theory of the holographic dark energy model is presented. It turns out the low energy effective theory of the holographic dark energy is the massive gravity theory whose graviton has 3 polarisations, including one scalar mode and two tensor modes. The Compton wavelength is the size of the future event horizon of the universe. The UV-IR correspondence in the holographic dark energy model stems from the scalar graviton's strong coupling at the energy scale that marks the breaking down of the effective field theory.


Introduction
The holographic dark energy model [1] [2] is based on the simple idea that the vacuum energy density arising from the quantum fluctuation of the UV-cut-off quantum field theory should relate to the boundary surface of a system in the way where M 2 p = 8πG is the reduced Planck mass, L is the size of a system (thus L 2 is essentially the area of the boundary surface). This relation between the energy density in a bulk and the area of its space-time boundary surface is rooted in the holographic principle [6]. The argument that leads to the above relation is the following. The zero-point energy diverges quartically and thus the energy density scales as ρ Λ ∼ Λ 4 if we simply cut-off the divergence at the UV scale Λ 1 . However, this simple scaling is violated when the total energy in the system with size L, i.e. L 3 Λ 4 , approaches to the mass of the black hole in the same size LM 2 p . In fact, the total energy must be bounded by the mass of the black hole from above, as the effective quantum field theory breaks down at the Schwarzschild radius scale. This energy bound is stronger than the Bekenstein entropy bound [7] which inspired the proposal of the holographic principle [6]. Nowadays it is widely believed that the holographic principle is one of the most important cornerstones of quantum gravity.
Assuming that the energy bound is saturated on the cosmological background, one obtains the important UV-IR correspondence Λ ∼ M p /L (1.2) and the vacuum energy ρ Λ comparable to our current critical energy density. The related idea was discussed in Ref. [8] [9] with the Hubble radius as the IR cut-off. However, as pointed out in Ref. [10], the resultant equation of state is greater than −1/3 and thus the vacuum energy fails to accelerate the cosmic expansion, if we simply take the Hubble radius L = 1/H as the IR cut-off. Later in the same year, by adopting the future event horizon as the IR cut-off, an accelerated expanding solution in the Friedman equation was eventually obtained in Ref. [2]. Interestingly, the cosmic coincidence problem can also be resolved by inflation in this scenario, provided the minimal number of inflationary e-foldings. Since then, the holographic dark energy has drawn a lot of attention and has been widely studied, see Ref. [11] for a comprehensive review on the topic. In passing, I shall mention that it was pointed out very recently that the holographic dark energy model may alleviate the Hubble tension problem [12]. Given the phenomenological success, however, one of key pieces for our holographic jigsaw puzzle is still missing. Namely we do not know how to write down the general covariant action for the holographic dark energy model. An attempt was made in the unpublished work [13] in which a mini superspace action was given. Nevertheless, due to the absence of the general covariant action, one may easily spot several conceptual problems. For instance, the evolution of our universe in the past is dependent on the one in the future. This apparent causality violation may not imply the pathology of the model, but rather that the low energy effective field theory is missing. This causality violation may be partially addressed at the cosmological background evolution level, given the mini superspace action [13]. However, it is still not quite clear whether the local physics violates the causality, as the perturbation theory is also missing 2 .
In the current work, I aim at finding this important missing piece for our holographic jigsaw puzzle. As we will see in the remain of this paper, the low energy effective field theory (EFT) for the holographic dark energy model is actually a massive gravity theory. A graviton has three polarisations, namely one scalar mode, and two tensor modes. The UV-IR correspondence eq. (1.2) stems from the strong coupling of the scalar mode above the energy scale where our effective field theory breaks down. One may get puzzled at this point as the Poincare symmetry in 4 dimensional space-time implies a massive spin-2 particle has 5 polarisations (including helicity 0, ±1, ±2). However, the massive gravity theory that we are going to engage is different from the conventional one [15] [16], as the global Lorentz invariance is broken on the cosmological background and thus a graviton does not necessary to be well equipped with 5 polarisations [17] [18]. See Ref. [19] [20] for two comprehensive reviews on the topic of massive gravity. This paper is organised in the following. I will start from the mini superspace action in the section 2, and then write down the covariant action by adopting the Stueckelberg trick. Upon writing down the general covariant action, an Ostrogradsky ghost is spotted in the decoupling limit. The ghost is eliminated in the section 3, in which the linear perturbation and Hamiltonian structure are analyzed in detail. The conclusive remarks and outlooks are given in the section 4.

From mini superspace to general covariance
We take the flat FLRW ansatz ds 2 = −N 2 dt 2 + a 2 dx 2 and the following mini superspace action [13] as our starting point, where M 2 p ≡ 8πG is the reduced Planck mass, R is the 4 dimensional Ricci scalar, and L is the variable with length dimension and subject to the constraint equation enforced by the Lagrangian multiplier λ, i.e.L = −N/a. We may integrate this constraint equation from the infinite past −∞ to nowadays, namely where L(−∞) is the initial condition in the infinite past, which requires the input from new physics (such as quantum gravity) to fully determine its value, in light of the Hawking-Penrose's singularity theorem [21]. Therefore L(−∞) remains unknown to me due to my ignorance and obtuseness. On the other hand, one may try to take the integration from the other side, namely given the asymptotic solution in the infinite future t = +∞. It turns out L(+∞) = 0 for the asymptotic solution derived in Ref. [13] and later is confirmed by the numerical computation in Ref. [22]. Therefore, R h ≡ aL is exactly the size of the future event horizon 3 . It follows the energy density of the dark sector where the first term is the holographic term, and the second term is the dark radiation as it scales as a −4 . In passing, I shall mention that in some cases, for instance, a cyclic universe realised by means of some exotic matter content, the future event horizon may not exist as the integral in the eq. (2.3) may not converge to a finite and definite value. The effective field theory of the holographic dark energy, which will be developed in the remains of the this paper, does not apply to these cases. Therefore, we will only focus on the universe models whose future event horizon has finite, definite and non-vanishing size. 3 Mathematically a more general solution toL = −N/a is We may input the value of L at the moment t 0 in the past by hand, and then treat it as the initial condition of the equations of motion, which is equivalent to introducing an assumption about the UV completion of the theory.
The general covariance can be recovered by Stueckelberging the mini superspace action [23]. We may follow the following dictionary, where δ ab is the metric of the Stueckelberg scalar field space, δ ia is the pullback mapping between the physical space-time and the scalar field space, and i, j, k... are adopted as 3 dimensional spatial coordinate indices, while a, b, c... are adopted as the field space indices. These 4 Stueckelberg fields are not yet canonically normalized and thus they are of the length dimension. The field space is flat, and respects the SO(3) rotational invariance. In passing, I shall mention that the global Lorentz invariance is broken in this scalar configuration. Nevertheless, one should not be worried about this Lorentz-violating vacuum expectation values (VEVs) of scalars as the Lorentz invariance is broken anyway on the cosmological background. For instance, CMB is not invariant under a Lorentz boost. Before writing down the general covariant action, which is pretty easy and straightforward at this point, I shall remind readers the hierarchical structure of the theory. We have the Planck scale as the genetic fundamental scale in the first place, and a secondary fundamental scale for our EFT of the dark sector, which is expected to be generated via the non-trivial VEV of the time-like Stueckelberg field ϕ, namely Λ ∼ M p /ϕ. The low energy effective field theory that I begin with is the following, where c is a constant, λ is a Lagrangian multiplier, and S m is the action of the matter sector which minimally couples to gravity. The theory respects the SO(3) rotational invariance, Z 2 symmetry, and the rescaling invariance where is a constant. The Lagrangian multiplier λ generates a constraint between the temporal Stueckelberg field and the spatial Stueckelberg fields. A similar example of this sort can be found in the generalized unimodular gravity [24]. The hierarchical structure of the theory is manifest at the action level. One would expect that the hierarchy disappears if we trace the cosmic evolution all the way backward in time to the infinite past, and it was generated later as the universe expands. The physical significance of the energy scale Λ ∼ M p /ϕ will be shown in the next section. I adopt the FLRW ansatz, and the four Stueckelberg scalar fields take the space-time VEVs in the eq. (2.5). The background equations of motion read, where ρ m and p m are the energy density and pressure of the matter sector, the lapse has been absorbed into the redefinition of time, namely Ndt → dt, and I have rescaled the Lagrangian multiplier λ → λ · ϕ 2 4a 2 M 2 p for convenience. All equations of motion presented in Ref. [13] have been reproduced 4 . Some detailed analyses about this set of equations, including the analytical one and the numerical one, have already been conducted in Ref. [13] and Ref. [22], and I shall not repeat it in my current work. Before moving on to study the local dynamics, I shall briefly comment on the causality. The apparent causality violation stems from our ignorance of the UV completion of the theory, which compels us to integrate the differential equation (2.9) backward from t = ∞ to the moment of interest, as the asymptotic value of ϕ at the infinite future is known to be ϕ(∞) = 0 [13] [22]. Assuming we had enough knowledge about the quantum gravity and we were able to specify the initial condition at t = −∞, we might just integrate the differential equation from t = −∞ to the moment of interest, as we normally do, and there is no apparent causality violation in this setup. The result obtained from this setup should be equivalent to the one integrated backward from t = ∞. In other words, the apparent causality violation appears as an artifactual and mathematical treatment for our equations of motion, rather than a real phenomenon.
It is quite remarkable to see that the holographic dark energy, which is seemingly non local and causality-violating, actually originates from a simple and well defined local field theory. The trick of the game is the spontaneous symmetry breaking, namely that I start from the general covariant theory, then 4 scalar fields take the non-trivial space-time VEVs and break the spatial and temporal diffeomorphism invariance around this vacuum solution. According the Goldstone theorem, there is a massless boson for each generator of the symmetry that is broken. I shall start to introduce the Goldstone excitations around this asymmetric state, where π a is the space-like Nambu-Goldstone boson, and π 0 is the time-like one. Due to the SO(3) rotational symmetry, we can decompose the helicity π a = δ i a ∂ i π +π a , where π is the longitudinal mode, andπ a are two transverse modes satisfying the transverse condition In a gauge field theory, the Goldstone pions decouple from the gauge bosons at the scale much shorter than the Compton wavelength. The physics becomes very transparent in the decoupling limit. A similar decoupling limit can also be adopted, where the Goldstone bosons decouple from the graviton, which is regarded as a gauge boson in this case. The limit is defined by setting M 2 p → ∞ while keeping Λ = M p /ϕ fixed (I neglect the time dependence of ϕ for the time being for the schematic analysis). It turns out this decoupling limit is a very illuminating perspective allowing us to take a quick peek at the microscopic dynamics of the dark sector. Perturbatively expanding the action eq. (2.6) up to quadratic order in the Goldstone excitations, I get the following Goldstone action in the decoupling limit, Seemingly so far so good, but unfortunately there is a pitfall lies in the Goldstone action. Taking variation w.r.t δλ we get the constraint ∂ 2 π = 3π 0 . Inserting it back into the eq. (2.11), we get a higher order temporal derivative term 5 This term yields to a fourth order equation of motion, which requires 4 initial conditions to fully determine its evolution. It implies there are 2 scalar degrees of freedom in the system, instead of one. One of them is actually the infamous Ostrogradsky ghost [25] [26]. The ghost instability spoils the validity of the low energy effective field theory, and renders the theory inconsistent. I will devote the next section to kill the ghost.

Stueckelberg fields and gauge transformation
I define the linear perturbation of our FLRW metric as follows, where α, β, ψ and E are the scalar perturbations, S i and F i are the vector perturbations satisfying the transverse condition ∂ i S i = ∂ i F i = 0, and γ ij is the tensor perturbation satisfying the transverse and traceless condition γ ii = ∂ i γ ij = 0. Under the linear gauge transformation these four Stueckelberg fields transform accordingly, On the other hand, the vector Z µ defined by [27] transforms in the same manner Therefore, the combination Z i − π i andφZ 0 − π 0 are gauge invariant. It is very convenient to transform to the unitary gauge by simply muting all Goldstone bosons, while keeping all perturbation variables in the eq. (3.1). It can be achieved by adopting a proper diffeomorphism ξ µ (t, x). In this gauge, the Goldstone bosons are eaten by the graviton and the graviton develops a mass gap at the low energy spectrum. The graviton turns into a massive spin-2 particle with at most 5 polarizations (it can be less), if the Ostrogradsky ghost had been eliminated (which will be done soon!).
On the other hand, we can also transform back to the Goldstone bosons' gauge by muting the vector Z µ , which can also be done by adopting a proper diffeomorphism ξ µ . Noted that the way back to the Goldstone bosons' gauge is not unique, we may choose to mute some other perturbation variables in the eq. (3.1), based on their transformations under the diffeomprhism.

Hamiltonian analysis in the unitary gauge
The Hamiltonian analysis was introduced by P. Dirac in the 50s and 60s [28], as a way of counting dynamical degrees of freedom, and quantizing mechanical systems such as gauge theories (for busy/lazy readers see the appendix A in the Ref. [29] for a digested version of the method).
It is more convenient to perform the Hamiltonian analysis in the unitary gauge where π 0 = π i = 0, as all terms introduced by Stueckelberg fields appear only in the potential sector in the action and thus the Legendre transformation can be easily done. In the unitary gauge, the action eq. (2.6) reduces to (I set M 2 p = 1 in this subsection) where Λ 1 = ϕ(t) −2 and Λ 2 (t) =φ(t) 2 , namely ϕ andφ should not be treated as a canonical variable and its velocity, and instead, they should be treated as two time-dependent functions subject to some constraint equations which will be derived later. The reason is that In the unitary gauge, all 4 Stueckelberg fields are eaten by the graviton, and turn themselves into ADM variables, or some space-time functions, and cease to be the independent dynamical degrees of freedom. The Hamiltonian can be obtained by performing the Legendre transformation. To this end, I need to adopt the ADM decomposition, the inverse of the metric reads, The conjugate momenta are defined in the following, The Hamiltonian is obtained by performing the Legendre transformation, where N , i , λ are Lagrangian multipliers, and are 5 primary constraints, and where R is the 3-dimensional Ricci scalar, ∇ i is the covariant derivative compatible with 3dimensional induced metric h ij , and ≈ denotes the "weak equivalence", namely equalities hold on the constraint surface. Caution should be paid to the terms with different types of indices. For instance, in the 3 dimensional hyperspace h ab ≡ h ij δ a i δ b j and N a ≡ N i δ a i are scalars (as their space-time indices are dummy indices), δ a i is a vector, while h ij and N i are tensor and vector respectively. The consistency conditions of these 5 primary constraints give rise to the following 5 secondary constraints, where H 0 is the Hamiltonian constraint, H i is the momentum constraint, and C λ is the constraint introduced by hand at the action level, (3.14) These 5 secondary constraints must be conserved in time, which yields to the following 5 consistency conditions, Whether these 5 consistency conditions generate some tertiary constraints, or they only fix the Lagrangian multipliers in the Hamiltonian, crucially depends on the rank of the matrix (3.16) where φ A is the whole set of constraints that we have so far, i.e.
The commutation relations among these constraints are showed in the following, The rank of the matrix M AB is 10. Therefore the consistency conditions eq. (3.15) only fix the Lagrangian multipliers s, instead of generating new tertiary constraints. Let's collect all of primary and secondary constraints in the total Hamiltonian and treat them on the same footing, where N , i , λ , 0 ,˜ i , λ are Lagrangian multipliers, and we have absorbed the shift vector into the˜ i . The algebra closes here. The rank of the matrix also implies all of these 10 constraints are second class. This is the consequence of that the gauge symmetries in GR, namely the space-time diffeomorphisms are all broken in the unitary gauge. Now let's count the number of degrees. In the phase space we have got 22 degrees in the first place, where 20 degrees are from the 10 independent components of g µν and their conjugate momenta, and 2 degrees are from Lagrangian multiplier λ and its conjugate momentum. On the other hand, those 10 second class constraints φ A remove 10 degrees in the phase space, and the leftover is 22 − 10 = 12 degrees in the phase space, (3.19) which corresponds to 6 degrees of freedom in the physical space-time. The sixth mode is the Ostrogradsky ghost.

Ghost elimination in the Goldstone action in the decoupling limit
The Goldstone action in the decoupling limit is a very convenient perspective for us to hunt and eventually execute the ghost. Taking a closer look at the Goldstone action eq. (2.11), I notice that the ghost arises from the kinetic term of the Goldsteon boson π, and the constraint equating the temporal derivative of π 0 and the gradient of the π, namely Therefore, a direct way to eliminate the higher order temporal derivative term is to eliminate the kinetic term ∂ iπ ∂ iπ , which can be achieved by introducing the symmetry [17] π i (t, x) → π i (t, x) + ξ i (t) , (3.21) where ξ i is an arbitrary function of time. This symmetry prohibits all temporal derivative terms of the Goldstone pion π i , including the ones of the longitudinal mode π and the ones of the transverse modeπ i . At the lowest dimensional operator level, the building block which respects this symmetry is the following, In the unitary gauge where all Goldstone bosons are muted, we have Z ab = 1 3 h ij δ a i δ b j , namely the residual symmetry eq. (3.21) strips the shift off the g ij , and leaves us with only 3 dimensional induced metric h ij . Now let's look at the leftover action, which reads where c and λ are assumed to be positive, to ensure that the energy density of dark energy and dark radiation are positive. Therefore, the kinetic term of the leftover scalar mode π 0 has a wrong sign. Unfortunately after removing the problematic higher order temporal derivative term, the theory still contains a ghost. This ghost arises from the gradient term of the Goldstone boson π, namely − 4c + λ 12 ∂ 2 π∂ 2 π → − 3c + 3λ 4 π 0π0 as ∂ 2 π = 3π 0 . (3.24) I have to flip the sign of the term ∂ 2 π∂ 2 π to cure this ghost pathology. The remedy is the operator introduced in Ref. [30],δ Z ab ≡ Z ab − 3 Z ac Z db δ cd Z cd δ cd , (3.25) which is traceless up to the linear perturbation level, and thus it does not contribute to the background evolution. Additionally I shall stick to the Z 2 symmetry and the global scaling invariance eq. (2.7) to tighten up the structure and reduce the arbitrariness of the theory. At the lowest dimensional operator level, I write down the general covariant and ghost free action which respects the residual symmetry eq. (3.21), the global symmetry eq. (2.7) and the SO(3) rotational symmetry in the field space as follows, where Z ≡ Z ab δ ab , the coefficient 3/8 is chosen for the later convenience and S m is the action of the matter sector which minimally couples to gravity. This is the main result of my current work. Compared with the original action eq. (2.6), in the unitary gauge the g ab , which is produced by the spatial Stueckelberg fields ∂ µ φ a ∂ µ φ b δ ab , is replaced by one of our building blocks h ab defined in the eq. (3.22). Take the note of that h ab differs from g ab only at the perturbation level. On the other hand, the operatorδZ ab only appears at the perturbation level too as it is traceless. Therefore, all modifications are operated at the perturbation level, and do not alter the background evolution which is subject to the set of equations in the eq. (2.8) (2.9). The Goldstone action for the scalar graviton reads with the ghost free condition, This ghost free condition will be reproduced later in the full perturbation analysis in the unitary gauge, where the Goldstone bosons are muted.

Full perturbation analysis in the unitary gauge
We adopt the metric perturbation decomposition in the eq. (3.1). Due to the background SO(3) rotational invariance, the scalar perturbations, the vector perturbations and the tensor perturbations completely decouple from each other at the linear perturbation level. We also adopt the unitary gauge and mute all Goldstone bosons. The massive gravity is an analog of the Higgs mechanism in the particle physics, where the Higgs field resides in an asymmetric state. In the unitary gauge, the Goldstone pions in this symmetry broken phase are eaten by the gauge bosons and consequently the gauge bosons become massive. We expect the same phenomena should also occur in the massive gravity, namely once we mute the Goldstone bosons, the scalar graviton and gravitational waves are massive in the unitary gauge.
I will only focus on the perturbation analysis of the pure gravity, while the cases with matter included will be discussed in the future work.

Scalar perturbation
The quadratic action of the scalar perturbation is obtained after a straightforward computation, (3.29) where L EH is the contribution from the Einstein-Hilbert action, (3.30) and the L mass is the contribution from the dark sector, or in other word the graviton mass term, (3.31) I have absorbed the lapse into the redefinition of time, and rescaled the Lagrangian multiplier λ → λ · ϕ 2 4a 2 M 2 p for convenience once again. Taking the variation of the action with respect to the non-dynamical variables α, δλ, and β, I get the following three constraint equations, Substituting the solution of the above equations into the action, the quadratic action of the scalar perturbation takes the following form, where the eclipse denotes the potential term and gradient term. The ghost free condition for the scalar perturbation reads We have reproduced the ghost free condition eq. The scalar graviton is canonically normalised as The canonically normalised scalar action reads where the sound speed reads where ρ dark is energy density of the dark sector defined in the eq. (2.4), w dark is its equation of state, and ρ hde = cM 2 p /a 2 ϕ 2 is the energy density of the holographic dark energy. It is clear that the c 2 s becomes negative when the equation of state of the dark sector drops below −1, this is precisely the case for c < 6 in the asymptotic future [13] [22]. Therefore, if we demand the freeness of the gradient instability throughout the whole cosmic expansion history, the parameter region with c < 6 needs to be excluded. However, c < 6 is still allowed, if we only demand the freeness of the gradient instability at present and in the past, and allow the instability to occur in the future. The mass of the scalar mode reads where Ω rad is the fraction of the dark radiation defined by Ω rad ≡ ρ rad /3M 2 p H 2 , and Ω hde is the fraction of the holographic dark energy defined by Ω hde ≡ ρ hde /3M 2 p H 2 . The Compton wavelength of the scalar mode is about Hubble radius size, up to a coefficient whose value is dependent on the parameters b and c.
One may wonder whether the mass of scalar graviton is always positive along the cosmic evolution to ensure the stability in the scalar sector. This problem may require a thorough numerical analysis, and thus it will not be carried out in the current work. Instead, I would like to show the tachyon freeness condition in the asymptotic de-Sitter phase where c = 6 [13] and w dark → −1. The scalar mass in the asymptotic phase reads The tachyon freeness condition translates to The stability analysis for general space-time solution will be performed in the future work.

Tensor perturbation
The quadratic action for the tensor perturbation reads where The Compton wavelength of the tensor graviton is about the size of the future event horizon.
The following tachyon freeness condition is required to ensure the stability of the tensor sector. In the asymptotic de-Sitter phase where c = 6, w dark = −1, the stability condition translates to 0 < b < 24, which is weaker than the stability condition eq. (3.41) derived in the preceding subsection of the scalar mode analysis. The non-vanishing mass modifies the dispersion relation of the gravitational waves, and the group velocity is constrained with high precision −3 · 10 −15 < c gw − 1 < 7 · 10 −16 at low redshift regime by the multi-messenger observation GW170817 [31]. It gives us an upper bound on the graviton mass around m g < 10 −23 eV. Our graviton mass is much lower than this upper bound by around 10 orders of magnitude. It remains challenging to directly probe the non-vanishing graviton mass effects at the late time epoch. However, the size of the future event horizon is much shorter during the early universe, and thus it leads to a much larger graviton mass, which may give rise to some interesting observational effects on the stochastic gravitational waves background.

Vector perturbation
I do not expect to find the dynamical degrees of freedom in the vector sector at the linear perturbation level, as the residual symmetry eq. (3.21) has projected out all temporal derivative of the Goldstone vector bosons. I will try to confirm it in this subsection. The quadratic action of the vector perturbation reads Taking the variation with respect to the Lagrangian multiplier S i , I get Substituting the solution of this constraint equation back to the action, I get The kinetic term of the vector mode has vanished, and the action is subject to a new constraint F i = 0. After substituting this solution back to the action, the whole action vanishes and I conclude that no dynamical degree of freedom is found in the vector sector at the linear perturbation level. Whether there exists dynamical degrees of freedom at the higher order level requires a fully non-linear analysis, which will be provided in the next subsection.

The Hamiltonian analysis confirms the total number of degrees
I have performed the Hamiltonian analysis in the subsection (3.2), and I found that the original theory contains 6 modes. Compared to the Hamiltonian of the original theory (2.6), the modified one (3.26) strips off the term that is non-linear in the shift vector. Consequently the primary constraint Π i ≈ 0 commutes with all primary and secondary constraints. The consistency conditions for the momentum constraints conserve in time yield to another 3 tertiary constraints, Now let's add these 3 new tertiary constraints into the constraint set φ A , and remove the Π i ≈ 0 as it is first class, we have the new set We can show that the rank of the following 10 × 10 matrix is 10, no new constraints are generated and the algebra closes here. Therefore, all constraints in the eq. (3.49) are second class, and they remove 10 degrees in the phase space. On the other hand, the constraints Π i ≈ 0 remove 6 degrees in the phase space as they are first class. The total number of the leftover dynamical degrees of freedom in the phase space is which corresponds to 3 dynamical degrees of freedom in the physical space-time. As we have shown in our perturbation analysis, 2 of them are tensor modes, and the rest one is a scalar graviton. The additional degrees of freedom are fewer than the broken symmetry generators, because the Lorentz-invariance of the background space-time is broken and therefore the Goldstone theorem needs to be understood in the extended sense.

Physical interpretation of the UV-IR correspondence
The UV-IR correspondence is based on the assumption that the vacuum energy arising from the quartic divergence should not exceed the energy density of a black hole in the same size, and the effective quantum field theory should break down at the Schwarzschild radius scale. Take the note of that the Hubble radius of our universe coincides with its Schwarzschild radius, and thus the bound should be saturated on the cosmological background and we have the UV-IR relation Λ ∼ M p /L. One may ask how this UV-IR correspondence appears naturally in our effective field theory framework? The answer is that it should be interpreted as the scale where the scalar graviton strongly couples to itself. Let's look at the 4-votex interaction in the decoupling limit, where π 0 c is the canonical normalised Goldstone boson, and ω is the frequency of π 0 c . The strength of the coupling exceeds unity if ω > M p /ϕ, and our effective field theory breaks down. As shown in our perturbation analysis, we have learnt that ϕ is the graviton's Compton's wavelength (remember ϕ has the length dimension, and we set the scale factor a = 1 in the decoupling limit). Therefore, a sensible and physical interpretation of this UV-IR correspondence is that the UV cut-off scale of our EFT is inversely proportional to the Compton wavelength of our graviton.

Conclusion and Discussion
In the current work, a general covariant local field theory is proposed for the holographic dark energy model. I started from the mini superspace action proposed a few years ago [13], showed that the direct, and perhaps the simplest way of covariantizing the action gives rise to additional 4 degrees of freedom (6 in total including 2 degrees in the gravitational wave sector), which leads to the Ostrogradsky ghost instability. To remedy the ghost pathology, I introduced a new symmetry to prohibit the problematic terms in the theory. It has turned out that the low energy effective field theory of the holographic dark energy model is actually the Lorentz-violating massive gravity theory, whose graviton has 3 polarisations including 1 helicity 0 mode and 2 helicity 2 modes. The Compton wavelength of the graviton is about the size of the future event horizon of our universe, as shown in the linear perturbation analysis in the unitary gauge. To confirm the total number of dynamical degrees of freedom of our theory at the fully non-linear level, I have performed the Hamiltonian analysis and found that the momentum constraints yield to another 3 tertiary constraints. These 6 constraints, which are all second class, eliminate 3 out of 4 additional degrees and therefore there are only 3 dynamical degrees of freedom in total in the gravity sector, including a scalar mode and two tensor modes. Our effective field theory breaks down at the scale Λ ∼ M p /L, where L is the graviton's Compton wavelength, due to the strong coupling of the scalar graviton above this scale, which offers a natural and physical interpretation for the UV-IR correspondence in the holographic dark energy model.
Our effective field theory provides a general framework in which the holographic dark energy model can be tested. For instance, the non-vanishing mass of the gravitational waves, which is small at late time epoch but sizeable during the early universe, leads to a modified stochastic gravitational waves background. On the other hand, the existence of the scalar graviton whose Compton wavelength is about the size of the Hubble radius (up to an order O(1) ∼ O(10) coefficient) during inflation may be tested by the cosmological collider physics [32]. It is also very intriguing to ask what would happen if the Stueckelberg field, say the time like one ϕ, couples to the standard model fields. Our EFT framework also provides a general setup where the perturbation theory of the holographic dark energy can be developed, which allows us to investigate its impacts on the cosmological structure formation in detail. All these possibilities warrant further scrutiny.
I shall end by commenting that the global Lorentz invariance is broken in our scalar field configuration. As an analog to the Higgs mechanism in the particle physics, one may expect a similar symmetry restoration to occur at high energy scale, where the Lorentz invariance is recovered and the graviton becomes massless again. However, this Higgs-like mechanism is till missing. In fact, to my best knowledge it is still an open question for all massive gravity theories.