Exiting Inflation with a Smooth Scale Factor

The expectation that the physical expansion of space occurs smoothly may be expressed mathematically as a requirement for continuity in the time derivative of the metric scale factor of the Friedmann-Robertson-Walker cosmology. We explore the consequences of imposing such a smoothness requirement, examining the forms of possible interpolating functions between the end of inflation and subsequent radiation- or matter-dominated eras, using a straightforward geometric model of the interpolating behavior. We quantify the magnitude of the cusp found in a direct transition from the end of slow roll inflation to the subsequent era, analyze the validity several smooth interpolator candidates, and investigate equation-of-state and thermodynamic constraints. We find an order-of-magnitude increase in the size of the universe at the end of the transition to a single-component radiation or matter era. We also evaluate the interpolating functions in terms of the standard theory of preheating and determine the effect on the number of bosons produced.


Introduction
The problem of trying to reconcile physical theories regarding the form and evolution of the primordial universe with modern cosmological observations has occupied researchers for decades.Providing a complete explanation for the origins of the characteristics that the universe exhibits today has been challenging.We witness extreme uniformity and flatness and the absence of certain particles predicted in some Grand Unified Theory (GUT) models.Many physicists have detailed these well-known difficulties in texts and expository papers; see, for example, Refs.[1][2][3][4][5][6][7][8][9].
The issue of uniformity has the name the Horizon Problem.We see a homogeneous, isotropic universe on large scales.Despite the almost incomprehensible longevity of the universe, it is simply too immense to have grown to be uniform on large scales.Causality demands that local homogeneity takes time to develop, with equilibrium conditions dispersing at a rate no greater than the speed of light.Two local homogeneous elements dispersed in time must retain a causal connection to remain in equilibrium with one another.Yet opposite sides of the universe appear nearly identical to us.Even if we assume they started that way, not enough time has passed for space to have expanded a great enough distance to maintain the equilibrium-at least not according to the physics we understand-with signaling bounded by the limit of the speed of light.Cosmologists assess uniformity primarily using the temperature of the Cosmic Microwave Background (CMB), which is the thermal radiation emitted as the matter in the universe was cooling and transitioning from a conductive, opaque plasma to a neutral, transparent gas.The widely-accepted standard is T CMB ≈ 2.7255K [10].An early CMB probe, the Cosmic Background Explorer, found the temperature variation from this mean to be on the order of 10 −5 K.Each causally connected portion of the CMB that should have been able to thermalize before the recombination photons began to stream freely along their paths toward the Earth covers a solid angle of approximately 0.013 sr on the sky, so that about 10 4 such solid angles make up the CMB.How then, cosmologists ask, did ten thousand discrete portions of the CMB, which do not appear to have been causally connected at the time of recombination, collectively equilibrate at a common temperature that is uniform to within the order of 10 −5 K [1]?A period of superluminal expansion provided by inflation could provide the missing causal connection.
Expansion itself does not change the inherent topology of the universe: A universe that is closed, open, or flat remains so.However, expansion makes any curvature appear locally more flat.The Planck Collaboration has measured the spatial curvature Ω k as 0.001 ± 0.002 [11], which means our nearly flat universe presents a second cosmological difficulty, the Flatness Problem.The scale factors a(t) of matter-dominated and radiationdominated cosmologies are proportional to t 2/3 and t 1/2 , respectively.According to the first Friedman equation [1,12], Ω k scales as 1/ ȧ2 , so that in the absence of any other influences, the longevity of the universe means consistency with the Planck measurement requires that the curvature at the beginning of the radiation-dominated era must have been extraordinarily small.Again, inflation offers a remedy: The inflationary exponential scale factor would tend to drive down the spatial curvature to a level that could support subsequent evolution to the value observed today.
The additional Monopole Problem arises from the predictions of some GUT models [13,14] that a phase transition breaks the symmetry between the strong and weak forces when the temperature of the universe drops to a level consistent with the energy scale 10 16 GeV.A result would be the formation of a dust of massive magnetic monopoles, with a density that is subsequently proportional to a −3 , potentially thereby blocking the radiation and matter eras from taking place [5].The Monopole Problem calls for a mechanism to reconcile the GUT prediction of the creation of these massive particles with our accepted understanding of the chronology of the early universe and current cosmological observation.Inflation could provide dilution that would make magnetic monopoles so few and far apart that finding them would be essentially impossible.

The Inflation Solution
In his groundbreaking paper in 1981 [15], Alan Guth introduced inflation as a theory to address the inexplicable Horizon, Flatness, and Monopole Problems.However, he also acknowledged the difficulty his mechanism created: An exit from the false vacuum that drives inflation involved quantum tunneling from a false to the true vacuum state, an effect that would occur primarily in localized bubbles-that is, discrete regions subsequently characterized by the Klein-Gordon scalar field that drives inflation (the inflaton ϕ) having settled into its true vacuum state.Meanwhile, expansion of space would continue between the bubbles (where such tunneling had not yet occurred), and as a result, we would expect to see parcels of non-uniform space today.Intersecting bubbles would have similar effects.This model of inflation thus predicted a universe inconsistent with observation; Guth's original theory lacks a graceful exit.
In 1982, inflation pioneer Andrei Linde sought to solve the graceful exit problem with a new theory, slow-roll inflation [16].Instead of starting in a false vacuum, the inflaton rolls down a potential energy plateau to a minimum where it oscillates around a true vacuum state, which graph (a) of figure 1 depicts schematically.The assumption that the potential energy of the inflaton dominates the kinetic energy for a sufficient time results in exponential inflation.
The Hubble parameter H is the expansion rate of the universe, V (ϕ) is the potential of the inflaton, and the comma denotes the partial derivative with respect to ϕ.The theory assumes that the magnitudes of both the density and pressure of the inflaton become approximately equal to the potential by treating the inflaton condensate as a perfect fluid: In this regime, the second Friedmann equation [1,12], commonly known as the acceleration equation, because it governs ä, is essentially In our notation, ℏ = 1 and c = 1 hereafter unless otherwise noted.Thus the expansion of space undergoes inflationary acceleration, ä a > 0, as a result.The first Friedman equation, with H 2 = ȧ a 2 ≫ k a 2 , yields the scale factor solution This is the exponential expansion of space predicted by the theory of slow-roll inflation [8,17].A period of superluminal expansion would explain the homogeneity and isotropy of the observable universe by providing the necessary causal connection to solve the Horizon Problem.Superluminal expansion would also flatten the spatial curvature and decrease the density of magnetic monopoles.Numerical analysis provides insight into the question of the number of e-folds of expansion necessary to resolve these problems.However, this solution comes with its own associated shortcoming: Producing an outcome consistent with modern observations demands very specific initial conditions.Although the underlying physics of inflation (such as the existence of the inflaton field) remains unsubstantiated experimentally, the framework of inflation is widely accepted among cosmologists as a way of providing an underlying solution to various cosmological problems.Over the years, researchers have revised the concept by devising a diverse body of new theories.For our purposes, we shall focus on Linde's solution to the problem of the requirement of specific initial conditions in the slow roll theory.In 1983, he published his theory of chaotic inflation [18].In its simplest version [17], the inflaton potential has the form V (ϕ) = 1 2 m 2 ϕ 2 .The plateau is absent, and in an expanding universe, the friction term 3H φ in the inflaton equation of motion (1) has the effect of restricting the motion of the inflaton, as the slow roll plateau does, resulting again in exponential expansion.Figure 1 shows the two contrasting potentials.).Slow roll inflation requires a plateau to generate enough e-folds of inflationary expansion to solve the Horizon, Flatness, and Monopole Problems.In chaotic inflation, the effects of the friction term in the equation of motion replace that of the plateau in keeping the inflaton from moving to the true vacuum too quickly.After inflation, both models involve the inflaton oscillating around a minimum potential during a period of reheating, which we review in section 1.2.
Cosmologists have applied a variety of approaches to estimating the amount of inflation necessary to solve the Horizon, Flatness, and Monopole Problems.The amount by which the cosmos expands is normally expressed in terms of the number of times the size has increased by a constant factor-in other words, the number of the e-folds (or nepers) N = log a(t f ) a(t i ) .Linde [17] reports that a quadratic inflaton potential creates a wavelength for the inflaton comparable in size to our observable universe after about 61 e-folds.In a detailed analysis, Lyth [19] finds that for a quartic inflaton potential, a range of e-fold values from 47 to 61 results in a universe on the present scale.He further explains that at a minimum, more than 14 e-folds are needed to generate perturbations leading to structure formation, and that an extended period of domination of the inflaton kinetic term could increase his estimates of N by nearly 14, therefore concluding with an estimate of 14 < N < 75.
Other researchers have performed analyses to determine the number of e-folds required to solve specific inflationary problems [20].Solving the Horizon Problem entails that the comoving Hubble radius at the beginning of inflation [a(t i )H(t i )] −1 must contain what has become the comoving Hubble radius today [a(t 0 )H(t 0 )] −1 , so that the comoving [a(t 0 )H(t 0 )] −1 could have thermalized before expanding through the post-inflationary epochs of the universe up to the present.The Hubble radius is the distance light travels in time t = H −1 .Thus we have In eq. ( 10) we use the inverse relation between the scale factor and temperature, which is derived in appendix A for reference.Parameter values T (t 0 ) ≈ 2.75 K and H(t 0 ) ≈ 100 km/s/Mpc lead to which indicates that N is at least 67, because the temperature H(t i ) represents is greater than T (t f ).

