Next Article in Journal
Polymer Dynamics of Isotropic Universe in Ashtekar and in Volume Variables
Next Article in Special Issue
Dynamical Signature: Complex Manifolds, Gauge Fields and Non-Flat Tangent Space
Previous Article in Journal
Prospects for Probing Axionlike Particles at a Future Hadron Collider through Top Quark Production
Previous Article in Special Issue
Is the Hubble Crisis Connected with the Extinction of Dinosaurs?
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A New Analytic Approximation of Luminosity Distance in Cosmology Using the Parker–Sochacki Method

Department of Mathematics, Faculty of Science, University of Malta, MSD2080 Msida, Malta
Universe 2022, 8(6), 300; https://doi.org/10.3390/universe8060300
Submission received: 10 March 2022 / Revised: 23 May 2022 / Accepted: 23 May 2022 / Published: 26 May 2022
(This article belongs to the Special Issue New Directions in Gravitational Physics and Cosmology)

Abstract

:
The luminosity distance d L is possibly the most important distance scale in cosmology and therefore accurate and efficient methods for its computation is paramount in modern precision cosmology. Yet in most cosmological models the luminosity distance cannot be expressed by a simple analytic function in terms of the redshift z and the cosmological parameters, and is instead represented in terms of an integral. Although one can revert to numerical integration techniques utilizing quadrature algorithms to evaluate such an integral, the high accuracy required in modern cosmology makes this a computationally demanding process. In this paper, we use the Parker–Sochacki method (PSM) to generate a series approximate solution for the luminosity distance in spatially flat Λ CDM cosmology by solving a polynomial system of nonlinear differential equations. When compared with other techniques proposed recently, which are mainly based on the Padé approximant, the expression for the luminosity distance obtained via the PSM leads to a significant improvement in the accuracy in the redshift range 0 z 2.5 . Moreover, we show that this technique can be easily applied to other more complicated cosmological models, and its multistage approach can be used to generate analytic approximations that are valid on a wider redshift range.

1. Introduction

The luminosity distance in the spatially flat Λ CDM cosmology is expressed exactly in terms of an elliptic integral. Evaluation of this integral via numerical quadratures tend to be computationally heavy when high accuracy is required in applications that involve ultra precise measurements of cosmological parameters such as the calculation of H 0 in view of the problem of the Hubble constant tension [1,2]. Different strategies for obtaining analytical approximations for this integral have been proposed during the last two decades, and these can be divided into three types. The first type involves obtaining an algebraic fitting formula expressed in terms of elementary functions for the integral which compares quite well with the numerical solution over a wide redshift range and for different values of the matter density parameter Ω m . This started with the work of Ue Li Pen [3] who obtained a formula which has a relative error of less than 0.4 % for 0.2 < Ω m < 1 . A similar approach was used by Wickramasinghe and Ukwatta [4] and Liu et al. [5] who obtained alternative fitting formulae that run faster and have smaller relative errors than that obtained by Pen [3]. Finally, the formula obtained by Adachi Kasai [6] which is based on the Padé approximation has even smaller relative error for the range of density parameter 0.3 Ω m 1 and redshift range 0.03 z 1000 . The second strategy is to identify exact analytic expressions for the luminosity distance which are often expressed in terms of transcendental functions, such as elliptic integrals [5,7,8,9] or hypergeometric functions [10]. Being exact expressions, this approach leads to a better accuracy, but the presence of such transcendental functions makes these more computationally demanding to evaluate. The third strategy involves using truncated Taylor series approximations [11,12,13] for the luminosity distance in a relatively small redshift range. This technique is easy to use but has divergence problems for higher values of z, and the fact that nowadays supernova data goes as far as z 2.35 [14] this can pose a significant problem. Thus several recent works suggest the use of the Padé approximant [15,16,17] to represent the luminosity distance, due to its better convergence over a larger redshift range. Instead of obtaining power series or Padé approximations to the luminosity integral directly, some authors have also considered the associated differential equation with corresponding initial conditions and obtained approximate analytical solutions for the luminosity distance using various methods, such as the homotopy perturbation method (HPM) [18] and its variants [19], the variational iteration method (VIM) [20] and the Daftardar-Jafari Method (DJM) [21]. Unlike the fitting formulae developed earlier, which are valid over a wide redshift range, the series approximations obtained using this strategy converge on a much limited redshift range and the relative error for the approximated luminosity can be larger for higher values of z meaning that one has to reach a compromise between accuracy and redshift convergence interval, although the accuracy and convergence can be somewhat improved by using a combination of the Padé approximant with the HPM [19]. However, these methods have the advantage that they can be easily applied to any cosmological model, unlike the earlier fitting formulae or exact expressions which only hold for Λ CDM cosmology. This is particularly important when comparing different dark energy cosmological models with observational data.
In this paper, we follow the latter strategy and apply the Parker–Sochacki method (PSM) [22,23] to obtain an approximate solution of the ordinary differential equation satisfied by the luminosity distance function. PSM is an extension of the Picard iteration [24] which is in turn an algorithm for solving linear systems of differential equations. It is a simple yet powerful technique for computing iteratively the coefficients of the Maclaurin’s series solution of a polynomial system derived from a nonlinear initial value problem (IVP) or a boundary value problem (BVP). Over the last two decades, interest in PSM has increased dramatically and it has been applied to various problems of technological and scientific importance (see for examples Refs. [25,26,27,28,29,30,31,32,33,34,35,36,37]) including stiff differential equations whose solution cannot be obtained efficiently with existing numerical solvers [38]. This would be the first example of the application of PSM in cosmology. The main advantage of the PSM over the other methods mentioned above, is that the radius of convergence R of the Maclaurin’s series solution to the polynomial system can be obtained easily and a priori in terms of the parameters in the system. This would allow the domain of the problem to be divided into subdomains over which the PSM can be used individually to obtain convergent solutions, such that the end result is a piecewise convergent solution. This allows one to extend the domain of convergence of the approximate analytic solution without compromising the accuracy of the solution. It has been shown [39] that such an approach compares relatively well with the Runge–Kutta method of order four.
This paper is structured as follows. In the next section, we introduce the Parker–Sochacki method. Then in Section 3 we obtain the differential equation satisfied by the luminosity distance function in Λ CDM cosmology and generate an approximate solution for d L ( z ) using the PSM. This will be followed by a comparison of this solution with the corresponding numerical solution and with the approximate solutions obtained by using the other methods mentioned above. The PSM is also used to obtain the luminosity distance in a dynamical dark energy model, which is chosen to be the Chevallier, Polarski and Linder (CPL) model [40,41,42] with the equation of state (EOS) function ω ( z ) = ω 0 + ω 1 z / ( 1 + z ) , where ω 0 and ω 1 are constants. This is followed with a discussion of the results and a conclusion. In this paper, we will use geometrized units such that 8 π g = c = 1 .

