A New Solution for the Friedmann Equations

Assuming a flat universe expanding under a constant pressure and combining the first and the second Friedmann equations, a new equation, describing the evolution of the scale factor, is derived. The equation is a general kinematic equation. It includes all the ingredients composing the universe. An exact closed form solution for this equation is presented. The solution shows remarkable agreement with available observational data for redshifts from a low of z = 0.0152 to as high as z = 8.68. As such, this solution provides an alternative way of describing the expansion of space without involving the controversial dark energy.


Introduction
The evolution of the universe has already been investigated through the first Friedmann equation.As discussed by Carrol [1], the first Friedmann equation is used because it only involves the first derivative of the scale factor.However, the resulting models such as ΛCDM-based models are in terms of limited numbers of parameters representing the ingredients of the universe.Two analytical solutions with restrictive assumptions are already available [2] [3], which will be presented in Section three.There exists no general analytical solution for the Friedmann equations.Here, through combining the first and the second Friedmann equations, a general equation is formulated.Assuming a flat universe, an exact closed form solution for this general equation is obtained.This solution is remarkably consistent with the observational data over a wide range of measured redshifts from a low of z = 0.0152 to the highest recently measured value of z = 8.68.
Except for the flatness assumption, the analytical solution is completely general.It includes all the ingredients forming the universe.The ΛCDM-based models only consider specific combinations of limited numbers of ingredients.Also the value of the cosmological density parameter is analytically estimated to be 0.685568 Λ Ω = .This is essentially identical to the 0.692 0.01 Λ Ω = ± value estimated in 2014 by the Plank Collaboration [4], based on combined data from "Plank + WP + high L + BAO".It is also shown that the contribution of the cosmological constant is to cancel the pressure term in the Friedmann acceleration equation.As a consequence, the expansion equation turns out to be a kinematic equation in terms of the scale factor and its rates of change.
In the next section we develop the new general equation and derive an analytical estimate of the cosmological density parameter.The new analytical solution together with the two existing analytical solutions is presented in Section 3. Comparison of the new analytical solution with the analytical solution involving matter and lamda is presented in Section 4.1.Comparisons of the new analytical solution with the ΛCDM-based models are presented in Section 4.2.In Section 5, the efficacy of the new analytical solution is shown through comparisons with two ΛCDM-based solutions and through comparisons with three sets of observational data.

The New General Equation
The first Friedmann equation, including the curvature, k, and the cosmological constant, Λ, is Alternatively, including the cosmological term in the total density, the above equation can be represented by where now the density, ( ) where m ρ represents the mass density; ( ) ρ represents the energy density; ( ) ρ represents the radiation energy density and ρ Λ is the intrinsic vacuum energy density, which is defined by The curvature can be represented by where 0 T represents the present time and ( ) a t is the scale factor.The Friedmann acceleration equation is given by Substitutions for ( ) 2) and for k from Equation ( 5) back into the above equation yield For the present time, 0 0 1 T H = , the scale factor, ( ) 0 1 a T = .Thus the above equation yields To evaluate Λ from Equation (8), we need to evaluate ( ) ( ) ( ) first.The value of ( ) at the present time is given by ( ) To evaluate the present time values of ( ) 0 R T and ( ) 0 R T , consider the conservation of energy or the alternative form of the first Friedman equation where t M is an equivalent mass.It represents the total mass-energy of the universe including the total mass of all forms of ordinary and non-ordinary masses as well as the total equivalent mass of all forms of energies; including the effect of the cosmological constant.Because of the conservation of the total mass-energy of the un iverse, t M has to be a constant.But ( ) ( ) Because at all points the expansion is taking place in all directions, to evaluate the rate of increase of space between any two galaxies, ( ) R t must be replaced by ( ) Correcting the velocity in the above relation for the effect of time dilation yields ( ) where c represents the speed of light.The value of ( ) 0 R T has been analytically determined [5] to be ( ) Therefore the corrected expansion velocity at the present time is given by Also the present time value of the expansion acceleration, according to Equations ( 13) and ( 14), is given by ( ) ( ) Substitutions for ( ) 0 R T from Equation ( 9), for the corrected velocity ( ) ( ) from Equation ( 15), and for ( ) 0 R T from Equation ( 16) back into Equation (8) yields Solving the above relation for Λ yields Therefore the cosmological density parameter can be represented by Now substitution for Λ from Equation ( 18) back into Equation ( 7) yields the Friedmann acceleration equation as In the next section, assuming a flat universe expanding under the constant pressure 0 P P = , we present an exact closed form solution for the above equation.Before getting to the next section, we will evaluate the value of the cosmological density parameter Λ Ω .The pressure, 0 P , in Equation ( 19) has been given by [5] the following relation where ( ) Assuming a flat universe, i.e., 0 k = , substitutions of the above values into Equation ( 19) yield the cosmological density parameter as This value of the energy density parameter is essentially identical to the 0.692 0.01 Λ Ω = ± value estimated in 2014 by the Plank Collaboration [4], based on data from "Plank + WP + high L + BAO".This remarkable agreement provides further evidence supporting the description of the pressure as given by Equation ( 21).