Reheating
The expansion of space by inflation dilutes the number densities of all particles and leaves the universe cold, with energy concentrated primarily in the inflaton.Following the end of inflation, reheating results in the transfer of energy from the inflaton to Standard Model particles or their precursors.Reheating has two stages, first the transfer of energy and then subsequent thermalization to a temperature sufficient to promote nucleosynthesis of light elements.A mechanism developed by Lev Kofman, Andrei Linde, and Alexei Starobinsky in their iconic 1997 paper, "Towards the Theory of Reheating after Inflation" [21], which they call preheating because it precedes thermalization, supersedes earlier explanations of reheating by way of perturbation theory and narrow parametric resonance.We reference a selection of the wide range of literature available on the subject of preheating [7,[22][23][24][25].
A preheating framework appears to be necessary, because perturbative processes prove too slow and inefficient to raise the reheating temperature enough to support nucleosynthesis.Also, the perturbative approach required certain conditions and treated inflatons collectively in a state of superposition of individual particles, each capable of decaying independently-rather than as coherent semiclassical fields.On the other hand, narrow parametric resonance models followed the approach that the inflatons formed a homogeneous, coherent, oscillating wave appropriate for classical treatment.In narrow parametric resonance, an inflaton wave interacts as a background source for a second scalar field χ.However, this theory itself can be problematic.Because the modes of the scalar field χ have physical wavelengths, the expansion of space redshifts modes outside the borders of the resonance band and also makes the band more narrow.In addition, the expansion and the decay of the inflaton into χ particles decrease the amplitude of the coherent inflaton wave.The number of particles being produced instantaneously is proportional to both the number of χ particles previously created and to the inflaton amplitude, so that the effects of expansion and decay lower the efficiency of the resonant conversion and tend to suppress the growth of the χ population.Narrow parametric resonance thus typically terminates well before reheating is complete.
The parametric resonance in preheating models is instead broad: All modes less than a specific momentum participate in the ϕ-χ coupling.A non-adiabatic transfer of energy leads to exponential growth in the number and number density of the χ quanta.Moreover, the expansion of space can actually make the resonance more effective by gradually redshifting additional modes down to below the maximum momentum, making them part of the process.The end of reheating depends on the possible range of values of parameters involved in preheating and the complex dynamics of backreaction and rescattering.However, preheating may still not be sufficient to complete reheating, and the reheating process may have to revert to a period of narrow parametric resonance, perturbative decay, or both to arrive at a temperature that is suitable for thermalization but not high enough to produce very massive particles like monopoles.
In section 2, the reader will find a description of the cusp discontinuity inherent in inflationary theories involving an exponential scale factor and our approach to quantifying the extent of the cusp.Section 3 introduces a method for finding an interpolating function to replace the cusp, by detailing the geometry of a simple circular model.Then in section 4, we focus on finding a more realistic interpolating function.We derive the formalism establishing smoothness in the expansion of space at the end of inflation and analyze the implications of the most straightforward interpolating candidates, power law functions.The equation-of-state and thermodynamic constraints provide additional means of restricting possible interpolating functions, and this is discussed in section 5. We analyze the effect on the size of the universe of a horizontal parabola-like power law serving as a transitional interpolating function in section 6.Finally, we look at further numerical analyses to determine the effect on the scalar χ number and number density predicted by the Kofman, Linde, and Starobinsky (KLS) model of preheating in section 7.

Period Scale Factors
The well-known expressions for the scale factor in the early universe include a curious unphysical approximation, a lack of smoothness at the end of the inflationary epoch.The inflationary and radiation-dominated scale factors, a 1 (t) and a 2 (t), respectively, follow [1] To demonstrate the discontinuity in ȧ, we assume contrariwise that the time derivatives of the scale factors are equal at the end of inflation, t f : By first expressing the Hubble parameter H in terms of N , the number of inflationary e-folds of expansion, and then equating derivatives, we find Continuity of the derivatives requires that N ≤ 1 2 , or else the time at the beginning of inflation is less than zero.Although much research into inflation has produced a wide range of proposed values for N , this result is particularly problematic.If taken literally, it would eliminate inflation as a solution to the kinds of problems the theory was designed to solve.

Quantifying the Discontinuity
Although early universe estimates are themselves quite problematic because of the uncertainty in the values of basic parameters, using reasonable values can provide some insight into the mathematical relationship between the scale factors in different periods of cosmological evolution.We can estimate a value for a(t f ) by taking advantage of the inverse relation between the scale factor and temperature, in conjunction with estimates of temperature then and now, T (t f ) and T (t 0 ), respectively, since by convention, a(t 0 ) = 1.The temperature of the CMB today is T (t 0 ) ≈ 2.73 K [10], and T (t f ) corresponds to the temperature equivalent to the value of H for a universe that supports the Standard Model, which is H ≈ 10 16 GeV [26].
Next we compare slopes at the end of inflation.The inflationary slope is For the radiation-era derivative, after solving eq. ( 16) for t f and substituting it into eq.( 15), we have The estimate by Liddle and Lyth that inflation began at t i = 10 −42 s [2] leads to Ht i ≈ 0.02.A reasonable assumption is that N ≈ 60 [17,19,20], which results in a measure of the discontinuity.The time derivative of the radiation-era scale factor is approximately 1 120 of the derivative of the inflationary scale factor.Graph (a) of figure 2 shows this change in the growth behavior qualitatively.We also note that these values, H = 10 16 GeV and N = 60, yield an estimate for the duration of inflation without the need to specify t i or t f : The circular arc defined by the transitional scale factor a cir (t) intersects tangentially with the inflationary scale factor a 1 (t) and the nowdisplaced radiation era scale factor noted with a prime, a ′ 2 (t).The parameter ∆ is the time period from the end of inflation until the time of continuity between a cir (t) and a ′ 2 (t).The a functions are not to scale.

The Transition
What kind of transitional function could provide continuity between the two period scale factors?At the end of inflation, the slope begins to decline.For simplicity, we require a steadily declining slope with no regions in which the universe undergoes contraction.A properly chosen intermediate power law would conform to our requirement, which we explore in section 4. Initially, as shown in graph (b) of figure 2, we shall use a circular arc a cir (t) to illustrate the geometry.The arc lies tangent to a 1 (t) at the end of inflation and tangent to the now-displaced radiation-era scale factor noted with a prime, a ′ 2 (t).The transitional duration ∆ remains to be determined.
In the more detailed view of figure 3, we see five unknown variables: • R -the radius of the circular arc • ∆ -the time between t f and the tangent point at which the circular arc a cir (t) meets the displaced radiation era a ′ 2 (t) • a(t f − δ) -the a-axis value at t f − δ, aligned with the center of the arc • δ -the measure of the t-axis displacement corresponding to the difference between a(t f ) and a(t f − δ) • a(t f + ∆) -the scale factor at t f + ∆, the t-axis point of tangency for a cir (t) and a ′ 2 (t).The transitional function is just the equation of the circle If a(t f ) and a(t f − δ) were to coincide, we would be introducing another discontinuity into the model, the change from the inflationary slope to the infinite slope of the circular arc.
With the displacement δ, the edge of the circular arc lies earlier on the timeline than t f , and δ sets that duration.We note that a(t f − δ) is therefore never the physical value of the scale factor and so has no direct effect on the expansion of space.We establish smoothness by equating the scale factors and their derivatives at the tangent points.Thus, we have five unknown parameters in the four matching conditions.In section 5.3, we invoke a fifth equation to specify the model completely.

Power Law Transitions
Imposing a circular arc is an unreasonably strict condition to use to define an interpolating function.We shall continue now by exploring more general power-law solutions for a(t) in the transition between the period scale factors.A properly chosen power-law section would provide continuity but would also potentially lengthen one or the other period, depending on its orientation.A power law with n < 1 has essentially no impact on the length of the radiation era.For a power law with n > 1, a very small increase in the inflationary period could have a substantial impact on the scale factor, as discussed further in section 4.2.

Interpolating Power Laws with n < 1
The power law takes the form We use the notation a p (t) for power laws with n < 1.The subscript p denotes the representative parabola for n = 1 2 , which opens to the right and has a horizontal axis.The unknown coefficient D is the analog of the unknown radius of the circular arc.Additional unknown parameters, a p (t f − δ), δ, ∆, and a(t f + ∆), correspond to the parameters displayed in figure 3 for the circular arc.We analyze the continuity of the forms of the scale factor at the two points of tangency.

Smoothness at t f
For the first matching condition-continuity of a(t)-at t f , a p (t) The second matching condition at t f equates the time derivatives, generating an expression for the a-axis vertex coordinate: where we have used the relations in eqs.( 14) and ( 18) to express these in terms of the current temperature T (t 0 ), since t = t 0 also provides the calibration scale for a.With the vertex coordinate from eq. ( 28), the scale factor is 4.1.2Smoothness at t f + ∆ The interpolating transition we have imposed between the end of inflation and the beginning of the radiation era shifts the eq.( 13) scale factor according to where we use the prime to distinguish this shifted expression.The vertex of the radiationera t 1/2 scale factor remains at (t = 0, a = 0).The third matching condition, in which the interpolating power-law equals a ′ 2 (t) at the point of tangency, yields the noninformative solution However, the final smoothness condition equates the derivatives of the scale factors a ′ 2 (t) and a p (t) at t f + ∆, so we have Thus we ultimately arrive at the condition, The formalism leaves us with the need to fix ∆ to evaluate the model; ∆ and a(t f + ∆) are physical but as yet unknown parameters.The others are mathematical constructs with no direct physical meanings.After substituting the early universe parameter values assumed in section 2.1, we find, for example, the solution for n = 1 2 at ∆ = 10 −35 s δ ≈ 2.65 × 10 −6 ∆ (38) (evaluated using Maple).Table 1 lists additional values of δ for a sample set of transition durations ∆.The purpose of having three significant figures listed in the table is to illustrate the relationship between the displacement and any changes to the transition scale.
We have accomplished the objective of parameterizing the transition from the inflationary scale factor with a family of power-law scale factors that ensure sufficient smoothness to have continuity of a and its first time derivative-although we have not imposed the condition of continuity on any higher-order derivatives.The slope of the inflation-era scale factor grows at a rate on the order of the Hubble parameter, and a requirement of continuity on the second derivative would effectively extend inflation into the subsequent period, rather than marking the physical end of inflation as the point at which the second derivative becomes negative.

Power Laws with n > 1
To continue the study of alternative transitions, we now examine power laws with n > 1, containing unknown parameters analogous to those analyzed in the previous section.The scale factor formula is The subscript v denotes the representative inverted parabola for n = 2 with the vertical axis parallel to the a-axis.The displacement of the vertex from a(t f ) now places δ at a time later than the tangent point at t f + ∆, as shown in figure 4.
Repeating the analysis of the matching conditions at t f and t f + ∆ yields the scale factor and a fourth matching condition As table 1 details for a sample set of transition durations, the vertex displacement δ for a power law with n = 2 and ∆ = 10 −35 s is almost six orders of magnitude farther from the point of tangency than that of the power law with n = 1 2 .The power law a p (t) ∝ t 1/2 can establish continuity with the slope of the inflationary scale factor with such a minute displacement, because the power law has an infinite slope at the vertex.However, the Figure 4: The displacement δ necessarily places the vertex of the representative inverted parabola with power law index n = 2 later than the end of inflation at t f and the tangency point at t f + ∆ (not to scale).Increasing the duration over which the the interpolating scale factor applies also shifts the vertex similarly.power law a v (t) ∝ t 2 has no such infinite slope, and the difficulty of establishing continuity with the large slope at the end of inflation, ∼ 10 11 s −1 , manifests itself in the displacement being many orders of magnitude greater than that of a p (t) ∝ t  Increasing the duration ∆ for the power law n = 1 2 tends to set the displacement of the vertex.However, because of the difficulty of establishing continuity with the inflationary slope at t f , the displacement δ for the n = 2 interpolator is many orders of magnitude greater.With its vertex located later on the timeline than t f , the table shows that increasing the transition has the effect of shifting the vertex of a v (t) farther away from the end of inflation.29) and (40), rescaled by a translation of the t-axis t → t ′ = t f + t.The timeline starts at the arbitrarily small initial value t = 10 −43 s.The graphs represent power laws for n < 1 and n > 1 with a representative sample of powers.However, the graphs of the scale factors themselves offer somewhat limited insight into the evaluation of the quality of the interpolating functions.For that, we now look instead to graphs of the Hubble parameter, In the transitional Hubble parameters below, H in is the constant inflationary Hubble parameter, taken to be 10 16 GeV.For the two classes of power laws, we find In figure 6, showing graphs of eqs.( 43) and (44), we note the requirement of smoothness at the beginning of the radiation era causes abrupt shifts downward and upward as t → ∆ for interpolating scale factors not proportional to t 1/2 .We note that the power laws with n > 1 that overlay each other in both figures 5 and 6 must exhibit the shift downward to establish continuity with the radiation era Hubble parameter.Unable to justify a physical Figure 5: The three a v (t) scale factors with power laws t 3/2 , t 2 , and t 5/2 essentially overlay each other, and the a p (t) scale factor proportional to t 1/4 approaches those of the power laws with n > 1.The graph also shows the unsmoothed scale factor defined by eq. ( 13), which exceeds that of a p (t) at t 3/4 .The t-axis timeline begins at t = 10 −43 s after inflation terminates, while the vertical dotted line at t = 10 −35 s marks the nominal start of the radiation era.basis for this behavior, we shall move forward in our analysis by eliminating these power laws as valid interpolating functions and focus on more specifically determining workable interpolating functions with n < 1.We also note that as n → 1  2 from above or below, the scale factor a p (t) transforms more seamlessly into the radiation era.So at this stage, we expect that the most suitable interpolating functions will correspond to the index value n = 1  2 , or something close to that.The power laws for the interpolating region and the subsequent radiation-dominated era are both horizontally-opening parabolas (or nearly so), which differ principally in their vertex placements and radii of curvature.Once again, the graphs based on the scale factors a v (t) for t 3/2 , t 2 , and t 5/2 essentially overlay each other.The discontinuity of the radiation-era Hubble parameter with H in exceeds two orders of magnitude.The graph based on scale factor a p (t) for t 3/4 fails to display asymptotic behavior with H in at small times, because the timeline has the same t-axis translation and does not start at t f .Again the vertical dotted line marks the start of the radiation era at t = 10 −35 s.

