#### Abstract

We investigate the ghost model of dark energy in the framework of DGP braneworld. We explore the cosmological consequences of this model by determining the equation of state parameter, , the deceleration, and the density parameters. We also examine the stability of this model by studying the squared of the sound speed in the presence/absence of interaction term between dark energy and dark matter. We find out that in the absence of interaction between two dark sectors of the universe we have in the late time, while in the presence of interaction can cross the phantom line . In both cases the squared of sound speed does not show any signal of stability. We also determine the statefinder diagnosis of this model as well as the plane and compare the results with the CDM model. We find that plane meets the freezing region in the absence of interaction between two dark sectors, while it meets both the thawing and the freezing regions in the interacting case.

#### 1. Introduction

The current acceleration of the universe expansion is strongly confirmed by the type Ia supernova observations [1] and also supported by the astrophysical data from WMAP [2]. A component which is responsible for this acceleration is called dark energy (DE) with negative pressure, which can overcome the gravity force between the galaxies and push them to accelerate, though its nature and origin is still an open question in the modern cosmology. There are two approaches for explanation of the cosmic acceleration. The first one is the modified gravity models such as gravity [3] and scalar-tensor theories [4], and the second approach is the idea of the existence of a strange type of energy whose gravity is repulsive such as the cosmological constant [5] and the dynamical DE models [6, 7]. Against the cosmological constant which has constant equation of state (EoS) parameter, , further observations detect a small variation in the EoS parameter of DE in favor of a dynamical DE with in the past and even in the late time [8].

An interesting model for probing the dynamical DE model is the ghost dark energy (GDE) model proposed in [9]. The advantages of this model are that it does not introduce any new degree of freedom in contrast to the most DE models which explain the accelerated expansion by introducing new degree(s) of freedom or by modifying the underlying theory of gravity. This is important because, with introducing new degrees of freedom, one needs to investigate the nature and new consequences in the universe so it seems to be impressive and economic if we can explain DE puzzle by using currently known fluids and fields of nature. Actually, GDE model which is based on the Veneziano ghost in Quantum Chromodynamics (QCD) can act as the source of DE [10] and its existence are required for resolution of the U problem in QCD [11]. Indeed, the ghosts are decoupled from the physical states and make no contribution in flat Minkowski space, but it produces a small vacuum energy density in a dynamic background or a curved spacetime proportional to , where is the Hubble parameter and is QCD mass scale of order a [12]. Different features of GDE have been studied in ample details [13].

Independent of the DE puzzle, for explanation of the cosmic acceleration, special attention is also paid to extra dimensional theories, in which our universe is realized as a -brane embedded in a higher dimensional spacetime. Based on the braneworld model, all the particle fields in the standard model are confined to a four-dimensional brane, while gravity is free to propagate in all dimension. One of the original models of braneworld is introduced by Dvali-Gabadadze-Porrati (DGP) [14], which describes our universe as a 4D brane embedded in a 5D Minkowskian bulk with infinite size. In this model the recovery of the usual gravitational laws on the brane is obtained by adding an Einstein-Hilbert term to the action of the brane computed with the brane intrinsic curvature. It is a well-known that the DGP model has two branches of solutions. The self-accelerating branch of DGP model can explain the late time cosmic speed-up without recourse to DE or other components of energy [15, 16]. However, the self-accelerating DGP branch has ghost instabilities and it cannot realize phantom divide crossing by itself. To realize phantom divide crossing it is necessary to add at least a component of energy on the brane. On the other hand, the normal DGP branch cannot explain acceleration but it has the potential to realize a phantom-like phase by dynamical screening on the brane. Adding a DE component to the normal branch solution brings new facilities to explain late time acceleration and also better matching with observations. These are the motivations to add DE to this braneworld setup [17, 18]. In this work we would like to investigate the GDE model in the framework of the DGP braneworld. This study is of great importance, since we can incorporate and disclose the effects of the extra dimension on the evolution of the cosmological parameters on the brane when the DE source is in the form of GDE.

This paper is organized as follows. In Section 2 we formulate the GDE model in the context of the DGP braneworld. We also consider both interacting and noninteracting cases and explore various cosmological parameters as well as cosmological planes. Besides the discussion of instability analysis, we study the plane and properties of statefinder parameters. We finish with closing remarks in Section 3.

#### 2. The GDE in the DGP Model

