Observational constraints on the jerk parameter with the data of the Hubble parameter

We study the accelerated expansion phase of the universe by using the kinematic approach. In particular, the deceleration parameter q is parametrized in a model-independent way. Considering a generalized parametrization for q, we first obtain the jerk parameter j (a dimensionless third time derivative of the scale factor) and then confront it with cosmic observations. We use the latest observational dataset of the Hubble parameter H(z) consisting of 41 data points in the redshift range of 0.07≤z≤2.36\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$0.07 \le z \le 2.36$$\end{document}, larger than the redshift range that covered by the Type Ia supernova. We also acquire the current values of the deceleration parameter q0\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$q_0$$\end{document}, jerk parameter j0\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$j_0$$\end{document} and transition redshift zt\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$z_t$$\end{document} (at which the expansion of the universe switches from being decelerated to accelerated) with 1σ\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$1\sigma $$\end{document} errors (68.3%\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$68.3\%$$\end{document} confidence level). As a result, it is demonstrate that the universe is indeed undergoing an accelerated expansion phase following the decelerated one. This is consistent with the present observations. Moreover, we find the departure for the present model from the standard Λ\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\Lambda $$\end{document}CDM model according to the evolution of j. Furthermore, the evolution of the normalized Hubble parameter is shown for the present model and it is compared with the dataset of H(z).


Introduction
Since the end of last century, various independent observations [1][2][3][4][5][6][7][8][9][10][11] have strongly suggested that the expansion of the universe is speeding up. But, understanding the fundamental physics behind this accelerated expansion is still an open question in modern cosmology. In order to give a reasonable explanation to this accelerating scenario, a large variety of attempts have been done. These attempts include the modified gravity models which relates to the changes of the geometry a e-mail: abdulla.physics@gmail.com b e-mail: bamba@sss.fukushima-u.ac.jp of the spacetime, and the dark energy models which involve the introduction of exotic matter sources (for a review, one can look into Refs. [12][13][14][15][16][17][18][19]). In general, these models correspond to the dynamics of the universe. Although these models fit the observational data, but they also have their own demerits. For example, the CDM model is the most natural one which shows very well consistence with the various observational data, however, it can not escape from the fine tuning and cosmological coincidence problems [20,21]. So, the study of explaining the cosmic acceleration is still continued.
In cosmology, another way to understand the cosmic acceleration is to analyze kinematic variables like the Hubble parameter (H ), the deceleration parameter (q), or the jerk parameter ( j), which are all derived from the derivatives of the scale factor (for details, see Sect. 2). The kinematic approach is advantageous since it does not need any model specific assumptions like the composition of the universe. It is described by a metric theory of gravity and is assumed that the universe is homogeneous and isotropic at cosmological scales (for review on this topic, see [22][23][24][25][26][27][28][29][30][31]). In the literature, there have been many attempts to constrain the present values of H , q and j by parametrizing q or j [31][32][33][34][35][36][37][38][39][40][41][42]. For example, Rapetti et al. [29] performed a systematic study of jerk parameter as the way towards building up a model in order to examine the expansion history of the universe. On the other hand, by setting the CDM model as the fiducial model and using the Type Ia Supernova and observational Hubble parameter data, recently Zhai et al.
[41] constrained four jerk models through different parametrizations of j ( j = CDM value + departure), as a function of the redshift z. In a pioneering work, Riess et al. [43] measured a transition from an early decelerating to present accelerating phase using a simple linear redshift parameterization of q (q(z) = q 0 + q 1 z). However, this model is not reliable at high redshift. In a recent work, Xu et al. [39] studied few kinematic models by considering the linear first order expansion of q (q(z) = q 0 + q 1 z 1+z ), constant jerk parameter and third order expansion of luminosity distance. Following the same line of thought, our main goal in this paper is to examine some simple kinematic model for the cosmic expansion based on a general parameterization for q. The features of this parametrization have been discussed in the next section. We have then derived the expression of j (z) for this specific choice of q. However, our work is more general and also different from other similar works [39,43] in different ways. Firstly, similar to the work of Zhai et al.
[41], we do not assume here a flat CDM model for the present universe a priori, but rather allow our model to behave in a more general way. Additionally, the present value of the jerk parameter is allowed to be fixed by the observational data. Secondly, in this work, we go one step further by studying the evolution of jerk parameter for a general deceleration parameter, which is independent of the matter content of the universe. Lastly, here, we employ the latest H (z) dataset as useful cosmic constraints. The paper is organized as the following. In Sect. 2, we have described the phenomenological model considered here. In Sect. 3, we have described the observational data used in the present work along with the statistical analysis and discussed the results in Sect. 4. Finally, we have summarized the main conclusions in Sect. 5.

