## Posing Constraints on the Free Parameters of a New Model of Dark Energy EoS : Responses Through Cosmological Behaviours

Promila Biswas ∗ 1 , Parthajit Roy ∗∗ 2 and Ritabrata Biswas ∗ 3

- ∗ Department of Mathematics, The University of Burdwan, Golapbag Academic Complex, City : Burdwan-713104, District : Purba Bardhaman, State : West Bengal, Country : India.
- ∗∗ Department of Computer Science, The University of Burdwan, Golapbag Academic Complex, City :

Burdwan-713104, District : Purba Bardhaman, State : West Bengal, Country : India.

## Abstract

Since the late 1990's observations of type Ia Supernova, our universe is predicted to experience a late time cosmic acceleration. Theoretical support to this observation were intended to be built via proposition of a hypothetical fluid which staying inside the universe exerts negative pressure. Cosmologists have prescribed many candidates for this exotic fluid so far. In this alley, a popular method is to choose time dependent equation of state parameter ω = p ρ and to parametrize it as a function of redshift. Again some common families of such parametrizations are constructed among which different members justify different properties of observed universe. Mainly, these were model dependent studies which comprise free parameters to be constrained by different observations. In this present article, a new expression for redshift parametrization is considered and we have constrain its free parameters for two Hubble parameter vs redshift data sets. These data sets are obtained depending on two basic methodologies known as different ages method and baryonic acoustic oscillation method. Different confidence contours for our model are located under the constraints of said data sets. Besides, different thermodynamic parameters related to the evolution of our universe are analysed. It is notified that our model indicates towards a delayed dark matter model which mimics EoS = -1(Λ) phenomena at the present epoch. Deceleration parameters behaviour's are studied. It is noticed that a possibility of future deceleration may occur for this new model. Outcomes for both the data sets are compared with each other.

Keywords : Dark Energy, Scale Factor, Redshift parametrization.

PACS Numbers : 98.80.-k, 95.35.+d, 95.36.+x, 98.80.Jk .

## 1 Introduction

Almost twenty years ago from now, late time cosmic acceleration was pointed out for the first time ever by the observations of two independent supernova collaborating teams (Riess, A.G. et al. 1998; Perlmutter, S. et al. 1999). This has speculated that the cosmic solution is comprised of time independent, spatially homogeneous hypothetical matter density and constant positive space curvature. This need of exotic matter/ energy density, at first, led us to the establishment of cosmological constant, 'Λ' model with Ω m ≈ 0 3, Ω . Λ ≈ 0 7 . . This model is preferred alternative to the Ω m = 1 scenario. Albert Einstein (1917), to build the first cosmological model, introduced cosmological constant term (Λ) to his field equations of general relativity (GR). The motivation was to make the related cosmic solution static. One can treat 'Cosmological Constant' as a constant valued energy density of the vacuum (Zeldovich, Y.B. 1968). As the dynamic cosmological models (Friedmann, A. 1922; Lemaitre, G. 1927) were noted and cosmic expansion was discovered (Hubble, E. 1929), this Λ term appeared as an unnecessary one. Cosmic Microwave Background (CMB) evidence for spatially flat universe (Bernardis, P. de et al. 2000; Hanany, S. et al. 2000), i.e., the proposition for Ω tot ≈ 1 was declared soon after. With Ω m &lt;&lt; 1 and Ω Λ = 0, this eliminated the free expansion theory. Scale factor a t ( ), which is governed by GR, grows at an accelerating rate if the pressure p &lt; -1 3 ρ . Another popular methodology to explain the late time cosmic acceleration was introduced where a hypothetical fluid/ energy component is assumed to exert negative pressure staying at the right hand side of Einstein's field equations. Names of different candidates

1 promilabiswas8@gmail.com

2 roy.parthajit@gmail.com

3 biswas.ritabrata@gmail.com

for such exotic matter are coined as quintessence, phantom or dark energy (DE). Value of the equation of state ω = p ρ is taken to be negative (Peebles, P. J. E. &amp; Ratra, B. 2003; Biswas, P. &amp; Biswas, R. 2018a). Most natural existence of such energy density can be simply obtained if the cosmological constant term is reintroduced into field equations. Cosmological upper bound ( ρ /lessorsimilar 10 71 GeV 4 ) is satisfied by more than the order of hundreds if DE interpolation of Λ is considered (Biswas, P. &amp; Biswas, R. 2018a). Besides this ΛCDM model, some ω CDM models are also famous ( ω is already mentioned DE EoS parameter). Huge discrepancies between observation and theory associated with the vacuum energy density enforces us to construct a better model. Time varying cosmological constant model (Sahni, V. &amp; Starobinsky, A. A. 2000), irreversible process of cosmological matter creation (Ozer, M. &amp; Taha, M. O. 1986), Chaplygin gas family (Biswas, P. &amp; Biswas, R. 2019a), redshift parametrization of the EoS parameters (Biswas, R. &amp; Debnath, U. 2013; Biswas, P. &amp; Biswas R. 2019b) etc are familiar ones. Scalar field φ with potential V ( φ ) (Ratra, B. &amp; Peebles, P. J. E. 1988) is proposed as another model. In the limit 1 2 ˙ φ 2 &lt;&lt; V φ | ( ) , | this model acts like the cosmological constant.

H z ( ) = 1 a t ( ) da t ( ) dt , the expansion rate of the universe is governed by the Friedmann equation. DE density can be written as a function of EoS ω z ( ) as

$$\frac { \rho _ { D E ( z ) } } { \rho _ { D E _ { 0 } } } = \frac { \Omega _ { D E } } { \Omega _ { D E _ { 0 } } } = e x p \left [ 3 \int _ { 0 } ^ { z } \left \{ \frac { 1 + \omega ( z ^ { \prime } ) } { 1 + z ^ { \prime } } \right \} d z ^ { \prime } \right ] \quad,$$

,

Where ρ DE 0 = ρ DE ( z = 0) and Ω DE 0 = Ω DE ( z = 0). In GR based linear perturbation theory, the pressureless matter's density contrast δ /vector , ( x t ) ≡ ρ /vector ,t ( x ) ρ t ( ) -1 grows in proportion to linear growth function G t ( ). This is followed by the differential equation

,

$$\ddot { G } + 2 H ( z ) \dot { G } - \frac { 3 } { 2 } \Omega _ { m } H _ { 0 } ^ { 2 } ( 1 + z ) ^ { 3 } G = 0$$

where '.' denotes differentiation with respect to t .

Upto a good approximation, logarithmic derivative of G z ( ) is

$$f ( z ) \equiv - \frac { d \{ l n G \} } { d \{ l n ( 1 + z ) \} } & \approx \left [ \Omega _ { m } ( 1 + z ) ^ { 3 } \frac { H _ { 0 } ^ { 2 } } { H ^ { 2 } ( z ) } \right ] ^ { \gamma } \quad, \\. \quad \ \. \ \. \ \. \ \. \ \. \ \. \ \. \ \. \ \. \ \. \ \. \ \. \ \. \$$

where relevant values of cosmological parameters (Linder, E.V. 2005) supports γ ≈ 0 55. .

Conventions exist there to phase constrain ω z ( ) in terms of linear evaluation model, ω a ( ) = ω 0 + ω a (1 -a ) = ω p + ω a ( a p -a ), where a = 1 1+ z , ω 0 is the value of ω at z = 0 and ω p is the value of ω at pivot redshift z p = a -1 p -1. More general set of cosmological parameters is constructed. Some necessary parameters are :-

- (i) Dimensionless Hubble's parameter h = H 0 100 Kms -1 Mpc -1 which determines the present day value of critical density and overall scaling of distances inferred from redshifts.
- (ii) Ω m and Ω tot affect expansion history and 'distance-redshift' relations
- (iii) The comoving distances through that the pressure wave is able to propagate between t = 0 and recombination, the sound horizon r s = ∫ t rec 0 c s ( ) t a t ( ) dt , c s is the sound speed through the exotic fluid, determines the physical scale for the Acoustic peaks in the CMB (Pan, Z. 2016) and BAO feature in low redshift matter clustering (Bassett, B. A. &amp; Hlozek, R. 2010).
- (iv) The quantity σ 8 ( z ) represents the amplitude of matter fluctuations and scales the overall amplitude of growth measurements. These may include weak density or redshift-space distortions. Time dependent modelling of DE EoS has a particular alley named as redshift parametrization which can not be obtained from the scale field dynamics. This is because of the fact that these parametrizations are neither limited functions nor they lie in the interval defined by ω = ˙ 2 φ 2 -V ( φ ) ˙2 φ 2 + V ( φ ) with scalar field φ and field potential V ( φ ). Two familiar series of redshift dependent EoS are found, viz,

$$\begin{array} { l } \text{Series I} \colon \omega ( z ) = \omega _ { 0 } + \omega _ { 1 } \left ( 1 - \frac { 1 } { 1 + z } \right ) ^ { n } \text{ and Series III} \colon \omega ( z ) = \omega _ { 0 } + \omega _ { 1 } \frac { z } { ( 1 + z ) ^ { n } }, \, n \in \mathbb { N }. \\ \text{Some other families such as Barboza-Alcaniz (2008), Efstathiou parametrization (Efstathiou, G. 199)!} \end{array}$$

Some other families such as Barboza-Alcaniz (2008), Efstathiou parametrization (Efstathiou, G. 1999; Silva, R. et al 2007), ASSS parametrization (Alam, U. et al 2004a, 2004b), Hannestad M¨rtsell Parametrization (Hannestad, S. o &amp; M¨ ortsell, E. 2004), Lee Parametrization (Lee, S. 2005), Feng Shen Li Li (FSLL) Parametrization (Feng, C., J. et al 2012), Polynomial Parametrization (Weller, J. &amp; Albrecht, A. 2002; Sendra, I. &amp; Lazkoz, R. 2012) can be found in literature.

In this present article, our motive is to propose a new EoS to replicate ω z ( ) = -1 epoch at z = 0. We will try to constrain the free parameters of the model under the observational data. While proposing the DE EoS we take care that it does not fall in any of Series I or Series II. The model has a completely new structure depending on redshift z . We will try to show whether our model generates ω z ( ) = -1 epoch in the neighbourhood of z = 0 or not. Tendency to occur a future deceleration will also be examined. We wish to study the nature of fractional dimensionless density parameters for our model. We plan to study the behaviour of deceleration parameter q z ( ) for the proposed model and whether phase transition-(s) from deceleration to acceleration or the converse take(s) place or not.

