Single scale factor for the universe from the creation of radiation and matter till the present

A scheme for incorporating the creation of radiation and matter into the cosmological evolution is introduced so that it becomes possible to merge the times before and after the creation of radiation and matter in a single scale factor in Robertson-Walker metric. This scheme is illustrated through a toy model that has the prospect of constituting a basis for a realistic model.


I. INTRODUCTION
The question of determining the model that best describes the universe is the ultimate goal of cosmology. The energy-momentum content of the present universe seems to be a perfect fluid mainly consisting of a dark sector (possibly consisting of a dark energy and a dark matter component), baryonic matter, and radiation [1]. In the standard model of cosmology (namely, ΛCDM) dark matter [2] and baryonic matter are considered to be dust, dark energy [3] is taken to be the Einstein's cosmological constant, and radiation is described by the usual energy-momentum term for radiation. Although the standard model seems to be compatible with observations yet it has some problems. The magnitudes of potential theoretical contributions to cosmological constant (CC) are extremely higher than the value of CC deduced from the energy density of the universe [4]. There are many attempts to solve this problem, namely, the CC problem. Nevertheless none is wholly satisfactory. The best option seems to employ a symmetry such as metric reversal symmetry [5] to cancel CC and then attribute the dark energy to something else e.g. to modified gravity [6], or to some scalar field such as quintessence [7]. Cold dark matter (i.e. dust-like dark matter with no or negligible interaction with itself and with baryonic matter and photons) scenario of ΛCDM as well suffers from some problems such as rotation curves of spiral-like galaxies i.e. cuspy halo problem, missing satellite galaxies problem [8,9]. There are many alternatives to cold dark matter (CDM) scenario including warm-dark matter [10], Bose-Einstein condensate dark matter [11,12], and scalar field dark matter [13].
The above considerations essentially hold for the time from the radiation dominated era till the present era. The standard paradigm for the era before the radiation dominated era is an inflationary era (that serves to solve the problems of the standard cosmology such as horizon, flatness, absence of monopoles problems) [14]. Usually the inflationary era and the epoch after this era are studied separately. This is not only due to the need to concentrate on each of these and to try to understand each epoch better before a possible unification. In fact the most serious problem in the direction of the unification 1 of the whole cosmic history is the difficulty of merging these two epochs because of the form of the dependence of the energy density of dust and radiation on scale factor (i.e. on redshift). In ΛCDM the energy density of radiation dominates over that of inflaton if one goes back to sufficiently large redshifts. This is due to the fact that the energy density of inflaton is essentially constant during inflationary era while the energy density of radiation scales like 1 a 4 where a is the scale factor. In other words, to have a true unification, the creation of radiation and matter after the inflationary era must be taken into account in the scale factor without destroying the standard cosmology before and after the inflation, and this is not an easy task. The models in literature that unify all eras of cosmological evolution in a single model [15,16] are not wholly realistic since they do not include baryonic matter although they are able to produce eras of cosmological evolution with correct equations of state in the corresponding eras, and some have graceful exit from inflationary era. The matter in these models must be identified with dark matter since the energy densities of these models do not contain energy components that scale proportional to 1 a 3 for all times (or at least for a sufficiently long time). The models in [15] use the energy densities expressed in terms of simple functions of Hubble parameter and/or scale parameter as the starting point rather than starting from the scale factor. Although one may, in principle, determine the scale factor from this information the form of scale factor may be rather complicated in some cases. On the other hand a relatively simple scale factor may result in a rather complicated and unmanageable functional form for the energy density when expressed in terms of the scale factor or the Hubble parameter.
Therefore in some cases it may be more suitable to consider a specific ansatz for scale factor such as in this study and in [16]. The same approach is adopted in this study. Moreover the present study introduces a general prescription to include dust and radiation into unification.
In this study, first, in Section II, I introduce a scheme to unify the cosmological evolution before and after the radiation dominated era. Then I give a concrete realization of this scheme in Section III. In Section IV I discuss the observational compatibility of this scheme in the context of the model introduced in Section III. Finally I conclude in Section V.
The scale factor in this model is a sum of two terms. The first term is a pure dark energy contribution. The second term is responsible for the baryonic matter and radiation terms and additional terms that may be mainly identified with dark matter. There is also an additional term due to coupling between these terms, and this term gives another contribution to the dark energy and dark matter. Some of the ideas employed here have been already studied in literature. In this study I do not make a sharp distinction between dark energy and dark matter because the dark energy and dark matter terms are coupled and the equation of state (EoS) of some terms e.g. EoS of the coupling term between dark matter and dark energy terms evolve with time. The superficiality of a distinction between dark energy and dark matter is considered in many studies in literature, either explicitly or implicitly [17][18][19][20][21].
This option is quite possible since dark energy and dark matter are not observed directly.
What we see observationally is a missing element in the energy-momentum tensor of the Einstein equations, other than baryonic matter and radiation, and this missing quantity may be described by two components; dark energy and dark matter. It is, in principle, equally possible that this quantity is composed of a single component, say, dark fluid. In [16] I had introduced a universe composed of a dark fluid (that may be written in terms of two scalar fields). In fact the scale factor in that study is essentially a 1 (t) in Eq.(2) of this paper. The present study, somewhat, may be considered as an extension of [16] where baryonic matter and radiation are included. However there are important differences as well. The main aim of this study is to introduce a scheme to merge the cosmological evolution of the time before and after the production of radiation into a single scale factor with the baryonic matter and the usual radiation terms included. The modified form of a 1 (t) in [16] only serves as a realization of this scheme. Furthermore I do not discuss the scalar field identification of the energy density due the part of the scale factor similar to a 1 (t) of [16] (although it can be easily done), and I do not consider the cosmological perturbations of these quantities, and the inflationary era in this study because these points would cause divergence of the main goal of the paper and would increase the volume of this study drastically. I leave these points to future studies.

