Comoving Distance in LCDM - Understanding an Approximation

Click For Summary
SUMMARY

The discussion focuses on calculating the comoving distance in the ΛCDM model using the integral $$\chi = c\int_0^z \frac{dz}{H(z)}$$. The user encounters difficulties with manual integration and opts for an approximation method as outlined in textbooks like Weinberg and Barbara. The approximation involves expanding the scale factor $$a(t)$$ around the present time $$t_0$$, leading to a series expansion for $$a(t)^{-1}$$. The discrepancy arises from the inclusion of the jerk parameter $$j_0$$ in the user's calculations, which is not present in Barbara's formulation.

PREREQUISITES
  • Understanding of the ΛCDM cosmological model
  • Familiarity with the Friedmann equations and Hubble parameter $$H(z)$$
  • Knowledge of Taylor series expansions and their applications in cosmology
  • Experience with numerical integration techniques
NEXT STEPS
  • Study the derivation of the Hubble parameter $$H(z)$$ in the ΛCDM model
  • Learn about Taylor series expansions in the context of cosmological scale factors
  • Explore numerical integration methods for cosmological distance calculations
  • Review the differences between jerk, snap, and other cosmological parameters in the context of cosmological models
USEFUL FOR

Astronomers, cosmologists, and physics students interested in advanced cosmological models and distance calculations in a flat universe.

Arman777
Insights Author
Gold Member
Messages
2,163
Reaction score
191
Homework Statement
Finding the Comoving Distance for LCDM universe.
Relevant Equations
$$\chi(t) = \int_{t_e}^{t_0}\frac{dt}{a(t)}$$
I am trying to find the comoving distance,

$$\chi = c\int_0^z \frac{dz}{H(z)}$$ for the ##\Lambda##CDM model (spatially flat universe, containing only matter and ##\Lambda##).

$$H^2 = H_0^2[\Omega_{m,0}(1+z)^3 + \Omega_{\Lambda, 0}]$$

When I put this into integral I am getting,

$$\chi = \frac{c}{H_0}\int_0^z \frac{dz}{\sqrt{\Omega_{m,0}(1+z)^3 + \Omega_{\Lambda, 0}}}$$

which seems that I cannot take the integral manually (only numeric integration is possible).

So I have decided to approximate the solution, as its done in many textbooks (including Weinberg and Barbara).

So I have written $$\chi(t) = \int_{t_e}^{t_0}\frac{dt}{a(t)}$$ and then I have expanded the ##a(t)## around ##t_0##.

$$a(t) \simeq 1 + H_0(t-t_0) - \frac{q_0H_0^2}{2}(t-t_0)^2 + \frac{j_0H_0^3}{6}(t-t_0)^3$$

from here I need to find $$a(t)^{-1}$$ so that I can put that inside the integral.
I have seen that $$1/x \simeq \sum_{n=0}^{\infty}(-1)^n(x-1)^n$$.

So by similar logic I should get

$$a(t)^{-1} \simeq 1 - H_0(t-t_0) + \frac{q_0H_0^2}{2}(t-t_0)^2 - \frac{j_0H_0^3}{6}(t-t_0)^3$$

However, it seems its not the case. In Barbara (the calculations do not include ##j_0##) its written as,

$$a(t)^{-1} \simeq 1 - H_0(t-t_0) + (1+q_0/2)\frac{H_0^2}{2}(t-t_0)^2$$

I just don't understand how this can be possible...I am doing a mistake somewhere but I don't know where.

Thanks
 
Physics news on Phys.org
a(t)^{-1}=\frac{1}{1+A}=1-A+A^2+..
where
A=H_0(t-t_0)-qH_0^2/2(t-t_0)^2+..
First order comes from A only . Second order comes from A and ##A^2## which is ##H_0^2(t-t_0)^2##. Similarly third order comes from A, A^2 and A^3.
 
Last edited:
  • Like
Likes   Reactions: Arman777

Similar threads

  • · Replies 2 ·
Replies
2
Views
2K
  • · Replies 21 ·
Replies
21
Views
3K
  • · Replies 2 ·
Replies
2
Views
755
  • · Replies 1 ·
Replies
1
Views
2K
  • · Replies 3 ·
Replies
3
Views
2K
  • · Replies 6 ·
Replies
6
Views
4K
Replies
4
Views
2K
  • · Replies 8 ·
Replies
8
Views
2K
  • · Replies 1 ·
Replies
1
Views
2K
  • · Replies 4 ·
Replies
4
Views
2K