The article is organised as follows : In section 2, we construct the cosmological model for our newly proposed parameterization. In section 3, we constrain the model under the data sets obtained with the help of differential ages method and Baryonic Acoustic Oscillation method. Section 4 comprises studies of different cosmological parameters related to our model. Finally, we briefly discuss our findings and conclude in the last section.

## 2 Mathematical Construction of a New Kind of Redshift Parametrization

Einstein's field equation in homogeneous and isotropic Friedmann-Lemaˆ ıtre-Robertson-Walker (FLRW) space time (with flat section) is written as

$$3 \frac { \dot { a } ^ { 2 } } { a ^ { 2 } } = \rho _ { d m } + \frac { 1 } { 2 } \dot { c } \dot { \phi } ^ { 2 } + V ( c \phi ) = \rho _ { d m } + \rho _ { c \phi } + \rho _ { r a d }$$

and

$$2 \frac { \ddot { a } } { a } + \frac { \dot { a } ^ { 2 } } { a ^ { 2 } } = - \frac { 1 } { 2 } \dot { c } \dot { \phi } ^ { 2 } + V ( c \phi ) = -$$

with Planck's units 8 πG = c = 1, as stated earlier, cφ is the scalar field in natural units, ρ cφ , p cφ and V ( cφ ) are the scalar field's density, pressure and potential respectively. ρ dm being the matter density.

Energy density and pressure of scalar field, ρ cφ and p cφ should have the structure as