2. The Parker–Sochacki Method

The Parker–Sochacki method is a modification of the conventional power series method for solving non-linear ordinary differential equations which may include transcendental functions. This involves expressing the given differential equation into a polynomial system of differential equations (also called polynomial projection) by introducing appropriate new variables to represent the non-polynomial functions in the original differential equation. The system is then solved using a standard power series method. The concept of polynomial projection may be better explained by referring to a simple example (taken from Ref. [43]). So consider the following first order nonlinear IVP
y ( t ) = sin y ( t ) e t 2 t , y ( 1 ) = 2 .
This can be expressed in a polynomial system of differential equations by introducing the variables
x 1 = y , x 2 = sin ( y e t 2 ) , x 3 = cos ( y e t 2 ) , x 4 = e t 2 , x 5 = t , x 6 = t 1 / 2 ,
such that (1) is recast in the form
x 1 = x 2 x 6 , x 2 = 2 x 1 x 3 x 4 x 5 + x 2 x 3 x 4 x 6 , x 3 = 2 x 1 x 2 x 4 x 5 x 2 2 x 4 x 6 , x 4 = 2 x 4 x 5 , x 5 = 1 , x 6 = 1 2 x 6 3 ,
with the initial conditions x 1 ( 1 ) = 2 , x 2 ( 1 ) = sin ( 2 e 1 ) , x 3 ( 1 ) = cos ( 2 e 1 ) , x 4 ( 1 ) = e 1 , x 5 ( 1 ) = 1 , x 6 ( 1 ) = 1 . So the solution y ( t ) to the given differential equation in (1) is now embedded in the polynomial system (3). This is called a projection for the solution y ( t ) . Now by expressing each of the six variables x i , i = 1 6 in terms of a truncated Maclaurin series in t and substituting in (3), one can apply the simple Picard iteration to solve the polynomial system and so get an approximate series solution for y ( t ) . In general, given a nonlinear IVP of the form
y = f ( t , y ) , y ( t 0 ) = y 0 , t 0 t t f ,
where t 0 R and y 0 R n , one can transform this to a projectively polynomial system by the introduction of auxiliary variables as seen in the above example, to get
x = g ( x ) , x ( t 0 ) = x 0 , t 0 t t f ,
where x 0 R m and g : R m R m . Taking t 0 = 0 without loss of generality and using the truncated Maclaurin’s expansions x ( t ) = i = 0 p x i t i and g ( x ) = i = 0 p g i t i and substituting in the above polynomial system, gives
x i + 1 = g i i + 1 ,
where the coefficients g i are obtained by the repeated use of the simple Cauchy product. By using (6), one can then obtain the truncated series solution
x ( t ) = i = 0 p x i t i ,
of the polynomial system in (18). We now introduce a few notations. The definition of the norm of the vector x is | | x | | = max 1 i m | x i | , where x R m . If g is a vector polynomial function and X α = x 1 α 1 x 2 α 2 x m α m , then the components of g are given by
g i ( x 1 , , x m ) = | α | k a α , i X α ,
where | α | = α 1 + + α m , and k = deg ( g i ) is the degree of g i . We also define
Σ g i = | α | k | a α , i | ,
and Σ g = max ( Σ g 1 , , Σ g m ) , deg ( g ) = max ( deg ( g 1 ) , , deg ( g m ) ) . Then the convergence for the above series solution in (7) obtained by the PSM is given by the following Taylor approximation theorem [44].
Theorem 1.
If x satisfies (5), k = deg ( g ) 2 , α = max ( 1 , | | x 0 | | ) , M = ( k 1 ) Σ g α k 1 , t 0 = 0 , and | t | < 1 / M then
| | x ( t ) i = 0 p x i t i | | α | M t | p + 1 1 | M t | .
In most cases the radius of convergence | t | < 1 / M of the series solution as stated in the above theorem is less than the interval of integration t 0 t t f . Therefore as proposed in [44] one can apply the PSM on smaller subintervals such that the convergence of the solution is preserved on each subinterval. In other words, the interval t 0 t t f is divided into subintervals by introducing the nodes t i = t 0 + i h where h = 1 / M and i = 0 , 1 , , N = t f t 0 h , such that the union of the non-overlapping sub intervals [ t k 1 , t k ] gives [ t k 1 , t k ] = [ t 0 , t f ] . Then, when applying the PSM on each subinterval the value of the function at the end-point of a given subinterval becomes the initial condition for the next subinterval. Combining the solutions obtained by the PSM for each subinterval yields a highly accurate piecewise smooth solution for the entire integration domain t 0 t t f . Moreover, by increasing the order of the approximation p in (7) one can also increase the accuracy of the approximation in each subinterval. This multistage approach for the implementation of the PSM has been used to obtain accurate series approximations of various non-linear and stiff ordinary differential equations (see for example Ref. [38]).
In the next section, we obtain the non-linear differential equation satisfied by the luminosity function and solve this for two different cosmological models by using the PSM via the multistage approach as described above.

