Evolution of spherical overdensity in Chaplygin gas model

Even though many scalar field models of dark energy have been considered in the literature, there is another interesting class of dark energy models involving a fluid known as a Chaplygin gas. In addition to describing the dark energy, both scalar-tensor model and the Chaplygin gas model are suitable candidates for explaining the spherical cosmological collapse. One of the most well-known scalar field models is the quintessence model, which was first introduced to explain an accelerating expanding universe. Using a special form of the quintessence model that is equivalent to Chaplygin gas, we describe evolution of a spherical collapse. We study the cosmological properties of the quintessence field with a special potential. In addition to the quintessence model, that can be converted into a Chaplygin gas model in a particular case, we claim that the fixed-potential tachyonic model is equivalent to the Chaplygin gas model. In this work, we obtain the spherical collapse parameters: the virialized over density parameters, radius, the energy density at the turnaround moment, etc. We compare the results of the proposed model with the standard model of cosmology and the Einstein–de Sitter model. We show that the formation of the large-scale structures within the framework of a Chaplygin gas model happens earlier than predicted in the standard model.


Introduction
Almost decades after the introduction of the Einstein-de sitter dust model, observations of the redshift in type Ia supernova a e-mail: am.rezaei@tabrizu.ac.ir (corresponding author) b e-mail: mohammad.ahmadii7373@gmail.com c e-mail: yousef_izadi@uml.edu d e-mail: smostafanazh2016@fau.edu e e-mail: wam@fau.edu suggested a positive acceleration of the universe [1]. These results led to the proposal of an unknown exotic matter which is called dark energy and it satisfies the condition ρ +3P < 0 where ρ represents the density of the matter and P denotes the pressure. There have been many candidates for the dark energy, the simplest model is based on the cosmological constant which is not dynamic. The cosmological constant faces two issues, "fine tuning" and "the coincidence problem" [2]. Models of interacting dark energy were introduced to solve the problem of coincidence. Quintessence [3] and K-essence [4] models were proposed to solve such issues, these models are based on either non-minimally or minimally coupled scalar fields [5,6]. Although these models face the challenge of fine-tuning, which states that, in terms of mathematical laws and probabilities, the universe could have been created in innumerable different ways. However, our universe is only a very special case. In recent literature some models have been investigated in which, in order to solve the fine tuning problem, it is necessary to consider dark energy in the initial stage of the universe [1]. Therefore, to explain the number of structures at present, the contribution of dark energy cannot be ignored. On the one hand, the process of structure formation should occur earlier and much slower than predicted in the standard cosmological model. Additionally, the analytical calculations show that by considering some assumptions for the collapse of a homogeneous and isotropic sphere, one can conclude that the number of structures formed within the scalar-tensor models is more than predicted in the ΛC DM model. Also, numerical analyses show that the classical spherical model is in good agreement with the simulations [7]. The dynamics of the overdensity parameter depend on the evolution of the flux, which causes the universe to expand. So, dark energy models, such as scalar-tensor models, can describe the growth of structures.
One of the most interesting ideas for the unification of dark matter and dark energy is the Chaplygin gas model in which the dark sector of the universe is expressed in terms of a single component that serves for both dark energy and dark matter. Chaplygin gas is a perfect fluid that behaves like a pressureless fluid in the early universe while it causes the accelerated expansion of the universe at the late time. It has been shown that the resulting evolution of the universe is compatible with the observational results [8]. The model states that the value of the effective cosmological constant is increasing. The simplest model of the Chaplygin gas satisfies the following equation of state in which A is a constant. In [9], the scenario based on the dynamics of a generalized D-brane in a space-time (D +1, 1) is investigated. They stated the conditions for homogeneity and demonstrated that the equation of state describes the evolution of a universe from a matter-dominated phase to the cosmological constant dominated phase. They wrote the effective equation of state as follows which is the state equation of the generalized Chaplygin gas with 0 < α ≤ 1. In [10], the constraints on this model of cosmic background radiation are discussed. In this paper, from the results of Archeops for the location of the first peak and, from the results of BOOMERANG for the location of the third peak, also, from the supernova, high-redshift observation, and gravitational lensing statistics [11], it is shown that this model is completely distinguishable from the standard model of cosmology. In line with previous works in [12], the WMAP temperature power spectrum and supernova data are compared with the generalized Chaplygin gas model, and it has been shown that the parameter α is in the range (0, 0.2) with 95% confidence. The authors have also shown that the generalized model at state α = 1 was ruled out as a candidate for dark energy with an accuracy of more than 99.99%. Furthermore, by examining this model in the non-flat background, one can see that the non-flat state corresponds to a flat case with an accuracy of 68% [13]. In [14,15] the modified Chaplygin gas model with the equation of state P = Aρ − B/ρ α in which A and B are constant, is considered. Assuming that the equation of state for the modified model of the Chaplygin gas is valid from the radiation period (A = 1/3, and for very large densities) to the current time (small densities), they show that their model can describe the accelerated expansion of the universe. In [16], it is shown that the inhomogeneous model of Chaplygin gas can explain dark matter and dark energy together in a geometric setting reminiscent of M-theory. Since the scalar field can explain both the holographic dark energy and the Chaplygin gas model, the relationship between the (interacting) holographic dark energy density and the energy density of the Chaplygin gas model is studied in [17,18]. Models that describe the dark energy issue not only explain the accelerated expansion of the universe, but they also play a fundamental role in the formation of cosmological structures. The large scale structures we see today are formed by the small perturbations in the inflationary phase of the early universe, and under the influence of gravity, these perturbations grow. In 1972 Gunn and Gott proposed the spherical collapse model [19], which led them to reach a simple explanation for the growth of the over dense structures. To understand the formation of the large scale structures and the growth of perturbations in a matter-dominated universe within the spherical collapse model, one can consider a spherical region with a radius that expands in time and with a non-uniform density. According to Birkhoff's theorem [20], the evolution of this radius depends only on the limited mass within the region. With the expansion of the cosmic background, this perturbed region begins to expand, and its radius reaches a maximum value. However, after a while, depending on its mass and due to its high density, it decouples from the Hubble background and starts evolving in a reverse process that decreases the radius or collapse in the non-linear region. Note that the collapse of this region does not proceed to a singular point. According to the virial theorem, the collapse of these perturbations ceases at half of the maximum radius, which leads to the formation of the structure. The spherical symmetry makes the problem easier to solve. Therefore, the spherical collapse model is a suitable model for understanding the evolution of the perturbations. Additionally, this model is successful in reproducing the simulation results [7,[21][22][23].
The studies and progresses mentioned in scalar-tensor theories were within the framework of Einstein's standard theory of general relativity. In 2004, Motta and Brock investigated the minimally coupled quintessence model under different potentials [24] and they concluded, in models with the assumption of non-minimal coupling of the scalar field, the evolution of the spherical overdensity parameter is quite different from the minimally coupled model. The non-linear evolution of structures within the framework of non-minimally coupled quintessence models was studied by Pace et al [7]. They conclude that the value of the virial overdensity and the linear density threshold for the spherical collapse are very close to the predictions of the standard cosmological model. Note that these two parameters play an essential role in the spherical collapse model. In the paper by Fan et al. [25], the non-minimally coupled model in the metric and Palatini formalisms is studied. It is shown that metric and Palatini formalisms lead to different results for linear growth rates. However, this difference is so small, and it is not observable now. Therefore, in this paper we study the non-linear evolution of structures in the framework of the spherical collapse model using the Chaplygin gas model. Section 2 presents an introduction to the quintessence model, and shows that by choosing an appropriate potential for the quintessence, one obtains the equation of state of the Chaplygin gas. In sect 3, this article studies the non-linear evolution of structures and shows that the quintessence model, in a special form, is equivalent to the Chaplygin gas model and plays an effective role in the spherical collapse model. We also discusse that the results of this model are the possibility for the formation of the structures in this framework occurs sooner than the cosmological standard model. Section 4 discusses that the tachyon field under a constant potential is equivalent to the Chaplygin gas model and examines the parameters of the spherical collapse model within this framework. Also, we show that the formation of structures under the Chaplygin gas model happens faster than what the standard cosmological model predicts. Section 5 contains a conclusion to this article.