$$\rho _ { c \phi } = \frac { 1 } { 2 } \dot { c } \phi ^ { 2 } + V ( c \phi ) \quad a n d \quad p _ { c \phi } = \frac { 1 } { 2 } \dot { c } \phi ^ { 2$$

We will write non interactive DE-dark matter conservation equations

$$\dot { \rho _ { c \phi } } + 3 H ( \rho _ { c \phi } + p _ { c \phi } ) = 0 \quad \ f o r \ e n e r g y$$

and

$$\rho _ { d m } + 3 H \rho _ { d m } = 0 \Rightarrow \rho _ { d m } = \rho _ { d m 0 } a ^ { - 3 } \quad f o r \ m a t t e r, \quad \quad \quad$$

where ρ dm 0 denotes the current time value of the matter field's energy density. Equation (4) can be rewritten as

$$\omega _ { c \phi } = \frac { p _ { c \phi } } { \rho _ { c \phi } } = - 1 - \frac { a } { 3 \rho _ { c \phi } } \frac { d \rho _ { c \phi } } { d$$

Among the above equations (1), (2), (4) and (5), only three equations are likely to be independent to each other. Bianchi identities can show the derivation of the fourth equation. So, we are going to solve four independent variables. Without additional input, it is impossible to find an exact solution. We propose an ansatz for the functional form of ρ cφ as

$$\frac { 1 } { \rho _ { c \phi } } \frac { d \rho _ { c \phi } } { d a } = - 3 \left [ \frac { \lambda _ { 1 } } { 1 + a k _ { 1 } } + \frac { \lambda _ { 2 } ( 1 - a ) } { ( 1 + a k _ { 2 } ) ^ { 2 } } \right ] \quad$$

where λ 1 , λ 2 , k 1 and k 2 are constants . 4

4 We propose this ansatz as we need four different equations to obtain an exact solution of four parameters H ρ , m , φ and V ( φ ) , by solving necessarily four equations of which only three we had. While proposing equation (7) , we have tried to increase a degree in 1 1+ z from the CPL parametrization (Chevallier, M., Polarski, D. 2001; Linder, E.V. 2003) retaining the original terms. This might be treated as a modification of JBP parameterization (Jassal, H. K., Bagla, J. S., Padmanabhan, T. 2005) only by changing the constant of the denominator by parameters. We were careful to reproduce the Λ = -1 era at the present time epoch. A particular case of

Integrating, we get

$$\rho _ { c \phi } = A \frac { ( 1 + a k _ { 2 } ) ^ { \frac { 3 \lambda _ { 2 } } { k _ { 1 } ^ { 2 } } } } { ( 1 + a k _ { 1 } ) ^ { \frac { 3 \lambda _ { 1 } } { k _ { 1 } } } } e x p \left \{ \frac { 3 \lambda _ { 2 } ( 1 + k _ { 2 } ) } { k _ { 2 } ^ { 2 } ( 1 + a k _ { 2 } ) } \right \} \quad,$$

where A = ρ cφ 0 (1+ k 1 ) 3 λ 1 k 1 (1+ k 2 ) 3 λ 2 k 2 2 exp { -3 λ 2 k 2 2 } and ρ cφ 0 is the present time (at z = 0) value of the scalar field density. We observe that the density depends on three distinct functions of a . If we make λ 2 = 0 and k 1 = 1, we see the solution will take a simple power law evolution of ρ cφ ( ∼ a -λ ), considered in many cosmological studies (Copeland, E. J. 2006). Equations (7) and (8) together give us the EoS parameter ω cφ as a function of redshift ( z = 1 a -1) as

$$\omega _ { c \phi ( z ) } & = - 1 + \frac { \lambda _ { 1 } } { ( 1 + k _ { 1 } ) + z } + \frac { \lambda _ { 2 } z } { \{ ( 1 + k _ { 2 } ) + z \} ^ { 2 } } \\. \. \. \. \. \. \.$$

This equation even can be treated to be same of

$$\omega _ { c \phi ( z ) } = \omega _ { 0 } + \frac { \omega _ { 1 } } { \omega _ { 2 } + z } + \frac { \omega _ { 3 } z } { ( \omega _ { 4 } + z ) ^ { 2 } } \quad \.$$

In the next section, we will constrain two of the free parameters of our model, namely, λ 1 &amp; λ 2 by the help of Hubble's parameter vs redshift data.

## 3 Constraining the Free Parameters for DA and BAO method :

Equation (10) depicts a new construction of DE EoS parameter. From (1), (5) and (9) we have

$$H ^ { 2 } = H _ { 0 } ^ { 2 } \left [ \Omega _ { r a d 0 } a ^ { - 4 } + \Omega _ { d m 0 } a ^ { - 3 } + \Omega _ { c \phi 0 } \beta \frac { ( 1 + a k _ { 2 } ) ^ { \frac { 3 \lambda _ { 2 } } { k _ { 2 } } } } { ( 1 + a k _ { 1 } ) ^ { \frac { 3 \lambda _ { 1 } } { k _ { 1 } } } } e x p \left \{ \frac { 3 \lambda _ { 2 } ( 1 + k _ { 2 } ) } { k _ { 2 } ^ { 2 } ( 1 + a k _ { 2 } ) } \right \} \right ] \quad,$$

$$\text{where} \, \beta = \frac { ( 1 + k _ { 1 } ) ^ { \frac { 3 \lambda _ { 1 } } { k _ { 1 } } } } { ( 1 + k _ { 2 } ) ^ { \frac { 3 \lambda _ { 2 } } { k _ { 2 } } } } e x p \{ - \frac { 3 \lambda _ { 2 } } { k _ { 2 } } \} \text{ is a constant and } \Omega _ { d m 0 } = \frac { \rho _ { d m 0 } } { 3 H _ { 0 } ^ { 2 } }, \Omega _ { r a d 0 } = \frac { \rho _ { r a d 0 } } { 3 H _ { 0 } ^ { 2 } } \text{ and } \Omega _ { c c 0 } = \frac { \rho _ { c c 0 } } { 3 H _ { 0 } ^ { 2 } } = 1 - \Omega _ { r a d 0 } - \Omega _ { d m 0 } \\ \text{convert} + \text{has a constant and to be a problem}$$

represent the current values of the dimensionless density parameters for the matter, radiation and the scalar field respectively. Now we will proceed for constraining the parametrization's free parameters with Hubble parameter vs redshift data. For this we must enlist the data first and mention the methods for collecting the data.

While we look at the sky, if the celestial distance through which we see some objects is shorter one, Cepheid variables are used as standard candles. When the turn comes for looking at distant galaxies, type Ia Supernova explosions (SNeIa) are taken as standard candles (Riess, A.G. et al. 1998; Perlmutter, S. et al. 1999). Fluctuations in visible baryonic matter's density which is caused by acoustic density waves in primordial plasma of the early universe is used as Standard rulers (Cole, S. et al. 2005; Eisenstein, D. J. et al. 2005). An acoustic wave can travel through primordial plasma until it is cooled to the point where it turns to neutral atoms and the maximum distance traversed by such a way is chosen to be the length of standard ruler. This oscillation's name is coined as baryon acoustic oscillation (BAO). Since the last twenty two years, studies of cosmic microwave background (CMB) (remnant electromagnetic radiation came out of early Big Bang cosmology) have enriched the understanding of expanding universe along with the already mentioned candles and rulers. These methodologies, however, do not directly constrain the Hubble's parameter. Another independent methodology is named as 'cosmic chronometer approach' constructed to constrain the history of the universe's expansion (Jimenez, R. &amp; Loeb, A. 2002; Valent, A. G. &amp; Amendola, L. 2018).

Moresco, M. et al. (2012) have analysed about ≈ 11,000 massive and passive galaxies and enlisted H z ( ) value in the range of 0 15 . &lt; z &lt; 1 1 after eight measurements with accuracy of 5 . -12%. When z &lt; 0 3, we can find the . most accurate constraints. Some authors (Moresco, M. et al 2012b; Wang, X. et al. 2012; Zhao, G. -B.et al. 2012; Riemer-Sorensen, S. et al. 2013) have comparatively discussed the standard probe's (like BAO and SNeIa) and cosmic a constant, The motivation is to examine how good This might be treated

our assumption will lead to a EoS ω a ( ) = ω 0 +( ω 1 + ω 2 ) a -ω a 2 2 ∼ ω Λ + ω a f + ω a s 2 , where ω Λ = -1 , ω f and ω s are two free parameters to be constrained by observation. our model fits with observational data and thermodynamic properties of the universe. to reside in the combination of the first three components of the parameterization 2 ∑ ( n =0 1 -a a 0 ) n .

chronometer. In the works of Moresco, M. et al. (2016), some more Hubble parameter values are enlisted in the range of 0 35 . &lt; z &lt; 0 5. .

A star's age can be calculated when we can analyse the spectra releasing out of it. We can take accumulation of stars, i.e., galaxies so that we can point out the ages. Thus we can observe a clock's behaviour and this 'so called clock' can be observed in archival data (Jimenez, R. et al 2003), GDDS- Gemini Deep Deep Survey (McCarthy, P. J. et al 2004). We use BAO signature density auto-correlation function in size method as the 'standard rod'. We believe that the entire theory is basically standing on the formation of a single 'burst' (Jimenez, R. et al 1999) of almost all stars of a galaxy.

We can notice more sensitive results for ω z ( ) in the differential ages (DA) case. In DA case, we must believe a clock and the dates of which may vary with the universe's age according to redshift. Spectroscopic dating of galaxy's ages provides this clock. We can write /triangle z /triangle t instead of dz dt based on /triangle t and /triangle z , where /triangle z is a small redshift interval and /triangle t is the age difference measurement between two passively reappear galaxies that created at the same time. So it is a more genuine method rather than age determination for galaxies (Dunlop, J. et al. 1996; Alcaniz, J. S. &amp; Lima, J. A. S. 2001; Stockton, A. 2001). The study of globular clusters, absolute stellar ages are more absorbent to well connected ages. Moreover, we can find the lower limit to the age of the universe with respect to absolute galaxy's ages and sort weak limitations of ω z ( ).

We can write Hubble parameter as

.

$$d t = - \frac { ( 1 + z ) ^ { - 1 } } { H ( z ) } d z$$

Applying this for old galaxies, one can obtain the value of H 0 , i.e., using this method to the elliptical galaxies (in the local universe), we can evaluate the current value of Hubble constant.

Now we will form two tables for two types of H z ( ) -z data sets and the corresponding error terms. First 46 data sets are provided by DA method and next, 26 data sets are provided by the standard ruler method, i.e., BAO method.

Table-I : Hubble parameter H z ( ) with redshift and errors σ H from DA method

| Sl No.   | z       | H ( z )   | σ ( z )   | Reference                                                  | Sl No.   |        | z H ( z )   | σ ( z )   | Reference                    |
|----------|---------|-----------|-----------|------------------------------------------------------------|----------|--------|-------------|-----------|------------------------------|
| 1        | 0       | 67.77     | 1.30      | Macaulay, E. et al. 2018                                   | 24       | 0.4783 | 80.9        | 9         | Moresco, M. et al. 2016      |
| 2        | 0.07    | 69        | 19.6      | Zhang, C. et al. 2014                                      | 25       | 0.48   | 97          | 60        | Stern, D. et al. 2010        |
| 3        | 0.09    | 69        | 12        | Simon, J. et al. 2005                                      |          |        |             |           | Alam, S., et al.             |
| 4        | 0.1     | 69        | 12        | Stern, D. et al. 2010                                      | 26       | 0.51   | 90.4        | 1.9       | 2016                         |
| 5        | 0.12    | 68.6      | 26.2      | Zhang, C. et al. 2014                                      |          | 0.57   | 96.8        | 3.4       | Anderson, L. et al.          |
| 6        | 0.17    | 83        | 8         | Stern, D. et al. 2010                                      | 27       |        |             |           | 2014                         |
| 7        | 0.179   | 75        | 4         | Moresco, M. et al. 2012a                                   | 28       | 0.593  | 104         | 13        | Moresco, M. et al. 2012a     |
| 8        | 0.1993  | 75        | 5         | Moresco, M. et al. 2012a; Valent, A. G. & Amendola, L.2018 | 29       | 0.6    | 87.9        | 6.1       | Blake, C. et al. 2012        |
| 9        | 0.2     | 72.9      | 29.6      | Zhang, C. et al. 2014                                      | 30       | 0.61   | 97.3        | 2.1       | Alam, S., et al. 2016        |
| 10       | 0.24    | 79.7      | 2.7       | Gazta˜ naga , E. et al. 2009                               | 31       | 0.68   | 92          | 8         | Moresco, M. et al. 2012a     |
| 11       | 0.27    | 77        | 14        | Stern, D. et al. 2010                                      | 32       | 0.73   | 97.3        | 7         | Blake, C. et al.             |
| 12       | 0.28    | 88.8      | 36.6      | Zhang, C. et al. 2014                                      |          |        |             |           | 2012                         |
| 13       | 0.35    | 82.7      | 8.4       | Chuang, C. H. & Wang, Y. 2013                              | 33       | 0.781  | 105         | 12        | Moresco, M. et al. 2012a     |
| 14       | 0.352   | 83        | 14        | Moresco, M. 2015                                           | 34       | 0.875  | 125         | 17        | Moresco, M. et al.           |
| 15       | 0.38    | 81.5      | 1.9       | Alam, S., et al. 2016                                      |          | 0.88   | 90          | 40        | 2012a Stern, D. et al. 2010  |
| 16       | 0.3802  | 83        | 13.5      | Moresco, M. et al. 2016                                    | 35 36    | 0.9    | 117         | 23        | Stern, D. et al. 2010        |
| 17       | 0.4     | 95        | 17        | Simon, J. et al. 2005                                      | 37       | 1.037  | 154         | 20        | Moresco, M. et al. 2012a     |
| 18       | 0.4004  | 77        | 10.2      | Moresco, M. et al. 2016                                    | 38       | 1.3    | 168         | 17        | Stern, D. et al. 2010        |
| 19       | 0.4247  | 87.1      | 11.2      | Moresco, M. et al.                                         | 39       | 1.363  | 160         | 33.6      | Moresco, M. 2015             |
| 19       | 0.4247  | 87.1      | 11.2      | 2016                                                       | 40       | 1.43   | 177         | 18        | Stern, D. et al. 2010        |
|          | 0.43    | 86.5      | 3.7       | Gazta˜ naga , E. et al. 2009                               | 41       | 1.53   | 140         | 14        | Stern, D. et al. 2010        |
| 20       |         |           |           |                                                            | 42       | 1.75   | 202         | 40        | Stern, D. et al. 2010        |
| 21       | 0.44    | 82.6      | 7.8       | Blake, C. et al. 2012                                      | 43       | 1.965  | 186.5       | 50.4      | Moresco, M. 2015             |
| 22       | 0.44497 | 92.8      | 12.9      | Moresco, M. et al. 2016                                    | 44       | 2.3    | 224.0       | 8.0       | Busca, N. G. et al. 2013     |
| 23       | 0.47    | 89        | 49.6      | Ratsimbazafy, A. L. et al. 2017; Valent,                   | 45       | 2.34   | 222         | 7         | Delubac , T. et al. 2015     |
| 23       | 0.47    | 89        | 49.6      | A. G. & Amendola, L. 2018                                  | 46       | 2.36   | 226         | 8         | Font-Ribera, A., et al. 2014 |

Table-II : Hubble parameter H z ( ) with redshift and errors σ H from BAO method

|   Sl No. |    z |   H ( z ) |   σ ( z ) | Ref. No.                     |
|----------|------|-----------|-----------|------------------------------|
|        1 | 0.24 |     79.69 |      2.99 | Gazta˜ naga , E. et al. 2009 |
|        2 | 0.3  |     81.7  |      6.22 | Oka A. et al. 2014           |
|        3 | 0.31 |     78.18 |      4.74 | Wang Y. et al. 2017          |
|        4 | 0.34 |     83.8  |      3.66 | Gazta˜ naga , E. et al. 2009 |
|        5 | 0.35 |     82.7  |      9.1  | Chuang, C.H. & Wang, Y. 2013 |
|        6 | 0.36 |     79.94 |      3.38 | Wang Y. et al. 2017          |
|        7 | 0.38 |     81.5  |      1.9  | Alam, S., et al. 2016        |
|        8 | 0.4  |     82.04 |      2.03 | Wang Y. et al. 2017          |
|        9 | 0.43 |     86.45 |      3.97 | Gazta˜ naga , E. et al. 2009 |
|       10 | 0.44 |     82.6  |      7.8  | Blake, C. et al. 2012        |
|       11 | 0.44 |     84.81 |      1.83 | Wang Y. et al. 2017          |
|       12 | 0.48 |     87.79 |      2.03 | Wang Y. et al. 2017          |
|       13 | 0.51 |     90.4  |      1.9  | Alam, S., et al. 2016        |

Fig. 1

|   Sl No. |    z |   H ( z ) | σ ( z )   | Ref. No.                     |
|----------|------|-----------|-----------|------------------------------|
|       14 | 0.52 |     94.35 | 2.64      | Wang Y. et al. 2017          |
|       15 | 0.56 |     93.34 | 2.3       | Wang Y. et al. 2017          |
|       16 | 0.57 |     87.6  | 7.8       | Chuang C-H. et al. 2013      |
|       17 | 0.57 |     96.8  | 3.4       | Anderson, L. et al. 2014     |
|       18 | 0.59 |     98.48 | 3.18      | Wang Y. et al. 2017          |
|       19 | 0.6  |     87.9  | 6.1       | Blake, C. et al. 2012        |
|       20 | 0.61 |     97.3  | 2.1       | Alam, S., et al. 2016        |
|       21 | 0.64 |     98.82 | 2.98      | Wang Y. et al. 2017          |
|       22 | 0.73 |     97.3  | 7.0       | Blake, C. et al. 2012        |
|       23 | 2.3  |    224    | 8.6       | Busca, N. G. et al. 2013     |
|       24 | 2.33 |    224    | 8 7       | Bautista, J. E. et al. 2017  |
|       25 | 2.34 |    222    | 8.5       | Delubac , T. et al. 2015     |
|       26 | 2.36 |    226    | 9.3       | Font-Ribera, A., et al. 2014 |

H(z)

z

Fig. 1 : H z ( ) -z graph using the data set from table I and II. Red circle and green rectangle symbol indicate data points of DA method and BAO method respectively.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00408v3-figure-55551181-cb8f-413a-b912-fd7bccffebdd.jpg)

Using the above data for both methods, we will plot H z ( ) -z graphs in figure 1.

We can see the values of H z ( ) corresponding to the values of redshift z are semi increasing. So the resultant graph increases with respect to z . For lower redshift, BAO method estimates the higher values of H z ( ) than DA method. Again, for higher redshift, DA method determines the higher values of H z ( ) as compared to BAO method. Now we will constrain our model's parameters given in equation (10) with the help of the data of table I and II.

In fig.2(a), we have plotted the best fit values for λ 1 and λ 2 along with their corresponding 1 σ , 2 σ , 3 σ confidence contours using the data sets from table I and table II respectively keeping k 1 = k 2 = 0 1. Figure 2( ) uses BAO with . b the data sets to constrain λ 1 and λ 2 . BAO and CMB along with the data sets are given figure 2(c). First we will analyse fig 2(a). For DA method, the point of best fit is denoted by ' P ' and 1 σ , 2 σ , 3 σ confidence regions are drawn in red, blue and black dashed lines respectively. For BAO method, ' Q ' is the best fit and the regions coloured as green, brown and red solid lines for 1 σ , 2 σ , 3 σ confidence contours respectively.

For both the data sets the confidence regions are of more or less elliptic structure and the semi major axis for the regions have the slope more than a right angle. However, slope of confidence contours for DA method is higher than that of BAO method. Tendency wise, high λ 1 is supported with low λ 2 to stay in to particular confidence region

Fig.2(a)-(b) : 1 σ , 2 σ , 3 σ confidence contours for H z ( ) -z data set and H z ( ) -z + BAO data set

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00408v3-figure-c17ed274-22db-4dd9-b83c-77539657441b.jpg)