3. Differential Equation for the Luminosity Distance in a Flat Universe

According to recent observations our universe is spatially flat [45,46] and is homogeneous and isotropic on the large scale, such that it can be described by the FLRW metric
d s 2 = d t 2 + a ( t ) 2 d r 2 + r 2 ( d θ 2 + sin 2 θ d ϕ 2 ) ,
where a ( t ) is the scale factor. One of the most fundamental distance scales in cosmology is the luminosity distance which is defined by d L = L / ( 4 π F ) , where L is the absolute luminosity of the source and F is the associated flux measured by the observer. In a spatially flat universe this is given by [47]
d L = c ( 1 + z ) H 0 0 z d z W ( z ) ,
where W ( z ) = ( H ( z ) / H 0 ) 2 depends on the Hubble parameter H ( z ) as a function of the redshift z and H 0 is the current value of the Hubble constant such that W ( 0 ) = 1 . Even in the simple case of the Λ CDM universe the above integral cannot be integrated explicitly and normally one has to revert to numerical integration or use any of the approximations mentioned in the Introduction. Introducing the function
u ( z ) = H 0 d L c ( 1 + z ) ,
the above integral can be written as
u ( z ) = 0 z d z W ( z ) .
The function u ( z ) then satisfies the non-linear differential equation
u ( z ) + 1 2 W ( z ) ( u ( z ) ) 3 = 0 ; u ( 0 ) = 0 , u ( 0 ) = 1 ,
where denotes differentiation with respect to z. In the next section, we obtain analytic approximations to this differential equation by using the multistage PSM described in the previous section for two cosmological models, namely the Λ CDM and the Chevallier–Polarski–Linder (CPL) models.

4. Luminosity Distance in Cosmology

4.1. Λ CDM Cosmological Model

For the Λ CDM model the function W ( z ) in the non-linear differential Equation (15) is given in terms of the cosmological parameter Ω m by
W ( z ) = Ω m ( 1 + z ) 3 + ( 1 Ω m ) .
By introducing the additional auxiliary variables
V ( z ) = u ( z ) , S ( z ) = W ( z ) = 3 Ω m ( 1 + z ) 2 , Q ( z ) = ( 1 + z ) ,
the differential equation in (15) can be incorporated in the following polynomial system which will be solved using the PSM as explained earlier
u = V V = 1 2 S V 3 S = 6 Ω m Q Q = 1 ,
where u ( 0 ) = 0 , V ( 0 ) = 1 , S ( 0 ) = 3 Ω m , Q ( 0 ) = 1 . Note that the required function u ( z ) which is related to the luminosity distance function d L ( z ) by (13) is contained in the above polynomial system. Using the series expansions for the auxiliary variables
u ( z ) = i = 0 p u i z i , V ( z ) = i = 0 p V i z i , S ( z ) = i = 0 p S i z i , Q ( z ) = i = 0 p Q i z i ,
and substituting in the above polynomial system, leads to the following recursive relations for the coefficients
u i + 1 = V i i + 1 V i + 1 = 1 2 ( i + 1 ) j = 0 i S i j p = 0 j k = 0 p V k V p k V j p S i + 1 = 6 Ω m Q i i + 1 Q i = 0 : i > 1 1 : i = 1 .
According to Theorem 1 the interval of convergence for the above series solution to the polynomial system (18) is given by | z | < 1 / M where M = ( k 1 ) Σ g α k 1 . In our case k = 4 , Σ g = 6 Ω m , and α = 3 Ω m such that M = 6 ( 3 Ω m ) 4 . So for example taking Ω m = 0.28 we get 1 / M = 0.33 and so we choose a slightly tighter bound | z | < 0.25 to make sure that the series approximations converge on this interval. Then computing the coefficients by using the above recursive relations gives a series approximation valid on the interval [ 0 , 0.25 ] . As described further in the Appendix A, this process can be carried on until the full redshift range is covered, yielding the piecewise convergent approximate solution shown in Table 1. For computational simplicity we have chosen to limit ourselves to third-order approximations and limit the full redshift range to 0 z 2.5 . The order of the approximation and the redshift range can be increased as required. Figure 1a shows the luminosity distance function H 0 d L ( z ) / c = ( 1 + z ) u ( z ) for the Λ CDM model obtained via the multistage PSM in the red-shift range 0 z 2.5 , for Ω m = 0.28 . The piecewise convergent luminosity distance has an excellent agreement with the solution obtained via numerical integration of the integral in (12).
To determine the accuracy of the approximate luminosity distance we define the absolute relative percentage error by
Δ E = d L approx d L num d L num × 100 % ,
where d L approx is the luminosity distance values calculated by the PSM and d L num are the values obtained via numerical integration of the integral in (12). This is shown in Figure 1b.

4.2. CPL Cosmological Model