II. OUTLINE OF THE MODEL
Consider the Robertson-Walker metric I take the 3-dimensional space be flat, i.e.g ij = δ ij for the sake of simplicity, which is an assumption consistent with cosmological observations [22,23]. I let the form of the scale factor be where t 0 denotes the present time. We will see that a 1 (t) is the part of the scale factor responsible for dark energy and dark matter, and a 2 (t) is the one mainly responsible for dust and radiation and additional contribution to dark matter-energy, and we shall see later that a mixing between the sectors due to a 1 and a 2 act as an additional source of dark energy. We assume that a 1 (t) and a 2 (t) are chosen in such a way that a(t) > 0 for all t.
In general one may identify the dust by a mixture of baryonic matter and dust-like dark matter. The best fit values that I could find by trial and error for the specific toy model considered in this study for implementation of the present scheme seem to prefer the case where the dust term is wholly or almost wholly due to baryonic matter.
We first focus on the a 2 (t) term and specify it as where x(t) is some function that its form will be specified later. Eqs. (2) and (3) may be used to relate a(t) and a 1 (t), a 2 (t) in a more applicable way, and to derive the corresponding Hubble parameter. We observe that In a similar way the Hubble parameter is found to bė where we have usedȧ Note that a(t 0 ) = 1 by convention.
We let x 1 (t) = exp t t 0H (1) 2 dt (8) where α o1 , α o2 , α b , α r , α x , α K are some other constant coefficients. In fact, in (7) we could take the simpler form where α c = 0, α o1 = α o2 = 1, c 1 = 1, c 2 = 0. This would be enough as long as we are concerned only with merging of the eras before and after the radiation domination, and the resulting model would be compatible with Union2 data set at an order of magnitude level. The more involved form in (7) is used to make the model phenomenologically more viable. This point will be discussed when we discuss the phenomenological viability of the model in Section IV. One may determineẋ in Eq.(5) by using Eq. (7),ẋ where Hence one may express (5) as where We let where . Because the three dimensional part of metric is taken to be flat the present energy density is equal to the critical energy density the above equations imply that Note that, at this pointΩ 1 ,Ω b ,Ω r ,Ω x ,Ω K cannot be identified as density parameters since density parameters should satisfy Ω 1 +Ω b +Ω r +Ω x +Ω K =1. In Chapter IV we will see that this condition is not satisfied for the phenomenologically viable sets of parameters, soΩ 1 ,Ω x , Ω K cannot be identified as density parameters separately, instead one must define the total density parameter for dark sector by Ω 1 2 D =Ω 1 −Ω x −Ω K rather than the separate contribution due to H 1n and H ∆ while I identifyΩ b ,Ω r as the density parameters corresponding dust and radiation. Therefore to retain the physical content of this paper more evident I will not make a distinction betweenΩ b ,Ω r and the density parameters for baryonic matter, radiation; Ω b , Ω r while I keep this distinction for the others i.e. for the ones due to H 1n and H ∆ terms.
and αr a 2 terms result in energy densities that are identified as the energy densities for baryonic matter and radiation. In principle, there may be also contributions due to the Ξ 1 a 3 and ψ 1 a . The sign of the Ξ 1 a 3 term is negative of the usual stiff matter. It may be identified as stiff matter under pressure so that it has a negative deceleration parameter. The main function of this term is to damp the energy densities of baryonic matter and radiation in the time before the radiation dominated era. The function of the 1 a term is similar. It ensures the behavior of the energy density in late times be well-behaved (i.e. preventing the energy density to grow too fast (through the 1 a term in x 2 (t) and x 3 (t))). Although the ψ 1 a term is similar to that of a negative curvature 3-space it is different from such a term since its origin is the Hubble parameter H while a usual 3-curvature term arises from the 3-curvature part of metric. Note that this term arises even in a flat 3-space in this construction. Therefore I identify the Ξ a 3 and ψ a terms in H as additional contributions to dark sector. Another point worth to mention is; It is evident that the square of (5) (in conjunction with (10)) results in an A 2H 2 2 term containing A 2 α 2 b a 3 and A 2 α 2 r a 4 terms that may be identified with the standard baryonic matter and radiation terms, respectively if A is taken to be constant while it depends on time in this scheme as it is evident from (7). In fact variation of A with time makes it possible to go to zero before the radiation dominated era as desired. Therefore, given the considerable success of the standard model at least in the observed relatively small redshifts, the variation in A after the matter -radiation decoupling time should be small so that this scheme mimics the standard model at relatively small redshifts where observational data is available. If one takes dA dt t≃ t 0 sufficiently small one may guarantee an almost constant value for A for a sufficiently long time (e.g. from the present time till the beginning of the radiation dominated era). We will see in Section IV that there exist such values of A with reasonable phenomenological viability. Another term arising fromH 2 2 is the cross term, A 2 α b αr a 7 2 . This term may be identified as the energy density term due to the transitory time where massive particles that act as radiation at high energies turn into more dust-like at intermediate energies. Another term in H 2 is H 2 1n . This term will be considered as a pure dark sector term. Finally the cross term 2H 1nH2 gives an additional contribution to the dark sector for the phenomenologically viable values of the parameters. It may be easily shown that this term does not necessarily imply strong interaction between the dark fluid and radiation and baryonic matter as its form may suggest if the parameters of the underlying physics at microscopic scale satisfy some restrictions. Otherwise one may use screening mechanisms such as [24][25][26] to explain the unobservablity of dark matter-energy.
Next we derive the general form of the equation of state for this model. We derive the explicit form of the equation of state after (EOS) after we give the explicit form of a 1 (t) in the section. However giving the general form of EOS in this scheme provides us a more model independent formula and may be useful for other choices of a 1 (t) in future. After using Eqs. (15,16,17,18) one obtains EOS, ω as The terms inside the first parenthesis in the second line correspond to the contribution of the dark sector term H 1n . The other terms in the same line correspond to the contributions of dust and radiation and their coupling with dark sector term H 1n . The remaining terms are the term corresponding to variation of A, the term corresponding to coupling of curvaturelike term and the stiff matter under negative pressure with dust and radiation, the term corresponding to coupling of curvature-like term and the stiff matter under negative pressure with H 1n , the term corresponding to coupling of curvature-like term and the stiff matter under negative pressure with the other terms, and the contribution of the curvature-like term and the stiff matter under negative pressure, respectively. It is evident from (21) that the pressure for baryonic matter is zero as should be, and the pressure for radiation is 1 3 as expected. A point worth to mention at this point is; The coupling term between baryonic matter and radiation in Eq. (21) has an equation of state 1 6 (that may be seen by considering the ratio of the α b αr a 7 2 in p by the corresponding term in ρ i.e. 2 α b αr a 7 2 ). The redshift dependence of this term is between that of baryonic matter and radiation. This time dependence is more natural than the standard picture where there is no such term. Massive particles at high energies act as radiation and at lower turns into dust. The coupling term accounts for the transitory time when massive particles pass from radiation to dust state.
In order to obtain the evolution of ω as a function redshift or time explicitly, H 1n must be specified. This will be done in the next section. However I give a ω versus redshift graph in Fig.1 for a 1n introduced in the next section for a phenomenologically viable set of parameters (i.e. those with small χ 2 values and with energy densities for recombination and nucleosynthesis as discussed in Section IV) to have an idea about the evolution of ω with redshift. To draw this graph I have converted time, t to redshift, z (for Union2.1 data) through the relation z = 1 a − 1, and then used Mathematica to use this relation to make the calculations (although the original quantities are expressed in terms of time). This procedure is applicable for small redshifts. However, in general, it becomes inapplicable due to highly non-linear form of scale factor and Hubble parameter since it requires huge RAM and CPU for computation, if it can be done at all, and hence requires a separate computational physics project by itself. Therefore I have used equation of state versus and energy density versus time graphs (instead of redshift) in Section IV. In fact, even that option required a long time of order of months to make the necessary computations.

