# Considering Late-Time Acceleration in Some Cosmological Models

Authors: S. Davood Sadatian
Journal: Advances in High Energy Physics
 (2013)

---

## Abstract

We study two cosmological models: a nonminimally coupled scalar field on brane world model and a minimally coupled scalar field on Lorentz invariance violation model. We compare some cosmological results in these scenarios. Also, we consider some types of Rip singularity solution in both models.


---

## Body

## 1. Introduction

In theories of extra spatial dimensions, ordinary matter is captured on the brane but gravitation extends through the entire space time [1–3]. In these scenarios, the cosmological evolution on the brane is taken by an effective Friedmann equation that combines the effects of the bulk in a nontrivial kind [4–6]. In other view point, the important result in brane models is an alternative scenario for late-time expansion of the universe. This result predicts deviations from the 4-dimensional gravity at limited distances. In other hand, the model considered by Dvali, Gabadadze, and Porrati (DGP) is different from previous models since it also predicts deviations from the standard 4-dimensional gravity in large distances [4, 7]. Generally one can study the effect of a caused gravity term as a quantum discipline in any brane world model. The existence of a higher dimensional embedding space lets for the bulk or brane matter that can evidently affect the cosmological evolution on the brane [8]. A special form of bulk or brane context is a scalar field. Scalar fields take an important key both in models of the early universe and late-time acceleration. These scalar fields give a dynamical model for matter fields in a brane world scenarios [9–12].In other hand, Lorentz invariance violation models (LIV) have been considered in the scalar-vector-tensor theories [13]. It has described that Lorentz violating vector fields influence the dynamics equations in the inflationary models. An interesting result of this model is that the exact Lorentz violating inflationary solutions are dependent on the absence of the inflation potential. In this case, the inflation is exactly collaborated with the Lorentz violation [14]. Therefore, we study this symmetry breaking on the dynamics of equation of state for some cosmological aspects.However, some evidence from supernova data [15, 16], (CMB) results [17–19] and (WMAP) data [20–24] points out an accelerating phase of cosmological expansion, and this characteristic shows that the picture of universe by pressureless fluid is not enough; the universe should contain some type of additional negative pressure known as dark energy (for brief introduction in this field see [14]). Also, the merged analysis of the WMAP data with the supernova Legacy survey (SNLS) [22] compels the equation of state wde, in accordance with 74% donation of dark energy in the currently accelerating universe. Moreover, observations show some sort of a dark energy equation of state, wde<-1 [23]. Hence, a practical cosmological model should accept a dynamical equation of state that may have crossed the value wde=-1, in late time of cosmological evolution.In the following, we consider cosmological results of a nonminimally coupled scalar field on the brane and a minimally coupled scaler field in LIV model. Also we determine late-time behavior of our equations and obtain some restriction on the parameters of models to have an accelerating universe.We understand that dark energy is an important problem in modern cosmology, especially, if the equation of state parameterω is less than -1. In this paper, we study basically solving Friedman equations and obtain the evolution of the effective equation of state parameter ω. The phantom crossing conditions that we obtained for these set-ups, are interesting at some degree themselves. Furthermore, we know that just considering background quantity is not enough in such kind of unconventional cosmological models. The current constraint on ω is obtained by combining the result of the observation of CMB, which means that the information of the evolution of perturbations should be also included. Also we understand in the brane world model that the perturbations of brane are coupled with these of bulk, which gives nontrivial effect. Therefore, in first stage we just study the background dynamics, and the evolution of perturbations will be studied in the future.In other view, there is a no-go theorem which explicitly points out that a conventional dark energy model involving single degree of freedom in the frame of standard Einstein gravity is forbidden to realize such a scenario. This no-go theorem was proven in the appendix of [25]. In this regards, we point out the original work in [26] and several papers addressing this scenario [27–29]. This theorem is exactly the reason why we study a number of nonconventional dark energy models in realizing the Quintom scenario [24], such as the nonminimally coupled on brane world model and the Lorentz-violating model considered in the present paper. However, a nonminimally coupling scalar field identified on the brane model and scalar field coupling minimally to gravity in LIV model have some similar cosmological results; this is an interesting theoretical motivation of these models. Also we study other solutions admitted as Rip singularity that occur in the condition ω<-1 increases rapidly. However, it possible different types of singularity, dependingenergy density and scale factor how increases with time [30, 31].


## 2. Nonminimally Coupled Scalar Field on the Brane

Here we study a brane world model where a scalar field is coupling nonminimally to the Ricci scalar of the brane. In the following we only consider a scalar field in the matter Lagrangian without taking into account baryons, cold dark matter, and radiation. This kind of analysis suffers from a potential that may cause that the model is not able to explain the evolution of a realistic universe at background level when confronting with observations. However, this problem can be solved by suitable fine-tuning parameters of model.The action in the absence of ordinary matter can be given as [12, 32]
(1)S=∫‍d4x-g[1k42α(ϕ)R[g]-12gμν∇μϕ∇νϕ-V(ϕ)],
where we have made a general nonminimal coupling α(ϕ). For simplicity, in the following, we set k42≡8πGN=1. We can obtain Einstein equations with variation of the action with respect to brane metric
(2)Rμν-12gμνR=α-1𝒯μν,
where 𝒯μν, energy-momentum tensor of the scalar field, nonminimally coupled to gravity is taken by
(3)𝒯μν=∇μϕ∇νϕ-12gμν(∇ϕ)2-gμνV(ϕ)+gμν□α(ϕ)-∇μ∇να(ϕ),
where □ shows 4-dimensional d’Alembertian. We study the FRW universe with line element as the following:
(4)ds2=-dt2+a2(t)dΣk2,
where dΣk2 is the line element for a constant curvature k=+1,0,-1. The Ricci scalar is obtained with the equation of motion for scalar field as
(5)∇μ∇μϕ=V′-α′R[g],
where a prime denotes the derivative of each parameter with respect to ϕ, that may be written by
(6)ϕ¨+3a˙aϕ˙+dVdϕ=α′R[g],
where a dot means the derivative of each parameter with respect to t. The intrinsic Ricci scalar for a FRW brane is given as
(7)R[g]=6(H˙+2H2+ka2),
and Friedmann’s equations are determined by
(8)a˙2a2=-ka2+ρ3,(9)a¨a=-16(ρ+3p).
We take a scalar field, ϕ, only dependent on time. So, with (3), we obtain
(10)ρ=α-1(12ϕ˙2+V(ϕ)-6α′Hϕ˙),p=α-1(12ϕ˙2-V(ϕ)+2(α′ϕ¨+2Hα′ϕ˙+α′′ϕ˙2)),
where H=a˙/a is Hubble parameter. Now equation of state has the following form:
(11)w≡pρ=ϕ˙2-2V(ϕ)+4(α′ϕ¨+2Hα′ϕ˙+α′′ϕ˙2)ϕ˙2+2V(ϕ)-12α′Hϕ˙.
When ϕ˙=0, we have p=-ρ. In this illustration ρ depended on a and V(ϕ), and it has the duty of a cosmological constant. In the minimal case in which ϕ˙2<V(ϕ), by (9), we take p<-ρ/3 that shows a late-time accelerating universe. In nonminimal case, the position depends on the choice of nonminimal coupling parameter. In the following we show how for a suitable range of coupling parameter, a late-time accelerating expansion can be described. We first consider a moving domain wall of brane world to discuss quintessence behavior, and then we study a special nonminimal coupling for late-time acceleration.


### 2.1. Late-Time Acceleration in a Brane World Model