For the CPL model [40,41,42] the function W ( z ) in (15) is given by
W ( z ) = Ω m ( 1 + z ) 3 + ( 1 Ω m ) ( 1 + z ) 3 ( 1 + ω 0 + ω 1 ) e 3 ω 1 z 1 + z ,
which now contains the additional parameters ω 0 and ω 1 . The case ω 0 = 1 , ω 1 = 0 corresponds to the Λ CDM. Introducing the auxiliary variables
V ( z ) = u ( z ) ; S ( z ) = W ( z ) ; P ( z ) = e 3 ω 1 z / ( 1 + z ) ; R ( z ) = ( 1 + z ) ( 3 ω 0 + 3 ω 1 ) ; Q ( z ) = ( 1 + z ) ; Z ( z ) = 1 / Q ( z ) ,
the differential equation in (15) can be projected into the polynomial system
u = V ( z ) V = 1 2 S ( z ) V 3 ( z ) S = 6 Ω m Q ( z ) + 3 ( 1 + ω 0 + ω 1 ) ( 2 + 3 ω 0 + 3 ω 1 ) ( 1 Ω m ) P ( z ) Q ( z ) R ( z ) 6 ω 1 ( 2 + 3 ω 0 + 3 ω 1 ) ( 1 Ω m ) P R + 9 ω 1 2 ( 1 Ω m ) P ( z ) R ( z ) Z ( z ) P = 3 ω 1 P ( z ) Z ( z ) 2 R = 3 ( ω 0 + ω 1 ) R ( z ) Z ( z ) Z = Z 2 ( z ) Q = 1 ,
with the initial conditions u ( 0 ) = 0 , V ( 0 ) = 1 , S ( 0 ) = 3 ( 1 + ω 0 ( 1 Ω m ) ) , P ( 0 ) = 1 , R ( 0 ) = 1 , Z ( 0 ) = 1 , Q ( 0 ) = 1 . Using the series expansions
u ( z ) = i = 0 p u i z i , V ( z ) = i = 0 p V i z i , S ( z ) = i = 0 p S i z i , P ( z ) = i = 0 p P i z i , R ( z ) = i = 0 p R i z i , Q ( z ) = i = 0 p Q i z i , Z ( z ) = i = 0 p Z i z i .
and substituting in the above polynomial system, gives the following recursive relations for the coefficients
u i + 1 = V i i + 1 V i + 1 = 1 2 ( i + 1 ) j = 0 i S i j p = 0 j k = 0 p V k V p k V j p S i + 1 = 1 i + 1 6 Ω m Q i + 3 ( 1 + ω 0 + ω 1 ) ( 2 + 3 ω 0 + 3 ω 1 ) ( 1 Ω m ) j = 0 i k = 0 j P k Q j k R i j 6 ω 1 ( 2 + 3 ω 0 + 3 ω 1 ) ( 1 Ω m ) j = 0 i P j R i j + 9 ω 1 2 ( 1 Ω m ) j = 0 i k = 0 j P k R j k Z i j P i + 1 = 3 ω 1 i + 1 j = 0 i k = 0 j Z k Z j k P i j R i + 1 = 3 ( ω 0 + ω 1 ) i + 1 j = 0 i R j Z i j Z i + 1 = 1 i + 1 j = 0 i Z j Z i j Q i = 0 : i > 1 1 : i = 1 .
Again by Theorem 1 the interval of convergence is given by | z | < 1 / M where M = ( k 1 ) Σ g α k 1 . We use the values of the parameters for the CPL model obtained in Ref. [48], namely Ω m = 0.311 , ω 0 = 0.937 , and ω 1 = 0.064 . Then k = 4 , Σ g = 2.619 , α = 1.06322 such that 1 / M = 0.106 . Hence we start by computing the coefficients of the series approximations by using the above recursive relations in the subinterval [ 0 , 0.1 ] and again we consider third order approximations. These are shown in the Appendix A. As in the case of the Λ CDM, the process can be continued to cover the required redshift range, thereby yielding the piecewise convergent approximation shown in Table 2. Figure 2a below shows the luminosity distance function for the CPL cosmological model obtained via the multistage PSM in the red-shift range 0 z 1.2 , for Ω m = 0.311 , ω 0 = 0.937 , and ω 1 = 0.064 . Again, as in the previous case the obtained piecewise convergent luminosity distance has an excellent agreement with the solution obtained via numerical integration of the integral in (12). This is evidenced by the very small relative percentage error obtained from (21) as shown in Figure 2b.

5. Discussion and Conclusions

In this paper, we have used the Parker–Sochacki method to obtain an analytical approximation of the luminosity distance for the Λ CDM and CPL cosmological models. This is done by solving the nonlinear differential equation satisfied by the luminosity distance function by first projecting it into a polynomial system and then implementing the PSM in a multistage fashion for the redshift subintervals so that the result is a piecewise convergent luminosity distance function on the entire required redshift range, which has an excellent agreement with the numerical solutions. The overall accuracy is controlled by the length of these subintervals and the order of the approximation used which can be increased for a better overall accuracy. As mentioned in the Introduction, there are other methods proposed in the literature for obtaining analytical approximations of the luminosity function by solving the differential Equation (15). For example, Shchigolev applied the variational iteration method (VIM) [18], the homotopy perturbation method (HPM) [20] and the Daftardar–Jafari method (DFM) [21] to obtain analytical approximations for d L ( z ) for the Λ CDM model. These methods yield fairly accurate approximations for low red-shifts z < 1 , but the accuracy is nowhere near that obtained by using the PSM as can be seen by comparing Figure 1 in Refs. [18,20,21] with Figure 1 in this paper. The relative percentage errors of the luminosity distance approximations for the PSM are still significantly less than those obtained by applying a variant of the HPM, called the HPM–Padé technique which was proposed recently by Yu et al. [19] as an improvement to the approximations generated by the HPM in [20]. In fact, for the Λ CDM model with the same value of the cosmological parameter Ω m = 0.28 and the same redshift range 0 z 2.5 , the HPM–Padé yields a maximum relative percentage error Δ E ∼4% for the best approximation over this redshift range. In our case Δ E < 0.1 % , and this is comparable or even better than some of the earlier fitting formulae (such as those developed by Pen [3], Wickramasinghe and Ukwatta [4], Liu et al. [5] and Baes et al. [10]). More importantly, unlike previous approximations obtained via the VIM, HPM, DFM, and HPM–Padé, the relative percentage error Δ E in the approximated luminosity distance does not increase with redshift. The multi-stage approach in the implementation of the PSM makes it possible to apply it for even higher redshift ranges without compromising accuracy.
One can argue that unlike the various analytic approximations obtained by solving (15), the direct fitting formulae mentioned above are much easier to apply and have a shorter computation time, and they are valid over a very wide range of redshift. This is undoubtedly correct. However, a major drawback of these simple formulae is that they were specifically developed for the Λ CDM model and cannot be used for other cosmological models. On the other hand, as shown in this paper and in previous articles, analytic approximations for the luminosity distance obtained by solving the differential Equation (15) can be easily applied to various cosmological models; even more complex models than the CPL model analyzed here. This is particularly important in modern cosmology when testing and investigating new dark energy models in alternative theories of gravity. The need for an efficient and precise calculation of the luminosity distance in this case is paramount to uncover any subtle differences with the standard Λ CDM cosmology [49].

