arxiv: v2 [gr-qc] 19 Oct 2018

Transcription

1 Observational constraints on the jerk parameter with the data of the Hubble parameter arxiv:85.85v [gr-qc] 9 Oct 8 Abdulla Al Mamon, and Kauharu Bamba, Department of Mathematics, Jadavpur University, Kolkata-7, India Division of Human Support System, Faculty of Symbiotic Systems Science, Fukushima University, Fukushima 96-96, Japan We study the accelerated expansion phase of the universe by using the kinematic approach. In particular, the deceleration parameter q is parametried in a model-independent way. Considering a generalied parametriation 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() consisting of data points in the redshift range of.7.6, larger than the redshift range that covered by the Type Ia supernova. We also acquire the current values of the deceleration parameter q, jerk parameter j and transition redshift t (at which the expansion of the universe switches from being decelerated to accelerated) with σ errors (68.% 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 ΛCDM model according to the evolution of j. Furthermore, the evolution of the normalied Hubble parameter is shown for the present model and it is compared with the dataset of H(). PACS numbers: 98.8.Hw Keywords: deceleration parameter, jerk parameter, kinematical approach, cosmic chronometer I. INTRODUCTION Since the end of last century, various independent observations [ ] 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 of the spacetime, and the dark energy models which involve the introduction of exotic matter sources (for a review, one can look into Refs. [ 9]). 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 [, ]. So, the study of explaining the cosmic acceleration is still continued. In cosmology, another way to understand the cosmic acceleration is to analye 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 section II). 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 [ ]). In the literature, there have been many attempts to constrain the present values of H, q and j by parametriing q or j [ ]. For example, Rapetti et al. [9] 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. [] constrained four jerk models through different parametriations of j (j = ΛCDM value + departure), as a function of the redshift. In a pioneering work, Riess et al. [] measured a transition from an early decelerating to present accelerating phase using a simple linear redshift parameteriation of q (q() = q +q ). However, this model is not reliable at high redshift. In a recent work, Xu et al. [9] studied few kinematic models by considering the linear first order expansion of q (q() = q + q + ), constant jerk parameter and third order Electronic address: abdulla.physics@gmail.com Electronic address: bamba@sss.fukushima-u.ac.jp