Quintessence filed as a Chaplygin gas
In 1988, Ratra and Peebles [26] proposed a scenario in which one can use a scalar field which is universal, rolling, self-interacting, and homogeneous to describe the dynamical energy density of the universe. Models based on this scenario are known as quintessence. Also, the quintessence fields are one of the best candidates for the scalar field to explain the formation of the large scale structures in the universe. Following Ratra and Peebles, quintessence fields were studied more by Robert R. Caldwel, Duhal Dave and Paul Steinhardt [27]. Some people even considered it as the fifth fundamental force of nature. The action for the quintessence field can be written as where V (φ) is the scalar field potential. We assume that the field is only a function of time φ = φ(t). We use the FLRW metric to describe the expanding universe, so we have √ −g = a(t) 3 where a(t) is the scale factor. Then the pointlike Lagrangian of this field becomes One can easily write the equation of motion as If we choose the potential in this model as follows one can conclude that the scalar quintessence model is equivalent to the Chaplygin gas model (Appendix A). For the aforementioned potential, A and κ are constants. Many observations place constraints on the model parameters, κ and A should approximately take the values k 0.06, 0.28< A <0.550. The energy density and the pressure are obtained as One can use the expressions in (7) and (8) to find the equation of state parameter (EoS) as In the slow-roll approximationφ 2 << V (φ), we have P −ρ. In Fig.1, we draw a graph of the EoS parameter for the quintessence verse the redshift, and we compare it with the results of the ΛC DM and non-minimally coupled tachyon models (NMCT). As we see from the figure, ω is a constant (ω = −1) for the ΛC DM model. However, for the quintessence and the NMCT models, at high redshifts the value of ω is slightly different from the standard model, while at low redshifts, it takes a value which completely agrees with the standard model. One can use the components of the stress tensor to find energy densities of the field and matter as which leads to the continuity equationṡ By taking the time derivative of Eq. (10), we find the variations in the field density aṡ where Using the following change of variable one can write the change of field density with respect to the redshift Note that prime denotes the derivative with respect to z. The value of Ω φ approaches to zero for the large redshifts (early times), while from the observations we know that the value of this parameter for small z (present time) approximately equals to 0.7. From Fig. 2, one can see that the value of Ω φ in the quintessence model is very close to the prediction of the standard model of cosmology. Now, using the above Fig. 3 The variations of the dimensionless Hubble parameter with respect to the redshift for the quintessence, minimally-coupled tachyon and the standard cosmological models equations, we find the relation between the Hubble parameter and the energy density of the scalar field as follows where in which H 0 is the Hubble parameter at the present time (In Fig. 3, the dimensionless Hubble parameter E(z) is plotted for different models in terms of the redshift.). One can see that for z = 0, we have E(0) = 1. Formation of large scale cosmological structures such as super clusters is one of the main questions in cosmology and astrophysics.
Probably some factors such as inhomogeneous distribution of matter are responsible for the growth of the perturbations in the early universe that caused the formation of the cosmological structures during the period of matterdominant. To understand the perturbations growth in the matter-dominant epoch which formed the structures that we have today, we consider a slice of the expanding universe with a radius of R and non-uniform density ρ =ρ(1 + δ) ( where δ denotes overdensity parameter in a matter-dominant universe). Thus, the linear evolution of overdensity parameter (Appendix B) can be written as Here, we assume that the collapse is spherical and homogeneous. We also ignore the effects of rotation and shear. Evidently, a better approximation for the spherical collapse can be obtained if one considers an angular momentum. However, note that we do not see shear effect for a sphere, and the shear tensor equals to zero. Nevertheless, we write an equation of evolution of overdensity parameter in a non-linear In the equations mentioned above, δ, E, and a denote density contrast, dimensionless Hubble parameter, and the scale factor; respectively. From WMAP and Planck's observations we approximately have Ω m = 0.3. We numerically solve the non-linear growth equations of the perturbations. We imposed the initial conditions in a way to have the growth of non-linear perturbations larger than 10 7 . Note that these initial conditions are necessary for the formation of the large scale structures. In Figs. 4 and 5, we plot linear and nonlinear density contrast versus the scale factor and compare the results of the quintessence model with ΛC DM and coupled tachyon model. Therefore, the results of the quintessence model are similar to the standard cosmological model.

