Cosmological screening and the phantom braneworld model

The scalar and vector cosmological perturbations at all length scales of our Universe are studied in the framework of the phantom braneworld model. The model is characterized by the parameter ΩM≡M3/2m2H0\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\Omega _M\equiv M^3/2m^2H_0$$\end{document}, with M and m the five- and four-dimensional Planck scales, respectively, and H0\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$H_0$$\end{document} the Hubble parameter today, while ΩM→0\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\Omega _M\rightarrow 0$$\end{document} recovers the ΛCDM\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\Lambda \mathrm CDM$$\end{document} model. Ignoring the backreaction due to the peculiar velocities and also the bulk cosmological constant, allows the explicit computation of the gravitational potentials, Φ\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\Phi $$\end{document} and Ψ\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\Psi $$\end{document}. They exhibit exponentially decreasing screening behaviour characterized by a screening length which is a function of the quasidensity parameter ΩM\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ \Omega _M$$\end{document}.


Introduction
In the braneworld (BW) model the 3 + 1-dimensional Universe we live in is a timelike hypersurface (the brane) of codimension one or more, embedded in a higher dimensional spacetime (the world), see [1,2] for a vast review and also references therein. Unlike the higher dimensional theories such as Gauss-Bonnet gravity, e.g. [3], in the BW model all standard model matter fields are confined on the brane whereas only gravity can propagate in the extra dimension(s).
The existence of the extra dimension implies departure from General Relativity. For example in the Randall-Sundrum model with a single extra dimension, the modification occurs at the small scales [4,5]. The extra dimension needs to neither be small nor compact and can even be infinite. Compact extra dimensions, on the other hand, imply an infinite and discrete Kaluza-Klein spectrum on the brane, see e.g. [6]. We further refer our reader to [7][8][9][10][11][12] for a description a e-mail: sbhatta@iitrpr.ac.in b e-mail: skousvos@physics.uoc.gr c e-mail: sromanop@physics.uoc.gr d e-mail: tomaras@physics.uoc.gr of fitting the galaxy rotation curves and the study of gravitational lensing in this model. While the extra dimension is usually taken to be spacelike, we refer our reader to [13] for a timelike extra dimension.
Discussions on static solutions such as a black hole in the BW model can be seen in [14][15][16][17] and references therein. For the so called two branch RS-I model, from the modification of Newton's law, the upper bound on the bulk anti-de Sitter radius turns out to be l 14 µm; whereas for the one branch RS-II model, the binary gravity wave data puts a bound : l 3.9 µm [18]. Probing the extra dimensional effects by studying the strong gravitational lensing can be seen in [19]. We refer our reader to [20] for a modification of the RS model with cosmological constants associated with both the bulk and the brane, fine tuned to make the bulk flat. This scenario is in particular helpful to estimate the energy lost by the brane via the Kaluza-Klein gravitons. In [21,22], the effect of brane -bulk energy exchange on cosmology was investigated and a model where our current universe is obtained as a late time attractor was proposed. We further refer our reader to [23] for a vast review and an exhaustive list of references pertaining to gravity and cosmology in the context of the braneworld model.
In this paper, we shall be interested in an extension of the Dvali-Gabadadze-Porrati braneworld (DGP) model [24][25][26][27] containing in the action, the 4-dimensional Ricci scalar on the brane, induced by the one loop correction due to the graviton-matter interaction, and the extrinsic curvature of the brane. This model, unlike the Randall-Sundrum case, modifies gravity only beyond a characteristic length scale, depending on the five-and four-dimensional Newton constants. The relevant equation of motion gives rise to two branches of cosmological solutions, both with flat spatial sections, one being self accelerated without requiring any dark energy/cosmological constant, whereas the other branch (the normal branch) requires at least one cosmological constant to accommodate for the current accelerated expansion [28][29][30]. However, the former was shown to have ghost instability in subsequent works [31,32], leaving only the "normal" branch to be a possible alternative to the CDM model. Furthermore, the equation of state parameter for the effective dark energy source is time dependent, w(t), and turns out to be less than minus one today [33][34][35][36][37]. For a certain range of parameter values, w(t) will reach asymptotically the value − 1 (the de Sitter phase). Otherwise, the universe can even leave at some stage the phase of accelerated expansion reentering matter domination, thus evading the so called phantom disaster [38]. Since w(t) < −1 in the current epoch, this model is often called "phantom braneworld model". Interestingly, this model indicates that the expansion of our universe was stopped at redshift z 6 and 'loitered' there for a long period of time favouring structure formation. Arguments supporting this, based on the observed data of population of the quasistellar objects and supermassive black holes in 6 z 20 can be found in [36]. Scalar cosmological perturbation theory in the phantom braneworld model and further details are studied in [39][40][41], while in [42] the stability analysis of large scale cosmic structures via their size-versus-mass study in the context of the present model and in the presence of a bulk cosmological constant, was performed. We also refer our reader to [43] for constraints on the braneworld model via gravity wave data. See also [44] for phenomenological arguments in the favour of w(t) < −1 via the neutrino mass higherarchy.
The braneworld model under study here is assumed to have 'zero thickness' in the extra dimension. Interesting effects however, may arise when one considers a thick brane [45,46]. In particular, in such a scenario, with a large extra dimension, one can have a new energy scale on the brane, determined by both brane thickness and the size of the extra dimension. For energies much larger than this new scale, the physics in the brane depends upon the position along the extra dimension, while for much smaller energies the equivalence principle may be violated, resulting in certain fine tuning to preserve it.
Given that the phantom braneworld model modifies gravity significantly at large scales, it becomes an interesting task to investigate this model's prediction at arbitrarily large distances. One such arena seems to be the study of screening effects, where certain terms in the scalar perturbation equation, which we can ignore at small scales, lead to modifications of the gravitational potential at large scales [47][48][49][50][51][52][53][54][55][56][57]. By approximating the inhomogeneities of our universe as delta function sources, a first order analytical formalism for the cosmological scalar and vector perturbations for the CDM model was developed recently in [47], where a Yukawa-like fall-off of the gravitational potential was derived at large scales. Various extensions of this work, including the case of interacting fluid sources, can be found in [48][49][50][51][52][53]. Discussions on the N -body simulations in the context of cosmic screening can be seen in [54,55]. We further refer our reader to [56,57] for second order computations on the scalar perturbation pertaining respectively to the CDM and the Einstein de Sitter models. The extra dimensional scenario is certainly not included in the above examples. Motivated by this, we shall study in this work the first order cosmological screening in the phantom braneworld model. Our chief goal would be, apart from casting the perturbation equations in a suitable form and solving them, to point out differences of this model from CDM, that can arise at very large scales.
The paper is organized as follows. In the next section we briefly review the phantom braneworld model. In Sect. 3 we develop the first order equations pertaining to the scalar and the vector perturbations with no bulk cosmological constant. In Sect. 4 we solve for the scalar perturbation ignoring the peculiar velocities, and compare it both analytically and numerically with the CDM model. We conclude with a Sect. 5.
We shall use mostly negative signature for the metric and will set c = 1 throughout.