The Equation of State
We continue with the evaluation of the usefulness of the possible interpolations by considering a parameter ϵ H , which is an alternative to the equation-of-state parameter ω that satisfies p = ωρ, according to Following a graphical technique described by Kaloian Lozanov [25], we shall interpret the formula as the slope in a plot of the evolution of the scale factor from inflation through reheating, matter domination, and finally the dark-energy-dominated era.Appendix B provides more information about this expression.Figure 7 shows our version of the Lozanov graphical approach for the power law a p (t) ∝ t n , with n ranging from 0.05 to 0.95.Table 2 lists statistics for some of the graphed power laws, as well as smaller and larger values of n.Aside from the footnoted observations in table 2, we note a further curious feature of figure 7. The graphs at the upper and lower extremes of n display almost cusp-like changes of slope at the tangent point between the transition and the radiation-dominated era at t = t f + ∆. Between inflation and the start of the radiation-dominated era, ϵ H changes from a value much greater than 2 to less than 2 for n > 1 2 .Conversely, for n < 1 2 , the parameter starts less than 2 and then becomes greater than 2. With our expectation that ϵ H = 2, only the power law n = 1 2 (the same power law index as in the radiation era itself) appears able to transition seamlessly to the radiation era, which suggests that all powers except n =   46), the expected ϵ H = 2 for a radiation-dominated scale factor occurs at n ≈ 0.5002.(b) Values of ϵ H → 2.6 − and ω ≈ 0.73 signify unphysical, exotic tachyon-like particles with velocities greater than the speed of light, which section 5.2 discusses in detail.(c) For ϵ H ≈ 1.50 and ω ≈ 0.00, we have a transition from inflation to an equation of state that would be consistent with a matter-dominated universe.We take up consideration of the single-component matter-dominated universe in section 5.3.(d) As ϵ H → 1.0 + and ω → − 1  3 , the scale factor remains inflationary, effectively eliminating the transition.
Figure 8 plots ϵ H versus the power law index n at time (t f + ∆): For n > 1 2 , the negative term in the numerator of eq. ( 46) is small relative to the other terms appearing in the fraction, which are on the order of a(t f ) and essentially cancel with the factor of the same characteristic size in the denominator.The linear relation ϵ H − 1 ∝ (1 − n) remains, as the graph and table 2 show.In contrast, for n < 1  2 , the scale factor a(t f ) is small relative to the other terms in the fraction.As n → 0, the fraction and (1 − n) are both increasing, and ϵ H increases to approximately 2.6.
Figure 8 aligns with the possibility raised by figure 5 that only n = 1 2 results in a transition to the radiation era without an awkward, cusp-like feature-whose very presence would seem to be contrary to the dictum we have adopted of modeling the transitions in a smooth fashion.However, precise calculations consistently indicate that a slightly different n ≈ 0.5002 actually produces the seamless transition.In the same way as the abrupt shift that we cannot explain in the graph of the Hubble parameters tends to disqualify all power laws except n = 1 2 , the cusp-like features again appear likely to signify unphysical, unexplained behavior.However, before we attempt to resolve these conflicts, we shall review additional constraints on the equation of state, starting with constraints related to the speed of sound.