whereas low value of λ 1 is accompanied with high λ 2 . This signifies that the DE model given in equation (9) allows either of the second or third term to dominate. At the best fits, if we take z = 0, then ω z ( ) = -0 553388 for DA . method and ω z ( ) = -0 891208 for BAO method. . The spans of different regions for both cases are enlisted in table III.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00408v3-figure-099e123f-3933-42f8-8d8a-457a598656ad.jpg)

Comparison shows that the length of major axis for confidence contours of BAO data (table I) is less than that for DA data (table II). So lesser amount of region is enclosed as 1 σ confidence if we consider BAO method. The maximum portion of 1 σ region of BAO is common to that of DA. So BAO method has more tendency to constrain the parameters.

In reference of Eisenstein, D. J. et al. (2005), BAO peak parameter is suggested. BAO signal at a scale of ∼ 100 MPc was detected. Using the BAO peak joint analysis, here we will look into the values of λ 1 and λ 2 for our redshift parametrization model. We will analyse our model in the range of 0 &lt; z &lt; z b , where z b = 2 36 which is called typical redshift (Do-. ran, M. et al 2007) while we use SDSS data samples. BAO peak parameter is defined as

Fig.2(c) : 1 σ , 2 σ , 3 σ confidence contours for H z ( ) -z data set + BAO + CMB

