Anisotropic evolution of D-dimensional FRW spacetime

We examine the time evolution of the D=d+4 dimensional Einstein field equations subjected to a flat Robertson-Walker metric where the 3D and higher-dimensional scale factors are allowed to evolve at different rates. We find the exact solution to these equations for a single fluid component, which yields two limiting regimes offering the 3D scale factor as a function of the time. The fluid regime solution closely mimics that described by 4D FRW cosmology, offering a late-time behavior for the 3D scale factor after becoming valid in the early universe, and can give rise to a late-time accelerated expansion driven by vacuum energy. This is shown to be preceded by an earlier volume regime solution, which offers a very early-time epoch of accelerated expansion for a radiation-dominated universe for d=1. The time scales describing these phenomena, including the transition from volume to fluid regime, are shown to fall within a small fraction of the first second when the fundamental constants of the theory are aligned with the Planck time. This model potentially offers a higher-dimensional alternative to scalar-field inflationary theory and a consistent cosmological theory, yielding a unified description of early- and late-time accelerated expansions via a 5D spacetime scenario.

We examine the time evolution of the D = d + 4 dimensional Einstein field equations subjected to a flat Robertson-Walker metric where the 3D and higher-dimensional scale factors are allowed to evolve at different rates. We find the exact solution to these equations for a single fluid component, which yields two limiting regimes offering the 3D scale factor as a function of the time. The fluid regime solution closely mimics that described by 4D FRW cosmology, offering a late-time behavior for the 3D scale factor after becoming valid in the early universe, and can give rise to a late-time accelerated expansion driven by vacuum energy. This is shown to be preceded by an earlier volume regime solution, which offers a very early-time epoch of accelerated expansion for a radiation-dominated universe for d = 1. The time scales describing these phenomena, including the transition from volume to fluid regime, are shown to fall within a small fraction of the first second when the fundamental constants of the theory are aligned with the Planck time. This model potentially offers a higher-dimensional alternative to scalar-field inflationary theory and a consistent cosmological theory, yielding a unified description of early-and late-time accelerated expansions via a 5D spacetime scenario.

I. INTRODUCTION
In an attempt at unifying gravity and electromagnetism, Theodor Kaluza and Oskar Klein extended general relativity to a 5D spacetime ultimately hypothesizing the existence of an extra spatial dimension. [1] It was later shown that when the components of the 5D metric are chosen to be independent of the extra spatial dimension, the 5D vacuum field equations yield the 4D Einstein field equations of gravitation and Maxwell's equations for electromagnetic theory, plus an additional equation that determines the dynamics of a scalar field. Since the advent of Kaluza-Klein theory, the notion of incorporating extra dimensions into 4D theories to potentially explain certain physical phenomena continues to garner attention and often falls under the designation of Modern Kaluza-Klein theories (for a review, see Applequist,Chodos,and Frend and/or Wesson [2]). Since the discovery of string theory in the late 1960s, research into extra dimensional scenarios has only increased in popularity, many years after Kaluza and Klein's influential pioneering work. String theory remains as a potentially viable candidate in offering a quantum theory of gravity, but is accompanied by the stringent requirement of residing in a higher-dimensional spacetime. The complete lack of observational evidence of these hypothesized extra dimensions has continued to remain a mystery. Traditionally, these extra dimensions have been hypothesized to be compactified objects, as is the extra dimension of Kaluza-Klein theory, and hence possibly elude observation due to having an incredibly tiny size. In the late 1980s, Dai et al. and Horǎva offered an alternative explanation for the hidden extra dimensions with the discovery of D(irichlet)-branes as fundamental extended objects, with the idea being that our universe may be a D3-brane residing in a higher-dimensional spacetime. [3] This discovery has lead to several different braneworld scenarios with large extra spatial dimensions. [4] Over the last few decades, cosmology has evolved into more of a precision science as recent observations have continued to test the predictions of theoretical cosmology. Observations of distant Type Ia supernovae, precise measurements of the size and spectral distribution of the cosmic microwave background, and the accurate mapping of galaxies in the night sky collectively indicate that the universe, on the largest of scales, is flat, homogenous and isotropic, and is currently undergoing an epoch of accelerated expansion. [5][6][7] These observations, when compared with standard 4D Friedmann-Robertson-Walker (FRW) cosmology, have constrained the six independent parameters of the Lambda-CDM model and, consequently, have upheld the model as a consistent cosmological theory. The current accelerated expansion of the universe can be understood through hypothesizing the existence of a vacuum * Electronic address: chmiddle@coloradomesa.edu † Electronic address: babrouse@mavs.coloradomesa.edu ‡ Electronic address: sdjackson@mavs.coloradomesa.edu energy density, which must comprise ∼ 70% of the total energy and matter content of the universe, as determined through this parameter fitting. According to standard FRW cosmology, a non-zero vacuum energy would dominate over pressureless matter and radiation in the late universe and give rise to an accelerated expansion via a negative pressure. Vacuum energy is predicted in quantum field theory as resulting through the continuous creation and annihilation of particle-antiparticle pairs. However, the ratio of the theoretically estimated value to that which is determined via observation and parameter fitting equates to a factor of ∼ 10 120 ; this ratio remains as the second worst discrepancy between theory and experiment in history, after the ultraviolet catastrophe of the blackbody spectrum of the late 19th/early 20th century. [8] Although remarkably successful at yielding a consistent cosmological picture of the complete time evolution of our universe, standard FRW cosmology leaves several mysteries associated with physical observations unanswered. Although a flat, homogeneous and isotropic universe is certainly allowed within the framework of general relativity, these features force a seemingly unnatural fine tuning of the initial conditions of the hot big bang model in the very early universe. Alternatively, a hypothesized epoch of early-time accelerated expansion, or inflation, neatly explains these phenomena, without the need for fine tuning. [9] Scalar-field inflationary theory, which generates the desired early-time inflationary epoch, equates as a supplement to standard 4D FRW cosmology, where the energy and matter content of the early universe is assumed to be dominated by a homogenous scalar field over the fluid contributions of radiation, matter, and vacuum energy. This hypothesized scalar field gives rise to a negative pressure, which drives an early-time inflationary epoch, when the kinetic energy of the scalar field is dominated by its associated potential energy. This inflationary epoch occurs when the scalar-field 'slowly rolls' down its potential energy curve and abruptly ends when the slow roll conditions cease to be met, leading to an era of reheating and particle production. Although highly successful at solving the horizon and flatness problem, and at giving an explanation for the origin of structure in the universe, the theory of inflation is incomplete as the potential energy function is left undetermined, although it should be noted that these features exist for a wide range of potentials. Time scales describing the beginning and end of inflation, each respectively defined by the domination of the potential or kinetic energy of the scalar field, are hence left unpredicted by the theory but are hypothesized to equate to when the universe was approximately 10 −36 to 10 −34 seconds old.
Shortly after the advent of scalar-field inflationary theory, higher-dimensional cosmological models were studied as to whether they could potentially account for the large observed entropy of the universe. [10] These models were shown to predict inflation of the 3D space during a radiation-dominated epoch of the early universe, hence possibly offering an alternative explanation for inflation to that of scalar-field inflationary theory. There, they considered a higherdimensional spacetime described by a Robertson-Walker metric, where the positively-curved, higher-dimensional compact spatial manifold was allowed to evolve at a different rate than that of the flat 3D noncompact space. There it was shown that when the higher-dimensional compact space undergoes an era of recollapse towards a (suspected) minimum value, the 3D scale factor rapidly expands. Since the publication of these original papers, several additional articles have addressed inflation within higher-dimensional spacetime scenarios, each subjected to varying assumptions and/or approximations. [11,12] In 2002, Mohammedi offered a possible higher-dimensional alternative explanation to that of vacuum energy driving a late-time accelerated expansion via dynamical compactification, where the higher-dimensional scale factor was assumed to evolve as the inverse of a power of the 3D scale factor. [13] There, the metric described in the aforementioned paragraph is subjected to a perfect fluid stress-energy tensor, where the pressure in the higher-dimensional space is allowed to differ from that of the 3D space. The D-dimensional FRW field equations were shown to reduce precisely to a 4D form once an effective pressure is defined. The solution for the 3D scale factor was found to undergo accelerated expansion in the late universe for positive 3D and higher-dimensional pressures under this assumption of dynamical compactification, where the higher-dimensional scale factor evolves through an inverse power law. Additional papers on anisotropic evolution have emerged in recent years, some including extensions of this dynamical compactification scenario. [14][15][16][17] Other papers have discussed observational constraints placed on large extra dimensions via the Large Hadron Collider and multi-messenger gravitational wave events. [18] What this manuscript offers beyond that already found in the literature begins with the discovery of the exact solution to the D-dimensional Einstein field equations subjected to a flat, anisotropic Robertson-Walker metric, where the 3D and higher-dimensional scale factors are allowed to evolve at different rates. Here, the energy-matter content of the D-dimensional universe is treated as a perfect fluid, where the pressures in the 3D and higher-dimensional spaces are allowed to differ in the general case. By adopting two equations of state linearly relating the 3D and higherdimensional pressures to the density, the D-dimensional FRW field equations are decoupled and an exact expression relating the higher-dimensional scale factor to a function of the 3D scale factor is obtained. This decoupling allows for the construction of the effective 4D FRW field equations, written solely in terms of the 3D scale factor, which ultimately can be solved exactly for a single fluid component. This exact treatment allows for the identification of two limiting regimes, which ultimately describes the behavior of the early-and late-time universe.
The fluid regime solution is found to closely mimic that of 4D FRW cosmology, offering a late-time behavior for the 3D scale factor after becoming valid in the early universe, and can give rise to a late-time accelerated expansion driven by vacuum energy. Interestingly, the lowest-order contribution of the series solution of this regime equates to a generalization of the dynamical compactification scenario of Mohammedi [13]. Here it is shown that the higherdimensional scale factor evolves as the inverse of a power of the 3D scale factor for a limited range of 3D and higher-dimensional equation of state parameters. We further show that accelerated expansion of the 3D scale factor is only obtained for a limited range of negative 3D and higher-dimensional equation of state parameters in this generalized treatment and that if the 3D scale factor is in fact undergoing accelerated expansion in this regime, then the higher-dimensional scale factor must be simultaneously expanding. The volume regime yields two branches of solutions, where two distinct cases arise for each branch. These correspond collectively to a generalization of the D-dimensional accelerating vacuum solutions of Levin and decelerating vacuum solutions of Chodos and Detweiler, where the time dependence of the 3D and higher-dimensional scale factors are merely functions of the number of spacetime dimensions. [11,19] In this manuscript, we show that both cases of each branch of solutions can in fact arise, but only for a limited range of 3D and higher-dimensional EoS parameters, w and v. For the unique case of d = 1, the volume regime solution is shown to offer a very early-time epoch of accelerated expansion, when radiation is the dominant energy component, for any negative value of the higherdimensional equation of state parameter. This volume regime solution is shown to turn on and then off, hence offering a natural entrance and exit from a possible inflationary epoch.
The fluid regime solution is then shown to turn on only after the volume regime solution turns off and consequently remains on indefinitely. This time ordering is obeyed so long as the predicted constant initial 3D scale factor of the volume regime is larger than a threshold value, where this lower bound is a function of the higher-dimensional equation of state parameter. Further, we show that the time scales marking the end of the volume regime and the beginning of the fluid regime are constrained to reside within a predicted range, with both times scales bounded from below and in some cases from above. By aligning the fundamental time constant of the theory with the Planck time and the corresponding constant initial 3D scale factor with the value predicted by standard 4D FRW cosmology at this time, we show that the aforementioned time scales marking the end of the volume regime and the beginning of the fluid regime take on maximum values of t vol,f,max ∼ 10 −35.5 s and t flu,i,max ∼ 10 −13 s. Curiously, these time scales match remarkably well with the times predicted for the unification of the strong and electroweak force and the unification of the electromagnetic and weak force, respectively. This paper is organized as follows. In Sec. II, we subject the D = d + 4 dimensional Einstein field equations to a flat, anisotropic Robertson-Walker metric. By adopting two equations of state, we arrive at an expression for the higher-dimensional scale factor as a function of the 3D scale factor. In Sec. III, we write the D-dimensional FRW field equations exclusively in terms of the 3D scale factor, arriving at a set of effective 4D field equations. These equations can be solved exactly, however, we withhold the presentation of the exact treatment until Appendix A. In Sec. IV, we present the fluid regime solution, which equates to a generalized treatment of the dynamical compactification scenario of Mohammedi [13]. In Sec. V, we present the volume regime solutions, which equate to a generalized treatment of the D-dimensional accelerating vacuum solutions of Levin and decelerating vacuum solutions of Chodos and Detweiler.[11,19] In Sec. VI, we study the strong inequalities that define the fluid and volume regimes for d = 1. We show that one of the two volume regime solutions gives rise to an epoch of accelerated expansion in the very early universe for the 5D case and offers a possible higher-dimensional alternative to scalar-field inflationary theory. Finally, in Sec. VII we summarize our results.

