Polytropic equation of state and primordial quantum fluctuations

We study the primordial Universe in a cosmological model where inflation is driven by a fluid with a polytropic equation of state $p = \alpha\rho + k\rho^{1 + 1/n}$. We calculate the dynamics of the scalar factor and build a Universe with constant density at the origin. We also find the equivalent scalar field that could create such equation of state and calculate the corresponding slow-roll parameters. We calculate the scalar perturbations, the scalar power spectrum and the spectral index.


Introduction
Recently the Bose-Einstein condensate (BEC) dark matter (DM), a cosmological model based on condensate state physics, has appeared in several works as an attempt to explain the origin and nature of DM [1,2,3,4,5,6,8] and was initially used to describe DM halos. The equation of state of the BEC can be found from the Gross-Pitaevskii (GP) equation [1,2] and is given by where l s is the scattering length, = h/2π is Planck's constant, ρ is the density distribution of the single component BEC DM and m is the mass of particles that have been condensed. Assuming the hypothesis that the cold dark matter in a galaxy is in form of BEC the density distribution of the static gravitationally bounded single component BEC DM is given by ρ(r) = ρ * sin kr/kr, where ρ * = ρ(0) is the density in the center of the condensate and k is a constant. Giving the conditions ρ(R) = 0 and kR = π, where R is the condensate radius, the condensate DM halo radius can be fixed as R = π( 2 l s /Gm 3 ) 1/2 . The calculated total mass of the condensate DM halo is M = 4π 2 ( 2 l s /Gm 3 ) 3/2 ρ * = 4R 3 ρ * /π. So the mass of the particles of the condensate is [2] m = π 2 2 l s GR 2 1/3 ≈ 6.73 × 10 −2 l s 1 fm The Bose-Einstein condensation process, that is a very well observed phenomenon in terrestrial experiments, occurs when a gas of bosons is cooled at very low temperatures, near absolute zero, what makes a large fraction of the particles occupy the same ground state. The BEC model can also be applied to cosmology in order to describe the evolution of the recent Universe. In this attempts it can be assumed that this kind of condensation could have occurred at some moment during the cosmic history of the Universe. The cosmic BEC mechanism was broadly discussed in [3,4]. In general the BEC takes place when the gas temperature is below the critical temperature T crt < 2π 2 n 2/3 /mk B , where n is the particles density, m is the particle mass and k B is the Boltzmann's constant. Since in an adiabatic process a matter dominated Universe behaves as ρ ∝ T 3/2 the cosmic dynamics has the same temperature dependence. Hence we will have the critical temperature at present T crt = 0.0027 K if the boson temperature was equal to the radiation temperature at the red-shift z = 1000. During the cosmic adiabatic evolution the ratio of the photon temperature and the matter temperature evolves as T r /T m ∝ a, where a is the scale factor of the Universe. Using as value for the present energy density of the Universe ρ = 9.44 × 10 −30 g/cm 3 BEC will happens if the boson mass satisfies m < 1.87 eV.
Recently the cosmological process of the condensation of DM was investigated [5,6] and in this model it is assumed that the condensation process is a phase transition that occurs at some time during the history of the Universe. In this case the normal bosonic DM cools below the critical condensation temperature, that turns to be memorable to form the condensate in which all particles occupy the same ground state. In this new period the two phases coexist for some time until all ordinary DM is converted into the condensed form, when the transition ends. The time evolution of cosmological parameters as the energy density, temperature, scale factor and both scalar and tensorial perturbations is changed during the phase transition process.
There are other cosmological scenarios with BEC. In [7], the author consider a Post-Newtonian cosmological approach to study the global cosmological evolution of gravitationally self-bound BEC dark matter and the evolution of the small cosmological perturbations. The author verifies that the presence of BEC dark matter changes the cosmological dynamics of the Universe. In [8] the author do the same as in [7] with the generalized EoS p = (αρ + kρ 2 ) c 2 . Optimal parameters in good agreement with the ΛCDM model are found. A natural justification and generalization of the adhesion model, the Burges equation and the cosmological Kardar-Parisi-Zhang equation, that describes the large-scale structure of the Universe but it is introduced heuristically, are studied in [9]. The possibility that a significant part of the compact astrophysical objects are made of BEC is also considered in [10]. The power spectrum for the BEC dark matter is studied in [11], where the authors limit the mass of the dark matter particle in the range 15 meV < m < 700 meV. The effects of a finite dark matter temperature on the cosmological evolution of the BEC dark mater systems is verified in [12]. In this paper the authors have shown that the presence of thermal excitations leads to an overall increase in the expansion rate of the Universe. With the aim of understanding whether the gravitational interactions of axions can generate entropy, the behavior of the axions during non-linear galaxy evolution is studied in [13], but the assumption that the axion field can form a Bose-Einstein condensate is not yet confirmed. In [14] the authors assume that the dark matter particles are described by a spin-0 scalar field, called scalar field dark matter. These bosonic particles are ultralight, with masses down to the order of 10 −33 eV/c 2 . An ultralight phase-space density is suggests and there is the possibility of formation of a Bose-Einstein condensate. They shown that the scalar field dark matter is compatible with observations of the cosmic microwave background and the abundance of the light elements produced by Big Bang nucleosynthesis. With a relativistic version of the Gross-Pitaevski equation, the authors in [3,15,16] propose i) a novel mechanism of inflation, ii) a natural solution for the cosmic coincidence problem with the transition from dark energy into dark matter, iii) a very early formation of highly non-linear objects like black holes and iv) log-z periodicity in the subsequent BEC collapsing time.
We can generalize the BEC equation o state (EoS) (1) [15,16] as follows where k > 0 represent repulsive and k < 0 attractive self-interaction and the linear term describe the well known radiation (α = 1/3), dust matter (α = 0) and cosmological constant (α = −1), and the less known stiff matter (α = 1). The stiff matter model is a specific cosmological model where the matter content of the Universe has an equation of state of the form, p = αρ, with α = 1, where ρ and p are, respectively, the fluid energy density and pressure [17]. This model can also be described by a massless free scalar field. The energy density of the stiff matter is proportional to 1/a(t) 6 and this result indicates that there may have existed a phase earlier than that of radiation, where α = 1/3 and ρ ∝ 1/a(t) 4 , and after inflation in our Universe, which was dominated by stiff matter. This peculiarity motivated us to investigate their behavior in the analyses made here in this work and to consider the implications of the presence of a stiff matter perfect fluid in FRW cosmological models. The EoS p = αρ + kρ 2 is the sum of the linear term and a quadratic term, that describes BECs. At late times, when the density is low, the BECs contribution to the EoS is negligible and the evolution is determined by the linear term. But in the early Universe, when the density is high, the term due to BECs in the EoS is dominant and modifies the dynamics of the Universe. Lately this model was used as a model of the early Universe. We can assume that this EoS holds before radiation era and for the repulsive self-interaction the Universe starts at t = 0 at a singularity with infinite density but finite radius. For the case of attractive self-interaction the Universe has always existed and for the non-physical limit t → −∞ the density tends to a constant value and the radius goes to zero, both exponentially [15,16].
In this letter we study the generalized EoS [18] to describe the physical state of the matter content of the Universe, where n = 1 and α = 0 describes cosmological BECs. With the generalized EoS this model can present a phase of early accelerated expansion. It can be also used to describe a phase of late accelerated expansion, depending on the choice of the parameters [15,16]. We calculate the primordial cosmological dynamics in this model for α = 1/3, α = 0 and α = 1 [15]. We also find a scalar potential that can generate this EoS and calculate both scalar and tensorial perturbations [16]. We study the corresponding slow-roll parameters and the power spectrum and spectral index are calculated. To motivate the model studied here we can see an analogy between this polytropic equation of state and a cosmological model where the fluid that fills the Universe has an effective bulk viscosity [19]. If we write p = αρ − 3Hη, where η is the viscous coefficient, we have exactly the generalized EoS, when η ∝ ρ and H ∝ ρ 1/n .
The present letter is organized as follows. In section 2 we introduce the generalized Gross-Pitaevskii equation used to describe the BEC in the short-ranged scale. In section 3 we study the evolution of the Universe filled with fluid described by a polytropic EoS. For the case of a non-singular inflationary Universe we find the scalar potential that could generate the polytropic EoS and calculate the slow-roll parameters. The primordial fluctuations, such as gravitational waves, perturbations for the gravitational potential and for the density contrast are calculate, and we find the quantum power spectrum and spectral index in section 4. We present our conclusions and discussions in section 5.