Funding

This research received no external funding.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A. Approximate Solutions

In this appendix, we obtain the series approximations to the polynomial system of equations for the Λ CDM and CPL cosmological models given by (18) and (24) respectively, by dividing the entire redshift range into subintervals as described in Section 4. The series expansions for the variable u ( z ) (which is directly related to the luminosity via (13)) in these subintervals gives the piecewise convergent approximation shown in Table 1 and Table 2.

Appendix A.1. Λ CDM Cosmological Model

From the recursive relations (20) for the first subinterval z [ 0 , 0.25 ] , it follows that
u ( 1 ) ( z ) = z 0.21 z 2 0.0518 z 3 V ( 1 ) ( z ) = 1 0.42 z 0.1554 z 2 + 0.2040 z 3 S ( 1 ) ( z ) = 0.84 ( 1 + z ) 2 Q ( 1 ) ( z ) = 1 + z ,
where the superscript (1) in the above expressions indicates that these series approximations are valid on the first subinterval. The PSM algorithm is applied in a multistage approach as explained in Section 2. Therefore on the second subinterval [ 0.25 , 0.50 ] the PSM algorithm is applied again by solving the same polynomial system in (18), but now with the initial conditions u ( 2 ) ( 0 ) = u ( 1 ) ( 0.25 ) , V ( 2 ) ( 0 ) = V ( 1 ) ( 0.25 ) , S ( 2 ) ( 0 ) = S ( 1 ) ( 0.25 ) and Q ( 2 ) ( 0 ) = Q ( 1 ) ( 0.25 ) . This yields:
u ( 2 ) ( z ) = 0.2361 + 0.8885 ( z 0.25 ) 0.2301 ( z 0.25 ) 2 0.0035 ( z 0.25 ) 3 V ( 2 ) ( z ) = 0.8885 0.4603 ( z 0.25 ) 0.0106 ( z 0.25 ) 2 + 0.1653 ( z 0.25 ) 3 S ( 2 ) ( z ) = 1.3125 + 2.1 ( z 0.25 ) + 0.84 ( z 0.25 ) 2 Q ( 2 ) ( z ) = 1 + z .
This process is repeated until the required full redshift range [ 0 , 2.5 ] is covered. Listing the obtained u ( i ) ( z ) as shown in Table 1 gives a piecewise convergent approximation for the luminosity distance function d L ( z ) valid on this redshift range.

Appendix A.2. CPL Cosmological Model

On the first subinterval z [ 0 , 0.1 ] , the recursive relations in (26) give
u ( 1 ) ( z ) = z 0.2658 z 2 0.0164 z 3 V ( 1 ) ( z ) = 1 0.5316 z 0.0493 z 2 + 0.2454 z 3 S ( 1 ) ( z ) = 1.0632 + 1.8927 z + 0.8016 z 2 + 0.2330 z 3 P ( 1 ) ( z ) = 1 0.192 z + 0.2104 z 2 0.2300 z 3 R ( 1 ) ( z ) = 1 2.619 z + 4.7391 z 2 7.2966 z 3 Z ( 1 ) ( z ) = 1 z + z 2 z 3 Q ( 1 ) ( z ) = 1 + z .
As was done in the previous case the PSM is applied in a multistage approach, and so for the next subinterval [ 0.1 , 0.2 ] , the polynomial system in (24) is solved with the initial conditions u ( 2 ) ( 0 ) = u ( 1 ) ( 0.1 ) , V ( 2 ) ( 0 ) = V ( 1 ) ( 0.1 ) , S ( 2 ) ( 0 ) = S ( 1 ) ( 0.1 ) , P ( 2 ) ( 0 ) = P ( 1 ) ( 0.1 ) , R ( 2 ) ( 0 ) = R ( 1 ) ( 0.1 ) , Z ( 2 ) ( 0 ) = Z ( 1 ) ( 0.1 ) and Q ( 2 ) ( 0 ) = Q ( 1 ) ( 0.1 ) . This gives
u ( 2 ) ( z ) = 0.0973 + 0.9466 ( z 0.1 ) 0.2673 ( z 0.1 ) 2 + 0.0055 ( z 0.1 ) 3 V ( 2 ) ( z ) = 0.9466 0.5347 ( z 0.1 ) + 0.0164 ( z 0.1 ) 2 + 0.1924 ( z 0.1 ) 3 S ( 2 ) ( z ) = 1.2607 + 2.0588 ( z 0.1 ) + 0.8550 ( z 0.1 ) 2 + 0.1328 ( z 0.1 ) 3 P ( 2 ) ( z ) = 0.9827 0.1559 ( z 0.1 ) + 0.1540 ( z 0.1 ) 2 0.1520 ( z 0.1 ) 3 R ( 2 ) ( z ) = 0.7782 1.8526 ( z 0.1 ) + 3.0473 ( z 0.1 ) 2 4.2648 ( z 0.1 ) 3 Z ( 2 ) ( z ) = 0.909 0.8263 ( z 0.1 ) + 0.7511 ( z 0.1 ) 2 0.6827 ( z 0.1 ) 3 Q ( 2 ) ( z ) = 1 + z .
The process is repeated until required redshift range is covered. Then the functions u ( i ) ( z ) represent the piecewise convergent approximation for the luminosity distance function d L ( z ) as shown in Table 2.