$$& = \sqrt { \Omega _ { m } } \Big ( \frac { \int _ { 0 } ^ { z _ { b } } E ( } { \sqrt { E } ) } \\ & \mathfrak { w } \mod e l, \ \mathcal { A } _ { \mathcal { B }. }$$

$$\mathcal { A } _ { \mathcal { B } A O } = \sqrt { \Omega _ { m } } \left ( \frac { \int _ { 0 } ^ { z _ { b } } E ( z ) ^ { - 1 } d z } { \sqrt { E ( z _ { b } ) } z _ { b } } \right ) ^ { \frac { 3 } { 3 } } \quad.$$

For flat FLRW model, A BAO = 0 469 . ± 0 017. . Hence, for our analysis, the χ 2 BAO can be written as,

$$\chi _ { \mathcal { B } \mathcal { A O } } ^ { 2 } = \frac { ( \mathcal { A } _ { \mathcal { B } \mathcal { A O } } - 0. 4 6 9 ) ^ { 2 } } { 0. 0 1 7 ^ { 2 } } \quad.$$

.

We have plotted the figure 2(b) for H z ( ) -z data sets along with BAO constraint and the addition of it does not change the basic nature of the confidence curves. The best fits and regions are enlisted are in table IV. In this case, ω z ( ) | z =0 is -0 554712 for DA method and . -0 892417 for BAO method. .

$$\mathcal { R } = \sqrt { \Omega _ { m } } \int _ { 0 } ^ { z _ { c } } \frac { d z ^ { \prime } } { E ( z ^ { \prime } ) }$$

CMB power spectrum's shift parameter peak is given by (Efstathiou, G. &amp; Bond J. R. 1999; Elgaroy, O. &amp; Multamaki, T. 2007)

,

where z c describes the value of the particular z at the last scattering surface. WMAP predicts the value of R = 1 726 . ± 0 018 at . z c = 1091 3. The . χ 2 CMB function (for CMB measurement) is defined as

$$\chi _ { e \mathcal { M B } } ^ { 2 } = \frac { ( \mathcal { R } - 1. 7 2 6 ) ^ { 2 } } { 0. 0 1 8 ^ { 2 } } \quad.$$

$$\overline { 0. 0 1 8 ^ { 2 } }$$

We have imposed this method and have plotted the confidence contours in fig 2(c).

The best fits and the spans of confidence contour curves are mentioned in the table V. Secondly, it is noticeable that the best fits are placed in the fourth quadrant. In this case at present epoch, i.e., at z = 0, ω z ( ) = -1 for DA method and -0 893717 for BAO method. .

As ω z ( ) | z =0 = -1 after using BAO and CMB in our model with H z ( ) -z data from DA method, it is significant that our cosmological model with BAO and CMB is the best measurement to support our result with ΛCDM model. We will vary k 1 and k 2 to draw λ 1 -λ 2 confidence contours in figures 2(d)-(g).

Data set obtained by DA method is used to plot fig 2(d) (varying k 1 , k 2 is fixed at 0.1) and fig. 2(e) (varying k 2 , k 1 is fixed at 0.1). If k 2 is fixed and k 1 is increased it is observed that the regions of the confidence contours increase along with a rotation in the semi major axis with positive angles. The corner of low λ 1 high λ 2 is fixed as the center of rotation. This signifies high λ 1 and low λ 2 values are permissible if k 1 is increased. If k 2 is increased and k 1 is kept to be fixed, angle of rotation of the semi major axis of the confidence regions. So increment in k 2 prefer low λ 1 with low λ 2 to stay in the same confidence. In this case, the center of rotation also stays in the low λ 1 -high λ 2 region.

Confidence contours obtained from BAO method's data set are plotted in fig 2(f) ( k 2 = 0 1, varying . k 1 ) and fig 2(g) ( k 2 varying, k 1 = 0 1). . Basic natures do match with 2(d) and 2(e) respectively. Only the difference is the center of rotation for major axis is near by best fit values. So if k 1 is increased BAO data support high λ 1 -low λ 2 regions along with low λ 1 - high λ 2 regions. On the otherhand increment of k 2 supports low λ 1 - low λ 2 or high λ 1 - high λ 2 regions. For both the data sets it is observed that if k 1 and k 2 are low, constraining is more stronger.

Table III :The best fit values of λ 1 , λ 2 , χ 2 and corresponding region of 1 σ , 2 σ and 3 σ for both DA and BAO method using H z ( ) -z data set with k 1 = k 2 = 0 1 .

| Tools            | stat. info   | Region of the contours         | Region of the contours      | Region of the contours        | Region of the contours            |
|------------------|--------------|--------------------------------|-----------------------------|-------------------------------|-----------------------------------|
| H ( z ) - z data | Best fits    | DA Method                      | DA Method                   | BAO Method                    | BAO Method                        |
| H ( z ) - z data | Best fits    | 2 = 44 . 3734 , ω ( z ) ∣ z =0 | ∣ - λ 2 =                   | = 20 . 0144 , ω ( z λ 1 =     | ∣ ∣ z =0 - λ 2 = - 0 . 297328 ( V |
| H ( z ) - z data | Best fits    | λ 1 = 0 . 491273 ( V 11 )      | - 2 . 75766 ( V 12 )        | 0 . 119671 ( V 21 )           | 22                                |
| H ( z ) - z data | 1 σ          | V 11 +0 . 704727 - 0 . 587388  | V 12 +3 . 38546 - 4 . 64334 | V 21 +0 . 712429 - 0 . 566171 | V 22 +3 . 392328 - 4 . 689672     |
| H ( z ) - z data | 2 σ          | V 11 +0 . 846727 - 0 . 684273  | V 12 +3 . 87666 - 5 . 62634 | V 21 +0 . 853029 - 0 . 646471 | V 22 +3 . 867328 - 5 . 600672     |
| H ( z ) - z data | 3 σ          | V 11 +1 . 021727 - 0 . 800673  | V 12 +4 . 49166 - 6 . 85434 | V 21 +1 . 033329 - 0 . 746871 | V 22 +4 . 461328 - 6 . 789672     |

Table IV :The best fit values of λ 1 , λ 2 , χ 2 and corresponding region of 1 σ , 2 σ and 3 σ for both DA and BAO method using H z ( ) -z + BAO data set k 1 = k 2 = 0 1 .

| Tools                  | stat. info   | Region of the contours        | Region of the contours        | Region of the contours        | Region of the contours          |
|------------------------|--------------|-------------------------------|-------------------------------|-------------------------------|---------------------------------|
| H ( z ) - z data + BAO | Best fits    | DA Method                     | DA Method                     | BAO Method                    | BAO Method                      |
| H ( z ) - z data + BAO | Best fits    | 2 = 804 . 712 , ω ( z )       | ∣ ∣ z =0 = - 0 . 554712 λ =   | 2 = 780 . 327 , ω ( z         | ) ∣ ∣ z =0 = - 0 . 892417 λ 2 = |
| H ( z ) - z data + BAO | Best fits    | λ 1 = 0 . 489817 ( V 11 )     | 2 - 2 . 75048 ( V 12 )        | λ 1 = 0 . 118341 ( V 21 )     | - 0 . 289766 ( V 22 )           |
| H ( z ) - z data + BAO | 1 σ          | V 11 +7 . 06183 - 0 . 590517  | V 12 +3 . 373787 - 4 . 614513 | V 21 +0 . 714359 - 0 . 564441 | V 22 +3 . 420766 - 4 . 667234   |
| H ( z ) - z data + BAO | 2 σ          | V 11 +0 . 846183 - 0 . 685317 | V 12 +3 . 914487 - 5 . 618513 | V 21 +0 . 852759 - 0 . 650141 | V 22 +3 . 900766 - 5 . 585234   |
| H ( z ) - z data + BAO | 3 σ          | V 11 +1 . 025183 - 0 . 794917 | V 12 +4 . 454487 - 6 . 853513 | V 21 +1 . 037659 - 0 . 742441 | V 22 +4 . 469766 - 6 . 809234   |

Table V :The best fit values of λ 1 , λ 2 , χ 2 and corresponding region of 1 σ , 2 σ and 3 σ for both DA and BAO method using H z ( ) -z + BAO + CMB data set k 1 = k 2 = 0 1 .

| Tools                        | stat. info   | Region of the contours                 | Region of the contours                 | Region of the contours                          | Region of the contours                          |
|------------------------------|--------------|----------------------------------------|----------------------------------------|-------------------------------------------------|-------------------------------------------------|
| H ( z ) - z data + BAO + CMB | Best fits    | DA Method                              | DA Method                              | BAO Method                                      | BAO Method                                      |
| H ( z ) - z data + BAO + CMB | Best fits    | χ 2 = 9999 . 38 , ω ( z ) ∣ z =0 = - 1 | χ 2 = 9999 . 38 , ω ( z ) ∣ z =0 = - 1 | χ 2 = 9974 . 99 , ω ( z ) ∣ z =0 = - 0 . 893717 | χ 2 = 9974 . 99 , ω ( z ) ∣ z =0 = - 0 . 893717 |
| H ( z ) - z data + BAO + CMB | Best fits    | λ 1 =                                  | ∣ λ 2 =                                | λ 1 =                                           | ∣ λ 2 = 0 . 289766 ( V )                        |
| H ( z ) - z data + BAO + CMB | Best fits    | 0 . 489817 ( V 11 )                    | - 2 . 75048 ( V 12 )                   | 0 . 118341 ( V 21 )                             | - 22                                            |
| H ( z ) - z data + BAO + CMB | 1 σ          | V 11 +0 . 700183 - 0 . 582677          | V 12 +3 . 38758 - 4 . 63252            | V 21 +0 . 720459 - 0 . 564841                   | V 22 +3 . 344766 - 4 . 658234                   |
| H ( z ) - z data + BAO + CMB | 2 σ          | V 11 +0 . 846183 - 0 . 683817          | V 12 +3 . 86148 - 5 . 61952            | V 21 +0 . 861059 - 0 . 645141                   | V 22 +3 . 899766 - 5 . 608234                   |
| H ( z ) - z data + BAO + CMB | 3 σ          | V 11 +1 . 016183 - 0 . 803917          | V 12 +4 . 41448 - 6 . 88352            | V 21 +1 . 041659 - 0 . 745541                   | V 22 +4 . 414766 - 6 . 797234                   |

Table VI :The best fit values of λ 1 , λ 2 , χ 2 from 1 σ , 2 σ and 3 σ confidence contours for both DA and BAO method using H z ( ) -z data set and different k 1 and k 2 values.

-

| DA method   | DA method    | DA method    | DA method   | DA method   | DA method    | DA method      | DA method   |
|-------------|--------------|--------------|-------------|-------------|--------------|----------------|-------------|
| k 2 = 0 . 1 | k 2 = 0 . 1  | k 2 = 0 . 1  | k 2 = 0 . 1 | k 1 = 0 . 1 | k 1 = 0 . 1  | k 1 = 0 . 1    | k 1 = 0 . 1 |
| k 1         | λ 1          | λ 2          | χ 2         | k 2         | λ 1          | λ 2            | χ 2         |
| - 0 . 99    | 0 . 0217069  | - 1 . 32615  | 46 . 0804   | - 0 . 99    | - 0 . 473434 | 0 . 0607605    | 47 . 9418   |
| - 0 . 5     | 0 . 229342   | - 2 . 36384  | 44 . 5134   | - 0 . 5     | 0 . 67951    | - 1 . 28682    | 46 . 7378   |
| 0 . 1       | 0 . 491273   | - 2 . 75766  | 44 . 3734   | 0 . 1       | 0 . 491273   | - 2 . 75766    | 44 . 3734   |
| 1           | 0 . 900393   | - 3 . 03312  | 44 . 4522   | 1           | 0 . 414204   | - 5 . 98958    | 42 . 9764   |
| 2           | 1 . 36386    | - 3 . 18998  | 44 . 577    | 2           | 0 . 379361   | - 10 . 9337    | 42 . 2397   |
| BAO method  | BAO method   | BAO method   | BAO method  | BAO method  | BAO method   | BAO method     | BAO method  |
| 2 = 0 . 1   | 2 = 0 . 1    | 2 = 0 . 1    | 2 = 0 . 1   | k 1 = 0 . 1 | k 1 = 0 . 1  | k 1 = 0 . 1    | k 1 = 0 . 1 |
| k 1         | λ 1          | λ 2          | χ 2         | k 2         | λ 1          | λ 2            | χ 2         |
| - 0 . 99    | 0 . 00478714 | 0 . 103182   | 20 . 2056   | - 0 . 99    | 0 . 113188   | - 0 . 00458441 | 20.0675     |
| - 0 . 5     | 0 . 0548766  | - 0 . 188829 | 20 . 0389   | - 0 . 5     | 0 . 0551597  | 0 . 0382559    | 20 . 0798   |
| 0 . 1       | 0 . 119671   | - 0 . 297328 | 20 . 0144   | 0 . 1       | 0 . 119671   | - 0 . 297328   | 20 . 0144   |
| 1           | 0 . 216863   | - 0 . 354743 | 20 . 0289   | 1           | 0 . 138063   | - 1 . 11864    | 19 . 8405   |
| 2           | 0 . 322417   | 0 . 376157   | 20 . 0515   | 2           | 0 . 143973   | 2 . 49855      | 19 . 694    |

-

## 4 Studies of Different Cosmological Parameters :

The deceleration parameter q is defined as

$$q = - \frac { \ddot { a } } { a H ^ { 2 } } = - \left ( 1 + \frac { \dot { H } } { H ^ { 2 } } \right ) \quad,$$

where ˙ H = dH dt = aH dH da .

From equations (11) and (12), the expression for q in terms of scale factor a can be written as,

$$q ( a ) = - 1 + \frac { 2 \Omega _ { r a d 0 } a ^ { - 4 } + \frac { 3 } { 2 } \Omega _ { d m 0 } a ^ { - 3 } - \frac { 3 a } { 2 } \beta } \, \Omega _ { c c g 0 } \, \frac { ( 1 + a k _ { 2 } ) ^ { \frac { 3 \lambda _ { 2 } } { 4 k _ { 1 } } } } { ( 1 + a k _ { 1 } ) ^ { \frac { 3 \lambda _ { 1 } } { 4 k _ { 2 } } } } \left \{ \frac { \lambda _ { 2 } } { 1 + a k _ { 2 } } - \frac { \lambda _ { 2 } \, ( 1 + k _ { 2 } ) } { ( 1 + a k _ { 2 } ) ^ { 2 } } - \frac { \lambda _ { 1 } } { 1 + a k _ { 1 } } \right \} e x p \left \{ \frac { 3 \lambda _ { 2 } ( 1 + k _ { 2 } ) } { k _ { 2 } ^ { \frac { 1 } { 4 k _ { 2 } } ( 1 + a k _ { 2 } ) } } \right \} } { \Omega _ { r a d 0 } a ^ { - 4 } + \Omega _ { d m 0 } a ^ { - 3 } + \beta \, \Omega _ { c c g 0 } \frac { ( 1 + a k _ { 2 } ) ^ { \frac { 3 \lambda _ { 2 } } { 4 k _ { 1 } } } } { ( 1 + a k _ { 1 } ) ^ { \frac { 3 \lambda _ { 1 } } { 4 k _ { 1 } } } } e x p \left \{ \frac { 3 \lambda _ { 2 } ( 1 + k _ { 2 } ) } { k _ { 2 } ^ { \frac { 1 } { 4 k _ { 2 } } ( 1 + a k _ { 2 } ) } } \right \} } \, \quad ( 1 3 )$$

Fig.2(d)-(g) : 1 σ , 2 σ , 3 σ confidence contours in λ 1 -λ 2 plane for H z ( ) -z data set obtained by DA method (fig.2(d) and fig.2(e)) and BAO method (fig.2(f) and fig.2(g)) for different k 1 (say fig.2(d) and fig.2(e)) and k 2 (say fig.2(f) and fig.2(g)) values.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00408v3-figure-bdfa8bf1-afac-491f-9960-16824a9f7183.jpg)

Now, the equation (13), in terms of redshift z is

$$& \text{Now, the equation (13), in terms of redshift } z \text{ is} \\ q ( z ) = - 1 + \left [ \Omega _ { r a d o } ( 1 + z ) ^ { 4 } + \Omega _ { d m o } ( 1 + z ) ^ { 3 } + \beta \Omega _ { c g o } \left ( \frac { 1 + z } { 1 + k _ { 1 } + z } \right ) ^ { \frac { 3 \lambda _ { 1 } } { k _ { 1 } } } \left ( \frac { 1 + k _ { 2 } + z } { 1 + z } \right ) ^ { \frac { 3 \lambda _ { 2 } } { k _ { 2 } } } e x p \left \{ \frac { 3 \lambda _ { 2 } ( 1 + k _ { 2 } ) ( 1 + z ) } { k _ { 2 } ^ { 2 } ( 1 + k _ { 2 } + z ) } \right \} \right ] ^ { - 1 } \\ & \times \left [ 2 \Omega _ { r a d o } ( 1 + z ) ^ { 4 } + \frac { 3 } { 2 } \Omega _ { d m o } ( 1 + z ) ^ { 3 } - \frac { 3 \beta } { 2 } \Omega _ { c g o } \left \{ \frac { \lambda _ { 2 } } { 1 + k _ { 2 } + z } - \frac { \lambda _ { 2 } } { ( 1 + k _ { 2 } ) ( 1 + z ) } { ( 1 + k _ { 2 } + z ) ^ { 2 } } - \frac { \lambda _ { 1 } } { 1 + k _ { 1 } + z } \right \} \\ & \times \left ( \frac { 1 + z } { 1 + k _ { 1 } + z } \right ) ^ { \frac { 3 \lambda _ { 1 } } { k _ { 1 } } } \left ( \frac { 1 + k _ { 2 } + z } { 1 + z } \right ) ^ { \frac { 3 \lambda _ { 2 } } { k _ { 2 } } } e x p \left \{ \frac { 3 \lambda _ { 2 } ( 1 + k _ { 2 } ) ( 1 + z ) } { k _ { 2 } ^ { 2 } ( 1 + k _ { 2 } + z ) } \right \} \right ] \text{ (14)} \\ & \text{To study the situation in every direction, we analyse the density parameters for the matter field } ( \Omega _ { d m } ) \text{ and scalar field} \\ ( \Omega _ { c o } ) \text{ as},$$

To study the situation in every direction, we analyse the density parameters for the matter field (Ω dm ) and scalar field (Ω cφ ) as,

$$\Omega _ { d m } ( z ) = \frac { \Omega _ { d m 0 } ( 1 + z ) ^ { 3 } } { \Omega _ { r a d 0 } ( 1 + z ) ^ { 4 } + \Omega _ { d m 0 } ( 1 + z ) ^ { 3 } + \beta \, \Omega _ { c c p 0 } \, \left ( \frac { 1 + z } { 1 + k _ { 1 } + z } \right ) ^ { \frac { 3 3 } { k _ { 1 } } } \left ( \frac { 1 + k _ { 2 } + z } { 1 + z } \right ) ^ { \frac { 3 3 } { k _ { 2 } } } \, \exp p \left \{ \frac { 3 3 2 ( 1 + k _ { 2 } ) ( 1 + z ) } { k _ { 2 } ^ { 2 } ( 1 + k _ { 2 } + z ) } \right \} } \, } \quad \ a n d \, \ ( 1 5 )$$