The generalized Gross-Pitaevskii equation
The Gross-Pitaevskii (GP) equation is a long-wavelength theory widely used to describe dilute BEC, but it fails [20] in the case of short-ranged repulsive interactions in low dimensions. Therefore the interparticle interaction term in the GP equation must be modified and in this model the ground state features of the BEC is described by the generalized GP equation [1,5] ı where φ(t, r) is the wave function of the condensate, m is the particles mass, V is the gravitational potential that satisfies the Poisson's equation ∇ 2 V ( r) = 4πGρ, g = dg/dn, n = |φ(t, r)| 2 is the BEC density and ρ = mn. To understand the physical properties of a BEC we can use the Madelung representation of the wave function [1,3,4], which is where S(t, r) has the dimension of an action. This transformation will make the generalized GP equation (5) breaks into two equations where V Q = −( 2 /2m)∇ 2 √ ρ/ρ is a quantum potential and v = ∇S/m is the velocity of the quantum fluid. The effective pressure of the condensate [3,4,5] is given by Now writing g ∝ ρ γ we find the generalized EoS where k is a proportionality constant that will be determined in the context of our model, or can be related to the mass and the scattering length of the boson in the case of the long-wavelength theory, and γ ≡ 1 + 1/n is the polytropic index.

Equation of state and cosmic dynamics
Following the present data [21,22,23] we assume a flat homogeneous and isotropic Universe, whose geometry is described by the Friedmann-Robertson-Walker metric, given by where a(t) is the scale factor of the Universe that describes the cosmic evolution, t is the cosmic time and we made the speed of light c = 1. The gravitational dynamics is given by the Einstein's field equations We also consider the Universe filled by a perfect fluid, described by the energy-momentum tensor where ρ is the density of the fluid, p is the pressure and g µν is the metric tensor. This perfect fluid has a general equation of state (EoS), presented in [15,16], that is a sum of a standard linear EoS and a polytropic term, where −1 ≤ α ≤ 1, k is the polytropic constant and 1 + 1/n is the polytropic index. In the linear term α = −1 represents vacuum energy, α = 1/3 is radiation, α = 0 is pressureless matter and α = 1 is stiff matter. The polytropic term may represent self-gravitating Bose-Einstein condensate (BEC) with repulsive (k > 0) or attractive (k < 0) self-interaction, where n = 1 corresponds to the standard BEC.
Here we will consider the high density case, (1 + α + kρ 1/n ) ≥ 0 and n > 0 to describe the primordial Universe which means that the density decreases with the radius. The case 1 + α + ρ 1/n ≤ 0 represents a phantom Universe, where the density increases with the radius [18]. In both cases the polytropic term in the EoS (14) dominates when the density is high and n > 0 and when the density is low and n < 0.
For the EoS (14) the energy conservation equation iṡ where overdot denote the derivative with respect to the cosmic time t and H =ȧ/a is the Hubble parameter. With α = −1 this equation is integrated to give with the minus sign corresponding to k > 0, the plus sign corresponding to k < 0, (a/a 0 ) = a/a 0 , where a 0 is a constant of integration, and ρ * = [(1 + α) /|k|] n . For the repulsive self-interaction (k > 0) the density is defined only for a 0 < a < ∞, where In the case of an attractive self-interaction (k < 0) the density is defined for 0 < a < ∞, and with, in the same limits, p = −ρ * and p → 0.
As we are analyzing the primordial Universe we can see that, when a → 0, the density ρ → ρ * , and we can identify the Planck density ρ * = ρ P = c 5 /G 2 ≈ 5.16 × 10 99 g/m 3 , where ρ P corresponds to a maximum value for the density in the limit a → 0. Likewise, the constant of integration a 0 can be considered a reference point for the transition between the vacuum energy era, when the polytropic component dominates the EoS (14), and the era dominated by linear term in (14). When a a 0 , the scale factor increases exponentially as where t P ≈ 5.39 × 10 −44 s is the Planck time and l P = c t P = (G /c 3 ) 1/2 ≈ 1.62 × 10 −35 m is the Planck length. From the mathematical point of view, in this scenario there is no primordial singularity since this Universe exists at any time in the past (a → 0 and ρ → ρ P for t → −∞). It is obvious that in the limit a → 0 a quantum theory of gravity is required. Even so, exponential solution (19) provides a semi-classical description of the early Universe (see [24] and (3.1) section for more detail). We can also write the equation (14) as where the effective EoS parameter ω(t) is With equation (20) we can calculate the sound speed in the fluid, which is Once again we can find the limits for both repulsive and attractive self-interaction. For k > 0 we have and for k < 0

