Extensions of modified Chaplygin gas from Geometrothermodynamics

We derive modified classes of Chaplygin gas by using the formalism of Geometrothermodynamics. In particular, our strategy gives us extended versions of Chaplygin gas, providing a novel thermodynamic explanation. Thus, we show that our models correspond to systems with internal thermodynamic interaction. Bearing this in mind, we find new free parameters which are derived from thermodynamics and we give them an interpretation. To this end, we predict the range of values that every term can take in the context of homogeneous and isotropic universe. We also show that our new versions of modified Chaplygin gas can be interpreted as unified dark energy models, independently from the introduction of our new additional terms. Finally, we compare our theoretical scenarios through a fit on a grid based on the Union 2.1 compilation and we evaluate the growth factor of small perturbations. In this respect, we show that our model better adapts to the theoretical Λ\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\varLambda $$\end{document}CDM value, namely γΛCDM=611\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\gamma _{\varLambda CDM}=\frac{6}{11}$$\end{document}, than previous versions of modified Cgaplygin gas. We show numerical constraints at late and early redshift domains, which turn out to be compatible with previous results on standard versions of Chaplygin gas models.


Introduction
A growing number of data surveys indicate that the observable universe is undergoing a phase of accelerated expansion [1][2][3][4][5][6][7][8][9][10]. The most popular explanation for the cosmic speed up states that the fluid pushing up the universe to accelerate is due to a vacuum energy cosmological constant [11]. a e-mail: hbenaoum@sharjah.ac.ae b e-mail: orlando.luongo@lnf.infn.it c e-mail: quevedo@nucleares.unam.mx Alternatively, non-equilibrium barotropic fluids may originate dark energy, representing an alternative paradigm which does not make use of the cosmological constant. Among all the possible mechanisms of dark energy production, unified dark energy models have reached a large consensus during last years [12]. The underlying philosophy of these approaches consists in assuming a single-component fluid, which behaves as dust and dark energy at different regimes. A famous prototype, which goes toward unified scenarios, is dubbed Chaplygin gas, which has definitively attracted great interest during last years [13][14][15][16][17][18][19][20][21][22][23][24][25][26][27][28][29][30][31][32]. Chaplygin proposed his equation of state to study the lifting forces on a plane wing in aerodynamics, which has been modified accordingly to recent variants like p = Aρ − B ρ α , where A, B and α are universal positive constants [33]. As A → 0 and α → 1, one recovers the original formulation made by Chaplygin. Alternative models with A = 0 and α > 0 are known as generalized Chaplygin gases [34,35]. Unified dark energy models and, in particular, the modified Chaplygin gas has the remarkable property of describing the dark sector of the universe through a single dark fluid, which passes from a matter-dominated era to a cosmological constant-dominated era [36][37][38]. Although several papers discuss the behavior of the modified Chaplygin gas to reconcile the standard cosmological model with observations [39][40][41][42], a full explanation of the mechanism behind the origin of unified dark energy model is still lacking. In other words, if dark energy is not under the form of a genuine vacuum energy fluid, what is the mechanism enabling matter to behave anti-gravitationally? 1 In this work, we demonstrate that a possible solution to the aforementioned issue may be found in Geometrothermodynamics (GTD), which represents a formalism that has been developed during the past few years to describe ordinary thermodynamics by using differential geometry [46]. This formalism has encoded diverse applications in black hole thermodynamics [47], relativistic cosmology [48], mathematical chemistry [49], and others. Here, we show that it is possible to get the modified Chaplygin models by means of GTD and we obtain a new modification of the Chaplygin gas, with additional free parameters. Each parameter holds its own physical significance since it comes from GTD considerations and describes different thermodynamic properties. Hence, we show that the class of models we get is mathematically equivalent to previous versions of modified Chaplygin gas, albeit physically different. In particular, we combine thermodynamic and geometric considerations in the GTD framework, and find out a new equation of state of the modified Chaplygin gas, which explains the recent cosmological observations. In so doing, we give an explanation to the mechanism which produces unified dark energy models by using a pure thermodynamic approach. Further, we analyze the kinematics of the model by means of cosmography and fix limits over the new free constants, providing flat priors over them. Thus, we compare our paradigm with the standard cosmological model and show the main differences.
The work is organized as follows. In Sect. 2, we briefly review fundamental topics on GTD. Later, we show how to build up our cosmological model and discuss its thermodynamic and kinematic properties. In Sect. 3, we develop a numerical analysis making use of cosmic data surveys and compare our results with the ones obtained from the standard cosmological model. Finally, in Sect. 4, we discuss the conclusions and perspectives of our work.

