New generalizations of cosmography inspired by the Padé approximant

The current accelerated expansion of the universe has been one of the most important fields in physics and astronomy since 1998. Many cosmological models have been proposed in the literature to explain this mysterious phenomenon. Since the nature and cause of the cosmic acceleration are still unknown, model-independent approaches to study the evolution of the universe are welcome. One of the powerful model-independent approaches is the so-called cosmography. It only relies on the cosmological principle, without postulating any underlying theoretical model. However, there are several shortcomings in the usual cosmography. For instance, it is plagued with the problem of divergence (or an unacceptably large error), and it fails to predict the future evolution of the universe. In the present work, we try to overcome or at least alleviate these problems, and we propose two new generalizations of cosmography inspired by the Padé approximant. One is to directly parameterize the luminosity distance based on the Padé approximant, while the other is to generalize cosmography with respect to a so-called yβ\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$y_\beta $$\end{document}-shift yβ=z/(1+βz)\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$y_\beta =z/(1+\beta z)$$\end{document}, which is also inspired by the Padé approximant. Then we confront them with the observational data with the help of the Markov chain Monte Carlo (MCMC) code emcee, and find that they work fairly well.


Introduction
From the observation of distant type Ia supernovae (SNIa) [1,2], it has been found in 1998 that the universe is experiencing an accelerated expansion. This amazing discovery was confirmed later by the observations of e.g. cosmic microwave background (CMB) [3][4][5][6], and large scale structure (LSS) [7,8]. In fact, this mysterious phenomenon has a Corresponding author: haowei@bit.edu.cn been one of the most important fields in physics and astronomy.
In the literature (see e.g. [9][10][11][12][13] for reviews), there are two representative categories of cosmological models accounting for the current cosmic acceleration. One is to introduce a new component with negative pressure, called "dark energy", in the right-hand side of the Einstein field equation in the framework of general relativity. The other is to modify the left-hand side of the Einstein field equation, namely to modify general relativity on a cosmological scale (known as "modified gravity theory"). We can constrain dark energy models and modified gravity theories by using the observational data. However, most of the observational constraints are modeldependent in fact. On the other hand, both dark energy models and modified gravity theories seem to be in agreement with the observational data; the physical mechanism to accelerate the cosmic expansion is still unclear by now [9][10][11][12][13][14]. Furthermore, it is argued in e.g. [15][16][17][18][19] that dark energy models cannot be distinguished from modified gravity theories even by using the observations of both the expansion and the growth histories. These confusions suggest that a more conservative approach to the problem of the cosmic acceleration, relying on as few model-dependent quantities as possible, is welcome. Thus, various model-independent approaches have been proposed in the literature [9][10][11][12][13][14]. A well-known one is the parameterization of equation-of-state parameter (EoS), such as w = w 0 + w 1 z [20,21], and w = w 0 + w a z/(1 + z) [22,23], where z is the redshift. Another powerful model-independent approach is cosmography [14,[24][25][26][27][28][29][30][31][32][33][34][35][36][76][77][78][79][80][81][82][83][84][85][86][87]. To the best of our knowledge, it was first discussed by Weinberg [27,28] and extended by Visser [29,30] recently. Using cosmography, one can analyze the evolution of the universe without assuming any underlying theoretical model. The only necessary assumption of cosmography is the cosmological principle, so that the spacetime metric is the one of the Friedmann-Robertson-Walker (FRW) universe, ds 2 = −c 2 dt 2 +a 2 (t) dr 2 1 − kr 2 + r 2 dθ 2 + sin 2 θ dφ 2 , (1) in terms of the comoving coordinates (t, r, θ, φ), where c is the speed of light, a is the scale factor, and k >, =, < 0 corresponds to a spatially close, flat, open universe, respectively. Introducing the so-called cosmographic parameters, namely the Hubble constant H 0 , the deceleration q 0 , the jerk j 0 , and the snap s 0 (defined below), one can expand the scale factor a in terms of a Taylor series with respect to cosmic time t [27][28][29][30], and also the luminosity distance d L with respect to redshift z [14,24,25,[27][28][29][30][31][32], So, one can study the universe in a model-independent way by using cosmography. It is easy to see that the key of cosmography is to expand the quantities under consideration as a Taylor series with respect to redshift z. However, it is well known that the Taylor series converges only for small z around 0, and it might diverge at high redshift (especially when z > 1). A possible remedy is to replace the redshift z with the so-called y-shift, y ≡ z/(1 + z) [26,[33][34][35]80,81]. In this case, y < 1 holds in the whole cosmic history 0 ≤ z < ∞, and hence the Taylor series with respect to y converges. However, there still exist several serious problems in the case of y = z/(1 + z). The first is that the error of a Taylor approximation throwing away the higher order terms will become unacceptably large when y is close to 1 (say, when z > 9). The second is that the cosmography in terms of y = z/(1 + z) cannot work well in the cosmic future −1 < z < 0. The Taylor series with respect to y = z/(1 + z) does not converge when y < −1 (namely z < −1/2), and it drastically diverges when z → −1 (it is easy to see that y → −∞ in this case). So, this y-shift cosmography fails to predict the future evolution of the universe. Note that there are other y-shifts considered in the literature, for instance, y 1 ≡ arctan (z/(1 + z)), y 2 ≡ z/(1 + z 2 ), and y 3 ≡ arctan z [36]. However, they are purely written by hand, without solid motivation. On the other hand, |y| > 1 at suit-able redshift z, since the function arctan x ∈ (−π/2, +π/2), and hence the Taylor series does not converge.
In the present work, we try to overcome the problems of cosmography mentioned above. We are mainly interested in the cosmography of the luminosity distance d L , since it can be confronted with the observational data directly. In fact, the new generalizations considered in the present work are inspired by the so-called Padé approximant. In Sect. 2, we briefly review the key points of the Padé approximant, and then we parameterize the luminosity distance d L based on the Padé approximant. We confront this parameterization of the luminosity distance d L with the observational data and see whether it works well. Note that in this work we use the Markov chain Monte Carlo (MCMC) code emcee [37,38] in the data fitting. In Sect. 3, inspired by the Padé approximant, we propose a new y β -shift and then derive the cosmography of the luminosity distance d L by expanding it as a Taylor series with respect to the new y β -shift. This cosmography is completely free from the problems mentioned above. We also confront it with the observational data. Finally, some brief concluding remarks are given in Sect. 4.

Padé parameterization of the luminosity distance
The so-called Padé approximant can be regarded as a generalization of the Taylor series. For any function f (x), its Padé approximant of order (m, n) is given by the rational function [39][40][41][42][43][44][45][46][47][48][49] f (x) = α 0 + α 1 x + · · · + α m x m 1 + β 1 x + · · · + β n x n , where m and n are both non-negative integers, and α i , β i are all constants. Obviously, it reduces to the Taylor series when all β i = 0. Actually in mathematics, a Padé approximant is the best approximation of a function by a rational function of given order [42]. In fact, the Padé approximant often gives a better approximation of the function than truncating its Taylor series, and it may still work where the Taylor series does not converge [42]. So, considering the Padé approximant in cosmology is well motivated.
Here, we directly parameterize the luminosity distance d L based on the Padé approximant, Note that the speed of light c and Hubble constant H 0 are introduced from dimensional point of view. How to choose the order (m, n) of the Padé approximant is important. If the order is too low, the error of the Padé approximant deviating from the real luminosity distance d L will be unacceptably Table 1 The best-fit model parameters with 1σ , 2σ , 3σ uncertainties. The corresponding χ 2 min and χ 2 min /dof are also given. These results are obtained by fitting the Padé parameterization (6)  large. If the order is too high, the number of free coefficients are too many and the uncertainties will be large. So, we choose a moderate order (2,2) in this work, and then Eq. (5) becomes Obviously, it can work well in the whole redshift range −1 < z < ∞, including not only the past but also the future of the universe. Especially, it is still finite even when z 1. Also, it is easy to ensure the denominator not equal to zero at any redshift z for suitable β 1 and β 2 . Thus, this parameterization based on the Padé approximant can easily avoid the problems of the usual cosmography mentioned above. It is worth noting that this parameterization is a generalization of cosmography, since the Padé approximant is a generalization of the Taylor series in fact.
Naturally, it is important to confront the Padé parameterization (6) with the observational data, and see whether this parameterization works well. Since SNIa data are directly related to the luminosity distance, we can use them to constrain the parameterization (6). Here, we consider the Union2.1 SNIa dataset [50,51] consisting of 580 data points, which are given in terms of the distance modulus μ obs (z i ). On the other hand, the theoretical distance modulus is given by [1,2,27,28,[52][53][54][55][56][57][58][59][60][61] where μ 0 ≡ 42.38-5 log 10 h, and h is the Hubble constant H 0 in units of 100 km/s/Mpc. In our case, D L has been given in Eq. (6). Correspondingly, the χ 2 from 580 Union2.1 SNIa is given by where σ is the corresponding 1σ error. The best-fit model parameters are determined by minimizing χ 2 . In this work, we use the Markov chain Monte Carlo (MCMC) code emcee [37,38] to find the best fits and the corresponding 68.3 %, 95.4 %, and 99.7 % confidence levels. We present the best-fit parameters with 1σ , 2σ , 3σ uncertainties and the corresponding χ 2 min in Table 1. The 1D marginalized distribution, and 1σ , 2σ , 3σ contours in the 2D model parameter spaces are also given in Fig. 1.
In addition to SNIa, the observation of cosmic microwave background (CMB) anisotropy [62,63] is another useful probe. However, using the full data of CMB to perform a global fitting consumes a large amount of computation time and power. As an alternative, one can instead use the shift parameter R [64,65] from CMB, which has been used extensively in the literature (including the works of the Planck and the WMAP Collaborations). It is argued in e.g. [66][67][68] that the shift parameter R is model-independent and contains the main information of the observation of CMB. As is well known, the shift parameter R is defined by [62][63][64][65][66][67][68][69] where d A (z) is the angular diameter distance to the redshift z, which can be related to the luminosity distance d L through (see e.g. the textbooks in Refs. [27,28]) and, in our case, D L is given in Eq. (6). The redshift of the recombination z * = 1089.90, which was determined by the latest Planck 2015 data [62]. m0 is the present frac- Fig. 1 The 1D marginalized distribution, and 1σ , 2σ , 3σ contours in the 2D model parameter spaces. These results are obtained by fitting the Padé parameterization (6) to SN data. See the text for details tional density of pressureless matter, and it was determined as m0 = 0.308 by the latest Planck 2015 data [62]. On the other hand, the observational value of R has also been determined to be R obs = 1.7382 ± 0.0088 by the latest Planck 2015 data [63]. So, the χ 2 from CMB is given by and then the total χ 2 from the combined SN + CMB data reads where χ 2 SN is given in Eq. (8). Again, we use the MCMC code emcee [37,38] to find the best fits and the corresponding 68.3 %, 95.4 %, and 99.7 % confidence levels to the combined SN + CMB data. We present the best-fit parameters with 1σ , 2σ , 3σ uncertainties and the corresponding χ 2 min in the last column of Table 1. The 1D marginalized distribution, and 1σ , 2σ , 3σ contours in the 2D model parameter spaces are also given in Fig. 2. Thanks to CMB data, it is easy to see that the constraints on all parameters are significantly tightened (nb. Table 1). By confronting the Padé parameterization (6) with the observational data, we see that this generalized cosmography works well.

y β -shift cosmography
In this section, we propose another generalization of cosmography inspired by the Padé approximant. We stress that it is completely independent of the one proposed in the previous section. At first, we propose a new y β -shift inspired by the Padé approximant, and then derive the cosmography of the luminosity distance d L by expanding it as a Taylor series with respect to this new y β -shift. We will also confront it with the observational data and see whether it works well.

The formalism of y β -shift cosmography
At first, we give the motivation to propose the new y β -shift. As is well known, the standard cosmography of the luminosity distance d L with respect to the redshift z is given by [14,24,25,[27][28][29][30][31][32] One can instead parameterize it based on the Padé approximant, Fig. 2 The same as in Fig. 1, except for SN + CMB data Requiring that Eq. (14) coincides with Eq. (13) at z → 0, we find that α 0 → 1. So, we consider the Padé approximant up to order (1, 1), When z 1, Eq. (15) can be expanded as Comparing it with Eq. (13), Putting Eqs. (13) and (17) together, it is quite interesting to see that the quantity z/(1 + βz) plays a role similar to the redshift z. This inspires us to propose a so-called y β -shift, where β is a dimensionless constant. Obviously, y β = z if β = 0, and y β = z/(1 + z) if β = 1. Noting that the redshift z [14,24,25,[27][28][29][30][31][32] and y-shift y = z/(1 + z) [26,[33][34][35] are extensively used in the usual cosmography, our y β can be regarded as their natural generalization. So, the cosmography with respect to y β -shift is also a natural generalization of the usual cosmography in fact. We stress that the above discussions are only arguments to justify y β , rather than strict derivations. Now, let us see how it might overcome the problems of the usual cosmography mentioned in Sect. 1. First, y β is inspired by the Padé approximant, and hence it is well motivated, not purely written by hand. Second, if 1 − β 1, we have y β < 1 even for z 1. So, the y β -shift cosmography can converge safely. Third, it can describe the future evolution of the universe, if 1 + βz = 0 in the redshift range −1 < z < 0. In fact, to ensure y β remains regular in not only the past but also in the future of the universe (−1 < z < ∞), it is required that In the following, let us derive the cosmography of the luminosity distance d L with respect to y β -shift. As in e.g. [14,[24][25][26][27][28][29][30][31][32][33][34][35], it is convenient to introduce the following functions: which are usually referred to as the Hubble, deceleration, jerk, snap, and lerk parameters. Using these definitions, the Taylor series expansion of scale factor a up to 5th order around the time t 0 reads where the subscript "0" indicates the value of the corresponding quantity evaluated at the time t 0 . The physical distance traveled by a photon that is emitted at the time t and absorbed at the current epoch t 0 is given by where the time difference δt = t 0 − t is called the "lookback time". So, we have Using Eq. (25), its right-hand side can be expanded as a Taylor series with respect to H 0 D/c , Therefore, from Eqs. (27) and (28), we obtain where Z (1) We can expand D as a Taylor series up to 5th order with respect to y β , After some algebra, Eq. (36) becomes in which we have used Z The role of the observable physical quantity is played by the luminosity distance. Let the photon be emitted at rcoordinate r = 0 at the time t, and absorbed at r -coordinate r = r 0 at the time t 0 . Then the luminosity distance d L is given by (see e.g. the textbooks in Refs. [27,28]) To calculate d L (D), we need r 0 (D). From Eq. (1), it is easy to obtain [27][28][29][30] Using Eq. (25), we can calculate the integration [14] t 0 where a 0 = a(t 0 ). At first glance, we should deal with three cases with space curvature k = +1, 0, −1 separately. Fortunately, we need not to do so. As is well known, the Taylor series expansions of sin x and sinh x are given by sin x = x − x 3 /3! + x 5 /5! − x 7 /7! + · · · and sinh x = x + x 3 /3! + x 5 /5! + x 7 /7! + · · · , respectively. Noting that r 0 (D) → sin( * ) for k = +1, and r 0 (D) → sinh( * ) for k = −1, we can write the Taylor series expansion of r 0 (D) in Eq. (44) as a uniform expression for k = +1, 0, −1, i.e. [14,29,30] Substituting Eq. (45) into Eq. (46), we have Substituting Eqs. (28), (47), (35) with (38)-(42) into Eq. (43), we finally obtain the cosmography of the luminosity distance d L with respect to y β -shift, where D (1) One can check that if β = 0 or β = 1, these results reduce to the one of cosmography with respect to redshift z (e.g. [14,29,30]) or y-shift y = z/(1 + z) (e.g. [14,26,[33][34][35]80,81]), respectively. The cosmography with respect to y β -shift y β = z/(1 + βz) obtained here is more general.

Cosmological constraints from the observational data
It is natural to confront our new cosmography with the observational data, and see whether this new cosmography works well. Here, we only consider a flat FRW universe with Similar to Sect. 2, we first consider the constraints from the Union2.1 SNIa dataset [50,51]. The χ 2 from 580 Union2.1 SNIa is given in Eq. (8), in which μ th is given by Eq. (7), D L ≡ H 0 d L /c, and d L is given by Eq. (48) in our case. Again, we use the MCMC code emcee [37,38] to find the best fits and the corresponding 68.3 %, 95.4 % and 99.7 % confidence levels to SN data. Note that the prior 0 < β < 1 is required in Eq. (19). We present the best-fit parameters with 1σ , 2σ , 3σ uncertainties and the corresponding χ 2 min in Table 2. The 1D marginalized distribution, and 1σ , 2σ , 3σ contours in the 2D model parameter spaces are also given in Fig. 3. It is easy to see that the parameter β cannot be well constrained, and its 1D marginalized distribution is not Gaussian. So, we temporarily relax the prior to 0 < β < 2. In this case, the best-fit parameters with 1σ , 2σ , 3σ uncertainties and the corresponding χ 2 min are given in the last column of Table 2, while the 1D marginalized distribution, and 1σ , 2σ , 3σ contours in the 2D model parameter spaces are also Table 2 The best-fit model parameters with 1σ , 2σ , 3σ uncertainties. The corresponding χ 2 min and χ 2 min /dof are also given. These results are obtained by fitting the y β -shift cosmography (48) to SN data with the priors 0 < β < 1 and 0 < β < 2, respectively. See the text for details  . 3 The 1D marginalized distribution, and 1σ , 2σ , 3σ contours in the 2D model parameter spaces. These results are obtained by fitting the y β -shift cosmography (48) to SN data with the prior 0 < β < 1. See the text for details given in Fig. 4. Now, the constraint on β looks well. The following discussions hold for both cases with 0 < β < 1 and 0 < β < 2 in fact. It is easy to see that SN data favor a best-fit β close to 1, and hence y β is close to the usual y-shift y = z/(1+z). However, β can still significantly deviate from 1 in the wide 1σ , 2σ , 3σ regions. In particular, β can even be close to 0 in the 3σ region. On the other hand, we find that q 0 < 0 at 1σ confidence level, which indicates the cosmic expansion is accelerating [nb. the definition in Eq. (21)]. Of course, it is not surprising that the constraints on l 0 and s 0 are fairly loose, since we use only SN data here. Let us further consider the constraints from the combined SN + CMB data. The corresponding χ 2 is given by Eq. (12), in which χ 2 SN is given by Eq. (8) and χ 2 CMB = (R − R obs ) 2 /σ 2 R , while D L ≡ H 0 d L /c and d L is given by Eq. (48) in our case. Note that the prior 0 < β < 1 is still required in Eq. (19). We present the best-fit parameters with 1σ , 2σ , 3σ uncertainties and the corresponding χ 2 min in Table 3. The 1D marginalized distribution, and 1σ , 2σ , 3σ contours in the 2D model parameter spaces are also given in Fig. 5. Obviously, the constraints on all parameters are significantly tightened, mainly thanks to the CMB data. It is interesting to see that the combined SN + CMB data favor a best-fit β close to 0, and hence y β is close to the usual redshift z. Noting that SN data favors a best-fit β close to 1 as mentioned above, this indicates that there is tension between the SN and the CMB data. However, we stress that the constraints from SN and SN + CMB data are still consistent within the 3σ confidence level. On the other hand, we see that q 0 < 0, j 0 > 0 beyond 3σ confidence level, which strongly indicates the cosmic expansion is accelerating, and the acceleration is increasing [nb. the definitions in Eqs. (21) and (22)].

Concluding remarks
The current accelerated expansion of the universe has been one of the most important fields in physics and astronomy since 1998. Many cosmological models have been proposed in the literature to explain this mysterious phenomenon. Since the nature and cause of the cosmic acceleration are still unknown, model-independent approaches to study the evolution of the universe are welcome. One of the powerful model-independent approaches is the so-called cosmography. It only relies on the cosmological principle, without postulating any underlying theoretical model. However, there are several shortcomings in the usual cosmography. In the present work, we try to overcome or at least alleviate these Fig. 4 The same as in Fig. 3, except for the prior 0 < β < 2 Table 3 The same as in Table 2 problems, and propose two new generalizations of cosmography inspired by the Padé approximant. We also confront them with the observational data with the help of the Markov chain Monte Carlo (MCMC) code emcee [37,38], and find that they work fairly well. In the literature (e.g. [14,[24][25][26][27][28][29][30][31][32][33][34][35][36][76][77][78][79][80][81][82][83][84][85][86][87][88]), there exist a number of works on cosmography. They focused on various issues in cosmology and made much significant progress. However, the previous works used ordinary cosmography mainly with respect to the redshift z or the so-called y-shift y = z/(1 + z). We stress that they are all plagued with the problem of divergence (or the unacceptably large error when y ∼ 1), and all fail to predict the future evolution of the universe (especially when z ∼ −1 or y < −1), as mentioned in Sect. 1. In the present work, the two new generalizations of cosmography proposed in Sects. 2 and 3 can instead avoid or at least alleviate these problems of ordinary cosmography. In addition, our new generalizations of cosmography are well motivated by the Padé approximant, rather than written purely by hand. In fact, it has been solidly proved by mathematicians that any (even unknown) function can be well approximated by a Padé approximant [39][40][41][42][43][44], Fig. 5 The same as in Fig. 3, except for SN + CMB data, and the prior 0 < β < 1 is required and the Padé approximant often gives a better approximation of the function than truncating its Taylor series, and it may still work where the Taylor series does not converge [42]. Noting that ordinary cosmography is actually based on a Taylor series, the advantages of the Padé approximant mentioned above make our new generalizations better than the ordinary cosmography used in the literature. Although the Padé approximant has been considered in cosmology (see e.g. [45][46][47][48][49][70][71][72][73][74][75]), here we instead use it in a fairly different issue. For example, in e.g. [48,49] the Padé approximant was used to study the issues of an analytical approximation of the luminosity distance in the XCDM model, EoS parameterizations, and gamma-ray burst cosmology, while we use it to generalize cosmography in this work. By all the above arguments, we stress that the present work is significantly different from the previous works on both cosmography and the Padé approximant in the literature.
The key to avoid or at least alleviate the problems of the ordinary cosmography is the denominator of the Padé approximant. Considering Eqs. (5) or (4), if the order n of the denominator is larger than or equal to the order m of the numerator, this Padé approximant with n ≥ m will not diverge even for z 1. In fact, this is just the case of our two generalizations [nb. Eqs. (6) and (18)]. On the other hand, for suitable parameters β i , it is easy to ensure the denominator not to equal zero for the very wide redshift range −1 < z < ∞, and hence the Padé approximant can avoid divergence. The shortcomings of ordinary cosmography mainly have their roots in the Taylor series. So, generalizing a Taylor series to the Padé approximant brings about a possible way out. If all β i = 0 in the denominator of the Padé approximant, it reduces to the usual Taylor series. However, a denominator not equal to 1 and 0 makes a big difference.
Here, we would like to clarify the main difference between the two new generalizations proposed in Sects. 2 and 3. Noting that ordinary cosmography is based on a Taylor series, the first one generalizes cosmography by directly generalizing the Taylor series to the Padé approximant when we expand the luminosity distance d L [nb. Eq. (5)]. On the other hand, the second one instead generalizes cosmography by generalizing the redshift z or the y-shift y = z/(1 + z) to the so-called y β -shift y β = z/(1 + βz) (which is inspired by the Padé approximant), but we still expand the luminosity distance d L in a Taylor series [nb. Eq. (48)], rather than the Padé approximant itself. This is the main difference. Noting that y β = z and z/(1 + z) when β = 0 and 1, respectively, we call y β = z/(1 + βz) the "y β -shift" in analogy to the well-known terminology of the "y-shift" y = z/(1 + z) used in the literature, while the terminology "y-shift" comes from the terminology "red-shift" in fact.
It is of interest to quantitatively compare our two new generalizations of cosmography, and also compare them with the ordinary cosmography (we thank the referee for pointing out this issue). Let us compare our y β -shift cosmography with the ordinary cosmography first. The observational constraints on our y β -shift cosmography from the SN + CMB data are presented in Table 3 and Fig. 5. It is easy to see that β = 0 and β = 1 deviate from the best fit far beyond the 3σ confidence level; see especially the leftmost column of Fig. 5. Note that our y β -shift reduces to the ordinary redshift z and the y-shift y = z/(1 + z) when β = 0 and 1, respectively. Therefore, our y β -shift cosmography can fit the SN + CMB data significantly better than ordinary cosmography. On the other hand, from the last column of Table 1 and Table 3, χ 2 min = 562.171, and χ 2 min /dof = 0.978 by fitting the Padé parameterization to the SN + CMB data, while χ 2 min = 562.654 and χ 2 min /dof = 0.979 by fitting the y βshift cosmography to the SN + CMB data. Therefore, our first generalization of cosmography is slightly better than the second one.
Other remarks are in order. First, the Padé approximant has been used previously in cosmology, e.g. in slow-roll inflation, the reconstruction of the scalar field potential, data fitting and analytical approximation of the luminosity distance, EoS parameterizations, gamma-ray burst cosmology, and cosmological perturbations in LSS (see e.g. [45][46][47][48][49][70][71][72][73][74][75]). We advocate further uses of the Padé approximant in cosmology. Second, in the present work we only derive the generalized cosmography of the luminosity distance d L . In fact, one can easily obtain the corresponding cosmography of other observable quantities, such as the angular diameter distance d A , the photon flux distance d F , the photon count distance d P , and the deceleration distance d Q , since they can be readily related to the luminosity distance d L (see e.g. [26]). Third, we confront the generalized cosmography with the observational data only for the spatially flat FRW universe (k = 0). In fact, one can do this for the k = 0 cases easily, and hence we do not present them here. Finally, the usual cosmography with respect to the redshift z and y-shift y = z/(1 + z) has been extensively applied to various cosmological issues, for instance, the EoS of dark energy, modified gravity theories like f (R) and f (T ) theories, gamma-ray burst cosmology, and so on (see e.g. [14,[24][25][26][27][28][29][30][31][32][33][34][35][36][76][77][78][79][80][81][82][83][84][85][86][87][88]). The new generalizations of cosmography proposed in the present work can also be used in these cosmological issues, and we leave it to the future works.