The phantom braneworld model
Let us first briefly review the basic features of the phantom braneworld model, details of which can be seen in e.g. [40] and references therein. The relevant action is given by, where R and R are the Ricci scalars corresponding to five (the bulk) and four dimensions (the brane) and M and m are the respective Planck masses. The quantity 5D is the cosmological constant in the bulk and σ is the brane tension, related to the brane cosmological constant by = σ/m 2 . K is the trace of the extrinsic curvature of the brane. L(g μν , φ) stands collectively for all matter fields, φ, confined to the brane and g μν is the induced metric on it. For our current purpose, φ would correspond only to the cold dark matter.
Being interested in the 3 + 1-dimensional physics, we choose to measure energies in units of the four-dimensional Planck mass m. So, we set m = 1 throughout.
Using the Gauss-Codacci relations, the Einstein equations on the brane become are convenient parameters, and The tensor C μν is traceless, coming from the projection of the five-dimensional Weyl tensor onto the brane. Taking the divergence of Eq. (2) yields the constraint equation, The spatially homogeneous Einstein equation reads (in conformal time, η) with the cold dark matter as the source, where a ≡ a(η) is the scale factor, H = a −1 da/dη is the Hubble rate andρ the time independent background homogeneous cold dark matter density in co-moving coordinates.
The constant C is due to the existence of the Weyl tensor in the bulk. Due to the radiation like behavior of the term containing C, it is often named "Weyl radiation". We shall ignore its backreaction effects onto the cosmological background, though we shall take into account the inhomogeneous perturbations of the projection of the Weyl tensor. We will also ignore the backreaction effects of 5D . Taking M → 0 in the above equation one recovers the CDM limit. Notice that Eq. (6) in the absence of 5D and C may be conveniently expressed as We will also need the derivative of this equation with respect to conformal time η Consequently, Eq. (6) takes the simple form in everything that follows we have replacedρ in favor of 3H 2 m a, M 3 in favor of 2 M H/a and in favor of 1 We would like to mention that the quantity M is often defined as √ l in the related literature. 3H 2 σ . It is very convenient since everything we derive may be expressed as functions of m and M , only ( σ is solved for from Eq. 10). The advantage of this procedure is twofold, firstly these parameters are dimensionless and we claim rather intuitive to handle, secondly these will make comparison to CDM trivial by simply taking M to zero.