Non-singular inflationary Universe
We assume that the Universe is filled by the fluid with EoS (14), with n > 0 and k < 0. With the metric (11), the Einstein's field equations (12) and the equation (14) we find the Friedmann equation For small values of scale factor a, i.e., when a a 0 we have (a/a 0 ) → 0 and we can expand the Friedmann equation (25), for x ≡ (a/a 0 ) 3(1+α)/n 1, as We impose the condition 3(1 + α)/n ≥ 1, which means that and we keep only the null order terms to find where H * = 8πGρ * 3 . This means that under these conditions the Universe is inflationary and the singularity can be found at the non-physical limit t → −∞ with a nearly constant finite density. This indicate that the Universe can start at any time t * , which we will define as t * = 0.
For the case of n > 0 and k < 0 and a a 0 the fluid with EoS (14) behaves like the vacuum energy, with constant density. The value of ρ * defines a maximum value for the density and it can be limited by the value of the Hubble parameter at the end of inflation [25]. With ρ * we fix [15] The Friedmann equation isȧ and we can calculate the scale factor for a a 0 a ≈ e tH , where H ≈ 8πG 3 ρ * . After the inflationary stage, for a a 0 , the linear term of the EoS (14) dominates, which means that at some point the inflationary stage will come to an end, and the density will behave as stiff matter, radiation or dust, depending on the value of α. If α = 1/3, for example, the Universe undergoes a radiation era, which would correspond to the standard Universe model. Here we are interested only in the inflationary phase, so we will deal only with solution (31), which is valid for any α and n, since condition (27) is true and k < 0, α = −1 and n > 0.
In case that the linear term has a proportionality constant with value α = −1 we can solve the conservation equation (15) in order to find a density ρ that decays with the scale factor a only if k > 0 and n > 0, and it can be found [15] that where ρ * ≡ (n/3k) n . In this case the density ρ(a) is defined for a ≥ a 0 , and we can't reproduce the behavior of an inflationary Universe. We will not study the perturbations for the case α = −1.