The New Analytical Solution for the Friedmann Equations
Already there exist two analytical solutions for the Friedmann equations.One is for the case of a flat universe containing only matter, 0 1 m Ω = , where the analytical solution is ( ) ( ) , and the other is for the case of a flat universe containing matter and lambda [2] [3].The analytical solution for the second case is given as ( ) where 0 t T τ = , ( ) a τ is the scale factor and 0 1 m Λ Ω = − Ω .The new analytical solution considers a flat universe containing all the ingredients including matter, energy, radiation, etc. Considering Equation (20), it is seen that as the time flows from the initiation of expansion to the present time, 0 T , the coefficient of the curvature term, k, reduces.Thus as time flows, the universe tends toward becoming flatter.Therefore, according to Equation (20), by the present time, ( ) 0 1 a T = , the universe has become completely flat.Assuming a flat universe eliminates the effects of curvature and simplifies the Friedmann acceleration equation to the following form Assuming the pressure  to be constant implies that in the above equation the term ( ) represents the contribution of the cosmological constant, ( ) being zero is consistent with the vacuum pressure being equal to the negative of the vacuum density as discussed by Carroll [6].Thus the Friedmann acceleration equation is further simplified to the following form: It is clear that the above equation satisfies the present time boundary conditions as given in Equations ( 9), ( 15) and ( 16).To non-dimensionalize the above equation, let where now dot denotes differentiation with respect to τ , and ( ) a τ represents the scale factor.The present time boundary conditions on the above equation, as derived from Equations ( 9) and ( 15), are ( ) ( ) The exact solution of the above nonlinear differential equation with the specified boundary conditions is The beauty of the above analytical solution is the fact that it does not involve the fractional components forming the mix of the universe.Using the Mathematica code [7], a plot of the scale factor as given by the above equation together with its first and second derivatives is presented in Figure 1.
As seen from the above figure, at the present time, 1 τ = , the scale factor, ( ) .Based on observational data [8]- [10], it has been concluded that the expansion initially decelerates but then continues to grow with an accelerating rate.
Considering Equation (18), it is clear that in Equation ( 20), the pressure 0 P represents the contribution of the cosmological constant.It is this pressure, 0 P , that cancels the constant pressure, P, allowing the Friedmann acceleration equation to be simplified to the form given by Equation (28).
In the following subsections the analytical solution given by Equation ( 29) is compared with the analytical solution for a universe containing only matter and lambda as given by Equation (25).It is also compared with other models based on ΛCDM parameterizations.

Comparison of Analytical Solutions
In order to compare the analytical solution given by Equation (29) with the one given by Equation (25), we need to first decide on the value of Λ Ω for substitution in Equation (25).We will use the analytically predicted value as given by Equation (24), i.e., 0.685568 Λ Ω = .The value of the constant C in Equation ( 25) is calculated by equating the present time value of the scale factor to unity.In this way, the value of C is determined to be 0.060274 . A plot of scale factors from Equations ( 25) and ( 29) is presented in Figure 2.  .Equation (29) describes the evolution of the scale factor for a flat universe including all its ingredients.The scale factor calculated from Equation (25), shown as a dashed line, is for a universe composed of matter and Λ only.Its zero time occurs at 0.0485307 τ = . The difference between these two analytical solutions is due to the fact that one of them in-cludes all the ingredients of the universe and the other one considers a universe composed of a specified mix of only two ingredients, matter and Λ.

Comparisons of the New Solution with ΛCDM-Based Models
The ΛCDM-based models characterize the universe with a limited number of energy density parameters as fractions of constituent ingredients.The values of these fractions are estimated through finding the optimum fit to the observationally measured data.The results are presented in terms of distance modulus versus redshift.Here, a comparison of the variation of scale factors versus redshift will be carried out first.The relation between the scale factor and the redshift, z, is defined by To express the scale factor given by Equation ( 29) in terms of redshift, substitution for the scale factor, a, from Equation (29) back into Equation (31) yields the relation between the red shift, z, and the time τ as ( ) ( ) ( ) For ΛCDM-based models, considering Equation (30), and rewriting Friedmann Equation ( 2) in terms of fractions of constituent densities, for a universe containing mass, energy, radiation and curvature, one obtains the ( ) Using the above equation, the values of ( ) are plotted versus z in Figure 4 for two sets of ΛCDM-based parameters, estimated by the Plank Collaboration [4], as presented in Table 1.As seen from Figure 4, the ΛCDM-based curves are consistent with the analytical curve.But they are not identical to the analytical curve.There are two reasons for not being identical.The first reason is the differences in the fractions of ingredients included in the models.The second reason is the fact that the initial times for the scale factors of the ΛCDMbased models are not defined.
In the next subsection the analytical curve and the curves based on ΛCDM parameterization are compared with the observational data.

