Access the full text.
Sign up today, get DeepDyve free for 14 days.
Loading next page...
/lp/arxiv-cornell-university/computing-the-luminosity-distance-via-optimal-homotopy-perturbation-xofZl8D8pN
References
References for this paper are not available at this time. We will be adding them shortly, thank you for your patience.
- ISSN
- 2212-6864
- eISSN
- ARCH-3330
- DOI
- 10.1016/j.dark.2020.100734
- Publisher site
- See Article on Publisher Site
Abstract
We propose a new algorithm for computing the luminosity distance in the at universe with a cosmological constant based on Shchigolev's homotopy perturba- tion method, where the optimization idea is applied to prevent the arbitrariness of initial value choice in Shchigolev's homotopy. Compared with the some exist- ing numerical methods, the result of numerical simulation shows that our algo- rithm is a very promising and powerful technique for computing the luminosity distance, which has obvious advantages in computational accuracy,computing eciency and robustness for a given Keywords: Distance scale, Numerical-cosmology, Optimal homotopy perturbation method 1. Introduction Considerable attentions have been paid to the numerical computation of cosmological distances in recent research[1,2,3,4,5,6,7]. One of the most funda- mental distances in cosmology is the luminosity distance d , which depends on the redshift and the cosmological parmeters. As is well known, the luminosity Fully documented templates are available in the elsarticle package on CTAN. Corresponding author Email address: tjzhang@bnu.edu.cn (Tong-Jie Zhang) Preprint submitted to Elsevier January 13, 2021 arXiv:2101.04110v1 [astro-ph.CO] 10 Jan 2021 distance can not be expressed in terms of a simple analytical function of redshift and the parameters of the underlying cosmological models. Taking the general lambda cold dark matter (CDM), the luminosity distance is only expressible in terms of transcendental functions or elliptic integral functions. Although the computing speed of the state-of-the-art computers is very fast, the distance anal- ysis involving in many studies (e.g. supernova cosmology) still require extensive computions of the numerical integrals, and is time-comsuming[2,6]. Therefore, accurate and ecient strategies of computing luminosity distance are vitial in modern cosmology. Up to now, many analytical and numerical approaches have been proposed for avoiding the diculty of heavy computation pressure. Pen(hereafter Pen99) proposed a simple algebraic tting formula that has a global relative error of less than 4% [1] . Liu et al. [2] presented two ecient numerical strategies of calculating elliptic integrals for luminosity distance in at CDM models. Wickramasinghe and Ukwatta(hereafter WU10) [3] obtained a di erent approx- imate expression of luminosity distance, which has a smaller relative error than Pen99. Hao et al. [4] used the Pad e approximant technique to obtain an di erent analytical expression. Maarten et al. used hypergeometric functions to derive an another analytical formula of luminosity distance for at CDM models [6] .By solving a certain di erential equation based on the homotopy perturba- tion method (HPM), a new way to calculate luminosity distance is proposed by Shchigolev (hereafter Shch17) [7], which is di erent from that of previous methods. In general, these evaluation methods of luminosity distance can be classi ed into two kinds; one is the method based on simpli cation of elliptic in- tegrals in luminositydistance and the other is the solution of a ceatain nonlinear di erential equation which the luminosity distance satis es to. Based on methods of simpli cation of elliptic integrals simple and ecient formulas over a large range of redshift usually can be obtained. However, meth- ods of simpli cation of elliptic integrals have no advantage in calculating pre- cision in certain small redshift range. In contrast, a formula for calculating the luminosity distance with less relative error over a small range of redshift 2 can be derived by solving a certain nonlinear di erential equation based on the homotopy perturbation method. The disadvantage of the solution of a ceatain nonlinear di erential equation is that the redshift range of calculating luminos- ity distance is slightly smaller than that of the method based on simpli cation of elliptic integrals. The homotopy perturbation method that takes full advantage of homotopy and perturbation in solving nonlinear di erential equations, was rst proposed by He [8,9]. It has been shown that a wide range of nonlinear di erential equa- tions can be solved to yield an highly accurate solution by homotopy perturba- tion method with one or two iterations.Later, several modi cations of homotopy perturbation method have been quickly proposed, such as optimal homotopy perturbation method[10,11], which can get a reliable approach to nonlinear problems,and optimal homotopy perturbation method coupled with the least squares technology[12,13],and so on.In addition, it is noted that the selection of initial value in Shchigolev's homotopy is arbitrary, meaning that optimization to the selection of initial value can be performed on his method. In this paper, a new algorithm is proposed for the purpose of getting a more accurate and ecient expression of the luminosity distance over a relatively large range of redshift. The rest of the paper is organized as follows. The di erential equa- tion which the luminosity distance in at CDM models satis es to is built in Section 2. In Section 3, the modi ed algorithm based on Shchigolev's method is presented. The performance of the proposed algorithm against some exist- ing exact methods are tested in Section 4. Finally, the conclusion is given in Section 5. 2. Di erential equation of luminosity distance for at Lambda cold dark matter(CDM) universe As mentioned in Liu2011, the luminosity distance d in the spatially at Lambda cold dark matter universe is given by c(1 + z) dt d (z) = p ; (1) 0 0 (1 + t) + 3 where are the energy densities corresponding to the matter and cosmo- logical constant, respectively: = 1. Following the notation in Shch16, we introduce the notation W (z) = (1 + z) + ; W (z)j = 1: (2) m z=0 Then we rewrite Eq.(1)as: c(1 + z) dt d (z) = p : (3) W (t) 0 0 d (z)H dt L 0 = : (4) c(1 + z) W (t) By di erentiating the Eq.(4), one can get d d (z)H 1 L 0 = p : (5) dz c(1 + z) W (z) For simplicity sake, we introduce: d (z)H L 0 1 + z = x; u(x) = ; (6) cx and rewrite equation (5) as du 1 j = p (7) x=1 dx W (x 1) Combining the Eq.(2), (4),(5),(6) and (7),we have u (x)j = 1; u(x)j = 0: (8) x=1 x=1 According to Eq.(7) and (8),the second derivative of u(x) is equal to d u(x) 1 3 dW (x 1) = W (x 1) : (9) dx 2 dx Combining the Eq.(8) and (9),one can obtain the Cauchy problem 00 1 0 0 3 0 u + W (x 1)u = 0; u (x)j = 1; u(x)j = 0: (10) x=1 x=1 where the superscript stands for the derivative with respect to x, and W (x 1)j = 1. x=1 4 3. Optimal homotopy perturbation method (OHPM)for calculation of luminosity distance Eq.(10) is a typical nonlinear second-order di erential equation, which can be accurately solved by integration method, and its exact result is formula (1). So we will use analytic method to obtain the approximate solution of this equation. In general,the nonlinear di erential equation can be solved very well to yield an highly accurate solution by homotopy perturbation method with one or two iterations.Using homotopy perturbation technique, nonlinear di erential equa- tions can be transformed into an in nite number of linear (simple)di erential equations. Now since the homotopy perturbation technique has become standard and concise, the reader can refer to[8,9] for its basic idea. Let us assume that a series in p can be used to represent the solution of Eq.(10).The speci c expression for the series inp is as follows: 2 3 u = u + pu + p u + p u + : : : ; (11) 0 1 2 3 where p 2 [0; 1] is an embedding parameter. Setting p = 1 results in the approximate solution of Eq.(10). In according with the procedure of the HPM, we build the homotopy as follows: 00 1 0 0 u + c + p W (x 1)u c = 0; (12) 1 1 where p 2 [0; 1] ,and the constant c is introduced into Eq.(10) with an initial condition u + c = 0 at p = 0. Substituting equation (11) into equation(12), and equating coecients with the identical powers of p, a set of equations are obtained as follows: p :u + c = 0; (13) 0 1 00 1 0 0 3 p :u + W (x 1)u c = 0; (14) 1 0 1 : : : : : : : : : : : : 5 According to Eq.(10), the initial conditions for u (x) can be set as follows: u j = 0; u j = 1; 0 x=1 0 x=1 u j = 0; u j = 1; (15) i x=1 i x=1 whrere i 1. By solving Eq.(13) with the initial conditions(15), one can obtain u = (x 1) (x 1) ; (16) Combining this equation and Eq.(14)with the initial conditions(15), we get 1 c 1 u = (x 1) + (x 1) W (t 1) 2 2 2 (17) [1 c (t 1)] f1 c [4(t 1) 3(x 1)]gdt; 1 1 Substituting Eq.(16) and Eq.(17) into Eq.(11), the approximate expression for luminosity distance with unknown constant c can be obtained c(1 + z) 3 7 3 2 6 d (z) = f0:01c z + (0:028c 0:042c )z 1 1 1 3 2 5 2 (18) + (0:021c 0:126c + 0:063c )z (0:105c 0:21c 1 1 1 1 1 4 3 2 + 0:035)z + (0:21c 0:14)z 0:21z + zg; where c is the speed of light,H is the Hubble constant, c is the unknown 0 1 constant, and z is redshift. Lastly and most importantly, the unknown constant c can be determined by minimizing the raltive error in the approximation of the luminosity distance: num d d 4E = (19) num num where d and d stand for the values of luminosity distance calculated form our approximate expression and the numerical method, respectively. The implementation of OHPM can be summarized as follows: Step1: According to the CDM model, obtain the di erential equation Eq.(10) that the luminosity distance d should satisfy to; 6 Step2: Build the homotopy Eq.(12); substitute Eq.(11) into Eq.(12), and then equate coecients with the identical powers of p; Step3: Solve the set of di erential equations consisting of Eq.(13) and Eq.(14) that obtained from Step2 with the initial conditions(15), and then get a set of equations with the unknown constant c ; Step4: According to Eq.(11),Eq.(16) and Eq.(17), one can obtain the ap- proximate expression for luminosity distance with an unknown constant c ; Step5: By minimizing Eq.(19) , the unknown constant c can be determined, and then we obtain the approximate expression for luminosity distance. 4. Performance of OHPM In this section, the performance of OHPM proposed in Sect.?? is assessed. The assessment is mainly carried out from two aspects: accuracy and eciency. In the at CDM models, we set = 0:28 and = 0:72,as an example. 4.1. Accuracy Figure 1 shows the comparisons of relative error percentages of approximate solutions to d (4E)for di erent values of c . By minimizing Eq.(19) , the L 1 unknown constant c is determined to be 0.44274 for the xed = 0:28, and 1 m then we obtain the approximate expression for luminosity distance. Using Python, one can obtain a series of numerical solutions to the Eq.(1),and the approximate solutions of Eq.(18). Compared with some existing methods, the relative error percentages of the numerical and approximate solutions of the same sample are given in Table 1. Seen from Table 1, a best approximation to the exact value of d for redshift range 0:1 < z 1:129 can be obtained by our method. For a xed = 0:28,it is clear from relative error percent- ages of d in Figure 2 that our method obviously outperforms some existing methods for 0 z 1:129. Relative error percentages of approximate solutions to d (4E)as a function of z for di erent is shown in Figure 3. From the L m Figure 4, for any redshift in 0 z 0:5 the relative error percentages of d are between-0.18% and 0.22% for within 0:26 0:30.The global error m m 7 Figure 1: The comparison of approximate solutions to d for di erent constants c . L 1 Table 1: Relative error percentages of approximate solutions to d (Errors %)in cases of OHPM, Shch17,Pen99 and WU10 z OHPM Shch17 Pen99 WU10 0.1 0.00351 0.00164 0.25934 0.25998 0.3 0.01067 0.04243 0.32274 0.16452 . 0.5 0.01192 0.17769 0.28412 0.11771 0.7 0.05888 0.41710 0.21090 0.91915 0.9 0.07540 0.71596 0.13634 0.07626 1.1 0.02791 0.98781 0.07311 0.06598 surface plot is shown in Figure 5.It shows that the error rst decreases and then increases when the variation of from 0.26 to 0.3, for 0:26 <0:28 the m m error rst decreases and then increases when the variation of from 0.26 to 0.3,and for 0:28 0:30 the error rst increases and then decreases when the variation of z from 0 to 0.5. 4.2. Eciency A comparison of the eciency of the some exact formulas to calculate the function of d is the main purpose of our numerical test. For this purpose, a sample of SN Ia redshifts is created based on the SNAP observation which has 1326 SN data points within 0.1<z <1.1 (Sha eloo et al)[14].We mock sample has the same redshift distribution as ducial SNAP, but is 100 times larger 8 Figure 2: Relative error percentages of approximate solutions to d (4E)as a function of z for = 0:72. Figure 3: Relative error percentages of approximate solutions to d (4E)as a function of z for di erent .The relative error percentages in d for 0 z 0:5 is ampli ed,which is shown m L in the inset. 9 Figure 4: The contour plot for the dietribution of relative error percentages using approximate solutions to d (4E) corresponding to within 0:26 0:30. The \pits" and \peaks"in m m the right region( 0:295) dominate the global error. Figure 5: The global error surface plot for approximate solutions to d (4E) corresponding to within 0:26 0:30. When the variation of from 0.26 to 0.3, the error rst m m m decreases and then increases. When the variation of z from 0 to 0.5, the error rst decreases and then increases for 0:26 <0:28, and the error rst increases and then decreases for 0:28 0:30. 10 Figure 6: A comparison of the running time of our method,Pen99 and WU10 for redshift range 0:1 < z 1:1. in data points than it.Because the calculation accuracy of Shch17 is relatively low, we no longer compare it when comparing the calculation eciency. So we conducted custom implementations of the methods from Pen99, WU10 and our method in the Python and used its time module. Each implementation of calculating d values from the created sample that contains 132600 SNe redshift points within0:1 z 1:1 is repeated 100 times. The histogram of the running time in cases of our method and two other methods is shown in Figure 6. Seen from the Figure 6, our method is obviously faster than that of Pen99 and WU10, and WU10 is slightly faster than Pen99.We note that the results of numerical test may vary depending on the compiler used and the hardware con guration. 5. Conclusions and discussion In this paper, a new algorithm for computing the luminositiy distance for at universes with a cosmological constan is proposed,which is named OHPM.The proposed algorithm integrates the optimization idea into homotopy perturbation method,where the modi ed method is applied to prevent the arbitrariness of initial value choice in Shchigolev's homotopy. 11 The results of numerical simulation indicate that OHPM has obvious advan- tages in computational accuracy. The relative error percentages is less than 0.08 percent error for redshit within 0 z 1:129 for the xed = 0:28. Figure 3,Figure 4 and Figure 5 indicate our algorithm has certain robustness for the dif- ferent . In the respects of enhancing computational eciency, our algorithm possesses great advantage. In addition, OHPM can be extended to other cos- mological models.Therefore, OHPM is a very promising and powerful technique to solve the calculation of luminosity distance in theoretical cosmology. Acknowledgments Bo Yu would like to thank Prof. Jin-Yu He for his kind help.This work was supported by National Key R&D Program of China (2017YFA0402600) and the National Science Foundation of China (Grants No. 11929301, 11573006). References [1] C. Clarkson and C. Zunckel, Direct Reconstruction of Dark Energy, Phys. Rev. Lett. 2010. [2] Pen, Ue-Li, et al., Analytical Fit to the Luminosity Distance for Flat Cos- mologies with a Cosmological Constant, Astron.Astrophys.Suppl.Ser. 1999. [3] M. Li, X.-D. Li, S. Wang, Y. Wang, Dark Energy, Communications in Theoretical Physics 56 (3) (2011) 525{604. [4] Liu, De-Zi and Ma, Cong and Zhang, Tong-Jie and Yang, Zhiliang, Numeri- cal strategies of computing the luminosity distance, Mon.Not.R.Astron.Soc. [5] Wickramasinghe, T. and Ukwatta, T. N., An analytical approach for the determination of the luminosity distance in a at universe with dark energy, Mon. Not. R.Astron.Soc.2010. [6] Wei, Hao and Yan, Xiao-Peng and Zhou, Ya-Nan, Cosmological applica- tions of Pad e approximant, J COSMOL ASTROPART P. 2014. 12 [7] Adachi, M. and Kasai, M.,An Analytical Approximation of the Luminosity Distance in Flat Cosmologies with a Cosmological Constant,Prog. Theor. Phys.2012. [8] Baes, Maarten and Camps, Peter and Van De Putte, Dries, Analytical ex- pressions and numerical Analytical expressions and numerical evaluation of the luminosity distance in a at cosmology, Mon. Not. R. Astron. Soc.2017. [9] Shchigolev, V. K., Calculating luminosity distance versus redshift in FLRW cosmology via homotopy perturbation method,Gravitation and Cosmol- ogy.2017. [10] He, J.,Homotopy perturbation technique,Computer Methods in Applied Mechanics and Engineering.1999. [11] He, Ji-Huan, Some Asymptotic Methods for Strongly Nonlinear Equations, International Journal of Modern Physics B.2006. [12] Sha eloo, Arman and Alam, Ujjaini and Sahni, Varun and Starobinsky, Alexei A., Smoothing supernova data to reconstruct the expansion history of the Universe and its age,Mon. Not. R. Astron. Soc.2006. [13] Nicolae Heris and Marinca, Vasile, Optimal Homotopy Perturbation Method for a Non-Conservative Dynamical System of a Rotating Electrical Machine,Ztschrift Fur Naturforschung A.2012. [14] Gupta, A. K. and Saha Ray, S., Comparison between homotopy pertur- bation method and optimal homotopy asymptotic method for the soliton solutions of Boussinesq{Burger equations, Computers and Fluids.2014. [15] Hayman Thabet and Subhash Kendre, Modi ed least squares homo- topy perturbation method for solving fractional partial di erential equa- tions,Malaya Journal of Matematik.2018. [16] Constantin Bota and Bogdan C~aruntu, Approximate analytical solutions of nonlinear di erential equations using the Least Squares Homotopy Pertur- bation Method, Journal of Mathematical Analysis and Applications.2016
Journal
Astrophysics – arXiv (Cornell University)
Published: Jan 10, 2021