Thanks for the comments. Another variable they deal with is efficiency---the fraction of gravitational energy that is converted to radiation or other output---and I think their model assumes efficiency increase over time. The model starts at some time far back in the past---large z---and extends up to the present---namely z = 0.
How could they go further than that? There is no negative z