Speed of Sound Constraints
Another tool for evaluating the interpolators is the application of constraints on the speed of sound to the equation of state.A speed of sound less than zero or greater than the speed of light would violate stability or causality, respectively [27,28].Stability requires that the speed must be real; imaginary phase speeds would correspond to imaginary frequencies, or modes that grow exponentially with time.At the other end, special relativity imposes the standard limitation that information carried by arbitrary quanta cannot propagate faster than the speed of light c = 1 in vacuum.Thus, we expect the sound speed of the transition waves to obey inequalities We assume, as is standard, that the inflaton wave oscillations are fast and thus adiabatic, so that a passing wave brings about temperature changes without conductive heat transfer.The thermodynamic behavior is reversible, and so the entropy per unit mass is constant [29] as an inflaton wave passes through.Pressure p = p(s, ρ) becomes a function of the density ρ only, p = p(ρ).We also make the assumption that the inflaton condensate at the end of inflation is a perfect fluid, allowing us to apply a linear, single-component equation of state, p = ωρ, expressing the dependence of pressure on density in terms of a ρ-independent equation-of-state parameter ω.This environment yields a sound speed and using ω = 2 3 ϵ H − 1, the condition 0 ≤ ω ≤ 1 implies 3 2 ≤ ϵ H ≤ 3.At the precise end of inflation, ∆ = 0 in eq. ( 46), leaving Figure 9: Values of the parameter ϵ H = 1 + (1−n) Hδ at the end of inflation, t f ≈ 4 × 10 −39 s, for a p (t) ∝ t n with n ≈ 1 2 .The gray band depicts the permissible values of ϵ H , 3 2 ≤ ϵ H ≤ 3, subject to the assumption that the inflaton condensate at the end of inflation is a singlecomponent perfect fluid with equation of state p = ωρ.
Figure 9 displays the effect of enforcing the sound speed restrictions from eq. (47) on ϵ H ; these conditions severely restrict the permissible range of power laws indices.The line plot is a Python cubic spline interpolation of Maple-generated solutions for the parameter ϵ H from eq. (46) in 0.001 increments around n = 1 2 .The section of the spline interpolation within the gray horizontal band contains valid values of ϵ H , corresponding to interpolating theories with stable, causal sound speeds.Reading off the graph, we see the permissible physical power law band for the continuous function transitioning from the end of inflation to the radiation era lies approximately between 0.4990 < n < 0.5005.This narrow band is consistent with the large separation in figure 7 between the power law n = 1 2 and the closely adjacent powers, n = 0.45 and n = 0.55.

Continuity of the Equation of State
Although we have found restrictions on the allowed range of power law values, a tension between the expected and derived equations of state for n = 1 2 actually remains.We do not find that ϵ H = 2 corresponds exactly to n = 1 2 .We instead recall the expected ϵ H = 2 was associated with n ≈ 0.5002 reported in table 2, and we now seek to understand the reason for the slight discrepancies between these values and the derived ϵ H ≈ 2.00039 that corresponds to the exact n = 1 2 .Comparison of the scale factor in eqs.( 28) and ( 29) of section 4.1 with the first two terms in the numerator of the fraction in the expression below, indicates that those two terms arise from the a-axis scale factor displacement a p (t f − δ) of the interpolating function's vertex.The third term in the numerator represents the functional dependence of the scale factor on time.Table 3 separates the difference between the expected and derived ϵ H into the relative contributions from the components of the numerator, ϵ H displacement and ϵ H time .For comparison, we repeat the analysis for a single-component, matter-dominated transition function establishing continuity between the end of inflation and a matter-dominated era (that is, with n = 2 3 power laws).We see the results are qualitatively the same.The displacement of the vertex of the scale factor along the a-axis is responsible for the discrepancies.
2.65 × 10 −41 2 2.00039 0.00039 1.00000 7.71 × 10 −42 1.5 1.50020 0.00020 0.50000 Table 3: These parameters correspond to power laws a p (t) with indices n = 1 2 and n = 2 3 transitioning to radiation-dominated and matter-dominated eras with scale factors likewise proportional to t 1/2 and t 2/3 .The column ϵ H derived reconstructs the parameter as the sum of 1 and the contribution from the displacement and the time components.We conclude that the displacement causes the difference from ϵ H expected .
If not for the contribution of the vertex displacement, we would have seamless transitions of the equation of state between the interpolating power laws and the radiation or matter eras.A first-order phase transition at (t f + ∆) might be responsible for the cusp, but we reason against that possibility.The dynamics of the expansion of space at the tangent point undergoes no change.Prior to and after (t f + ∆), the power law index n governing expansion remains approximately the same for each single-component era.Also, the transition precedes the period of preheating described in section 7 and subsequent thermalization, so that we expect temperature to evolve smoothly at t f + ∆.
Instead, invoking continuity of the equation of state at (t f +∆) and noting ϵ H expected = 1 n for a single-component universe, as appendix C shows, we solve eq. ( 50) for the displacement δ and find Substituting a(t f ) = T (t 0 )/H from eq. ( 18) with H = H in yields Analysis of eq. ( 50) demonstrated that the displacement caused the cusp-like feature.With the eq.( 53) result, and recalling a(t f ) = T (t 0 ) H , we instead have ϵ H displacement = 0, and the bump on the curve is gone.
However, making the assumption of continuity of the equation of state at t f + ∆ destroys the smoothness of the scale factor that we imposed at both t f + ∆ and t f .So we must reexamine the matching condition at t f + ∆; returning to the fourth matching condition and trying to solve for ∆, we see that With δ = n H , this simplifies to So the transition period ∆ is undefined for n = 1 2 , which invalidates the claim of firstderivative smoothness imposed by the eq.( 38) parameters, δ ≈ 2.65 × 10 −6 ∆.For n ≈ 0.5002, associated with ϵ H = 2, the new formula's value of ∆ is in fact less than zero.Since ∆ is supposed to represent the length of time over which the interpolating function applies, this value is manifestly unphysical.
Furthermore, substituting δ = n/H in the interpolating scale factor, eliminates the a-axis displacement of a p (t).We introduced the displacement of the powerlaw vertex in section 4.1.1 in order to enforce the smoothness condition at t f , but this would be undone by the assumption of exact continuity of the equation of state.The Lozanov graphical approach to analyzing the equation of state suggests a range of power law indices, 0.5000 ≲ n ≲ 0.5002, are reasonable, and this is supported by the values the are permitted by the speed-of-sound constraints, 0.4990 ≲ n ≲ 0.5005.Having tried unsuccessfully to establish continuity with the equation of state at t f + ∆, we now seek an explanation of the discontinuity.If the transition represents a continuation and ultimately a termination of inflation, a local discontinuity might result from a weak phase transition of unknown character.A second explanation may be that a power law index not equal to 0.5 in the transition signals that the composition of the universe is not strictly radiation-dominated as the transition ends, and so a single-component model is not sufficient to describe the dynamics.