In the DGP cosmology, a homogeneous, spatially flat, and isotropic -dimensional brane which is embedded in a -dimensional Minkowskian bulk can be described by the following Friedmann equation [19]or equivalentlywhere is the Hubble parameter, [20] is the crossover length scale reflecting the competition between 4D and 5D effects of gravity, and corresponds to the two branches of solutions of the DGP model. Before going any further, it is worth noting that if (early times) the 4D general relativity is recovered; otherwise the 5D effect becomes significant. Also corresponds to the self-accelerating solution where the universe may accelerate in the late time purely due to modification of gravity [15, 16], while can produce the acceleration only if a DE component is included on the brane. Here, to accommodate GDE into the formalism we take .

The fractional energy density parameters are defined aswhere is the Hubble parameter at redshift . The Friedmann equation (2) can be rewritten in terms (3) asWe introduce , which comes from the extra dimension. Thus the Friedmann equation (4) can be reexpressed asFor the GDE density we have where is a constant of order and is the QCD mass scale [12]. Taking the time derivative of the energy density and using (6) we obtain For the FRW universe filled with DE and DM, with mutual interaction, the energy-momentum conservation law can be written aswhere is considered as the interaction term between DE and DM also is the coupling constant of interaction .

We know that (i) our universe is in a DE dominated phase and (ii) our universe that is our habitat is stable. These imply that any variable DE model should result in a stable DE dominated universe. So it is worth investigating the stability of the GDE in DGP braneworld against perturbation. The intended indicator for checking the stability of a proposed DE model is to study the behavior of the squared sound speed () [21]. If we have the classical instability of a given perturbation because the perturbation of the background energy density is an oscillatory function and may grow or decay with time. When , we expect a stable universe against perturbations because the perturbation in the energy density propagates in the environment. We continue discussion of stability in the linear perturbation regime where the perturbed energy density of the background can be written as where is unperturbed background energy density. For the energy conservation equation () which yields [21] we encounter two cases. In the first case where , we observe an ordinary wave equation which has a wave solution in the form (stable universe). In the second case where , the frequency of the oscillations becomes pure imaginary and the density perturbations will grow with time as (unstable universe). Since plays a crucial role in determining the stability of DE model, we rewrite it in terms of EoS parameter as where is the pressure of DE, is the total energy density of DE and DM, and is the energy density ration.

On the other sides, Sahni et al. [22] proposed new geometrical diagnostic pair parameter , known as statefinder parameter, for checking the viability of newly introduced DE models. Unlike the physical variables which depend on the properties of physical fields describing DE models, the statefinder pair primarily depends on the scale factor and hence it depends on the metric of the spacetime. The and parameters are defined as [22] where can rewrite as and then Let us note that in the plane, corresponds to a quintessence-like model of DE and corresponds to a phantom-like model of DE. Also the studies on a flat CDM model and matter dominated universe (SCDM) show that for these models and , respectively. In above equations is the deceleration parameter which is given by In what follows we discuss the plane which introduced by Caldwell and Linder [23] for analyzing the dynamical property of various DE models and distinguish these models ( represents the evolution of ). The models can be categorized into two different classes: (i) and which present the thawing region. (ii) and which present the freezing region. It should be noted that the CDM model corresponds to a fixed point in the plane. We shall consider the noninteracting and interacting cases, separatively.** It deserves to mention here that in all figures we selected values for the fixed parameters, which is custom and the observational evidence confirms it, as****[24] or****that can be between 0.7 and .75 and****that can be between 60 and 100. It should be noted that the results and conclusions are not sensitive to these choices and also do not do any observational work on this title; thus we by helping papers and experience selected these values for the fixed parameters.**

##### 2.1. Noninteracting Case

We start to obtain the cosmological parameters for GDE in the DGP braneworld by ignoring the interaction term . The deceleration parameter can be obtained by taking the time derivative of (2), which lead to Using relation (16), we find Inserting (7) in (9) we have where by replacing (17) in it, we get Also we obtain from the above equation as Note that, in order to find the evolution of density parameter , we take the time derivative of relation , after combining the result with (7) and (17), which yieldsIn Figure 1, we plot the evolution of versus redshift parameter . It is obvious that tends to in the early universe where , while at the late time where , we have . Clearly, (20) for the EoS parameter shows that, at the late time where , the EoS parameter mimics the cosmological constant, namely, ,** which is showed in **Figure 2. In Figure 3, the behavior of the deceleration parameter is plotted and indicates that indeed there is a decelerated expansion at the early stage of the universe followed by an accelerated expansion. The energy density ratio is defined as , which by using (5) can be writtenDifferentiating (23) and (20) and then substituting the results in (12) we get the squared of sound speed as