$$\Omega _ { c \phi } ( z ) = \frac { \beta \, \Omega _ { c \phi 0 } \, ( \frac { 1 + z } { 1 + k _ { 1 } + z } ) ^ { \frac { 3 \lambda _ { \alpha } } { k _ { 1 } } } ( \frac { 1 + k _ { 2 } + z } { 1 + z } ) ^ { \frac { 3 \lambda _ { 2 } } { k _ { 2 } } } \ e x p \left \{ \frac { 3 \lambda _ { 2 } ( 1 + k _ { 2 } ) ( 1 + z ) } { k _ { 2 } ^ { 2 } ( 1 + k _ { 2 } + z ) } \right \} } { \Omega _ { r a d 0 } ( 1 + z ) ^ { 4 } + \Omega _ { d m 0 } ( 1 + z ) ^ { 3 } + \beta \, \Omega _ { c \phi 0 } \, \left ( \frac { 1 + z } { 1 + k _ { 1 } + z } \right ) ^ { \frac { 3 \lambda _ { \alpha } } { k _ { 1 } } } \left ( \frac { 1 + k _ { 2 } + z } { 1 + z } \right ) ^ { \frac { 3 \lambda _ { 2 } } { k _ { 2 } } } \ e x p \left \{ \frac { 3 \lambda _ { 2 } ( 1 + k _ { 2 } ) ( 1 + z ) } { k _ { 2 } ^ { 2 } ( 1 + k _ { 2 } + z ) } \right \} } \\ \text{Now add/inc} \, ( 3 ) \text{ and } ( 4 ) \text{, we can obtain.}$$

Now adding (3) and (4), we can obtain,

$$c \dot { \phi } ^ { 2 } = ( 1 + z ) ^ { 2 } H ^ { 2 } \left ( \frac { d c \phi } { d z } \right ) ^ { 2 } \Rightarrow \frac { d c \phi ( z ) } { d z } = \pm \sqrt { 3 } ( 1 + z ) ^ { - 1 } \left [ \frac { \lambda _ { 1 } } { 1 + k _ { 1 } + z } + \frac { \lambda _ { 2 } z } { ( 1 + k _ { 2 } + z ) ^ { 2 } } \right ] ^ { \frac { 1 } { 2 } }$$

$$\times \left [ 1 + \frac { \Omega _ { r a d 0 } ( 1 + z ) ^ { 4 } + \Omega _ { d m 0 } ( 1 + z ) ^ { 3 } } { \Omega _ { c \phi 0 } \beta \left ( \frac { 1 + z } { 1 + k _ { 1 } + z } \right ) ^ { \frac { 3 \lambda _ { 1 } } { k _ { 1 } } } \left ( \frac { 1 + k _ { 2 } + z } { 1 + z } \right ) ^ { \frac { 3 \lambda _ { 2 } } { k _ { 2 } ^ { 2 } } } e x p \left \{ \frac { 3 \lambda _ { 2 } ( 1 + k _ { 2 } ) ( 1 + z ) } { k _ { 2 } ^ { 2 } ( 1 + k _ { 2 } + z ) } \right \} } \right ] ^ { - \frac { 1 } { 2 } } \quad. ( 1 7 ) \\ \intertext { n d } ( 4 ), we can rewrite the potential in terms of that scalar field as,$$

Again from equations (3) and (4), we can rewrite the potential in terms of that scalar field as,

$$V ( c \phi ) = \frac { 1 } { 2 } \rho _ { c \phi } ( 1 - \omega _ { c \phi } ) \quad \.$$

In terms of z , we can express V ( z ) as

$$V ( z ) & = V _ { 0 } \left ( \frac { 1 + z } { 1 + k _ { 1 } + z } \right ) ^ { \frac { 3 \lambda _ { 1 } } { k _ { 1 } } } \left ( \frac { 1 + k _ { 2 } + z } { 1 + z } \right ) ^ { \frac { 3 \lambda _ { 2 } } { k _ { 2 } } } e x p \left \{ \frac { 3 \lambda _ { 2 } ( 1 + k _ { 2 } ) ( 1 + z ) } { k _ { 2 } ^ { 2 } ( 1 + k _ { 2 } + z ) } \right \} \left \{ 1 - \frac { \lambda _ { 1 } } { 2 ( 1 + k _ { 1 } + z ) } - \frac { \lambda _ { 2 } z } { 2 ( 1 + k _ { 2 } + z ) ^ { 2 } } \right \} \quad, \\ \text{show} \quad V _ { 1 } - 2 \, \hat { H } ^ { 2 } \Omega \.. \, \hat { R }$$

where , V 0 = 3 H 2 0 Ω cφ 0 β .

Fig.3(a)-(c) : Ω dm , Ω cφ vs z graphs, where dotted lines show DA method and solid lines show BAO method in the case of H z ( ) -z data, H z ( ) -z data with BAO, H z ( ) -z data + BAO + CMB respectively.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00408v3-figure-6276234b-d0ed-4776-baf5-916cce453a02.jpg)

In figure 3(a)-3(c) we have plotted the fractional dimensionless densities Ω dm for matter and Ω cφ for the exotic matter with respect to redshift z for three of the cases, H z ( ) -z data, H z ( ) -z data+ BAO, H z ( ) -z data + BAO + CMB respectively. We observe that the fractional densities are of increasing nature in past with more or less same slope. Ω dm &gt; Ω cφ for high z . But as time grows, Ω cφ increases and after a certain point, z = z 1 (say), Ω cφ turns greater than Ω dm . These graphs, to some extents, support the theory that in extreme part, the universe was matter dominated with Λ = 0. But as time grows, a delayed decay in matter world took place. This converted matter into relativistic hypothetical energy counterpart and finally Λ = -1 epoch came to exist at the present time. This theory (Turner, M. S. et al 1984; Mathews, G. J. 2008) even helped a lot to bypass different theoretical discrepancies faced by Λ = -1 model alone. This is consistent with the results of Pan, S. &amp; Chakraborty, S. (2013) and Chakraborty, S. et al. (2014). Fig 3(a)-(c) also depict that DA method is more appropriate than BAO to explain the transit from Λ = 0 to -1.

Fig.4(a)-(c) : q vs z graphs, where dotted lines are drawn for DA method and solid lines are drawn for BAO method. From fig 4(a) to 4(c) the graphs are for H z ( ) -z data, H z ( ) -z data + BAO and H z ( ) -z data + BAO + CMB respectively.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00408v3-figure-7e80bc6e-98c7-4c1d-913f-3de67b40f4c3.jpg)

We plot deceleration parameter, q , as a function of z in fig 4(a)-(c). q is found to be a decreasing function of time. The rate of q 's contraction is low at high z and high at low z . We do not find any value z where q changes its sign. So a transition from deceleration to acceleration or the opposite is not allowed for our model. We find at least two values of ( z = z 2 and z 3 , say) where, as the time grows, i.e., z decreases q BAO ( z ) &lt; q DA ( z ) as z &gt; z 3 , then q BAO ( z ) &gt; q DA ( z ) for the interval z 2 &lt; z &lt; z 3 and finally q BAO ( z ) &lt; q DA ( z ) when z &lt; z 2 . So z 2 and z 3 are such two points for which the same deceleration is permitted by both the data sets. This signifies high accelerated expansion is supported by BAO method than DA method in the interval z 2 &lt; z &lt; z 3 . The opposite happens the rest of the z values.

Fig.5(a)-(c) : ω cφ vs z graphs using H z ( ) -z data, H z ( ) -z data + BAO and H z ( ) -z data with BAO +CMB where dotted lines shows DA method and solid line shows BAO method.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00408v3-figure-77cf30f2-ed0b-49c6-984a-7af4ea749a27.jpg)

In figure 5(a)-(c), we plot ω cφ vs z. The whole curve stays in negative zone/ fourth quadrant of z -ω cφ plane. The value of ω cφ is decreasing for a region of high z . As z turns low, ω cφ starts to increase and its value becomes almost equal to -1 at a little past or a small neighbourhood of present time or z = 0. Tendency of ω cφ points that a future deceleration might be possible where the cosmos will transit from phantom to quintessence era again. This result does match with the work of Mamon, A. A. (2018). We find several works Shafieloo, A. et al. (2009), Li, Z. (2010), Magana, J. et al. (2014) where authors have presented that cosmic acceleration is presently witnessing its slowing down. Distinct methods have been used to establish the results. In our model we also find that universe has accelerated its expansion untill recent past. However, in present epoch and preferably in near future, the cosmos is likely to slow down its acceleration. This result is consistent with the above mentioned articles.

Fig 6(a)-(c) are plots of dω cφ dz vs z . The rate of changes of ω cφ with respect to z increases in high z and then it

Fig.6(a)-(c) : Using H z ( ) -z data, D = dω cφ dz vs z graphs are drawn for H z ( ) -z data, H z ( ) -z data + BAO and H z ( ) -z data + BAO + CMB respectively, where dotted lines describes DA method and solid line describes BAO method.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00408v3-figure-80da70b0-1cb3-4217-bfbe-e7e8fb210c51.jpg)

falls near the present time.

## 5 Brief Discussions and Conclusions :