Summary of Numerical Results
We have found that a transition function with power-law index n can provide seamless first-derivative smoothness over the period between the end of inflation and the development of a single-component n-era universe, while obeying fundamental stability and causality constraints.During the transition, the scale factor increases by approximately an order of magnitude more, compared with what it would have been in a model with a sharp cusp dividing the inflationary from post-inflationary functional forms; and that additional accumulated expansion factor remains as time progresses.Figure 10 shows the key comparisons.
Figure 11 depicts the increases toward asymptotic limits more clearly.Both figures 10 and 11 also reveal that these increases occur primarily in the vicinity of t = 10 −37 s and do not particularly depend on the duration of the transition ∆.Table 4 contains further data, including how much larger, relatively speaking, the universes with the smooth interpolations are than the models without smoothing.The underlying numbers show that  after 10 −34 s, the asymptotic values 9.8 and 11.2 have completely stabilized (to over 12decimal-place precision).Even by 10 −37 s, the increased expansion factors have already grown to be within 2% of their asymptotic values.
Figure 11: The scale factor ratios ap(t) a ′ 2 (t) .At approximately t = 10 −37 s, the ratios reach greater than 98% of the asymptotic values of 9.8 and 11.2 for a p (t) ∝ t 1/2 and a p (t) ∝ t 2/3 , respectively.We are left with the interesting result about what happens when we insert an interpolating function after the end of inflation to smooth out the dynamics.Compared with the models with discontinuous derivatives-signifying abrupt transitions between the inflationary period and a period with a different equation of state-the total expansion of the scale factor is greater by about an order of magnitude (or between 2 and 3 e-folds).In a way, this is unsurprising, since the interpolating function allows the inflationary expansion to tail off a bit more gradually, and so the net result is always a larger universe at later times.This kind of increase in the scale factor will form the basis for our analysis of the effect of continuity in the numerical analysis going forward.

The Smooth Scale Factor in the Preheating Model
Having concluded that enforcing a smooth transition results in an order-of-magnitude increase in the ultimate scale factor of the subsequent single-component universe, we shall now examine the effects of this change on reheating, based on the preheating model of KLS [21], in which the inflaton couples to a second scalar field χ in the era following after inflation, which is taken to be a matter-dominated universe.We shall evaluate preheating effects using a smooth interpolating power law with n = 2 3 , as described previously in section 4.1 as an example of a power law with n < 1.We compare our results to those of the KLS model, which employs the scale factor a(t) ≈ a f (t/t f ) 2/3 with a discontinuous slope.Our numerical analysis shows that the larger scale factor in the smooth model decreases the χ occupation numbers n k and dilutes the total number density n χ .The dilution arises naturally out of the volume increase due to the greater expansion of space-although the broad parametric resonance during preheating partially offsets the effect.Broad parametric resonance involves all modes of the scalar field χ less than a specific maximum being involved in quasi-resonant interactions with the inflaton, and it causes an exponential increase in the number of χ particles created.