II. ANISOTROPIC D-DIMENSIONAL FRW COSMOLOGY
We begin with the Einstein field equations in D = d + 4 spacetime dimensions of the form where A, B are indices that run over all spacetime dimensions and G D is the higher-dimensional Newtonian constant. For notational simplicity, we set this coupling constant 8πG D /c 2 equal to one. The higher-dimensional stress-energy tensor will be assumed to be that of a perfect fluid and of the form where p(t) and p d (t) are the pressures of the 3D and higher-dimensional spaces, respectively. As is obvious from Eq. (2.2), we are allowing the pressure in the higher-dimensional space to be different, in general, from the pressure in the 3D space. Hence, this stress-energy tensor describes an anisotropic perfect fluid in d + 4 spacetime dimensions.
We choose a metric ansatz of the form ds 2 = −dt 2 + a 2 (t) dr 2 + r 2 dθ 2 + sin 2 θdφ 2 +b 2 (t)(dy 2 1 + dy 2 2 + ... + dy 2 d ), (2.3) where we allow the scale factor of the higher-dimensional manifold, b(t), to evolve at a different rate, in general, than the 3D scale factor, a(t). This metric ansatz describes flat, homogeneous and isotropic 3D and higher-dimensional spaces, where we work in units where the speed of light is set equal to unity. By adopting the above perfect fluid stress-energy tensor and metric ansatz, the D-dimensional FRW field equations and the D-dimensional conservation equation are of the form linearly relating the 3D and higher-dimensional pressures to the density. Notice that the 3D and higher-dimensional EoS parameters, w and v, can in general be time-dependent, in this manuscript they are simply treated as constants. Also notice that here we limit the EoS parameters to the domain −1 ≤ w, v < 1. Lastly, notice that Eqs. (2.2), (2.3) and (2.8) equate to the main assumptions of this manuscript. Remarkably, the scale factors a(t) and b(t) can be decoupled from one another. Using Eq. (2.8) to eliminate ρ, p, and p d from Eqs. (2.4)-(2.6) and performing some algebra, we obtain an exact differential equation of the form where we defined the power Notice that ε lies within the range 0 < ε < ∞ and hence is strictly positive as here we consider EoS parameters within the domain −1 ≤ w, v < 1. Also notice that when the 3D and higher-dimensional EoS parameters, w and v, are set equal to one another, the aforementioned power takes on the value n = −1. See Fig. 1 for a plot of the power n versus ε, where both are defined in Eq. (2.10). Notice that n takes on positive values in two distinct regions specified by ε < 2/3 and ε > d/(d − 1) and negative values in the region defined by 2/3 < ε < d/(d − 1) for all d > 1. Also notice that n diverges as ε → d/(d − 1), approaching either negative or positive infinity as ε approaches d/(d − 1) from the left or right, respectively. Interestingly, it is noted that n only takes on values n ≤ 2/d or n > 3/(d − 1) for our allowed EoS parameter space. Lastly, for the special case of d = 1, n simplifies to a linear expression in ε as the denominator of Eq. (2.10) reduces to one and, consequently, n remains positive only for ε < 2/3. This special case of d = 1 was studied by one of the authors in a previous work. [14] Equation (2.9) can be integrated twice where the higher-dimensional scale factor can be written as a function of 3D scale factor. Performing the aforementioned integrations, one finds an expression for the higher-dimensional scale FIG. 1: Plot of n versus ε for d > 1. Notice that n takes on positive values in two distinct regions (where ε < 2/3 and ε > d/(d − 1)), negative values in the region where 2/3 < ε < d/(d − 1), and diverges as ε → d/(d − 1). Also notice that n only takes on values given by n ≤ 2/d or n > 3/(d − 1) for 0 < ε < ∞. The hashed region indicates the forbidden range of values for n as ε > 0. This plot was generated for d = 6, however, the same basic features persist for all d > 1.
factor of the form b(t) = 1 a n (t) where γ 0 and γ 1 are constants of integration and n is defined in Eq. (2.10). Notice that although γ 1 is a unitless numerical scale, γ 0 has units of inverse seconds. It should be emphasized that Eq. (2.11) is an exact expression for one fluid component, relating the higherdimensional scale factor to a function of the 3D scale factor. As is obvious from Eq. (2.11), the higher-dimensional scale factor is determined by two competing terms, each characterized by a corresponding integration constant. Notice that if γ 0 is set equal to zero, we obtain a generalization of the dynamical compactification scenario of Mohammedi where, there, n was an arbitrary power that was held positive to ensure contraction of the higher-dimensional scale factor. [13] Here, dynamical compactification emerges naturally as a special case of an exact treatment for a limited range of 3D and higher-dimensional EoS parameters as the value of n is determined by w and v via Eq. (2.10).
The discovery of Eq. (2.11) allows us to decouple the D-dimensional FRW equations and obtain a set of differential equations written exclusively in terms of the 3D scale factor. In the next section, we first present this set of effective 4D FRW field equations prior to obtaining the exact solution for the 3D scale factor and, more illuminating, the approximate solutions in two distinct regimes.