Slow-roll formalism
Here we will represent our fluid as a scalar-field φ and we find the scalar potential V (φ) [26,27] that generates the EoS (14). The scalar field representation can more conveniently retain the features we could expect from fluids with negative pressure, responsible for the inflationary phase, mainly for those that are interesting for cosmology, as the scenarios resulting from phase transitions [28]. The scalar field must obey the Klein-Gordon equationφ where V ,φ = dV /dφ, and we define ρ =φ p =φ Inflation will only occur [26,27] ifφ We combine equations (34) and (35) to finḋ In order to find how the scale factor a varies with the scalar field φ we use the chain rule and combine equation (37) with the Friedmann equation (30) to have With help the equation (21) we can invert the above solution to give us the solution (a/a 0 ) 3(1+α)/n = sinh 2 (ψ) , with ψ defined as With equations (38) and (40) combined we have [16] V (φ) = ρ * 2 1 − α 1 + (a/a 0 ) 3(1+α)/n n + The inflationary expansion of the Universe will happen while the ψ 1. If we expand the scalar potential (42) for small values of the scalar field ψ we get which resembles a symmetry breaking scalar field potential. The first slow-roll parameter, that is related to the measure of accelerated expansion during inflation, is The accelerated expansion occurs while < 1, and in our model, we have 1 for a a 0 . Inflation ends when ≈ 1, and from the evolution of the slow-roll parameter in Figure 2 we can clearly see that the inflationary regime will end as the evolution approaches a 0 .
With the background equations, the EoS parameter (21) and the Klein-Gordon equation (33) we can show that (45) See Figure 2 for the behavior of the equation (45) as a function of the scale factor (a/a 0 ) in two different situations, with n = 1 and n = 2.
The second slow-roll parameter tells us how long the accelerated expansion will be sustained. It is related to the smallness of the second time derivative of the scalar field, and we can write where dN = Hdt and V ,φφ = d 2 V /dφ 2 . During inflation we have that |η| < 1. The slow-roll conditions are , |η| 1 .
See [29] for more discussions on the slow roll parameters. As already mentioned before, after the inflationary stage, for a a 0 , the linear term of the EoS (14) dominates, which means inflation will come to an end, and the Universe will be dominated by the linear component of the EoS (14). If α = 1/3, for example, the Universe undergoes a radiation era, which would correspond to the standard Big Bang model, in which the Universe undergoes a radiation dominated era. Although it is possible to reheat the Universe using the scalar-field formalism and the potential (42), we can not say the same about the fluid formalism. A complete reheating analysis should be made in this model in order to see if the Universe can be reheated.