The Kinematic model
In what follows, we have assumed a homogeneous and spatially flat Friedmann-Robertson-Walker (FRW) metric described by the line element: where a(t) is the cosmic scale factor (which is scaled to be unity at the current epoch, i.e., a 0 = 1) and t is the cosmic time.
As discussed in the previous section, the Hubble parameter, deceleration parameter and jerk parameter are purely kinematical, since they are independent of any gravity theory, and all of them are only related to scale factor a or redshift z (since, a = 1 1+z ). In particular, the jerk parameter, a dimensionless third derivative of the scale factor a(t) with respect to cosmic time t, can provide us the simplest approach to search for departures from the concordance CDM model. It is defined as [28,29,[44][45][46].
where H =ȧ a is the Hubble parameter and the "dot" implies derivative with respect to t. In terms of the deceleration parameter q (a dimensionless second derivative of a(t) with respect to t), the jerk parameter j can be written as where q = −ä a H 2 . Blandford et al. [28] described how the jerk parameterization provides an alternative and a convenient method to describe cosmological models close to concordance CDM model. A powerful feature of j is that for the CDM model j = 1 (constant) always. It should be noted here that Sahni et al. [47,48] drew attention to the importance of j for discriminating different dark energy models, because any deviation from the value of j = 1 (just as deviations from the equation of state parameter ω = −1 do in more standard dynamical approaches) would favor a non-CDM model. The simplicity of the jerk formalism thus enables us to constrain the departure from the CDM value in an effective manner. Also, the Eq. (3) is useful when the parametric form of the deceleration parameter q(z) is given. In fact, a variety of q-paramertized models have been proposed in the literature (for details see [30,31,43,[49][50][51][52][53][54][55][56][57][58][59][60]). Following this line of thought, in the present work, we are interested to investigate the evolution of j for a general q-parametrized model given in [55]. It is given by where q 0 , q 1 and α are arbitrary model parameters. In Eq. (4), q 0 indicates the present value of q, and q 1 indicates the derivative of q(z) with respect to the redshift z. From Eq. (4), one can easily recover few popular q-parametrized models in the following limits [43,49,50]: It is worth noting here that the generalized parametrization is not valid at α = 0. The interesting cosmological characteristics of the parametrization, as given in Eq. (4), are extensively discussed in [55]. Our main goal in this paper is to examine some simple kinematic model for the cosmic expansion based on the parameterization for q(z) in Eq. (4). With this choice of q(z), the expression for the Hubble parameter is obtained as where H 0 denotes the present value of the Hubble parameter. In this case, the transition redshift (whereä(t) vanishes) can be obtained as With the corresponding q-parametrization, the jerk parameter j (z), defined in Eq. (3), is obtained as . Since the expression for the jerk parameter is explicit, so we can think that we are actually parameterizing j (z) instead of q(z). The advantage of this type of jerk parametrization is that it incorporates a wide class of viable models of cosmic evolution based on the choice of the parameter α. Similar to the generalized qparametrization, the above jerk parametrization is also not valid at α = 0. As mentioned in the previous section, Zhai et al.
[41] parameterized j (z) phenomenologically aiming at measuring the departure of j from the CDM value. An important difference with the work of Zhai et al.
[41] is that we do not assume a flat CDM model for the present universe a priori, but rather allow our model to behave in a more general way. Also, the value of j 0 is allowed to be fixed by the observational dataset. Obviously, the cosmological characteristics of the model given in Eq. (4) (or Eq. (8)) strongly depend on values of the parameters q 0 , q 1 and α. Using the latest H (z) data, in the next section, we have constrained the parameters q 0 and q 1 for some specific values of α.

Observational constraints on the model parameters
In this section, we have described the latest observational data used in our analysis and the method used to analyze them.
It is known that the Type Ia Supernova or CMB or BAO dataset is powerful in constraining the cosmological models. However, the integration in its formula makes it hard to reflect the precise measurement of the expansion rate of the universe as a function of redshift, i.e., H (z) [61]. This is the most direct and model independent observable of the dynamics of the universe. Therefore, the fine structure of the expansion history of the universe can be well indicated by the H (z) dataset. From the observational point of view, the ages of the most massive and passively evolving galaxies will provide direct measurements of H (z) at different redshifts, which develop another type of standard probe (namely, standard clocks) in cosmology [62]. It should be noted that H (z) measurements are always obtained from two different techniques: galaxy differential age (also known as cosmic chronometer) and radial BAO size methods. The Hubble parameter depending on the differential ages as a function of redshift z can be written in the form of So, H (z) can be obtained directly if dz dt is known [63]. Also, the apparently small uncertainty of this measurement naturally increases its weight in the χ 2 statistics. For this dataset, the χ 2 is defined as where H obs is the observed Hubble parameter at z i and H th is the corresponding theoretical value given by Eq. (6) [67]. In our analysis, we have also used the present value of Hubble parameter H 0 , determined from the combined analysis with Planck+highL+WP+BAO [68]. The best fit values of the model parameters (say, θ * ) from the Hubble data are estimated by minimizing the χ 2 function in equation (10). It should be noted that the confidence levels 1σ (68.3%) and 2σ (95.4%) are taken proportional to χ 2 = 2.3 and 6.17 respectively, where χ 2 = χ 2 (θ )−χ 2 (θ * ) and χ 2 m is the minimum value of χ 2 . An important quantity which is used for data fitting process is where subscript "dof" is the abbreviation of degree of freedom, and it is defined as the difference between all observational data points and the number of free parameters. If χ 2 m dof ≤ 1, then the fit is good and the observed data are consistent with the considered model.  Table 1 Best fit values of q 0 and q 1 with 1σ error bars obtained via χ 2 minimization method. Also, the values of z t and j 0 are given for the best fit model

Results
Following the χ 2 analysis (as described in the previous section), we have obtained the constraints on the model parameters q 0 and q 1 by fixing the parameter α to some constant values (1, 0.5, 0.3, 0.007, − 0.5 and − 1) for the latest H (z) dataset. The 1σ and 2σ contours in q 0 − q 1 plane for the proposed model is shown in Fig. 1. The best fit values with 1σ errors of the parameters q 0 , q 1 , z t and j 0 are displayed in Table 1. It is clear from Table 1 that the model has almost same goodness in the viewpoint ofχ 2 for different values of α. But, there the values of q 0 , z t and j 0 are different. The plot of the deceleration parameter q(z), as given in Fig. 2, clearly shows that our model successfully generates late time cosmic acceleration along with a decelerated expansion in the past for all values of α. Figure 2 also shows that the transition from decelerated to accelerated expansion phase took place in the redshift range 0.685 ≤ z ≤ 0.974 within 1σ errors for different values of α (for details, see fifth column of Table 1). This range is well consistent with those previous imposed the jerk parameter to mimic as the flat CDM model via their parametrization at z = 0, but our work relaxes that requirement. In a word, the constraint results of j 0 tend to favor the dynamical jerk parameter.
In Fig. 4, we have shown the evolution of the normalized Hubble parameter h(z) = H (z) H 0 for our model and have compared that with the latest 41 points of H (z) dataset [64,65]. We have also plotted data points for h(z) with 1σ error bars which have been obtained from the H (z) dataset using the current value of H (z) given by Planck observations [68]. The corresponding error in h can be estimated as [53] where σ H 0 and σ H are the errors in H 0 and H measurements respectively. We have observed from Fig. 4 that the our model is well consistent with the H (z) data against redshift parameter for different values of α (except α = −1 case). The reason is simple that the model with α = −1 [i.e., q(z) = q 0 + q 1 z, see Eq. (5)] is not reliable at high redshift.

Conclusions
In this paper, we have studied the accelerated expansion phase of the universe by using the kinematic approach. In this context, we have parameterized the deceleration parameter q in a model independent manner. The functional form of q is chosen in such a way that it reproduces three popular q-parametrized models, such as q ∝ z, q ∝ ln(1 + z) and q ∝ z 1+z for α = −1, α → 0 and α = 1 respectively. Consequently, the jerk parameter j also incorporates a wide class of viable models of cosmic evolution based on the choice of the parameter α. We have also constrained the model parameters by χ 2 minimization technique using the latest 41 points of H (z) dataset. Figure 1 shows the 1σ and 2σ confidence level contours in the q 0 -q 1 parametric plane for different choices of α. The best fit values of the model parameters, transition redshift z t and j 0 within 1σ error regions for different values of α are displayed in the Table 1. Theχ 2 implies the same goodness of the model for all the values of α considered here. We have also shown the evolution of the normalized Hubble parameter for our model and have compared that with the latest H (z) dataset. In what follows, we have summarized our main conclusions in more detail: For all the values of α considered here, the kinematic model shows a smooth transition from the deceleration (q > 0) phase to acceleration (q < 0) phase of the universe in the recent past. It has been found that values of the transition redshift z t (from decelerated to accelerated expansion) and q 0 depend upon the choice of α . However, the changes in the values of z t and q 0 do not differ by very large values. It has also been found that the values of z t obtained in this work for a wide range in the values of α, are in good agreement with the previous results as reported in [29,31,53,54,[69][70][71][72][73]. We have found from Fig. 3 and Table 1 that the CDM model is not well supported within 1σ confidence level at the present epoch, except for the cases α = 0.5 and 0.3. As discussed earlier, the present model is allowed to pick up any values of j depending on the parameters to be fixed by the observed data as contrary to the work of Zhai et al. [41], where j is constrained to mimic a flat CDM model at z = 0. It has also been found that our model is well consistent with the H (z) data at the low redshifts for different choices of α (Fig. 4). Therefore, we conclude the present H (z) data provides well constrained values of j and our model remains at a very close proximity of the standard CDM model. However, it is natural to extend the present work with addition of Hubble parameter dataset from the GW standard sirens of neutron star binary system [74].