III. AN EXPLICIT REALIZATION OF THE MODEL
Now we focus on the a 1 (t) term. We take where A 1 < 1, p 1 , p 2 , b 2 , b 1 are some constants that to be fixed or bounded by consistency arguments or cosmological observations. This scale factor is a generalization of the scale factor in [16] where r = 1, s = 6. A similar scale factor is considered in [27] as well. One of the shortcomings of [16] is that the present value of the equation of state parameter in that model (for phenomenologically relevant choices of parameters where the model mimics [22,23]. In the present study there is an additional contribution due to mixing of the terms due to a 1 and a 2 and hence there is less need to modify the scale factor in [16]. However I prefer to adopt the more general form in (23) to seek a greater parameter space and to insure the correct equation of state parameter.
We had shown in Eq. (5) We let and where t 0 is the present age of the universe. One observes from (19) and the above expression We will see in the next section thatΩ 1 cannot be identified as the density parameter corresponding to H 1n . Instead one must define an overall density parameter for the dark sector by Ω Observational values of H 0 = ȧ a t=t 0 and 1 t 0 are almost the same. Therefore ξ 2 ξ 2 1 ≃ ξ 2 1 . After determining the H 1n we are almost ready to find the explicit values of the energy density and the equation of state. The only missing element for calculation of these quantities is to find A, Ξ, ψ in (15,16,17,18). Another point to be addressed is to show that there exist sets of A whose variation with time are small for low redshifts so that the terms that are proportional to 1 a 3 2 and 1 a 2 in AH 2 term may be identified with dust and usual radiation terms, respectively.
In order to determine A, Ξ, ψ (and to determine the rate of variation of A with time) one should derive an approximation scheme for the evaluation of these quantities because these quantities depend on x 1 (t), x 2 (t), x 3 (t) (that are defined in (8) and (9)), and these quantities, in turn, are defined in a recursive way since , and a(t), in turn, depends on x i (t) through Eq. (4). In other words, in order to determine the approximate values of x i (t) one must identify the zeroth order approximation and a method how to obtain the higher order approximations in an iterative way. One may use the following observations to obtain the zeroth order approximation; This implies that the zeroth order approximation for the scale factor a(t) should be taken as a (0) (t) = a 1n (t) Hence for phenomenologically viable cases (whereȦ ∼ 0 for small redshifts) one may take the zeroth order approximations as whereH (10)  2 ,H 2 ,H Then One may get the next order approximation by using The next order quantities A (1) , x (1) may be obtained from (32) and (33) by replacing the superindices (0) by (1) where In principle this may be done up to arbitrarily higher order approximations but it is quite difficult to calculate even A (1) even with the help of computers. In fact I have divided the interval t − t 0 in to coarser subintervals to decrease the CPU time and have used the approximate numerical values in the i'th interval (by assuming A (0) to be almost constant in those intervals) by using the formula to find A (1) . I have seen (by trial and error) that it is possible to find almost constant A (0) and A (1) values for many relevant (i.e. of small χ 2 values considered in the next section) For example the variations of A (0) and A (1) with time for one of the phenemonologically viable sets in Table III is given in Table I. 2-In searching for the phenomenologically viable parameter space I start from the values of the parameters in [16] i.e. r = 1, s = 6, ξ = 1, and β ∼ O(1) since the universe studied in [16] mimics the true universe roughly.
3-Due to the highly non-linear relation between the Hubble parameter and the scale factor I seek the relevant parameter space usually by trial and error rather than a continuous scan of the parameter space. Therefore the optimum values obtained here most probably may not correspond to the best possible optimization. Rather they hopefully correspond to a good approximation to the best optimal values.
A. Compatibility with Union2.1 Data In this subsection we use the Union2.1 compilation data set to find the optimal values of β, r, s starting from β = 3, r = 1, s = 6. We find the theoretical values of distance moduli, µ for the redshift values of Union2.1 and calculate the corresponding χ 2 value by using the measured values of µ and their errors.
The expression for distance modulus is where for small redshifts reduces to where we have used the requirement that A 1 c 0 1+c 0 x ≃ 1 at small redshifts as discussed in the preceding section (see Table II), and a 0 = a(0) = 1. In ΛCDM dt a(t) is usually expressed in terms of redshift, z and Hubble parameter H, and then the results for different z's are compared with the data directly. This is not possible in this model because H cannot be expressed in terms of a(t) in a simple way. Therefore in this study first we convert redshift values of Union2 to time values by using z = 1 a(γ) − 1 ≃ 1 a 1n (γ) − 1 and then solve it for γ. The corresponding expression for the theoretical value of the luminosity distance d L in this case (i.e. in terms of γ) is where a 1n (t) is expressed in terms of β, r, s, γ = t t 0 by using the parameterization given in the preceding section. Eq.(40) may be written in more standard form in terms of H 0 by using H 0 t 0 = ξ. Then we find (39) numerically for each of the γ corresponding to observational redshifts. Finally we find the corresponding χ 2 0 values by using the formula where the subscript 0 in χ 0 and the superscript (0) in µ th(0) stands for the fact that a(t) is approximated by its zero'th order approximation i.e. by a 1n , the superindices th and obs stand for the theoretical and observational values of µ, and the subindices i denote the values of the corresponding quantity for the i'th data point in Union2 data set.
One may try a better approximation by replacing a 1n (t) in (39) by a better approximation of a(t) i.e. by c 0 A 1 1+c 0 x (0) (t) a 1n (t) where x (0) (t) is defined by Eq.(31). In principle, then, one may evaluate the integral (38) after replacing a 1n (t) by c 0 A 1 1+c 0 x (0) (t) a 1n (t). However this seems to be inapplicable for standard computers because of the complicated form of the integral. One needs a separate computational physics project for this aim. Instead one may try a rough approximation (hopefully better than a 1n ); we take the 1+c 0 x (0) c 0 A 1 term in the integral to outside of the integral with its γ value being the bound of the integral. This approximation is a good approximation provided that c 0 A 1 1+c 0 x(t) does not vary much in the time interval between t 0 and the time corresponding to the given redshift value. Otherwise the higher order approximation may worsen the approximation rather than improving. The corresponding formulas (in the first order approximation) become After trial and error I have found many sets of parameters with relatively small χ 2 0 , χ 2 values. For example the χ 2 0 , χ 2 values for two phenemonologically viable sets of parameters are given in Table III where the reduced χ 2 0 , χ 2 red 0 = The sets of parameters (that I could by trial and error) with relatively small χ 2 values satisfy c 1 ≃ c 2 ≃ 1, α c ≪ 1. By using this information one may check the validity of (20) and determine if one may identifyΩ 1 ,Ω,Ω K by the corresponding density parameters; Ω 1 , Ω x , Ω K for the phenomenologically relevant parameters by using Eq. (19) and Eq. (17) and Eq. (18). We observe that x(0)=A 2 , x 1 (0)=x 2 (0)=x 3 (0)=1, c 1 = 1 and for relevant values of the parameters Hence, after using (19), we obtaiñ Ω 1 2 x We observe that for phenemonologically viable sets of parameters, for example, for those in x +Ω 1 2 K ) 2 = 1. However one may define a total density parameter for the dark sector by Then the density parameters satisfies the necessary condition, (Ω r ) 2 = 1. In other words, H 1n and H ∆ terms can not be identified as separate contributions to dark sector, rather they must be considered as just a single object in order not to introduce an ambiguity in their identification.