2 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 parameteriation for q. The features of this parametriation have been discussed in the next section. We have then derived the expression of j() for this specific choice of q. However, our work is more general and also different from other similar works [9, ] in different ways. Firstly, similar to the work of Zhai et al. [], 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() dataset as useful cosmic constraints. The paper is organied as the following. In section II, we have described the phenomenological model considered here. In section III, we have described the observational data used in the present work along with the statistical analysis and discussed the results in section IV. Finally, we have summaried the main conclusions in section V. II. THE KINEMATIC MODEL In what follows, we have assumed a homogeneous and spatially flat Friedmann-Robertson-Walker (FRW) metric described by the line element: ds = dt a (t)[dr +r dω ] () where a(t) is the cosmic scale factor (which is scaled to be unity at the current epoch, i.e., a = ) 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 (since, a = + ). 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 [8, 9, 6] j(a) = d a dt ah () 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 [ j(q) = q(q +)+(+) dq ] () d where, q = ä ah. Blandford et al. [8] described how the jerk parameteriation 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 = (constant) always. It should be noted here that Sahni et al. [7, 8] drew attention to the importance of j for discriminating different dark energy models, because any deviation from the value of j = (just as deviations from the equation of state parameter ω Λ = 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 equation () is useful when the parametric form of the deceleration parameter q() is given. In fact, a variety of q-paramertied models have been proposed in the literature (for details see [,,, 9 6]). Following this line of thought, in the present work, we are interested to investigate the evolution of j for a general q-parametried model given in [57]. It is given by [ (+) α ] q() = q q α () where q, q and α are arbitrary model parameters. In equation (), q indicates the present value of q, and q indicates the derivative of q() with respect to the redshift. From equation (), one can easily recover few popular

3 q-parametried models in the following limits [, 9, 5]: q +q, for α = q() = q +q ln(+), for α. (5) ( ) q +q, for α = + + It is worth noting here that the generalied parametriation is not valid at α =. The interesting cosmological characteristics of the parametriation, as given in equation (), are extensively discussed in [57]. Our main goal in this paper is to examine some simple kinematic model for the cosmic expansion based on the parameteriation for q() in equation (). With this choice of q(), the expression for the Hubble parameter is obtained as [ +q(x) ] [ H() = H exp +x dx = H (+) (+q+q ) q {(+) α } ] α exp α (6) where H denotes the present value ofthe Hubble parameter. In this case, the transition redshift (where ä(t) vanishes) can be obtained as ( q ) α t = (7) αq +q With the corresponding q-parametriation, the jerk parameter j(), defined in equation (), is obtained as j() = q (+) α +α (αq +q q (+) α )(α+αq +q q (+) α ) (8) with j = j( = ) = q +q (+q ). Since the expression for the jerk parameter is explicit, so we can think that we are actually parameteriing j() instead of q(). The advantage of this type of jerk parametriation is that it incorporates a wide class of viable models of cosmic evolution based on the choice of the parameter α. Similar to the generalied q-parametriation, the above jerk parametriation is also not valid at α =. As mentioned in the previous section, Zhai et al. [] parameteried j() phenomenologically aiming at measuring the departure of j from the ΛCDM value. An important difference with the work of Zhai et al. [] 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 is allowed to be fixed by the observational dataset. Obviously, the cosmological characteristics of the model given in equation () (or equation (8)) strongly depend on values of the parameters q, q and α. Using the latest H() data, in the next section, we have constrained the parameters q and q for some specific values of α. III. 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 analye them. It isknownthat thetypeiasupernovaorcmborbaodatasetispowerfulinconstrainingthecosmologicalmodels. 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() [6]. 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() dataset. From the observational point of view, the ages of the most massive and passively evolving galaxies will provide direct measurements of H() at different redshifts, which develop another type of standard probe (namely, standard clocks) in cosmology [6]. It should be noted that H() measurements are always obtained from two different techniques: galaxy differential age (also known as cosmic chronometer) and radial BAO sie methods. The Hubble parameter depending on the differential ages as a function of redshift can be written in the form of H() = d (+) dt So, H() can be obtained directly if d dt is known [65]. Also, the apparently small uncertainty of this measurement naturally increases its weight in the χ statistics. For this dataset, the χ is defined as χ H = i= [H obs ( i ) H th ( i,h,θ)] σh (, () i) (9)

4 where H obs is the observed Hubble parameter at i and H th is the corresponding theoretical value given by equation (6). Also, σ H ( i ) represents the uncertainty for the i th data point in the sample and θ denotes the model parameter. In this work, we have used the latest observational H() dataset consisting of data points in the redshift range,.7.6, larger than the redshift range that covered by the Type Ia supernova. Among them, 5 new data points of H() are obtained from the differential age method by Moresco et al. [66] and 6 data points ( data points are deduced from the radial BAO sie method and 6 data points are obtained from the galaxy differential age method) are compiled by Meng et al. [67], Cao et al. [68] and Zhang et al. [69]. In our analysis, we have also used the present value of Hubble parameter H, determined from the combined analysis with Planck+highL+WP+BAO [7]. The best fit values of the model parameters (say, θ ) from the Hubble data are estimated by minimiing the χ function in equation (). It should be noted that the confidence levels σ(68.%) and σ(95.%) are taken proportional to χ =. and 6.7 respectively, where χ = χ (θ) χ (θ ) and χ m is the minimum value of χ. An important quantity which is used for data fitting process is χ = χ m dof 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 χ m dof, then the fit is good and the observed data are consistent with the considered model. () IV. RESULTS Following the χ analysis (as described in the previous section), we have obtained the constraints on the model parameters q and q by fixing the parameter α to some constant values (,.5,.,.7,.5 and ) for the latest H() dataset. The σ and σ contours in q q plane for the proposed model is shown in figure. The best fit values with σ errors of the parameters q, q, t and j are displayed in table I. It is clear from table I that the model has almost same goodness in the viewpoint of χ for different values of α. But, there the values of q, t and j are different. The plot of the deceleration parameter q(), as given in figure, clearly shows that our model successfully generates late time cosmic acceleration along with a decelerated expansion in the past for all values of α. Figure also shows that the transition from decelerated to accelerated expansion phase took place in the redshift range within σ errors for different values of α (for details, see fifth column of table I). This range is well consistent with those previous results given in [9,, 55, 56, 7 75], though the transition redshift t slightly increases as the value of α decreases. Similarly, the evolutions of the jerk parameter j() within σ error regions are shown in figure for different values of α. It has been found from figure that the concordance ΛCDM model (dashed line) is not compatible within σ confidence level at the present epoch (i.e., = ) for α =,.7,.5 and. It has also been found that the model (for α =,.7,.5 and ) does not deviate very far from the ΛCDM model at the current epoch. These deviations of j from the ΛCDM value also need attention as the real cause behind the cosmic acceleration is still unknown. However, it is also observed from figure that for α =., the ΛCDM model is just marginally consistent within σ confidence level at =. On the other hand, the ΛCDM model is well consistent within σ confidence level at the present epoch for α =.5. Hence, our results (for α =.5 and.) also incorporate the flat ΛCDM model well within the σ error region like the work on the reconstruction of j() by Zhai et al. []. The main difference is that Zhai et al. [] imposed the jerk parameter to mimic as the flat ΛCDM model via their parametriation at =, but our work relaxes that requirement. In a word, the constraint results of j tend to favor the dynamical jerk parameter. In figure, we have shown the evolution of the normalied Hubble parameter h() = H() H for our model and have compared that with the latest points of H() dataset [66, 67]. We have also plotted data points for h() with σ error bars which have been obtained from the H() dataset using the current value of H() given by Planck observations [7]. The corresponding error in h can be estimated as [55] σ h = h σh H + σ H H ()

5 5. Α Α.5 Α q. q q q q q. Α.7.9 Α.5.7 Α q.8 q q q q q FIG. : This figure shows the σ (yellow region) and σ (red region) confidence contours in the q -q plane using the latest H() dataset and for different choices of α, as indicated in each panel. In each panel, the black dot represents the best-fit values of the pair (q,q ). where, σ H and σ H are the errors in H and H measurements respectively. We have observed from figure that the our model is well consistent with the H() data against redshift parameter for different values of α (except α = case). The reason is simple that the model with α = (i.e., q() = q +q, see equation (5)) is not reliable at high redshift. V. CONCLUSIONS In this paper, we have studied the accelerated expansion phase of the universe by using the kinematic approach. In this context, we have parameteried 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-parametried models, such as q, q ln( + ) and q + for α =, α and α = 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 χ minimiation technique using the latest points of H() dataset. Figure shows the σ and σ confidence level contours in the q -q parametric plane for different choices of α. The best fit values of the model parameters, transition redshift t and j within σ error regions for different values of α are displayed in the table I. The χ implies the same goodness of the model for all the values of α considered here. We have also shown the evolution of the normalied Hubble parameter for our model and have compared that with the latest H() dataset. In what follows, we have summaried our main conclusions in more detail: For all the values of α considered here, the kinematic model shows a smooth transition from the deceleration (q > ) phase to acceleration (q < ) phase of the universe in the recent past. It has been found that values of the transition redshift t (from decelerated to accelerated expansion) and q depend upon the choice of α. However, the changes in the values of t and q do not differ by very large values. It has also been found that the values of t obtained in this work for a wide range in the values of α, are in good agreement with the previous results as reported in [9,, 55, 56, 7 75]. We have found from figure and table I that the ΛCDM model is not well supported within

6 6 TABLE I: Best fit values of q and q with σ error bars obtained via χ minimiation method. Also, the values of t and j are given for the best fit model. α q q χ m t j χ (σ) (σ) (σ) (σ) q q... Α... Α.7 q q Α.5 Α.5 q q Α. Α FIG. : Plots of the deceleration parameter q() as a function of redshift are shown in σ error regions by considering different values of α. In each panel, the central dark line denotes the best fit curve, while the horiontal line denotes q() =. σ confidence level at the present epoch, except for the cases α =.5 and.. 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. [], where j is constrained to mimic a flat ΛCDM model at =. It has also been found that our model is well consistent with the H() data at the low redshifts for different choices of α (figure ). Therefore, we conclude the present H() 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 [76].

7 7 Α Α.5 Α. j...8 j j..5 j.5..5 Α.7 j Α.5 j 6 5 Α FIG. : Plots of the jerk parameter j() as a function of redshift are shown in σ error regions by considering different values of α. In each panel, the central dark line denotes the best fit curve, while the horiontal dashed line represents the concordance ΛCDM (j = ) model. Α Α.5 Α. h h h Α.7 Α.5 Α 6 5 h h h FIG. : Figure shows the evolution of the normalied Hubble parameter h() = H() H, as a function of, for the present model. We have also plotted h() data (red dots) with σ error bars calculated from the compilation of points H() data [66, 67]. In each panel, the central dark line denotes the best fit curve.

8 8 VI. ACKNOWLEDGMENTS AAM acknowledges the financial support from the Science and Engineering Research Board (SERB), Government of India through National Post-Doctoral Fellowship Scheme (File No: PDF/7/8). AAM also wishes to thank the Inter University Center for Astronomy and Astrophysics (IUCAA), Pune for their warm hospitality as a part of the work was done during a visit. The work of KB was supported in part by the JSPS KAKENHI Grant Number JP586 and Competitive Research Funds for Fukushima University Faculty (7RI7). [] A. G. Riess et al., Astron. J., 6, 9 (998). [] S. Perlmutter et al., Astrophys. J., 57, 565 (999). [] M. Tegmark et al., Phys. Rev. D., 69, 5 (). [] U. Seljak et al., Phys. Rev. D., 7, 55 (5). [5] D. J. Eisenstein, Astrophys. J., 6, 56 (5). [6] E. Komatsu et al., Astrophys. J. Suppl., 9, 8 (). [7] G. Hinshaw et al., Astrophys. J. Suppl., 8, 9 (). [8] P. A. R. Ade et al. [Planck Collaboration], Phys. Rev. Lett., (). [9] P. A. R. Ade et al. [Planck Collaboration], Phys. Rev. Lett., (5). [] P. A. R. Ade et al. [Planck Collaboration], Phys. Rev. Lett. 6, (6). [] P. A. R. Ade et al. [Planck Collaboration], A&A, 59, A (6). [] V. Faraoni and S. Capoiello, Fundam. Theor. Phys., 7 (). [] S. Nojiri and S. D. Odintsov, Phys. Rept., 55, 59 (). [] S. Capoiello and M.De Laurentis, Phys. Rept., 59, 67 (). [5] K. Bamba and S. D. Odintsov, Symmetry, 7, (5). [6] S. Nojiri, S. D. Odintsov and V. K. Oikonomou, Phys. Rept., 69, (7). [7] V. Sahni and A. A. Starobinsky, Int. J. Mod. Phys. D, 9, 7 (). [8] E. J. Copeland, M. Sami and S. Tsujikawa, Int. J. Mod. Phys. D, 5, 75 (6). [9] K. Bamba, S. Capoiello, S. Nojiri, S. D. Odintsov, Astrophys. Space Sci.,, 55 (). [] S. Weinberg, Rev. Mod. Phys., 6, (989). [] P. J. Steinhardt et al., Phys. Rev. Lett., 59, 5 (999). [] S. Dodelson et al., Phys. Rev. Lett., 85, 576 (). [] G. German, A. de la Macorra, Lect. Notes Phys, 66, 59 (). [] C. Shapiro, M. Turner, Astrophys. J., 69, 56 (6). [5] S. D. P. Vitenti, M. Penna-Lima, JCAP, 9, 5 (5). [6] A. Aviles, C. Gruber, O. Luongo, H. Quevedo, Phys. Rev. D, 86, 56, (). [7] S. Capoiello, R. D Agostino, O. Luongo, MNRAS, 76,, (8). [8] R. D. Blandford et al., ASP Conf. Ser., 9, 7 () [astro-ph/879]. [9] D. Rapetti, S. W. Allen, M. A. Amin and R. D. Blandford, Mon. Not. Roy. Astron. Soc., 75, 5 (7). [] M. S. Turner, A. G. Riess, Astrophys. J., 569, 8 (). [] R. Nair et al., JCAP,, 8 (). [] J. Lu, L. Xu, M. Liu, Phys. Lett. B, 699, 6 (). [] A. R. Neben, M. S. Turner, Astrophys. J., 769, (). [] O. Luongo, G. B. Pisani, A. Troisi, Int. J. Mod. Phys. D, 6, 755, (6) [5] I. Sendra, R. Lako, Mon. Not. Roy. Astron. Soc.,, 776 (). [6] S. del Campo, I. Duran, R. Herrera, D. Pavon, Phys. Rev. D, 86, 859 (). [7] L. Xu, Y. Wang, Phys. Lett. B, 7, (). [8] R.-G. Cai, Q. Su, H.-B. Zhang, JCAP,, (). [9] L. Xu, W. Li, J. Lu, JCAP, 7, (9). [] A. Mukherjee, N. Banerjee, Class. Quant. Grav.,, 56 (7). [] Z.-X. Zhai et al., Phys. Lett. B 77, 8 (). [] L. Tedesco, arxiv:8. [gr-qc] (8). [] A. G. Riess et al., Astrophys. J., 67, 665 (). [] T. Chiba, T. Nakamura, Prog. Theor. Phys.,, 77 (998). [5] M. Visser, Class. Quant. Grav.,, 6 (). [6] O. Luongo, Mod. Phys. Lett. A, 8, 58 (). [7] V. Sahni, T. D. Saini, A. A. Starobinsky, U. Alam, JETP Lett., 77, (). [8] U. Alam, V. Sahni, T. D. Saini, A. A. Starobinsky, MNRAS,, 57 (). [9] J. V. Cunha, J. A. S. Lima, MNRAS, 9, (8). [5] J. V. Cunha, Phys. Rev. D, 79, 7 (9). [5] O. Aksaru et al., EPJ Plus, 9, ().

9 [5] J. V. Cunha, J. A. S. Lima, MNRAS, 9, (8). [5] J. V. Cunha, Phys. Rev. D, 79, 7 (9). [5] S. del Campo et al., Phys. Rev. D, 86, 859 (). [55] A. A. Mamon, S. Das, Int. J. Mod. Phys. D., 5, 65 (6). [56] A. A. Mamon, S. Das, Eur. Phys. J. C, 77, 95 (7). [57] A. A. Mamon, Mod. Phys. Lett. A,, 8556 (8). [58] B. Santos, J. C. Carvalho, J. S. Alcani, Astropart. Phys., 5, 7 (). [59] Y. Gong, A. Wang, Phys. Rev. D., 7, 856 (6). [6] Y. Gong, A. Wang, Phys. Rev. D., 76, 5 (7). [6] L. Xu, H. Liu, Mod. Phys. Lett. A,, 99 (8). [6] L. Xu, J. Lu, Mod. Phys. Lett. A,, 69 (9). [6] M. Seikel, S. Yahya, R. Maartens, C. Clarkson, Phys. Rev. D 86, 8 (). [6] R. Jimene, A. Loeb, Astrophys. J., 57, 7 (). [65] J. Simon, L. Verde and R. Jimene, Phys. Rev. D., 7, (5). [66] M. Moresco et al., JCAP, 5, (6). [67] X.-L. Meng et al., arxiv:57.57 (5). [68] S.-L. Cao et al., Eur. Phys. J. C, 78, (8). [69] C. Zhang et al., A & A,, (). [7] P. A. R. Ade et al. [Planck Collaboration], A & A, 57, A6 (). [7] O. Farooq, B. Ratra, Astrophys. J., 766, L7 (). [7] S. Capoiello, O. Farooq, O. Luongo, B. Ratra, Phys. Rev. D., 9, 6 (). [7] J. Magana et al., JCAP,, 7 (). [7] A. A. Mamon, K. Bamba, S. Das, Eur. Phys. J. C, 77, 9 (7). [75] O. Farooq, F. R. Madiyar, S. Crandall, B. Ratra, ApJ, 85, 6 (7). [76] L. Yang et al., arxiv:7.79 (7). 9