Spherical collapse via quintessence field
In this section, our goal is to understand how the perturbations in the universe grow within the framework of the quintessence model and form a dense structure. The intended area expands within an expanding background up to a maximum radius R max . Then the expansion of the perturbed region stops, and the process of collapse continues. However, it is important to note that the intended area does not collapse to a singular point, and the expanding process ends at half of the maximum radius, which is called virialized radius. Hence, it is essential to realize that the intended area does not collapse to a singular point, and in half of the maximum radius, it is virialized, making the collapse to stop [28,29].
To find the growth of the perturbations, we numerically solve the Eq. (22). As mentioned before, solving this equation is very sensitive to the choice of the initial conditions. Thus, to solve the equation and find the correct value for the linear overdensity parameter δ c . The value of this parameter in the Eds model is constant and equals to 1.68 [7]. However, one can see from Fig. 6 that overdensity parameter changes with the redshift, and eventually, both models share the value of the Eds model. To determine the virialized overdensity parameter, let us consider the sphere with density in the background of ρ b as defined below in which R is the radius of the structure, a c denotes the virialized scale factor and ρ is the density of the structure. Equation (24) can be rewritten as where x c is defined as a c /a ta . a c and a ta denote the scale factor in the virialization and the turn-around epochs; respectively. In the Eq. (25), λ and ξ present the radius of structures and perturbations in the turn-around time. In the Eds model, The value of Δ V is constant and equals to 178. However, in the other two models, it varies with the redshift. But for the large redshifts, Δ V has the same value as the Eds model. In Fig. 6, we plot the virialized overdensity parameter versus z c for the quintessence, ΛC DM, NMCT, and Eds models. Virialized overdensity is a constant for the Eds model, and it equals to 178. However, for the other models, it varies with the redshift. For the quintessence model and the standard model, the value of this parameter is approximately the same. But in the tachyonic model, the slope of the graph approaches the Eds model slightly faster, which means that the perturbations in the tachyonic model turn to virialized overdensity faster. Thus, the possibility for the formation of the structures in this model occurs sooner than the other two models. Fig. 6 The evolution of virialized overdensity parameter Δ V versus z c for the ΛC DM, quintessence, NMCT and Eds models Fig. 7 The behavior of λ versus z c for the ΛC DM, quintessence, Eds and NMCT models As we mentioned earlier, the growth of the perturbations during the collapse phase is virilized at half of the maximum radius, and so the structures can be formed. At the turnaround moment, we plot the radius of the cosmological structures λ for the four models, as shown in Fig. 7. One can see that λ takes its maximum value for the Eds model, which is 0.5. However, in the other models, the value of λ depends on the redshift and for higher redshifts, it tends to the Eds values.
The ξ parameter that appears in Eq. (25) is as follows It can be seen from Fig. 8 that the parameter ξ for the Eds model is a constant and equal to 5.6. But, for the other models, it varies with respect to virialized redshift, and for large redshifts takes a value which is close to what it has in the Eds model. But the coupled tachyonic model tends to the Eds model with a steeper slope, which means that within the framework of the tachyonic model, during virialization, the structures take the return path faster than other models. Therefore, it can be concluded that in this model the structures are formed faster than the other models.