We study a bulk configured by two 5-dimensional anti de Sitter-Schwarzschild (AdSS5) black hole spaces combined on a moving domain wall. To insert this moving domain wall into 5-dimensional bulk, it is required to indicate normal and tangent to the domain wall by determination of normal instructing to the brane. We take that domain wall is identified at coordinate r=a(τ) where a(τ) is considered by Israel junction conditions [33]. Here we study the following line element [32]:
(12)ds5±2=-(k-η±r2+r2ℓ2)dt2+1k-η±/r2+r2/ℓ2dr2+r2γijdxidxj,
where ± is for left (−) and right (+) side of the moving domain wall, also ℓ is curvature radius of AdS5 manifold, and γij is the horizon metric. η±≠0 creates the electric part of the Weyl tensor on two sides [32].We assume that usual mater on the brane has a perfect fluid form,Tμν=(ρ+p)uμuν+phμν where ρ=ρm+σ and p=pm-σ. Energy density of the confined matter on the brane and analogous pressure is shown with ρm and pm, while σ is the brane tension. with Israel junction conditions [33] and Gauss-Codazzi equations, we obtain generalized Friedmann equations [32] as the following:
(13)a˙2a2+ka2=ρm3+ηa4+ℓ236ρm2,(14)a¨a=-ρm6(1+3w)-ηa4-ℓ236ρm2(2+3w),
where we require a Z2-symmetry by η+=η-≡η, and we have supposed pm=wρm. Following we take case: η=0. For η=0, each submanifold of bulk space time is accurate AdS5 space time. Now we study a confined nonminimally coupled scalar field on the brane and consider cosmological aspects. We use energy density and pressure of scalar field (10), and by (14), hence for cosmic expansion we have(15)a¨a=-16α(12ϕ˙2+V(ϕ)-6α′Hϕ˙)×(1+3ϕ˙2-2V(ϕ)+4(α′ϕ¨+2Hα′ϕ˙+α′′ϕ˙2)ϕ˙2+2V(ϕ)-12α′Hϕ˙)-ℓ236α2(12ϕ˙2+V(ϕ)-6α′Hϕ˙)2×(2+3ϕ˙2-2V(ϕ)+4(α′ϕ¨+2Hα′ϕ˙+α′′ϕ˙2)ϕ˙2+2V(ϕ)-12α′Hϕ˙),
where H=a˙/a is Hubble parameter. This is a complex equation; to describe cosmological aspects, we must study some limiting cases or specify α(ϕ), V(ϕ), and ϕ.Equation (15) shows the condition for an accelerating universe depending on the choice of nonminimal coupling and the scalar field potential. Hence we determine the following nonminimal coupling:
(16)α(ϕ)=12(1-ξϕ2).
For an accelerating universe we should have ρm+3pm<0 [12, 34]. If we take V>0, we obtain
(17)(1-3ξ)ϕ˙2-V(ϕ)-3ξϕ(ϕ¨+Hϕ˙)<0.
Dynamics of this scalar field are given as
(18)ϕ¨+3Hϕ˙+ξRϕ+dVdϕ=0.
We obtain
(19)(1-3ξ)ϕ˙2-V(ϕ)+3ξ2Rϕ2+6ξHϕϕ˙+3ξϕdVdϕ<0.
With ρm in this case, we determine
(20)y≡(1-ξϕ2)ρm-2V(ϕ)+ϕ˙2(12-3ξ)+3ξ2Rϕ2+3ξϕdVdϕ<0,
that is, the requirement to have an accelerated universe. In the following, we study the weak energy condition ρm≥0 and limit our consideration to the case by ξ≤1/6. Therefore we obtain(21)-2V+3ξϕdVdϕ<y<0,
and a main condition for cosmic acceleration is
(22)V-3ξ2ϕdVdϕ>0,ξ≤16.
For the model to have a cosmic acceleration by ξ>0, the potential V(ϕ) must change with ϕ  slower than power-law potential Vc(ϕ)=V0(ϕ/ϕ0)2/3ξ (Figure 1). In other hand, when ξ<0, the main condition for cosmic acceleration needs V to grow faster than Vc as ϕ increases [34]. Now we take a simple special example to show how this model processes. By potential of the form V(ϕ)=λϕn,  condition (22) takes λ(1-3nξ/2)>0. In this case just for ξ≤2/3n we have the accelerated expansion.Figure 1
Vvary for different values of the ϕ and nonminimal coupling parameter ξ for V=ϕ2/3ξ.More general case isη≠0 in (13) and (14). In the following we convert the generalized Friedmann equation as
(23)χ˙+4a˙a[χ+112α(12(1-2ξ)ϕ˙2-V(ϕ)-2ξϕ(ϕ¨+3a˙aϕ˙))]=0,
where χ=(a˙/a)2 is dark energy field. Then we take ϕ(t)≈A/tβ  where A is a constant and β an unspecified power will be studied. For scale factor we take the ansatz a(t)≈Btν where B is a constant and ν is to be determined. For scalar field potential we take V(ϕ)=λ(ϕ2/2). Therefore we obtain
(24)a¨a+(a˙a)2≈σ212αu+σ24αuρ(1-3w)+16α(V(ϕ)+Λ)-ϕ˙212α(1+4α′′)-α′3α(ϕ¨+3a˙aϕ˙),ϕ¨+3a˙aϕ˙+u′2uϕ˙2≈-6αuV′+10α′u(V(ϕ)+Λ)+4α′(u+α′u′)σu3ρ(1-3w)+8α′(u+α′u′)u3σ2,
where u=6[1-ξ(1-(16/3)ξ)(ϕ2/2)]. We take that usual matter on the brane has an equation of state such as  ρ=Dt-3ν(1+w). Hence we have
(25)ν(2ν-1)t2≈124{σ23+4Λ+(σ(1-3w)D6)1t3ν(1+w)+A2t2β[2λ+2ξ(Λ+4(3/8-ξ)9σ2)+4ξD(1-3w)(3/8-ξ)σ9×1t3ν(1+w)2βt2×(β(1-8ξ)+4ξ(3ν-1))(Λ+4(3/8-ξ)9σ2)](σ(1-3w)D6)1t3ν(1+w)},
and scalar field equation is given as
(26)β(β+1-3ν)Atβ+2≈-λAtβ[1-83ξ2A2t2β]-5ξ3Atβ[Λ+(λ2+ξ2(1-163ξ)Λ)A2t2β]-ξσ9(2σ+Dt2(1-3w))Atβ×[1+ξ(1-163ξ)(1+ξ)A2t2β],
where β>1. Two conditions must be solved numerically to limit the values of ν and β for an accelerating expansion. After some calculations we obtain
(27)-(1-3w)ξσ9D=Aβ(β+1-3ν),λ=-ξσ212,
where we passed over terms of order 𝒪(t-3β) and higher. Now, if we arrange λ, since ν(2ν-1)=(σ(1-3w)/144)D,  we find
(28)β=3ν-12±[(3ν-1)22-16ξν(2ν-1)]1/2,
where ν=2/3(1+w). For ν<1 the parameter β is real for any value; for ν>1 the exponent β is real only if we have
(29)ξ≤116(3ν-1)22ν(2ν-1).
These equations limited the values of nonminimal coupling ξ to obtain an accelerating expansion (Figure 2). Hence, it seems that nonminimally coupled scalar field identified on the brane world supplies natural candidate for late-time expansion.Variation ofξ with respect to ν with (29) (a) and variation of β with respect to ν and ξ with (28) by sign +, for sign − in (28) β has negative value. We choose special range of parameter ξ with (29) in (b) figure.
(a)
(b)Also we can use (11) to determine dynamical equation of state ω(t).Figure3 shows possible crossing of phantom divided barrier in this framework; this result is important because previous considerations have shown that crossing phantom divided line with a scalar field nonminimally coupled to gravity can exist [35, 36].Figure 3
Dynamical equation of stateω(t) varies for different values of the t and nonminimal coupling parameter ξ by (11) and specific choosing for parameters space.In the following, we explain the physical reasons why the crossing of the phantom divide can occur in this model. Here, we have a scalar field that is nonminimally coupled in gravity and an ordinary matter. Generally, according to [35], these sources of energy momentum can switch together and describe evolution of universe. In this regard, a promising candidate for the dark matter is either a positive cosmological constant or a slowly evolving scalar known as “quintessence.” Also the quintessence can be considered as an alternative method to resolve the cosmological constant problem. This is possible because instead of the fine tuning, the quintessence provides a model of slowly decaying cosmological constant. However, there exists another question for the present-day quintessence. This is to ask whether the expansion will keep accelerating forever or it will decelerate again after some time. This is very similar to the exit problem in the inflationary cosmology. According to the theory of quintessence, the dark energy of the universe is dominated by the potential of a scalar field ϕ which is still rolling to its minimum at V=0.