This article comprises of the construction of a redshift dependent model of dark energy and this model's behaviour under the constraints given by two particular redshift-Hubble parameter data sets. We have started with a cosmological model which is mainly governed by two independent components of Einstein's field equations for FLRW metric and equation of continuity for energy and matter. We have noticed that only three among these four governing equations can be independent of each other. But we were to solve four different quantities. This is why the requirement to consider a fourth is followed. This we have done by a process which gave birth of a new equation of state for the corresponding dark energy present in the current cosmos. As we have introduced a new dark energy representative, we require to specify the values of its different parameters. To do so we have motivated ourselves to constrain the parameters for two H z ( ) -z data sets : derived from namely the differential ages method and Baryonic Acoustic Oscillation method. Firstly we plot these two data with each other and observe that the H z ( ) graph is almost increasing with respect to corresponding z . The interesting part of this graph is that the higher values of H z ( ) for low redshift can be seen in BAO case rather than DA method. The similar opposite phenomena happens for higher values of both redshift and H z ( ) in DA method compared to BAO case. Next, we locate the best fit values of two parameters of our model under the data sets obtained by DA method and BAO method (along with BAO scaling and CMB constraints). We have plotted the 1 σ , 2 σ and 3 σ confidence contours for both data sets. We have observed that the contours are elliptic type, the semi major axis of which is inclined with a slope greater than one right angle. Confidence contours for DA method is almost a superset of that for BAO method. So BAO method constrains the model more than DA method does. Increment in the values of k 1 and k 2 increases the span of λ 1 and λ 2 causes the increment of the span of λ 1 -λ 2 confidence contours in different ways. We have noted down the best its of the parameters λ 1 and λ 2 along with the span of the confidence contours in different tables. Fractional dimensionless densities for our model lies in the interval [0 , 1] and fractional density for dark energy increases with time. On the other hand, deceleration parameter's value decreases with time.

The interesting result is found when we check the variation of the equation of state parameter with redshift. EoS parameter decreases with time and then increases again to become equal to almost -1 at the past neighbourhood of present time. Tendency of ω z ( ) shows a possibility future deceleration. From our model, we can theoretically construct the algebraic structure of the deceleration parameter, fraction dimensionless density, dω cφ ( z ) dz etc. We have plotted them as well to understand the deeper insight.

Variations of fractional dimensional densities show quite interesting phenomena. We observe the fractional density of dark energy to start almost from zero at high z and to grow gradually. The same parameter for matter shows completely the opposite behaviour. For low redshift, the fractional density for dark energy grows high and almost becomes asymptotic to unity. This matches with a delayed decay of dark matter into dark energy with time.

Study of deceleration parameter vs redshift does not show any transition from deceleration to acceleration or converse. For present time neighbourhood q falls abruptly. Variation of ω z ( ) shows that at present epoch ω z ( ) is converging to -1, especially when H z ( ) -z data + BAO + CMB is applied as constraining tool.

To conclude in brief, our model which is of inverse quadratic nature fits with H z ( ) vs z data sets derived with the help of different ages method and Baryonic Acoustic Oscillation method by giving the best fits of the ( λ , λ 1 2 ) type as (0.489817, -2.75048) and (0.118341, -0.289766) respectively. Both the data sets indicates to a ω ∼ -1 cosmology in z = 0 epoch. DA does it with a prompt jump than a slower slope of BAO. Probable decay from matter to energy in late time universe is noted.

Acknowledgments: This research is supported by the project grant of Government of West Bengal, Department of Higher Education, Science and Technology and Biotechnology (File no:ST/P/S &amp; T/ 16 G -19 / 2017). PB thanks Department of Higher Education, Science &amp; Technology and Biotechnology, West Bengal for Swami Vivekananda Merit-Cum-Means Scholarship. RB thanks IUCAA, Pune for Visiting Associateship.

## References :

Alam, S., et al. :-'The clustering of galaxies in the completed SDSS-III Baryon Oscillation Spectroscopic Survey: cosmological analysis of the DR 12 galaxy sample' , (2016), [arXiv:1607.03155]

Alam, U., Sahni, V., Saini, T. D. &amp; Starobinski, A. A. :'Is there supernova evidence for dark energy metamorphosis?' , Mon. Not. R. Astron. Soc. , 354 (2004a) 275 [arXiv:0311364].

Alam, U., Sahni, V. &amp; Starobinski, A. A. :'The case for dynamical dark energy revisited' , JCAP , 0406 , (2004b) 008 [arXiv:0403687v2].

Alcaniz, J. S. &amp; Lima, J. A. S. :'Dark Energy And The Epoch Of Galaxy Formation' ApJ , 550 , L133 (2001).

Anderson, L. et al. [BOSS Collaboration] :'The clustering of galaxies in the SDSS-III Baryon Oscillation Spectroscopic Survey: baryon acoustic oscillations in the Data Releases 10 and 11 Galaxy samples.' , Mon. Not. R. Astron. Soc , 441(1) (2014) 24 [arXiv:1312.4877].

Barboza Jr., E. M. &amp; Alcaniz, J. S. :'A parametric model for dark energy' , Phys. Lett. B , 666 , (2008) 415 [arXiv:0805.1713v1].

Bassett, B. A. &amp; Hlozek, R. :'Baryon Acoustic Oscillations' , Dark Energy, Ed. P. Ruiz-Lapuente , (2010) [arXiv:0910.5224v1].

Bautista, J. E. et al. :-'Measurement of baryon acoustic oscillation correlations at z = 2 3 . with SDSSDR12 L yα -Forests' , Astron. Astrophys. , 603 , (2017) 23 [arXiv:1702.00176].

Bernardis, P. de et al.[Boomerang Collaboration] :A Flat Universe from High-Resolution Maps of the Cosmic Microwave Background Radiation' , Nature 404 , (2000) 955 [arXiv:0004404v1].

Biswas, R. &amp; Debnath, U. :'Constraining redshift parametrization parameters of dark energy: loop quantum gravity as background' , Eur.Phys.J.C , 73 (2013) no.5, 2424.

Biswas, P. &amp; Biswas, R. :'Evolution of universe as a homogeneous system: Changes of scale factors with different dark energy Equation of States' , Modern Physics Letters A , 33 , No. 19, (2018) 1850106 [arXiv:1710.06307].

Biswas, P. &amp; Biswas, R. :'Interacting Models of Generalized Chaplygin Gas and Modified Chaplygin Gas with Barotropic Fluid' , Modern Physics Letters A , 34 9(2019a)1950064[arXiv:1805.03962].

Biswas, P. &amp; Biswas, R :'Barboza-Alcaniz Equation of State Parametrization : Constraining the Parameters in Different Gravity Theories' , Modern Physics Letters A , 34 21(2019b)1950163[arXiv:1807.10608].

Blake, C. et al. :-'The WiggleZ Dark Energy Survey: joint measurements of the expansion and growth history at z ≤ 1 .' , Mon. Not. R. Astron. Soc , 425 (2012) 405 [arXiv:1204.3674].

Busca, N. G. et al. :'Baryon acoustic oscillations in the Ly alpha forest of BOSS quasars' , Astron. &amp; Astrophys. , 552 (2013) 18 [arXiv:1211.2616].

Chakraborty, S. et al. :-'A third alternative to explain recent observations: Future deceleration' , Phys. Lett. B , 738 (2014) 424 [arXiv:1411.0941].

Chevallier, M., Polarski, D. :'Acclerating Universes with Scaling Dark Matter ' , Int. J. Mod. Phys. D , 10 , (2001) 213 [arXiv:gr-qc/0009008].

Chuang, C-H. et al. :-'The clustering of galaxies in the SDSS-III Baryon Oscillation Spectroscopic Survey : single-probe measurements and the strong power of f ( z σ ) 8( z ) on constraining dark energy' , Mon. Not.Roy. Astron. Soc. , 433 , (2013). 3559-3571 [arXiv:1303.4486].

Chuang, C.H. &amp; Wang, Y. :'Modeling the Anisotropic Two-Point Galaxy Correlation Function on SmallScales and Improved Measurements of H z ( ) , DA z ( ) , and f ( z ) σ 8( z ) from the Sloan Digital Sky Survey DR7 Luminous Red Galaxies' , Mon. Not. Roy. Astron. Soc. , 435 , (2013) 255-262 [arXiv:1209.0210].

Cole, S. et al. :-'The 2dF Galaxy Redshift Survey: power-spectrum analysis of the final data set and cosmological implications' , Mon. Not. Roy. Astron. Soc. , 362 (2005) 505 [arXiv:0501174].

Copeland, E. J., Sami, M. &amp; Tsujikawa, S. :'Dynamics of dark energy' , Int. J. Mod. Phys. D , 15 (2006) 1753 [arXiv:0603057v3].

Delubac , T. et al. [BOSS Collaboration], 'Baryon acoustic oscillations in the L y forest of BOSS DR11 quasars' , Astron. Astrophys. , 574 (2015) A59 [arXiv:1404.1801].

Doran, M., Stern. S, &amp; Thommes, E. :'Baryon acoustic oscillations and dynamical dark energy' , JCAP , 0704 , (2007) 015 [arXiv:0609075].

Dunlop, J. et al. :-'A 3 5 . -Gyr-old galaxy at redshift 1.55' , Nature , 381 , 581 (1996).

Efstathiou, G. :'Constraining the equation of state of the Universe from Distant Type Ia Supernovae and Cosmic Microwave Background Anisotropies' , Mon. Not. R. Astron. Soc. , 310 , (1999) 842 [arXiv:9904356v1].

Efstathiou, G. &amp; Bond J. R. :'Cosmic confusion: degeneracies among cosmological parameters derived from measurements of microwave background anisotropies' , MNRAS , 304 , (1999) 75 [arXiv:9807103].

Einstein, A. :'Cosmological Considerations in the General Theory of Relativity.' ; Sitzungsber. Preuss. Akad. Wiss. Berlin (Math. Phys.), 142 (1917).

Eisenstein, D. J. et al. (SDSS collaboration):'Detection of the baryon acoustic peak in the large-scale correlation function of SDSS luminous red galaxies' , Astrophys. J. , 633 (2005) 560 [arXiv:0501171].

Elgaroy, O. &amp; Multamaki, T. :'On using the cosmic microwave background shift parameter in tests of models of dark energy' , Astron. Astrophys , 471 (2007) 65E [arXiv:0702343].

Feng, C., J., Shen, X. -Y., Li, P. &amp; Li, X. -Z. :-A New Class of Parametrization for Dark Energy without Divergence' , JCAP , 1209 , (2012) 023, [arXiv:1206.0063].

Font-Ribera, A., et al. :-'Quasar-Lyman α Forest Cross-Correlation from BOSS DR11 : Baryon Acoustic Oscillations' , JCAP , 1405 , (2014) 027 [arXiv:1311.1767].

Friedmann, A. :'On the curvature of space.' ; Z. Phys. 10 , (1922) 377.

Gazta˜ naga , E., Cabr e /primereverse , A. &amp; Hui, L. :'Clustering of luminous red galaxies - IV. Baryonacoustic peak in the lineof-sight direction and a direct measurement of H z ( ) ' , Mon. Not. R.Astron. Soc. , 399

(2009) 1663 [arXiv:0807.3551] Hanany, S. et al. :-'MAXIMA-1: A Measurement of the Cosmic Microwave Background Anisotropy on angular scales of 10 arcminutes to 5 degrees' , Astrophys. J. 545 , (2000) L5 [arXiv:0005123].