B. Compatibility with Recombination and Nucleosynthesis
In this subsection we investigate if this model is compatible with the cosmological depiction of the recombination and nucleosynthesis, at least, at the order of magnitude level.
In a similar vein as the preceding subsection we require this model mimic the standard model, ΛCDM, as much as possible. We assume that the radiation and the baryonic matter are in thermal equilibrium in the eras of recombination and nucleosynthesis since we adopt the same equations of thermal equilibrium as ΛCDM. Therefore, in the following, first we drive the condition for thermal equilibrium for this model. Then we find the sets of parameters with least χ 2 values that may produce successful recombination and nucleosynthesis eras. The correct choices should have sufficient radiation energy densities in these eras. In other words the redshift at the recombination time, z re should be in the order of (1 + z re ) 4 > (1 + z * ) 4 ≃ (1100) 4 where * denotes time of last scattering surface; and in the nucleosynthesis era the energy density of neutrinos should reach energy densities of the order of (1 MeV ) 4 . We seek an approximate, rough agreement with ΛCDM since the search of the parameter space is done by trial and error rather than a systematic search of the whole parameter space. Therefore a detailed, thorough analysis and compatibility survey would be too ambitious especially considering this is a toy model. Before checking if there exist a set of parameters compatible with recombination and nucleosynthesis we should check if the thermal equilibrium is maintained in these eras in for the given set of parameters because we adopt the standard analysis in ΛCDM, and that analysis assumes existence of thermal equilibrium. As is well known, if there is thermal equilibrium then we should have Γ > H where Γ is the rate of the interaction between radiation and the matter and H is the Hubble parameter. However the implementation of this condition in this model is not exactly the same as in ΛCDM. In the case of recombination era the implementation of this condition does not give exactly the same result as ΛCDM since , in ΛCDM the recombination takes place in radiation dominated era and the total energy density is almost wholly due to radiation while, in this model, the total energy density of the universe at this era is not almost wholly due to radiation although the equation of state parameter for phenomenologically relevant cases is similar that of radiation dominated universe at the time of recombination and we require the radiation energy density to be the same or almost the same as ΛCDM. In the case of nucleosynthesis, even the equation state parameter in this model does not mimic that of a radiation dominated universe. Therefore we should derive the corresponding conditions for thermal equilibrium for this model.
The condition for thermal equilibrium in the recombination era is Here we have used the identities, where α 1 ≃ 1 is the ΛCDM value and α 2 1 ≤ 1 at the time of recombination and is not constant in this model, and I have used the PDG values, H 0 = 72 km Mpc −1 s −1 , Ω ph = 4.8 × 10 −5 . Note that Γ in (46) is the same as the ΛCDM value while H is different from the ΛCDM value.
Next consider the condition on thermal equilibrium at and before the time of nucleosynthesis. In thermal equilibrium we have where we have used the identity similar to (47), where α 1 and the subindex ph are replaced by α 2 and r, respectively and the ratio is evaluated at the time of nucleosynthesis. In this case, as well, Γ ν is the same as its ΛCDM value while the expression for H in terms of temperature is different since α 2 ≤ 1 and is not a constant (i.e. it gives a different value when evaluated at different time during nucleosynthesis) in this model while α 2 = 1 in ΛCDM. During thermal equilibrium the ratio of neutrinos to all nucleons, X n is given by where Q is the rest mass energy difference between a neutron and a proton, Q = m n − m p = 1.239 MeV . After the thermal equilibrium between the neutrinos and the nucleons are lost i.e. after decoupling the value of X n further decreases due to decay of free neutrons as where X n0 is the X n of Eq.(49) at the time of decoupling, and τ 0 = 885.7 seconds is the lifetime of a free neutron. Therefore the effect of this model is to change the value of X n0 (that depends on α 2 ) and probably the value of X n as well.
Now we are ready to check the viability of this model. I could give only four graphs and three tables that partially summarize the results of my calculations related to this and the next paragraphs in order not to expand the size of the paper too much. Otherwise the size of the manuscript would be almost doubled. First we check the viability of the model for recombination and nucleosynthesis eras. To this end I have used the equations (21,28,17,18) in the zeroth order approximation where a(t) ≃ a 1n (t) (as discussed before Eq.(30)) to draw ω, ρr ρ r0 , versus time graphs by using a Mathematica code that I have prepared for this aim for the sets of the parameters, r,s,β,ξξ 1 , , that correspond to some relatively small χ 2 values obtained in preceding subsection. Then I have tried to find at least one set of parameters with phenomenologically viable ω 0 , ρr ρ r0 , ρ ρ 0 values i.e. ω 0 , in the range −0.68 -− 0.74; ρr ρ r0 > (1100) 2 ≃ 10 12 (in the range of redshifts z ∼ 800 − 3000) at the time of recombination, and ρr ρ r0 > (1 M eV ) 4 5× 10 −5 (2.5× 10 −3 eV ) 4 > 10 38 at the time of nucleosynthesis where I have approximated 3 (t) (that are defined in Eq.(30)), and a(t) by a 1n (t) (that is defined in (23)) as discussed in the preceding section. I have found two sets of parameters given in Table III that satisfy these conditions. A comment is in order at this point. The zeroth order approximation is reliable only for small redshifts. However this approximation is reliable at any redshift if one is only interested in the energy density -redshift relation. This may be seen as follows: Assume that the energy density ρ is related to redshift z by ρ = f (z) in the zeroth order approximation (where f (z) is an arbitrary function), and in an approximation better than the zeroth order . Then the energy density after the correction is ρ ′ = f (z ′ ). If one rescales z ′ as 1 x z ′ = z then one obtains the same redshift and energy density values. In other words the redshift -energy density relation is invariant under such corrections. However this is not true for the redshift -time relation. If the approximation is not a good approximation to the true value then the redshift -time relation will be distorted. This, in turn, may cause the distortion of the value of the equation of state and the distortion of the variation of the energy densities with time in an amount depending on the reliability of zeroth order approximation. Keeping these observations in mind I content to use zeroth order approximation for the times of recombination and nucleosynthesis because even employing zeroth order approximation needs a lot of computer CPU and RAM, and in many cases the use of first order approximation as well does not improve the situation. We will come back to these points when discuss the times of recombination and nucleosynthesis.
Next I have checked if thermal equilibrium is maintained at the times of recombination and nucleosynthesis and if recombination and nucleosynthesis are realized in this model.
One may have an idea on thermal equilibrium at the time of recombination by using the values of Table III at z ≃ 1100 and Eq.(46). However a more rigorous way is to draw (that may be obtained from Eq.(46)) versus time graphs to determine the time intervals (and then the corresponding redshift intervals) where ρ ρ 0 ρr ρ r0 for each of the sets A and B. In fact I have used T T ph0 = 1 + z for the relevant redshifts. The resulting intervals are the intervals where thermal equilibrium is maintained as shown in Figure 4 for the set B in Table III. The smallest redshifts where the thermal equilibrium is lost are z = 2317 (γ = 1.615× 10 −10 with ρr ρ r0 ≃ 5.52× 10 13 ) and z = 1625 (γ = 2.82× 10 −10 ρr ρ r0 ≃ 3 × 10 13 ) for the sets A and B, respectively. This implies that the photon electron decoupling takes place before the time of last scattering at an energy of ≃ 2317 × 6 × 10 −4 A detailed comprehensive separate study is need to see if these imply some interesting phenomenologically viable alternatives or just an artifact of the toy model and/or the sets of parameters considered. This may also be due to the limitation of the applicability of zeroth order approximation that we have discussed above. a(t) ≃ a 1n (t) is not violated badly at the time of recombination for the most of the relevant sets of parameters. For example for the sets of parameters given in Table III the first order approximation results in a(t) ≃ 0.4 a 1n (t) i.e. c 0 A 1 A 1 −A 2 ≃ 0.4 and does not vary much at the time of recombination. Therefore it seems that the effect of the limitation of the applicability of zeroth order approximation to the time of recombination must be limited. However this shift does not introduce a major problem since the redshift values, hence the photon energy density at recombination remains almost the same and thermal equilibrium is maintained. Next I have checked if thermal equilibrium is maintained at the peaks in Table III where the energy densities are sufficient for nucleosynthesis. I have used Eq.(48) to find the range of temperatures where thermal equilibrium is maintained. I have found that this condition is satisfied for T > 3 × 10 10 K (provided that Ω r ≃ 5 × 10 −5 ) for the second peaks. This value gives us X n0 in (50) by using (49) as X n0 ≃ 0.39 which is quite large compared to the ΛCDM value of ≃ 0.25.
The time that takes 3 × 10 10 K ≃ 1 MeV drop to 0.07 MeV (that is when ρr ρ r0 ∼ 10 32 ) in this model is something like ∼ 2 × 10 −16 × t 0 ≃ 90 seconds. Therefore X n0 does not drop significantly through Eq.(50). In other words the final result X n ≃ 0.35 is much larger than the ΛCDM value ≃ 0.13 (which agrees well with observations). Probably the main source of this discrepancy is inapplicability of zeroth order approximations to redhifts and energy densities to this era to obtain correct energy density -time relations. The variations of c 0 A 1 A 1 −A 2 and A are quite large and their values are quite different than those at z ∼ 0 at the time of nucleosynthesis that makes the applicability of the zeroth order approximation extremely difficult to obtain correct energy density -time relation. In other words the main source of the discrepancy may be due to the fact that the real time of free decay may be of the order of ∼ 1000 seconds in this model instead of 90 seconds. The use of first order approximation does not improve the situation because in the calculation of first order approximations to c 0 A 1 A 1 −A 2 and A(t) one uses the zeroth order approximation a(t) ≃ a 1n (t) in the integrals for x i , i = 1, 2, 3. This, in turn, results in over contribution of large redshifts and hence larger and more varying x i with respect to their true values since x i < 1 and they get smaller i.e. c 0 A 1 A 1 −A 2 gets larger at larger redshifts. Therefore the energy density versus time graphs in the figures 2, 3 and 4 must be considered with some care: The time values in those graphs should be taken with utmost care especially in the case of nucleosynthesis while the magnitudes of energy densities and the corresponding redshifts are expected to be the same as the exact values. All these points must be studied in more detail in future studies. However I have been able to show that this scheme can produce a model that mimics the standard model: There is a current accelerated epoch whose present equation of state (for the whole universe) is -0.7 (that is, at least, roughly in agreement with observations e.g. see the value in Table   III for a phenomenologically relevant set of parameters). Before this epoch ω changes sign and the time near this sign change may be considered as the matter dominated era. Although the sign change of ω occurs at a later time in this model compared to ΛCDM the time and the redshift of onset of the accelerated era (i.e. ω ≃ − 1 3 ) are comparable with those of ΛCDM. There is an epoch before the matter dominated era where ω is on average close to 1 3 , and may be identified by radiation dominated era, and the time of the maximum value of ω may be considered as the time when the universe was like stiff matter or denser (as in the cores of stars). Then ω changes sign again reaches to two minima peaks as mentioned before and eventually approaches to -1 as time goes to zero (due to the H 1n , in particular the first part of it) and this epoch probably may be considered as the inflationary era. Moreover the model is able to give relatively small reduced χ 2 0 and χ 2 values for Union2.1 data set, and it can, at least roughly, account for recombination and nucleosynthesis times. I think this is a sufficiently well starting point for a toy model whose main aim is to embody the creation of matter and radiation in the scale factor of Robertson-Walker metric. However there is a great deal of points to be clarified and addressed in future studies such as checking the whole parameter space of this model by using a more elaborate software and to use more powerful computers that may give scan the whole parameter space in a better approximation than the one given here, and considering a more detailed analysis of the recombination and nucleosynthesis epochs, studying the evolution of cosmological perturbations in this model, and considering possible extensions of this model towards a more realistic model.