### 2.2. Rip Singularity

The additional interest to the models with the crossing phantom divideω<-1 is caused by their prediction of a Rip singularity [37]. Theoretically, the scale factor of the universe becomes infinite at a finite time in the future which was dubbed Big Rip singularity. There were proposed several scenarios to cure the Big Rip singularity: (I) to consider phantom acceleration as transient phenomenon. This is possible for a number of scalar potentials; (II) to account for quantum effects which may delay/stop the singularity occurrence [38]; (III) to modify the gravitation itself in such a way that it appears to be observationally-friendly from one side but it cures singularity; (IV) To couple dark energy with dark matter in the special way or to use special (artificial) form for dark energy equation of state [39]. Note that for quintessence dark energy, other (milder) finite-time singularities may occur. For instance, type II (sudden) singularity or type III singularity occurs with finite scale factor but infinite energy and/or pressure. The closer examination shows that the condition ω<-1 is not sufficient for a singularity occurrence [40]. First of all, a transient phantom cosmology is quite possible. Moreover, one can construct such models that ω asymptotically tends to -1 and the energy density increases with time or remains constant, but there is no finite-time future singularity [40]. Of course, most evident case is when Hubble rate tends to be constant (cosmological constant or asymptotically de Sitter space), which may also correspond to the pseudo rip [31]. Very interesting situation is related with Little Rip cosmology [30] where Hubble rate tends to infinity in the infinite future. The key point is that if ω approaches -1 sufficiently rapidly, then it is possible to have a model in which the time required for singularity is infinite. However, this position does not appear in our model, because H is finite.In the following, if we want to consider some types of Rip singularity, we should use Hubble parameter, energy density, and pressure relations that are determined in the previous equations. Therefore, after calculations withξ=0.14 (Figure 4), we obtain results as the following.Variation of pressure (a) and energy density (b) with respect to time for a nonminimal coupling model (whereξ=0.14).
(a)
(b)Case 1 (type II singularity (sudden)).
Conditions: scale factor and energy density reach finite value butp→∞.
In this model, we find that scale factor, pressure, and energy density are finite. Therefore, the model does not predict a sudden singularity.Case 2 (type III singularity).
Conditions, scale factor reaches finite value but (ρ→∞) and |p|→∞.
This means that the energy density grows so rapidly with time and scale factor does not reach the infinite value. The key difference between Cases1 and 2 for this model is that the potential of scalar field has a pole in a case of singularity. Here, according to numerical calculations, if scale factor reaches finite value, energy density and pressure are finite too. Hence, type III singularity does not occur in this setup.Case 3 (pseudo-rip singularity).
Conditions: pressurep→-ρ in infinite time and leads to H˙→0.
This means that the expansion of the universe asymptotically approaches the exponential regime. Obviously in our model with the previous setup, pseudo-rip singularity predict, because in infinite time, whenH˙ reach to zero, p→-ρ.


## 2.1. Late-Time Acceleration in a Brane World Model

