View Single Post
  #1  
Old January 20th 20, 12:05 AM posted to sci.astro.research
Phillip Helbig (undress to reply)[_2_]
external usenet poster
 
Posts: 273
Default easy(?) cosmology question: what to integrate to get R(t)

When discussing Friedmann models, the usual approach is to start with
the Friedmann equation then express dr/dt as a function of R, Omega,
lambda, and so on, then re-arranging it to get an expression for dt as a
function of R (scale factor), Omega, and lambda, which can be integrated
to give the lookback time (or time since the big bang) as a function of
the scale factor. This can be re-arranged to express it as a function
of redshift, then one can compute distances as a function of redshift
and so on.

This is standard stuff.

But what about an expression which one can integrate to get R(t), the
scale factor as a function of time?

For special case, one can invert t(R) to get R(t). Of course, if I
numerically calculate t(R), I can invert it to get R(t). But what I
would like is an expression for dR as a function of t which I can
integrate from 0 to t to get R(t), just like I have an expression for dt
as a function of R which I can integrate from 0 to R to get t(R).

Does such a thing exist?

Formally, dR/dt = f(R;Omega,lambda) (with H as a scale factor). The
usual approach is then dt = dr/f, which I can integrate from 0 to R to
get t, which is fine because f = f(R). Algebraically, I can also write
dR = f(R)dt, but I can't integrate it from 0 to t to get R, since f is a
function of R, not of t.

Another way to ask the question: I give you arbitrary Omega
(Omega_matter) and lambda (Omega_Lambda) and ask you to make a plot of
R(t). How would you do it?

I'm pretty sure that one can do it with Jacobian elliptic functions
(essentially the inverse functions of Legendre elliptic integrals).
This is more elegant, and faster, but more difficult to implement.
However, I want something simpler: a function which I can integrate
numerically to get R(t) for arbitrary lambda and Omega. Any ideas?