This study considers the specific case of a flat minimally coupled to gravity quintessence cosmology with a dark energy quartic polynomial potential that has the same mathematical form as the Higgs potential. Previous work on this case determined that the scalar field is given by a simple expression of the Lambert W function in terms of the easily observable scale factor. This expression provides analytic equations for the evolution of cosmological dark energy parameters as a function of the scale factor for all points on the Lambert W function principal branch. The Lambert W function is zero at a scale factor of zero that marks the big bang. The evolutionary equations beyond the big bang describe a canonical universe that is similar to inlineformulammlmath idmm1mmlsemanticsmmlmommlmommlsemanticsmmlmathinlineformulaCDM making it an excellent dynamical template to compare with observational data. The portion of the W function principal before the big bang extends to the infinite prebang past. It describes a noncanonical universe with an initially very low mass density that contracts by rolling down the dark energy potential to a singularity big bang at the scale factor zero point. This provides a natural origin for the big bang. It also raises the possibility that the universe existed before the big bang and is far older and that it was once far larger than its current size. The recent increasing interest in the possibility of a dynamical universe instead of inlineformulammlmath idmm2mmlsemanticsmmlmommlmommlsemanticsmmlmathinlineformulaCDM makes the exploration of the nature of such universes particularly relevant.
['cosmological constraint', 'dark energy', 'theoretical model', 'General Relativity and Quantum Cosmology', 'Astrophysics - Cosmology and Nongalactic Astrophysics']
NonCanonical Dark Energy Parameter Evolution in a Canonical Quintessence Cosmology
{'No Header': 'Citation: Thompson,R.I.Special Issue for the Journal Universe "[Dark Energy andDarkMatter]". Journal Not Specified 2024 , 1 , 0. https://doi.org/  \nReceived: Aug. 12 2024  \nAccepted: Sept. 2, 2024  \nPublished: Sept. 5. 2024  \n<!-- image -->  \nArticle', 'Non-Canonical Dark Energy Parameter Evolution in a Canonical Quintessence Cosmology': "Rodger I. Thompson 1,†,‡ *  \nPublisher's Note: MDPIstaysneutral with regard to jurisdictional claims in published maps and institutional affiliations.  \nCopyright: © 2024 by the author. Submitted to Journal Not Specified for possible open access publication under the terms and conditions of the Creative Commons Attribution (CC BY) license (https:// creativecommons.org/licenses/by/ 4.0/).  \n- 1 Department of Astronomy and Steward Observatory University of Arizona; Tucson, AZ. 85721,USA; rit@arizona.edu  \nAbstract: This study considers the specific case of a flat, minimally coupled to gravity, quintessence cosmology with a dark energy quartic polynomial potential that has the same mathematical form as the Higgs potential. Previous work on this case determined that the scalar field is given by a simple expression of the Lambert W function in terms of the easily observable scale factor. This expression provides analytic equations for the evolution of cosmological dark energy parameters as a function of the scale factor for all points on the Lambert W function principal branch. The Lambert W function is zero at a scale factor of zero that marks the big bang. The evolutionary equations beyond the big bang describe a canonical universe that is similar to Λ CDMmaking it an excellent dynamical template to compare with observational data. The portion of the W function principal before the big bang extends to the infinite pre bang past. It describes a noncanonical universe with an initially very low mass density that is contracting by rolling down the dark energy potential to a singularity, big bang, at the scale factor zero point. This provides a natural origin of the big bang. It also raises the possibility that the universe existed before the big bang and is far older and was once far larger than its current size. The recent increasing interest in the possibility of a dynamical universe instead of Λ CDMmakes the exploration of the nature of such universes particularly relevant.  \nKeywords: cosmological constraint; dark energy; theoretical model", '1. Introduction': 'Flat, minimally coupled to gravity quintessence was one of the first and is probably one of the simplest rolling scalar field cosmologies. Even so quintessence has a complex range of properties that depend on the mathematical form of the dark energy potential and its arguments, as well as the boundary conditions on the cosmological parameters such as the current dark energy equation of state w 0 . The dark energy potential in this study is a quartic polynomial with the same mathematical form as the Higgs potential that is hence referred to as the Higgs Inspired, HI potential. In a Ratra Peebles format [1,2] the HI potential is  \nV ( κθ ) = M 4 (( κθ ) 2 -( κδ ) 2 ) 2 = M 4 (( κθ ) 4 -2 ( κθ ) 2 ( κδ ) 2 +( κδ ) 4 ) (1)  \nwhere M is a constant with units of mass in terms of the reduced Planck mass Mp , and κ is the inverse reduced mass 1 Mp The scalar field is θ and δ is a constant, both with units of mass, which makes κθ and κδ dimensionless. All of the HI potential dimensionality is in the M 4 term. Previous work [3] on this combination of cosmology and potential determined that the scalar field θ ( a ) has a simple solution in terms of the scale factor a and the Lambert Wfunction [4].  \nκθ ( a ) = κδ √ -W ( χ ( a )) . (2)  \nThe argument χ ( a ) of the W function is a power law function of the scale factor a .  \nχ ( a ) = qa p . (3)  \nFigure 1. The upper panel shows portions of the principal, solid line, and minus, dashed line, branches of the Lambert W function W ( χ ) . The big bang lies at the 0,0 origin. Evolution is from right to left as indicated by the arrow in the figure. The dotted vertical line in the pre bang region marks the pre bang W ( χ ( a )) = 1 point. The slightly thicker portion of the χ = 0 axis indicates the width of the post bang universe. The lower panel shows the post bang universe. The extent of the solid line shows the present universe for w 0 = -0.999 and the extent of the dashed line for w 0 = -0.99. The terminator of each line is at a scale factor of one.  \n<!-- image -->  \nBoth q and p are constants determined by the boundary conditions and the HI dark energy potential. The expression ( κθ ) 2 -( κδ ) 2 appears in the HI potential and several other places in the evolutionary functions. In terms of the W function it is.  \n( κθ ) 2 -( κδ ) 2 = -( κδ ) 2 ( W ( χ ( a )) + 1 ) . (4)  \nUsing Eqn. 4 the HI potential in the W function format is  \nV ( a ) = ( M κδ ) 4 ( W ( χ ( a )) + 1 ) 2 = ( M κδ ) 4 ( W ( χ ( a )) 2 + 2 W ( χ ( a )) + 1 ) (5)  \nThe combination of the Lambert W function scalar, HI potential and the Quintessence cosmology is hereinafter referred to as the LHQ cosmology and universe.  \nFigure 1 shows small portions of the W function principal and minus branch on both sides of the big bang. Although the Lambert W function minus branch is not considered in this work it is included to show that the maximum negative value of χ ( a ) lies at the end of the principal branch and beginning of the minus branch. At the terminal point of the principal branch χ = -1 e and W ( χ ) = -1.  \nFigure 1 displays many of the unique aspects of the LHQ cosmology and its evolution. In keeping with its definition as the ratio of the size of the universe relative to the current size set as one, the scale factor a is positive on either side of the big bang. As shown later in section 3.2.1 q is negative in the post bang region and positive in the pre bang region, making χ negative in the post bang and positive in the pre bang regions. Figure 1  \nshows that The Lambert W function has the same sign as χ . This, in turn, indicates that the scalar field is real and positive in the post bang region and imaginary in the pre bang region. All of the pre bang cosmological parameters are, however, real. This divides the LHQ cosmology universe into a canonical post bang era containing the canonical universe and the noncanonical pre bang era. Although only plotted out to a pre bang χ of 5 the W function principal branch extends to infinity and the equation for the scalar is mathematically valid at all points on principal branch.  \nArecent publication by the DESI collaboration [5] has claimed evidence for a preference of a dynamical cosmology over Λ CDM, drawing more interest to dynamical universe models. Several publications, both positive and negative, have appeared relative to the DESI results, eg.[6-12], employing parametric analysis, except for [12]. This highlights the need for accurate, analytic, solutions based on physics, rather than parameterizations, for the evolution of dark energy cosmological parameters.  \nThe investigation centers on two regions of the principal branch. The first is the far past in the pre bang epoch that extends to the infinite past. The study adopts a more practical starting point at a pre bang scale factor of 1000. The region between a pre bang scale factor of one and 1000 is the far past region. The region between a post bang scale factor of one and a pre bang scale factor of one, straddling the big bang is the second region, designated the transition zone. In comparison to the far past the transition zone is miniscule, however, it includes all of the known post bang universe, the big bang, and the pre bang approach to the big bang.  \nThe existence of analytic simple equations for the dark energy scalar and potential provides the primary motivation for this investigation. These equations along with the quintessence cosmology, reviewed in Sect. 2 and the Friedmann constraints provide the tools to calculate mathematically deterministic solutions for cosmological parameters as functions of the scale factor at all points on the W function principal branch. The evolution of the dark energy parameters is determined by the principal branch track of the Lambert Wfunction illustrated in Fig. 1. The matter and radiation average densities are strictly functions of their boundary condition densities and the scale factor.  \nThe remainder of the manuscript is arranged as follows. After this introduction the basic quintessence equations are reviewed in Sect. 2. The mathematical properties of the cosmology are then established in terms of the Lambert W function in Sect. 3. These properties are the equations for the evolution of cosmological parameters as analytic functions of the scale factor. The boundary conditions are set in Sect. 4. The physical properties for the far past and transition zone are established in Sec.5. These contain both mathematical and graphical representations of the evolution of cosmological parameters. A set of conclusions are presented in Sect. 6. The study uses natural units with c , ¯ h and 8 π G set to one. The unit of mass is the reduced Planck mass Mp . In these units κ = 1 but is retained in equations to show the proper units.', '2. Flat Minimally Coupled Quintessence': 'Even though they are well known, some of the more useful quintessence equations are presented here for easy reference. The equations are expressed in the familiar form with ϕ as the true scalar. The true scalar is ϕ = M κθ for the Ratra Peebles scalar κθ used in this work. Equations 2 and 1 plus the equations in this section, the Friedmann constraints in Sec. 3.1, along with the boundary conditions in Sec 4 form a complete set of tools to produce deterministic solutions for the evolution of the primary cosmological parameters.  \nThe action is  \nS = ∫ d 4 x √ -g [ R 2 -1 2 g µν ∂µ∂νϕ -V ( ϕ )] + Sm + Sr (6)  \nwhere R is the Ricci scalar, g is the determinant of the metric g µν and V ( ϕ ) is the dark energy potential. Sm and Sr are the actions of the matter and radiation fluids. By definition the big bang has not occurred in the pre bang epoch so the radiation term is zero in the  \npre bang action. The dark energy equation of state is P ϕ ρϕ with a dark energy density and pressure of  \nρϕ ≡ ˙ ϕ 2 2 + V ( ϕ ) , P ϕ ≡ ˙ ϕ 2 2 -V ( ϕ ) . (7)  \nThe kinetic term ˙ ϕ 2 2 is often referred to as X . It follows that  \nP ϕ + ρϕ = ˙ ϕ 2 , P ϕ + ρϕ ρϕ = ( w + 1 ) (8)  \nwhich provides the useful relation  \n˙ ϕ 2 = ρϕ ( w + 1 ) . (9)  \nThe average matter and radiation densities are simply  \nρ m = ρ m 0 a 3 , ρ r = ρ r 0 a 4 (10)  \nwhere ρ m 0 and ρ r 0 are the matter and radiation densities at a scale factor of one. Only these average densities are considered in this study.', '3. The Mathematical Properties of The LHQ Universe': 'This section develops the analytic functions for the evolution of cosmological parameters in both the pre and post bang epochs. Some duplication of derivations in [3] occurs but are included in order to have most derivations in this manuscript rather than having to refer back to the previous publication.', '3.1. Friedmann constraints': 'The first and second Friedmann constraints, given below, are important tools in deriving the cosmological parameter evolution equations.  \n3 H 2 = ρ de + ρ m + ρ r . (11)  \n3 ( ˙ H + H 2 ) = -ρ de + ρ m + ρ r + 3 P 2 . (12)  \nThe radiation density term ρ r is absent in the pre bang constraints. The cosmological parameter evolutionary equations are expressed in terms of the Lambert W function when possible.', '3.2. Mathematical properties, equations and functions': 'The equations and functions needed to describe the physical properties of the LHQ universe are developed in the following sections.', '3.2.1. The dark energy scalar constants and current values': 'From [3] the constants p and q in χ are  \np = 8 ( κδ ) 2 q = -e c ( κδ ) 2 ( κδ ) 2 c = 2 ( κδ ) 2 ln ( κθ 0 ) -( κθ 0 ) 2 . (13)  \nand the present day scalar κθ 0 is  \nκθ 0 = -4 -√ 16 + 12 Ω θ 0 ( w 0 + 1 )( κδ ) 2 2 √ 3 Ω θ 0 ( w 0 + 1 ) (14)  \nFigure 2. The value of q is plotted as a function of w 0 for pre and post bang w 0 values between -0.99 and -1.01. The sign of q transitions from negative for w 0 less negative than minus one to positive for w 0 values more negative than minus one.  \n<!-- image -->  \nwhere Ω θ 0 is the ratio of the present dark energy density to the critical density.  \nFigure 2 shows the solution for q in Eqn. 13 in terms of w 0 using Eqns. 13 and 14. The post bang values of q are negative which produce a negative value of χ and a real scalar. The pre bang values are positive making the scalar imaginary. Although both the pre bang scalar and its time derivative are imaginary all of the observable cosmological parameters are real. Positive values of the W function require being on the positive portion of the principal branch. Figure 2 shows this is achieved with w 0 values more negative than minus one. This requirement is addressed in the boundary conditions for the pre bang values of w 0 in Sect. 4.', '3.2.2. The time derivative of the scalar': 'The time derivative of the scalar is  \n˙ θ = d θ da da dt = d θ da Ha (15)  \nwhich requires an equation for d θ da .  \nThe derivative of the W function W ( χ ) with respect to its argument χ is  \ndW ( χ ) d χ = W ( χ ) χ ( W ( χ ) + 1 ) . (16)  \nThe derivative of the scalar with respect to χ is then  \nd θ d χ = δ 2 ( -W ( χ )) -1 2 ( -dW ( χ ) d χ ) (17)  \nwhich results in  \nd θ d χ = δ 2 √ -W ( χ ) χ ( W ( χ ) + 1 ) . (18)  \nThe derivative with respect to the scale factor, using Eqn. 3, is  \nd θ da = d θ d χ d χ da = p δ √ -W ( χ ) 2 a ( W ( χ ) + 1 ) (19)  \nwhere p is the power of the scale factor in Eqn. 3. The time derivative of the scalar is then  \n˙ θ = p δ √ -W ( χ ) 2 a ( W ( χ ) + 1 ) Ha = p δ √ -W ( χ ) 2 ( W ( χ ) + 1 ) H = 4 √ -W ( χ ) ( κδ )( W ( χ ) + 1 ) H (20)  \nusing Eqn. 13 for the constant p . It shows that both the scalar and its time derivative are zero at the singularity where the W function is zero.', '3.2.3. The HI potential constants δ and M': 'The kinetic term is zero at the singularity making the dark energy EoS -V V = -1. This requires a thawing post bang w that starts at minus one and thaws to less negative than minus one at larger scale factors. The previous study showed values of δ near one put the zero point in the future and produced thawing evolutions. The constant δ is therefore set at 0.95 to satisfy the requirements.  \nM is a constant therefore its value can be calculated at any point on the principal branch which is taken to be a post bang scale factor of one. The time derivative of the scalar at a = 1 is  \n˙ θ 0 = p ( κδ ) √ -W ( χ ( 1 )) 2 ( W ( χ ( 1 )) + 1 ) H 0 = 4 √ -W ( χ ( 1 )) ( κδ )( W ( χ ( 1 )) + 1 ) H 0 . (21)  \nThe first Friedmann constraint at the present time is.  \n3 H 2 0 = -8 W ( χ ( 1 )) ( κδ ) 2 ( W ( χ ( 1 )) + 1 ) 2 H 2 0 +( M κδ ) 4 ( W ( χ ( 1 )) + 1 ) 2 + ρ m 0 + ρ r 0 . (22)  \nRearranging Eqn. 22 gives  \n3 H 2 0 ( 1 + 8 W ( χ ( 1 )) 3 ( κδ ) 2 ( W ( χ ( 1 )) + 1 ) 2 ) = ( M κδ ) 4 ( W ( χ ( 1 )) + 1 ) 2 + ρ m 0 + ρ r 0 . (23)  \nSolving Eqn. 23 for M 4 by replacing ˙ θ 0 with Eqn. 21 and ( κθ ) 2 -( κδ ) 2 with Eqn. 4 yields  \nM 4 = 3 H 2 0 ( 1 + 8 W ( χ ( 1 )) 3 ( κδ ) 2 ( W ( χ ( 1 ))+ 1 ) 2 ) -ρ m 0 -ρ r 0 ( κδ ) 4 ( W ( χ ( 1 )) + 1 ) 2 . (24)  \nSince κθ 0 is a function of w 0 there is a different M value for each of the two post bang values of w 0 . The calculated values of M are listed in Table 1 along with the choice of δ .', '3.2.4. The Hubble parameter': 'Rearranging Eqn. 23 and replacing H 0 with H and ( 1 ) with ( a ) provides a solution for the post bang Hubble parameter.  \nH ( a ) = √ √ √ √ √ ( M κδ ) 4 ( W ( χ ( a )) + 1 ) 2 + ρ m 0 a 3 + ρ r 0 a 4 3 + 8 W ( χ ( a )) ( κδ ) 2 ( W ( χ ( a ))+ 1 ) 2 (25)  \nThe pre bang Hubble parameter is the same as the post bang except for the absence of the radiation density and a minus sign indicating contraction.', '3.2.5. The time derivative of the Hubble parameter': 'From [3] the time derivative of the Hubble parameter, needed for the second Friedmann constraint, is  \n˙ H ( a ) = -1 2 (( M κ ) 4 ˙ θ 2 + ρ m 0 a 3 + ρ r 0 a 4 ) (26)  \nwhere the radiation density term was added to the previous study result for the post bang epoch. Equation 20 provides the equation for ˙ θ in terms of the W function.', '3.2.6. The kinetic term X': 'Equation 25 for the Hubble parameter provides the complete equation for ˙ θ in terms of the W function.  \nκ 2 ˙ θ = 4 √ -W ( χ ) κδ ( W ( χ ) + 1 ) √ √ √ √ √ ( M κδ ) 4 ( W ( χ ( a )) + 1 ) 2 + ρ m 0 a 3 + ρ r 0 a 4 3 + 8 W ( χ ( a )) ( κδ ) 2 ( W ( χ ( a ))+ 1 ) 2 . (27)  \nMultiplying through by 1 κδ ( W ( χ ( a ))+ 1 ) gives  \nκ 2 ˙ θ = 4 √ -W ( χ ) √ ( M κδ ) 4 ( W ( χ ( a )) + 1 ) 2 + ρ m 0 a 3 + ρ r 0 a 4 3 ( κδ ) 2 ( W ( χ ( a )) + 1 ) 2 + 8 W ( χ ( a )) . (28)  \nUsing Eqn. 28 the post bang kinetic term is  \nX ( a ) = -8 W ( χ ( a )) ( ( M κδ ) 4 ( W ( χ ( a )) + 1 ) 2 + ρ m 0 a 3 + ρ r 0 a 4 3 ( κδ ) 2 ( W ( χ ( a )) + 1 ) 2 + 8 W ( χ ( a )) ) . (29)  \nThe pre bang kinetic term is identical except for the absence of the radiation term. The value of ˙ θ and the kinetic term are real and positive in the post bang epoch. In the pre bang era ˙ θ is imaginary making the kinetic term negative and real.', '3.2.7. The dark energy density and pressure': 'The quintessence dark energy density and pressure evolutionary functions are given in Eqns. 7. Using equation Eqn. 29 the post bang dark energy density is  \nρ de ( a ) = -8 W ( χ ) ( ( M κδ ) 4 ( W ( χ ) + 1 ) 2 + ρ m 0 a 3 + ρ r 0 a 4 3 ( κδ ) 2 ( W ( χ ) + 1 ) 2 + 8 W ( χ ) ) +( M κδ ) 4 ( W ( χ ( a )) + 1 ) 2 . (30)  \nThe post bang dark energy pressure is  \nP de ( a ) = -8 W ( χ ) ( ( M κδ ) 4 ( W ( χ ) + 1 ) 2 + ρ m 0 a 3 + ρ r 0 a 4 3 ( κδ ) 2 ( W ( χ ) + 1 ) 2 + 8 W ( χ ) ) -( M κδ ) 4 ( W ( χ ) + 1 ) 2 (31)  \nThe pre bang functions are identical except for the absence of the radiative term.  \n3.2.8. The dark energy equation of state  \nFrom the definition of the dark energy EoS w is the ratio of the dark energy pressure and density. The dark energy EoS is  \nw ( a ) = X ( a ) -( M κδ ) 4 ( W ( χ ( a )) + 1 ) 2 X ( a ) + ( M κδ ) 4 ( W ( χ ( a )) + 1 ) 2 (32)  \nwhere X ( a ) are the pre and post bang kinetic terms.', '4. Boundary Conditions': 'Table 1 shows the boundary conditions for the pre and post bang epochs. Most are post bang values where observational evidence exists. Post bang H 0 = 73 (km/sec)/Mpc. The pre bang scale factor one Hubble parameter values are near, but not exactly 73. They do not appear in Table 1 because they are not boundary conditions. The two post bang w 0 values are -0.99 and -0.999, purposely set close to Λ CDM. Pre bang are -1.01 and -1.001 symmetric to post bang. Ω θ 0 is 0.7 at a = 1 in both epochs. The matter and radiation ratios are 0.2999 and 0.0001 post bang respectively with the pre bang matter at 0.3.  \nTable 1. The post and pre bang boundary conditions  \nThe dimensionfull units are: H 0 , ( km / sec ) / Mpc , ρ θ 0 , ρ M 0 and ρ M 0 , M 4 p , M and δ , Mp . The other boundary conditions are dimensionless.', '5. The Physical Properties of the LHQ Universe': 'This section examines the evolution of cosmological parameters in the far past region and the transition zone. The physical properties of the pre bang universe calculated by the mathematically deterministic functions from Sec.-3.2 are presented in Sect. 5.1 followed by the transition zone physical properties in Sect. 5.2.', '5.1. The physical properties of the pre bang epoch': 'Although the equations and mathematical properties of the pre bang epoch are identical to the post bang the physical properties are significantly different. Some differences are due to the purely imaginary scalar and its derivatives. A major difference is a contracting universe as the scalar rolls down the HI potential with a negative Hubble parameter. At the beginning of the far past, a pre bang scale factor of 1000, the universe is cold and starless. It remains so for the entire pre bang region that terminates at pre bang scale factor of one. Non-dark matter is most likely elementary particles, quarks, leptons, gauge bosons and the Higgs boson. The density of matter is 10 -9 of the current post bang density matter density. Dark energy dominates with a density more than 2000 times greater than the current post bang density. In the equation for the dark energy density the HI potential dominates the much smaller kinetic term. This condition continues up to the present time. The far past start scale factor of 1000 is deep into the region where both the W function and the scalar are significantly greater than one. This means that the several ( W + 1 ) terms are equal to Wand the first term of the HI potential is dominant. In the following the abscissas of the plots are in Log 10 of the scale factor to properly display the evolutions at all scale factors between one and 1000. The ordinates are linear unless stated otherwise such as in the top panel of Fig. 3.', '5.1.1. Chi, the W function and the scalar': 'Figure 3 shows the evolution of χ ( a ) , the Lambert W function, and the scalar κθ from a scale factor of 1000 to 1. The top χ ( a ) panel is a log log plot but the bottom two panels are the standard linear log plots. Time evolution is from right to left with the scale factor contracting. The Lambert W function is positive and real while the scalar is positive and imaginary. The nature of the early time W function and scalar evolution persists until it nears the transition zone where it flattens out. The evolution of χ , W ( χ ) and κθ are very similar for the two pre bang w 0 boundary conditions of -1.01 and -1.001. All three of the functions are monotonically decreasing as the scalar rolls down the HI potential. The imaginary scalar makes the square of the scalar negative instead of the post bang positive values. The flattening of the evolution of the W function and the scalar near the  \ntransition zone is an indicator that the transition zone evolution is atypical of the general LHQ evolution.', '5.1.2. The HI potential': 'Figure 4 shows the pre bang evolution of the HI potential in ( M P ) 4 . The imaginary scalar renders all terms of HI potential positive. The upper panel shows the evolution in the standard log linear format along with the evolution of the kinetic term X demonstrating the dominance of the HI potential in the pre bang region. The evolution is a smooth monotonic decline. The w 0 = -1.01 case has a slightly higher value than the w 0 = -1.001 case. The lower panel shows that the W 0 = -1.01 case crosses below the w 0 = -1.001 case as it enters the transition zone. Comparison of the two panels indicates that the HI potential decreased by roughly a factor of 2300 from a pre bang scale factor of 1000 to one.', '5.1.3. The dark energy density and pressure': 'Figure 5 shows the pre bang evolution of the dark energy density and pressure. As expected, the range of evolution is similar to that of the HI potential since both the density and pressure are potential dominated in the pre bang epoch. The magnitude of the pressure exceeds that of the density since both the kinetic and potential terms in the pressure are negative while in the density the potential term is positive but the kinetic term is negative. Both the dark energy density and pressure have almost flat evolutions as they approach the transition zone. The dark energy density and pressure drop by approximately a factor of 10 4 between a scale factor of 1000 and one.', '5.1.4. The Hubble parameter': 'Figure 6 displays the evolution of the Hubble parameter which is negative in all parts of the far past. The Hubble parameter is dark energy dominated and becomes less negative with time following the decrease in dark energy density shown in the upper panel of Fig. 5. The evolution is linear in the log linear plot from a scale factor of 1000 to a scale factor of 3. The Hubble parameter evolution is essentially the negative of the W function evolution in Fig. 3. This is due to the dominance of W ( χ ) over one in the ( W ( χ ) + 1 ) terms in Eqn. 25. The pre bang Hubble parameter evolution flattens at a scale factor of 3 as it approaches the transition zone.', '5.1.5. The time derivative of the scalar and the kinetic term': 'Figure 7 shows the kinetic term and the time derivative of the pre bang scalar. The pre bang time derivative of the scalar is negative and purely imaginary since the scalar is decreasing and purely imaginary. It has a similar evolution as previous parameters with a steep decrease in magnitude and a flattening of evolution as it approaches the transition zone. The kinetic term is negative but real since it is proportional to the square of an imaginary number.', '5.1.6. The dark energy equation of state': 'The upper panel of Fig. 8 shows the pre bang evolution of the dark energy EoS w . The most prominent feature is the negative dip at scale factors between two and three reaching w values near -2.5. The lower panel in Fig. 8 shows how the dip occurs. At scale factors near the transition zone the kinetic term plays a larger, but not dominant, role as the HI potential declines. The ratio of the pressure to the density reaches a maximum at a scale factor of three producing the dip. At smaller scale factors the kinetic term is approaching zero making w ≈ -V / V = -1 as it enters the transition zone. The overall evolution of w is from a value almost minus one at a scale factor of 1000, becoming more negative as it evolves to a minimum near -2.5 near a scale factor of 3 then returning toward minus one as it approaches the transition zone.  \nFigure 3. The w 0 = -1.01 cases are plotted with dashed lines and the w 0 = -1.001 cases with solid lines. These line styles are used in the remainder of the pre bang plots unless stated otherwise. The top panel shows the evolution of χ ( a ) in a log, log plot. The other panels are log linear plots. Time evolution is from right to left.  \n<!-- image -->  \nFigure 4. The figures show the distant past evolution of the HI dark energy potential in the pre bang epoch. The upper panel shows the evolution of the HI potential in the standard log linear plot. The smaller kinetic term is at the bottom of the panel. The lower panel shows the detail of the cross over of the potentials near the pre bang scale factor of one.  \n<!-- image -->  \nFigure 5. The figures show the evolution of the HI dark energy density and pressure in the pre bang epoch. The upper panel shows the evolution of the dark energy density and the lower panel shows the pressure.  \n<!-- image -->  \nFigure 6. The distant past evolution of the Hubble parameter in the pre bang epoch.  \n<!-- image -->  \nFigure 7. The top panel shows the evolution of the time derivative of the scalar ˙ θ . The time derivative is purely imaginary. The bottom panel show the kinetic term X = ˙ θ 2 2  \n<!-- image -->  \nFigure 8. The upper panel shows the distant past evolution of the dark energy equation of state w . The evolution of the dark energy density and pressure creating the w dip is shown in the lower panel.  \n<!-- image -->', '5.1.7. Friedmann constraints in the pre bang era': 'This section examines the compliance of the derived parameters to the two Friedmann constraints, Eqns. 11 and 12 in the pre bang zone. The pre and post bang equations and deviations are identical except the pre bang lacks the radiation density in the first constraint. The equations for fractional deviations from the first, F1, and second, F2, Friedman constraints are  \ndev = 3 H 2 -( ρ de + ρ m + ρ r ) 3 H 2 and 3 ( ˙ H + H 2 ) + ( ρ de + ρ m + ρ r + 3 P 2 ) 3 ( ˙ H + H 2 ) . (33)  \nFigure 9 shows the constraint deviations for the first and second Friedmann constraints with w 0 = -1.001. Their deviations are all less than one part in 10 15 . The nature of the fractional errors suggest that they are the digital limitations of the Mathematica calculations which are on the order of a few times 10 -16 . The plots of the w 0 = -1.01 case are similar to the -1.001 case satisfying the Friedmann constraints at a high level.', '5.1.8. Pre bang chronology': 'This section examines the predicted pre bang chronology. The time derivative of the scale factor is simply H ( a ) a for the Hubble parameter H ( a ) and the scale factor a . The time as a function of a is  \ndt = da H ( a ) a (34)  \nwhere H ( a ) is the pre bang Hubble parameter given by the negative of Eqn. 3.2.4 without the radiation term. The time in gigayears from a scalar factor of 1000 to a scale factor of a is calculated by numerically integrating Eqn. 34. Figure 10 shows the time derivative of the scale factor as a function of the scale factor. The time since a scale factor of 1000 is shown at  \nFigure 9. The panels show the error level of the first and second Friedmann constraints for w 0 = -1.001 case . The Friedmann constraints are satisfied to better than one part in 10 15 at all points in the pre bang zone.  \n<!-- image -->  \nthe bottom for scale factors of 1000, 100, 10 and 1. The time derivative of the scale factor is of course negative in the pre bang epoch and is quite large at the pre bang scale factor of 1000. The contraction from a scale factor of 1000 to 100 takes only 1.1 gigayears. The contraction to a pre bang scale factor of 1 takes 32.6 gigayears, less than three times the age of the post bang universe. The time to contract from a pre bang scale factor of one to the big bang is 14.9 gigayears, slightly longer than the age of the post bang universe. These numbers illustrate the significant difference between the rapid contraction in the very early pre bang universe relative to the contraction and expansion rates in the transition zone. It further emphasizes that the transition zone is a remarkable time in the evolution of the LHQ universe.  \nFigure 10. The time derivative of the scale factor is plotted for pre bang scale factors between 1000 and 1. The time in gigayears for the scale factors of 1000, 100, 10, and 1 is given at the bottom of the figure.  \n<!-- image -->', '5.1.9. Still a Quintessence cosmology': "The pre bang epoch has two properties associated with a phantom cosmology: a dark energy EoS more negative than minus one and a negative kinetic term. It is, however, still a quintessence cosmology rather than a phantom cosmology with well known instability issues [13-18]. In a phantom cosmology the kinetic term is negative because it has been changed to negative in Eqns. 7 to give the 'wrong sign' in the definitions of the dark energy density and pressure [19]. In this study the equations retain their quintessence form. The negative sign for the kinetic term is due to the imaginary value of ˙ θ . Since the equations retain their quintessence form the sound speed, c 2 s = P , X ρ , X , is one in both the pre and post bang epochs. From [20] this indicates that the LHQ universe is stable in both epochs.", '5.2. The physical properties of the transition zone': 'The transition zone contains the end of the pre bang epoch, the big bang and the entirety of the post bang observable universe between post bang scale factors from zero to one. It is a region of relatively little dark energy evolution but of a highly evolving matter density. The big bang is a natural consequence of the contraction of the pre bang universe to a singularity. The big bang is, therefore, no longer the origin of time and the universe. This is counter to the canonical view of the big bang being the origin of both time and the universe, eg. [21-25] and many others that invoke a universe that is the result of a quantum entanglement collapsing to a given quantum state. In LHQ the universe is in a given state far before the big bang. In this work our chosen scale factor of 1000 for the beginning of the evolutionary calculations is presumed to be after the state has been established. A quantum analysis of the establishment of the LHQ state is beyond the scope of this investigation. Two possibilities, however, are worth consideration. One is that the universe has always been in the LHQ state and does not involve a quantum entanglement. The other possibility invokes a multiverse option that there is a quantum entanglement that exists in the distant past before the big bang. New universes form from the collapse of part of these entanglements at regular intervals and we are in the one that collapsed to the LHQ state. Neither of these possibilities have been examined in this work beyond pure speculation.  \nThe big bang is, however, an essential event required to produce favorable conditions for the formation of stars and galaxies that is not possible in the pre bang epoch. Contraction to a singularity and transition to a post bang universe is not a new concept [26-30]. Several studies of cyclic universes [31-34] proposed universes that cycled through numerous big bangs, however, with some entropy issues due to the multiple occurrences. The LHQ universe is not cyclic. It contains only one big bang.', '5.2.1. Expansion normalized variable examination of the transition zone': 'Asuccessful transition from the pre bang to post bang epoch requires the universe to evolve through the big bang rather than bouncing back into the pre bang epoch such as described in [35]. Expansion Normalized, EN, variables [36] are employed to determine whether a transition is possible. The calculation is made near either side of the singularity. The quintessence kinetic x and potential y variables are defined in [36] as  \nx = κ ˙ θ √ 6 H y = κ √ V √ 3 H . (35)  \nIn the pre bang epoch x is imaginary since ˙ θ is imaginary but x is real in the post bang epoch. The value of y is real in both epochs.  \nFigure 11 shows the evolution of the EN variables in the transition zone. It indicates that a transition through the singularity from the pre to post bang epoch is possible. The figure also shows the dominance of the y potential variable over the kinetic x variable. Even though the dark energy potential V does not go to zero the y value is zero at the singularity due to the infinite Hubble parameter there. The kinetic evolution is near zero on both sides  \nFigure 11. The evolution of the EN variables x and y at the pre and post bang boundary. The arrows indicated the direction of the evolution. The numbers at the ends of the track indicate the appropriate w 0 values for the track.  \n<!-- image -->  \nof the singularity for a significant portion of the plot, particularly for the w 0 value closest to minus one.', '5.2.2. χ ( a ) , the Lambert W function and the scalar': 'The flatness of the kinetic term evolution in the EN plot is an indicator of a region of suppressed evolution of dark energy on either side of the singularity. It is due to the slow evolution of χ ( a ) for scale factors less than one. Equation 13 for χ ( qa p ) shows that the power of the scale factor is larger than 8 which dampens the evolution for scale factors less than one.  \nThe dampening of the evolution of χ ( a ) is apparent in the upper panel of Fig. 12. The middle panel of Fig. 12 shows the Lambert W function evolution and the bottom panel the evolution of the scalar κθ . The scalar has a flat evolution between scale factors of 0.5 on either side of the singularity. The w 0 = -0.99 case has significantly more evolution than the w 0 = -0.999 case showing the effect of w 0 on evolution. All three of the parameters in Fig. 12 are zero at the singularity. The change of the scalar from imaginary to real is an important aspect of the transition. The scalar is much smaller than one throughout the transition zone, particularly for the w 0 case closest to minus one. This is a significant factor in the evolution of the HI potential discussed in the next section.', '5.2.3. The HI dark energy potential': 'At the singularity the HI dark energy potential is equal to the constant term ( M κδ ) 4 since the scalar κθ = 0. Figure 13 shows that the evolution of the HI potential is essentially flat in the transition zone. The w 0 = -0.999 case deviation is only 0.1% at a scale factor of one. The slow evolution occurs because the scalar remains much less than one in the transition zone as shown in the lower panel of Fig. 12. This makes the first two terms of the HI potential in Eqn. 1 much smaller than the constant third term. Figure 13 shows that the HI potential continues to decrease as it exits the transition zone. The decrease is due to the second term of HI potential changing from positive in the pre bang epoch to negative in the post bang era as the scalar changes from imaginary to real. The flatness of the dark energy potential evolution in the transition zone is a combination of the scalar being almost  \nFigure 12. The upper panel shows the evolution of χ ( a ) through the transition zone. As was the case in the pre bang plots the w 0 = -0.999 case is the solid line and the w 0 = -0.99 case is the dashed line. The evolution is from the pre to post bang epoch, right to left. The middle panel shows the evolution of the W function and the bottom panel the evolution of the scalar.  \n<!-- image -->  \nFigure 13. The figure shows the evolution of the HI potential, Eqn. 1, in the transition zone for the two different w 0 cases. Compared to their absolute values there is relatively insignificant evolution of the HI potential in the transition zone.  \n<!-- image -->  \nFigure 14. The upper panel shows the transition zone evolution of the Hubble parameter from scale factors between 0.01 and 1.0. in units of M 2 P on either side of the singularity. The lower panel is for scale factors between 0.4 and 1.0. The tracks for the two w 0 cases are the same to the width of the plot line but not identical.  \n<!-- image -->  \nFigure 15. The figure shows the fractional deviation of the two LHQ Hubble parameter evolutions from the Λ CDMevolution for the post bang era. It is plotted with time evolution from right to left to be consistent with the previous plots.  \n<!-- image -->  \nzero near the singularity due to the W function, the presence of the constant term in the HI potential, and the dampening described in Sect. 5.2.2. As discussed in Sect. 5.2.6 this leads to an evolution of cosmological parameters in the transition zone very close to Λ CDMwith no fine tuning.', '5.2.4. The Hubble parameter': 'The negative pre bang and positive post bang Hubble parameters create the expected discontinuity marking the transition from a contracting to an expanding universe at the singularity. Figure 14 shows the transition zone Hubble parameter evolution near the singularity. The upper panel tracks terminate at a scale factors of 0.01 to show the discontinuity. The lower panel tracks terminate at scale factors of 0.4 to better show the evolution at larger scale factors The tracks for the two different w 0 values overlap at the resolution of the plot. This demonstrates the insensitivity of the Hubble parameter to w 0 , consistent with previous studies [37,38].  \nFigure 15 shows the fractional deviation of the post bang Hubble parameter from Λ CDM, defined as ( H LHQ -H Λ ) H Λ , for the two w 0 values. The maximum fractional deviation  \nFigure 16. The top panel shows the transition zone evolution of the time derivative of the scalar ˙ θ . It is negative and imaginary in the pre bang epoch and transitions to positive and real in the post bang epoch. The bottom panel shows the evolution of the kinetic term X which is negative in the pre bang epoch and positive in the post bang epoch.  \n<!-- image -->  \nis 0.006 for the w 0 = -0.99 case and only 0.00005 for the w 0 = -0.999 fiducial case. The Hubble parameter is not a sensitive discriminator between static and dynamical cosmologies for w 0 values close to minus one. This makes it hard to verify Λ CDMor falsify LHQ based solely on observations of the Hubble parameter.', '5.2.5. The time derivative of the scalar and the kinetic term.': 'The upper panel of Fig. 16 shows the evolution of the time derivative of the scalar ˙ θ in units of M 2 P and the lower panel shows the evolution of the kinetic term X = ˙ θ 2 2 in units of M 4 P . Like the scalar, ˙ θ is zero at the transition and stays near zero for scale factors of ≈ 0.5 on either side of the transition, particularly for the w 0 = -0.999 case. ˙ θ is negative and imaginary in the pre bang epoch. It becomes real and positive at the big bang transition point.  \nThe kinetic term X evolution is quite flat in the transition zone. This is consistent with the dominance of the potential in the EN variables in Fig. 11. Its value is real and negative in the pre bang epoch since it is the square of an imaginary number.', '5.2.6. The dark energy density and pressure': 'Figure 17 shows the evolution of the dark energy density and pressure. The dark energy densities for both w 0 values are identical at a post bang scale factor of one since the density is a boundary condition. Both the dark energy density and dark energy pressure have smooth transitions from the pre bang to post bang era. Their evolutions are almost flat through the transition zone but are not zero due to the constant in the HI potential. The slight increase in the dark energy density in the post bang era is primarily due to an increase in the kinetic term X as shown in fig. 16.  \nFigure 17. The top panel shows the evolution of the dark energy density in the transition zone. The bottom panel shows the evolution of the dark energy pressure. The post bang dark energy scale factor one densities are identical because it is a boundary condition.  \n<!-- image -->  \nFigure 18. The evolution of the dark energy equation of state w in the transition zone.  \n<!-- image -->  \nThe almost constant dark energy density in the transition zone, particularly for the w 0 = -0.999 case, acts like a cosmological constant. It means that LHQ has essentially the same success in matching the observations as Λ CDMin the post bang observable universe. Alarge part of this success is due to the w 0 values very close to minus one. Much larger deviations from minus one would not be as successful as shown by the higher deviation of the w 0 = -0.99 case. On the other hand if observational constraints become much more rigorous favoring, Λ CDM,the LHQ cosmology can match them by decreasing the deviation of w 0 from minus one. Even though the match to the observations may be the same the scalar field source of dark energy in LHQ is profoundly different than the cosmological constant source in Λ CDM. The very different properties of the two universes outside of the transition zone are also profound.', '5.2.7. The dark energy equation of state': 'Figure 18 shows the evolution of w in the transition zone. It shows that w is or is very close to minus one between the pre and post bang scale factors of 0.5 on either side of the singularity. This corresponds to the flat evolution of the dark energy density and pressure between the same scale factors in Fig. 17. In the post bang era w demonstrates a thawing evolution. The w 0 = -0.999 case w does not deviate from minus one until very recent times and even then it would be difficult to reliably distinguish it from minus one. Confirmation  \nFigure 19. The panels show the error level of the first Friedmann constraint for w 0 = -0.99 case in the top panel and the second Friedmann constraint in the bottom panel.  \n<!-- image -->  \nof such a deviation of w from minus one would negate Λ CDMand be significant evidence for a dynamical cosmology such as LHQ. The LHQ predictions for the evolutions of the observables such as the Hubble Parameter and the dark energy EoS meet all of the current observational constraints.', '5.2.8. Friedmann constraints in the transition zone': 'Figure 19 shows the fractional deviation for the w 0 = -0.99 case of the first Friedmann constraint in the upper panel and the second Friedmann constraint in the lower panel. All of the fractional deviations are less than 10 -15 for the first Friedmann constraint indicating that it is satisfied. The spikes in the second Friedmann constraint plots are caused by the denominator in Eqn. 33 becoming zero when ˙ H and H 2 are equal but with opposite signs. The plots are discrete points therefore the height of the spikes depends on how close the discrete plot point is to the zero point. The plot away from the spikes is similar to the first Friedmann constraint plot. Even with the spikes the second Friedmann constraint is also well satisfied. The plots for the w 0 = -0.999 case are similar.', '6. Conclusions': 'LHQ is a self consistent, analytic set of deterministic cosmological parameter evolutions. The evolutions are mathematically valid at all points on the Lambert W function principal branch. The principal branch spans from the infinite past to a terminal point in the future. LHQ has two distinct epochs, the pre big bang epoch and the post big bang epoch. The post bang epoch, from the big bang to the present time is a canonical quintessence universe with a quartic polynomial dark energy potential. As shown in fig. 1 it is after the big bang at χ and W ( χ ) of zero where both χ and W ( χ ) are negative resulting in a real and positive scalar. The bottom panel of fig. 1 shows that the entire canonical universe only occupies a minute region compared to the total span of the W function principal branch.  \nThe pre bang epoch spans from the infinite past to the big bang singularity, a vast space on the principal branch relative to the post bang space. It raises the possibility of the universe existing before the big bang. The big bang is, in fact, a natural consequence of the contraction of the pre bang universe to the big bang singularity. In order to explore a possible pre bang universe this work assumes that the mathematically valid evolutions  \nin the pre bang epoch are physically valid as well. Under this assumption the calculated cosmological parameter evolutions reveal a noncanonical universe that is greatly, perhaps infinitely, older and at early times larger, than the canonical post bang universe. The evolution of the total LHQ universe starts with a very low matter density, high dark energy density and vastly larger universe that is contracting rather than expanding. Its basic equations are the same quintessence equations of the post bang universe but it is cold, starless, and has a scalar whose value is imaginary rather than real. It is rolling down the HI potential reaching a matter singularity at a scale factor of zero. The values of χ , W ( χ ) and the scalar κθ are also zero but the HI potential is not. Since the scalar is zero the HI potential is equal to its constant term M 4 ( κδ ) 4 producing a constant dark energy density similar to Λ CDMat the instant of the big bang. The HI potential is proportional to ( W ( χ ( a )) + 1 ) 2 , eqn. 5, therefore it takes a long time for W ( χ ( a ))) to increase enough to significantly alter the HI potential. This makes the LHQ universe similar to Λ CDMin most of the transition zone.  \nThe Λ CDM like behavior in the transition zone is not achieved by fine tuning but occurs due to the selection of two important arguments, the HI potential constant term κδ and the dark energy EoS boundary values w 0 . The Λ CDMdark energy EoS is a constant w = -1. In [3] is was found that κδ values near one produced thawing dark energy EoS evolutions with w = -1 at a scale factor of zero, hence the choice of κδ = -0.95. A goal of the study was to provide post bang analytic evolutions close to Λ CDMfor use in likelihood comparisons between Λ CDMand LHQ therefore the w 0 values are set close to minus one, -0.99, -0.999, -1.001 and -1.01. This limits the evolution of w in the transition zone to be between minus one and the boundary condition w 0 values.  \nThe canonical post bang LHQ universe has significant flexibility as a dynamical cosmology template for comparison to observations. The investigation in [3] determined that dark energy EoS evolution can be thawing, as in the case here, thawing but transition to freezing, or freezing only, by changing the value of the constant κδ in the HI potential to be near one, two or three respectively. This provides the means to determine the relative likelihoods over a wide range of EoS evolutions. Varying the present day boundary conditions on w 0 increases that range. The same is true for all of the boundary conditions listed in Table 1. This is one of the significant advantages of having analytic equations as a function of the scale factor for the evolutions of cosmological parameters as opposed to numerical solutions that do not provide the same insights as analytic equations.  \nAlthough not a subject of the present investigation it is worthwhile to mention one aspect of the future evolution of the LHQ universe. The termination of the Lambert W function principal branch at χ ( a ) = -1 e is the minimum and zero point of the HI potential. Any evolution onto the minus branch requires a reduction of the scale factor and an increase in the magnitude of the HI potential. As such the end of the principal branch may be a benign stable equilibrium point for the LHQ universe. This would provide a future with no big rip or crunch. This possibility will be the subject of a subsequent investigation.  \nThe LHQ universe is an important addition to a list of cosmologies for comparison with current and the expected future cosmological observational data. It provides a noncanonical possibility of a universe that exists before the big bang which is far older and once much larger than the canonical Λ CDMor post bang LHQ universes. Exciting new observations deserve investigations that are outside of the box if we are to discover new aspects of cosmological physics not present in our current canonical universe models.  \nAcknowledgments: The author acknowledges the very useful and informative discussions with Sergey Cherkis on the mathematical properties of the Lambert W function.  \nFunding: This research received no external funding  \nData Availability Statement: No new data was produced in this study.  \nConflicts of Interest: The author declares no conflicts of interest.', 'References': '- 1. Ratra, B. and Peebles, P. 1988, Physical Review D, 37, 12\n- 2. Peebles, P. and Ratra, B. 1988, The Astrophysical Journal, 325, L17\n- 3. Thompson, R. I. 2023, Universe, 9(4), 172\n- 4. Olver, F.W.J., Lozier, D. W., Boisvert, R.F., and Clark, C.W., 2010, in NIST Handbook of Mathematical Functions, Chap. 4, p. 111, 1st ed, (Cambridge University Press, New York)\n- 5. Desi Collaboration 2024 arXiv:2024.03002 V1\n- 6. Mukherjee, P. and Sen, A. A. 2024 arXiv.2405.19178v1\n- 7. Roy, N. 2024 arXiv:2406.00634v1\n- 8. Jia, X.D., Hu, J.P., and Wang, F.Y. 2024 arXiv:2406.02019v1\n- 9. Orchard, L. and Cardenas, V. H. 2024 arXiv:2407.05579v1\n- 10. Pourojaghi, S., Mlekjani, M. and Davari Z. arXiv:2407.09767v1\n- 11. Giare, W., Najafi, M., Pan, S., Di Valentino, E., and Firouzjaee, J. T. arXiv:2407.16689v1\n- 12. Dinda, B. R., and Maartens, R. arXiv:2407.17252v1\n- 13. Caldwll, R. R. 2002, Phys. Lett. B, 545, 23\n- 14. Carroll, S. M., Hoffman, M. and Trodden, M. 2003, Phys. Rev. D, 68, 023509\n- 15. Cline, J. M., Jeon, S. Y. and Moore, G. D. 2004, Phys. Rev. D, 70, 043543\n- 16. Scherrer, R. J. 2005, Phys. Rev. D, 71, 063519\n- 17. Vikman, A. 2005, Phys. Rev. D, 71, 023525\n- 18. Ludwick, K. J. 2017, Modern Physics Letters A, 32, 1730025\n- 19. Barenboim, G., Kinney, W. H. and Morse, M. J. P. 2018, Phys. Rev. D. 98, 083531\n- 20. Erickson, J. K., Caldwell, R., Steinhardt, P. J., Armendariz-Picon, C. and Mukhanov, V. 2002, Physical Review Letters, 88, 121301\n- 21. Hartle, J. B. and Hawking, S. W. 1983, Phys. Rev. D, 28, 2960\n- 22. Hawking, S. W. 1984, Nucl. Phys. B, 239, 257\n- 23. Hartle, J. B., Hawking, S. W. and Hertog, T. 2008, Phys. Rev. D, 77, 123537\n- 24. Hartle and Hertog, T. 2013, Phys. Rev. D, 88, 123516\n- 25. Hawking, S. W. 2014, European Physics H, 39, 413\n- 26. Penrose, R., 1965, Phys. Rev. Lett., 14, 57\n- 27. Penrose, R., 1969a, Contemporary Physics, 1, 545.\n- 28. Penrose, R., 1969b, Rivista del Nuovo Dimento Numero Speziale, 1, 252\n- 29. Gasperini, M. and Veneziano, G., 1993, Astropart Phys. 1, 317\n- 30. Khoury, J., Ovrut, B. A., Seiberg, N., Steinhardt, P. J. and Turok, N., 2002, Phys. Rev. D, 65, 086007-1,8\n- 31. Linde, A. D., 1986, Physics Letters B, 175, 395\n- 32. Steinhardt, P. J. and Turok, N. 2002, Physical Review D, 65, 126003\n- 33. Steinhardt, P. J. and Turok, N. 2002, Science, 296, 1436\n- 34. Gao, C., Lu, Y. and Shen, Y.-G. 2014, General Relativity and Gravity, 46, 1791\n- 35. Veneziano, G. 2004, Journal of Cosmology and Astroparitcle Physics, 03, 004\n- 36. Bahamonde, S., Bohmer, C., Carloni, S.,Copeland, E., Fang, W., and Tamanini, N. 2018, Physics Reports, 775-777, 1-122\n- 37. Thompson, R.I., 2018, MNRAS, 477, 4104.\n- 38. Thompson, R.I., 2019, MNRAS, 482, 5448.'}