We study a bulk configured by two 5-dimensional anti de Sitter-Schwarzschild (AdSS5) black hole spaces combined on a moving domain wall. To insert this moving domain wall into 5-dimensional bulk, it is required to indicate normal and tangent to the domain wall by determination of normal instructing to the brane. We take that domain wall is identified at coordinate r=a(τ) where a(τ) is considered by Israel junction conditions [33]. Here we study the following line element [32]:
(12)ds5±2=-(k-η±r2+r2ℓ2)dt2+1k-η±/r2+r2/ℓ2dr2+r2γijdxidxj,
where ± is for left (−) and right (+) side of the moving domain wall, also ℓ is curvature radius of AdS5 manifold, and γij is the horizon metric. η±≠0 creates the electric part of the Weyl tensor on two sides [32].We assume that usual mater on the brane has a perfect fluid form,Tμν=(ρ+p)uμuν+phμν where ρ=ρm+σ and p=pm-σ. Energy density of the confined matter on the brane and analogous pressure is shown with ρm and pm, while σ is the brane tension. with Israel junction conditions [33] and Gauss-Codazzi equations, we obtain generalized Friedmann equations [32] as the following:
(13)a˙2a2+ka2=ρm3+ηa4+ℓ236ρm2,(14)a¨a=-ρm6(1+3w)-ηa4-ℓ236ρm2(2+3w),
where we require a Z2-symmetry by η+=η-≡η, and we have supposed pm=wρm. Following we take case: η=0. For η=0, each submanifold of bulk space time is accurate AdS5 space time. Now we study a confined nonminimally coupled scalar field on the brane and consider cosmological aspects. We use energy density and pressure of scalar field (10), and by (14), hence for cosmic expansion we have(15)a¨a=-16α(12ϕ˙2+V(ϕ)-6α′Hϕ˙)×(1+3ϕ˙2-2V(ϕ)+4(α′ϕ¨+2Hα′ϕ˙+α′′ϕ˙2)ϕ˙2+2V(ϕ)-12α′Hϕ˙)-ℓ236α2(12ϕ˙2+V(ϕ)-6α′Hϕ˙)2×(2+3ϕ˙2-2V(ϕ)+4(α′ϕ¨+2Hα′ϕ˙+α′′ϕ˙2)ϕ˙2+2V(ϕ)-12α′Hϕ˙),
where H=a˙/a is Hubble parameter. This is a complex equation; to describe cosmological aspects, we must study some limiting cases or specify α(ϕ), V(ϕ), and ϕ.Equation (15) shows the condition for an accelerating universe depending on the choice of nonminimal coupling and the scalar field potential. Hence we determine the following nonminimal coupling:
(16)α(ϕ)=12(1-ξϕ2).
For an accelerating universe we should have ρm+3pm<0 [12, 34]. If we take V>0, we obtain
(17)(1-3ξ)ϕ˙2-V(ϕ)-3ξϕ(ϕ¨+Hϕ˙)<0.
Dynamics of this scalar field are given as
(18)ϕ¨+3Hϕ˙+ξRϕ+dVdϕ=0.
We obtain
(19)(1-3ξ)ϕ˙2-V(ϕ)+3ξ2Rϕ2+6ξHϕϕ˙+3ξϕdVdϕ<0.
With ρm in this case, we determine
(20)y≡(1-ξϕ2)ρm-2V(ϕ)+ϕ˙2(12-3ξ)+3ξ2Rϕ2+3ξϕdVdϕ<0,
that is, the requirement to have an accelerated universe. In the following, we study the weak energy condition ρm≥0 and limit our consideration to the case by ξ≤1/6. Therefore we obtain(21)-2V+3ξϕdVdϕ<y<0,
and a main condition for cosmic acceleration is
(22)V-3ξ2ϕdVdϕ>0,ξ≤16.
For the model to have a cosmic acceleration by ξ>0, the potential V(ϕ) must change with ϕ  slower than power-law potential Vc(ϕ)=V0(ϕ/ϕ0)2/3ξ (Figure 1). In other hand, when ξ<0, the main condition for cosmic acceleration needs V to grow faster than Vc as ϕ increases [34]. Now we take a simple special example to show how this model processes. By potential of the form V(ϕ)=λϕn,  condition (22) takes λ(1-3nξ/2)>0. In this case just for ξ≤2/3n we have the accelerated expansion.Figure 1
Vvary for different values of the ϕ and nonminimal coupling parameter ξ for V=ϕ2/3ξ.More general case isη≠0 in (13) and (14). In the following we convert the generalized Friedmann equation as
(23)χ˙+4a˙a[χ+112α(12(1-2ξ)ϕ˙2-V(ϕ)-2ξϕ(ϕ¨+3a˙aϕ˙))]=0,
where χ=(a˙/a)2 is dark energy field. Then we take ϕ(t)≈A/tβ  where A is a constant and β an unspecified power will be studied. For scale factor we take the ansatz a(t)≈Btν where B is a constant and ν is to be determined. For scalar field potential we take V(ϕ)=λ(ϕ2/2). Therefore we obtain
(24)a¨a+(a˙a)2≈σ212αu+σ24αuρ(1-3w)+16α(V(ϕ)+Λ)-ϕ˙212α(1+4α′′)-α′3α(ϕ¨+3a˙aϕ˙),ϕ¨+3a˙aϕ˙+u′2uϕ˙2≈-6αuV′+10α′u(V(ϕ)+Λ)+4α′(u+α′u′)σu3ρ(1-3w)+8α′(u+α′u′)u3σ2,
where u=6[1-ξ(1-(16/3)ξ)(ϕ2/2)]. We take that usual matter on the brane has an equation of state such as  ρ=Dt-3ν(1+w). Hence we have
(25)ν(2ν-1)t2≈124{σ23+4Λ+(σ(1-3w)D6)1t3ν(1+w)+A2t2β[2λ+2ξ(Λ+4(3/8-ξ)9σ2)+4ξD(1-3w)(3/8-ξ)σ9×1t3ν(1+w)2βt2×(β(1-8ξ)+4ξ(3ν-1))(Λ+4(3/8-ξ)9σ2)](σ(1-3w)D6)1t3ν(1+w)},
and scalar field equation is given as
(26)β(β+1-3ν)Atβ+2≈-λAtβ[1-83ξ2A2t2β]-5ξ3Atβ[Λ+(λ2+ξ2(1-163ξ)Λ)A2t2β]-ξσ9(2σ+Dt2(1-3w))Atβ×[1+ξ(1-163ξ)(1+ξ)A2t2β],
where β>1. Two conditions must be solved numerically to limit the values of ν and β for an accelerating expansion. After some calculations we obtain
(27)-(1-3w)ξσ9D=Aβ(β+1-3ν),λ=-ξσ212,
where we passed over terms of order 𝒪(t-3β) and higher. Now, if we arrange λ, since ν(2ν-1)=(σ(1-3w)/144)D,  we find
(28)β=3ν-12±[(3ν-1)22-16ξν(2ν-1)]1/2,
where ν=2/3(1+w). For ν<1 the parameter β is real for any value; for ν>1 the exponent β is real only if we have
(29)ξ≤116(3ν-1)22ν(2ν-1).
These equations limited the values of nonminimal coupling ξ to obtain an accelerating expansion (Figure 2). Hence, it seems that nonminimally coupled scalar field identified on the brane world supplies natural candidate for late-time expansion.Variation ofξ with respect to ν with (29) (a) and variation of β with respect to ν and ξ with (28) by sign +, for sign − in (28) β has negative value. We choose special range of parameter ξ with (29) in (b) figure.
(a)
(b)Also we can use (11) to determine dynamical equation of state ω(t).Figure3 shows possible crossing of phantom divided barrier in this framework; this result is important because previous considerations have shown that crossing phantom divided line with a scalar field nonminimally coupled to gravity can exist [35, 36].Figure 3
Dynamical equation of stateω(t) varies for different values of the t and nonminimal coupling parameter ξ by (11) and specific choosing for parameters space.In the following, we explain the physical reasons why the crossing of the phantom divide can occur in this model. Here, we have a scalar field that is nonminimally coupled in gravity and an ordinary matter. Generally, according to [35], these sources of energy momentum can switch together and describe evolution of universe. In this regard, a promising candidate for the dark matter is either a positive cosmological constant or a slowly evolving scalar known as “quintessence.” Also the quintessence can be considered as an alternative method to resolve the cosmological constant problem. This is possible because instead of the fine tuning, the quintessence provides a model of slowly decaying cosmological constant. However, there exists another question for the present-day quintessence. This is to ask whether the expansion will keep accelerating forever or it will decelerate again after some time. This is very similar to the exit problem in the inflationary cosmology. According to the theory of quintessence, the dark energy of the universe is dominated by the potential of a scalar field ϕ which is still rolling to its minimum at V=0.


## 2.2. Rip Singularity

The additional interest to the models with the crossing phantom divideω<-1 is caused by their prediction of a Rip singularity [37]. Theoretically, the scale factor of the universe becomes infinite at a finite time in the future which was dubbed Big Rip singularity. There were proposed several scenarios to cure the Big Rip singularity: (I) to consider phantom acceleration as transient phenomenon. This is possible for a number of scalar potentials; (II) to account for quantum effects which may delay/stop the singularity occurrence [38]; (III) to modify the gravitation itself in such a way that it appears to be observationally-friendly from one side but it cures singularity; (IV) To couple dark energy with dark matter in the special way or to use special (artificial) form for dark energy equation of state [39]. Note that for quintessence dark energy, other (milder) finite-time singularities may occur. For instance, type II (sudden) singularity or type III singularity occurs with finite scale factor but infinite energy and/or pressure. The closer examination shows that the condition ω<-1 is not sufficient for a singularity occurrence [40]. First of all, a transient phantom cosmology is quite possible. Moreover, one can construct such models that ω asymptotically tends to -1 and the energy density increases with time or remains constant, but there is no finite-time future singularity [40]. Of course, most evident case is when Hubble rate tends to be constant (cosmological constant or asymptotically de Sitter space), which may also correspond to the pseudo rip [31]. Very interesting situation is related with Little Rip cosmology [30] where Hubble rate tends to infinity in the infinite future. The key point is that if ω approaches -1 sufficiently rapidly, then it is possible to have a model in which the time required for singularity is infinite. However, this position does not appear in our model, because H is finite.In the following, if we want to consider some types of Rip singularity, we should use Hubble parameter, energy density, and pressure relations that are determined in the previous equations. Therefore, after calculations withξ=0.14 (Figure 4), we obtain results as the following.Variation of pressure (a) and energy density (b) with respect to time for a nonminimal coupling model (whereξ=0.14).
(a)
(b)Case 1 (type II singularity (sudden)).
Conditions: scale factor and energy density reach finite value butp→∞.
In this model, we find that scale factor, pressure, and energy density are finite. Therefore, the model does not predict a sudden singularity.Case 2 (type III singularity).
Conditions, scale factor reaches finite value but (ρ→∞) and |p|→∞.
This means that the energy density grows so rapidly with time and scale factor does not reach the infinite value. The key difference between Cases1 and 2 for this model is that the potential of scalar field has a pole in a case of singularity. Here, according to numerical calculations, if scale factor reaches finite value, energy density and pressure are finite too. Hence, type III singularity does not occur in this setup.Case 3 (pseudo-rip singularity).
Conditions: pressurep→-ρ in infinite time and leads to H˙→0.
This means that the expansion of the universe asymptotically approaches the exponential regime. Obviously in our model with the previous setup, pseudo-rip singularity predict, because in infinite time, whenH˙ reach to zero, p→-ρ.


## 3. A Lorentz Violating Model