III. EFFECTIVE 4D FRW FIELD EQUATIONS
Now that we have an expression for the higher-dimensional scale factor in terms of the 3D scale factor, we can write the D-dimensional FRW field equations solely in terms of a(t). After performing some algebra, we find that Eqs. (2.4)-(2.7) can be written in the form where we defined the parameters where we used Eq. (2.10) and defined the roots of η 1 as It is noted that these roots are positive for all d. As is evident from Eq. (3.5), the parameter η 1 takes on positive values in two distinct regions, where n > α + and n < α − , and negative values in the region where α − < n < α + . Notice that for the special case of d = 1, η 1 becomes linear in n as the α + root diverges whereas α − becomes unity. For d > 1, the α ± roots take on a limited range of values constrained by 3/(d − 1) < α + and 0 < α − < 2/d. Notice that Eq. (3.2) represents an effective pressure and equates to a linear combination of Eqs. (2.4)-(2.6), which can be shown with the help of Eq. (2.11). This effective pressure presented in Eq. (3.2) is generated through a combination of the density and pressures given bỹ (3.8) wherew represents an effective EoS parameter. By defining this effective pressure, notice that the D-dimensional conservation equation given by Eq. (3.4) reduces to a form that precisely mimics the 4D conservation equation with w replaced withw once the integration constant γ 0 is set equal to zero or for the special case of v = −1. This effective pressure was originally defined elsewhere, but proves useful in our analysis of the following sections. [13] Also notice that if one redefines the coupling constant (which was previously set equal to one for notational simplicity) to absorb the parameter η 1 , Eqs. (3.1) and (3.2) also mimic that of 4D FRW cosmology with p replaced withp, when γ 0 is set equal to zero. In arriving at the effective 4D FRW field equations of Eqs. (3.1)-(3.3), we defined the higher-dimensional volume element, x, as where we again used Eq. (2.11) to express this volume element solely in terms of 3D scale factor, a(t). Hence, the effective 4D FRW field equations are non-linear ODEs that contain integrals of the scale factor via the γ 0 /x terms (and squares of the integral of the scale factor via the γ 2 0 /x 2 terms for d = 1) and therefore effectively equate to thirdorder differential equations for the integral. It is interesting to note that the γ 2 0 /x 2 terms are somewhat analogous to the curvature terms of standard non-flat 4D FRW cosmology. These volume terms compete for dominance with the density and pressure terms as they collectively determine the time evolution of the scale factor. It is further noted that these curvature-like γ 2 0 /x 2 terms are only present for d > 1, which add to the intricacy of the field equations.
As a check of the effective 4D field equations, one can solve Eqs. (3.1) and (3.2) forȧ 2 /a 2 andä/a in terms of the density, pressure and the higher-dimensional volume element and then substitute these expressions into Eq. (3.3). Using Eq. (2.8), the parameter relations and performing some algebra, one can in fact verify that Eq. (3.3) is a redundant expression, as should be expected. Similar to 4D FRW cosmology, the D-dimensional conservation equation yields an expression for the density as a function of the 3D scale factor. Integrating Eq. (2.7) and using Eq. (2.11), we arrive at an expression for the higher-dimensional density of the form , (3.11) where ρ 0 is another constant of integration. Notice that in order to obtain a constant energy density for our allowed domain of EoS parameters, we must set w = v = −1; this becomes evident upon close examination of Eqs. (2.10), (3.8) and (3.11). Also notice that when the higher-dimensional EoS parameter is set to v = −1, Eqs. (2.7) and (3.4), and consequently Eq. (3.11), reduce precisely to the 4D form. This special case equates to the first law of thermodynamics for cosmology in 4D. Hence, for the special case of v = −1, the dynamical evolution of the density is completely independent of the extra dimensions. Before proceeding to a general treatment of the effective 4D field equations, we first discuss the conditions necessary for accelerated expansion. Combining Eqs. (3.1) and (3.2), we arrive at an expression for the acceleration of the 3D scale factor of the formä As is obvious from Eq. (3.12), the value of this acceleration is determined by three competing terms. When the first term on the right-hand side dominates, one obtains accelerated expansion for −1 ≤w < −1/3, η 1 > 0, and a positive energy density. In the next section, we present the series solution for the 3D scale factor in this aptly named fluid regime, which yields this accelerated expansion scenario. We show that accelerated expansion in this regime is only obtained for a limited range of negative 3D and higher-dimensional EoS parameters, which corresponds to higher-dimensional dark energy or vacuum energy when w = v = −1. We also show that if the 3D scale factor is in fact undergoing accelerated expansion in this regime, then the higher-dimensional scale factor must be simultaneously expanding.
For the unique case of d = 1, the third term on the right-hand side of Eq. (3.12) vanishes and the value of the acceleration is consequently determined by the two remaining terms. Interestingly, we note that the second term on the right-hand side is positive, yielding accelerated expansion if dominant, for η 1 > 0, γ 0 > 0, and an expanding 3D space. In Sec. V, we find the series solution for the 3D scale factor in this volume regime, which yields this aforementioned scenario of second-term domination and accelerated expansion when d = 1. This corresponding solution equates to one of two branches of solutions and will be referred to as the α + solution, where γ 0 > 0 for all allowed values of the EoS parameters.
Astonishingly, the effective 4D field equations of Eqs. (3.1)-(3.3) can be solved exactly for a single fluid component, meaning they can be integrated thrice for an integral of the 3D scale factor raised to a power. The resultant expression relates a function of the 3D scale factor, namely a product of a power of an integral of the 3D scale factor and a hypergeometric function, to the time and proves to be rather cumbersome in a general analytical treatment. As hypergeometric functions contain singularities, and these singularities can be expanded about by writing in terms of a hypergeometric series, the exact solution reduces to more convenient series solutions in ultimately two distinct regimes. This exact treatment is withheld until Appendix A.
Equivalently, one can solve the effective 4D FRW field equations perturbatively in each of these two distinct regimes, each characterized by a set of strong inequalities that ultimately define their relevant time regimes. In the following sections, we present the series solutions for the 3D and higher-dimensional scale factors in the fluid and volume regimes. We show that the fluid regime solution closely mimics that of standard 4D FRW cosmology, gives rise to a late-time accelerated expansion driven by a higher-dimensional vacuum energy, and can become valid in the early universe following an even earlier epoch of volume regime solution validity. We then show that the volume regime yields two branches of solutions, where two distinct cases arise for each branch, which collectively correspond to a generalization of the D-dimensional accelerating vacuum solutions of Levin and decelerating vacuum solutions of Chodos and Detweiler. [11,19] We show that both cases of each branch of solutions can arise, but only for a limited range of EoS parameters. For the unique case of d = 1, we show that one of the two volume regime solutions offers a very early epoch of accelerated expansion for the 3D scale factor for any 3D EoS parameter value, w, as long as v < 0. This volume regime solution is valid in an early-time radiation-dominated epoch and gives rise to a potential higher-dimensional alternative to scalar-field inflation, as the solution is shown to naturally turn on and off, hence offering a graceful exit from an inflationary epoch. This solution is followed by the relatively late-time fluid regime solution. The time scales for this transition are shown to be within a small fraction of the first second when the fundamental constants of the theory are aligned with the Planck scale, with the bounds determined by the values of the constant initial 3D scale factor and the higher-dimensional EoS parameter.

IV. THE FLUID REGIME
Here we are interested in the approximate solution to the effective 4D FRW field equations of Eqs. (3.1)-(3.3) for the 3D and higher-dimensional scale factors in the regime when In this regime, the density and pressure components dominate over terms involving the inverse of the higherdimensional volume element. The effective field equations reduce to a form similar to that of standard 4D FRW cosmology at lowest order. Here, the time evolution of the 3D scale factor closely mimics that predicted by standard 4D cosmology and the solution in this regime is shown to become valid in the early universe and remain so indefinitely. It should be noted that the above strong inequalities are in reference to the magnitude of the terms only; at this point the constants of integration γ 0 , γ 1 can be of either sign. In this regime, the higher-dimensional scale factor of Eq. (2.11) takes on the approximate form b(t) ∼ 1 a n (t) (4.2) at lowest order. Hence, this regime generalizes the dynamical compactification scenario first investigated by Mohammedi[13], where here the value of n is determined by the 3D and higher-dimensional EoS parameters as demanded by Eq. (2.10). The effective 4D FRW field equations can be solved perturbatively by writing the 3D scale factor as a series solution of the form where a 0 (t) is the zeroth-order solution, a 1 (t) is the first-order correction term, etc. Notice that in this regime, the γ 0 terms in the effective field equations of Eqs. , grouping by order of γ 0 /γ 1 and solving, we obtain the desired series solution for the 3D scale factor, which takes the form where the integration constantã 0 > 0 and should not be confused with the zeroth-order perturbative solution, a 0 (t), as defined in Eq. (4.3). Notice thatã 0 can be fixed by normalization so that a(t 0 ) = 1 when evaluated at the present time. Also notice that to lowest order, this solution for the 3D scale factor has the same functional form as that of standard 4D FRW cosmology with w replaced with the effective EoS parameter,w, which was previously defined in Eq. (3.8). This effective EoS parameter can potentially lie within the range −1 ≤w < −1/3, hence yielding accelerated expansion for the 3D scale factor, for positive values of w and v. We will explore this possibility later in this section. For Eq. (4.4) to be a valid and consistent solution to all three of the effective 4D FRW field equations, the integration constantã 0 is found to obey the algebraic relation Notice that Eq. (4.5) agrees with the relation one obtains from standard 4D FRW cosmology with the exception of the additional factors γ 1 and η 1 . Now, inserting Eq. (4.4) into Eq. (2.11), the higher-dimensional scale factor takes the form where the coefficient of the higher-dimensional scale factor in this fluid regime is determined bỹ Notice that if one demands that the higher-dimensional scale factor be positive and real, then we arrive at the requirement that γ 1 > 0. Also notice that γ 1 plays the interesting role of an overall numerical factor that scales the time-dependent expression for the higher-dimensional scale factor given by Eq. (4.6). Although the sign of n determines whether the higher-dimensional scale factor expands or contracts at lowest order, a vanishingly small value of γ 1 can scale away its overall significance in the line element, given by Eq. (2.3), and possibly offers an explanation for its absence from observation.