The evolution of against for the noninteracting GDE in the framework of DGP braneworld is plotted in Figure 4. From graphical analysis of one concludes that this model does not indicate any signal of stability, that is, during the history of the universe. We can also find the statefinder parameters and by taking derivative of (18) and using (13) and (15). The results are

The graphical behavior of the statefinder parameters given in (25) and (26) is plotted in Figures 5 and 6, showing that, at late time where , we have which implies that GDE mimics the cosmological constant at the late time, as expected.

Let us study the trajectory in the statefinder plane and analyze this model from the statefinder viewpoint. For this purpose, we plot the statefinder diagram in the in Figure 7, which shows the cure gets to the point in the end, which implies that the model corresponds to the CDM model at the late time. For complementarity of the diagnostic, we also plot the trajectories of statefinder pair in Figure 8 which ends in the future to , corresponding to the de Sitter expansion.

The plane for the noninteracting GDE in the DGP scenario is shown in Figure 9. Again, we see that this plane corresponds to CDM model, i.e., , and the trajectory meets the freezing region as well.

##### 2.2. Interacting Case

Differentiating the modified Friedmann equation (2) and using (7) and (8) we reach Next, the EoS parameter can be determined by substituting (7) into the semiconservation law equation (9) and using (27). We findTaking differentiation with respect to from above equation we get where the prime indicates derivative with respect to . We can obtain the equation of motion for as

To illustrate the cosmological consequences of the interacting GDE in the DGP braneworld, we plot their evolution in terms of redshift parameter . In Figure 10, we present the graphical of versus for the different values of the coupling constant . As expected, we see both and for late time and early time, respectively. The graphical behavior of the EoS parameter for the different values of shows crossing of phantom line as plotted in Figure 11. The stability of interacting GDE in DGP model can be obtained by differentiating with respect to time of (23) and (29)

The evolution of the deceleration parameter and the squared of sound speed versus redshift parameter are plotted in Figures 12 and 13, respectively. In Figure 12, we see for different values of with the interacting GDE in DGP model, our universe has a phase transition from deceleration to an acceleration, while by keeping the same situation in Figure 13, this universe cannot be stable. As the value of decreases the severity of instability also decreases. Like previous section, the statefinder parameters are obtained by taking derivative of (28) and using (13) and (15)

We obtain for CDM model from (33) and (34) in the limiting case where , , and (in the late time). Also, Figures 14 and 15 show that and are positive through the entire life of the universe and turn to 1 and 0 at the late time, respectively. The evolutionary trajectories for the interacting GDE in the framework of the DGP braneworld for different values of are shown in Figure 16. From Figure 16, we can see that at the late time all curves tend to the CDM fixed point ; also different results in different evolution trajectories of statefinder which states that is smaller when is larger. The diagrams are plotted for different values of in Figure 17 which mimics the de Sitter expansion, namely, , in the far future where . In Figure 18, we plot the plane for different values of which show that the trajectories meet both the thawing and the freezing regions as well.

#### 3. Closing Remarks

We have made a versatile study on both noninteracting and interacting GDE in the framework DGP model through well-known cosmological parameters as well as planes. We summarize our results as follows. For noninteracting case, we have found that the density parameter tends to zero at the early universe while at the late time we have . Meanwhile the EoS parameter cannot cross the phantom line and mimics the cosmological constant at the late time (Figure 2). We have shown that our universe has a phase transition from deceleration to an acceleration, though we do not receive any signal of stability. The statefinder plane shows that the trajectory corresponds to quintessence model ( and ) while at late time we have for CDM model as expected. The plane in Figure 9 meets the freezing region as well.

For interacting case, we find that the density and the deceleration parameters as well as the EoS parameter are consistent with observational data. We have seem that as the value of decreases the severity of instability decreases. From plane, we can see that at the late time all cures tend to the CDM fixed point . Besides, for different values of , the different evolution trajectories of statefinder are shown which indicates that is smaller when is larger. The plane is plotted in Figure 17 which mimics the de Sitter expansion, namely, , in the far future where . In the end, the plane exhibits both freezing and thawing regions of the universe for all values of . Again, in this case which implies that interacting GDE in the DGP braneworld is not stable against perturbation.

#### Data Availability

No data were used to support this study.

#### Conflicts of Interest

The authors declare that they have no conflicts of interest.

#### Acknowledgments

The authors thank Shiraz University Research Council. This work has been supported financially by Research Institute for Astronomy & Astrophysics of Maragha (RIAAM), Iran.