Tachyon field with a constant potential
Let us assume that the Chaplygin gas acts as a tachyon field, so its Lagrangian becomes The potential in this model is a constant and equals to V 0 . Using the Euler-Lagrange method, one can obtain the equation of motion as Using the energy-momentum tensor, we find the energy density and pressure for the tachyon field The EoS parameter corresponding to the tachyon field is So, the energy density of the tachyon and matter field are as follows with the continuity equationṡ anḋ We know that the tachyon energy density iṡ Using the Friedman and the continuity equations, we find the expression foṙ Using the change of variable d d(ln a) = −(1 + z) d dz , the Eq. (35) can be rewritten as where prime represents the derivative with respect to z.
According to the Planck's recent results, this parameter is approximately 0.7, but tends to zero for the very large redshifts [7]. As discussed in the previous section, the spherical collapse model for a fixed-potential tachyonic model (FPTM) which is equivalent to a Chaplygin gas can be studied. In Fig. 9, the EoS parameter ω is plotted for both the FPTM and the ΛC DM models. One can see that, for the small redshifts, the value of ω φ approaches to −1. In Fig. 10, the behaviour of energy density of the FPTM is compared with the standard model of cosmology. The plot shows that the variation of Ω(z) for the different redshifts is similar for both of the models. The Hubble dimensionless parameter E(z) is plotted in Fig. 11. This graph shows the compatibility of the fixedpotential tachyon model with the standard cosmological model. We solve Eqs. (22) and (23) for the FPTM and obtained the growth of linear and non-linear tachyonic perturbations. In Figs. 12 and 13, we plot the growth of linear and non-linear perturbations in terms of a function of the scale factor. Then to describe the spherical collapse model in the framework of the fixed-potential tachyonic model, we need to solve Eqs. (25) and (26) for the FPTM state. We numerically solved these equations and plot the evolution of virialized overdensity Δ V , linear overdensity δ c , radius of the structure λ and ξ parameters in Figs. 14, 15, 16 and 17; respectively. These graphs show that in the tachyonic model, the δ c and ξ curves tend to Eds model faster than the Fig. 9 The behaviour of Eos parameter ω versus z for the ΛC DM and FPTM models Fig. 10 The variation of energy density parameter Ω(z) versus z for the ΛC DM and FPTM models Fig. 11 The variation of the Hubble dimensionless parameter E(z) versus z for the ΛC DM and FPTM models Fig. 12 The behaviour of the non-linear tachyonic perturbations δ l versus a for the ΛC DM and FPTM models quintessence and standard models. Therefore, we conclude that within the framework of the tachyon model, the pertur-