A. Positive energy density
The requirement that the higher-dimensional energy density be positive is equivalent to the parameter constraint η 1 > 0. This can be easily witnessed through Eq. (4.5) as both γ 1 andã 0 are positive as we require the 3D and higher-dimensional scale factors to be both positive and real. Employing Eqs. (2.10) and (3.5), we find two possible cases of EoS parameter inequalities given by where the parameters α ± were defined in Eq. (3.7). Notice that these two distinct cases arise from the fact that η 1 > 0 if n < α − or if n > α + . Special care must be taken in arriving at Eqs. (4.8) and (4.9) as the denominator of n can be positive or negative, which can be easily seen through examination of Eq. (2.10). Figure 2 shows a plot of the allowed parameter space subjected to the two cases of dual inequalities of Eqs. (4.8) and (4.9). The hashed regions equate to η 1 < 0, hence a negative energy density, and are not considered in this manuscript. Although this EoS parameter space plot was generated for d = 6, the same basic features remain for all values of d > 1. Notice that for the unique case of d = 1, the upper wedge-like hashed region vanishes and all values of the EoS parameters in the second quadrant are consequently allowed.
B. Accelerated expansion of the 3D scale factor Accelerated expansion occurs for the perturbative 3D solution analyzed in this section when the effective EoS parameter resides in the range −1 ≤w < −1/3, which can be seen by examining the power of the zeroth-order solution given in Eq. (4.4) or Eq. (3.12). Employing Eqs. (2.10) and (3.8), we find accelerated expansion for the 3D scale factor when the set of inequalities are collectively satisfied. Figure 2 shows a plot of the EoS parameter space that yields −1 ≤w < −1/3. It is interesting to note that the entire parameter space considered in the manuscript obeys the first inequality, which equates to −1 ≤w. Hence, the 3D scale factor expands for all values of w and v whereas the higher-dimensional scale factor can either expand or contract, depending on the values of w and v, as discussed in the previous subsection.

C. Dynamical compactification of the higher-dimensional scale factor
We now wish to examine the EoS parameter space for regions where dynamical compactification occurs. Dynamical compactification equates to n > 0 in the fluid regime, where the higher-dimensional scale factor contracts as the 3D scale factor expands, which can be seen via Eq. (4.2). Employing Eq. (2.10), we find two possible cases of EoS parameter inequalities given by the relations These two distinct cases arise from the fact that n > 0 when the numerator and denominator of n share the same sign. Figure 2 shows a plot of the EoS parameter space subjected to both cases of the dual inequalities of Eqs. (4.11) and (4.12). It is interesting to note that the area of the blue lower-right wedge-like region is independent of the number of extra spatial dimensions whereas the area of the blue upper-left wedge-like region varies with d. For the unique case of d = 1, the blue upper-left wedge-like region vanishes and all values of EoS parameters in the second quadrant of this parameter space give rise to n < 0.
Finally, we explore the EoS parameter space constrained by the requirement that the energy density remains positive and present the range of EoS parameters that yields dynamical compactification of the higher-dimensional scale factor and accelerated expansion of the 3D scale factor. Figure 2 shows a plot of the EoS parameter space. Notice that neither of the two regions of the EoS parameter space that yield dynamical compactification overlap with the region that yields accelerated expansion of the 3D scale factor. Thus, one can conclude that if the 3D scale factor is undergoing accelerated expansion when the perturbative solution of this section is valid, then the higher-dimensional scale factor cannot be simultaneously dynamically compactifying. A similar result was found in the five dimensional case studied elsewhere. [14,15] FIG. 2: EoS parameter space plot of v versus w that indicates the parameter space values that yield accelerated expansion for the 3D scale factor or dynamical compactification of the higher-dimensional scale factor in the fluid regime (or neither). Notice that no allowed parameter values yield both. The hashed regions indicate the values of EoS parameters that yield a negative energy density and are not considered here. This plot was generated for d = 6, however, the same basic features persist for all d > 1.

V. THE VOLUME REGIME
Here we are interested in arriving at the approximate solutions to the effective 4D FRW field equations of Eqs.
hold, with the unique exception of one of two solutions for d = 1. This special case in the volume regime for d = 1 obeys a modified set of strong inequalities defining its regime of validity and will be further discussed later in this section. It is important to note that the above strong inequalities are obtained only after securing the corresponding volume regime solutions. This was the case for the fluid regime solution, where the strong inequalities of Eq. (4.1) were found to hold for the series solution of Eq. (4.4), which was realized via the method of Eq. (4.3).
In this regime, the terms involving the inverse of the higher-dimensional volume element dominate over the density and pressure components for a negative higher-dimensional EoS parameter, v < 0. Here, two distinct branches of solutions are found that give rise to 3D and higher-dimensional scale factors that undergo either decelerated or accelerated expansion or contraction. For the aforementioned unique case of d = 1, the 3D scale factor is found to be initially of constant finite size, however, a vanishingly small higher-dimensional volume element and density remain possible as the higher-dimensional scale factor can evolve from zero initial size, hence offering an early-time big bang scenario.
It is noted that the strong inequalities of Eq. (5.1) are in reference to the magnitude of the terms only. Whereas in the previous section it was found that γ 1 > 0, in this section we find that the integration constant γ 0 , with units of inverse seconds, can be of either sign for d > 1. We find that the allowed EoS parameter space is restricted in the volume regime, with the allowed regions determined by the signs of γ 0 and the 3D Hubble parameter.
In this regime, when the aforementioned dual strong inequalities are valid, the higher-dimensional scale factor of Eq. (2.11) takes on the approximate form The effective 4D FRW field equations can again be solved perturbatively, this time by writing the 3D scale factor as a series solution of the form where, as was the case in the previous section, a 0 (t) is the zeroth-order solution, a 1 (t) is the first-order correction term, etc. Notice that when Eq. (5.3) is inserted into Eqs. (3.1)-(3.3), the density and pressure terms on the left-hand side do not yield a contribution to the zeroth-order solution as these terms' lowest-order contribution is to the first-order correction term. Additionally, it is noted that the γ 1 terms in Eqs. (3.9) and (3.11) are of order (γ 1 /γ 0 ) and will not contribute to the first-order solution as these are of higher-order than the (γ 1 /γ 0 ) (1+v) terms, so long as v < 0. This contrasts with the perturbative treatment of the previous section where the γ 0 terms of Eqs. (3.9) and (3.11) contributed to the first-order solution as they are of order (γ 0 /γ 1 ) and therefore of the same order as that of the expansion parameter. Following this method, we obtain the series solutions for the 3D scale factor, which take the general form to first-order in the expansion parameter, where a 0 , c, and k are constants of integration. The coefficient of the first-order correction term takes the form where we defined the parameter combinations Notice that in arriving at Eqs. (5.4) and (5.5) we used Eq. (4.5) to exchange ρ 0 in favor of the constantã 0 , which sets the scale for the 3D scale factor in the fluid regime (see Eq. (4.4)). Also notice that the coefficient for the 3D scale factor, a 0 , of Eq. (5.4) is left undetermined by the series solution method of this section but must be held positive for a positive 3D scale factor. The series solutions of Eq. (5.4) equate to two branches of solutions in the volume regime, where the power 1/(3 − dα ± ) takes on a negative or positive value for the top or bottom solution, respectively, where two distinct cases arise for each branch of solution. The first case is defined by k > 0, where the integration constant c is arbitrary and only constrained by c ≥ 0 to ensure a real 3D scale factor for all t ≥ 0 and can therefore be set equal to zero. For this case, the zeroth-order solutions of Eq. (5.4) equate to the D-dimensional decelerating vacuum solutions of Chodos and Detweiler. [19] The 'α + solution', which corresponds to the solution of the top sign, yields decelerated contraction for d > 1 whereas the 'α − solution', which corresponds to the bottom-sign solution, yields decelerated expansion for d ≥ 1 for the 3D scale factor.
Notice that another distinct class of solutions exist if the integration constant k < 0 and c > 0 for a real 3D scale factor for all time. For this case, one can set c = 1 without loss of generality. It is noted that this zeroth-order α + solution equates to the D-dimensional accelerating vacuum solution initially realized by Levin, where the negative integration constant k is (5.8) where H 0± (0) is the corresponding zeroth-order Hubble parameter evaluated at t = 0. [11] Considering only the zerothorder vacuum contribution, the time-contant t vol,f marks the temporal end of the validity of this inflationary solution for negative k.
Notice that both the decelerating and accelerating classes of solutions yield either expansion or contraction for the 3D scale factor and one must chose the expanding solution as the physically relevant one. Although the α + solution yields an epoch of accelerated expansion for k < 0, the zeroth-order vacuum solution suffers from a graceful exit problem as it becomes singular at a finite time. [11] This problem appears nonexistent with the series solution of this section, so long as the strong inequalities of Eq. (5.1) render this regime no longer valid near t = t vol,f .
Lastly, notice that for d = 1, the zeroth-order α + solution of Eq. (5.4) becomes constant as α + diverges for an integration constant k of either sign. To understand the dynamics of the 3D scale factor for this branch of solution for d = 1, one must examine the first-order correction term. This uniquely special case is closely examined later in this section and the next. Now, inserting Eq. (5.4) into Eq. (2.11), the higher-dimensional scale factor takes the form where κ ± and δ ± were defined in Eqs. (5.5) and (5.7). Notice that the coefficient of the higher-dimensional scale factor in the volume regime is determined by the expression The series solutions of Eq. (5.9) also equate to two branches of solutions, with two distinct cases arising for each branch. These solutions yield either decelerated or accelerated expansion or contraction for the higher-dimensional scale factor in the volume regime. Figures 3(a) and 3(b) show plots of the powers of the lowest-order solutions of a(t) and b(t), namely 1/(3 − dα ± ) and −α ± /(3 − dα ± ), versus the integer number of higher dimensions, d. Notice that for the unique case of d = 1, the power of the lowest-order α + solution is zero, which corresponds to a constant initial scale factor. The series solutions for the 3D and higher-dimensional scale factors, given by Eqs. (5.4) and (5.9), respectively, generalize the D-dimensional decelerating vacuum solutions of Chodos and Detweiler and accelerating vacuum solutions of Levin, where the time-dependence of the zeroth-order solutions are independent of the EoS parameters w, v and merely functions of the number of spacetime dimensions. [11,19] In both of these papers, the vacuum field equations were solved without the exact relation between the 3D and higher-dimensional scale factors, hence, b 0 was simply an arbitrary integration constant. In this manuscript, Eq. (2.11) yields a coupling between a 0 and b 0 as witnessed via Eq. (5.10). By demanding that the 3D and higher-dimensional scale factors be positive and real for all values of d, we find two distinct possibilities for both positive and negative γ 0 . In the next subsection, we find the relevant constraints and the allowed regions of the EoS parameter space.
Similar to that of the fluid regime, the coefficient for the higher-dimensional scale factor, given by Eq. (5.10), is determined collectively by the coefficient of the 3D scale factor and the integration constant γ 0 . Notice that a vanishingly small value of γ 0 can scale away the significance of Eq. (5.9) in the line element, possibly offering an alternative explanation to compactification for the absence of detecting the hypothesized extra dimensions.