From geometrothermodynamics to extended modified Chaplygin gas
The main purpose of this section is to show how to reproduce the dark energy contribution by means of GTD and how to handle it in view of the universe expansion history. In particular, we study the possibility that GTD can be seen as a method to generate unified dark energy models. To do so, we first analyze the whole framework of GTD as a differential geometric formalism for thermodynamics. Then, we show how to get from GTD the corresponding equation of state for the modified Chaplygin gas.

Basics of Geometrothermodynamics
The starting point of Geometrothermodynamics is a (2n+1)dimensional contact manifold T , called phase space, which is endowed with a Riemanian metric G and a contact 1-form Θ.
A set of coordinates {Z A } A=1,...,2n+1 = {Φ, E a , I a } a=1,...,n is introduced, where Φ represents the thermodynamic potential and E a and I a represent extensive and intensive thermodynamic variables, respectively. An important property of GTD is that all its geometric objects are constructed such that they are invariant with respect to Legendre transformations. In classical thermodynamics, this is equivalent to saying that the physical properties of a particular thermodynamic system do not depend on the choice of thermodynamic potential used for its description. A particular GTD metric G, which is invariant with respect to partial and total Legendre transformations, can be written as (summation over all repeated indices is implied) where the fundamental Gibbs 1-form is Θ = dΦ−δ ab I a d E b and Λ is a real constant, which can be set equal to one without loss of generality. The equilibrium n-dimensional submanifold E ∈ T is defined by the smooth map, such that the condition ϕ (Θ) = 0 is satisfied, implying that the first law of thermodynamics is satisfied on E. Applying the pullback ϕ to the metric G, we get the induced thermodynamic metric g given by: According to the GTD prescription, one only needs to specify the fundamental equation Φ = Φ(E a ) in order to find explicitly the metric g of the equilibrium submanifold E.

A new derivation of modified Chaplygin gas from GTD
To derive the modified Chaplygin gas, we choose the entropy S = S (U, V ) to be the thermodynamic potential and U , V to be the extensive variables. Then, the corresponding thermodynamic phase space T is a five dimensional space endowed with the set of independent coordinates Z A = {S, U, V, 1 T , P T }. The Gibbs' fundamental 1-form Θ S is given by: By defining the space of equilibrium states E by ϕ S (Θ S ) = 0, we obtain both the first law of thermodynamics, and the equilibrium conditions, To derive the modified Chaplygin gas, we consider the Nambu-Goto system of differential equations where is the d'Alembert operator and Γ A BC are the Christoffel symbols associated with the metric G AB of the phase space. The above system of differential equations follows from the variational principle δ E √ det(g)d n E = 0, which implies that the equilibrium space E with metric g constitutes an extremal subspace of the phase space T . A relevant solution of the Nambu-Goto differential equations is the fundamental equation where c 0 , c 1 , c 2 , α and β are real constants. Now, according to Eq. (2), the induced metric in the space of equilibrium states E is given as follows where the components of the thermodynamic metric g can be expressed as in which we used the auxiliary Using the induced metric, we obtain the scalar curvature for the particular case β = α as: where N (U, V ) and D (U, V ) are given in "Appendix A". The Legendre invariant scalar curvature is in general nonvanishing, which indicates the presence of internal thermodynamic interaction.

Extending the modified Chaplygin gas
From the fundamental Eq. (6), we obtain the the equilibrium conditions which lead to the equation of state This is a generalization of the equation of state proposed by hand by Chaplygin and of the other modifications investigated so far in the literature. Consequently, the GTD formulation is able to generate models of unified dark energy in terms of an equation of state, which contains the results of previous approaches and extensions of them. We will see that the generality of this new fundamental equation is constrained by a degeneration between parameters. In particular, it is not possible to constrain arbitrarily c 0 and c 1 , but the ratio c 0 c 1 only. Further, a multiplicative degeneration occurs among c 0 , c 1 , α and β. These degenerations suggest the existence of theoretical limits on the applicability of each particular case, as we will discuss below.