Conclusion
In this paper, we used the Chaplygin gas model to describe the spherical collapse of the large-scale structures. We solved the quintessence model with potential V (φ) = √ A(cosh 2 ( √ 3φ) + 1)/(2 cosh( √ 3φ)). Numerical analysis of this model showed that for 0.28 ≤ A ≤ 0.55, our model agrees with the current observational results. Then, we obtained the equation of state parameter for this model. As we expected, the EoS parameter of our model changed over time, and in the case of small redshifts, the standard cosmological model was favored. Also, using the parameter of the equation of state and the continuity equation, we obtained the evolution of the Hubble dimensionless parameter and the energy density of the scalar field. The results showed that within the framework of the Chaplygin gas model, during virialization, the structures take the return path faster than other models. Therefore, it can be concluded that in this model, the structures are formed faster than the other models. Then, using the numerical solution of the differential equations of growth of non-linear and linear perturbation, we investigated the structure formation parameters. In studying the structure formation of the spherical collapse model, we obtained parameters such as structure radius, virialized overdensity Δ V and ξ parameters. By considering the Chaplygin gas as an ordinary tachyonic scalar field, φ, one can obtain a corresponding potential for it. Finally, we considered the fixed-potential tachyon model as another candidate for the Chaplygin gas and calculated all the parameters of the spherical collapse for it. During virialization, the results confirmed that in the framework of the fixed potential tachyonic model, the structures take the return path faster than other models. Therefore, it can be concluded that in this model, the structures are formed faster than the other models. Therefore, we conclude that the Chaplygin gas model can properly explain the spherical collapse of the large-scale structures.
Data Availability Statement This manuscript has no associated data or the data will not be deposited. [Authors' comment: This is a theoretical study and no experimental data has been listed.] Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecomm ons.org/licenses/by/4.0/. Funded by SCOAP 3 .

Appendix A: Quintessence and tachyon equivalence to the Chaplygin gas
The simplest form of the equation of state for the Chaplygin gas is where A is a positive constant. Equation (38) is a particular form for the tachyon's equation of state with a constant potential. The general form of the equation of state for the tachyonic model is where φ denotes the tachyonic field. In Eq. (39), one can treat the potential as a constant, and find the EoS for the Chaplygin gas. Now, using Eq. (38) and the continuity equationρ + 3(ȧ/a)(ρ + P) = 0, we find where A and B are positive constants and a is the scale factor. One can see that for a << (B/A) 1/6 , the energy density ρ will be in the order of ∼ √ B/a 3 and when a is much larger than (B/A) 1/6 , ρ tends to −P √ A. Thus, one can conclude that in the past, when the scale factor was smaller, the Chaplygin gas can be considered as a pressureless matter. As we get closer to the present time, the scale factor grows, and the density of the Chaplygin gas becomes the same as the density of cosmological constant. Therefore, this fluid can be the source of undergoing accelerated expansion at the current time. Now, we will show that by choosing a particular potential in the quintessence model [30], we obtain the equation of state of the Chaplygin gas. In the quintessence model, the energy density of the scalar field and its pressure are written as ρ =φ 2 /2 + V (φ) and P =φ 2 /2 − V (φ). Using these equations, we finḋ Using Eqs. (38) and (40) we have: Applying the change of variable ∂ ∂t = −a H ∂ ∂a and using where κ equals to √ 8π G. Integrating the above equation, we obtain substituting Eq. (46) into Eq. (44), we find the proper form of the potential function as follows

Appendix B: The evolution of overdensity parameter
To investigate the cosmic fluid, we use the continuity, Euler, and Poisson equations, which describe the dynamic of the anḋ ρ + 3H (P +ρ) = 0, where v is velocity, and φ is gravitational potential. The physical coordinate r that is related to the comoving coordinate x via r = ax in which a is the scale factor. Since it is not possible to find a general solution to the non-linear equations, we write the quantities such as density, pressure, and gravitational potential as two parts of the background field and the perturbations around it. Writing these quantities in the comoving coordinates, the energy density, the equation of state, velocity, and the gravitational potential take the following forms and Using the quantities defined above, the equations of the dynamics of the fluid can be written aṡ and ∇ 2 φ − 4π G(1 + 3ω)a 2ρ δ = 0.
Taking the divergence from (57), we have where θ = ∇ x .u. Also, the stress tensor σ i j and the rotational tensor ω i j are defined as which can be used to define ω 2 = ω i j ω i j and σ 2 = σ i j σ i j . Now if we take derivative from Eq. (56) with respect to time and use Eqs. (56, 57, 58 and 59) , we obtain the following equation for the evolution of the spherical overdensitÿ δ + 2H −ω 1 + ω δ − 4 + 3ω 3(1 + ω)δ 2 1 + δ − 4π Gρ(1 + ω)(1 + 3ω)δ(1 + δ) Using the change of variable ∂ ∂t = −ȧ ∂ ∂a , we rewrite Eq. (62) as follows where E(a) = H/H 0 and ω m0 denotes the density parameter of matter in present. In the spherical collapse model, we can ignore the stress and rotational tensors. Also, we assume the matter is pressureless. So for the nonlinear regime, we obtain For the linear regime, the equation for the evolution of the spherical overdensity becomes