A. Positive and real 3D and higher-dimensional scale factors
The requirement that the 3D and higher-dimensional scale factors remain positive and real for all d is equivalent to the constraints γ 0 > 0, n > α + or γ 0 < 0, n < α − for the case of an expanding 3D scale factor in d > 1, where we used Eqs. (5.6), (5.8), and (5.10) and demanded that the energy density be positive in arriving at these results. Using Eq. (2.10), we find two possible cases of EoS parameter inequalities for both branches of expanding 3D solutions, which , which correspond to the powers of the lowest-order solutions for a(t) and b(t) in the volume regime, respectively, versus d. Notice that the α+ solution exhibits decelerated contraction for k > 0 and accelerated expansion for k < 0 for the 3D scale factor for d > 1, whereas the α− solution exhibits decelerated expansion for k > 0 and accelerated contraction for k < 0 for the 3D scale factor for d ≥ 1. Also notice that for the unique case of d = 1, the zeroth-order α+ solution for the 3D scale factor is constant as α+ diverges and the corresponding power is zero. The time evolution of the 3D scale factor for this unique case is thus determined by the first-order correction term and is found to exhibit accelerated expansion for any 3D EoS parameter, w. Interestingly, this special case of accelerated expansion for the 3D scale factor for d = 1 is a member of the class of solutions that exhibit decelerated contraction for d > 1, defined by k > 0. are given by of Eqs. (5.11) and (5.12) for both branches of expanding 3D solutions, with these allowed regions shaded green. In contrast, the regions indicating contracting 3D solutions are shaded red. These regions have also been constrained by the additional requirement of a positive higher-dimensional energy density. Notice that although the series solution method of this section generates k and a 0 as arbitrary constants of integration, the general treatment of the effective 4D FRW field equations, which is presented in Appendix A, yields a prediction for their values in terms of constants and parameters (see Eqs. (A8) and (A15)). There we find that the α + (α − ) branch of solutions must have γ 0 > 0 (γ 0 < 0) for a real 3D scale factor. We also find that the regions of the EoS parameter space that are shaded green yield a 3D scale factor that undergoes accelerated expansion with k < 0 for the α + solution and decelerated expansion with k > 0 for the α − solution.
As previously mentioned, the α + solution for the 3D scale factor for d = 1 is constant to lowest order, hence, one must consider the first-order correction term to understand its time evolution. In the next subsection we show that this special case of solution is a member of the class of solutions that exhibits decelerated contraction for d > 1, defined by k > 0, yet gives rise to accelerated expansion for d = 1. For this unique case, the entire EoS parameter space with positive energy density is allowed for γ 0 > 0.

B. The uniquely special case of the α+ solution in 5D
In this subsection, we wish to closely examine the α + solution in the volume regime for the special case of d = 1. Although this solution is found via the method outlined in Eq. (5.3) and is accurately described by Eq. (5.4) when k > 0 for d = 1, its regime of validity is not dictated by Eq. (5.1) but rather when the strong inequalitieṡ a 2 a 2 ρ ∼ γ 0 xȧ a asȧ a γ 0 x and γ 1 γ 0 a(t) (n−3) dt (5.13) are satisfied. In this regime, the terms in the effective 4D FRW field equations involving the inverse of the higherdimensional volume element are found to be of the same order as the density and pressure components, but only for the α + solution when d = 1. Here, terms involving the square of the Hubble parameter are found to be much smaller than those of the density and pressure at lowest order in this regime. It is interesting to note that inflation in 4D scalar field inflationary theory arises in the 'slow-roll regime' when the time evolution of the scalar field is sufficiently gradual, namely when theφ 2 term is small when compared with the potential energy. [20] In what follows, we find an early-time inflationary epoch in a regime when the Hubble parameter is small when compared with the inverse of the higher-dimensional volume element, without the need of a scalar field. Notice that by setting d = 1, the zeroth-order solution for the 3D scale factor, given by Eq. (5.4), becomes constant as α + diverges for an arbitrary integration constant k. Hence, on first glance it appears that k can take on either sign. However, notice that by setting d = 1, the zeroth-order solution for the higher-dimensional scale factor, given in Eq. (5.9), becomes linear in time and exhibits either coasting expansion or contraction, depending on the sign of k. To gain further insight into the dynamics of the higher-dimensional scale factor for d = 1, one needs to look no further than the first D-dimensional FRW field equation, given by Eq. (2.4). In this volume regime where the strong inequalities of Eq. (5.13) hold, notice that this first field equation demands that the 3D and higher-dimensional Hubble parameters share the same sign, if one demands a positive higher-dimensional density. In this section we show that the 3D scale factor expands and so must the higher-dimensional scale factor. Hence, we find that the α + solution for d = 1 has the additional constraint that k > 0. This solution is exclusively a member of the class of solutions that exhibit decelerated contraction for d > 1, which will soon be shown to exhibit accelerated expansion for any 3D EoS parameter, w.
Setting d = 1 and choosing the top sign solution, Eqs. (5.4) and (5.9) take on the slightly simplified form where the coefficients of the first-order correction term and the higher-dimensional scale factor, originally presented in Eqs. (5.5) and (5.10), now take the form (5.16) for d = 1 where we used Eq. (3.10). Notice that in arriving at Eqs. (5.14), (5.15), and (5.16), we set c = 0 and absorbed k into our definition of κ + and b 0 . Also notice that for this α + solution in d = 1, the higher-dimensional scale factor is positive and real for γ 0 > 0, as can easily be seen by Eq. (5.16). This α + solution for d = 1 in the volume regime is special for two reasons. As discussed earlier in this section, this solution for the 3D scale factor is constant to lowest order, where a 0 corresponds to its size at t = 0. This approximate solution was first found in an earlier work of one the authors where the exact solution for the 3D scale factor in the 5D treatment, which was described via a product of a power of the 3D scale factor and a hypergeometric function, was expanded around a branch point. [14] There it was found that the higher-dimensional EoS parameter was constrained by v ≤ 0 when one demands that the 3D scale factor remain both positive and real. In this work, this approximate solution in the volume regime was found through a series solution, which is valid when v < 0.
Using the perturbative solution of Eq. (5.14), the Hubble parameter and acceleration for the 3D scale factor can easily be calculated and take the formȧ (5.18) to lowest order, where we used Eq. (4.5) in simplifying the expressions. As seen in Eq. (5.16), γ 0 > 0 for the entire EoS parameter space for this α + solution when d = 1. As the constants ρ 0 and a 0 are both held positive, and this approximate solution is only valid for v < 0, Eq. (5.18) shows that the 3D scale factor exhibits accelerated expansion in this regime for any negative value of the higher-dimensional EoS parameter, v, regardless of the value of the 3D EoS parameter, w. Notice that the 3D Hubble parameter increases with time, hence a positive time rate of change, whenever v < 0 and drives the earliest expansion of the 3D scale factor from an initial constant value. This phenomenon of accelerated expansion in the volume regime for the case of k > 0 is unique to the 5D case as all other approximate α + solutions exhibit decelerated contraction via the lowest order term (see Fig. 3(a)). Interestingly, this solution evades the question of naturalness, unlike the case of solutions that yield accelerated expansion for k < 0, as there is no need to choose the expanding solution for the 3D space and the contracting solution for the higher-dimensional manifold. For this case, both the 3D and higher-dimensional scale factors expand for v < 1.
In the next section, we explore when the fluid and volume regime solutions are valid. As each regime is characterized by a set of strong inequalities, and the time evolution of the 3D scale factor is known perturbatively, we obtain time scales that define each regime of validity. We then show that the approximate solution for the 3D scale factor in the fluid regime can describe the early universe, and remains valid indefinitely, shortly after the even earlier volume regime solution comes into and fades out of existence. The approximate α + solution for the 3D scale factor in the volume regime for d = 1 describes a very early-time epoch of accelerated expansion and potentially offers an alternative to scalar-field inflationary theory.