In this section with using [13, 14, 41], we consider the cosmological aspects of Lorentz invariance violating model. We emphasize that these model parameters are tightly bounded by both astronomical and cosmological experiments; some works such as [14, 42–44] explained that the model parameters are able to be bounded by comparing the model with cosmological observations. Also a global analysis and a group of canonical values for the model parameters to constrain these parameters are given in [45], so that the model is expected to be viable and background evolution is reasonable.We want to obtain a connection between Lorentz Invariance violation parameter and dynamics of scalar field. First we define an action for a representative scalar-vector-tensor theory which allows Lorentz invariance violation as(30)S=Sg+Su+Sϕ,
where the actions for the vector field Su, the tensor field Sg, and the scalar field Sϕ are assumed as
(31)Sg=∫‍d4x-g116πGR,Su=∫‍d4x-g[-β1∇μuν∇μuν-β2∇μuν∇νuμSu=∫‍d4x-g-β3(∇μuμ)2-β4uμuν∇μuα∇νuαSu=∫‍d4x-g+λ(uμuμ+1)]Sϕ=∫‍d4x-gℒϕ.
This action has a nongravitational degrees of freedom in the structure of Lorentz violating scalar-tensor-vector theory. We take uμuμ=-1 and the expected value of vector field uμ is 〈0|uμuμ|0〉=-1 [45]. βi(ϕ) (i=1,2,3,4) are unrestricted parameters by dimension of mass squared, and also ℒϕ is the Lagrangian density for the scalar field in this model [13, 45]. If universe is homogeneous and isotropic, we consider the universe by metric as
(32)ds2=-𝒩2(t)dt2+e2α(t)δijdxidxj,
where 𝒩 is a lapse function and the universe scale is given by α [13, 14]. If the action varies with respect to metric and selecting a fitting gauge, field equations are obtained as
(33)Rμν-12gμνR=8πGTμν,
where Tμν=Tμν(u)+Tμν(ϕ) is the total energy-momentum tensor, and energy-momentum tensors of vector and scalar fields are given by Tμν(ϕ). The time and space elements of the total energy-momentum tensor are obtained as [14]
(34)T00=-ρu-ρϕ,Tii=pu+pϕ,
where the energy density and pressure of the vector field are determined as
(35)ρu=-3βH2,pu=(3+2H′H+2β′β)βH2,β≡β1+3β2+β3,
a prime point to the derivative of any parameters with respect to α   and   H≡dα/dt=α˙ is the Hubble parameter. The energy equations for scalar field ϕ and the vector field u are given as
(36)ρu′+3(ρu+pu)=+3H2β′,ρϕ′+3(ρϕ+pϕ)=-3H2β′.
So, the total energy equation for the vector and the scalar fields is written as
(37)ρ′+3(ρ+p)=0,(ρ=ρu+ρϕ).
Now, dynamics of the model are determined with the following Friedmann equations [13, 14]:
(38)(1+18πGβ)H2=13βρϕ,(1+18πGβ)(HH′+H2)=-16(ρϕβ+3pϕβ)-H2β′β.
For the scalar section of this model we assume the following Lagrangian:
(39)ℒϕ=-η2(∇ϕ)2-V(ϕ),
where (∇ϕ)2=gμν∂μϕ∂νϕ. Usual scalar fields are matched to η=1 while η=-1 describes phantom fields. The homogeneous scalar field has the density ρϕ and pressure pϕ as
(40)ρϕ=η2H2ϕ′2+V(ϕ),pϕ=η2H2ϕ′2-V(ϕ),
and equation of state parameter is obtained as
(41)ωϕ=pϕρϕ=-1-ηH2ϕ′2/2V1+ηH2ϕ′2/2V.
Now the Friedmann equation takes a form as [14]
(42)H2=13β-[η2H2ϕ′2+V(ϕ)],
where β-=β+1/8πG. With this equation we obtain(43)ϕ′=-2ηβ-(H,ϕH+β-,ϕβ-).
If this equation is substituted into the Friedmann equation, we can find the potential of the scalar field as [14]
(44)V=3β-H2[1-23ηβ-(β-,ϕβ-+H,ϕH)2],
where H=H(ϕ(t)). The equation of state is found in following form:
(45)ωϕ=-1+43ηβ-(H,ϕH+β-,ϕβ-)2=-1+13ηϕ′2β-.
Relations (43) and (45) have main role in following considerations. We can see violation of the Lorentz invariance which has been described by existence of a vector field in the action; now it has been combined in the dynamics of scalar field and equation of state by β- parameter. This situation allows us to consider phantom divided line crossing in LIV model. Now we should solve (43) and (45), to obtain dynamics of scalar field ϕ and the equation of state ωϕ. Therefore, we must define the Hubble parameter H(ϕ) and the vector field, β-(ϕ). In the following, we choose a general case of the Hubble parameter H(ϕ) and the vector field β-(ϕ) and then consider crossing of phantom divided line, late-time acceleration, and Rip singularity scenario.


### 3.1. Considering Late-Time Acceleration in LIV Model

For an acceleration universe that should bea¨>0, we rewrite it to this form H′/H>-1. To consider this case, we must determine some equations. We study a general case for the vector field and the Hubble parameter; note that these equations are a function of ϕ [46]
(46)H=H0ϕζ,β-(ϕ)=mϕn,n>2.
Following, we just study a quintessence scalar field by η=1. By using relation (43), we calculate
(47)ϕ(t)=1A,
where
(48)A=[H0(t-t0)(-4ζm+4ζmn+2ζ2mA=-4mn+2mn2)+ϕ0](1/(n+ζ-2)),
and by using relation (45), we can obtain
(49)ωϕ(t)=-1+43mϕn-2(t)(ζ+n)2.
Following these relations (46), (45), and (49), we find
(50)m2<14(-1)nϕn-2(t)(ζ+n)2,n>2.
This equation describes a restriction in LIV parameters for a late-time acceleration.Now with using (46) and (47) to obtain scale factor as
(51)a(t)=a0(t0)×e(2ϕ0ζ/(n+ζ-2)(m(n+ζ)(n+ζ-2)H0(t-t0)+(1/2)ϕ0)eF-ϕ0)/2ϕ0ζ/(n+ζ-2)(n-2)m(n+ζ),
where(52)F=ζln(1eln(2m(n+ζ)(n+ζ-2)H0(t-t0)+ϕ0)/(n+ζ-2)),
we should select a special space parameter; therefore we approach (51) by a Taylor series in special space parameter
(53)a(t)=.286504+.447663t+.237821t2+.0352651t3.
We can see in Figure 5 acceleration evolution in late time. Now, we study equation of state to determine crossing phantom divided line. By ϕ according to (47), the equation of state is obtained as
(54)ωϕ(t)=-1+43m(ζ+n)2E,
where
(55)E=[H0(t-t0)(-4ζm+4ζmn+2ζ2m-4mn+2mn2)E=+ϕ0]((n-2)/(n+ζ-2)),
which implicitly has a dynamical behavior. This model lets us choose a suitable parameter space to crossing phantom divided line. Moreover, this parameter space must be compared by observational data. We insist in this model that a scalar field and a vector field together can explain crossing phantom divided barrier and late-time acceleration. We can see in Figure 6 the crossing of phantom divided line for a dynamical equation of state. Figure 6 perhaps explains why we are living in an era of   ω<-1. Note that β-(t) has the main role of Lorentz invariance violation model. The dynamical equation for β-(t) has an interested mean: with a suitable fine tuning we can obtain a Lorentz violating cosmology correspond to observational data.Figure 5
Variation of scale factora(t) for different values of t for n=3, ζ=-2 and m=-0.1. The values of ζ are determined by relation (50).Figure 6
ω
ϕ vary for different values of the vector field by m and t for n=3 and ζ=-2. Positive values of ζ do not show a phantom divided line crossing. The values of ζ are obtained with (50).We know many different models that explain phantom divided line crossing, but this model is special because it contains only a scalar field and a Lorentz violating vector field that checks the crossing [47–49]. However, two notices must be explained in this paper: first, in Figures 3 and 6 we can see that there are some sudden jerks of the equation of state. In some models that equation of state crossing the phantom divided line, ω surge around -1 (see [50] and references therein). These jerks are certainly a signature of chaotic manner of equation of state. Second, in those figures we see that crossing of the phantom divided line appears at late time. This means that second cosmological coincidence problem requires extra fine tuning in model parameters.For describing the physical reasons why the crossing of the phantom divide can occur in this model, we can state by a suitable coupling between a quintessence scalar fields and other matter content can lead to a constant ratio of the energy densities of both components which is compatible with an accelerated expansion of the universe or crossing of phantom divide line. In our model, we have three sources of energy momentum:(1) standard ordinary matter, (2) scalar field as a candidate of dark energy, and (3) energy-momentum content depending on Lorentz violating vector field. Here we assume that standard matter has negligible contribution on the total energy-momentum content of the universe. For other two energy-momentum contents, it is possible to use the “trigger mechanism” to explain dynamical equation of state. This means that we assume that scalar-vector-tensor theory contains Lorentz invariance violation which acts like the hybrid inflation models. In this situation, vector and scaler fields play the roles of inflaton and the “waterfall” fields, respectively. Therefore it is reasonable to expect that one of them will eventually dominate to explain inflation or accelerating phase and crossing of phantom divided line.