Comparison of the New Solution with Observational Data
In this part the curve of the analytical scale factor given by Equation (29), transformed through Equation (32), is compared with the curves based on ΛCDM, as given by Equation (33), for the two different sets of ingredients presented in Table 1.To check how well these curves represent the reality, the following three sets of observational data will be used: 1) A set of 557 SNe data with redshifts from a low of 0.0152 z = to a maximum of 1.4 z = , as reported in 2010 in the Union2 Compilation [11]; 2) A set of 394 extragalactic distances to 349 galaxies at cosmological redshifts significantly higher than the Union2 Compilation with redshifts from a low of 0.133 z = to a maximum of 6.6 z = , as reported in 2008 by Mador and Steer [12]; 3) A set of data for a quasar and the three most distant recently confirmed galaxies, as presented in Table 2.
To compare with the aforementioned observational data, the scale factors have to be presented in terms of distance modulus and redshift.The SNe and the Union2 data are already available in terms of distance modulus and redshift.The data for the galaxies and for the quasar are listed in Table 2, and, for the new analytical solution, the data for the scale factors in terms of distance modulus and redshift are represented through the following relation where µ represents distance modulus; a represents the scale factor; 0 is in megaparsecs and the factor ( ) K z represents the effects of observational data such as source luminosity, and data processing correc- tions including the instrument corrections and the K-correction.These are well known corrections and they are considered in various ways [17]- [20].We evaluate the factor ( ) .But most recently (Mostaghel [5], 2015), the value of the Hubble constant is analytically estimated to be Because this value is remarkably consistent with the observationally estimated values, we will use this value and substitute it together with the values of ( ) 1 a − , as obtained through Equations ( 29) and (32), into Equation (34).Through matching the curve of Equation (34) with the first set of observational data presented in Figure 5, the factor ( ) K z is found to be given by  ( ) ( ) As can be seen from )( ) ( ) K z , we include the second and the third sets of observational data.Equation (37) will be used to compare the analytical solution with the ΛCDM-based models and the observational data.The parameters for the ΛCDM curves are given in Table 1.For the analytical solution, as mentioned above, we use

km s Mpc H
. The three sets of observational data, the analytical curve based on Equations ( 29) and (32), and the two ΛCDM curves based on Equation (33), are presented in Figures 6-9.As seen from these figures, in all cases, the analytical curve is remarkably consistent with the observational data as well as with the ΛCDM-based curves.The log-linear plots and the linear plots show how well the curves represent the observational data at the low and high values of the redshifts respectively.
The excellent match of the analytical curve and the ΛCDM curves with the second and the third set of the observational data validates Equation (36) representing the factor ( ) K z .It also confirms the analytically eva- luated value for the Hubble constant as given in Equation (35).It should be noted that, except for the flatness assumption, the analytical solution is completely general.It includes all the ingredients forming the universe.The ΛCDM-based solutions only consider specific combinations of limited numbers of ingredients.

Summary and Remarks
The value of the energy density parameter was analytically estimated to be 0.685568 . This value is essentially identical to the estimated value based on the observational data.This fact and the remarkable consis-  tency of the analytical solution with the observational data, as well as with the ΛCDM-based models, provide the necessary confidence in the fidelity of the analytical solution in the representation of reality.
The pressure is cancelled from the Friedmann acceleration equation through the contribution of the cosmological constant.As the result, Equation (28) may be interpreted as a kinematic equation.Its solution, Equation  (29) describes the evolution of the expansion of space.As such, this equation provides an alternative way of describing the expansion of space without involving the controversial dark energy.
and its rates of change are,

Figure 1 .
Figure 1.Evolution of the scale factor.

Figure 2 . 1 τ = in Figure 2
Figure 2. Evolution of the scale factor with time.The time 1 τ = in Figure 2 represents the present time.For the new analytical solution, shown as a solid line in this figure, the zero time of the scale factor occurs at the time 0.0641186 τ =

Figure 3 . 4 .
Figure 3. Variation of time with the redshift.

Parameters ΛCDM- 1 2 WMAP
Plank + WP + high L + BAO ΛCDM- z through matching of the analytical curve with the first set of the observational data.Then we check the validity of its value through comparisons with the second and third sets of observational data as well as with the ΛCDM-based curves.To evaluate ( ) K z , using the analytical curve, we only need to have the value for the Hubble constant.The recent estimated values of the Hubble constant based on observational data are: (the Seven-Year Wilkinson Microwave Anisotropy Probe [

Figure 5 ,
with this ( ) K z the analytically derived scale factor fits the first set of observational data remarkably well.Substitution for ( ) K z from the above equation back into Equation (35) yields the distance modulus as ( ) ( relation directly gives the distance modulus in terms of the scale factor and the redshift.Now, to check the validity of the factor ( )

Table 1 .
Cosmological parameters used for ΛCDM based models.

Table 2 .
Data for recently confirmed galaxies.