Comoving Distance in LCDM - Understanding an Approximation

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:
Hi, I had an exam and I completely messed up a problem. Especially one part which was necessary for the rest of the problem. Basically, I have a wormhole metric: $$(ds)^2 = -(dt)^2 + (dr)^2 + (r^2 + b^2)( (d\theta)^2 + sin^2 \theta (d\phi)^2 )$$ Where ##b=1## with an orbit only in the equatorial plane. We also know from the question that the orbit must satisfy this relationship: $$\varepsilon = \frac{1}{2} (\frac{dr}{d\tau})^2 + V_{eff}(r)$$ Ultimately, I was tasked to find the initial...
The value of H equals ## 10^{3}## in natural units, According to : https://en.wikipedia.org/wiki/Natural_units, ## t \sim 10^{-21} sec = 10^{21} Hz ##, and since ## \text{GeV} \sim 10^{24} \text{Hz } ##, ## GeV \sim 10^{24} \times 10^{-21} = 10^3 ## in natural units. So is this conversion correct? Also in the above formula, can I convert H to that natural units , since it’s a constant, while keeping k in Hz ?
Back
Top