Primordial quantum perturbations: scalar perturbations
In this section we calculate scalar perturbations generated in the early Universe. Scalar quantum fluctuations can be the source of the seeds that originated the large scale structures we see today. First we introduce the conformal time τ , such that dt = a(t)dτ .
During the inflation we have where a e is the scale factor at the end of inflation. As H is approximately constant during this period we can consider that The scale factor at the end of inflation is much large than in the middle, (a e a). So, we find that In order to find the fluctuations that originated the large scale structures we introduce the perturbed metric where Φ(τ, x) is the gauge-invariant Bardeen's potential [26]. We substitute the metric (52) in the Einstein's field equations (12), and keeping only the first order terms we find where δT µ ν is the gauge-invariant perturbed stress-energy tensor and H = a /a is the Hubble parameter in terms of the conformal time.
Using the hydrodynamics description of the polytropic fluid we first perturb the density ρ → ρ + δρ and equations (53) and (55) will be It is easy to show that δp = c 2 s δρ. We join equations (56) and (57) to find where k is the modulus of the wave-number and we made ∇ 2 Φ = −k 2 Φ. To obtain Φ during inflation we use the scale factor (51) and the transformations and we find With help of µ = |τ |g we will get the Bessel differential equation with z = k|τ |c s . The solution is where H (i) ν (z) are the Hankel functions of order ν = β 2 + β + 1/4. The same procedure we used to quantize the primordial gravitational waves [26,30,31] can be applied to the wave equation (61). This will transform the classical field Φ into the quantum field which allow us, using the two-points correlation function Φ (τ, x)Φ † (τ, y) , to find the quantum power spectrum . We can also write the quantum power spectrum as where A Φ is the power spectrum amplitude, k p is some fixed pivot scale and n Φ is spectral index, which can be calculated following the relation where aH = c s k is the horizon crossing condition, is given by In order to find the evolution of the density perturbation δρ we only need to perturb the conservation equations T µν ;ν [26] to find the hydrodynamics perturbations In the regime a a 0 the equation (69) will become where δ = δρ/ρ is the density contrast, and the classical solution is To find the quantum evolution of the density contrast δ we need to find a wave equation. This can be done if we substitute the conservation equation (69) into (70) and consider the appropriate approximations to have Making the transformation δ = a γ µ, where γ = − 1 2 5 + c 2 s = 1 2 α+1 n − 4 we find the wave equation and finally, with the transformations we will find which is the Bessel differential equation with the same solution (63), with the order ν = γ 2 − 9γ + 81/4. We can use again the already discussed quantization process to make δ →δ to find both power spectrum and spectral index where µ is described in terms of the Hankel function as µ = √ π 2 |τ |H (2) ν (z), and

Conclusions
In this work we assumed that the primordial Universe was filled with a fluid described by the equation of state p = αρ + kρ 1+1/n , (14), that is the sum of a standard linear equation of state and a polytropic term. The polytropic term, with n = 1, can be considered as a generalization of the standard Bose-Einstein condensate dark matter equation of state. Following [15,16], but letting the parameters α and n free, we show that the EoS (14) can describe a inflationary Universe in the case of attractive self-interaction. We found the slow-roll parameters and η.
In Figure 1 we plotted the scalar field's potential as a function of the scale factor ratio a/a 0 . We can see that the difference between the panel with the curves representing dust (α = 0) and radiation (α = 1/3) are the smallest. The behavior follows a pattern where the case with n = 1 the concavity of the curve is downward while for n = 2 and n = 3 the concavity is upward. On the other hand, the panel with figures representing the behavior of rigid material, either for n = 1, n = 2 and n = 3 have the same concavity within the limits of the scale factor ratio a/a 0 . Moreover, in this same figure we note that in a model of Universe filled only with dust and n = 1, which represents the BEC model, the slow-roll period ends earlier when compared with the stiff matter or radiation cases. As can be seen in Figure 1, the case where n = 1 the slow-roll period lasts longer than the cases with n > 1. In all cases the inflationary era is longer for stiff matter (α = 1).
We calculated the slow-roll conditions for the scalar field during inflation (see Figure 2). We can see that in the model with stiff matter the slow-roll period, compared with the scale factor ratio a/a 0 , is longer than the scenarios with dust (α = 0) and radiation (α = 1/3) for both n = 1 and n = 2. This situation indicates that the accelerated expansion of the Universe with stiff matter is slower than with dust and radiation. This characteristic may have important consequences in the process of evolution of the Universe, seen that the presence of stiff matter in FRW cosmological models produces an abundance of relic species of particles after the Big Bang due to the expansion and cooling of the Universe [32]. The presence of stiff matter in FRW cosmological models may also help explaining the baryon asymmetry and the density perturbations of the right amplitude for the large scale structure formation in our Universe [33], and may also play a important role in the spectrum of relic gravity waves created during inflation [34]. These two important consequences may be changed due to the behavior of the slow-roll parameter for the model of stiff matter. Figure 3 shows the unnormalized power spectrum of Φ for various values of the parameters α and n as a function of wavenumber k. We can clearly see that the curvature power spectrum is not scale invariant, and we see more power for small scales. It's spectral index (68) is proportional to the fluid's speed of sound, which, during inflation, can be approximated to Although we are dealing with an inflationary phase, which prevents agglomeration of matter, we also study the density contrast of the polytropic fluid, that has the quantum spectral index (which is also a function of the fluid's speed of sound) give by equation (79). The power spectrum behavior can be seen in Figure 4. In the left upper panel and in the right bottom panel of Figure 4 we see that the power spectrum is scale invariant for α = 1 and n = 1. For others combinations of α and n we see more power for big scales.
It is known that when we consider a fluid with negative pressure, the equivalence between hydrodynamical and field representation exists only at the background level: at perturbative level, the model behave in a complete different way [28]. Hence, in situations where negative pressures are concerned, a field representation leads to a much more complete scenario, being closer to a realistic model. We hope to present a more general analysis involving the comparison between the hydrodynamical model and the scalar representation of this polytropic equation of state in a future study, both in terms of background level and in the pertubative level in order to verify the behavior of the power spectrum in the early Universe.
To have a more robust analysis of the behavior of the cosmological models with the polytropic equation of state we can use here the Bayesian chi-square χ 2 minimization technique to limit the different parameters of the EoS for a viable cosmological model considering the observational data available today, but we have here a more complex situation: one primordial inflationary phase described by the polytropic equation of state with k < 0 and n > 0 and a current phase of accelerated expansion described by the same equation of state but with k > 0 and n < 0. In practice they are two different cosmological models. The idea was already developed in [9,15,16]. We intend to use a combination of these two models in order to apply the statistical techniques mentioned above. The best-fit values of the model parameters are then determined from the chi-square function to study the evolution of the Universe. We plan to show this comparison with observations in a future work.
To summarize, the polytropic equation of state represents a interesting scenario to study the evolution of the Universe. The similarities with the models that are described by a linear equation of state, more than a simple coincidence, should be investigated with others kind of representations that not just the hydrodynamical representation and with statistical methods to verify the feasibility of the model.