Occupation Numbers
In this section and section 7.2, we briefly summarize the foundations of the detailed, extensive case that KLS present in support of their theory.The Lagrange density for the scalar field χ coupled to the inflaton, in expanding space with vanishing mass parameter m χ = 0, generates the equation of motion where k = √ k 2 , for the Fourier mode χ k in momentum space.The inflaton at the end of inflation is a coherently oscillating field of form ϕ(t) = Φ(t) sin(mt), with amplitude envelope Φ(t) = M P √ 3πmt [23], so that In slow roll inflation, chaotic inflation, and other inflationary models in which the friction term 3H φ in the equation of motion (1) becomes negligible, the inflaton exhibits sinusoidal oscillating behavior around ϕ = 0. (Here the argument of the sine function has time t in units of m −1 , which the KLS model uses throughout.)The appearance of the Planck mass M P in Φ(t) derives from the Hubble parameter expressed in terms of the gravitational constant.The units of k are m, and the scale factor, normalized in the Robertson-Walker metric with a(t 0 ) = 1 today, remains dimensionless.Broad parametric resonance consists of non-adiabatic oscillation of the χ field in Fourier-space regions where the equation of motion is unstable.The character of the instability is revealed by converting eq. ( 64) into the standard Mathieu equation.Rescaling the scalar field, eliminates the friction-like term and so yields Now, recasting the argument of the oscillating term by setting z = mt completes the conversion into the Mathieu equation: The prime represents the derivative with respect to the argument z, and the two parameters in the equation are The resonance behavior of solutions to the Mathieu equation depends on the values of these A k and q, which determine the stable and unstable regions.Appendix D reproduces the standard plot depicting the stability and instability regions in the q-A k plane with a graph of the Mathieu equation parameters.
The oscillations of the scalar field exhibit adiabatic instability when and energy transfer occurs between the inflaton and the scalar field χ.Trial solutions of the Mathieu equation, are unstable for real values of the Floquet characteristic exponent µ k [30,31].Section 7.2 discusses µ k in more detail.The mode occupation number n k is the energy of the mode in question, divided by the single-particle energy ω k : (The adjustment −n k 0 to account for the zero-point energy density is effectively negligible.) Figure 12 reproduces the results of the discontinuous scale factor of the KLS model, for the scalar field mode amplitude X k and the exponential increase in the corresponding occupation number log n k .The t-axis timeline of both graphs becomes a count of the number of oscillations of the inflaton after t is expressed in units of 2π m , with which the revised equation of motion (66) is Broad parametric resonance preheating requires certain preconditions on Φ(t) and the Mathieu equation parameter q, and it begins shortly after the end of inflation, after approximately one quarter of an oscillation of the inflaton.(KLS use this approximation to advance their analysis.)With time defined in terms of the number of oscillation cycles, t f = π 2m ≈ 10 −37 s, which makes the timeline consistent with that which we have found for the continuous scale factor; our order-of-magnitude increase in the size of the cosmos also appears at around 10 −37 s.
In graph (b) of figure 12, the scalar field spans many instability bands in the first ∼ 10 oscillations, as q decreases substantially, and the resonances cause exponential growth in the occupation number.From about 12 to 17 oscillations, the growth flattens as q lessens while crossing the stability region corresponding to q values decreasing from about 2 to 1. Broad resonance and growth resume in the next 10 oscillations in the instability band for q ≲ 1 and A k ≈ 1, before ultimately terminating after ∼ 34 oscillations.Appendix D also shows graph (b) of figure 12 superimposed on the final three instability regions of figure 16 (corresponding to decreasing q as time progresses).In figure 13, we repeat the presentation from figure 12 using the smooth transitional scale factor in place of the kinked scale factor of the KLS model.The scalar field and occupation number show sharp decreases from figure 12 to figure 13.We can examine the effect of the continuous scale factor more precisely by analyzing the root mean squares of X k and n k averaged over the 10 oscillations following the end of preheating, which occurs after approximately 34 oscillations.With time t in units of m −1 according to the KLS formalism, we can convert the scale factor units of time in seconds to oscillations, We have also used the assumption for t f that broad parametric preheating begins after inflation ends, at one fourth of an oscillation.Then we apply a factor of 10 for the approximate order-of-magnitude increase in the continuous a(t), For X k (in the k = 4m mode), we find a modest decline of ∼ 0.03 in the root mean square, due to the order-of-magnitude increase in the scale factor.Figure 14 shows log n k for both forms of the scale factor for 10 oscillations following the end of broad resonance.The decrease in log n k because of the effect of the larger scale factor causes a reduction of just ∼ 0.002 in the root mean square of the occupation number n k at 10 oscillations after broad resonance terminates.Local maxima in log n k for the smooth scale factor in figure 14 occur at every half oscillation of ϕ at t = 1 2 , 1, 3 2 , . ... At these points, where sin(2πt) = 0 in eq. ( 72), the frequency reduces to with values less than one, 0.080 ≲ ω k ≲ 0.095.For the 10 oscillation periods under consideration with the smooth scale factor model, this range of fractional frequencies has the effect of increasing the contribution of the term containing the kinetic energy Ẋk in the occupation number, even as it tends to suppress the contribution of the potential-like X k term.Thus, the small fractional frequency generates the local maxima.The range of larger frequencies with the cusped scale factor following the end of resonance, 0.80 ≲ ω k ≲ 0.95, has less of an effect and intersperses some local minima, depending on the relative values of Ẋk and X k at the half-oscillation times.
We are able to provide some understanding of the differences in appearance of log n kthat is, the greater degree of dispersion of the amplitudes above the average occupation number in graph (b) of figure 13 in comparison with figure 12-by examining in detail the effect of the fractional frequency.At oscillation 36, for example, the occupation numbers log n k (36) are approximately 45.3 and 40.2 for the cusped and smooth scale factor models, respectively.The kinetic term in the energy, amplified by the frequency, for the most part determines the occupation number in both models.The average occupation numbers over 4 oscillations from oscillation 34 to 38 are approximately 43.9 and 36.8,respectivelyyielding an increase during this period of ∼ 0.03 with the cusped scale factor and ∼ 0.09 with the smooth model.The lower level of the scalar field in the smooth model and (more importantly) its time derivative moderate what would otherwise be an approximately 10-fold difference in the increases based on the values of ω k alone.Thus, we see the greater dispersion of amplitudes above the average log n k in figure 13.Appendix E contains a table that lists some of the supporting data associated with the behavior around oscillation 36, as well as related graphs.

Number Density
The number density of the scalar field quanta has its basis in the process of broad parametric resonance KLS characterize in their paper as stochastic-that is, random.They show that the variation in the phase θ k of the scalar field χ in the course of semiclassical interactions between the χ-particles and the oscillating inflaton field is very much greater than π, which makes successive phases effectively random.However, this does not mean that the there is no net energy flow from one sector to the other.In fact, a growth in the number of particles between classical scattering events can be as much as 3 times as probable as a decrease, based on the numerical effect of possible values for the phase angle in the recurrence relation governing resonance.KLS also separate preheating into two time periods.The first period precedes all backreaction and rescattering, and the second period involves the effect of those interactions on number density, which can be significant.Backreaction and rescattering are quantum effects in which the created χ-particles interact with the background inflaton field.In backreaction, interactions can alter the effective masses of the particles and the frequency of the inflaton oscillations.Rescattering involves a created particle scattering again, either off an inflaton or another χ-particle.However, KLS conclude that the duration of the second period is so brief that during it they can safely neglect the expansion of the universe, and their analysis of that part does not depend on the scale factor.Therefore, here we shall determine the effect of the continuous scale factor on number density conversely without including backreaction and rescattering.
Semiclassical scattering leading to quantum-mechanical χ-particle production involves the interaction of the scalar field χ and the background inflaton field oscillating around zero.KLS derive the number density of the χ field from the adiabatic approximation solution to eq. ( 66), with the scalar field phase θ j k = t j dt ω k and t j representing the time at end of the j th oscillation-such that as time t → t j , the inflaton field is oscillating around its minimum, ϕ → 0. The functions α k (t) and β k (t) are time-dependent Bogoliubov transformation coefficients [32].Around ϕ ≈ 0, eq.(66) becomes The scalar field χ with an effectively random phase θ j k completes a half-oscillation at time t → t j for j = 1, 2, 3, . ... As t → t j for each half-oscillation of χ, the inflaton field concurrently oscillates near zero, creating a period of non-adiabatic energy transfer, which leads to exponential growth in the number of χ-quanta according to eq. (69).At other times, the number density n χ remains stable.Introduction of parameters τ = k * (t − t j ) and κ = k ak * (79) recasts eq. ( 78) as a differential equation with a parabolic cylinder function solution, which is also the Schrödinger equation with an unstable quadratic potential, V (ϕ) ∝ −τ 2 .Appendix F derives the largest mode to participate in the broad parametric resonance, k * = √ gmΦ.The scattering of solutions X k of eq. ( 66) leads to a recurrence relation for the Bogoliubov coefficients, which may be represented by transfer matrix, KLS provide the reflection R k and transmission D k amplitudes from the solutions of the parabolic cylinder equation and also the phase angle φ k , which is a complicated function of the parameter κ, With these, the recurrence relation becomes Noting that the occupation number n k in eq. ( 71) just depends on the Bogoliubov coefficient β k [33], and that for for a coherent process n k ≫ 1, leads to the recurrence relation with the accumulated phase Because the variation in the phases θ j k is very much greater than π, the randomness of θ j kand by extension the randomness of α j k and β j k as functions of θ j k -make θ j tot stochastic.Noting that resonance begins to be suppressed unless πκ 2 ≲ 1, KLS find that for πκ 2 ≪ 1, a growth in the number of particles is three times as likely as a decrease.Within the range 0 < θ j tot ≤ 2π, values of 0 < θ j tot < π 4 and 3π 4 < θ j tot ≤ 2π cause an increase in the number of particles according to eq. ( 85); only over one quarter of the possible range of phases, π 4 < θ j tot ≤ 3π 4 , does the number of χ-particles decrease, as energy flows (incoherently) back to the inflaton field.A second recurrence relation also obtainable [30] from the Mathieu equation (67), in combination with eq. ( 85), yields the Floquet characteristic exponent, Integration of n k for all modes that participate in broad parametric resonance gives rise to the total number density of χ-quanta, The units of number density n χ are the expected m −3 , since occupation number n k (t) is dimensionless.KLS evaluate the integral on the far-right-hand side of eq. ( 89) by the steepest descent method and estimate the number density to be They also determine the maximum Floquet characteristic exponent µ associated with an unknown maximum k max , estimated as k max ≈ k * 2 .We use the proportionality to perform a numerical analysis of the effect of the continuous scale factor by examining the ratio The terms n χ a(t) smooth and n χ a(t)cusp represent the number densities of the smooth and cusped scale factor models, respectively.We anticipate a decrease in the number density due to the increase in volume, moderated to a certain amount by the dependence of the proportionality in eq. ( 91) on µ.The use of the proportionality eliminates the dependence on the unknown mode k * , which KLS estimate as k * (t) ≈ gmΦ(t), as detailed in Appendix F.
In the absence of e 2µmt √ µmt , the greater time allowed for the expansion of space in the smooth model would on its own cause dilution-that is, a decrease in the number density.The order-of-magnitude increase in the smooth scale factor alone would reduce the number density by the cube of the scale factor increase, ∼ 10 −3 .However, the effect of the broad parametric resonance in preheating-in particular, the term e 2µmt √ µmt > 1-may modestly offset the mere increase in the volume of space.The extent of the offset is dependent on the stochastic sin(θ j tot ) in eq.(88).Figure 15 displays the ratio R χ of number density of the smooth scale factor to the discontinuous scale factor as a function of time (again expressed as the number of oscillations).The value of the ratio at the start of preheating, ∼ 10 −3 , Figure 15: The ratio R χ of n χ for the smooth scale factor, to that with a(t) with a discontinuous derivative at the end of inflation, as used by KLS.The increase in volume in the smooth model, resulting from the extra time given for space to expand as inflation tails off, dilutes n χ by ∼ 10 −3 , which the broad parametric resonance term e 2µmt / √ µmt partially tends to offset.Maximizing the offset with total phase sin(θ j tot ) = −1 in eq. ( 88) minimizes the dilution, and that is what is shown in this figure.So with the smooth scale factor, n χ should be diluted by at least the ratios shown here.
reflects the effect only of the expansion of space.As preheating progresses, however, R χ rises to a level slightly greater than 1.7 × 10 −3 at the end of broad parametric resonance, at around 34 oscillations, in the limiting case in which sin(θ j tot ) is consistently equal to −1.In contrast, as sin(θ j tot ) increases toward 1, R χ decreases.For example, at sin(θ j tot ) = 0, R χ ∼ 1.6 × 10 −3 , and R χ is about 1.3 × 10 −3 at sin(θ j tot ) = 0.65.With a slightly larger stochastic value, n χ a(t)cusp is not directly calculable via this method at lower oscillations, and with a stochastic phase of 0.8, the calculations of both n χ a(t)cusp and n χ a(t) smooth , even at the end of 34 oscillations, because that would require the Floquet index µ in eq. ( 88) to be negative.The negative Floquet index signals an essentially unphysical solution, which the model formalism does not support; physically this scenario would describe a net energy flowing back into the inflaton field, while mathematically the formalism breaks down because the saddle point integration method is no longer usable.Thus, examination of the ratio R χ places a bound on the effect of the continuous scale factor.The reduction of the number density due to the expansion of space alone, ∼ 10 −3 , increases only slightly, by at most about 1.7 × 10 −3 after preheating, depending on the values of the stochastic sin(θ j tot ) angles.