Limits on the modified Chaplygin gas and variable Chaplygin gas
Through the following recasting Equation (18) can be written as Consequently, we have two possibilities: β = α. This corresponds to a genuine modified Chaplygin gas. -β = α. This corresponds to the variable modified Chaplygin gas.
Assuming that ρ, as a function of time, can be represented as a function of the scale factor ρ(a), the conservation law can be integrated and so we assume the GTD equation of state to get the final solution for the Hubble rate. The models are formally equivalent to the ones provided in the literature. However, since each new parameter has its own physical interpretation, as we underline later in the manuscript, our approach differs from the ones of standard Chaplygin models. For each model, we can find its physical properties, investigating the corresponding effective thermodynamics and kinematics to check its goodness.

Thermodynamics
One of the advantages of knowing the fundamental Eq. (8) explicitly is that it can be used to derive the complete thermodynamic information of the system. In particular, two thermodynamic quantities are important for the cosmological model under consideration, namely, the temperature and specific heats.
Dark temperature: Equation (16) takes the form: Heat capacity: The arbitrary parameters which enter Eq. (22) must be chosen such that the temperature is always positive definite. Moreover, from the heat capacity one can infer the phase transition structure of the system. Indeed, for a phase transition to take place (C V → ∞), the condition must be satisfied.

Cosmological dynamics
We consider a homogeneous and isotropic universe pictured by the cosmological principle. The maximally symmetric spacetime which accounts for the cosmological principle is described by the Friedmann-Lemaître-Robertson-Walker (FLRW) line element where the whole dynamics is encapsuled into a(t), i.r. the scale factor. Spatial curvature is described by k = 0, ±1, respectively for spatially flat, closed or open universes. In agreement with recent Planck measurements, we assume k = 0 and we take into Einstein's equations the onecomponent perfect fluid through the requirement T μν = Pg μν + (ρ + P) u μ u ν . This leads to the most general Friedmann equations: and to the conservation law: where H ≡ȧ a is the Hubble parameter and c = 8π G = 1.
Under the hypothesis D = 4, we get: The final solution for the normalized Hubble rate thus becomes where E ≡ H/H 0 . It is evident from the above relation that one can assume c 1 = 1, without losing generality. Moreover, the limit to the concordance paradigm is found assuming the set: α → 0, β → 0, c 2 → Ω m,0 − 1, c 0 → 0.

Kinematics of our model
The kinematics of the universe consists in writing the scale factor in terms of observable quantities and to expand it around time, making use of the cosmological principle only. The corresponding definition of kinematic quantities is model-independent and furnishes a set of quantities which can be fixed by using cosmic data directly. Likely, the most suitable kinematic definitions are the deceleration and jerk parameters: or in terms of the Hubble rate At current time, these parameters enter the a(t) expansion in the form in which the subscript 0 indicates the value of the coefficients at present time, i.e., at z = 0. The deceleration parameter q provides the sign of the cosmic expansion, i.e., if the universe accelerates or decelerates. The sign of the jerk parameter j indicates whether the universe passes through a decelerated phase before acceleration or not. Current limits on q 0 and j 0 are essentially given by which are bounds useful for fixing constraints over our free constants defined in our extensions of the modified Chaplygin gas.
We compute the deceleration at present time and obtain and the jerk parameter which is reported in "Appendix A". We use the intervals defined in Eq. (37) to combine the constants c 0 , c 2 , α and β to show which limits might be imposed on them.