### 3.2. Rip Singularity

According to Rip singularity scenario that is explained in previous Section2.2, we consider some types of Rip singularity in this LIV model. As shown in Figures 5 and 6, we provide suitable conditions for considering Rip singularity scenario. Hence, corresponding to equation of state and using equations such as energy density, pressure, and Hubble parameter, we can study some Rip singularity solutions.According our setup in this paper that selectedζ=-2 and n=3 in anzats (46) we have
(56)H=H0ϕ-2,β=mϕ3.
Now with substituting previous equations in energy density and pressure relations, we can consider dynamical behaviors in these relations. Note that we selected for ordinary matter ρm=0 and pm=0. Now after numerical calculations (Figure 7), we obtain results as the following.Variation of pressure (a) and energy density (b) with respect to time for a minimal coupling model (whereζ=-2 and n=3).
(a)
(b)In finite time, scale factor and energy density increase, but pressure decreases. Therefore type II singularity (sudden) does not predict. But type III singularity appears because|p| decrease so rapidly. In infinite time, Hubble rate tends to zero; therefore Little Rip singularity is impossible, but pseudo-rip singularity can occur. We emphasize this relations to determine original setup in this paper. This means by different selection in parameters of  ζ and m, all dynamics of relations change. Therefore, one can change setup and consider other options. However, in this situation, our work is limited.


## 3.1. Considering Late-Time Acceleration in LIV Model

For an acceleration universe that should bea¨>0, we rewrite it to this form H′/H>-1. To consider this case, we must determine some equations. We study a general case for the vector field and the Hubble parameter; note that these equations are a function of ϕ [46]
(46)H=H0ϕζ,β-(ϕ)=mϕn,n>2.
Following, we just study a quintessence scalar field by η=1. By using relation (43), we calculate
(47)ϕ(t)=1A,
where
(48)A=[H0(t-t0)(-4ζm+4ζmn+2ζ2mA=-4mn+2mn2)+ϕ0](1/(n+ζ-2)),
and by using relation (45), we can obtain
(49)ωϕ(t)=-1+43mϕn-2(t)(ζ+n)2.
Following these relations (46), (45), and (49), we find
(50)m2<14(-1)nϕn-2(t)(ζ+n)2,n>2.
This equation describes a restriction in LIV parameters for a late-time acceleration.Now with using (46) and (47) to obtain scale factor as
(51)a(t)=a0(t0)×e(2ϕ0ζ/(n+ζ-2)(m(n+ζ)(n+ζ-2)H0(t-t0)+(1/2)ϕ0)eF-ϕ0)/2ϕ0ζ/(n+ζ-2)(n-2)m(n+ζ),
where(52)F=ζln(1eln(2m(n+ζ)(n+ζ-2)H0(t-t0)+ϕ0)/(n+ζ-2)),
we should select a special space parameter; therefore we approach (51) by a Taylor series in special space parameter
(53)a(t)=.286504+.447663t+.237821t2+.0352651t3.
We can see in Figure 5 acceleration evolution in late time. Now, we study equation of state to determine crossing phantom divided line. By ϕ according to (47), the equation of state is obtained as
(54)ωϕ(t)=-1+43m(ζ+n)2E,
where
(55)E=[H0(t-t0)(-4ζm+4ζmn+2ζ2m-4mn+2mn2)E=+ϕ0]((n-2)/(n+ζ-2)),
which implicitly has a dynamical behavior. This model lets us choose a suitable parameter space to crossing phantom divided line. Moreover, this parameter space must be compared by observational data. We insist in this model that a scalar field and a vector field together can explain crossing phantom divided barrier and late-time acceleration. We can see in Figure 6 the crossing of phantom divided line for a dynamical equation of state. Figure 6 perhaps explains why we are living in an era of   ω<-1. Note that β-(t) has the main role of Lorentz invariance violation model. The dynamical equation for β-(t) has an interested mean: with a suitable fine tuning we can obtain a Lorentz violating cosmology correspond to observational data.Figure 5
Variation of scale factora(t) for different values of t for n=3, ζ=-2 and m=-0.1. The values of ζ are determined by relation (50).Figure 6
ω
ϕ vary for different values of the vector field by m and t for n=3 and ζ=-2. Positive values of ζ do not show a phantom divided line crossing. The values of ζ are obtained with (50).We know many different models that explain phantom divided line crossing, but this model is special because it contains only a scalar field and a Lorentz violating vector field that checks the crossing [47–49]. However, two notices must be explained in this paper: first, in Figures 3 and 6 we can see that there are some sudden jerks of the equation of state. In some models that equation of state crossing the phantom divided line, ω surge around -1 (see [50] and references therein). These jerks are certainly a signature of chaotic manner of equation of state. Second, in those figures we see that crossing of the phantom divided line appears at late time. This means that second cosmological coincidence problem requires extra fine tuning in model parameters.For describing the physical reasons why the crossing of the phantom divide can occur in this model, we can state by a suitable coupling between a quintessence scalar fields and other matter content can lead to a constant ratio of the energy densities of both components which is compatible with an accelerated expansion of the universe or crossing of phantom divide line. In our model, we have three sources of energy momentum:(1) standard ordinary matter, (2) scalar field as a candidate of dark energy, and (3) energy-momentum content depending on Lorentz violating vector field. Here we assume that standard matter has negligible contribution on the total energy-momentum content of the universe. For other two energy-momentum contents, it is possible to use the “trigger mechanism” to explain dynamical equation of state. This means that we assume that scalar-vector-tensor theory contains Lorentz invariance violation which acts like the hybrid inflation models. In this situation, vector and scaler fields play the roles of inflaton and the “waterfall” fields, respectively. Therefore it is reasonable to expect that one of them will eventually dominate to explain inflation or accelerating phase and crossing of phantom divided line.


## 3.2. Rip Singularity