Derivation of scalar and vector perturbation equations
We shall extend below the linear perturbation scheme developed for the CDM model in [47] to the phantom braneworld model described in the preceeding section. We start with the ansatz for the first order McVittie metric on the brane in the Cartesian coordinates, where , and B i 's are respectively the scalar and vector perturbations and the bold font is used to indicate a vector, which determines the position in space where the potentials are evaluated at. Note that unlike the CDM, = here, owing to the anisotropic stresses originating from the bulk, e.g. [40]. We shall consider the backreaction effects due to N self gravitating moving point masses. Following [47], we define the proper interval for the n-th mass, The peculiar velocities appearing above can be evaluated by subtracting from the observed velocity of the mass, the velocity due to the Hubble flow, e.g. [38]. The energy momentum tensor for these point masses is then given by where x n is the value of the x coordinate (as defined in the metric Eq. 11) where the nth particle is located at. Existing data shows that the peculiar velocities are in general rather small or non-relativistic, at most of the order of 10 6 ms −1 [58]. Putting these all in together, we find from Eq. (13) the energy momentum tensor up to the first order, where each ρ n corresponds to a delta function point mass located at r n , We decompose the total energy density ρ in Eq. (14) as, where δρ(η, x) stands for the contribution of the inhomogeneities. The index n runs over all N particles in the Universe. Note here that δρ is not treated as a perturbation, due to the fact that it is dominant at small scales (see [60]). Since we must have | |, | | 1 in Eq. (11), we write from Eq. (14) at first order, the geodesic equation for the nth particle in Eq. (12) also reads, where the 'prime' denotes differentiation once with respect to the conformal time η and the variations δT μν , δC μν and δ Q μν depend on both space and time. Since we wish to build a perturbation scheme valid all the way to superhorizon scales, we cannot assume that the perturbations' spatial variations dominate over the temporal ones, unlike the case of the study of cosmic structures, e.g. [38]. Finally, we come to the perturbation of the Weyl tensor's projection onto the brane, δC μν . Its most generic form is given by, e.g. [40], where δπ i j = (∇ i ∇ j − g i j /3)δπ C ( stands for the Euclidean 3-Laplacian) is trace free and δρ C , v C and δπ C are scalars. In particular, v C can be regarded as a momentum potential, whose backreaction effects will also be ignored, while considering its time evolution also negligible. The Einstein equations on the brane (Eq. 2), at first order read, after using Eqs. (17,19), which is the 00 component, and for We also have for the vector perturbation, (22) where as earlier is the Laplacian on the Euclidean 3-space and also the function m eff ≡ m eff (η) has been introduced The CDM limit in the above equations is obtained by letting M → 0 in which case we recover the results of [47]. At small length scales relevant to cosmic structures, the spatial derivatives of the potential in Eq. (20) dominate over its temporal derivatives and the other effective mass-like terms appearing on the left hand side. Accordingly, at such small scales, Eq. (20) reduces to the Poisson equation, yielding a gravitational potential falling off as 1/r , along with a modified Newton's constant [42]. For CDM in particular, we have δρ C = 0, yielding Newton's potential. However, at length scales much larger than those of cosmic structures, the temporal derivative and the effective mass terms can be comparable and, as we will show in Sect. 4, this leads to a significant modification in the behavior of the solution of Eq. (20), as is expected due to the presence of the mass-like term on its left-hand side.
In this work, we are chiefly interested in distinguishing the phantom braneworld model from CDM with respect to the cosmological screening, which is certainly impossible unless we go to very large length scales. Note that at such scales, the backreaction effects due to the peculiar velocities, which are essentially non-relativistic, would be negligible, e.g. [58]. Thus for our current purpose, we shall from now on ignore the peculiar velocities (and hence the vector perturbation) throughout.