Conclusion
This work has explored the consequences of applying the reasonable expectation of smoothness to the physical expansion of space, as expressed by the characteristic scale factors defining the early universe evolving through its generally-accepted, broadly-defined epochs.We focused on the nearly instantaneous slice of time separating the inflationary era and the subsequent era in which the stress-energy tensor was assumed to be dominated by a single component, either radiation or matter.We focused on the transition out of inflation specifically because it is where we inevitably expect to find the sharpest change in the behavior of the scale factor; assuming some realistic values for primordial parameters reveals that the time derivative of the scale factor can decrease by a factor of 1 120 between inflation and the radiation era.Rather than being guided by a specific equation of state model, we imposed a first-derivative smoothness requirement upon the scale factor and looked at phenomenalistic interpolating functions that could connect the inflation and subsequent eras.The assumption of a continuously, steadily declining (but not contracting) slope after the end of inflation led to an in-depth examination of families of interpolating candidates with shifted power-law dependencies on time.We imposed the same requirements of smoothness at the beginning and at the end of the brief interpolating transition period.
From these matching conditions, we uncovered that it was necessary to place the vertices of power law interpolating functions with indices n < 1 prior to the end of inflation at t f and the vertices of functions with n > 1 subsequent to t f , with the displacement in either case parameterized by δ.Also initially unknown was the transition period ∆-the duration of the period between the end of inflation and the single-component universe (whether modeled as composed of radiation or matter).However, implicit in our transition model was a remaining uncertainty in the parameters of the model.We cannot find specific expressions for all of them without imposing additional conditions, and we can do no better than finally expressing the displacement δ in terms of the transition period ∆.Graphical analysis of the Hubble parameter and the equation-of-state and speedof-sound stability and causality constraints allowed us to identify physically reasonable interpolating power-law functions as those having indices that approached the powerlaw indices 1  2 and 2 3 for radiation and matter single-component universes, respectively.Numerical analyses demonstrate the remarkable result that the actual transition lasts approximately 10 −37 s, essentially regardless of the composition of the single-component universe that follows the transition and the duration ∆.In addition, the universe enters the single-component era about an order of magnitude (or approximately 2-3 e-folds) larger than it would have been if subject to a scale factor with a discontinuous slope, which switched instantaneously to t 1/2 or t 2/3 behavior at the end of the inflationary epoch.Although the form of the interpolating function is not exponential, the increase in the lifespan of the universe, 10 −37 s, is not inconsequential compared to the assumption for the inflationary expansion of the universe, N ≈ 60 e-folds.We understand the outcome to be a universe given an additional short sliver of time in which to grow larger simply because we have imposed a condition of smoothness on the physical expansion of space.The numerical analysis adds precision to this result.For a radiation-dominated era following the transition, at 10 −37 s the increase in the size of the universe has attained 98.5% of its asymptotic value, and the corresponding figure for a subsequent matter-dominated era is 98.3%.Generalizing the approach we have used to a multiple-component universe would also be interesting, as would considering the high-scale physics of inflation that might provide the friction needed to end inflation in a smoother way.
We proceeded to examine the effect of the theoretical changes we had described to the dynamic expansion of space (characterized by a smooth scale factor and the resulting predicted increase in the size of the universe) on a subsequent preheating era.The evolution of the universe after inflation remains highly speculative because of the challenges implicit in experimental confirmation.A period of reheating appears to be required in order to be consistent with the later stages of cosmological development, but the details of the reheating dynamics can depend sensitively on the nature of the particle species available to be excited-including as-yet unobserved high-mass species that may not be accessible at standard model scales but could nonetheless have been active participants in the dynamics of the hot, dense early universe.However, we have also discussed the intricate, highly-technical theory of preheating developed by Kofman, Linde, and Starobinsky to address some generic problems with reheating.We applied the KLS formalism to our model with a smooth interpolating scale factor leading into a matter-dominated universe, in order to gauge the effects of the smoothing on the most sensitive χ-particle occupation number n k and the corresponding number density n χ .We were able to estimate the numerical changes compared with the results obtained using the standard cusped scale factor, and we concluded that the differences are not necessarily numerically significant, apart from a dilution in the total particle density that should be common to all models that predict somewhat larger universes after the end of inflation.Specifically, for the occupation number of the most aggressively growing χ mode, we find a modest decline in log n k = 4m of ∼ 2 × 10 −3 in the root mean square for 10 oscillations following the end of broad parametric resonance, which is a consequence of a decrease of just ∼ 3 × 10 −2 in the root mean square of the scalar field χ k over the same period.In addition, by constructing a relation consisting of the ratio of the number density in the cosmology with the smooth scale factor to that with the cusped scale factor, we determine a partial offset to the expected dilution of the quantity of bosons produced by broad parametric resonance due to the approximate 10 3 increase in the unit volume of space caused by the larger smooth scale factor.The stochastic nature of broad parametric resonance precludes a specific prediction, but we find an additional modest increase in the proportion, with an upper bound of ≲ 1.7 × 10 −3 .
It may be somewhat surprising that the effect of a proposed smoothing of the scale factor is so minor-mostly limited to the natural rarefaction of the χ particles that comes with a spatially larger universe.Regarding the possibility (in a case of optimal phase alignment) of, at most, an additional near doubling per unit volume of number density, we note that a doubling of a small number of something in a unit volume may easily be thought of as not negligible.However, in terms of the many, many orders of magnitude of primal particles in a unit volume of early space, we consider the outcome of having, at most, close to twice as many as not of substance.Thus, we view the result of the numerical analysis of the effect of a not insignificant increase in the size of the universe to represent confirmation of the comparative invariance of the KLS preheating model to these kinds of modifications.We are satisfied that result should represent a modestly useful contribution to the body of work in support of this iconic theory.