According to Rip singularity scenario that is explained in previous Section2.2, we consider some types of Rip singularity in this LIV model. As shown in Figures 5 and 6, we provide suitable conditions for considering Rip singularity scenario. Hence, corresponding to equation of state and using equations such as energy density, pressure, and Hubble parameter, we can study some Rip singularity solutions.According our setup in this paper that selectedζ=-2 and n=3 in anzats (46) we have
(56)H=H0ϕ-2,β=mϕ3.
Now with substituting previous equations in energy density and pressure relations, we can consider dynamical behaviors in these relations. Note that we selected for ordinary matter ρm=0 and pm=0. Now after numerical calculations (Figure 7), we obtain results as the following.Variation of pressure (a) and energy density (b) with respect to time for a minimal coupling model (whereζ=-2 and n=3).
(a)
(b)In finite time, scale factor and energy density increase, but pressure decreases. Therefore type II singularity (sudden) does not predict. But type III singularity appears because|p| decrease so rapidly. In infinite time, Hubble rate tends to zero; therefore Little Rip singularity is impossible, but pseudo-rip singularity can occur. We emphasize this relations to determine original setup in this paper. This means by different selection in parameters of  ζ and m, all dynamics of relations change. Therefore, one can change setup and consider other options. However, in this situation, our work is limited.


## 4. Comparison with Observational Data and Other Dark Energy Models

### 4.1. Analyzing SNe Data

The parameters of the cosmological models can be determined from a strict comparison of their predictions with accurate observational data. Here we consider the data coming from SNe observations, the evolution of the Hubble parameter in the brane model.The modulusμ versus redshift z=a0/a-1 relation corresponding to type Ia supernovae from the Supernova Cosmology Project ([51], Uni 2013 [52]) is, as well known,
(57)μ(z)=μ0+5lgDL(z).
The relation for the luminosity distance DL(z) as a function of the redshift in the FRW cosmology (FC) is
(58)DLFC=cH0(1+z)∫0z‍h-1(z)dz,
where h(z)=[Ωm0(1+z)3+ΩD0F(z)]1/2.Here,Ωm0 is the total fraction of matter density, ΩD0 is the fraction of DE energy density, and H0 is the current Hubble parameter. The constant value μ0 depends on the chosen Hubble parameter:
(59)μ0=42.384-5logh,h=H0/100km/s/Mpc.
The function F(z)=ρD(z)/ρD0 can be determined from the continuity equation
(60)ρ˙D-3a˙ag(ρD)=0,
which can be rewritten as
(61)∫ρD0ρD(z)‍dyg(y)=-3ln(1+z).
For simplicity, we neglect the contribution of radiation. For cosmology on the brane (BC), (58) can be rewritten as
(62)DLBC=cH0(1+z)∫0z‍h-1(z)[1+δh2(z)]-1/2(1+δ)1/2dz,
where for convenience the parameter δ=ρ0/2 has been introduced. For the analysis of the SNe data one needs to calculate the parameter χ2, which is defined by
(63)χSN2=∑i‍[μobs(zi)-μth(zi)]2σi2,
where σi is the corresponding 1σ error. The parameter μ0 is independent of the data points, and, therefore, one has to perform a uniform marginalization over μ0. Minimization with respect to μ0 can be done by simply expanding the χSN2 with respect to μ0,
(64)χSN2=A-2μ0B+μ02C,
where
(65)A=∑i‍[μobs(zi)-μth(zi;μ0=0)]2σi2,B=∑i‍μobs(zi)-μth(zi)σi2,C=∑i‍1σi2.
The expression has a minimum for μ0=B/C at
(66)χ-SN2=A-B2C.
One can minimize χ-SN2 instead of χSN2. With [53] and Table 1, one determines the 56.01% confidence level by Δχ2=χ2-χmin2<1.01 for the one-parameter or 1.9 for the two-parameter model. Similarly, the 91.2% confidence level is determined by Δχ2=χ2-χmin2<3.78 or 5.63 for the one- and two-parameter models, respectively.Table 1
Hubble parameter versus redshift data.
z
H
obs
(
z
)
σ
H
km s−1 Mpc−1
km s−1 Mpc−1
0.170
83
8
1.530
140
14
1.750
202
40


### 4.2. Dark Energy Cosmology

In this paper, we consider two models that explain dark energy in content of the universe, but according to [54] and references therein, there are many different models in this subject. Therefore, in this section for a worthy review, we briefly point out them.As we know a number of popular dark energy models, such as theΛCDM model, Little Rip and Pseudo-Rip, scenarios, the phantom and quintessence cosmologies with the four types (I, II, III, and IV) of the finite-time future singularities and nonsingular universes filled with dark energy have been studied.Investigations show theΛCDM model and recent various cosmological observations to give the bounds on the late-time acceleration of the universe. Furthermore, researchers have considered a fluid description of the universe in which the dark fluid has a general form of the EoS covering the inhomogeneous and imperfect EoS. They have shown that all the dark energy cosmologies can be realized by different fluids and also considered their properties. It has also been demonstrated that at the current stage the cosmological evolutions of all the dark energy universes may be similar to those of the ΛCDM model, and hence these models are compatible with the cosmological observations. In other hand, they have studied the equivalence of different dark energy models, that is, single and multiple scalar field theories, tachyon, scalar theory and holographic dark energy, in which the quintessence/phantom cosmology with the current cosmic acceleration can be represented, and eventually verified those equivalence to the corresponding fluid descriptions. In this regards, as another equivalent class of dark energy models, in which dark energy has its geometrical origins, namely, modifications of gravitational theories, one has examined F(R) gravity including its extension to F(R) Hor˘ava-Lifshitz gravity and f(T) gravity. It has clearly been explored that in these models, the ΛCDM model or the late-time cosmic acceleration with the quintessence/phantom behavior can be performed. It is shown that all the previous dark energy universes may mimic the ΛCDM model currently, consistent with the recent observational data. Furthermore, special attention is paid to the equivalence of different dark energy models.However, it is significant to remark that there are a number of various dark energy models such asF(𝒢) gravity, where 𝒢≡R2-4RμνRμν+RμνρσRμνρσ with Rμν and Rμνξσ being the Ricci tensor Riemann tensors, respectively, is the Gauss-Bonnet invariant, F(R,𝒢) gravity, scalar-Gauss-Bonnet dark energy, k-essence dark energy models, ghost condensates scenario, viscous dark energy, nonminimal derivative dark energy models, G-essence dark energy models, nonlocal gravity produced by quantum effect, which is investigated to account for the coincidence problem of dark energy and dark matter, and galilean dark energy models. In particular, galilean gravity has recently been studied very extensively in the literature. The most important feature of the Lagrangian for the galilean scalar field is that the equation of motion derived from the Lagrangian is up to the second order, so that the appearance of an extra degree of freedom with the existence of a ghost can be avoided. The galilean field originates from a brane bending mode in the Dvali-Gabadadze-Porrati (DGP) brane world scenario, and therefore galilean gravity might be regarded as an indirect resolution for the issue of a ghost in the self-accelerating branch of the DGP model. In all of the previous models, in each model, the difference is only the forms of the energy density ρDE and pressure PDE of dark energy. Hence, the expression of the Hubble parameter H to describe the concrete cosmology, for example, the ΛCDM, quintessence, and phantom cosmologies, can be reconstructed by using gravitational field equations. Similarly, by applying ρDE and PDE, the EoS wDE≡DE/ρDE in the fluid description can also be presented.Finally, we emphasize the role of cosmography in this discussion. As shown, it is a fundamental tool because it allows, in principle, to discriminate among models without a priori assumptions but just laying on constraints coming from data. However, the main criticism to this approach is related to the extension of the Hubble series, the quality and the richness of data samples. In particular, observations cannot be extended at any redshift and, in most of cases, are not suitable to track models up to early epochs. However, the forthcoming observational campaigns should ameliorate the situation removing the degeneration emerging at low redshifts and allowing a deeper insight of models.


## 4.1. Analyzing SNe Data

