A. Neumaier said:
I am looking for reliable information about the functional dependence of the diameter ##d(t)## of the visible universe on the time ##t## since the big bang singularity, based on the different hypotheses currently deemed competitive.
One caveat: few cosmologists use the diameter of the universe. Sometimes the radius of the observable universe is used. There's no measure for the size of the universe beyond the observable patch.
That said, there is no closed-form expression. It has to be estimated numerically.
The equation of interest for determining this is the first Friedmann equation:
##H^2 = {8\pi G \over 3}\rho - {kc^2 \over a^2}##
Where:
##H = {1 \over a}{da \over dt}##
Here ##a(t)## is the scale factor, usually defined so that ##a = 1## at the current time. If the diameter of the universe at the current time is ##d_0##, then the diameter of the universe at any other time will be ##d_0 a(t)##. ##\rho## is the energy density of the universe divided by ##c^2##. ##k## is the spatial curvature.
What is usually done is to make use of stress-energy conservation to show how the energy density of each component of the universe changes over time. For example, the energy density of matter scales as ##1/a^3##, the energy density of radiation scales as ##1/a^4##, and the energy density from dark energy is constant (in the simplest case). We then parameterize the equation in terms of the current density fraction, like so:
##H^2 = H_0^2\left({\Omega_r \over a^4} + {\Omega_m \over a^3} + {\Omega_k \over a^2} + \Omega_\Lambda \right)##
Here each ##\Omega## is a dimensionless number, and ##H_0## is the current expansion rate. Since we have defined ##a = 1## at the current time, when ##a = 1##, ##H^2 = H_0^2##. Thus for the above equation to be valid, ##\Omega_r + \Omega_m + \Omega_k + \Omega_\Lambda = 1##. This is why the ##\Omega## terms are called the density fractions for each component. These four density fractions and the current Hubble expansion rate ##H_0## must be measured experimentally. Once you have values for the five parameters, it's possible to use a differential equation solver to get ##a(t)##, which you can use to get the diameter as a function of time.