VI. A POSSIBLE 5D ALTERNATIVE TO SCALAR-FIELD INFLATIONARY THEORY
The series solutions presented in Eqs. (4.4) and (5.14) are valid only when the respective set of strong inequalities of Eqs. (4.1) and (5.13) are satisfied. 1 Inserting these perturbative solutions into the corresponding strong inequalities, we obtain time-dependent expressions of the form , (6.3) where the first two inequalities collectively define the regime of validity for the d = 1 volume regime α + solution.
Notice that the set of strong inequalities of Eq. (4.1) can be shown to be effectively equivalent, with the help of the fluid regime solution in Eq. (4.4), and reduce to the lone expression given by Eq. (6.3), which solely defines the regime of validity for the fluid solution. Also notice that in arriving at Eqs. (6.2) and (6.3) we dropped overall multiplicative factors that depend on EoS parameters as here we're interested in obtaining the order of magnitude of the relevant time scales. In arriving at the above strong inequalities, we defined the fundamental time scale (6.4) which is present in each expression and proves to be a useful definition. Notice that once the time becomes much larger than this fundamental time scale, the strong inequality of Eq. (6.1) is consequently satisfied and the approximate solution of the volume regime, given by Eq. (5.14), becomes valid and is consequently 'turned on'. Eqs. (6.2) and (6.3) feature two competing time-dependent quantities, with the second term a ratio of the fluid to volume regime solutions at lowest order. It is interesting to note that the regimes of validity of the volume and fluid regime solutions are ultimately determined by the values of the constants t vol,i and a 0 , in addition to the EoS parameters, w and v. As the coefficient for the 3D scale factor,ã 0 , can be chosen to normalize the fluid regime solution to yield a(t 0 ) = 1 today, we setã a 0 ≡ 10 m t vol,i ≡ 10 τ s, (6.5) where t 0 ∼ 10 17 s and we parameterized the constant initial 3D scale factor and the fundamental time constant. As the time constant t vol,i sets the scale for when the volume regime solution becomes valid, we will eventually set this constant equal to the Planck time, t vol,i = t P ∼ 10 −43 s, for now we leave it as a free parameter. Notice that at the Planck time, the 3D scale factor is predicted by flat 4D FRW cosmology to have a value a 4D (t P ) ∼ 10 −30 for w = 1/3 in a radiation-dominated epoch. Therefore, we anticipate that m should have a value near m ∼ −30 if the constant initial scale factor a 0 is to align with the value predicted by standard 4D FRW cosmology for the 3D scale factor when evaluated at the Planck time.
Upon adopting the parameterizations of Eq. (6.5), Eqs. (6.1)-(6.3) can be rewritten into the form where t flu,i ≡ 10 p s, (6.8) where the powers of the aforementioned time scales are determined by the algebraic expressions (6.10) with τ a free parameter. Again, once the time becomes much larger than the fundamental time constant defined in Eq. (6.4), Eq. (6.6) becomes satisfied and the volume regime is turned on. Notice that the right-hand side of Eq. (6.7) also increases with time as v < 1. Hence, the α + volume regime solution is 'turned off' once the strong inequality of Eq. (6.7) is no longer satisfied. As the α + solution of the volume regime exhibits accelerated expansion for any value of the 3D EoS parameter, w, as long as the higher-dimensional EoS parameter is negative, and is found to turn on and off as dictated by Eqs. (6.6) and (6.7), this solution potentially describes an early-time inflationary epoch for the 3D scale factor when radiation is the dominant energy density component. Additionally, we find that the right-hand side of Eq. (6.8) increases with time as α < −1 for d = 1 for the EoS parameters in the domain −1 ≤ w, v < 1 and a positive energy density, where α was defined in Eq. (3.10). Hence, the series solution of the fluid regime, given by Eq. (4.4), is only 'turned on' once the strong inequality of Eq. (6.8) becomes satisfied. As the fluid regime is defined by only one unique strong inequality, the fluid regime solution is never turned off once turned on and offers a relatively late-time behavior for the 3D scale factor. It is interesting to note that the special case τ = q = p occurs if 11) which is equivalent to the more revealing equality a 0 =ã 0 t 2/3(1+w) vol,i . Hence, if Eq. (6.11) is satisfied there is no regime of validity for the 5D volume solution as the volume regime turns off when turning on as the relevant time scales marking these phenomena are made equivalent, which is easily witnessed by examination of Eqs. (6.6) and (6.7) with q = τ . Additionally, the time scale indicating the beginning of the fluid regime equals that of the volume regime as p = τ . For this case, the scale factor evolving as the fluid regime solution at t vol,i must equate to the value of the constant initial scale factor of the volume regime, if Eq. (6.11) is to be satisfied. Equation (6.11) equates to a temporal boundary condition between solutions when the time interval for the volume regime is collapsed zero.
FIG. 5: Plot of the threshold power m th versus v for w = 1/3, which corresponds to a radiation-dominated universe for a traceless 4D stress energy tensor. This threshold power defines a lower bound on the constant initial scale factor, a0 = 10 m , for the volume regime solution as m > m th if time ordering is to be obeyed. Notice that if the fundamental time constant is set equal to the Planck time, which is equivalent to setting τ = −43, then m th = −30 at both v = 0 and v = −1. This equates to the lower bound on the constant initial scale factor matching the value predicted by standard 4D FRW cosmology at the Planck time for a radiation-dominated universe. For −1 < v < 0, the threshold power is smaller than this value with an absolute minimum of m th,min −32.8 at v min −0.41, nearly three orders of magnitude smaller than the value predicted by 4D FRW cosmology.