References

  1. Lin, W.; Mack, K.J.; Hou, L. Investigating the Hubble constant tension: Two numbers in the standard cosmological model. Astrophys. J. Lett. 2020, 904, L22. [Google Scholar] [CrossRef]
  2. Fanizza, G. Precision Cosmology and Hubble tension in the era of LSS surveys. arXiv 2021, arXiv:2110.15272v2. [Google Scholar]
  3. Pen, U.-L. Analytical fit to the luminosity distance for flat cosmologies with a cosmological constant. Astrophys. J. Supp. Ser. 1999, 120, 49–50. [Google Scholar] [CrossRef] [Green Version]
  4. Wickramasinghe, T.; Ukwatta, T.N. An analytical approach for the determination of the luminosity distance in a flat universe with dark energy. Mon. Not. R. Astron. Soc. 2010, 206, 548–550. [Google Scholar] [CrossRef] [Green Version]
  5. Liu, D.-Z.; Ma, C.; Zhang, T.-J.; Yang, Z. Numerical strategies of computing the luminosity distance. Mon. Not. R. Astron. Soc. 2011, 412, 2685–2688. [Google Scholar] [CrossRef] [Green Version]
  6. Adachi, M.; Kasai, M. An analytical approximation of the luminosity distance in flat cosmologies with a cosmological constant. Prog. Theor. Phys. 2012, 127, 145–152. [Google Scholar] [CrossRef] [Green Version]
  7. Eisenstein, D.J. An analytic expression for the growth function in a flat universe with a cosmological constant. arXiv 1997, arXiv:astro-ph/9709054. [Google Scholar]
  8. Mészáros, A.; Řípa, J. A curious relation between the flat cosmological model and the elliptic integral of the first kind. Astron. Astrophys. 2013, 556, A13. [Google Scholar] [CrossRef] [Green Version]
  9. Zaninetti, L. An analytical solution in the complex plane for the luminosity distance in flat cosmology. J. High. Energy Phys. 2016, 2, 581–586. [Google Scholar] [CrossRef] [Green Version]
  10. Baes, M.; Camps, P.; van de Putte, D. Analytical expressions and numerical evaluation of the luminosity distance in a flat cosmology. Mon. Not. R. Astron. Soc. 2017, 468, 927–930. [Google Scholar] [CrossRef] [Green Version]
  11. Clarkson, C.; Zunckel, C. Direct reconstruction of dark energy. Phys. Rev. Lett. 2010, 104, 211301. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  12. Capozziello, S.; D’Agostino, R.; Luongo, O. High-redshift cosmography: Auxiliary variables versus Pad polynomials. Mon. Not. R. Astron. Soc. 2020, 494, 2576–2590. [Google Scholar] [CrossRef]
  13. Cattoën, C.; Visser, M. The Hubble series: Convergence properties and redshift variables. Class. Quantum Grav. 2007, 24, 5985–5998. [Google Scholar] [CrossRef]
  14. Scolnic, D.M.; Jones, D.O.; Rest, A.; Pan, Y.C.; Chornock, R.; Foley, R.J.; Huber, M.E.; Kessler, R.; Narayan, G.; Riess, A.G.; et al. The complete light-curve sample of spectroscopically confirmed SNe Ia from Pan-STARRS1 and cosmological constraints from the combined Pantheon sample. Astrophys. J. 2018, 859, 101. [Google Scholar] [CrossRef]
  15. Zhou, Y.-N.; Liu, D.-Z.; Zou, X.-B.; Wei, H. New generalizations of cosmography inspired by the Padé approximant. Eur. Phys. J. C 2016, 76, 281. [Google Scholar] [CrossRef] [Green Version]
  16. Li, S.-Y.; Li, Y.-L.; Zhang, T.-J.; Zhang, T. Model-independent determination of cosmic curvature based on the Padé approximation. Astrophys. J. 2019, 887, 36. [Google Scholar] [CrossRef]
  17. Wei, H.; Yan, X.-P.; Zhou, Y.-N. Cosmological applications of Padé approximant. J. Cosmol. Astropart. Phys. 2014, 1401, 045. [Google Scholar] [CrossRef] [Green Version]
  18. Shchigolev, V.K. Analytic approximation of luminosity distance in cosmology via Variational Iteration Method. Univ. J. Comp. Math. 2017, 5, 68–74. [Google Scholar] [CrossRef]
  19. Yu, B.; Zhang, J.-C.; Zhang, T.-J.; Zhang, T. A new analytical approximation of luminosity distance by optimal HPM-Padé technique. Phys. Dark Univ. 2021, 31, 100772. [Google Scholar] [CrossRef]
  20. Shchigolev, V.K. Calculating luminosity distance versus redshift in FLRW cosmology via Homotopy Perturbation Method. Grav. Cosmo. 2017, 23, 142–148. [Google Scholar] [CrossRef] [Green Version]
  21. Shchigolev, V.K. A simple analytic approximation of luminosity distance in FLRW cosmology using Daftardar-Jafari Method. Comp. Sci. Info. Tech. 2020, 8, 43–49. [Google Scholar] [CrossRef]
  22. Parker, G.E.; Sochacki, J.S. Implementing the Picard iteration. Neural Parallel Sci. Comput. 1996, 4, 97–112. [Google Scholar]
  23. Parker, G.E.; Sochacki, J.S. A Picard-McLaurin theorem for initial value PDEs. Abstr. Appl. Anal. 2000, 5, 47–63. [Google Scholar] [CrossRef] [Green Version]
  24. Picard, E. Traite D’Analyse; Gauthier-Villars: Paris, France, 2013; Volume 3, pp. 1922–1928. [Google Scholar]
  25. Abdelrazik, I.M.; Elkaranshawy, H.A. Extended Parker Sochacki method for Michaelis-Menten enzymatic reaction model. Anal. Biochem. 2016, 496, 50–54. [Google Scholar] [CrossRef] [PubMed]
  26. Abdelrazik, I.M.; Elkaranshawy, H.A.; Abdelrazek, A.M. Modified Parker-Sochacki method for solving nonlinear oscillators. Mech. Based Des. Struct. Mach. 2017, 45, 239–252. [Google Scholar] [CrossRef]
  27. Akindeinde, S.O. Parker-Sochacki method for the solution of convective straight fins problem with temperature-dependent thermal conductivity. Int. J. Nonlinear Sci. 2018, 25, 119–128. [Google Scholar]
  28. Akindeinde, S.O.; Okyere, E. New analytic technique for the solution of nth order nonlinear two-point boundary value problems. Br. J. Math. Comput. Sci. 2016, 15, 1–11. [Google Scholar] [CrossRef]
  29. Elkaranshawy, H.A.; Abdelrazek, A.M.; Ezzat, H.M. Power series solution to sliding velocity in three-dimensional multibody systems with impact and friction. Int. J. Math. Comput. Phys. Electr. Comput. Eng. 2015, 9, 585–588. [Google Scholar]
  30. Nakhjiri, N.; Villac, B. Modified Picard integrator for spacefight mechanics. J. Guid. Control Dyn. 2014, 37, 1625–1637. [Google Scholar] [CrossRef]
  31. Nurminskii, E.; Buryi, A. Parker-Sochacki method for solving systems of ordinary differential equations using graphics processors. Numer. Anal. Appl. 2011, 4, 223. [Google Scholar] [CrossRef]
  32. Ogundare, B.S.; Akindeinde, S.O.; Adewumi, A.O.; Aderogba, A.A. Improved Parker- Sochacki approach for closed form solution of enzyme catalyzed reaction model. J. Mod. Methods Numer. Math. 2017, 8, 90–98. [Google Scholar] [CrossRef] [Green Version]
  33. Rudmin, J.W. Application of the Parker-Sochacki Method to Celestial Mechanics; Technical Report; James Madison University: Harrisonburg, VA, USA, 1998. [Google Scholar]
  34. Singarimbun, A.; Fujimitsu, Y.; Djamal, M.; Dewi, R. Pressure transient modeling in geothermal reservoir by using Picard-Mclaurin iteration. Adv. Mater. Res. 2014, 1025–1026, 959–973. [Google Scholar]
  35. Stewart, R.D.; Bair, W. Spiking neural network simulation: Numerical integration with the Parker-Sochacki method. J. Comput. Neurosci. 2009, 27, 115–133. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  36. Stewart, R.D.; Bair, W. Spiking neural network simulation: Memory-optimal synaptic event scheduling. J. Comput. Neurosci. 2011, 30, 721–728. [Google Scholar] [CrossRef]
  37. Szynkiewicz, P. A novel GPU-enabled simulator for large scale spiking neural networks. J. Telecommun. Inf. Technol. 2016, 2, 34–42. [Google Scholar]
  38. Akindeinde, S.O. A new multistage technique for approximate analytical solution of nonlinear differential equations. Heliyon 2020, 6, e05188. [Google Scholar] [CrossRef]
  39. Carothers, D.C.; Parker, G.E.; Sochacki, J.S.; Warne, P.G. Some properties of solutions to polynomial systems of differential equations. Electron. J. Differ. Equ. 2005, 2005, 1–17. [Google Scholar]
  40. Chevallier, M.; Polarski, D. Accelerating universes with scaling dark matter. Int. J. Mod. Phys. D 2001, 10, 213–223. [Google Scholar] [CrossRef] [Green Version]
  41. Linder, E.V. Exploring the expansion history of the Universe. Phys. Rev. Lett. 2003, 90, 091301. [Google Scholar] [CrossRef] [Green Version]
  42. Linder, E.V. Probing gravitation, dark energy, and acceleration. Phys. Rev. D 2004, 70, 023511. [Google Scholar] [CrossRef] [Green Version]
  43. Guenther, J.; Wolf, M. An adaptive, highly accurate and efficient, Parker-Sochacki algorithm for numerical solutions to initial value ordinary differential equation systems. SIAM Undergrad. Res. Online 2019, 12, 257–281. [Google Scholar] [CrossRef]
  44. Pouly, A.; Graça, D.S. Computational complexity of solving polynomial differential equations over unbounded domains. Theor. Comp. Sci. 2016, 626, 67–82. [Google Scholar] [CrossRef] [Green Version]
  45. Riess, A.G.; Filippenko, A.V.; Challis, P.; Clocchiatti, A.; Diercks, A.; Garnavich, P.M.; Gilliland, R.L.; Hogan, C.J.; Jha, S.; Kirshner, R.P.; et al. Observational evidence from supernovae for an accelerating universe and a cosmological constant. Astron. J. 1998, 116, 1009. [Google Scholar] [CrossRef] [Green Version]
  46. Perlmutter, S.; Aldering, G.; Goldhaber, G.; Knop, R.A.; Nugent, P.; Castro, P.G.; Deustua, S.; Fabbro, S.; Goobar, A.; Groom, D.E.; et al. Measurements of omega and lambda from 42 high-redshift supernovae. Astrophys. J. 1999, 517, 565. [Google Scholar] [CrossRef]
  47. Weinberg, S. Gravitation and Cosmology: Principles and Applications of The General Theory of Relativity; John Wiley Press: New York, NY, USA, 1972. [Google Scholar]
  48. Solà, J.; Gómez-Valent, A.; de Pérez, J. Dynamical dark energy: Scalar fields and running vacuum. Mod. Phys. Lett. A 2017, 32, 1750054. [Google Scholar] [CrossRef] [Green Version]
  49. Nunes, R.C.; Pan, S.; Abreu, E.N.S.E.M.C. New observational constraints on f(R) gravity from cosmic chronometers. J. Cosmol. Astropart. Phys. 2017, 2017, 005. [Google Scholar] [CrossRef] [Green Version]