Appendices Appendix A
In this appendix, we review the derivation of the inverse relation between the radiation temperature of the universe and the scale factor, following the approach outlined by Ryden [1].In an isothermal environment, the first law of thermodynamics dQ = dE + p dV After substituting the pressure of a relativistic gas, p γ = ργ 3 , and the CMB black body energy density, ρ γ = 4σT 4 , we get In an expanding universe with volume element V = a 3 (t)L 3 , the relation becomes which is an elementary separable differential equation, satisfied for Therefore, There are some interesting and noteworthy differences between the behavior in the three broad-resonance growth regions for n k (at least for the particular, rapidly-growing value k = 4m we have selected).There are small oscillations visible, in addition to the secular growth in log n k .During periods when the parameters make the Mathieu equation stable (the white bands in figure 17), the oscillations are comparatively chaotic; this is also what was seen in figure 14 after the last resonant growth period has ending.There is a certain amount of approximately periodic behavior, due to the driving by the amplitude squared of the inflaton field, so there are fairly stark features every half an inflaton oscillation period.However, underneath these is a chaotically varying baseline.During the periods of resonance (the gray bands), the baseline behavior is different, with approximately exponential growth in the occupation number, as is typical in an unstable driven system.On top of this are additional oscillations, qualitatively similar in some ways to those in the stability regions.However, there are also clear manifestations of the nonlinearity of the Mathieu equation, in the form of period doubling or tripling.When the exponential growth is subtracted, the residual still has, on average, one peak per half oscillation of the inflaton field.However, these peaks are not evenly placed or of equal amplitude.During the second shown resonance region, the oscillating residuals have periods equal to the full inflaton oscillation period-a period doubling phenomenon.Within each full oscillation are two dissimilar up-and-down cycles.Moreover, in the vicinity of and during the first, shortest resonant period there is period tripling, with the periodic residuals taking one and half inflaton oscillation cycles to return fully to their original phase space positions.

Appendix
Taking k → 0 in eq. ( 112) generates an expression for the inflaton associated with the minimum-range of mode k: g 2 ϕmΦ Figure 20 shows a standard graphical representation of the bands of ϕ(t) associated with the minimum and maximum ranges of k.We note that k 2 max applies to a band of ϕ(t) for which |ϕ| ≤ 2ϕ * .Thus, we find k * = gmΦ. (119)

Figure 1 :
Figure 1: The contrast between the potentials of slow roll inflation, shown in graph (a), and chaotic inflation, in graph (b).Slow roll inflation requires a plateau to generate enough e-folds of inflationary expansion to solve the Horizon, Flatness, and Monopole Problems.In chaotic inflation, the effects of the friction term in the equation of motion replace that of the plateau in keeping the inflaton from moving to the true vacuum too quickly.After inflation, both models involve the inflaton oscillating around a minimum potential during a period of reheating, which we review in section 1.2.

Figure 2 :
Figure 2: (a) The graph depicts the discontinuity between inflationary and radiation era scale factors at the end of inflation.(b)The circular arc defined by the transitional scale factor a cir (t) intersects tangentially with the inflationary scale factor a 1 (t) and the nowdisplaced radiation era scale factor noted with a prime, a ′ 2 (t).The parameter ∆ is the time period from the end of inflation until the time of continuity between a cir (t) and a ′ 2 (t).The a functions are not to scale.

Figure 3 :
Figure 3: Five unknown parameters characterize the two points of tangency of a cir (t) with the period scale factors (not to scale).

Figure 5
Figure 5 displays the transitional scale factors of eqs.(29) and (40), rescaled by a translation of the t-axis t → t ′ = t f + t.The timeline starts at the arbitrarily small initial value t = 10 −43 s.The graphs represent power laws for n < 1 and n > 1 with a representative sample of powers.However, the graphs of the scale factors themselves offer somewhat limited insight into the evaluation of the quality of the interpolating functions.For that, we now look instead to graphs of the Hubble parameter,

Figure 6 :
Figure 6: The Hubble parameters corresponding to the scale factors shown in figure 5.Once again, the graphs based on the scale factors a v (t) for t 3/2 , t 2 , and t 5/2 essentially overlay each other.The discontinuity of the radiation-era Hubble parameter with H in exceeds two orders of magnitude.The graph based on scale factor a p (t) for t 3/4 fails to display asymptotic behavior with H in at small times, because the timeline has the same t-axis translation and does not start at t f .Again the vertical dotted line marks the start of the radiation era at t = 10 −35 s.

Figure 7 :
Figure 7: The slopes of the graphs equal the parameter ϵ H = d log(H −1 ) d log a .

Figure 8 :
Figure 8: The parameter ϵ H = − Ḣ H 2 with power law a p (t) ∝ t n for n < 1 at the start of the radiation era, with the transition period ∆ = 10 −35 s.

Figure 10 :
Figure 10: The growth in scale factors for single-component universes with smoothness enforced at t f + ∆ with ∆ = 10 −35 s and 10 −22 s in (a) and (b), respectively.We note that the approximate order-of-magnitude increases in the power law scale factors occur at around 10 −37 s in all cases.

Figure 12 :
Figure12:The scalar field and occupation number for the first 60 oscillations, in the model with a discontinuous scale factor, and for the inflaton mass m = 10 −6 M P .The taxis is in units of the number of oscillations, 2π/m.We have selected the specific mode of the KLS model with the wave number k = 4m to maximize the growth of the occupation number.To reproduce the broad-resonance exponential growth, we have used parameters g = 6.25 × 10 −4 , Ẋ(t f ) = 0.045, and Ẍ(t f ) ≈ 0; these were identified empirically, and varying the parameter values away from these will decrease the observable resonance effect.The scalar field derivative Ẋ(t f ) approximates what KLS advise-namely, that the positive-frequency solution X k (t) ≈ exp −iω k t/ √ 2ω k be applied as an initial condition.

Figure 13 :
Figure 13: The scalar field and occupation number for the first 60 oscillations in a model with the continuous scale factor.The graphs show reduced values of X k and log n k compared to figure 12 because of the effect of the order-of-magnitude increase in a(t).

Figure 14 :
Figure 14: These plots compare 10 oscillations of log n k after the end of broad resonance at ∼ 34 oscillations for the two functional forms of the scale factor.

Figure 18 :
Figure 18: Plots of scalar field X and its time derivative Ẋ with the cusped scale factor for 10 oscillations following the end of broad parametric resonance.

Table 1 :
The vertex displacements δ for power laws a(t) ∝ t 1/2 and t 2 for a sample set of transition durations ∆ between the end of inflation and the beginning of the radiation era.
12 result in a cusp in the evolution of ϵ H .This motivates a closer inspection.

Table 2 :
The displacement δ, parameter ϵ H , and equation-of-state parameter ω = 2 3 ϵ H − 1 for power laws a p (t) ∝ t n , with values of n ranging between 0 and 1.(a) Computation sets this value more precisely at ≈ 2.00039.After interpolation using the linear relation associated with eq. (

Table 4 :
The ratios of a p (t) to a ′ 2 (t) at different times during the interpolation period.The percentages represent the degree to which the ratios have approached the asymptotic values reached at 10 −34 s.

Table 5 :
E Data in support of the differences in appearance between graphs (b) of figures 12 and 13.The last column represents the increase in occupation number log n k (36) compared to the average value over 4 oscillations from oscillation 34 to 38, which are 43.9 and 36.8 for the KLS and smooth scale factor models, respectively.