A. Lower bound on the constant initial 3D scale factor
Here we are interested in the case where the fluid regime begins after the volume regime ends with the volume regime ending only after initially beginning, which is equivalent to the set of inequalities p > q > τ . These time ordering constraints are collectively met when the power of the initial constant scale factor of the volume regime is larger than a threshold value, given by m > m th ≡ 2(τ − 17) 3(1 +w) , (6.12) which is equivalent to the inequality a 0 > a 0,th ≡ã 0 t 2/3(1+w) vol,i . Namely, the value of the constant initial scale factor of the volume regime, a 0 , must be larger than what the scale factor would be at t vol,i if it were evolving as a fluid regime solution, if time ordering is obeyed. Figure 5 shows a plot of the threshold power, m th , defined in Eq. (6.12), versus the higher-dimensional EoS parameter, v, for several values of τ , where we set w = 1/3 for a radiation-dominated universe with a traceless 4D stress-energy tensor. Notice that the shape of the m th curve is independent of the chosen value of τ . Also notice that this threshold power matches the value predicted by standard 4D FRW cosmology when evaluated at the beginning of the volume regime, at both v = 0 and v = −1 for a radiation-dominated epoch described by w = 1/3. This is due to the fact thatw | v=0 = 1/3 for all w (6.13) w| v=−1 = w, (6.14) hence, a 0,th =ã 0 t 2/3(1+w) vol,i at v = 0 and v = −1. Notice that if the higher-dimensional EoS parameter takes on values in the range −1 < v < 0, the threshold power is smaller than that predicted by 4D FRW cosmology when evaluated at the beginning of the volume regime. The minimum of the threshold power, m th,min , can be found by setting the v-derivative of m th equal to zero, which yields For a radiation-dominated universe, we find m th,min = (τ − 17)(1 + 2 √ 2)/7 at v min = 1 − √ 2 −0.41 for all τ . It is noted that for the remainder of this section we limit our discussion to w = 1/3, describing a radiation-dominated universe with a traceless 4D stress-energy tensor.
B. Lower and upper bounds on the end of the volume regime and the beginning of the fluid regime Although the power of the constant initial scale factor of the volume regime is only bounded from below as m > m th , the power of the time scale describing the end of the volume regime, q, is bounded from both below and above for a given m value. Figure 6(a) shows a plot of this power, as defined in Eqs. (6.7) and (6.9), as a function of the higher-dimensional EoS parameter, v, for several values of m, where we set τ = −43 aligning the fundamental time constant with the Planck time. The hashed region equates to q < τ , which indicates the region where time ordering is not obeyed and is therefore not considered in this manuscript. Notice that the lower bound on the power of the time scale describing the end of the volume regime, q min , has two distinct possibilities depending upon the chosen value of m for a given τ , which are characterized by . (6.16) Notice that for the unique case of m = (τ − 17)/2, the lower bound occurs at both v = −1 and v = 0. Additionally, q min for m = (τ − 17)/2 satisfies both expressions in Eq. (6.16) as m = (τ − 17)/2 marks the boundary between these two distinct cases of q min . The upper bound on the power of the time scale describing the end of the volume regime, q max , can be found by setting the v-derivative of q equal to zero, which yields Hence, for a chosen value of m for a given τ , the time scale describing the end of the volume regime is constrained with its power residing within the range q min ≤ q ≤ q max , where q min and q max are determined by Eqs. (6.16) and (6.17), respectively. The power of the time scale describing the beginning of the fluid regime, p, is bounded from below and in some cases from above, depending on the chosen value of m for a given τ . Figure 6(b) shows a plot of this power, as defined in Eqs. (6.8) and (6.10), as a function of the higher-dimensional EoS parameter, v, for several values of m where we again set τ = −43. The hashed region equates to p < τ , which is not considered here as time ordering is not obeyed in this region. The lower bound on the power of the time scale describing the beginning of the fluid regime, p min , also has two distinct possibilities depending upon the chosen value of m for a given τ , which are characterized by . (6.18) Notice that for the unique case of m = (τ − 17)/2, p min satisfies both expressions in Eqs. (6.18) as m = (τ − 17)/2 marks the boundary between these two distinct cases. The upper bound on the power of the time scale describing the beginning of the fluid regime, p max also has two distinct possibilities depending upon the chosen value of m for a given τ , where only one of these possibilities is finite. FIG. 6: Plots of the powers q vs v and p vs v for w = 1/3, which corresponds to a radiation-dominated universe for a traceless 4D stress energy tensor. Here we set τ = −43, which sets the time constant that determines when the volume regime solution becomes valid, t vol,i , equal to the Planck time. These powers define the time scales t vol,f ≡ 10 q s and t flu,i ≡ 10 p s characterizing when the α+ volume regime solution 'turns off' and when the fluid regime solution 'turns on', respectively. Notice that for a traceless 5D stress energy tensor, equating to the EoS parameter relation v = 1 − 3w, similar ranges for q and p are found, however, the shape of the p vs v curve is markedly different.
The finite upper bound is found by setting the v-derivative of p equal to zero, which yields (6.19) Notice that the above set of algebraic equations are coupled, where v max must satisfy both expressions when yielding the local maximum of p. Interestingly, the largest finite p max occurs at v = −1 for m = (τ − 17)/2 with p max diverging to either positive or negative infinity at v = −1 for m > (τ − 17)/2 or m < (τ − 17)/2, respectively.
Hence, for a chosen value of m for a given τ , the time scale describing the beginning of the fluid regime is constrained with its power residing within the range p min ≤ p ≤ p max for (τ − 17)(1 + 2 √ 2)/7 < m ≤ (τ − 17)/2, where p min is determined by the first expression of Eq. (6.18) and p max is determined by Eq. (6.19). For the distinct case of m > (τ − 17)/2, the time scale describing the beginning of the fluid regime is only constrained from below with its power p ≥ p min , where p min is determined by the second expression of Eq. (6.18).
It is interesting to note that for τ = −43 and m = −30, which correspond to the constant initial scale factor, a 0 , aligning with the value of the 3D scale factor predicted by standard 4D FRW cosmology at the Planck time, the time scales determining the end of the volume regime and the beginning of the fluid regime take on maximum values given by t vol,f,max ∼ 10 −35.5 s t flu,i,max ∼ 10 −13 s, (6.20) at v max = −1/3 and v max = −1, respectively. These time scales, curiously, match remarkably well with the times predicted for the unification of the electromagnetic and weak force, or electroweak (ew) force, and the unification of the strong and electroweak force in the Grand Unified Theory (GUT) given by t ew ∼ 10 −12 s and t GUT ∼ 10 −36 s, respectively. [21] C. Lower and upper bounds on γ1/γ0 Finally, we discuss the behavior of the higher-dimensional scale factor, b(t), in 5D in the volume and fluid regimes. Notice that as the α + solution for the 3D scale factor in the volume regime undergoes accelerated expansion for −1 ≤ v < 0 and w = 1/3, which corresponds to a radiation-dominated epoch in the very early universe with a traceless 4D stress energy tensor, the higher-dimensional scale factor expands as b(t) ∼ t to lowest order. After the volume regime ends and the fluid regime begins, notice that n > 0 for w = 1/3 and −1 ≤ v < 0 (see Eq.(2.10) or Fig.2), which corresponds to a scenario where the higher-dimensional scale factor undergoes dynamical compactification (see Eq. (4.2)). Later, during a pressureless matter-dominated epoch of the fluid regime, w = 0, the higher-dimensional scale factor undergoes either dynamical compactification or decelerated expansion, depending on whether −1 ≤ v < −1/2 or v > −1/2. In the late universe, when dark energy is the dominant fluid component and the 3D scale factor undergoes accelerated expansion with −1 ≤w < −1/3, the higher-dimensional scale factor once again experiences decelerated expansion as n < 0. It is currently unclear whether the dynamical compactification of the higher-dimensional scale factor during the early fluid regime is enough to dominate over the expansion it experiences during the early volume and late fluid regimes, when the 3D scale factor exhibits accelerated expansion, to yield an extra dimension that has contracted to a small enough size that the universe becomes effectively four dimensional.
Vanishingly small values of γ 1 and γ 0 can scale away the overall significance of the higher-dimensional scale factor in the line element, given by Eq. (2.3), offering an alternative explanation for the extra dimension's absence from observation, regardless of whether the higher-dimensional scale factor undergoes dynamical compactification or expansion. Notice that the integration constants γ 1 and γ 0 act as overall scales for the higher-dimensional scale factor, b(t), in a given regime (see Eqs. (4.7) and (5.16)), and seem to be left undetermined by the theory. However, the ratio of the integration constants γ 1 /γ 0 is constrained and is bounded from below and above for a chosen value of m for a given τ . Equation (6.4) yields an expression for this ratio. Setting τ = −43 and m = −30, which corresponds to the constant initial scale factor, a 0 , aligning with the value of the 3D scale factor predicted by standard 4D FRW cosmology at the Planck time, the ratio of these integration constants is constrained to reside within the range 10 17 s ≤ γ 1 γ 0 ≤ 10 47 s with γ 1 γ 0 ∼ 10 32 s at q max , (6.21) where, again, γ 1 is a unitless numerical scale whereas γ 0 has units of inverse seconds. Again, although this ratio is constrained via Eq. (6.21), each integration constant is left undetermined by the theory and can seemingly be chosen to take on vanishingly small values, scaling away the significance of the higher-dimensional scale factor and offering an alternative explanation for its absence from observation.

