Accelerated cosmos in a non-extensive setup

Accelerated cosmos in a non-extensive setup

H., Alexander, Everton M. C. and Jorge Ananias Research Institute for Astronomy and Astrophysics of Maragha (RIAAM), Maragha 55134-441, Iran
Departamento de Física, Universidade Federal de Juiz de Fora, 36036-330, Juiz de Fora, MG, Brazil
Grupo de Física Teórica e Matemática Física, Departamento de Física, Universidade Federal Rural do Rio de Janeiro, 23890-971, Seropédica, RJ, Brazil

Here, we consider a flat FRW universe whose its horizon entropy meets the Rényi entropy of non-extensive systems. In our model, the ordinary energy-momentum conservation law is not always valid. By applying the Clausius relation as well as the Cai-Kim temperature to the apparent horizon of a flat FRW universe, we obtain modified Friedmann equations. Fitting the model to the observational data on current accelerated universe, some values for the model parameters are also addressed. Our study shows that the current accelerating phase of universe expansion may be described by a geometrical fluid, originated from the non-extensive aspects of geometry, which models a varying dark energy source interacting with matter field in the Rastall way. Moreover, our results indicate that the probable non-extensive features of spacetime may also be used to model a varying dark energy source which does not interact with matter field, and is compatible with the current accelerated phase of universe.

I Introduction

The violation of energy-momentum conservation law in curved spacetime has firstly been proposed by P. Rastall to modify general relativity theory of Einstein (GR) rastall (). After his pioneering work, various type of modified gravity in which matter fields are non-minimally coupled to geometry have been proposed cmc (); cmc1 (); cmc2 (); genras (); PRL (); shabprl (); shabprl1 (). It has been shown the Rastall correction term to the Einstein field equations can not describe dark energy meaning that a dark energy-like source is needed to model the current phase of universe expansion in this framework prd (). But, if one generalizes this theory in a suitable manner, then the mutual non-minimal coupling between geometry and matter field may be considered as the origin of the current accelerating phase and the inflation era genras (). More studies on Rastall theory can be found in Refs. gc1 (); al1 (); al2 (); al3 (); smal (); Rlag (); obs1 (); rastbr (); rascos1 (); rasch (); more1 (); more2 (); more3 (); more4 (); more5 (); hm (); msal (); plb1c (); msh (); clm (); net ().

In Einstein gravity, horizons may meet the Bekenstein-Hawking entropy-area law which is a non-extensive entropy rn1 (); rn2 (); rn4 (); rn5 (); rn6 (); rn7 (); rn8 (); rn9 (); rn10 (); rn11 (); rn12 (); rn13 (); rn14 (); rn15 (); rn16 (); rey2 (). Moreover, it has recently been argued that a deep connection between dark energy and horizon entropy may exist in gravitational theories cana (); cana1 (); mmg (); em (); mitra (); mitra1 (); mitra2 (); mitra3 (); mms (); md (); msgj (); mswr (). Indeed, although extensive statistical mechanics and its corresponding thermodynamics lead to interesting results about the universe expansion history mrej (), the mentioned points encourage physicists to use non-extensive statistical mechanics reyo1 (); tsallis () in order to study the thermodynamic properties of spacetime and its related subjects rn3 (); rey1 (); rsal1 (); rsal2 (); rsal (); anp (); ijtp (); cite5 (); Kom1 (); Kom2 (); Kom3 (); Kom (); EPJC (); MSK ().

Recently, applying Rényi entropy to the horizon of FRW universe and considering a varying dark energy source interacting with matter field, N. Komatsu found out modified Friedmann equations in agreement with observational data on the current phase of universe expansion EPJC (). Therefore, in his model, total energy-momentum tensor including the matter field and varying dark energy-like source is conserved. Combining this entropy with entropic force scenario, one can also obtain a theoretical basis for the MOND theory MSK (). In fact, the probable non-extensive features of spacetime may be considered as an origin for both the MOND theory and the current accelerated expansion phase in a universe filled by a pressureless source satisfying ordinary conservation law MSK (). Finally, it is useful to note here that both mentioned attempts EPJC (); MSK () used the Padmanabhan holographic approach Padd () in getting their models of the universe expansion.