The parameters of the cosmological models can be determined from a strict comparison of their predictions with accurate observational data. Here we consider the data coming from SNe observations, the evolution of the Hubble parameter in the brane model.The modulusμ versus redshift z=a0/a-1 relation corresponding to type Ia supernovae from the Supernova Cosmology Project ([51], Uni 2013 [52]) is, as well known,
(57)μ(z)=μ0+5lgDL(z).
The relation for the luminosity distance DL(z) as a function of the redshift in the FRW cosmology (FC) is
(58)DLFC=cH0(1+z)∫0z‍h-1(z)dz,
where h(z)=[Ωm0(1+z)3+ΩD0F(z)]1/2.Here,Ωm0 is the total fraction of matter density, ΩD0 is the fraction of DE energy density, and H0 is the current Hubble parameter. The constant value μ0 depends on the chosen Hubble parameter:
(59)μ0=42.384-5logh,h=H0/100km/s/Mpc.
The function F(z)=ρD(z)/ρD0 can be determined from the continuity equation
(60)ρ˙D-3a˙ag(ρD)=0,
which can be rewritten as
(61)∫ρD0ρD(z)‍dyg(y)=-3ln(1+z).
For simplicity, we neglect the contribution of radiation. For cosmology on the brane (BC), (58) can be rewritten as
(62)DLBC=cH0(1+z)∫0z‍h-1(z)[1+δh2(z)]-1/2(1+δ)1/2dz,
where for convenience the parameter δ=ρ0/2 has been introduced. For the analysis of the SNe data one needs to calculate the parameter χ2, which is defined by
(63)χSN2=∑i‍[μobs(zi)-μth(zi)]2σi2,
where σi is the corresponding 1σ error. The parameter μ0 is independent of the data points, and, therefore, one has to perform a uniform marginalization over μ0. Minimization with respect to μ0 can be done by simply expanding the χSN2 with respect to μ0,
(64)χSN2=A-2μ0B+μ02C,
where
(65)A=∑i‍[μobs(zi)-μth(zi;μ0=0)]2σi2,B=∑i‍μobs(zi)-μth(zi)σi2,C=∑i‍1σi2.
The expression has a minimum for μ0=B/C at
(66)χ-SN2=A-B2C.
One can minimize χ-SN2 instead of χSN2. With [53] and Table 1, one determines the 56.01% confidence level by Δχ2=χ2-χmin2<1.01 for the one-parameter or 1.9 for the two-parameter model. Similarly, the 91.2% confidence level is determined by Δχ2=χ2-χmin2<3.78 or 5.63 for the one- and two-parameter models, respectively.Table 1
Hubble parameter versus redshift data.
z
H
obs
(
z
)
σ
H
km s−1 Mpc−1
km s−1 Mpc−1
0.170
83
8
1.530
140
14
1.750
202
40


## 4.2. Dark Energy Cosmology

In this paper, we consider two models that explain dark energy in content of the universe, but according to [54] and references therein, there are many different models in this subject. Therefore, in this section for a worthy review, we briefly point out them.As we know a number of popular dark energy models, such as theΛCDM model, Little Rip and Pseudo-Rip, scenarios, the phantom and quintessence cosmologies with the four types (I, II, III, and IV) of the finite-time future singularities and nonsingular universes filled with dark energy have been studied.Investigations show theΛCDM model and recent various cosmological observations to give the bounds on the late-time acceleration of the universe. Furthermore, researchers have considered a fluid description of the universe in which the dark fluid has a general form of the EoS covering the inhomogeneous and imperfect EoS. They have shown that all the dark energy cosmologies can be realized by different fluids and also considered their properties. It has also been demonstrated that at the current stage the cosmological evolutions of all the dark energy universes may be similar to those of the ΛCDM model, and hence these models are compatible with the cosmological observations. In other hand, they have studied the equivalence of different dark energy models, that is, single and multiple scalar field theories, tachyon, scalar theory and holographic dark energy, in which the quintessence/phantom cosmology with the current cosmic acceleration can be represented, and eventually verified those equivalence to the corresponding fluid descriptions. In this regards, as another equivalent class of dark energy models, in which dark energy has its geometrical origins, namely, modifications of gravitational theories, one has examined F(R) gravity including its extension to F(R) Hor˘ava-Lifshitz gravity and f(T) gravity. It has clearly been explored that in these models, the ΛCDM model or the late-time cosmic acceleration with the quintessence/phantom behavior can be performed. It is shown that all the previous dark energy universes may mimic the ΛCDM model currently, consistent with the recent observational data. Furthermore, special attention is paid to the equivalence of different dark energy models.However, it is significant to remark that there are a number of various dark energy models such asF(𝒢) gravity, where 𝒢≡R2-4RμνRμν+RμνρσRμνρσ with Rμν and Rμνξσ being the Ricci tensor Riemann tensors, respectively, is the Gauss-Bonnet invariant, F(R,𝒢) gravity, scalar-Gauss-Bonnet dark energy, k-essence dark energy models, ghost condensates scenario, viscous dark energy, nonminimal derivative dark energy models, G-essence dark energy models, nonlocal gravity produced by quantum effect, which is investigated to account for the coincidence problem of dark energy and dark matter, and galilean dark energy models. In particular, galilean gravity has recently been studied very extensively in the literature. The most important feature of the Lagrangian for the galilean scalar field is that the equation of motion derived from the Lagrangian is up to the second order, so that the appearance of an extra degree of freedom with the existence of a ghost can be avoided. The galilean field originates from a brane bending mode in the Dvali-Gabadadze-Porrati (DGP) brane world scenario, and therefore galilean gravity might be regarded as an indirect resolution for the issue of a ghost in the self-accelerating branch of the DGP model. In all of the previous models, in each model, the difference is only the forms of the energy density ρDE and pressure PDE of dark energy. Hence, the expression of the Hubble parameter H to describe the concrete cosmology, for example, the ΛCDM, quintessence, and phantom cosmologies, can be reconstructed by using gravitational field equations. Similarly, by applying ρDE and PDE, the EoS wDE≡DE/ρDE in the fluid description can also be presented.Finally, we emphasize the role of cosmography in this discussion. As shown, it is a fundamental tool because it allows, in principle, to discriminate among models without a priori assumptions but just laying on constraints coming from data. However, the main criticism to this approach is related to the extension of the Hubble series, the quality and the richness of data samples. In particular, observations cannot be extended at any redshift and, in most of cases, are not suitable to track models up to early epochs. However, the forthcoming observational campaigns should ameliorate the situation removing the degeneration emerging at low redshifts and allowing a deeper insight of models.


## 5. Summary

Several pieces of evidence from observational data have shown that our universe is in the accelerated expansion era. In this paper, we have shown this reality in two different models: a nonminimally coupling scalar field identified on the brane and a scalar field coupling minimally to gravity in LIV model. According to complex of dynamical equations, we have limited our research to some special form of nonminimal coupling and scalar field potential. Then we have studied special form of time evolution for the scale factor and a scalar field. In brane world model, we have determined an equation that explains situations where an accelerating expansion of the universe is possible. As a consequence, we have obtained which nonminimally coupling scalar field identified on the brane is a suitable candidate for late-time expansion of universe. In LIV model, we have considered a new framework for crossing of phantom divided line with equation of state that combined a violation of Lorentz invariance in a cosmological model. However, we have explained with a suitable choice of parameter space that it is possible to have phantom divided line crossing just with a Lorentz invariance violating vector field and a single scalar field [35, 36]. In this regard, existence of a Lorentz invariance violating vector field prepares a setup for crossing phantom divided line just with a minimally coupling scalar field. In other view point, there is the possibility of a Rip singularity by suitable tuning in the parameters of models.Comparison of nonminimal model with WMAP data gives more accurate constraints on the values of nonminimal coupling. A detailed comparison between our results and some results of previous studies shows that our constraints for nonminimal coupling with exponential potential are consistent with holographic dark energy framework and also with the result of warped DGP brane model [35].



---

*Source: 802685-2013-04-14.xml*