VII. CONCLUSION
Here we examined the time evolution of the D = d + 4 dimensional Einstein field equations subjected to a flat Robertson-Walker metric, where the scale factors for the 3D and higher-dimensional spaces were allowed to evolve at different rates in the general case. We chose the higher-dimensional stress energy tensor to be that of a perfect fluid, where the 3D and higher-dimensional pressures were allowed to differ from one another. Following in a manner analogous to that of standard 4D FRW cosmology, we adopted two equations of state linearly relating the 3D and higher-dimensional pressures to the density. We arrived at an exact expression for the higher-dimensional scale factor, written solely in terms of a function of the 3D scale factor. This expression allowed us to decouple the higherdimensional field equations and arrive at a set of effective 4D FRW field equations, written exclusively in terms of the 3D scale factor. The higher-dimensional scale factor is determined by two competing terms, each characterized by an associated integration constant, whose dominance partially defines either regime of validity for the two series solutions.
We then examined the effective 4D FRW field equations in the general case and arrived at the exact solution, which equates to the third integral of these non-linear differential equations. This resultant expression relates a function of the 3D scale factor, namely a product of an integral function of the 3D scale factor and a hypergeometric function, to the time. This expression cannot, in general, be inverted to yield an analytical expression for the 3D scale factor as a function of time. This exact expression can, however, be expanded near the three singularities of the hypergeometric function, which ultimately yields three series solutions in two unique regimes. Equivalently, one can solve the effective 4D FRW field equations perturbatively in each of these two distinct volume and fluid regimes. Each regime is characterized by a set of strong inequalities that ultimately define the time interval when the corresponding series solution is valid. We withheld the exact treatment until Appendix A and presented the perturbative treatment within the body of this paper as the latter proves to yield a simpler path to our main results.
The fluid regime solution was found to be valid when a set of strong inequalities are obeyed; these are fundamentally characterized by when the Hubble parameter is much larger than the inverse of the higher-dimensional volume element. This regime emerges when the higher-dimensional scale factor is dominated by the first of its two terms and is consequently proportional to a power of the 3D scale factor. This, interestingly, equates to a generalized treatment of dynamical compactification first studied by Mohammedi.[13] There, n was an arbitrary power held constant to ensure contraction of the higher-dimensional scale factor, whereas here dynamical compactification occurs naturally for a limited range of 3D and higher-dimensional EoS parameters, w and v. The solution for the 3D scale factor in this fluid regime was found to have precisely the same functional form as that of standard 4D FRW cosmology to lowest order, with the EoS parameter, w, replaced with the effective EoS parameter,w. It was found that this solution offers a latetime epoch of accelerated expansion, but only for a limited range of the 3D and higher-dimensional EoS parameters, w and v. Interestingly, neither of the regions of the EoS parameter space plot that yield dynamical compactification overlap with the region that yields accelerated expansion for the 3D scale factor. Thus, one concludes that if the 3D scale factor is undergoing accelerated expansion in the fluid regime, then the higher-dimensional scale factor cannot be simultaneously undergoing dynamical compactification.
The volume regime solutions were found in the regime where the terms involving the inverse of the higherdimensional volume element dominate over the density and pressure components for d > 1. This regime emerges when the higher-dimensional scale factor is dominated by the second of its two terms and is proportional to an integral function of the 3D scale factor. The 3D scale factor in this volume regime is described by series solutions where the lowest order contribution equates to a generalized treatment of the D-dimensional accelerating vacuum solutions of Levin and decelerating vacuum solutions of Chodos and Detweiler. [11,19] These series solutions equate to two distinct branches of solutions, where the 3D and higher-dimensional scale factors can exhibit decelerated or accelerated expansion or contraction. We showed that either case of each branch of solutions arises for a limited range of 3D and higher-dimensional EoS parameters, w and v. For the uniquely special case of d = 1, the 3D scale factor takes the form of a constant to lowest order for one branch of solution; the time evolution of the 3D scale factor is hence determined by the first-order correction term at lowest order. This solution was found to obey a different set of strong inequalities, which equate to a time when the Hubble parameter was much smaller than the inverse of the higher-dimensional volume element. This solution was found to give rise to an early-time accelerated expansion for any value of the 3D EoS parameter, w, including during a radiation-dominated epoch, so long as the higher-dimensional EoS parameter is negative, v < 0. Hence, this model offers a possible alternative to scalar-field inflationary theory.
Lastly, we examined the strong inequalities that define the regimes of validity for the series solutions of the volume and fluid regimes for d = 1, which corresponds to a 5D spacetime. The volume regime is ultimately determined by two unique strong inequalities, where the corresponding solution was shown to turn on and then off in the very early universe, hence possibly offering a natural exit from an inflationary epoch. This very early-time solution is proceeded by the relatively late-time fluid regime solution, whose regime of validity is ultimately determined by only one unique strong inequality. This fluid regime solution was shown to turn on after the volume regime solution turns off and consequently remains on indefinitely, yielding a late-time solution that closely mimics that predicted by 4D FRW cosmology. We showed that this time ordering is obeyed so long as the constant initial 3D scale factor of the volume regime is larger than a predicted threshold value. This threshold value matches the value of the 3D scale factor of the fluid regime, which matches that predicted by standard 4D FRW cosmology at v = 0 and v = −1, when evaluated at the start of the volume regime. Further, we found that the time scales marking the end of the volume regime and the beginning of the fluid regime are constrained to reside within a predicted range, with both time scales bounded from below and above. By aligning the fundamental time constant, which marks the beginning of the volume regime, with the Planck time and the corresponding constant initial 3D scale factor with the value predicted by standard 4D FRW cosmology at this time, we showed that the time scales determining the end of the volume regime and the beginning of fluid regime take on maximum values of t vol,f,max ∼ 10 −35.5 s and t flu,i,max ∼ 10 −13 s. It was noted that these time scales match remarkably well with the times predicted for the unification of the strong and electroweak force and the unification of the electromagnetic and weak force, respectively, when the universe was a small fraction of a second old.
As the solutions presented in Sections IV and V describing the aforementioned phenomena are of a perturbative nature, it is left to show that the exact solution of the effective 4D FRW field equations, presented in Appendix A via Eq. (A7), can in fact generate the number of e-folds necessary to solve the horizon and flatness problem, if the model presented in this manuscript is to truly offer a higher-dimensional alternative to scalar-field inflation theory. This examination surely can be done through numerical methods. Additional work is also required to see if this model is capable of generating the irregularities necessary to lead to the formation of structure. Additionally, this model offers a possible alternative explanation for the lack of observation of the higher dimensions from that offered by dynamical compactification. The overall size of the higher-dimensional scale factor is dictated by, in part, the constants γ 0 and γ 1 for the volume and fluid regimes, respectfully. Although the ratio of these constants is constrained, their individual numeric values appear to be left undetermined by the theory. As previously suggested, vanishingly small values for γ 0 and γ 1 can scale away the significance of the higher-dimensional scale factor from the line element.
Appendix A: General treatment of the effective 4D FRW field equations Finally, we present an exact treatment of the 4D effective FRW field equations for a single fluid component. Remarkably, these equations can be integrated thrice to yield an exact solution involving the 3D scale factor. This resulting solution relates a product of a power of an integral of the 3D scale factor and a hypergeometric function to the time, and cannot, in general, be inverted to obtain a = a(t). Rather, this exact expression can be expanded in different regimes to yield approximate solutions for the 3D scale factor, which equate to the series solutions found in the fluid and volume regimes presented earlier in this manuscript.
Using Eq. (3.8) to eliminate ρ andp from Eqs. (3.1) and (3.2), factoring and performing some algebra, we obtain an expression dictating the behavior of the 3D scale factor of the form where α and β ± were previously defined in Eqs. (3.10) and (5.6) and we used the parameter relationship of Eq. (3.10). Eq. (A1) can be integrated with the help of the relation Notice that when (γ 1 /f ) δ± is small one can write the hypergeometric function of Eq. (A7) as a hypergeometric series, which is useful for obtaining the approximate behavior of Eq. (A7) in this limiting regime. Applying the series expansion of Eq. (A9) to Eq. (A7) and performing some algebra, we obtain a non-linear expression of the form where we kept terms up to order o(h 2 ) in generating the above expression and defined the quantities As we are interested in the time evolution of the 3D scale factor for small h, Eq. (A10) can be solved perturbatively by writing h as a series solution of the form where h 0 is the solution to the linear approximation of Eq. (A10), h 1 is the first-order correction term, etc. Using this method, we obtain an approximate solution for h of the form h = (c + kt) 1 − κ 2 (c + kt) + ... .
Now, using the definition of h given in Eq. (A11), we obtain the approximate solution for f , which can then be differentiated to arrive at a perturbative solution for the 3D scale factor. We obtain an approximate solution of the form a(t) = a 0 (c + kt) 1/(3−dα±) where we used Eq. (A4) in obtaining this result. The coefficient for the 3D scale factor is determined by where we used Eqs. (A4), (A8), and (A11) in arriving at these expressions. The series solutions of Eq. (A14) equate to two distinct branches of solutions in this small (γ 1 /f ) δ± regime. It is interesting to note that k and a 0 were found to be arbitrary constants of integration via the series solution method of Section V, but are determined by constants and parameters via Eqs. (A8) and (A15) through the general treatment presented here. Now, as c 0 and γ 1 have each been found to be strictly positive, we arrive at the additional constraint (1 + α)/γ 0 δ ± > 0 when one demands a real 3D scale factor. Employing Eqs. (3.7), (3.10), and (5.7), one can easily show that (1 + α)/δ + > 0 and (1 + α)/δ − < 0 for all allowed values of the EoS parameters that yield a positive energy density as −1 <w for −1 < w, v < 1. Hence, we find two possible cases that yield a real 3D scale factor described by γ 0 > 0 and (1 + α) δ + > 0 or γ 0 < 0 and (1 + α) δ − < 0.
For the α + (α − ) solution of Eq. (A14) we find that γ 0 > 0 (γ 0 < 0). Notice that when either of the above cases are met we obtain not only a real but also positive value for a 0 . This contrasts with the fact that the parameter k can be positive or negative, whose sign depends on the values of the EoS parameters w, v. Upon examination of Eq. (A8), we find that for an expanding 3D scale factor, where we used Eq. (5.6) and demanded a positive energy density. Figures 4(a) and where k and were defined in Eq. (A22). Using the solution for h, which was presented in Eq. (A13) with the parameters specified in Eqs. (A22) and (A23), we obtain an expression for f and then differentiate. We obtain an approximate solution for the 3D scale factor of the form a(t) = a 0 (c + kt) 1/(3−dα∓) · 1 − κ 2 (1 + 1/δ ± ) (dn − 3) (c + kt) −β∓δ∓ + ... , where the coefficient for the 3D scale factor is determined by where we used Eqs. (A4) and (A22) in arriving at these expressions. Notice that the series solutions of Eq. (A24) equate to two distinct branches of solutions in this small (f /γ 1 ) δ± regime and are equivalent to those presented in Eq. (A14). Upon close examination of the approximate solutions of this subsection and the last, one finds that one can arrive directly at Eqs. (A22)-(A25) from Eqs. (A8), (A11), (A14) and (A15) by making the parameter substitutions and by setting γ 1 = 1. Hence, the top sign solution of Eq. (A14) is equivalent to the bottom sign solution of Eq. (A24) and vice versa when one sets γ 1 = 1. The approximate solutions for the 3D scale factor near the z = 0 and z = ∞ singularities equate to two different branches of equivalent solutions.

f ∼ γ1
Lastly, we explore the behavior of Eq. (A7) in the vicinity of the f /γ 1 ∼ 1 branch point. In this section, we employ another hypergeometric transformation formula of the form where we again redefined the parameters k, , and h of the previous subsections as in this subsection, where we used Eq. (A26) in arriving at these expressions. As was the case in the previous subsection, the second term of the transformation equation when multiplied by f 1/β± reduces to a function of the EoS parameters and the integration constant γ 1 , which can be shown through use of the hypergeometric series of Eq. (A9). This term, which is void of the 3D scale factor, can be absorbed into the integration constant c, which effectively shifts the size of the scale factor at a given t. Now applying the series expansion of Eq. (A9) for small h to the hypergeometric function in Eq. (A28) and performing some algebra, we again arrive at an expression of the form given by Eq. (A10). In arriving at this expression we again kept terms up to order o(h 2 ) and redefined the parameter where k, , and h were previously defined in Eq. (A29). Using the solution for h, which was presented in Eq. (A13) with the parameters specified in Eqs. (A29) and (A30), we obtain an expression for f and then differentiate. We obtain an approximate solution of the form a(t) =ã 0 (c + kt) 2/3(1+w) (A31) where the coefficient for the 3D scale factor in this limiting regime is determined bỹ Notice that in the fluid regime, the parameters k andã 0 , given by Eqs. (A29) and (A32), are real when either of the cases of Eq. (A16) are satisfied. This result is precisely equivalent to that for real values of k and a 0 in the volume regime. However, in the volume regime we found that k could take on positive or negative values, with the sign of k determined by the inequalities of Eq. (A17). In the fluid regime, we find that k is strictly positive, which can easily be witnessed through examination of Eq. (A29) as γ 1 > 0. Notice that for k > 0, the integration constant c is arbitrary, only constrained by c ≥ 0 to ensure a real 3D scale factor for all t ≥ 0, and can therefore be set equal to zero. Now, inserting Eq. (A31) into Eq. (2.11), the higher-dimensional scale factor takes the form b(t) =b 0 (c + kt) −2n/3(1+w) (A33) · 1 + 1 dδ ± 1 + 2v · dn 3(1 +w)(2 + α) (c + kt) (1+α) + ... , where the coefficient for the higher-dimensional scale factor in this limiting regime is determined bỹ Notice that in calculating Eq. (A33), the integral term in Eq. (2.11) contributes to the first-order correction term and must be included. Also notice that the approximate solutions for the 3D and higher-dimensional scale factors given by Eqs. (A31) and (A33) can be compared to the series solutions obtained in the fluid regime, given by Eqs. (4.4) and (4.6). One finds perfect agreement between these solutions when the integration constant c = 0, which can be done without loss of generality as k > 0 in the fluid regime.