Here, we are interested in obtaining a model for the universe dynamics by applying the Clausius relation as well as the Rényi entropy to the horizon of FRW universe which has non-minimally been coupled to matter field. Therefore, total energy-momentum tensor does not necessarily satisfy the ordinary conservation law, and in fact, it follows the Rastall hypothesis in our setup.

The paper is organized as follows. In the next section, introducing our approach, we present a thermodynamic description for Friedmann equations in Rastall theory. Using the Rényi entropy, our model of universe is obtained in Sec. (III). In the fourth section, we consider a universe filled by a pressureless source, and show that, in our formalism, it can experience an accelerated expansion. Sec. (V) includes the observational constraints of model. The last section is devoted to a summary and concluding remarks. We use the unit of in our calculations.

Ii Thermodynamic description of Friedmann equations in Rastall theory

Based on the Rastall hypothesis rastall ()


where denotes the Rastall constant parameter, and is the energy-momentum tensor of source which fills background. Moreover, is the Ricci scalar of spacetime. This equation says that there is an energy exchange between spacetime and cosmic fluids due to the tendency of geometry to couple with matter fields in a non-minimal way genras (); clm (). For example, the case, where is called the Rastall gravitational coupling constant, can support the primary inflationary era in an empty FRW universe genras (). In this manner, the ability of geometry to couple with matter fields in a non-minimal way generates a constant energy density equal to genras (). Some features of this non-minimal mutual interaction and its corresponding energy flux as well as its applications in cosmic eras have been studied in Ref. genras ().

Bearing the Bianchi identity in mind and integrating Eq. (1), one can reach at rastall ()


The Newtonian limit of the Rastall theory implies msal (); msh ()


where . Applying the thermodynamic laws to the horizon of spacetime and using Eqs. (2) and (3), it is shown that the horizon entropy is achieved as hm (); msal (); plb1c ()


in which is the horizon area. Combining this result with Eq. (3), one can easily find that entropy is positive whenever either meets the or condition msal (). This equation can also be written as , where is the Bekenstein-Hawking entropy msal (), meaning that the second law of thermodynamics is obtained for the mentioned values of , if it is satisfied by the Bekenstein-Hawking entropy. In addition, simple calculation lead to


for the Rastall constant parameter msal (). Indeed, Eqs. (3) and (5) are the direct results of imposing the Newtonian limit on Eq. (2), indicating that only for we have rastall (); msh (). In Ref. net (), assuming and studying Neutron stars in Rastall gravity, authors found out that is so close to zero. Therefore, we see that since they consider , their result does not reject the Rastall hypothesis. Finally, it is worth to remind here that the Bekenstein-Hawking entropy () is also obtainable at the appropriate limit of .

For a flat FRW universe with scale factor and line element


the apparent horizon, equal to the Hubble horizon, is located at


and therefore . Now, if the geometry is filled by a prefect fluid with energy density and pressure (), then Eq. (1) leads to


where dot denotes derivative with respect to time. It is also useful to remind here that for the flat FRW universe


Moreover, the use of Eq. (2) yields hm ()


The evolution of density perturbation in this model has been studied in Refs. prd (); rascos1 (); obs1 (). It has been shown that the story is the same as those of the standard cosmology at the background and linear perturbation level prd (). Finally, one can use Eqs. (10) to get hm ()


From the standpoint of tensor calculus, Eq. (2) is a solution for Eq. (1) leading to the above results. But, does thermodynamics lead to the same solutions? Indeed, since entropy is the backbone of thermodynamic approach, it is expected that Eq. (2) and thus the above results are available only whenever the horizon entropy meets Eq. (4). In order to find the Friedmann equations corresponding on Eq. (1) from the thermodynamics point of view, we consider the general form of entropy as . Additionally, one can use


to evaluate the energy flux crossing the apparent horizon Cai2 (); CaiKim (). Here, we also focus on an energy-momentum source as which yields for the work density. Finally, we see hm (); plb1c ()


Now, applying the Clausius relation () to horizon abe (), and using the Cai-Kim temperature () CaiKimt (), one can easily find


where pwang (). Now, combining this result with Eq. (8), we obtain


as the differential form of the first Friedmann equation. The result of this equation can be combined with Eq. (14) to get the second Friedmann equation.