Hannestad, S. &amp; M¨rtsell, E. :o 'Cosmological constraints on the dark energy equation of state and its evolution' , JCAP , 0409 , (2004) 001, [arXiv:0407259].

Hubble, E. :'A relation between distance and radial velocity among extra-galactic nebulae.' Proc. Nat. Acad. Sci. 15 , (1929) 168.

Jassal, H. K., Bagla, J. S., Padmanabhan, T. :'WMAP constraints on low redshift evolution of dark energy' , Mon. Not. R. Astron. Soc. , 356 , (2005) L11 333 [arXiv:astro-ph/0404378].

Jimenez, R. et al. :-'Premature dismissal of high-redshift elliptical galaxies' , Mon. Not. R. Astron. Soc. , 305 , (1999) L16-L20 , [arXiv:9812222].

Jimenez, R. et al. :-'Constraints on the equation of state of dark energy and the Hubble constant from stellar ages and the cosmic microwave background' , Astrophys. J , 593 , (2003) 622-629, [arXiv:0302560].

Jimenez, R. &amp; Loeb, A. :'Constraining cosmological parameters based on relative galaxy ages' , Astrophys. J. , 573 (2002) 37 [arXiv:astro-ph/0106145].

Lee, S. :'Constraints on the dark energy equation of state from the separation of CMB peaks and the evolution of alpha ' , Phys. Rev. D , 71 , (2005) 123528 [arXiv:0504650].

Lemaitre, G. :'Un Univers homog` ene de masse constante et de rayon croissant rendant compte de la vitesse radiale des n´buleuses extra-galactiques. e ; Annales de la Soci´t´ Scientifique de Bruxelles e e 47 , (1927) 49.

Li, Z. et al. :-'Probing the course of cosmic expansion with a combination of observational data' , JCAP , 11 (2010) 31 [arXiv:1011.2036].

Lima, J. A. S. &amp; Alcaniz, J. S. :'Flat FRW Cosmologies with Adiabatic Matter Creation: Kinematic tests' , Astron. Astrophys. , 348 , (1999) 1 [arxiv:9902337].

Linder, E. V. :-'Exploring the Expansion History of the Universe' , Phys. Rev. Lett. , 90 (2003) 091301 [arXiv:0208512v1].

Linder, E. V. :'Cosmic growth history and expansion history' , Phys. Rev. D 72 , (2005) 043529 [arXiv:0507263]. Macaulay, E. et. al. (DES collaboration) :'First Cosmological Results using Type Ia Supernovae from the Dark Energy Survey: Measurement of the Hubble Constant' , [arXiv:1811.02376].

Magana, J. et al. :'Cosmic slowing down of acceleration for several dark energy parametrizations' , JCAP , 017 (2014) 1410 [arXiv:1407.1632].

Mathews, G. J., Lan, N. Q. &amp; Kolda, C. :'Late decaying dark matter, bulk viscosity, and the cosmic acceleration' , Phys. Rev. D , 78 , (2008) 043525 [arXiv:0801.0853].

McCarthy, P. J. et al. :-'Evolved galaxies at z &gt; 1 5 . fromthe Gemini Deep Deep Survey: The formation epoch of massive stellar systems' , Astrophys. J , 614 ,(2004) L9-L12, [arXiv:0408367].

Moresco, M. et al. :-'Improved constraints on the expansion rate of the Universe up to z ∼ 1 1 . from the spectroscopic evolution of cosmic chronometers' , JCAP , 08 (2012a) 006 [arXiv:1201.3609].

Moresco, M. et al. :-'New constraints on cosmological parameters and neutrino properties using the expansion rate of the Universe to z ∼ 1 75 . ' , JCAP , 07 (2012b) 053 [arXiv:1201.6658].

Moresco, M. :'Raising the bar: new constraints on the Hubble parameter with cosmic chronometers at z ∼ 2 ' , Mon. Not. Roy. Astron. Soc. , 450 (2015) L16 [arXiv:1503.01116].

Oka, A. et al. :-'Simultaneous constraints on the growth of structure and cosmic expansion from the multipole power spectra of the SDSS DR7 LRG sample' , Mon. Not. Roy. Astron. Soc. , 439 , (2014) 2515-2530 [arXiv: 1310.2820].

Ozer, M. &amp; Taha, M. O. :'A possible solution to the main cosmological problems' , Phys. Lett. B , 171 (1986) 363.

Pan, S. &amp; Chakraborty, S. :'Will there be again a transition from acceleration to deceleration in course of the dark energy evolution of the universe?' , Eur. Phys. J. C , 73 (2013) 2575 [arXiv:1303.5602].

Pan, Z., Knox, L., Mulroe, B. &amp; Narimani, A. :'Cosmic Microwave Background Acoustic Peak Locations' , MNRAS , 459 (2016) 2515 [arXiv:1603.03091v2].

Peebles, P. J. E. &amp; Ratra, B. :- 'The Cosmological Constant and Dark Energy' , Rev. Mod. Phys. , 75 (2003) 559 [arXiv:0207347v2].

Perlmutter, S. et al. (Supernova Cosmology Project collaboration):' Measurements of Ω and Λ from 42 high redshift supernovae' Astrophys. J. 517 (1999) 565 [astro-ph/9812133].

Ratra, B. &amp; Peebles, P. J. E. :'Cosmological consequences of a rolling homogeneous scalar field' , Phys. Rev. D 37 , (1988) 3406.

Ratsimbazafy, A. L. et al. :-'Age-dating Luminous Red Galaxies observed with the Southern African Large Telescope' , Mon. Not. Roy. Astron. Soc. , 467 (2017) 3239 [arXiv:1702.00418].

Riemer-Sorensen, S. et al. :'Simultaneous constraints on the number and mass of relativistic species' , Astrophys. J. , 763 (2013) 89 [arXiv:1210.2131].

Riess, A. G. et al. (Supernova Search Team collaboration) :' Observational evidence from supernovae for an accelerating universe and a cosmological constant' Astron. J. 116 (1998) 1009 [astro-ph/9805201].

Sahni, V. &amp; Starobinsky, A. A. :'The Case for a Positive Cosmological Lambda-term' , Int. J. Mod. Phys. D , 9 (2000) 373 [arXiv:9904398v2].

Sendra, I. &amp; Lazkoz, R. :'SN and BAO constraints on (new) polynomial dark energy parametrizations: current results and forecasts' , Mon. Not. Roy. Astron. Soc. , 422 , (2012) 776 [arXiv:1105.4943].

Shafieloo, A. et al. :-'Is cosmic acceleration slowing down?' , Phys. Rev. D , 80 , (2009) 101301(R) [arXiv:0903.5141].

Silva, R., Alcaniz J. S. &amp; Lima, J. A. S. :- it 'On the thermodynamics of dark energy' , Int. J. Mod. Phys. D 16 , (2007) 469.

Simon, J. et al. :-'Constraints on the redshift dependence of the dark energy potential' , Phys. Rev. D , 71 123001 (2005) [arXiv:0412269].

Stern, D. et al. :-'Cosmic chronometers: constraining the equation of state of dark energy. I: H z ( ) measurements' , JCAP , 1002 (2010) 008 [arXiv:0907.3149].

Stockton, A. :'The Oldest Stellar Populations at z ∼ 1 5 . ' , to be published in Astrophysical Ages and Time Scales, ASP Conference Series, (2001) [arXiv:0104191].

Turner, M. S., Steigman, G. &amp; Krauss, L. M. :'Flatness of the Universe: Reconciling Theoretical Prejudices with Observational Data' , Phys. Rev.Lett. , 52 , (1984) 2090.

Valent, A. G. &amp; Amendola, L. :-' H 0 from cosmic chronometers and Type Ia supernovae, with Gaussian Processes and the novel Weighted Polynomial Regression method' , JCAP , 1804 (2018) 051, [arXiv:2.01505v3].

Wang, X. et al. :-' Observational constraints on cosmic neutrinos and dark energy revisited' , JCAP 11 (2012) 018 [arXiv:1210.2136].

Wang, Y. et al. :-'The clustering of galaxies in the completed SDSS-III BaryonOscillation SpectroscopicSurvey: tomographic BAO analysis of DR12 combined sample configuration space , Mon. Not. Roy.Astron. Soc. , 469 , (2017) 3762-3774 [arXiv:1607.03154].

Weller, J. &amp; Albrecht, A. :'Future supernova observations as a probe of dark energy' , Phys. Rev. D , 65 , (2002) 103512 [arXiv:0106079].

Zeldovich, Y. B. :' Special Issue: the Cosmological Constant and the Theory of Elementary Particles.' ; Soviet Physics Uspekhi 11 , (1968) 381.

Zhang, C. et al. :-'Four new observational H z ( ) data from luminous red galaxies in the Sloan Digital Sky Survey data release seven' , Res. Astron. Astrophys. , 14 , 1221 (2014) [arXiv:1207.4541].

Zhao, G. B. et al. :-'Examining the evidence for dynamical dark energy' , Phys. Rev. Lett. , 109 (2012) 171301 [arXiv:1207.3804].
Promila Biswas Parthajit Roy Ritabrata Biswas
Posing Constraints on the Free Parameters of a New Model of Dark Energy EoS : Responses Through Cosmological Behaviours
Since the late 1990's observations of type Ia Supernova, our universe is
predicted to experience a late time cosmic acceleration. Theoretical support to
this observation were intended to be built via proposition of a hypothetical
fluid which staying inside the universe exerts negative pressure. Cosmologists
have prescribed many candidates for this exotic fluid so far. In this alley, a
popular method is to choose time dependent equation of state parameter $\omega
= \frac{p}{\rho}$ and to parametrize it as a function of redshift. Again some
common families of such parametrizations are constructed among which different
members justify different properties of observed universe. Mainly, these were
model dependent studies which comprise free parameters to be constrained by
different observations. In this present article, a new expression for redshift
parametrization is considered and we have constrain its free parameters for two
Hubble parameter vs redshift data sets. These data sets are obtained depending
on two basic methodologies known as different ages method and baryonic acoustic
oscillation method. Different confidence contours for our model are located
under the constraints of said data sets. Besides, different thermodynamic
parameters related to the evolution of our universe are analysed. It is
notified that our model indicates towards a delayed dark matter model which
mimics EoS = $-1(\Lambda)$ phenomena at the present epoch. Deceleration
parameters behaviour's are studied. It is noticed that a possibility of future
deceleration may occur for this new model. Outcomes for both the data sets are
compared with each other.