Bounds on the extra parameters
The aforementioned theoretical predictions define suitable bounds over the involved new coefficients, namely α, β, c 0 , c 1 and c 2 . In particular, we discuss below the physical priors which come from their definitions.
α > β As the universe increases its size, namely at current time, the right part of P(ρ) turns out to be negligible with respect to the left one. This implies that α − β > 0, i.e., α is always larger than β, although the sign of α will be fixed later. 1 + α > 0 The requirement that 1 + α > 0 guarantees that the asymptotic behavior of the specific heat at constant volume is not extensive with respect to the density ρ. Again, this does not fix the sign of α which is arbitrary inside the interval α > −1. Signs of α and c 2 The two parameters, α and c 2 have the same sign. This is a direct consequence of Eq. (24). In fact, since ρ and V are positive definite the signs of α and c 2 is the same. c 1 is arbitrary The sign and value of c 1 is unfixed. This is due to the degeneracy between c 1 and c 0 . Hence, for simplicity hereafter we fix c 1 = 1.
c 0 is negative definite The limit to the ΛCDM model guarantees that c 0 1+β 1+α < 0. In fact, the concordance paradigm is recover as α → 0, β → 0 with c 0 → −1 and B → 0. As the volume increases, however, the dominant term continues being the one ∝ ρ. Thus, conventionally one can assume −1 < c 0 < 0 without losing generality. Arbitrariness on c 2 The term c 2 A + 1+β 1+α is not bounded a priori. The only speculation that can be done is over the c 2 sign, which is arbitrary. As a consequence, we here take c 2 < 0, whose choice limits the intervals over α and β. In particular, by means of Eq. (37), we take α < 0 and β < 0, implying (1 + c 0 ) 1+β 1+α > 0. Bearing this in mind, we take 1 + β > 0, showing that α and β are closed to each other, but different. Positive T The temperature T is positive definite. This implies that c 2 is limited inside a particular set of values, albeit it is unbounded.
The above discussion means that only for this particular choices of the free parameters, a physical phase transition can occur and the Chaplygin scenarios are reproduced. Below we fix limits over the form of the free coefficients. To do so, we focus on H 0 , c 0 and c 2 .

Numerical constraints
We here test through our model with the Union 2.1 SNe of type Ia compilation. To do so, we employ the theoretical distance modulus, often parameterized by μ obs = m B − (M B −α X 1 +βC) ,, with absolute magnitude: Here, M host is the host stellar mass and {M, Δ M ,α,β} are the nuisance parameters. We here assume that the nuisance parameters do not influence the analysis and we neglect their effects. Thus, to perform our analysis we consider a fit on a grid, performed over a discretization of phase-space, having as likelihood the following in which we assume that all probabilities are independent. We thus assume uniform priors as given by (see 1), showing our numerical outcomes 2. For our joint analysis, we produce  in Fig. 1 the marginalized 2D 1σ and 2σ regions for the cosmological parameters.
The fitting procedure shows that at least at 1σ our model extends Chaplygin gas, obeying the properties emphasizing above for each free parameters. One may notice that the whole test, involving all free constants free to vary, suffers from a deep degeneration over the coefficients. Thus, a refined analysis would be useful to better underline the bounds over the entire sets of coefficients, in order to check the departures of our model from a genuine Chaplygin gas.