V. CONCLUSION
In this study a scheme for obtaining a scale factor (in Robertson-Walker metric) that may account for the times before, during, and after the radiation dominated eras is introduced.
The prescription to obtain the scale factor in this model is quite simple; First one introduces a scale factor for the pure dark sector, and then the full scale factor is obtained by a relation between these two scale factors. The result is a scheme to produce the scale factor for the whole universe, including baryonic matter, radiation, and dark energy-matter (i.e. dark sector) in such a way that the times before, during, and after radiation dominated era are expressed by a single scale factor in Robertson-Walker metric. Different choices of the pure dark sector scale factor (denoted by a 1 in this paper) and different choices of the relation between a 1 (t) and the scale factor of the full universe, a(t) give different models. As an illustration of this scheme a model with a specific scale factor for the pure dark sector and a specific relation between a 1 (t) and a(t) is considered. The phenomenological viability of this model is checked through its compatibility with Union2.1 data set, and with recombination and nucleosynthesis by using trial error and Mathematica software for almost randomly chosen sets of parameters. Two sets of parameters with relatively small χ 2 values for Union2.1 data set, and that are compatible with successful recombination and nucleosynthesis at an order of magnitude level are found. These results are encouraging in view of the fact that only a tiny portion of the whole parameter space could be considered in this way. A separate, detailed, and comprehensive computational project with more advanced software codes and/or powerful computing facilities that may scan the full parameter space and may employ better approximation schemes is needed to reach a definite view on the observational viability of this scheme and/or this model. Moreover the effect of this model on cosmological perturbations should be considered and possible implications and extensions of this scheme to inflationary era should be studied in future. Furthermore different pure dark sector scale factors and different options to relate pure dark sector and the full universe scale factors may be considered in future to see the full range of possibilities that this scheme may offer.