Now, if the system entropy meets Eq. (4), then Eq. (14) leads to . It is also easy to insert Eq. (4) into Eq. (15) to get


where, is the integration constant of Eq. (15). Now, adding and subtracting from the LHS this equation, and using the relation, we can reach at


It is apparent that, for the case, the original Friedmann equations in the Rastall framework are reproduced (16). Indeed, since divergence of is zero, one may add the term to the RHS of Rastall field equations (2) to directly get Eq. (16) instead of Eq. (10). Anyway, we know that this term represents an unusual fluid in the context of ordinary physics leading to dark energy concept and thus its problems.

Now, Let us consider a flat FRW universe filled by a fluid with constant state parameter defined as . In this manner, calculations lead to hm ()


where is constant, for the energy density profile. For a universe filled by a pressureless component, where , combining this equation with Eqs. (16), one reaches at




Now, it is natural expectation that the matter density should be diluted during the expansion of universe. This limits us to the Rastall theories with and . Applying this result to Eq. (19), we find out that, at the long run limit (), we have , and therefore, universe may experience an accelerating phase. In this situation, Eq. (20) implies meaning that a non-minimal coupling between geometry and energy-momentum source cannot describe the current accelerating phase of universe in the Rastall framework. In fact, as it is apparent, we should have to get , and thus, an accelerating universe. Therefore, the same as the standard Friedmann equations, a dark energy-like source is needed to model the accelerating universe in the Rastall theory, a result in agreement with recent study by Batista et al. prd ().

In summary, our thermodynamic based study shows the dark energy problem is also valid in this formalism prd (), unless one generalizes the Rastall theory in a suitable manner genras ().

Iii Rényi entropy and Friedmann equations

Recently, Rényi entropy has been used in order to study the effects of probable non-extensive aspects of spacetime which led to interesting results in both cosmological and gravitational setups rey1 (); rey2 (); EPJC (); MSK (); rey03 (); rey04 (); rey05 (); rey06 (). For a non-extensive system including discrete states, the Rényi entropy is defined as reyo1 ()


in which and denote the probability of state and the non-extensive parameter, respectively. Moreover, the Tsallis entropy of this system is as follows tsallis ()


The linear relation between entropy and area is the key point of the Bekenstein-Hawking entropy (), which can also be obtained from the Tsallis’ non-additive entropy definition rn3 (). As it is apparent from Eq. (4), the functionality of with respect to the horizon area () is the same as that of the Bekenstein-Hawking entropy () meaning that may be considered as a special case of Eq. (22). More detailed studies on gravitational and cosmological implications of Tsallis entropy (22) can be found in Refs. rn3 (); rsal1 (); rsal2 (); rsal (); anp (); ijtp (); cite5 (); Kom1 (); Kom2 (); Kom3 () and references therein. Eq. (22) can be combined with Eq. (21) to show that


where , and we used to obtain this equation rey1 (); MSK (). It has frequently been argued that the Bekenstein-Hawking entropy is not an extensive entropy rn1 (); rn2 (); rn4 (); rn5 (); rn6 (); rn7 (); rn8 (); rn9 (); rn10 (); rn11 (); rn12 (); rn13 (); rn14 (); rn15 (); rn16 (); rey2 (), and in fact, the Bekenstein-Hawking entropy can be considered as a proper candidate for in gravitational and cosmological setups rey1 (); rey2 (); EPJC (); MSK (), a choice in full agreement with Ref. rn3 (). Here, following the above arguments and recent studies rn3 (); rey1 (); rey2 (); EPJC (); MSK (), we use Eq. (4) as the Tsallis entropy candidate in our calculations leading to


for the Rényi entropy of horizon () and its derivative with respect to (), respectively. Here, , and it is easy to check that whenever the non-extensive features of system approach zero (or equally ), the relation recovers Eq. (4). Now, inserting Eq. (25) into Eqs. (14) and (15), one can obtain




where is the integration constant, respectively. Defining a new constant , one can rewrite the last equation as


Since , Eq. (9) can be rewritten as


Finally, combining this equation with Eq. (28), and inserting the result into Eq. (26), we find




for the first and second Friedmann equations in our model, respectively. In fact, one should combine Eqs. (30) and (26) with each other, and then add and subtract the term to the result to find the last equation. In the above equations,

denote effective energy density and effective pressure in this framework, respectively. It is easy to see that, at the limit (or equally ), the results of Rastall theory are recovered. Additionally, we have whenever . Now, since , one can rewrite Eq. (31) as


where its LHS has the same form as that of the standard Friedmann equation roos (). Besides, simple calculations reach


for the acceleration equation. Therefore, we deal with two fluids. The first fluid is the ordinary energy-momentum tensor corresponding to the real fluid with energy density and pressure . The second fluid, which has geometrical origin, is called the effective energy-momentum tensor, and it is defined as


In fact, the obtained effective fluid consists of two parts: the non-extensive aspects of spacetime, and the non-minimal coupling between geometry and matter fields which follows the Rastall hypothesis. Therefore, the limit of only includes the non-extensive effects. We show it by in which


where and . In fact, it is a geometrical fluid originated from the non-extensive aspects of spacetime, and recovers the ordinary cosmological constant model of dark energy at the appropriate limit of . It is also easy to check that this source satisfies the conservation law i.e.


Therefore, acts as a time-varying dark energy model which satisfies the conservation law only for .

Bearing the Bianchi identity in mind, since the LHS of Eqs. (30) and (33) are compatible with the Einstein tensor, we should have leading to


and thus


In fact, the above results would also be obtained by writing Einstein field equations as . In addition, Eqs. (37) and (39) tell us that there is no energy flux between geometry and matter fields at the appropriate limit of (or equally ), a result in full agreement with Eq. (1). Indeed, although at the limit, since is a divergence-less tensor (37), we have and thus the ordinary energy-momentum conservation law is met by the source. Applying the and limits to the above equations, one can easily reach at the standard Friedmann equations compatible with the Bekenstein-Hawking entropy of horizon ijtp (); md (); em (). Therefore, the limit helps us in obtaining the modification of considering Rényi entropy to the standard Friedmann equations as


where and follow Eq. (III). It is worth to mention that, independent of the values of and , we have whenever . As a check, one can also insert in Eqs. (3) and (III) to get these results. In this manner, the acceleration equation is


and the second line of Eq. (40) can also be written as


where its LHS is in the form of the standard Friedmann equation roos (). Bearing Eq. (37) as well as the argument after Eq. (39) in mind, it is apparent that, for , the source respects the continuity equation i.e.


Although the same as Refs. EPJC (); MSK (), we used the Rényi entropy to get the modified Friedmann equations, these equations differ from those of recent studies EPJC (); MSK (). It has three reasons. i) While the entropy expression appears in acceleration equations obtained in Refs. EPJC (); MSK (), which use the Padmanabhan approach, its derivative is the backbone of getting the acceleration equation in our model based on applying thermodynamics laws to horizon (see Eq. (14)). ii) The Komar mass definition has been used by authors in Refs. EPJC (); MSK () only for the source. As we saw, our results are also obtainable if one writes the Einstein field equations as meaning that the Komar mass should be written for the modified energy-momentum tensor instead of . iii) In Ref. EPJC (), the non-extensive features of spacetime has been introduced as an origin for a time-varying dark energy () which interacts with matter fields and does not meet Eq. (37). This is while time-varying dark energy candidate of our model interacts with matter fields only in the Rastall way, and meets Eq. (37) in the absence of the Rastall hypothesis ().

Iv A Universe filled by a pressureless fluid

In order to study a universe filled by a pressureless source, we insert into Eq. (26) and use Eq. (30) to reach at


where . It is clear that, for a Rastall theory of or , the RHS of this equation and thus its LHS are vanished at long run limit () meaning that . Now, can be evaluated from


where is a constant. It is also apparent that, depending on the value of , this equation may be solvable even if we have or . Additionally, since is unknown parameter, this equation helps us in finding its possible values as a function of . One can also use this result in order to apply the limit to Eq. (26) to see that whenever meets either or . In summary, based on our results, the dark energy problem in Rastall theory can be overcame by considering the probable non-extensive features of spacetime.

Now, we consider a universe filled by a pressureless source () whenever . In this manner, both of Eqs. (18) and (43) lead to for energy density. Inserting into Eqs. (40), and following the recipe which led to Eqs. (44) and (45), we reach at




where , and . It also means that, whenever the divergence of is zero, the probable non-extensive features of spacetime, which behave as a conserved fluid, can be considered as the nature of current accelerating phase of universe if and meet the above equation. It is worth to note here that this equation is solvable even if . This result (the case) is in agreement with our previous results, where we found out the parameter should either meet or .

Bearing the results addressed after Eqs. (20) and (45) in mind, it is worth to mention that from the view point of dynamics, the and intervals are permissible for . On the other hand, thermodynamic considerations (the results of Eq. (4)), insist only the and intervals are admissible. Comparing these results with each other, one can easily find that and are common intervals. This means that the values of obtained from observations are allowed, if they be within in these ranges.

V Observational Constraints

In what follows, let us discuss the observational constraints on the scenarios presented above. In order to constrains the free parameters of the models we use: the Union sample 2012ApJ…746…85S (), which contains Supernovae type Ia (SNIa) in the redshift range , 36 Observational Hubble Data () in the range () compiled in 2015arXiv150702517M () and the Baryon Acoustic Oscillations (BAO) distance measurements at different redshift, in order to diminish the degeneracy between the free parameters.

v.1 Supernovae type Ia

In order to study the constraints applied to a cosmological model by the SNIa data, one can use the distance modulus defined as


where . Moreover, is the dimensionless Hubble parameter, and is the luminosity distance calculated as


in the flat FRW universe. Here, , and thus, using Eqs. (30) and (40), we can easily reach

whenever and (Model I), and

while and (Model II), respectively. In the above results, and is the critical density defined as . Now, applying the and conditions (the usual normalization conditions at ) to Eqs. (V.1) and (V.1), we get




in Model I and Model II, respectively. Therefore, Model I has three free parameters as , and Model II has two free parameters including . It is worthwhile to remind here that if and , then the model is recovered for .

Observational constraints on cosmological model can be obtained by minimizing given by metodo (); Nesseris & Perivolaropoulos (2005)




and we have marginalized over the nuisance parameter and .

v.2 Baryon Acoustic Oscillations (BAO)

The expanding spherical wave of baryonic perturbations, which comes from acoustic oscillations at recombination and co-moving scale of about , helps us in identifying the peak of large scale correlation function measured from SDSS (Sloan Digital Sky Survey). It is worth to note that the BAO scale depends on i) the scale of sound horizon at recombination, ii) the transverse and radial scales at the mean redshift of galaxies in the survey. In order to obtain the corresponding constraints on the cosmological models, we begin with for the WiggleZ BAO data 2011MNRAS.415.2892B () given as


where is data vector at , denotes the ordinary transpose, and is 2005ApJ…633..560E ()


in which


is the distance scale. Here, denotes the angular diameter distance defined as . Additionally, is the inverse covariance matrix for the WiggleZ data set given by


For the SDSS DR7 BAO distance measurements, can similarly be expressed as 2010MNRAS.401.2148P ()


where is the data points at and . is also defined as


in which is the radius of co-moving sound horizon given by




is the sound speed. Here, and = 2.726K. at the baryon drag epoch fitted with the formula, proposed in 1998ApJ…496..605E (),








is the inverse covariance matrix for the SDSS data set. Additionally, we use the Six Degree Field Galaxy Survey (6dF) measurement bao1 (), the Main Galaxy Sample of Data Release of Sloan Digital Sky Survey (SDSS-MGS) bao3 (), the LOWZ and CMASS galaxy samples of the Baryon Oscillation Spectroscopic Survey (BOSS-LOWZ) bao3 () and the distribution of the LymanForest in BOSS (BOSS - ) bao4 (). These measurements and their corresponding effective redshifts () are summarized in Table 1. Therefore, the total includes data point (for all the BAO data sets)

Survey z Parameter Measurement Reference

0.106 bao1 ()

0.57 bao3 ()

0.32 bao3 ()

2.36 bao4 ()

Table 1: Baryon acoustic oscillations (BAO) data measurements used in our statistical analysis.

v.3 History of the Hubble parameter

The differential evolution of early type passive galaxies provides direct information about the Hubble parameter . We adopt Observational Hubble Data (OHD) at different redshifts () obtained from 2015arXiv150702517M (), where data are deduced from the differential age method, and the remaining data belong to the radial BAO method. Here, we use these data to constrain the cosmological free parameters of the models under consideration. The corresponding can be defined as metodo ()


where is the theoretical value of the Hubble parameter at the redshift . This equation can be re-written as metodo ()


in which


The function depends on the model parameters. To marginalize over , we assume that the distribution of is a Gaussian function with standard deviation width and mean . Then, we build the posterior likelihood function that depends just on the free parameters , as




is a prior probability function widely used in the literature. Finally, we minimize with respect to the free parameters to obtain the best-fit parameters values.

v.4 Statistic analysis and results

Maximum likelihood , is the procedure of finding the value of one or more parameters for a given statistic which maximizes the known likelihood distribution. The maximum likelihood estimate for the best fit parameters is


and therefore, 2010arXiv1012.3754A (). In order to find the best values of the free parameters of the models, we consider


Moreover, the Fisher matrix is widely used in analyzing the constraints of cosmological parameters from different observational data sets 2009arXiv0901.0721A (); 2012JCAP…09..009W (). Having the best fit , the Fisher matrix can be calculated as


where depends on the uncertainties of the parameters for a given model. The inverse of the Fisher matrix also provides an estimate of the covariance matrix through . Its diagonal elements are the squares of uncertainties in each parameter marginalizing over the others, while the off-diagonal terms yield the correlation coefficients between parameters. The uncertainties obtained in the propagation of errors are also given by . Note that the marginalized uncertainty is always greater than (or at most equal to) the non-marginalized one. In fact, marginalization cannot decrease the error, and it has no effect, if all other parameters are uncorrelated with it. Previously known uncertainties of parameters, known as priors, can trivially be added to the calculated Fisher matrix.

Figure 1: Contour plots for the free parameter at and CL for Model I, from the joint analysis (left panel). Additionally, we present the corresponding marginalized one-dimensional posterior distributions for parameter (right panel).
Figure 2: Marginalized one-dimensional posterior distributions for parameter at and CL for Model II, from the joint analysis .

Table 2 summarizes the main results of the statistical analysis carried out by using the data sets SNIa, SNIa + BAO, and SNIa + BAO + for two scenarios including i) the Rényi entropy is taking into account within Rastall framework (Model I), and ii) the particular case of (Model II). The parameter takes into account the content of cold dark matter plus baryons to the present. It is useful to note that SNIa does not constrain very well, and in fact, results are improved by introducing the other observational tests including and . We can also see that the sign change of does not affect the main thermodynamic consideration obtained from Eq. (4). Indeed, since its obtained values meet the condition, entropy is always positive and dynamics, i.e. the acceleration of the universe, is in agreement with the observational data, an outcome in agreement with the results of previous section. For Model I, the likelihood contours arisen from the fitting analysis for the set of free parameters (, ), and marginalized one-dimensional posterior distributions for (PDF), considering the best fit values for used data sets , are presented in Fig. 1. Moreover, Fig. 2. includes marginalized one-dimensional posterior distributions for parameter in Model II. Here, we can appreciate slight deviations from the standard model (or equally the limit). This possibility does not rule out the standard model, and may in principle be used to distinguish between the and our models.

The equation of state (EoS) considering a given cosmological model can be written as


which has been derived from the combination of Eqs. (30) and (31), for the general case and Eq. (40) for the particular case of . Fig. 3. shows the behavior of the total EoS for both cases, with error propagation at 68.27% CL regarding the best fit values presented in table 2. We note that the total EoS, due to the mechanism presented in this paper, does not cross the phantom division line for the best fit of parameters. To high redshift approach asymptotically to a value of zero, that is, behaving like a fluid without pressure. In general, we note the behavior from the best fit values. Similar behavior are found in unification models in the dark sector of the universe.

Figure 3: Reconstruction of the EoS at 68% CL (gray region) from our joint analysis Model I (top panel) and Model II () (bottom panel). The blue line represents the best fit value for all data set .

On the other hand, it is natural to describe the kinematics of the cosmic expansion through the Hubble parameter , and its dependence on time, i.e. the deceleration parameter . The deceleration parameter is defined as combined with the relation, where , to get


From Eq. (41) and by considering along with , it follows that


where its RHS has to be evaluated at , and we have defined . In general, if is sufficiently large (i.e. ), then , which corresponds to an accelerated expanding universe.