Perturbations and growth index
On sub-horizon scales, the first order perturbations of matter δ = δρ m ρ m are given by 2 This defines the evolution of δρ m , i.e. matter fluctuations in linear approximation. Immediately, it is possible to compute the evolution equation for the growth index. It is thus defined by: which has been obtained by replacing the time derivative by  a(t). Furthermore, we assume that a solution of Eq. (43) is given by: whose expression is given by: in which we have used the usual assumption: To insert the information coming from the Chaplygin equation of state, we need to consider the normalized matter, Ω m (a), and the corresponding equation of state. For the standard modified Chaplygin gas, with equation of state ω, we simply obtain [50]: Thus, in the case of a genuine modified Chaplygin gas, one argues that In our version of the Chaplygin gas, we can use the assumptions that we have made on the free constants coming from GTD. To do so, let us soon notice that in the case of our approach, the pressure becomes: in which we have plugged the behavior of the density in terms of the volume and the information on α and β on the exponents of the above relation. Approximatively this can be viewed as two evolving fluids whose net pressure evolves like P ∼ a −9 + a −6 . Clearly, as the volume becomes smaller, i.e. at the regime of small perturbations and high redshifts, we have that the term ∼ V −3 becomes dominant. This implies that our fluid behaves like the ωCDM model as the redshift increases. Plugging this expression in (44), we immediately get in which we have considered the zero-order term of γ in approximating it with our equation of state for the modified version of Chaplygin gas inferred from GTD. In order to check the goodness of our approach with previous paradigms, let us conventionally define a discrepancy as: whose aim is to evaluate how different models 3 depart from an arbitrary fiducial scenario, that here we assume to be the theoretical expectation coming from the standard cosmological model, i.e. γ ΛC DM . The discrepancies have been reported in Table 3.
The results have been obtained using Eq. (49) and γ ΛC DM = 6 11 , by means of numerical values of γ i respectively for the modified and variable modified Chaplygin gas got from [51], i.e. γ MCG = 0.562 and γ V MCG = 0.564. It soon comes evident that our model generalizes the standard Chaplygin approaches even at the level of small perturbations. The discrepancy between our approach and the fiducial one are mostly smaller than modified scenarios of Chaplygin fluids. The model is thus expected to work better at the level of small perturbations than previous Chaplygin scenarios. This is prevalently due to the fact that our model has been built up in terms of two counterparts: the first proportional to a quintessence-like fluid, namely P ∼ ωρ, whereas the second proportional to a variable modified Chaplygin gas. Thus the fluid has two possible regimes, at separate redshift domains, behaving differently at distinct epochs of the universe expansion. The second advantage of our approach is clearly due to the fact that the free constants entering the model are intertwined among them. By construction, in fact, A and B depend upon α, β and on c 0 , c 1 , c 2 . Even though this leads to a degeneracy problem as one fits the free constants, it "calibrates" the corresponding values to precise outcomes which turn out to be compatible with current numerical bounds on the growth index. It is, however, evident that a more refined analysis, which combines small and large scales, will better focus on the arbitrariness of c 2 . This would heal the degeneracy on coefficients, increasing the precision on the whole set of free constants, without fixing α and β to any theoretical values.

Final outlooks and perspectives
In this work, we employed the formalism of GTD to build up a new version of the modified Chaplygin gas describing the dark sector. The model shows up new physical constants which cannot escape a physical interpretation, since they represent a byproduct of our method. The advantage consists in handling pure extensive thermodynamic quantities to provide the equation of state, pushing up the universe to accelerate. Our scenario defines all the physical properties of the system by using thermodynamic standard laws, providing a fundamental GTD equation. Naturally, we obtain an equation of state which is then used to integrate the Friedmann equations of relativistic cosmology. By analyzing the physical properties of the resulting model, we conclude that from GTD it is possible to obtain fundamental equations for thermodynamic systems that can be used to develop physically reasonable cosmological models.
In particular, under certain conditions our formalism constructs a generalization of the modified Chaplygin gas cosmological model. To get this result, we inferred a particular fundamental equation that determines an equilibrium space embedded in the phase space by means of a map. This procedure satisfies the Nambu-Goto variational principle and relates the entropy of a thermodynamic system with its internal energy and volume. It turned out that our fluid corresponds to an equilibrium space with non-zero thermodynamic curvature, indicating the presence of internal thermodynamic interaction. The so-obtained pressure depends explicitly on the volume and has a non-vanishing small temperature associated to the dark fluid itself.
We showed the limiting cases of our approach with respect to the previous versions of modified Chaplygin gas and we performed a detailed analysis of its behavior. Dark energy and dark matter have been unified under our scheme by properly choosing the constants involved in the model. We discussed the role played by each constant and we finally gave priors over their evolution. As a final check, we forecasted a simple numerical analysis and showed that the limits over the free parameters of our model fulfill the basic requirements imposed by our theoretical priors. To do so, we assume the SNe data based on Union 2.1 catalog and baryonic acoustic observations. We thus discussed the main properties of our model in view of the cosmic evolution of each thermodynamic quantity.
Last but not least, we provided refined results on the growth factor γ in the regime of linear perturbations. We compared our expectations, got from the numerical values that we obtained through our fit, with the ones obtained from previous approaches. We demonstrated that, under the numerical bounds found in this work, the γ factor matches well the theoretical predictions of the standard ΛCDM model up to a ∼ 1% discrepancy. Future developments will focus on refined numerical analyses with additional data sets of SNe Ia and also with joint analyses making use of early time data sets. Moreover, we will clarify the role of each parameter even at the level of small perturbations, in order to check the goodness of our approach in view of structure formation. In particular, this will alleviate the arbitrariness of c 2 and will refine the degeneracy among all the other coefficients.