Figure 1. (a) The luminosity distance function ( H 0 / c ) d L ( z ) (solid curve) vs. the redshift z for the Λ CDM model with Ω m = 0.28 . The corresponding numerical solution obtained by integrating (12) is also shown as plotted points for discrete values of z. (b) The absolute relative percentage error of the luminosity distance Δ E vs. the redshift z.
Figure 1. (a) The luminosity distance function ( H 0 / c ) d L ( z ) (solid curve) vs. the redshift z for the Λ CDM model with Ω m = 0.28 . The corresponding numerical solution obtained by integrating (12) is also shown as plotted points for discrete values of z. (b) The absolute relative percentage error of the luminosity distance Δ E vs. the redshift z.
Universe 08 00300 g001
Figure 2. (a) The luminosity distance function ( H 0 / c ) d L ( z ) (solid curve) vs. the redshift z for the CPL model with Ω m = 0.311 , ω 0 = 0.937 , and ω 1 = 0.064 . The corresponding numerical solution obtained by integrating (12) is also shown as plotted points for discrete values of z. (b) The absolute relative percentage error of the luminosity distance Δ E vs. the redshift z for the same redshift range.
Figure 2. (a) The luminosity distance function ( H 0 / c ) d L ( z ) (solid curve) vs. the redshift z for the CPL model with Ω m = 0.311 , ω 0 = 0.937 , and ω 1 = 0.064 . The corresponding numerical solution obtained by integrating (12) is also shown as plotted points for discrete values of z. (b) The absolute relative percentage error of the luminosity distance Δ E vs. the redshift z for the same redshift range.
Universe 08 00300 g002
Table 1. Piecewise convergent approximate solution for u ( z ) with Ω m = 0.28 .
Table 1. Piecewise convergent approximate solution for u ( z ) with Ω m = 0.28 .
z  u ( z )
0 z 0.25   z 0.21 z 2 0.0518 z 3
0.25 z 0.50   0.2361 + 0.8885 ( z 0.25 ) 0.2301 ( z 0.25 ) 2 0.0035 ( z 0.25 ) 3
0.50 z 0.75   0.4438 + 0.7753 ( z 0.5 ) 0.2202 ( z 0.5 ) 2 + 0.0272 ( z 0.5 ) 3
2.25 z 2.5   1.3133 + 0.3112 ( z 2.25 ) 0.0245 ( z 2.25 ) 2 0.0336 ( z 2.25 ) 3
Table 2. Piecewise convergent approximate solution for u ( z ) with Ω m = 0.311 , ω 0 = 0.937 , and ω 1 = 0.064 .
Table 2. Piecewise convergent approximate solution for u ( z ) with Ω m = 0.311 , ω 0 = 0.937 , and ω 1 = 0.064 .
z  u ( z )
0 z 0.1   z 0.2658 z 2 0.0164 z 3
0.1 z 0.2   0.0973 + 0.9466 ( z 0.1 ) 0.2673 ( z 0.1 ) 2 + 0.0055 ( z 0.1 ) 3
0.2 z 0.3   0.1893 + 0.8935 ( z 0.2 ) 0.2631 ( z 0.2 ) 2 + 0.0222 ( z 0.2 ) 3
1.1 z 1.2   0.8088 + 0.5194 ( z 1.1 ) 0.1481 ( z 1.1 ) 2 + 0.0390 ( z 1.1 ) 3
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Sultana, J. A New Analytic Approximation of Luminosity Distance in Cosmology Using the Parker–Sochacki Method. Universe 2022, 8, 300. https://doi.org/10.3390/universe8060300

AMA Style

Sultana J. A New Analytic Approximation of Luminosity Distance in Cosmology Using the Parker–Sochacki Method. Universe. 2022; 8(6):300. https://doi.org/10.3390/universe8060300

Chicago/Turabian Style

Sultana, Joseph. 2022. "A New Analytic Approximation of Luminosity Distance in Cosmology Using the Parker–Sochacki Method" Universe 8, no. 6: 300. https://doi.org/10.3390/universe8060300

APA Style

Sultana, J. (2022). A New Analytic Approximation of Luminosity Distance in Cosmology Using the Parker–Sochacki Method. Universe, 8(6), 300. https://doi.org/10.3390/universe8060300

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop