Dark matter as the Bose–Einstein condensation in loop quantum cosmology

We consider the FLRW universe in a loop quantum cosmological model filled with radiation, baryonic matter (with negligible pressure), dark energy, and dark matter. The dark matter sector is supposed to be of Bose–Einstein condensate type. The Bose–Einstein condensation process in a cosmological context by supposing it as an approximate first-order phase transition, has already been studied in the literature. Here, we study the evolution of the physical quantities related to the early universe description such as the energy density, temperature, and scale factor of the universe, before, during, and after the condensation process. We also consider in detail the evolution era of the universe in a mixed normal–condensate dark matter phase. The behavior and time evolution of the condensate dark matter fraction is also analyzed.


Introduction
When one considers a universe following the standard Einstein cosmology, that is, when the dynamics of the universe is described by the general relativity equations, one deduces that there exists a mysterious singularity at the beginning of time, the so-called big bang. This singularity may be considered as a deficiency of Einstein cosmology at high energies [1]. Actually, one can assert that the big bang implies the breakdown of general relativity at scales with high energies, whereas we know from the observational evidence, such as the existence of the cosmic microwave background, that the big bang model works at scales lower than Planck energy scale. At these scales, the universe was full of a hot photonbaryon combination almost described by a radiation fluid. This hot combination was cooling down as the universe was a e-mail: atazadeh@azaruniv.ac.ir b e-mail: f.darabi@azaruniv.edu c e-mail: mousavi@azaruniv.ac.ir experiencing the expansion, and due to the different scaling behaviors of pure radiation and pure matter, the energy density of non-relativistic matter started to dominate over the energy density of radiation and led to the formation of structures.
Loop quantum gravity (LQG) is one of the potential candidates for the study of quantum gravity [2]. Homogeneous and isotropic space time, reduces loop quantum gravity LQG to loop quantum cosmology (LQC) [3]. According to this theory, due to the effective quantum gravitational effects, the big bounce happens to remove the big bang singularity, after which a super-inflation phase occurs, and then the universe enters a normal inflation regime [4,5]. The loop quantum cosmology effects manifest themselves in a new effective form of the modified Friedmann equation. These effective equations have been derived in the literature for FLRW loop quantum cosmologies considering different matter sectors in the early universe [6][7][8][9][10][11][12].
An extended scenario of the matter bounce cosmology, considering a single scalar field with an approximately exponential potential, was proposed as an alternative to the slow roll inflation, in the context of LQC, in which the universe has experienced a quasi-matter contracting phase with a variable equation of state parameter [13]. This matter bounce scenario, in the teleparallel version of LQC, was shown to be in good agreement with the new BICEP2 data [14]. A Gauss-Bonnet extension of LQC, by introducing holonomy corrections in modified f (G) theories of gravity was also developed, where the authors have provided a perturbative expansion in the critical density as well as a parameter characteristic of LQG, and one obtained leading-order corrections to the classical f (G) theories of gravity. They also presented a reconstruction method which makes possible to find the LQC corrected f (G) theory capable of realizing various cosmological scenarios [15]. In another work, in order to avoid singularities of f (R) = R + a R 2 model, holonomy corrections to this model were introduced in Einstein frame and a detailed analytical and numerical study was performed when holonomy corrections are taken into account in both Jordan and Einstein frames. They obtained, in Jordan frame, a dynamics which is different qualitatively, from the one of the original model, at early times. According to this dynamics, the universe is not singular, neither at early times in the contracting phase, nor at bouncing to enter in the new expanding inflationary phase. This dynamics may lead to better predictions for the inflationary phase in comparison to the current observations [16].
Motivated by the above mentioned success of LQC in describing the very early stage of inflationary universe, we are encouraged to study the evolution of the universe after the inflationary era and investigate the possible establishment of a specific kind of dark matter (see below) and study its contribution in the subsequent evolution of the universe. Explicitly, we assume that after inflation the universe includes ordinary dark matter (bosonic particles), dark energy, radiation, and baryonic matter, and that the ordinary dark matter component in the sufficiently cooled universe has gone through a phase transition provided by the Bose-Einstein condensation, which might have happened during the early stages of cosmological evolution of the universe with a low temperature comparable to the critical temperature for Bose-Einstein condensation T cr ∼ 2πh 2 n 2/3 /mk B where m is the particle mass, n is the particle density and k B is Boltzmann's constant [17][18][19][20][21][22][23][24][25][26][27][28][29][30][31][32].
Particularly, in the Bose-Einstein Condensation (BEC) model, the dark matter can be described as a non-relativistic Newtonian gravitational condensate, whose pressure and energy density are related by a barotropic equation of state. Based on the above arguments, we assume the possibility that the Bose-Einstein condensation might have happened during the early stages of cosmological evolution of the universe with a temperature comparable to the critical temperature for Bose-Einstein condensation T cr ∼ 2πh 2 n 2/3 /mk B where m is the particle mass, n is the particle density, and k B is Boltzmann's constant [17][18][19][20][21][22][23][24][25][26][27][28][29][30][31][32].
The modification of cosmological evolution of the early universe, in the presence of condensate dark matter in standard cosmology, has already been studied in [33]. In this paper, we generalize this model from standard cosmology to LQC. The organization of the present paper is as follows. In Sect. 2, we outline the modification of Friedmann equations in LQC model and the basic properties of the normal and BEC dark matter in addition to the presentation of the relevant physical quantities and the equation of state. In Sect. 3, we give a brief explanation about the cosmological dynamics of the Bose-Einstein condensate. Next, we introduce the numerical values of the cosmological parameters at the condensation point, and then obtain the important h(t) parameter as a volume fraction of matter in the Bose-Einstein condensed phase which explains the evolution of the dark matter energy density during the transition process, and extract the order of BEC time interval in LQC model. In the following, we study the post condensation phase in LQC model describing the evolution of the scale factor parameter under the effect of the existence of the BEC dark matter in LQC. The paper ends with a conclusion.

Bose-Einstein condensate dark matter in the LQC universe
In this section, we intend to review the Bose-Einstein condensation and explore the dynamics of the condensation and some of the relevant cosmological implications. We consider the line element of the flat Friedmann-Robertson-Walker metric where a is the scale factor presenting the cosmological expansion. We describe the matter by the perfect fluid energymomentum tensor, where (ρ, p) stands for (ρ b , p b = 0), (ρ rad , p rad ), and (ρ χ , p χ ) for baryonic matter, radiation and dark matter, respectively. Note that we do not consider the interaction between these energy terms, in other words the energy of these components are individually satisfying the conservation equation.
Taking into account the well-known loop quantum gravity constraints and coupling the matter to classical phase space of FLRW universe, we obtain [34][35][36][37][38][39] where H matt is the matter Hamiltonian, γ ≈ 0.24 [40,41] represents the Barbero-Immirzi parameter, and λ ≈ 2.27 is the length gap. Meanwhile, b is the conjugate momentum and V is proportional to the physical volume of a cubical cell, with unit comoving volume V = a 3 2πγ . A dot means the derivative with respect to the time. By rewriting Eqs. (3) and (4), we have where H is the Hubble parameter and the matter density is limited as ρ < ρ c where ρ c = As we mentioned, the effective equations have been derived for FLRW loop quantum cosmologies considering different matter sectors in the early universe [6][7][8][9][10][11][12]. So, we may use the modified FLRW equation in the loop quantum cosmology model, namely Eq. (7), and assume the matter density ρ as the combination of different matter components ρ b , ρ rad , ρ χ , and indicating the baryonic matter, radiation, dark matter, and cosmological constant, respectively, so that The second modified FLRW equation is also given by and also for the energy density conservation equation we havė Now, regarding the cosmological evolution of radiation and baryonic matter we consider the relations ρ rad = ρ rad,0 /(a/a 0 ) 4 and ρ b = ρ b,0 /(a/a 0 ) 3 , respectively; meanwhile ρ b,0 and ρ rad,0 are the energy densities corresponding to their values at a = a 0 , and also for the dark matter we work with a general form of density as ρ χ = ρ χ,0 / f (a/a 0 ). It should be noted that f (a/a 0 ) is an arbitrary function which depends on the special dark matter model. As we know we have the critical density ρ cr,0 = 3H 2 0 , where H 0 is the value of the Hubble parameter at a = a 0 . In addition, we can write the dimensionless parameters known as density parameters as i,0 = ρ i,0 /ρ cr,0 , with i = b, rad, χ. Using these relations helps us to write the following new form of the modified FLRW equation: where is the dark energy density parameter. Thus, we can write the constraint b,0 + rad,0 + χ,0 + = 1. Now, we suppose that in the early universe the bosonic particles, with mass m χ and a very high temperature T , were in equilibrium state with other particles in a hot relativistic plasma. After the expansion of the universe and rapid decrease of the temperature it was decoupled from this equilibrium state at a high chemical potential μ m χ or low decoupling temperature T D T . At this decoupling temperature, we assume that the bosonic particles are still in kinetic equilibrium among themselves in an isotropic gas with low temperature. Consequently, for the spatial number density we have where g is the helicity state number, h is Planck's constant, and f ( p) is defined as follows: where E is the energy E = p 2 + m 2 χ c 4 and p is the momentum of the particle. Due to the redshift in the momentum of particles, we have p(t)a(t) = p D a D where p D and a D denote for the momentum of particles and the scale factor of the universe, respectively, at the temperature T D denoting for the decoupling temperature from the rest of plasma. Besides, we have the scaling evolution relation for the number density given by n χ ∼ a −3 [47,48]. For the extreme-relativistic case E ≈ pc (μ = μ D a D /a and T = T D a D /a), the distribution function is derived as [47,48]. For a frozen distribution of dark matter, the associated energy density and kinetic energymomentum tensor T μ ν are as follows, respectively: The mentioned pressure is described by where we have used v = pc 2 /E [49]. The density of dark matter ρ χ for the non-relativistic case with E = m χ c 2 and P ≈ m χ v χ is given by ρ χ = m χ n χ whereas the dark matter pressure is given by which leads to where we have used σ 2 =< v x 2 > /3c 2 , σ being the onedimensional velocity dispersion. By means of the dark matter conservation equatioṅ we can find the following general solution: where ρ χ,0 is the density of dark matter at the present value of scale factor a = a 0 . In the standard model, for the description of the dynamics of our universe with a normal dark matter we work with the following equation: However, in the LQC model we will have the following modification: Since we aim to extract numerical results for dynamical BEC quantities, we have to mention the values of several parameters and quantities such as the Hubble constant H 0 = 70 km/s/Mpc = 2.27 × 10 −18 s −1 , the Hubble time t H = H −1 0 = 4.39 × 10 17 s, the critical density ρ cr,0 = 9.24 × 10 −30 gr/cm 3 , and the present day dark matter densities b,0 = 0.045, χ,0 ≈ 0.228, rad,0 = 8.24 × 10 −5 , = 0.73, respectively [50]. It is worth to recall that the following dark matter case is non-relativistic, so the global cosmological evolution of the universe will not be affected by the variation of the numerical values of σ 2 .
As is explained in the introduction, at very low temperature all dilute Bose gas particles are condensed to the same quantum ground state which leads to the BEC formation. When the particle's wavelengths overlap, they will experience correlation with each other and the mean inter-particles distance l becomes definitely smaller than the thermal wavelength λ T . This is possible when the temperature is almost T cr ≈ 2πh 2 ρ 2/3 /m 5/3 k B , where m is the particle mass in the condensate, k B is Boltzmann's constant [51] and ρ is the density. Under the condition that the universe is sufficiently large and the temperature is low enough, there will be a coherent state in progress due to the overlap of the particles wavelengths. Our main assumption is that the dark matter halos may be constructed by a strongly coupled dilute BEC at absolute zero. As a result, all the particles are in the condensate. Therefore in this situation, just binary collisions at low energy are considerable. These collisions are specified by just one parameter l a , s-wave scattering length, which is thoroughly independent of the two-body potential. Consequently, the interaction potential can be replaced by an effective interaction, [51]. The Gross-Pitaevskii (GP) equation, describing the ground state properties of the dark matter, for the dark matter halos can be obtained from the GP energy functional as follows: where U 0 = 4πh 2 l a /m χ [51] and obviously we have applied ( r ) as the wave function of the condensate. In Eq. (22), the first and second terms are the quantum pressure and the interaction energy, respectively, and the third term is the gravitational potential energy. It should be noted that the condensate dark matter mass density is and that the wave function satisfies the normalization condition N = | ( r )| 2 d r , where N is the total number of particles. Eventually, by imposing the variational procedure while μ plays the role of the chemical potential. In addition the gravitational potential V follows the well-known Poisson equation, ∇ 2 V = 4π Gρ. Having considered the time dependency, the generalized GP equation, showing a gravitationally trapped rotating BEC, is found as By means of the well-known Madelung representation of the wave function, we can reach our purpose more easily than using the last generalized GP equation. In order to go through this, we should consider the following representation of : where S( r , t) represents the classical action. Inserting the last expression of ( r , t) into Eq. (25) leads to the following differential equations: where v = ∇S/m χ is the velocity of the quantum fluid and V Q is the quantum potential, given by By defining u 0 = 2πh 2 l a m 3 χ , the condensate effective pressure is obtained The quantum pressure term makes a great contribution just close to the boundary of the condensate provided that the number of particles in the gravitationally bounded BEC becomes large enough. As a result, the quantum stress term in the condensate equation of motion can be ignored, which is the so-called Thomas-Fermi approximation. When the number of particles in the condensate goes to infinity, the mentioned Thomas-Fermi approximation becomes accurate. It is also related to the classical limit of the theory. According to its definition, the velocity field is irrotational, following the condition ∇ × v = 0. In accordance with Eqs. (10) and (30) we can writė which has the following general solution with the integration constant C χ : Having imposed the condition ρ χ = ρ χ,0 for a = a 0 , we can find the density of the condensate as where

Loop quantum cosmological dynamics of the Bose-Einstein condensation
For an interacting Bose system case, where all dark matter particles are supposed to be in a single-particle state, the transition procedure from the normal to the condensed phase, and the order of the phase transition has been extensively studied in the recent literature [33,52,53]. In principle, BEC represents a spontaneous breaking of U (1) gauge symmetry, with the condensate fraction showing the order parameter, leading to the second-order phase transition [54,55]. Nevertheless, it has been shown that BEC is a first-order phase transition [56,57]. From another point of view, in thermodynamics the first-order phase transitions result in a genuine mathematical singularity. The question "can finite systems in nature explicitly show such a behavior?" is a mysterious controversial question in physics [58]. Nonetheless, a significant study about the thermodynamic instability in an ideal Bose gas, confined in a cubic box, has shown that a system including a finite number of particles can exhibit an interrupted phase transition, describing a real mathematical singularity, with constant pressure [57]. Of course the mentioned result was obtained without considering a continuous approximation or a thermodynamic limit. Therefore, according to above discussion, we will assume that the intrinsic dynamics of BEC is described in the framework of a first-order phase transition.

Cosmological parameters at the condensation point
As we know, in general, the chemical potential μ in a physical system is a function which depends on the temperature and particle density n = N /V (n is the total particle number and V is the volume). By considering the property of extensivity of the Helmholtz free energy F = F(N , V, T ), we are allowed to write not only F = V f (n, T ) but also F = Nf (v, T ), where v = V /N = n −1 and f = nf . Using the Helmholtz free energy we can obtain μ(n, T ) = (∂ f /∂n) T and p(n, T ) = −(∂f /∂v) T , having the same information. Imposing the thermodynamical requirements, the chemical potential and the pressure are supposed to be singlevalued which means that for arbitrary values of n, v, and T , there must be just one choice for μ or p, respectively [56]. Hence, the pressure at the transition moment during the BEC process must be continuous. Imposing this thermodynamic condition on Eqs. (18) and (30) helps us to fix the transition density ρ cr χ as follows: Apparently, in order to have the numerical value of the transition density, we should know the scattering length, the dark matter particle mass and the particle velocity disper-sion. Supposing typically the mass of dark matter particle, the mean velocity square, and the scattering length to be of the order of 1 eV, 10 15 cm 2 /s 2 , and 10 −10 cm, respectively, we can write Therefore, the critical temperature at the BEC moment is obtained by where ζ(3/2) is the Riemann zeta function, or equivalently The critical pressure of the dark matter fluid at the critical point is given by Consequently, we can find the critical redshift as follows: The obtained critical values of density, temperature, pressure, scale factor, and redshift are important in the study of the subsequent LQC evolution of the universe during the Bose-Einstein condensation toward the end time of the phase transition when the universe enters in the complete Bose-Einstein condensed dark matter phase. We will use these critical values in the study of this evolution in the following subsection.

Loop quantum cosmological evolution during the Bose-Einstein condensation phase
In the first-order phase transition, we have constant temperature and pressure, respectively as T = T cr and p = p cr . Moreover, the enthalpy W = (ρ + p)a 3 and entropy S = sa 3 are also conserved quantities. The phase transition begins through a decrease in the density of dark matter and increase in the density of Bose-Einstein condensed state. This is because the density ρ cr χ (T cr ) ≡ ρ nor χ is converted to ρ χ (T cr ) ≡ ρ BEC χ during the phase transition. One may define a useful time dependent parameter h(t), playing the role of ρ χ (t), which is called the volume fraction of the matter in the Bose-Einstein condensed phase By simplifying (44), we can extract the time evolving dark matter energy density during the transition process as where we have defined By looking at Eq. (44), we can determine the beginning and final points as h(t cr ) = 0 (t cr being the beginning time of transition in which ρ χ (t cr ) ≡ ρ nor χ ) and h(t BEC ) = 1 (t BEC being the end time of transition in which ρ χ (t BEC ) ≡ ρ BEC χ ), respectively. At the end time of phase transition the universe enters in the Bose-Einstein condensed dark matter phase. As mentioned, by using h(t) instead of ρ χ (t) in Eq. (10), we can find the following relation: where The last expression implies that in general we have r < 0, r ∈ (−1, 0), and also n χ < 0, because ρ BEC χ < ρ nor χ . Clearly, imposing the condition h(t cr ) = 0, we can obtain the relation between the scale factor and h(t) from Eq. (47) as where a cr = a(t cr ). The value of the scale factor at the end of phase transition is as follows: This is an important result, since this condensation process has modified the expansion rate of the universe. Subsequently, we can find the evolution equations for the baryonic matter and the radiation density during the phase transition, respectively, as and Using all these equations, we consider the time evolution of the volume fraction of the condensed matter dark energy density during the BEC process, which describes the dynamics of this process in the cosmic history in the LQC model, Here, we have used τ as a dimensionless time variable τ = H 0 t and χ,nor = ρ nor χ /ρ cr,0 . Considering Eq. (48) and p cr /ρ nor χ = σ 2 1, we deduce the approximation r ≈ n χ in the following calculations. Another approximation comes from neglecting the radiation energy density contribution, so Eq. (53) can be integrated to extract h(t) where we have defined and In order to obtain the necessary time interval for the entire converting process of the total normal dark matter into the Bose-Einstein condensed phase, we put h(t = t BEC ) = 1 in our last calculation, By inserting the value of dark energy density parameter as 0.68 and also considering the standard values ρ c 0.82M 4 p 1.84 × 10 170 eV/cm 3 , a = 10 −10 cm, m = 1 eV, and tr = 5.02 × 10 8 , we can obtain the following expression for the conversion time of the dark matter to the Bose-Einstein condensed phase: For instance, for r = −0.2 and t H = 4.39 × 10 17 s, we will have t cond 3.82 × 10 16 s. Since it is thoroughly useful to depict the time evolution of h(t), we have represented (54) for different values of r in Fig. 1.

The post-condensation phase in LQC model
After the complete conversion to the BEC phase, the postcondensation phase begins, which starts at t = t BEC and a BEC = a cr (1+r ) −1/3 . Similar to the era, during the condensation phase, in which the cosmological dynamics is affected by the BEC (see (50)), the Bose-Einstein condensed dark matter, in the post-condensation phase, may change the cosmological dynamics of the universe in LQC model. Moreover, the Bose-Einstein condensed dark matter may have considerable impact on the subsequent structure formation in the universe. Therefore, the post-condensation phase and the corresponding cosmological values in LQC model are important and deserve more scrutiny.
According to (33), we can obtain the cosmological density of the dark matter as follows: where the numerical value of ρ 0χ is obtained from (34), ρ 0χ = 1.63 × 10 −15 × (l a /10 −10 cm)(m/10 −33 gr) −3 1 + 1.63 × 10 −15 × (l a /10 −10 cm)(m/10 −33 gr) −3 . (60) Clearly, the density (59) should be positive, which leads to the condition (a cr /a 0 ) > ρ0χ 1/3 (1 + r ) 1/3 . Simplifying Eq. (21) and defining BE = χ,0 1+ρ χ,0 u 0 /c 2 gives us the time evolution of the scale factor in the BEC phase as follows: By integration of this relation, we obtain exactly an expression of the following form: where C is an integration constant, determined from the constraint a = a BEC at t = t BEC as For the universe presented here, filled with radiation, dark energy, baryonic matter, and Bose-Einstein condensed dark matter, the equation describing time evolution of the scale factor is given by which can be integrated by considering the initial conditions a(t BEC ) = a(0) = a BEC . To study the time evolution of the scale factor in such a universe which contains BEC dark matter, it will be useful to depict it for different values of the ρ 0χ in Fig. 2. It is remarkable that the existence of the condensed dark matter changes the cosmological dynamics of the universe in the post-condensation phase, and the magnitude of this change decreases with the increase of the Bose-Einstein condensate parameter ρ 0χ . This behavior is expected because as ρ 0χ increases, the corresponding gravity contributes stronger to the dynamics and decreases the expansion rate. On the other hand, it is seen from Fig. 2 that as ρ 0χ increases, the Bose-Einstein condensation phase (where a = a BEC ) occurs at earlier times. This is also justified because as ρ 0χ increases, the ρ BEC χ increases as well and so the scale factor a BEC corresponding to this increased density tends toward the scale factors having higher density, which naturally occurs at earlier times.

Conclusions
When the Bose-Einstein condensation happens, dark matter contains two phases, the normal and the condensed phase, respectively. Under a thermodynamic condition, we expect to have the continuity of the pressure profile of these two phases right at the beginning of the condensation process. This leads to uniquely determination of dark matter density at the condensation time as well as other thermodynamical quantities, like temperature and pressure. The explicit numerical values of the condensation quantities is related to mean square velocity of the normal dark matter σ 2 and the mass m χ of the dark matter particle and the scattering length l a . As a result of the uncertainty of the values of σ 2 , m χ , and l a , apparently it is difficult to predict the precise and exact cosmological moment of the Bose-Einstein condensation and consequently the corresponding cosmological quantities. Nevertheless, by considering some "standard" numerical values, it will be possible to obtain a total qualitative picture of the transition. Consequently, we estimate the values of the dark matter mass and the mean velocity of the non-relativistic dark matter particles, of the order of 1 eV [19,59] and 900 km/s, respectively. According to the general analysis of the condensation process, the first phase of condensation belongs to the standard cold dark matter ( CDM) model. In this paper, we have generalized this analysis and have gone through the loop quantum cosmology context. With the standard numerical values of the dark matter quantities, the condensation was begun at a redshift around z = 1200. Here, we have studied the time evolution of the fraction factor h(t) of the condensed dark matter through the loop quantum cosmology. In the next step, we have extracted the modified phase transition time interval in LQC, which is approximately of the order of Hubble time and shows the viability of our results. We have studied the evolution during the Bose-Einstein condensation phase and found that this condensation process modifies the expansion rate of the universe. Also, we have studied the post-condensation phase in LQC model, describing the time evolution of the universe right after the end of phase transition, which leads to (64). The cosmic time dependence of the global formation of BEC dark matter (at t BEC ) is an advantage, compared to the ordinary dark matter, which makes BEC dark matter to be preferred than the ordinary dark matter. This is an important feature which makes BEC dark matter to be similar to the global dark matter candidates like neutrinos, WIMPs, axions and LSP. The ordinary dark matter has not such capability and cannot form global sources of dark matter at a definite cosmic time, and that is why we have favored BEC dark matter in the present study. We have plotted the time behavior of Eqs. (53) and (64) in Figs. 1 and 2 for different values of r and ρ 0χ , showing the definite effect of the phase transition and Bose-Einstein condensed dark matter density in the modified cosmological history in the loop quantum cosmology. A major point regarding the present study is that assuming two phases of dark matter, one as the bosonic particles before BEC, and the other one as BEC dark matter, one may ask: "is there any experimental indication for the occurrence of BEC, the existence of these two phases, and the preference of one to the other one?" and "why is it necessary to distinguish ordinary dark matter from BEC dark matter?". Of course, at present there is no direct evidence for the occurrence of BEC, the existence of these two phases, and the preference of one to the other one. However, in order to answer properly these questions, we may try to answer the alternative question: "is there any indication for distinguishing between the dynamics of the universe having bosonic particles without the occurrence of BEC, and the one with BEC dark matter". In this regard, we may point out that the time evolution of the scale factor with BEC is different from that without BEC for many possible reasons, such as: (1) the equation of state parameter of BEC dark matter is different from that of bosonic particles without BEC, due to the first-order phase transition, (2) in the absence of BEC, the bosonic particles may interact or participate in some possible transformations into other matter components, like visible matter or dark energy, each resulting in different dynamical history of the universe; however, with BEC the condensed bosonic particles cannot decay into other forms of matter and they remain as the BEC dark matter permanently, with a definite dynamical history distinguishable from the variety of dynamics corresponding to the regular bosonic particles having different kinds of interaction with other matter components.