Solutions ignoring peculiar velocities
On the other hand we can write Eq. (21) and recall that in a marginally closed universe with a vanishing bulk cosmological constant, one has [39], Combining Eqs. (28) and (29) to eliminate δπ C we get In order to solve for and we want one more equation. This comes from the spatial component of Eq. (5), after using Eqs. (18) and (27), we obtain We can substitute Eq. (30) into Eqs. (27) and 31 to obtain a system of two equations with only two unknowns, the perturbations and . The constant in Eqs. (28) and (31) has to be zero in order for the potential to be vanishing at infinity. The solution of the system is straightforward and − 9H 2 where For a single particle -a single central over-density -the solution for the potential , valid for all length scales is where m 0 is the mass of the central overdensity and r = |x|.
We have dropped the 1 3 which is generated by the existence of an infinite number of point particles. We will prove that this occurs naturally when considering an isolated sub-region of the Universe (e.g. the observable Universe) at the end of this section. The exponential appearing above clearly indicates the suppression of Newton's potential at large scales, originating from the term present in the perturbation equation behaving as an effective mass. Thus the length scale, λ, should be interpreted as a screening length.
In every case we can find solving Eq. (32) = I  [59] and we examine the gravitational potentials for one particle with mass M = 1.989 · 10 30 kg. Figure 1 depicts the behavior of the effective mass density parameter and the screening length versus M . The CDM limit is obtained by letting M → 0.
We also note that since the screening length is typically of the order of O(10 3 ) Mpc (Fig. 1), at length scales comparable of the size of a typical cosmic structure i.e. O(100) Mpc, Eq. (37) recovers the 1/r fall-off of the gravitational potentials. However, the 1/m 2 eff, ≡ I /m 2 eff, and 1/m 2 eff, terms present modify Newton's 'constant' in and respectively and make it time dependent, as discussed in [42]. Figure 2 depicts the behavior of the effective Newton's constant for and . In the M → 0 limit both of them aproach 1 recovering the CDM limit. Note also that in this limit setting furtherρ → 0 ( m → 0) removes the exponential fall off since then λ → ∞ (cf., Eqs. 34, 36), yielding Newton's potential for a point mass located in a de Sitter universe. It is easy to verify that, as expected, this is the linearized approximation of the Schwarzschild-de Sitter metric in the McVittie coordinate frame. Similar conclusions hold  (r ) r/Mpc and (r ) r/Mpc versus M and log(r /Mpc) for one particle with mass equal to one solar mass for the potential one particle . Finally, we depict the potentials in Fig. 3.
We would now like to show that with respect to the universe visible to an observer located at some point x, we can actually get rid of the first term in Eq. (35). Indeed, let N be the total number of point sources in Eq. (35) and let N be the number located within the Hubble horizon radius of an observer located at x. Clearly, we may expect that only these N particles would contribute significantly into Eq. (35). On the other hand, since we should have N → ∞ in order to obtain a non-vanishingρ, it is natural to consider N N . We next split the summations in Eq. Since the second summation gets contributions from all particles outside the Hubble horizon of the observer, we can average the potential of this part following [47] and using ρ = ∞ 1 m n /V ≈ ∞ N +1 m n /V . It is easy to see that this average cancels-out the constant (1/3) term in Eq. (35), leading to many particle; average = − 1 8π m 2 eff, a N n=1 m n |x − x n | e − |x−xn | λ where "average" in the subscript refers to the aforementioned averaging over sources located outside the observer's Hubble horizon. Note that the above formula has a smooth one particle ( N = 1) limit, recovering Eq. (37).

Discussion
At very large length scales of our universe not decoupled from the cosmic expansion, one might expect the gravitational potential to be modified from that of Newton's, as has explicitly been demonstrated for the CDM model in [47]. It is an interesting task to investigate the same for other viable gravity models as well. Being motivated by this, we have investigated the cosmological screening at such large length scales for the phantom braneworld model described in Sect. 2, with the expectation that the qualitative differences of this model compared to CDM should be maximum at the (super-)horizon scales of our universe. We have presented the equations governing the first order scalar and vector perturbations in Sect. 3. Finally, by ignoring the backreaction effects due to the bulk cosmological constant and the vector perturbation, we have demonstrated analytically and numerically, the behaviour of the two potentials up to the superhorizon length scale in Sect. 4. It seems to be an interesting task to investigate the tensor perturbation for this model in an early universe scenario. We hope to address this issue in future work.